Convergence of dissolving and melting at the nanoscale †

Phase transitions of water and its mixtures are of fundamental importance in physical chemistry, the pharmaceutical industry, materials sciences, and atmospheric sciences. However, current understanding remains elusive to explain relevant observations, especially at the nanoscale. Here, by using molecular dynamics simulations, we investigate the dissolution of sodium chloride (NaCl) nanocrystals with volume-equivalent diameters from 0.51 to 1.75 nm. Our results show that the dissolution of NaCl in aqueous nanodroplets show a strong size dependence, and its solubility can be predicted by the Ostwald – Freundlich equation and Gibbs – Duhem equation after considering a size-dependent solid – liquid surface tension. We ﬁ nd that the structure of dissolved ions in the saturated aqueous nanodropplet resembles the structure of a molten NaCl nanoparticle. With decreasing nanodroplet size, this similarity grows and the average potential energy of NaCl in solution, the molten phase and the crystal phase converges.


Introduction
5][6][7] For example, size-dependent solubilities have been found for sodium chloride (NaCl) and ammonium sulfate (AS) nanoparticles, which can be well explained by the Ostwald-Freundlich effect. 3As the particle size of NaCl decreases, increased solubility and dissolution rate have also been reported based on MD simulations. 8,9hen comes the question of what will happen if the particle size keeps decreasing, e.g., down to 1 nm or smaller?Will the solubility become innitely large?What's the state of these nanoparticles?So far, there are very few quantitative studies on size-dependent solubility of nanoparticles with a diameter of 2 nm or smaller.][12] Here, the size-dependent solubility of NaCl nanoparticles in water droplets is investigated for a diameter range from 0.507 to 1.745 nm by classical MD simulations.The size-dependent solid-liquid surface tension and structures under the corresponding critical conditions were also analyzed.In the MD simulations, an improved detection approach is developed to directly distinguish solid-like and liquid-like structures of NaCl nanoparticles in water droplets.With this approach, the solubility of a certain size of NaCl nanoparticles is determined.Then, solidliquid surface tension and solubility of NaCl nanoparticles at other sizes are calculated by the Gibbs-Duhem equation, the Ostwald-Freundlich equation and the Tolman equation.We also analyzed the structure of saturated aqueous nanodroplets of different sizes.In the end, the connection of dissolution and melting processes and its implications in the phase transition of aerosol particles is also discussed.

Results and discussion
Size-dependent solubility of NaCl nanoparticles Fig. 1 shows the solubility of NaCl nanoparticles with 7 to 108 pairs of NaCl ions calculated by the MD simulations.Here, we used the q 8 detection method, a local order parameter approach 29 to determine the solubility of NaCl (see Methods section).A strong size dependence of solubility is found for NaCl particles with diameters larger than 0.665 nm (13 pairs of NaCl ions) (red and orange bars in Fig. 1) according to the improved q 8 detection approach.4][15][16] The trend is qualitatively consistent with the experimental results (green triangles in Fig. 1). 3 The simulated solubility data are a bit lower compared to the experimental data, where the experimental saturation mass fraction of bulk NaCl (brown triangle in Fig. 1) is ∼1.5 times higher than the simulated saturation mass fraction of bulk NaCl.
Size-dependent solubility based on the simulations cannot be explained by the combination of the Ostwald-Freundlich and Gibbs-Duhem equations with a single bulk solid-liquid surface tension (s sl,bulk ) (eqn (3)).The prediction curve obtained based on the single s sl,bulk = 102 mN m −1 and the bulk saturation mass fraction ðc * s;bulk Þ calculated by MD simulations 17 (pink line in Fig. 1) does not t the data well.This indicates that the simulated solubility data become sensitive to s sl as the crystal gets smaller.Therefore, the size-dependent solid-liquid surface tension of NaCl nanoparticles should be considered to reproduce the sizedependent solubility of NaCl nanoparticles.
The saturation mass fraction ðc * s Þ of NaCl particles with diameters less than 0.665 nm does not vary monotonically.The saturation ratio ðc * s =c * s;bulk Þ remains similar when the diameter of the crystal decreases from 0.665 to 0.530 nm, with 13 to 8 pairs of NaCl ions.The saturation mass fraction suddenly jumps to ∼1 when the diameter of the crystal drops to 0.507 nm (7 pairs of NaCl ions), suggesting an innite solubility.This result is consistent with that in the melting process simulations (Fig. S9 †), where the NaCl particle with 7 pairs of NaCl ions is considered as molten at 300 K.
With a size-dependent solid-liquid surface tension (s sl ) based on the secondorder Tolman equation (ESI Note 3 †), [40][41][42] the size-dependent solubility from MD simulations can be well explained by the combination of the Ostwald-Freundlich and Gibbs-Duhem equations (blue shaded curves in Fig. 1).Compared to the tting curves from the rst-order Tolman equation (light pink shaded curves in Fig. 1), the blue tting curves from the second-order Tolman equation better capture the variation of NaCl solubilities.In addition, the size-dependent surface tension competes with the effect of crystal size on the solubility of nanoparticles.The reduction in particle size contributes to an increase in solubility, but this gain   is inhibited by the reduced surface tension, particularly when the particles are smaller than 0.665 nm with 13 pairs of NaCl ions.

