Grain Boundary Structures and Collective Dynamics of Inversion Domains in Binary Two-Dimensional Materials

Phase-field crystal modeling of graphene/hexagonal boron nitride interfaces

Two-dimensional (2D) materials such as graphene and hexagonal boron nitride (h-BN) are an essential class of materials with enhanced structural and electronic properties compared to their bulk counterparts. The phase-field crystal (PFC) model can reach diffusive time scales to study nucleation, growth of crystallites, and relaxation of strain-driven 2D monolayers that are much larger in comparison to molecular dynamics (MD) and quantum mechanical density functional theory (QMDFT) methods while retaining atomic resolution. The model also naturally incorporates an atomic length scale and elastic and plastic deformations. We simulate the morphological transition of the crystal growth of various equilibrium crystal shapes. In this work, we generalize the one-mode PFC model to study the graphene/h-BN heterostructure interface by using conserved dynamics to describe the dynamics of the model. The model was used to find the equilibrium shape of the crystal of the h-BN crystal embedded in a graphene monolayer.

Growth mechanisms of monolayer hexagonal boron nitride (h-BN) on metal surfaces: theoretical perspectives

Two-dimensional hexagonal boron nitride (h-BN) has appeared as a promising material in diverse areas of applications, including as an excellent substrate for graphene devices, deep-ultraviolet emitters, and tunneling barriers, thanks to its outstanding stability, flat surface, and wide-bandgap. However, for achieving such exciting applications, controllable mass synthesis of high-quality and large-scale h-BN is a precondition. The synthesis of h-BN on metal surfaces using chemical vapor deposition (CVD) has been extensively studied, aiming to obtain large-scale and high-quality materials. The atomic-scale growth process, which is a prerequisite for rationally optimizing growth circumstances, is a key topic in these investigations. Although theoretical investigations on h-BN growth mechanisms are expected to reveal numerous new insights and understandings, different growth methods have completely dissimilar mechanisms, making theoretical research extremely challenging. In this article, we have summarized the recent cutting-edge theoretical research on the growth mechanisms of h-BN on different metal substrates. On the frequently utilized Cu substrate, h-BN development was shown to be more challenging than a simple adsorption-dehydrogenation-growth scenario. Controlling the number of surface layers is also an important challenge. Growth on the Ni surface is controlled by precipitation. An unusual reaction-limited aggregation growth behavior has been seen on interfaces having a significant lattice mismatch to h-BN. With intensive theoretical investigations employing advanced simulation approaches, further progress in understanding h-BN growth processes is predicted, paving the way for guided growth protocol design.

Control of phase ordering and elastic properties in phase field crystals through three-point direct correlation

Effects of three-point direct correlation on properties of the phase field crystal (PFC) modeling are examined for the control of various ordered and disordered phases and their coexistence in both three-dimensional and two-dimensional systems. Such effects are manifested via the corresponding gradient nonlinearity in the PFC free-energy functional that is derived from classical density functional theory. Their significant impacts on the stability regimes of ordered phases, phase diagrams, and elastic properties of the system, as compared to those of the original PFC model, are revealed through systematic analyses and simulations. The nontrivial contribution from three-point direct correlation leads to the variation of the critical point of order-disorder transition to which all the phase boundaries in the temperature-density phase diagram converge. It also enables the variation and control of system elastic constants over a substantial range as needed in modeling different types of materials with the same crystalline structure but different elastic properties. The capability of this PFC approach in modeling both solid and soft matter systems is further demonstrated through the effect of three-point correlation on controlling the vapor-liquid-solid coexistence and transitions for body-centered cubic phase and on achieving the liquid-stripe or liquid-lamellar phase coexistence. All these provide a valuable and efficient method for the study of structural ordering and evolution in various types of material systems.

Grain boundary and misorientation angle-dependent thermal transport in single-layer MoS2

Grain boundaries (GBs) are inevitable defects in large-area MoS₂ samples but they play a key role in their properties, however, the influence of grain misorientation on thermal transport has largely remained unknown. Here, the critical role of misorientation angle in thermal transport characteristics across 5|7 polar dislocation-dominated GBs in monolayer MoS₂ is explored using nonequilibrium molecular dynamics simulations. Results show that thermal transport characteristics of defective GBs are greatly dictated by the misorientation angle, with "U"-shaped thermal conductance as misorientation angle varying from around 5.06-52.26°, as well as by GB energy, 5|7 dislocation type and the grain size. Such unique thermal transport across GBs is primarily attributed to rising phonon-boundary softening and scattering with increasing dislocation density at GBs or GB energy, as well as an increase in localized phonon modes. The study establishes the fundamental relationship between GB and the thermal properties of single-layer MoS₂ and highlights the vital role of GBs in designing efficient thermoelectric and thermal management transition metal dichalcogenides.

