skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Friday, May 17 until 8:00 AM ET on Saturday, May 18 due to maintenance. We apologize for the inconvenience.

Title: Atomistic mechanisms of phase nucleation and propagation in a model two-dimensional system
We present a computational study on the solid–solid phase transition of a model two-dimensional system between hexagonal and square phases under pressure. The atomistic mechanism of phase nucleation and propagation are determined using solid-state Dimer and nudged elastic band (NEB) methods. The Dimer is applied to identify the saddle configurations and NEB is applied to generate the transition minimum energy path (MEP) using the outputs of Dimer. Both the atomic and cell degrees of freedom are used in saddle search, allowing us to capture the critical nuclei with relatively small supercells. It is found that the phase nucleation in the model material is triggered by the localized shear deformation that comes from the relative shift between two adjacent atomic layers. In addition to the conventional layer-by-layer phase propagation, an interesting defect-assisted low barrier propagation path is identified in the hexagonal to square phase transition. The study demonstrates the significance of using the Dimer method in exploring unknown transition paths without a priori assumption. The results of this study also shed light on phase transition mechanisms of other solid-state and colloidal systems.  more » « less
Award ID(s):
2308163 2305529 2046218 1930783
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Saddle point search schemes are widely used to identify the transition state of different processes, like chemical reactions, surface and bulk diffusion, surface adsorption, and many more. In solid-state materials with relatively large numbers of atoms, the minimum mode following schemes such as dimer are commonly used because they alleviate the calculation of the Hessian on the high-dimensional potential energy surface. Here, we show that the dimer search can be further accelerated by leveraging Gaussian process regression (GPR). The GPR serves as a surrogate model to feed the dimer with the required energy and force input. We test the GPR-accelerated dimer method for predicting the diffusion coefficient of vacancy-mediated self-diffusion in body-centered cubic molybdenum and sulfur diffusion in hexagonal molybdenum disulfide. We use a multitask learning approach that utilizes a shared covariance function between energy and force input, and we show that the multitask learning significantly improves the performance of the GPR surrogate model compared to previously used learning approaches. Additionally, we demonstrate that a translation-hop sampling approach is necessary to avoid overfitting the GPR surrogate model to the minimum-mode-following pathway and thus succeeding in locating the saddle point. We show that our method reduces the number of evaluations compared to a conventional dimer method. 
    more » « less
  2. Dimerized valence bond solids appear naturally in spin-1/2 systems on bipartite lattices, with the geometric frustrations playing a key role both in their stability and the eventual `melting' due to quantum fluctuations. Here, we ask the question of the stability of such dimerized solids in spin-1 systems, taking the anisotropic square lattice with bilinear and biquadratic spin-spin interactions as a paradigmatic model. The lattice can be viewed as a set of coupled spin-1 chains, which in the limit of vanishing inter-chain coupling are known to possess a stable dimer phase. We study this model using the density matrix renormalization group (DMRG) and infinite projected entangled-pair states (iPEPS) techniques, supplemented by the analytical mean-field and linear flavor wave theory calculations. While the latter predicts the dimer phase to remain stable up to a reasonably large interchain-to-intrachain coupling ratio r≲0.6, the DMRG and iPEPS find that the dimer solid melts for much weaker interchain coupling not exceeding r≲0.15. We find the transition into a magnetically ordered state to be first order, manifested by a hysteresis and order parameter jump, precluding the deconfined quantum critical scenario. The apparent lack of stability of dimerized phases in 2D spin-1 systems is indicative of strong quantum fluctuations that melt the dimer solid. 
    more » « less
  3. Abstract

    Reliable and controllable growth of two-dimensional (2D) hexagonal boron nitride (h-BN) is essential for its wide range of applications. Substrate engineering is one of the critical factors that influence the growth of the epitaxial h-BN films. Here, we report the growth of monolayer h-BN on Ni (111) substrates incorporated with oxygen atoms via molecular beam epitaxy. It was found that the increase of incorporated oxygen concentration in the Ni substrate through a pretreatment process prior to the h-BN growth step would have an adverse effect on the morphology and growth rate of 2D h-BN. Under the same growth condition, h-BN monolayer coverage decreases exponentially as the amount of oxygen incorporated into Ni (111) increases. Density functional theory calculations and climbing image nudged elastic band (CI-NEB) method reveal that the substitutional oxygen atoms can increase the diffusion energy barrier of B and N atoms on Ni (111) thereby inhibiting the growth of h-BN films. As-grown large-area h-BN monolayer films and fabricated Al/h-BN/Ni (MIM) nanodevices were comprehensively characterized to evaluate the structural, optical and electrical properties of high-quality monolayers. Direct tunneling mechanism and high breakdown strength of ∼11.2 MV cm−1are demonstrated for the h-BN monolayers grown on oxygen-incorporated Ni (111) substrates, indicating that these films have high quality. This study provides a unique example that heterogeneous catalysis principles can be applied to the epitaxy of 2D crystals in solid state field. Similar strategies can be used to grow other 2D crystalline materials, and are expected to facilitate the development of next generation devices based on 2D crystals.

    more » « less
  4. null (Ed.)
    In systems with many local degrees of freedom, high-symmetry points in the phase diagram can provide an important starting point for the investigation of their properties throughout the phase diagram. In systems with both spin and orbital (or valley) degrees of freedom such a starting point gives rise to SU(4)-symmetric models.Here we consider SU(4)-symmetric "spin'' models, corresponding to Mott phases at half-filling, i.e. the six-dimensional representation of SU(4). This may be relevant to twisted multilayer graphene.In particular, we study the SU(4) antiferromagnetic "Heisenberg'' model on the triangular lattice, both in the classical limit and in the quantum regime. Carrying out a numerical study using the density matrix renormalization group (DMRG), we argue that the ground state is non-magnetic.We then derive a dimer expansion of the SU(4) spin model. An exact diagonalization (ED) study of the effective dimer model suggests that the ground state breaks translation invariance, forming a valence bond solid (VBS) with a 12-site unit cell.Finally, we consider the effect of SU(4)-symmetry breaking interactions due to Hund's coupling, and argue for a possible phase transition between a VBS and a magnetically ordered state. 
    more » « less
  5. TiO 2 has been extensively studied in many fields including photocatalysis, electrochemistry, optics, etc. Understanding the mechanism of the anatase–rutile phase transition (ART) process is critical for the design of TiO 2 -based high-activity photocatalysts and tuning its properties for other applications. In this work, the ART process using individual anatase micro-particles with a large percentage of (001) facets was monitored and studied. Phase concentration evolution obtained via Raman microscopy was correlated with the morphological evolution observed in scanning electron microscope (SEM) images. The ART of anatase microcrystals is dominated by surface nucleation and growth, but the ART processes of individual anatase particles are distinctive and depend on the various rutile nucleation sites. Two types of transformation pathways are observed. In one type of ART pathway, the rutile phase nucleated at a corner of an anatase microcrystal and grew in one direction along the edge of the crystal firstly followed by propagation over the rest of the microcrystal in the orthogonal direction on the surface and to the bulk of the crystal. The kinetics of the ART follows the first-order model with two distinct rate constants. The fast reaction rate is from the surface nucleation and growth, and the slow rate is from the bulk nucleation and growth. In the other type of ART pathway, multiple rutile nucleation sites formed simultaneously on different edges and corners of the microcrystal. The rutile phase spread over the whole crystal from these nucleation sites with a small contribution of bulk nucleation. Our study on the ART of individual micro-sized crystals bridges the material gap between bulk crystals and nano-sized TiO 2 particles. The anatase/rutile co-existing particle will provide a perfect platform to study the synergistic effect between the anatase phase and the rutile phase in their catalytic performances. 
    more » « less