Structures of NaCl nanoparticles at different dissolution states
Here, we examined the structures of NaCl nanoparticles at three states of dissolution: the incomplete, critical, and complete dissolutions.By incrementally adding water molecules and the q 8 criteria, we determined the range of water molecule quantities and corresponding solubility range for crystal dissolution.When the number of water molecules surrounding the crystal falls within the solubility range (between the upper and lower limits in Table S1 †), we consider it to be in a critical dissolution state.When the crystal is surrounded by fewer water molecules than the solubility range, we consider it to be in an incomplete dissolution state.Conversely, when the crystal is surrounded by more water molecules than the solubility range, we consider it to be in a complete dissolution state.At this stage, we can compare the distribution of NaCl ions and water molecules within the nanoscale particle under the three different states.We nd that the distribution of ions and water molecules is inhomogeneous in the incomplete and critical dissolution states, but close to homogeneous in the complete dissolution state.The droplet composed of 32 NaCl pairs (corresponding to 1 nm diameter) with 100 to 160 water molecules is taken as an example, where the particle containing 110-130 water molecules is in the critical dissolution state.
As shown by the light red lines in Fig. 2, in the incomplete dissolution state, both Na + and Cl − ions are primarily distributed within ∼1 nm distance from the center of mass (COM) of the particle, and more concentrated within distances of ∼0.4 nm.Na + ions have the highest density distribution at ∼0.3 nm distance from the COM, while Cl − ions are concentrated at the COM and ∼0.4 nm distance from the COM.Meanwhile, water molecules are distributed at distances >0.4 nm and mostly concentrated at distances of ∼1 nm from the COM.
As shown by the red and blue lines in Fig. 2, in the critical dissolution state, the concentration distribution of Na + and Cl − ions becomes similar but not homogeneous, and the concentration of ions decreases gradually from the COM to the particle surface.Water molecules are concentrated at a position closer to the COM compared to that in the incomplete dissolution state, and individual water molecules enter the interior of the particle (<0.4 nm distance from the COM).This shows that only a few water molecules are required to destroy the crystal structure inside the NaCl nanoparticle, where saturation can be much higher than that on the surface of the particle.
As shown by the light blue lines in Fig. 2, in the complete dissolution state, the distribution of Na + and Cl − ions is even more similar.This indicates that the distribution of ions and water molecules would be homogeneous when the mass saturation of NaCl particle is close to the bulk mass saturation.Note that the xed COM position of the total system in the MD algorithm might cause small uctuations of the nanocrystal position relative to the total COM, which may give rise to some smearing of the proles of the ion and water molecule densities near the COM of the total system.
The structure of different-size droplets in the critical dissolution state is further analyzed to quantify the characteristics of the inhomogeneous distribution of ions and water.The saturated droplet composed of 32 pairs of NaCl ions (corresponding to 1 nm diameter) with 130 water molecules is taken as an example.As shown in Fig. 3a, the prole of number density of ions and water molecules at different distances from the center of mass (COM) indicates ions are highly concentrated in the area around the COM of the droplet and are more diluted in the area near the surface.For example, the prole of the number of ions and water molecules at different distances from the COM in this droplet suggests only ∼4 water molecules exist in the area around the COM of the droplet with a radius of 0.5 nm, while there are 16 ions in the same area (Fig. 3b).
The two-dimensional density map of the relative fraction of ions shows more visually the gradient distribution of the ions (Fig. 3c).Except for a thin interface where no ion exists, the bulk solution with a concentration of 14.4 mol kg −1 is homogeneous.This is similar to the result from the previous simulation, where the smooth density prole of ions in the bulk solution with a concentration of 11.48 mol kg −1 indicates the homogeneous structure of the system. 18In contrast, the droplet composed of 32 NaCl and 130 water molecules shows inhomogeneity of the ion/water distribution with a gradient.The concentration becomes larger when the position is closer to the center.Furthermore, the distribution of Na + and Cl − ions is rather stratied.The bidimensional density map of the difference of the relative fraction of Na + and Cl − ions (Fig. 3d) shows that Na + ions are more concentrated in the center of the droplet, while Cl − ions are more concentrated at the surface of the droplet.

