Gabriele C. Sosso

Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations

The nucleation of crystals in liquids is one of nature’s most ubiquitous phenomena, playing an important role in areas such as climate change and the production of drugs. As the early stages of nucleation involve exceedingly small time and length scales, atomistic computer simulations can provide unique insights into the microscopic aspects of crystallization. In this review, we take stock of the numerous molecular dynamics simulations that, in the past few decades, have unraveled crucial aspects of crystal nucleation in liquids. We put into context the theoretical framework of classical nucleation theory and the state-of-the-art computational methods by reviewing simulations of such processes as ice nucleation and the crystallization of molecules in solutions. We shall see that molecular dynamics simulations have provided key insights into diverse nucleation scenarios, ranging from colloidal particles to natural gas hydrates, and that, as a result, the general applicability of classical nucleation theory has been repeatedly called into question. We have attempted to identify the most pressing open questions in the field. We believe that, by improving (i) existing interatomic potentials and (ii) currently available enhanced sampling methods, the community can move toward accurate investigations of realistic systems of practical interest, thus bringing simulations a step closer to experiments.

Fast Crystallization of the Phase Change Compound GeTe by Large-Scale Molecular Dynamics Simulations

Phase change materials are of great interest as active layers in rewritable optical disks and novel electronic nonvolatile memories. These applications rest on a fast and reversible transformation between the amorphous and crystalline phases upon heating, taking place on the nanosecond time scale. In this work, we investigate the microscopic origin of the fast crystallization process by means of large-scale molecular dynamics simulations of the phase change compound GeTe. To this end, we use an interatomic potential generated from a Neural Network fitting of a large database of ab initio energies. We demonstrate that in the temperature range of the programming protocols of the electronic memories (500-700 K), nucleation of the crystal in the supercooled liquid is not rate-limiting. In this temperature range, the growth of supercritical nuclei is very fast because of a large atomic mobility, which is, in turn, the consequence of the high fragility of the supercooled liquid and the associated breakdown of the Stokes-Einstein relation between viscosity and diffusivity.

The Many Faces of Heterogeneous Ice Nucleation: Interplay Between Surface Morphology and Hydrophobicity

What makes a material a good ice nucleating agent? Despite the importance of heterogeneous ice nucleation to a variety of fields, from cloud science to microbiology, major gaps in our understanding of this ubiquitous process still prevent us from answering this question. In this work, we have examined the ability of generic crystalline substrates to promote ice nucleation as a function of the hydrophobicity and the morphology of the surface. Nucleation rates have been obtained by brute-force molecular dynamics simulations of coarse-grained water on top of different surfaces of a model fcc crystal, varying the water-surface interaction and the surface lattice parameter. It turns out that the lattice mismatch of the surface with respect to ice, customarily regarded as the most important requirement for a good ice nucleating agent, is at most desirable but not a requirement. On the other hand, the balance between the morphology of the surface and its hydrophobicity can significantly alter the ice nucleation rate and can also lead to the formation of up to three different faces of ice on the same substrate. We have pinpointed three circumstances where heterogeneous ice nucleation can be promoted by the crystalline surface: (i) the formation of a water overlayer that acts as an in-plane template; (ii) the emergence of a contact layer buckled in an ice-like manner; and (iii) nucleation on compact surfaces with very high interaction strength. We hope that this extensive systematic study will foster future experimental work aimed at testing the physiochemical understanding presented herein.

Vibrational properties of crystalline Sb2Te3 from first principles

Phonon dispersion relations and infrared and Raman spectra of crystalline Sb(2)Te(3) were computed within density functional perturbation theory. Overall good agreement with experiments is obtained, which allows us to assign the Raman and IR peaks to specific phonons.

Neural network interatomic potential for the phase change material GeTe

GeTe is a prototypical phase change material of high interest for applications in optical and electronic nonvolatile memories. We present an interatomic potential for the bulk phases of GeTe, which is created using a neural network (NN) representation of the potential-energy surface obtained from reference calculations based on density functional theory. It is demonstrated that the NN potential provides a close to ab initio quality description of a number of properties of liquid, crystalline, and amorphous GeTe. The availability of a reliable classical potential allows addressing a number of issues of interest for the technological applications of phase change materials, which are presently beyond the capability of first-principles molecular dynamics simulations.

Breakdown of Stokes–Einstein relation in the supercooled liquid state of phase change materials†

