Qi-Jun Hong

Solid-liquid coexistence in small systems: A statistical method to calculate melting temperatures

We propose an efficient and accurate scheme to calculate the melting point (MP) of materials. This method is based on the statistical analysis of small-size coexistence molecular dynamics simulations. It eliminates the risk of metastable superheated solid in the fast-heating method, while also significantly reducing the computer cost relative to the traditional large-scale coexistence method. Using empirical potentials, we validate the method and systematically study the finite-size effect on the calculated MPs. The method converges to the exact result in the limit of large system size. An accuracy within 100 K in MP is usually achieved when simulation contains more than 100 atoms. Density functional theory examples of tantalum, high-pressure sodium, and ionic material NaCl are shown to demonstrate the accuracy and flexibility of the method in its practical applications. The method serves as a promising approach for large-scale automated material screening in which the MP is a design criterion.

Direct first-principles chemical potential calculations of liquids.

We propose a scheme that drastically improves the efficiency of Widoms particle insertion method by efficiently sampling cavities while calculating the integrals providing the chemical potentials of a physical system. This idea enables us to calculate chemical potentials of liquids directly from first-principles without the help of any reference system, which is necessary in the commonly used thermodynamic integration method. As an example, we apply our scheme, combined with the density functional formalism, to the calculation of the chemical potential of liquid copper. The calculated chemical potential is further used to locate the melting temperature. The calculated results closely agree with experiments.

Epicycle method for elasticity limit calculations

The task of finding the smallest energy needed to bring a solid to its onset of mechanical instability arises in many problems in materials science, from the determination of the elasticity limit to the consistent assignment of free energies to mechanically unstable phases. However, unless the space of possible deformations is low-dimensional and a priori known, this problem is numerically difficult, as it involves minimizing a function under a constraint on its Hessian, which is computionally prohibitive to obtain in low symmetry systems, especially if electronic structure calculations are used. We propose a method that is inspired by the well-known dimer method for saddle point searches but that adds the necessary ingredients to solve for the lowest onset of mechanical instability. The method consists of two nested optimization problems. The inner one involves a dimer-like construction to find the direction of smallest curvature as well as the gradient of this curvature function. The outer optimization then minimizes energy using the result of the inner optimization problem to constrain the search to the hypersurface enclosing all points of zero minimum curvature. Example applications to both model systems and electronic structure calculations are given.

Tetrahedron-tiling method for crystal structure prediction

Reliable and robust methods of predicting the crystal structure of a compound, based only on its chemical composition, is crucial to the study of materials and their applications. Despite considerable ongoing research efforts, crystal structure prediction remains a challenging problem that demands large computational resources. Here we propose an efficient approach for first-principles crystal structure prediction. The new method explores and finds crystal structures by tiling together elementary tetrahedra that are energetically favorable and geometrically matching each other. This approach has three distinguishing features: a favorable building unit, an efficient calculation of local energy, and a stochastic Monte Carlo simulation of crystal growth. By applying the method to the crystal structure prediction of various materials, we demonstrate its validity and potential as a promising alternative to current methods.

Combined computational and experimental investigation of high temperature thermodynamics and structure of cubic ZrO2 and HfO2

Structure and thermodynamics of pure cubic ZrO2 and HfO2 were studied computationally and experimentally from their tetragonal to cubic transition temperatures (2311 and 2530 °C) to their melting points (2710 and 2800 °C). Computations were performed using automated ab initio molecular dynamics techniques. High temperature synchrotron X-ray diffraction on laser heated aerodynamically levitated samples provided experimental data on volume change during tetragonal-to-cubic phase transformation (0.55 ± 0.09% for ZrO2 and 0.87 ± 0.08% for HfO2), density and thermal expansion. Fusion enthalpies were measured using drop and catch calorimetry on laser heated levitated samples as 55 ± 7 kJ/mol for ZrO2 and 61 ± 10 kJ/mol for HfO2, compared with 54 ± 2 and 52 ± 2 kJ/mol from computation. Volumetric thermal expansion for cubic ZrO2 and HfO2 are similar and reach (4 ± 1)·10−5/K from experiment and (5 ± 1)·10−5/K from computation. An agreement with experiment renders confidence in values obtained exclusively from computation: namely heat capacity of cubic HfO2 and ZrO2, volume change on melting, and thermal expansion of the liquid to 3127 °C. Computed oxygen diffusion coefficients indicate that above 2400 °C pure ZrO2 is an excellent oxygen conductor, perhaps even better than YSZ.