Xiao-Jie Zhang

Energy Landscape of Zirconia Phase Transitions

The solid-phase transitions of zirconia are important phenomena for many industrial applications. Because of the lack of tools for resolving the atom displacement pattern, the transition kinetics has been disputed for over 60 years. Here, first-principles-based stochastic surface walking (SSW) pathway sampling is utilized for resolving the mechanism of ZrO2 tetragonal-to-monoclinic solid-phase transition. Two types of lattice and atom correspondence allowed in phase transition are determined for the first time from energy criterion, which are originated from two nearly energy-degenerate lowest-transition pathways and one stress-induced ferroelastic transition channel of tetragonal phase. An orthorhombic crystal phase (Pbc2/1) is discovered to be a trapping state at low temperatures in phase transition, the presence of which does not create new orientation relation but deters transformation toughening significantly. This new finding may facilitate the design of new functional oxide materials in ceramic industry.

Double-Ended Surface Walking Method for Pathway Building and Transition State Location of Complex Reactions

Toward the activity prediction with large-scale computations, here a double-ended surface walking (DESW) method is developed for connecting two minima on a potential energy surface (PES) and locating the associated transition state (TS) using only the first derivatives. The method operates two images starting from the initial and the final states, respectively, to walk in a stepwise manner toward each other. The surface walking involves repeated bias potential addition and local relaxation with the constrained Broyden dimer method to correct the walking direction. We apply the method to a model PES, a large set of gas phase Baker reactions, and complex surface catalytic reactions, which demonstrates that the DESW method can establish a low energy pathway linking two minima even without iterative optimization of the pathway, from which the TS can be located readily. By comparing the efficiency of the new method with the existing methods, we show that the DESW method is much less computationally demanding and is applicable for reactions with complex PESs. We hope that the DESW method may be integrated with the PES sampling methods for automated reaction prediction.

From Atoms to Fullerene: Stochastic Surface Walking Solution for Automated Structure Prediction of Complex Material.

It is of general concern whether the automated structure prediction of unknown material without recourse to any knowledge from experiment is ever possible considering the daunting complexity of potential energy surface (PES) of material. Here we demonstrate that the stochastic surface walking (SSW) method can be a general and promising solution to this ultimate goal, which is applied to assemble carbon fullerenes containing up to 100 atoms (including 60, 70, 76, 78, 80, 84, 90, 96, and 100 atoms) from randomly distributed atoms, a long-standing challenge in global optimization. Combining the SSW method with a parallel replica exchange algorithm, we can locate the global minima (GM) of these large fullerenes efficiently without being trapped in numerous energy-nearly degenerate isomers. Detailed analyses on the SSW trajectories allow us to rationalize how and why the SSW method is able to explore the highly complex PES, which highlights the abilities of SSW method for surmounting the high barrier and the preference of SSW trajectories to the low energy pathways. The work demonstrates that the parallel SSW method is a practical tool for predicting unknown materials.

Graphite to Diamond: Origin for Kinetics Selectivity

Under mild static compression (15 GPa), graphite preferentially turns into hexagonal diamond, not cubic diamond, the selectivity of which is against thermodynamics. Here we, via novel potential energy surface global exploration, report seven types low energy intermediate structures at the atomic level that are key to the kinetics of graphite to diamond solid phase transition. On the basis of quantitative kinetics data, we show that hexagonal diamond has a facile initial nucleation mechanism inside graphite matrix and faster propagation kinetics owing to the presence of three coherent graphite/hexagonal diamond interfaces, forming coherent nuclei in graphite matrix. By contrast, for the lack of coherent nucleus core, the growth of cubic diamond is at least 40 times slower and its growth is inevitably mixing with that of hexagonal diamond.

Stochastic surface walking reaction sampling for resolving heterogeneous catalytic reaction network: A revisit to the mechanism of water-gas shift reaction on Cu

Heterogeneous catalytic reactions on surface and interfaces are renowned for ample intermediate adsorbates and complex reaction networks. The common practice to reveal the reaction mechanism is via theoretical computation, which locates all likely transition states based on the pre-guessed reaction mechanism. Here we develop a new theoretical method, namely, stochastic surface walking (SSW)-Cat method, to resolve the lowest energy reaction pathway of heterogeneous catalytic reactions, which combines our recently developed SSW global structure optimization and SSW reaction sampling. The SSW-Cat is automated and massively parallel, taking a rough reaction pattern as input to guide reaction search. We present the detailed algorithm, discuss the key features, and demonstrate the efficiency in a model catalytic reaction, water-gas shift reaction on Cu(111) (CO + H2O → CO2 + H2). The SSW-Cat simulation shows that water dissociation is the rate-determining step and formic acid (HCOOH) is the kinetically favorable product, instead of the observed final products, CO2 and H2. It implies that CO2 and H2 are secondary products from further decomposition of HCOOH at high temperatures. Being a general purpose tool for reaction prediction, the SSW-Cat may be utilized for rational catalyst design via large-scale computations.