No Arabic abstract
We extend the nested sampling algorithm to simulate materials under periodic boundary and constant pressure conditions, and show how it can be used to determine the complete equilibrium phase diagram, for a given potential energy function, efficiently and in a highly automated fashion. The only inputs required are the composition and the desired pressure and temperature ranges, in particular, solid-solid phase transitions are recovered without any a priori knowledge about the structure of solid phases. We benchmark and showcase the algorithm on the periodic Lennard-Jones system, aluminium and NiTi.
We introduce a massively parallel replica-exchange grand-canonical sampling algorithm to simulate materials at realistic conditions, in particular surfaces and clusters in reactive atmospheres. Its purpose is to determine in an automated fashion equilibrium phase diagrams for a given potential-energy surface (PES) and for any observable sampled in the grand-canonical ensemble. The approach enables an unbiased sampling of the phase space and is embarrassingly parallel. It is demonstrated for a model of Lennard-Jones system describing a surface in contact with a gas phase. Furthermore, the algorithm is applied to Si$_M$ clusters ($M=2, 4$) in contact with an H$_{2}$ atmosphere, with all interactions described at the textit{ab initio} level, i.e., via density-functional theory, with the PBE gradient-corrected exchange-correlation functional. We identify the most thermodynamically stable phases at finite $T, p$(H$_{2}$) conditions.
We develop a method to efficiently construct phase diagrams using machine learning. Uncertainty sampling (US) in active learning is utilized to intensively sample around phase boundaries. Here, we demonstrate constructions of three known experimental phase diagrams by the US approach. Compared with random sampling, the US approach decreases the number of sampling points to about 20%. In particular, the reduction rate is pronounced in more complicated phase diagrams. Furthermore, we show that using the US approach, undetected new phase can be rapidly found, and smaller number of initial sampling points are sufficient. Thus, we conclude that the US approach is useful to construct complicated phase diagrams from scratch and will be an essential tool in materials science.
Supersolid is a mysterious and puzzling state of matter whose possible existence has stirred a vigorous debate among physicists for over 60 years. Its elusive nature stems from the coexistence of two seemingly contradicting properties, long-range order and superfluidity. We report computational evidence of a supersolid phase of deuterium under high pressure ($p >800$ GPa) and low temperature (T $<$ 1.0 K). In our simulations, that are based on bosonic path integral molecular dynamics, we observe a highly concerted exchange of atoms while the system preserves its crystalline order. The exchange processes are favoured by the soft core interactions between deuterium atoms that form a densely packed metallic solid. At the zero temperature limit, Bose-Einstein condensation is observed as the permutation probability of $N$ deuterium atoms approaches $1/N$ with a finite superfluid fraction. Our study provides concrete evidence for the existence of a supersolid phase in high-pressure deuterium and could provide insights on the future investigation of supersolid phases in real materials.
The complexity of strongly correlated electron physics in vanadium dioxide is exemplified as its rich phase diagrams of all kinds, which in turn shed light on the mechanisms behind its various phase transitions. In this work, we map out the hydrostatic pressure - temperature phase diagram of vanadium dioxide nanobeams by independently varying pressure and temperature with a diamond anvil cell. In addition to the well-known insulating M1 (monoclinic) and metallic R (tetragonal) phases, the diagram identifies the existence at high pressures of the insulating M1 (monoclinic, more conductive than M1) phase, and two metallic phases of X (monoclinic) and O (orthorhombic, at high temperature only). Systematic optical and electrical measurements combined with density functional calculations allow us to delineate their phase boundaries as well as reveal some basic features of the transitions.
The phase diagram of Zn has been explored up to 140 GPa and 6000 K, by combining optical observations, x-ray diffraction, and ab-initio calculations. In the pressure range covered by this study, Zn is found to retain a hexagonal close-packed crystal symmetry up to the melting temperature. The known decrease of the axial ratio of the hcp phase of Zn under compression is observed in x-ray diffraction experiments from 300 K up to the melting temperature. The pressure at which the axial ratio reaches the square root of 3 value, around 10 GPa, is slightly affected by temperature. When this axial ratio is reached, we observed that single crystals of Zn, formed at high temperature, break into multiple polycrystals. In addition, a noticeable change in the pressure dependence of the axial ratio takes place at the same pressure. Both phenomena could be caused by an isomorphic second-order phase transition induced by pressure in Zn. The reported melt curve extends previous results from 24 to 135 GPa. The pressure dependence obtained for the melting temperature is accurately described up to 135 GPa by using a Simon-Glatzel equation. The determined melt curve agrees with previous low-pressure studies and with shock-wave experiments, with a melting temperature of 5060 K at 135 GPa. Finally, a thermal equation of state is reported, which at room-temperature agrees with the literature.