Water nanoparticles play an important role in atmospheric processes, yet their equilibrium and nonequilibrium liquid-ice phase transitions and the structures they form on freezing are not yet fully elucidated. Here we use molecular dynamics simulations with the mW water model to investigate the nonequilibrium freezing and equilibrium melting of water nanoparticles with radii R between 1 and 4.7 nm and the structure of the ice formed by crystallization at temperatures between 150 and 200 K. The ice crystallized in the particles is a hybrid form of ice I with stacked layers of the cubic and hexagonal ice polymorphs in a ratio approximately 2:1. The ratio of cubic ice to hexagonal ice is insensitive to the radius of the water particle and is comparable to that found in simulations of bulk water around the same temperature. Heating frozen particles that contain multiple crystallites leads to Ostwald ripening and annealing of the ice structures, accompanied by an increase in the amount of ice at the expense of the liquid water, before the particles finally melt from the hybrid ice I to liquid, without a transition to hexagonal ice. The melting temperatures T(m) of the nanoparticles are not affected by the ratio of cubic to hexagonal layers in the crystal. T(m) of the ice particles decreases from 255 to 170 K with the particle size and is well described by the Gibbs-Thomson equation, T(m)(R) = T(m)(bulk) - K(GT)/(R - d), with constant K(GT) = 82 ± 5 K·nm and a premelted liquid of width d = 0.26 ± 0.05 nm, about one monolayer. The freezing temperatures also decrease with the particles' radii. These results are important for understanding the composition, freezing, and melting properties of ice and liquid water particles under atmospheric conditions.

The phase behavior of confined water is a topic of intense and current interest due to its relevance in biology, geology, and materials science. Nevertheless, little is known about the phases that water forms even when confined in the simplest geometries, such as water confined between parallel surfaces. Here we use molecular dynamics simulations to compute the phase diagram of two layers of water confined between parallel non hydrogen bonding walls. This study shows that the water bilayer forms a dodecagonal quasicrystal, as well as two previously unreported bilayer crystals, one tiled exclusively by pentagonal rings. Quasicrystals, structures with long-range order but without periodicity, have never before been reported for water. The dodecagonal quasicrystal is obtained from the bilayer liquid through a reversible first-order phase transition and has diffusivity intermediate between that of the bilayer liquid and ice phases. The water quasicrystal and the ice polymorphs based on pentagons are stabilized by compression of the bilayer and are not templated by the confining surfaces, which are smooth. This demonstrates that these novel phases are intrinsically favored in bilayer water and suggests that these structures could be relevant not only for confined water but also for the wetting and properties of water at interfaces.

In order to predict accurately the chemical shift of NMR-active nuclei in solid phase systems, magnetic shielding calculations must be capable of considering the complete lattice structure. Here we assess the accuracy of the density functional theory gauge-including projector augmented wave method, which uses pseudopotentials to approximate the nodal structure of the core electrons, to determine the magnetic properties of crystals by predicting the full chemical-shift tensors of all 13 C nuclides in 14 organic single crystals from which experimental tensors have previously been reported. Plane-wave methods use periodic boundary conditions to incorporate the lattice structure, providing a substantial improvement for modeling the chemical shifts in hydrogen-bonded systems. Principal tensor components can now be predicted to an accuracy that approaches the typical experimental uncertainty. Moreover, methods that include the full solid-phase structure enable geometry optimizations to be performed on the input structures prior to calculation of the shielding. Improvement after optimization is noted here even when neutron diffraction data are used for determining the initial structures. After geometry optimization, the isotropic shift can be predicted to within 1 ppm.

Confinement of water usually depresses its melting temperature. Here we use molecular dynamics simulations to determine the liquid-crystal equilibrium temperature for water confined between parallel hydrophobic or mildly hydrophilic plates as a function of the distance between the surfaces. We find that bilayer ice, an ice polymorph in which the local environment of each water molecule strongly departs from the most stable tetrahedral structure, has the highest melting temperature (T(m)) of the series of l-layer ices. The melting temperature of bilayer ice is not only unusually high compared to the other confined ices, but also above the melting point of bulk hexagonal ice. Recent force microscopy experiments of water confined between graphite and a tungsten tip reveal the formation of ice at room temperature [K. B. Jinesh and J. W. M. Frenken, Phys. Rev. Lett. 101, 036101 (2008)]. Our results suggest that bilayer ice, for which we compute a T(m) as high as 310 K in hydrophobic confinement, is the crystal formed in those experiments.

Quasicrystals are structures with long-range order and no translational periodicity. Monatomic quasicrystals were predicted for model potentials, but no single-component atomic quasicrystal of an actual element has been reported to date. A dodecagonal quasicrystal was recently predicted to form in bilayer water. Water and silicon present striking similarities in their phase behavior, raising the question of whether quasicrystals may occur in silicon. Here, we show, using molecular simulations, that a confined silicon bilayer forms a quasicrystal upon compression between smooth surfaces. The quasicrystal is stable in a narrow region of the phase diagram and forms spontaneously upon cooling the liquid bilayer in a wide range of pressures. Cooling the liquid between atomically detailed plates incommensurate with the quasicrystal leads to its spontaneous formation at 1 atm of lateral pressure. This suggests that the silicon quasicrystal could be obtained in experiments at room pressure by tuning the structure and interactions of the surfaces.