Grain Boundary Motion in Two-Dimensional Hexagonal Boron Nitride

An in-depth understanding and precise controlling of grain boundary (GB) motion at the atomic scale are crucial for grain growth and recrystallization in polycrystalline materials. So far, the reported studies mainly focus on the GB motion in the ideal bicrystal system, while the atomic mechanisms of GB motion in polycrystals remain poorly understood. Herein, taking two-dimensional (2D) hexagonal boron nitride (h-BN) as a model system, we experimentally investigated the atomic-scale mechanisms of the GB motion in 2D polycrystals. Since GB motion is directly related to the GB structures, this article is organized following the configurations of GBs, which can be divided into straight (including symmetric and asymmetric GBs) and curved GBs. We revealed that (I) for symmetric GBs, the shear-coupled motion alone is insufficient to drive the continuous GB motion in polycrystalline materials, and GB sliding is also needed. (II) For asymmetric GBs, GB motion follows a defaceting-faceting process, in which dislocation reactions are crucial. (III) For curved GBs, shear-coupled GB motion (during grain shrinking) leads to grain rotation, and the rotation direction highly depends on the misorientation angles. (IV) Finally, we will discuss the characteristics of binary lattice h-BN and find that partial dislocations participate in the GB motion at high misorientation angles (>38°). Our results build up the framework of the atomic-scale mechanisms of the GB motion in 2D polycrystalline materials and will be instructive for technological applications such as grain growth and GB engineering.

Revealing roles of competing local structural orderings in crystallization of polymorphic systems

Most systems have more than two stable crystalline states in the phase diagram, which is known as polymorphism. Crystallization in such a system is often under strong influence of competing orderings linked to those crystals. However, how such competition affects crystal nucleation and ordering toward the final crystalline state is largely unknown. This is primarily because the competition takes place locally and thus is masked by large positional fluctuations. We develop a unique method to correctly identify local symmetries by removing their distortions due to positional fluctuations. This allows us to experimentally access the spatiotemporal fluctuations of local symmetries at a single-particle level in crystallization of a charged colloidal system near the body-centered cubic-face-centered cubic border. Thus, we successfully reveal the crucial roles of competing ordering in the initial selection of polymorphs and the final grain boundary motion toward the most stable state from a microscopic perspective.

Exploring the structure-property relationship of three-dimensional hexagonal boron nitride aerogels with gyroid surfaces

Three-dimensional hexagonal boron nitride aerogels (hBNAGs) are novel porous materials with many promising applications such as energy storage, thermal insulation and sensing. However, the structure-property relationships of hBNAGs in complicated thermo-mechanical coupled environments are still not clear. In this study, we employed a binary phase-field crystal (PFC) model to construct the atomic structures of hBNAGs, upon which the mechanical and thermal behaviors of hBNAGs were systematically investigated using large-scale atomistic simulations. It is found that the hBNAG geometry and topological defects strongly affect the mechanical and thermal properties. For example, the Young's modulus and tensile strength follow the scaling laws of mass density with a power factor of about 1.4 and 1.2, respectively, indicating that the stretching and bending combine toward tensile deformation. In addition, cracks nucleate around the octagon defects, indicating that the tensile strength is also influenced by the topological defects. Under compression, complicated crumpled deformations and ridges in the entire region are observed and the compression strength follows the scaling law of mass density with a power factor above 2.0, which means that a large portion of the hBNAGs do not contribute to the compression load bearing. We find that hBNAGs have a very low thermal conductivity of about two orders of magnitude lower than that of a hBN sheet. Also, the thermal conductivity of hBNAGs increases with increasing mass density, which also follows a scaling law. The power of the scaling law is about 0.5, indicating that the thermal conductivity has a strong nonlinear dependence on the mass density. Our work provides a deep understanding of the structure-property relationships of hBNAGs, which is useful for the engineering applications of hBNAGs.

Electron-Driven In Situ Transmission Electron Microscopy of 2D Transition Metal Dichalcogenides and Their 2D Heterostructures