Dissolution and melting of NaCl nanoparticles
Following the approach of Qi et al., 19 we determined the melting temperature for the nanocrystals above the size of 15 pairs and the melting points of the NaCl nanocrystals show a clear size dependence (ESI Note 1).As shown in Fig. S8, † the heat capacity (blue circles) is obtained from the derivative of the average potential energy (red circles) with respect to the temperature (see eqn (S1) †), and the melting temperature of each crystal is at the temperature with the maximum apparent heat capacity.
Above the size of 15 pairs of NaCl nanocrystals, the depression of the melting point of the nanocrystals in comparison to the bulk is quite proportional to N (−1/ 3) , showing a signicant size dependence (red dashed line in Fig. S9 † and N is the number of ions).0][21] The dependence can be derived from the decrease in phase stability of the nanoparticles due to surface effects. 22By using the q 8 calculation, the melting point of the NaCl nanocrystals is also determined (blue diamonds in Fig. S9 †).The results show that above the size of 8 pairs of NaCl nanocrystals, the melting point depression of the nanocrystals is proportional to N (−1/3) (blue dashed line in Fig. S9 †).The approach of Qi et al. 19 (2001) cannot be applied to determine the melting temperature for NaCl nanocrystals below the size of 15 pairs.Note that the melting temperature for the NaCl nanocrystals below the size of 15 pairs is already less than the Debye temperature, and the heat capacity no longer follows the Dulong-Petit law.The vibrations of the atoms must be considered for accurate heat capacity data.
By comparing the average potential energy of dissolved and molten NaCl nanoparticles with diameters below 6 nm, we nd that the structure of dissolved NaCl in saturated solution becomes similar to the structure of molten NaCl when the NaCl particle size decreases.The average potential energy of NaCl (m NaCl ) only depends on the relative position of the ions, 14 so that it can reect the spatial distribution of the ions.The average potential energies in different simulated systems are compared to analyze the similarity of NaCl in the dissolved state and molten state quantitatively (Fig. 4).Dark circles represent m NaCl in the molten NaCl particles at the respective melting points, and the values are all close to −700 kJ mol −1 (grey line).This means that the spatial distribution of ions in a molten NaCl particle at the respective melting point is almost the same and independent of size.This phenomenon was also found for other substances, such as copper. 23However, m NaCl in the dissolved droplets at the critical dissolution state (blue circles) decreases with the decrease of size.The blue circles approach −700 kJ mol −1 and deviate from the blue dashed line (z−358.35± 20.95 kJ mol −1 ), which represents m NaCl in bulk saturated solution.It suggests the structural similarity of dissolved NaCl in saturated nanodroplets to molten NaCl at the melting point increases when the size of the NaCl particle decreases.
In Fig. 4, the red circles represent m NaCl of the NaCl nanoparticles at 300 K, which deviate from the red line that represents m NaCl of the bulk NaCl crystal (z−780 kJ mol −1 ).The increase of m NaCl with the decrease of particle size is because the bigger particle is more stable and thus has a higher melting point temperature.As the particle size decreases to ∼7 pairs of NaCl ions, the red circles and blue circles almost converge at around −700 kJ mol −1 .The convergence of m NaCl suggests that 7 pairs of NaCl ions at 300 K is in a molten state and can be considered as a saturated droplet of NaCl with innite solubility.Due to the sizedependent solubility, smaller molten NaCl particles at 300 K (#7 pairs of NaCl ions) can also be considered as a saturated droplet of NaCl with innite solubility.
The radius of gyration of the particles and the radial distribution function of ions to ions (RDF ion-ion ) also supports the analysis of m NaCl .The similarity of the radius of gyration of the NaCl particles in the dissolution and melting processes becomes more pronounced when the particle size decreases (Fig. 5).The similarity of the relative positions of ions in the dissolution and melting processes increases when the particle size decreases (Fig. S12 †).Harbury 4 gave a similar assumption that there is a signicant comparability between a supersaturated solution and undercooled melts.In summary, the dissolution of bulk NaCl solids is a traditional dissolution, while the dissolution of NaCl particles at the nanoscale is analogous to the melting process of NaCl particles at the nanoscale.These two different types of phase transition, dissolving and melting, "meet" at #7 pairs of NaCl ions at 300 K.
The radial distribution functions of Na + to Na + (RDF Na-Na ), Cl − to Cl − (RDF Cl- Cl ), Na + to Cl − (RDF Na-Cl ), and Cl − to Na + (RDF Cl-Na ) give valuable additional information that the q 8 -determined solubility and melting point of NaCl nanocrystals are on rm ground.The RDF Na-Na and RDF Cl-Cl of the q 8 -determined unmelted NaCl nanocrystals and undissolved NaCl in saturated nanodroplets show rather sharp peaks at the rst to fourth nearest neighbor distances of ∼0.40 nm, ∼0.56 nm, ∼0.69 nm, and ∼0.79 nm (Fig. S13 and S15 †).The RDF Na-Cl and RDF Cl-Na of the q 8 -determined unmelted NaCl nanocrystals and undissolved NaCl in saturated nanodroplets show the same sharp peaks at the rst to fourth nearest neighbor distances of ∼0.28 nm, ∼0.49nm, ∼0.63 nm, and ∼0.84 nm (Fig. S17 and S19 †).This is well compatible with the simple predictions from the snapshot of 9 pairs in a NaCl crystal shown in Fig. S20, † indicating that the nanocrystals are still crystalline despite their smallness, and the crystal structure is not disturbed by the condensation of a few water molecules.Meanwhile, the RDF Na-Na and RDF Cl-Cl as well as the RDF Na-Cl and RDF Cl-Na of the q 8 -determined molten NaCl nanocrystals at the melting points and dissolved NaCl in saturated nanodroplets show atter peaks compared to the ones of the unmelted and undissolved NaCl nanoparticles (Fig. S12, S14, S16 and S18 †).This captures well the different structural features between the solid and liquid states of NaCl, thus supporting the results of the solubility and melting point of NaCl nanoparticles determined by the q 8 method.Snapshots of NaCl nanoparticles during the melting and dissolution processes in Fig. S20 † support the solubility and melting point of NaCl nanoparticles determined by the q 8 method as well.The snapshots of NaCl nanoparticles at a temperature lower than the q 8 -determined melting points show that the typical shape of nano-sized NaCl nanocrystals is not spherical.Both in the case of melting and the dissolution process, a spherical shape occurs only when the transition to the uid phase takes place.

Implications
Our MD simulations demonstrate that the solubility of NaCl at the nanoscale is size-dependent and can be much higher than that of a bulk solution.The melting point of NaCl is also found to be size-dependent.To transform a NaCl particle with a temperature of 300 K into the liquid state, the temperature of the environment should be increased or water should be added to the NaCl.In these processes, heat from a higher-temperature environment or the interaction with water molecules changes the NaCl crystal to the liquid state.In other words, the heat, and the interaction with water for NaCl nanocrystals are equivalent.As the size decreases, a lower temperature or less surrounding water molecules are needed to change the phase of NaCl.The nanoparticle with only 7 pairs of NaCl ions cannot exist in the solid state at 300 K, thus no heat or interaction with a few water molecules is needed to transform it into liquid.
Overall, this study provides a new perspective to connect the dissolution and melting processes of NaCl nanoparticles, and this connection may also be suitable for other solids.For a substance whose melting point is higher than 300 K at the macroscopic level, there exists a critical diameter (D c ) at which the melting point is 300 K. Naturally, when the diameter of the nanoparticle composed by this substance is below D c , the particle is expected to be liquid at room temperature.The nanoparticle whose diameter is larger than D c is solid at room temperature, while a small amount of adsorbed water molecules can turn it into a molten-like liquid with a thin solution coating.5][26] This study mainly focuses on the thermodynamics of dissolution.A deep understanding of the NaCl dissolution mechanism also requires knowledge of the kinetics, and we would like to direct the readers to other related research. 27,28

