No Arabic abstract
We introduce a local order metric (LOM) that measures the degree of order in the neighborhood of an atomic or molecular site in a condensed medium. The LOM maximizes the overlap between the spatial distribution of sites belonging to that neighborhood and the corresponding distribution in a suitable reference system. The LOM takes a value tending to zero for completely disordered environments and tending to one for environments that match perfectly the reference. The site averaged LOM and its standard deviation define two scalar order parameters, $S$ and $delta S$, that characterize with excellent resolution crystals, liquids, and amorphous materials. We show with molecular dynamics simulations that $S$, $delta S$ and the LOM provide very insightful information in the study of structural transformations, such as those occurring when ice spontaneously nucleates from supercooled water or when a supercooled water sample becomes amorphous upon progressive cooling.
For a long time, there have been huge discrepancies between different models and experiments concerning the liquid-liquid phase transition (LLPT) in dense hydrogen. In this work, we present the results of extensive calculations of the LLPT in dense hydrogen using the most expensive first-principle path-integral molecular dynamics simulations available. The nonlocal density functional rVV10 and hybrid functional PBE0 are used to improve the description of the electronic structure of hydrogen. Of all the density functional theory calculations available, we report the most consistent results through quantum Monte Carlo simulations and coupled electron-ion Monte Carlo simulations of the LLPT in dense hydrogen. The critical point of the first-order LLPT is estimated above 2000 K according to the equation of state. Moreover, the metallization pressure obtained from the jump of dc electrical conductivity almost coincides with the plateau of equation of state.
Electron tomography (ET) has been demonstrated to be a powerful tool in addressing challenging problems, such as understanding 3D interactions among various microstructures. Advancing ET to broader applications requires novel instrumentation design to break the bottlenecks both in theory and in practice. In this work, we built a compact four-degree-of-freedom (three-directional positionings plus self-rotation) nano-manipulator dedicated to ET applications, which is called X-Nano transmission electron microscope (TEM) holder. All the movements of the four degrees of freedom are precisely driven by built-in piezoelectric actuators, minimizing the artefacts due to the vibration and drifting of the TEM stage. Full 360o rotation is realized with an accuracy of 0.05o in the whole range, which solves the missing wedge problem. Meanwhile, the specimen can move to the rotation axis with an integrated 3D nano-manipulator, greatly reducing the effort in tracking sample locations during tilting. Meanwhile, in-situ stimulation function can be seamlessly integrated into the X-Nano TEM holder so that dynamic information can be uncovered. We expect that more delicate researches, such as those about 3D microstructural evolution, can be carried out extensively by means of this holder in the near future.
By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semi-local density functional theory (DFT), and thereby furnish a more accurate and reliable description of the electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we have devised a linear-scaling yet formally exact approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions, MLWFs) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems at the hybrid DFT level. We provide a comprehensive description of the exx algorithm, which is currently implemented in the Quantum ESPRESSO program and employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize high-performance computing (HPC) resources. This is followed by a critical assessment of the accuracy and parallel performance of this approach when performing AIMD simulations of liquid water in the canonical ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT based AIMD simulations of condensed-phase systems containing 500-1000 atoms with a walltime cost that is comparable to semi-local DFT. In doing so, exx takes us closer to routinely performing AIMD simulations of large-scale condensed-phase systems for sufficiently long timescales at the hybrid DFT level of theory.
FLAME is a software package to perform a wide range of atomistic simulations for exploring the potential energy surfaces (PES) of complex condensed matter systems. The range of methods include molecular dynamics simulations to sample free energy landscapes, saddle point searches to identify transition states, and gradient relaxations to find dynamically stable geometries. In addition to such common tasks, FLAME implements a structure prediction algorithm based on the minima hopping method (MHM) to identify the ground state structure of any system given solely the chemical composition, and a framework to train a neural network potential to reproduce the PES from $textit{ab initio}$ calculations. The combination of neural network potentials with the MHM in FLAME allows a highly efficient and reliable identification of the ground state as well as metastable structures of molecules and crystals, as well as of nano structures, including surfaces, interfaces, and two-dimensional materials. In this manuscript, we provide detailed descriptions of the methods implemented in the FLAME code and its capabilities, together with several illustrative examples.
We introduce a local machine-learning method for predicting the electron densities of periodic systems. The framework is based on a numerical, atom-centred auxiliary basis, which enables an accurate expansion of the all-electron density in a form suitable for learning isolated and periodic systems alike. We show that using this formulation the electron densities of metals, semiconductors and molecular crystals can all be accurately predicted using a symmetry-adapted Gaussian process regression model, properly adjusted for the non-orthogonal nature of the basis. These predicted densities enable the efficient calculation of electronic properties which present errors on the order of tens of meV/atom when compared to ab initio density-functional calculations. We demonstrate the key power of this approach by using a model trained on ice unit cells containing only 4 water molecules to predict the electron densities of cells containing up to 512 molecules, and see no increase in the magnitude of the errors of derived electronic properties when increasing the system size. Indeed, we find that these extrapolated derived energies are more accurate than those predicted using a direct machine-learning model.