Investigations on monolayered transition metal dichalcogenides (TMDs) and TMD heterostructures have been steadily increasing over the past years due to their potential application in a wide variety of fields such as microelectronics, sensors, batteries, solar cells, and supercapacitors, among others. The present work focuses on the characterization of TMDs using transmission electron microscopy, which allows not only static atomic resolution but also investigations into the dynamic behavior of atoms within such materials. Herein, we present a body of recent research from the various techniques available in the transmission electron microscope to structurally and analytically characterize layered TMDs and briefly compare the advantages of TEM with other characterization techniques. Whereas both static and dynamic aspects are presented, special emphasis is given to studies on the electron-driven in situ dynamic aspects of these materials while under investigation in a transmission electron microscope. The collection of the presented results points to a future prospect where electron-driven nanomanipulation may be routinely used not only in the understanding of fundamental properties of TMDs but also in the electron beam engineering of nanocircuits and nanodevices.

Mechanical relaxation and fracture of phase field crystals

A computational method is developed for the study of mechanical response and fracture behavior of phase field crystals (PFC), to overcome a limitation of the PFC dynamics which lacks an effective mechanism for describing fast mechanical relaxation of the material system. The method is based on a simple interpolation scheme for PFC (IPFC) making use of a condition of the displacement field to satisfy local elastic equilibration, while preserving key characteristics of the original PFC model. We conduct a systematic study on the mechanical properties of a sample nanoribbon system with honeycomb lattice symmetry subjected to uniaxial tension, for numerical validation of the IPFC scheme and the comparison with the original PFC and modified PFC methods. Results of mechanical response, in both elasticity and fracture regimes, show the advantage and efficiency of the IPFC method across different system sizes and applied strain rates, due to its effective process of mechanical equilibration. A brittle fracture behavior is obtained in IPFC calculations, where effects of system temperature and chirality on the fracture strength and Young's modulus are also identified, with results agreeing with those found in previous atomistic simulations of graphene. The IPFC scheme developed here is generic and applicable to the mechanical studies using different types of PFC free-energy functionals designed for various material systems.

Moiré and honeycomb lattices through self-assembly of hard-core/soft-shell microgels: experiment and simulation

Moiré and honeycomb lattices result from the sequential double deposition of monolayers of core/shell microgels in dependence of the drying conditions.

Separation of Elastic and Plastic Timescales in a Phase Field Crystal Model

A consistent small-scale description of plasticity and dislocation motion in a crystalline solid is presented based on the phase field crystal description. By allowing for independent mass motion and lattice distortion, the crystal can maintain elastic equilibrium on the timescale of plastic motion. We show that the singular (incompatible) strains are determined by the phase field crystal density, while the smooth distortions are constrained to satisfy elastic equilibrium. A numerical implementation of the model is presented and used to study a benchmark problem: the motion of an edge dislocation dipole in a triangular lattice. The time dependence of the dipole separation agrees with continuum elasticity with no adjustable parameters.

Size-Dependent Grain-Boundary Structure with Improved Conductive and Mechanical Stabilities in Sub-10-nm Gold Crystals

Low-angle grain boundaries generally exist in the form of dislocation arrays, while high-angle grain boundaries (misorientation angle >15°) exist in the form of structural units in bulk metals. Here, through in situ atomic resolution aberration corrected electron microscopy observations, we report size-dependent grain-boundary structures improving both stabilities of electrical conductivity and mechanical properties in sub-10-nm-sized gold crystals. With the diameter of a nanocrystal decreasing below 10 nm, the high-angle grain boundary in the crystal exists as an array of dislocations. This size effect may be of importance to a new generation of interconnects applications.

Two-component structural phase-field crystal models for graphene symmetries

We extend the three-point XPFC model of Seymour & Provatas (Seymour & Provatas 2016 Phys. Rev. B93, 035447 (doi:10.1103/PhysRevB.93.035447)) to two components to capture chemical vapour deposition-grown graphene, and adapt a previous two-point XPFC model of Greenwood et al. (Greenwood et al. 2011 Phys. Rev. B84, 064104 (doi:10.1103/PhysRevB.84.064104)) into a simple model of two-component graphene. The equilibrium properties of these models are examined and the two models are compared and contrasted. The first model is used to study the possible roles of hydrogen in graphene grain boundaries. The second model is used to study the role of hydrogen in the dendritic growth morphologies of graphene. The latter results are compared with new experiments.This article is part of the theme issue 'From atomistic interfaces to dendritic patterns'.

Heat transport in pristine and polycrystalline single-layer hexagonal boron nitride

Unusual thermal transport in polycrystalline h-BN prepared by phase field crystal model is revealed by large-scale molecular dynamics simulations.