MD simulations
4][15][16] Therefore, this force eld combination is applied in the study to simulate the NaCl-water mixture.The MD simulations are carried out with the GROMACS 2016.6 package. 29,30NaCl nanocrystals with different diameters are surrounded by certain water molecules initially.The systems are then energetically minimized by the steepest-descent method and equilibrated for 100 ps at 300 K.The MD simulations are further performed in the NVT ensemble with periodic boundary conditions and velocity-rescaling thermostat.Electrostatic interactions are calculated using the particle mesh Ewald (PME) algorithm, and van der Waals interactions are accounted for up to a cutoff distance of 10 Å.All simulations are carried out for at least 600 ns using a 1 fs time step, and conformations for analysis are saved every 2 ps.To determine the solubility, we kept increasing the number of water molecules in the system until we found a critical number that allows NaCl to dissolve completely based on the q 8 algorithm.Ideally, we would add water molecules one by one, but in practice, due to the limitations of the computational resources, we increased the step size by 5 to 20 molecules.Thus, instead of an exact number, we determine a narrow range of water molecules required for crystal dissolution and Table S1 † summarizes the upper and lower limits of this range.We chose a simulation duration of at least 600 ns, ensuring that even the largest simulated NaCl nanoparticle (108 pairs of NaCl ions with 650 water molecules, ∼1.745 nm diameter) would completely dissolve (Fig. S10 †).Each simulation box has a uniform size of 10 × 10 × 10 nm 3 , ensuring it can accommodate the largest simulated NaCl nanoparticle (108 pairs of NaCl ions with 650 water molecules, ∼1.745 nm diameter).More detailed information about the number of ions/ water molecules, the initial structure setup, and simulation times are summarized in Table S1 and Fig. S11.† In this study, the melting processes of NaCl nanocrystals is also simulated.The JC force eld is adopted to keep synchronous with the simulations of the NaCl dissolution process.The particles with different diameters are simulated at different temperatures to determine the melting point.More details about the setup are shown in ESI Note 1. † Solubility determination: q 8 detection approach To detect the dissolution and melting progress of NaCl nanocrystals, it is important to distinguish ions that are in the solid-like and liquid-like structures.To do this, the local order parameter approach of Steinhardt et al. 31 is adopted.For each ion (i), the local bond-orientational order parameter q l (i) is calculated as: and where l is a free integer parameter and m is an integer that runs from m = −l to l. Y l m is the spherical harmonic, r ij is the position vector of the neighbor ion (j) with respect to the central ion (i), q(r ij ) and B(r ij ) are the polar and azimuthal angles with respect to the reference coordinate of r ij .N is the number of neighboring ions.Note that in the perfect bulk NaCl crystal structure, each ion has 6 nearest neighbors of different species (e.g., Na + has 6 neighboring Cl − ions and vice versa) and 12 next-nearest neighbors of the same kind (e.g., Na + has 12 neighboring Na + ions).However, the NaCl nanoparticle structures during the melting and dissolution processes in the simulations are mostly not perfect crystal structures.Therefore, here the ions are considered to be identical and the number of neighboring ions is equal to 12. Lanaro and Patey 32 utilized a similar method to track the production of NaCl nuclei in aqueous solutions.Jiang et al. 33 also used this approach to investigate the mechanism of NaCl crystal nucleation from solutions with high supersaturations.The parameter q l (i) is sensitive to different crystal symmetries depending on the choice of l. q l (i) distributions for l = 4, 6, and 8 were examined, using the cubic nanocrystal with 108 pairs of NaCl in the solid phase and the dissolved state as representative of the liquid phase at 300 K.As shown in Fig. S3, † the local bond-orientational order parameter q 8 (i) provides the best separation between solid and liquid phase distribution for NaCl nanoparticles.
The NaCl nanocrystal is considered as dissolved or molten when the time average of the q 8 value of each ion in the crystal is less than 0.35.q 8 is the local bond-orientational order parameter to provide the best separation between solidlike (>0.35) and liquid-like (#0.35) states of each NaCl ion in the nanoparticles (Fig. S5 †).We added 5 water molecules each time for the certain NaCl nanocrystal in the simulation.The crystal with N + 5 water molecules is considered as dissolved when the time average of the q 8 value of each ion in the crystal is less than 0.35.At this point, we consider the solubility range of the crystal to be within the range of water molecule quantities of N to N + 5. Furthermore, based on the known molar masses of NaCl ions and water molecules, we can calculate the corresponding range of saturated mass fraction values, which represents the solubility range.
Since at 300 K the crystal structure of NaCl is already close to perfectly rigid, a simple study of the radial pair distribution functions between the different pairs of ions can clearly identify which NaCl clusters are still in a nanocrystal state and which are already in a liquid state (Fig. S6 †).

Size-dependent solubility estimation
Size-dependent solubility can be determined by the combination of the Ostwald-Freundlich 6,7 and Gibbs-Duhem equations [34][35][36] as: where R is the universal gas constant, T is the temperature, n is the molar volume of solid phase NaCl (atomic volume) given as 2.7 × 10 −5 m 3 mol −1 , and M is the molar weight of solute.x * s;bulk is the mass fraction of solute in a saturated bulk solution, and a * s;bulk is the solute activity at x * s;bulk .x * s is the mass fraction of solute in a saturated droplet of diameter D. a w is the water activity retrieved from the modied Tang-Munkelwitz (TM) model (ESI Note 2 †). 3,38s sl is the size-dependent solid-liquid surface tension estimated by the Tolman equation (ESI Note 3 †). 39In this study, D is the volume equivalent diameter of a NaCl crystal at the saturation dissolution point, which is equal to the diameter of the NaCl droplet of equivalent volume, and the volume of the NaCl crystal (V) is measured based on the initial structure in the MD simulations, where the distance between Na + and Cl − is 0.28 nm.

Structure analysis
The proles of number density and number of ions and water molecules at different distances from the center of mass (COM) are used to analyze the structure of saturated droplets.The relative fraction of ions is calculated by the following equation: RFI ¼ r ion:area =r water:area r ion:total =r water:total (5)   where r ion.area represents the number density of NaCl in a spherical shell area, r water.arearepresents the number density of water in a spherical shell area, the thickness of the spherical shell is 0.02 nm.r ion.total represents the number density of NaCl in the entire droplet, r water.totalrepresents the number density of water in the entire droplet.The radial distribution functions (RDFs) of ions with respect to other ions can reect the spatial distribution of ions. 29,37They are calculated as: Here, g ion-ion (r) represents the RDF for one conformation: r represents a given distance and the interval value is the thickness of a spherical shell equal to 0.02 nm, N ion represents the number of ions for the given conformation, i and j represent the ions and i s j, r ij represents the distance from ions i and ions j, n(r ij ) represents the number of ions j locates in the distance of r − 0.02 nm to r from the given reference ion i. g ion-ion (r,t o ) represents the average value of RDF from t o to t o + 1 ns.As the simulation time step is 0.002 ns, g ion- ion (r,t o ) can be considered as the average value of 500 conformations, g ion-ion (r,t o + t) is the RDF at the time t o + t.In this study, RDF ion-ion represents g ion-ion (r,t o ).The prole of RDF ion-ion over time can reect the change of relative positions of ions during the process of dissolving.

Average potential energy calculation
To calculate the average potential energy of dissolved or molten NaCl, the trajectories only containing the Na + and Cl − ions aer the MD simulations of a given NaCl-water system are extracted and then rerun for the new trajectories.
In the rerun simulation, energy is calculated every 2 ps, and the average value over the time aer the equilibrium of a system is taken to describe the system.

Radius of gyration
The radius of gyration is dened as: where N is the number of ions in the cluster, R i is the position of ion (i) and R COM is the position of the center of mass of the cluster.

Conflicts of interest
Authors declare that they have no competing interests.

Fig. 1
Fig. 1 Size-dependent solubility (saturation mass fraction) of NaCl nanoparticles at 300 K.The solubility range of NaCl nanoparticles is obtained from the MD simulations for particles with 8 to 108 pairs of NaCl ions (red and orange bars) and 7 pairs of NaCl ions (red circle).a The size-dependent solubility of NaCl nanoparticles is taken from experiments for particles with dry diameters of 5.18 to 54.2 nm (green triangles).b,c The NaCl bulk solubility from MD simulations (grey circle) and from experiment (brown triangle) is taken from ref. 9.The size-dependent solubility curves are obtained based on the combination of the Ostwald-Freundlich and the Gibbs-Duhem equation with and without considering size effects on surface tension.The size effect on solid-liquid surface tension is considered based on the 2nd order Tolman equations for the light blue shaded curves.d The size effect on solid-liquid surface tension is considered based on the 1st order Tolman equation with the bulk solid-liquid surface tension from ref. 17 for the light pink shaded curves.e The green curve without considering the size effect of surface tension is based on the NaCl bulk solubility obtained experimentally and the bulk solid-liquid surface tension from ref. 3. f The pink curve without considering the size effect of surface tension is based on the NaCl bulk solubility from MD simulations and the bulk solid-liquid surface tension from ref. 17.

17
Fig. 1 Size-dependent solubility (saturation mass fraction) of NaCl nanoparticles at 300 K.The solubility range of NaCl nanoparticles is obtained from the MD simulations for particles with 8 to 108 pairs of NaCl ions (red and orange bars) and 7 pairs of NaCl ions (red circle).a The size-dependent solubility of NaCl nanoparticles is taken from experiments for particles with dry diameters of 5.18 to 54.2 nm (green triangles).b,c The NaCl bulk solubility from MD simulations (grey circle) and from experiment (brown triangle) is taken from ref. 9.The size-dependent solubility curves are obtained based on the combination of the Ostwald-Freundlich and the Gibbs-Duhem equation with and without considering size effects on surface tension.The size effect on solid-liquid surface tension is considered based on the 2nd order Tolman equations for the light blue shaded curves.d The size effect on solid-liquid surface tension is considered based on the 1st order Tolman equation with the bulk solid-liquid surface tension from ref. 17 for the light pink shaded curves.e The green curve without considering the size effect of surface tension is based on the NaCl bulk solubility obtained experimentally and the bulk solid-liquid surface tension from ref. 3. f The pink curve without considering the size effect of surface tension is based on the NaCl bulk solubility from MD simulations and the bulk solid-liquid surface tension from ref. 17.
Fig. 1 Size-dependent solubility (saturation mass fraction) of NaCl nanoparticles at 300 K.The solubility range of NaCl nanoparticles is obtained from the MD simulations for particles with 8 to 108 pairs of NaCl ions (red and orange bars) and 7 pairs of NaCl ions (red circle).a The size-dependent solubility of NaCl nanoparticles is taken from experiments for particles with dry diameters of 5.18 to 54.2 nm (green triangles).b,c The NaCl bulk solubility from MD simulations (grey circle) and from experiment (brown triangle) is taken from ref. 9.The size-dependent solubility curves are obtained based on the combination of the Ostwald-Freundlich and the Gibbs-Duhem equation with and without considering size effects on surface tension.The size effect on solid-liquid surface tension is considered based on the 2nd order Tolman equations for the light blue shaded curves.d The size effect on solid-liquid surface tension is considered based on the 1st order Tolman equation with the bulk solid-liquid surface tension from ref. 17 for the light pink shaded curves.e The green curve without considering the size effect of surface tension is based on the NaCl bulk solubility obtained experimentally and the bulk solid-liquid surface tension from ref. 3. f The pink curve without considering the size effect of surface tension is based on the NaCl bulk solubility from MD simulations and the bulk solid-liquid surface tension from ref. 17.

Fig. 2
Fig. 2 The number density profiles of ions and water molecules in particles composed of 1 nm NaCl (32 pairs) and different numbers of water molecules for (a) Na + ions, (b) Cl − ions, (c) ions, and (d) water molecules.Particles with less than 110 water molecules are in the incompletely dissolved state (light red lines).Particles with 110-130 water molecules are in the critically dissolved state (red and blue lines), and particles with more than 130 water molecules are saturated NaCl aqueous nanodroplets (light blue lines).

Fig. 3
Fig. 3 Distribution of ions and water molecules in the droplet composed of 1 nm NaCl and 130 water molecules.(a) The number density of ions and water molecules at different distances from the COM in the droplet.(b) The cumulative number of ions and water molecules at different distances from the COM in the droplet.(c) Two-dimensional density map of the relative fraction of ions.The contour lines represent the relative fraction of water molecules.(d) Two-dimensional density map of the difference of the relative fraction of Na + and Cl − ions.The contour lines represent the relative fraction of water molecules.

Fig. 4
Fig. 4 Average potential energy of NaCl (m NaCl ) in solution, molten NaCl and crystal.Blue circles represent m NaCl in the saturated droplets at the critical dissolution state.The blue dashed line represents m NaCl in bulk saturated solution (z−358.35± 20.95 kJ mol −1 ).Black circles represent m NaCl in the molten NaCl at the respective melting points, and the values are all near −700 kJ mol −1 (gray dashed line).Red circles represent m NaCl in the NaCl crystals at 300 K, and the dashed red line is m NaCl in the bulk NaCl crystal (z−780 kJ mol −1 ).