The application of amorphous chalcogenide alloys as data-storage media relies on their ability to undergo an extremely fast (10-100 ns) crystallisation once heated at sufficiently high temperature. However, the peculiar features that make these materials so attractive for memory devices still lack a comprehensive microscopic understanding. By means of large scale molecular dynamics simulations, we demonstrate that the supercooled liquid of the prototypical compound GeTe shows a very high atomic mobility (D \sim 10^(-6) cm2/s) down to temperatures close to the glass transition temperatures. This behaviour leads to a breakdown of the Stokes-Einstein relation between the self- diffusion coefficient and the viscosity in the supercooled liquid. The results suggest that the fragility of the supercooled liquid is the key to understand the fast crystallisation process in this class of materials.

Vibrational properties of hexagonal Ge2Sb2Te5 from first principles

Phonons at the Γ point and the Raman spectrum of the hexagonal Ge(2)Sb(2)Te(5) were computed within density functional perturbation theory. The three different stackings of the Ge/Sb planes proposed in the experimental literature were considered. The theoretical Raman spectrum is similar for the three stackings with a marginally better agreement with experiments for the structure proposed by Matsunaga et al (2004 Acta Crystallogr. B 60 685) which assumes a disorder in Ge/Sb site occupation. Although the large broadening of the experimental Raman peaks prevents discriminating among the different stackings, the assignment of the Raman peaks to specific phonons is possible because the main features of the spectrum are rather insensitive to the actual distribution of atoms in the Sb/Ge sublattices. On the basis of the energetics (including configurational entropy) two stackings seem plausible candidates for GST, but only the mixed stacking by Matsunaga et al reproduces the spread of Ge/Sb-Te bond lengths measured experimentally.

Dynamical Heterogeneity in the Supercooled Liquid State of the Phase Change Material GeTe

A contending technology for nonvolatile memories of the next generation is based on a remarkable property of chalcogenide alloys known as phase change materials, namely their ability to undergo a fast and reversible transition between the amorphous and crystalline phases upon heating. The fast crystallization has been ascribed to the persistence of a high atomic mobility in the supercooled liquid phase, down to temperatures close to the glass transition. In this work we unravel the atomistic, structural origin of this feature in the supercooled liquid state of GeTe, a prototypical phase change compound, by means of molecular dynamic simulations. To this end, we employed an interatomic potential based on a neural network framework, which allows simulating thousands of atoms for tens of ns by keeping an accuracy close to that of the underlying first-principles framework. Our findings demonstrate that the high atomic mobility is related to the presence of clusters of slow and fast moving atoms. The latter contain a large fraction of chains of homopolar Ge-Ge bonds, which at low temperatures have a tendency to move by discontinuous cage-jump rearrangements. This structural fingerprint of dynamical heterogeneity provides an explanation of the breakdown of the Stokes-Einstein relation in GeTe, which is the ultimate origin of the fast crystallization of phase change materials exploited in the devices.

Microscopic Mechanism and Kinetics of Ice Formation at Complex Interfaces: Zooming in on Kaolinite

Most ice in nature forms because of impurities which boost the exceedingly low nucleation rate of pure supercooled water. However, the microscopic details of ice nucleation on these substances remain largely unknown. Here, we have unraveled the molecular mechanism and the kinetics of ice formation on kaolinite, a clay mineral playing a key role in climate science. We find that the formation of ice at strong supercooling in the presence of this clay is about 20 orders of magnitude faster than homogeneous freezing. The critical nucleus is substantially smaller than that found for homogeneous nucleation and, in contrast to the predictions of classical nucleation theory (CNT), it has a strong two-dimensional character. Nonetheless, we show that CNT describes correctly the formation of ice at this complex interface. Kaolinite also promotes the exclusive nucleation of hexagonal ice, as opposed to homogeneous freezing where a mixture of cubic and hexagonal polytypes is observed.

Electron-phonon interaction and thermal boundary resistance at the crystal-amorphous interface of the phase change compound GeTe

Phonon dispersion relations and electron-phonon coupling of hole-doped trigonal GeTe have been computed by density functional perturbation theory. This compound is a prototypical phase change material of interest for applications in phase change non-volatile memories. The calculations allowed us to estimate the electron-phonon contribution to the thermal boundary resistance at the interface between the crystalline and amorphous phases present in the device. The lattice contribution to the thermal boundary resistance has been computed by non-equilibrium molecular dynamics simulations with an interatomic potential based on a neural network scheme. We find that the electron-phonon term contributes to the thermal boundary resistance to an extent which is strongly dependent on the concentration and mobility of the holes. Further, for measured values of the holes concentration and electrical conductivity, the electron-phonon term is larger than the contribution from the lattice. It is also shown that the presence of Ge vacancies, responsible for the p-type degenerate character of the semiconductor, strongly affects the lattice thermal conductivity of the crystal.