Direct determination of Lennard-Jones crystal surface free energy by a computational cleavage method.

The surface free energy of solids, γ, plays a crucial role in all physical and chemical processes involving material surfaces. For the first time, we obtained γ directly from molecular dynamics simulations using a crystal cleavage method. The approach was successfully realized in a Lennard-Jones system by inserting two movable external walls, each consisting of a single crystal layer, into a bulk crystal to create flat, defect-free surfaces. The cleavage technique designed allowed us to calculate the surface free energy according to its definition and avoid surface premelting. The temperature dependence of γ was determined for the (100) and (110) crystal planes along the whole sublimation line and its metastable extension, up to T = 1.02 · Tm, where Tm is the melting point. Good agreement with indirect values of γ(T) was found. The proposed computational cleavage method can be applied to other solids of interest, providing valuable insight into the understanding of chemical and physical surface processes, and demonstrates the successful import of the cleavage method, traditionally used in technical preparation and study of crystal surfaces, into a modern atomistic simulation.

Neural Network Quantum Molecular Dynamics, Intermediate Range Order in GeSe2, and Neutron Scattering Experiments.

A remarkable property of certain covalent glasses and their melts is intermediate range order, manifested as the first sharp diffraction peak (FSDP) in neutron-scattering experiments, as was exhaustively investigated by Price, Saboungi, and collaborators. Atomistic simulations thus far have relied on either quantum molecular dynamics (QMD), with systems too small to resolve FSDP, or classical molecular dynamics, without quantum-mechanical accuracy. We investigate prototypical FSDP in GeSe2 glass and melt using neural-network quantum molecular dynamics (NNQMD) based on machine learning, which allows large simulation sizes with validated quantum mechanical accuracy to make quantitative comparisons with neutron data. The system-size dependence of the FSDP height is determined by comparing QMD and NNQMD simulations with experimental data. Partial pair distribution functions, bond-angle distributions, partial and neutron structure factors, and ring-size distributions are presented. Calculated FSDP heights agree quantitatively with neutron scattering data for GeSe2 glass at 10 K and melt at 1100 K.

Crystal Nucleation Kinetics in Supercooled Germanium: MD Simulations versus Experimental Data.

The validity of the classical nucleation theory (CNT), the most important tool to describe and predict nucleation kinetics in supercooled liquids, has been at stake for almost a century. Here, we carried out comprehensive molecular dynamics simulations of the nucleation kinetics of a fast quenched supercooled germanium using the Stillinger-Weber potential at six temperatures, covering a supercooling range of T/Tm = 0.70-0.86, where Tm is the equilibrium melting temperature. We used the seeding method to determine the number of particles in the critical crystal nuclei at each supercooling, which yielded n* = 150-1300 atoms. The transport coefficient at the liquid/nucleus interface and the melting point were also obtained from the simulations. Using the parameters resulting directly from the simulations, the CNT embraces the experimental nucleation rates, J(T), with the following fitted (average) values of the nucleus/liquid interfacial free energy: γ = 0.244 and 0.201 J/m2, for the experimental and calculated values of thermodynamic driving force, Δµ(T), respectively, which are close to the value obtained from n*(T). Without using any fit parameter, the calculated nucleation rates for the experimental and calculated values of Δµ(T) embrace the experimental J(T) curve. Therefore, this finding favors the validity of the CNT.

On the use of atomistic simulations to aid bulk metallic glasses structural elucidation with solid-state NMR.

Solid-state nuclear magnetic resonance (ssNMR) experimental 27Al metallic shifts reported in the literature for bulk metallic glasses (BMGs) were revisited in the light of state-of-the-art atomistic simulations. In a consistent way, the Gauge-Including Projector Augmented-Wave (GIPAW) method was applied in conjunction with classical molecular dynamics (CMD). A series of Zr-Cu-Al alloys with low Al concentrations were selected as case study systems, for which realistic CMD derived structural models were used for a short- and medium-range order mining. That initial procedure allowed the detection of trends describing changes on the microstructure of the material upon Al alloying, which in turn were used to guide GIPAW calculations with a set of abstract systems in the context of ssNMR. With essential precision and accuracy, the ab initio simulations also yielded valuable trends from the electronic structure point of view, which enabled an overview of the bonding nature of Al-centered clusters as well as its influence on the experimental ssNMR outcomes. The approach described in this work might promote the use of ssNMR spectroscopy in research on glassy metals. Moreover, the results presented demonstrate the possibility to expand the applications of this technique, with deeper insight into nuclear interactions and less speculative assignments.

Glass forming ability and alloying effect of a noble-metal-based glass former.

This work addresses the question on how the glass-forming ability (GFA) of a binary Pd-Ni metallic glass can be enhanced by the alloying effect of Pt. The structural features and slow dynamics of liquid and glassy states on both alloys are investigated by molecular dynamics simulations. Both alloys show typical features of glassy dynamics, namely, the non-Arrhenian behavior of diffusion and relaxation and the fractional Stokes-Einstein relation validity at low temperatures. On the basis of the analysis of the dynamical susceptibilities, we demonstrate that there is a strong influence of the alloying effect on the collective motion of the species, revealing that the GFA of the binary liquid increases with Pt alloying.