Search

Proximity ferroelectricity has recently been reported as a new design paradigm for inducing ferroelectricity, where a non-ferroelectric polar material becomes a ferroelectric by interfacing with a thin ferroelectric layer. Strongly polar materials, such as AlN and ZnO, which were previously unswitchable with an external field below their dielectric breakdown fields, can now be switched with practical coercive fields when they are in intimate proximity to a switchable ferroelectric. Here, we develop a general Landau-Ginzburg theory of proximity ferroelectricity in multilayers of non-ferroelectrics and ferroelectrics to analyze their switchability and coercive fields. The theory predicts regimes of both "proximity switching" where the multilayers collectively switch, as well as "proximity suppression" where they collectively do not switch. The mechanism of the proximity ferroelectricity is an internal electric field determined by the polarization of the layers and their relative thickness in a self-consistent manner that renormalizes the double-well ferroelectric potential to lower the steepness of the switching barrier. Further reduction in the coercive field emerges from charged defects in the bulk that act as nucleation centers. The application of the theory to proximity ferroelectricity in Alx-1ScxN/AlN and Zn1-xMgxO/ZnO bilayers is demonstrated. The theory further predicts that multilayers of dielectric/ferroelectric and paraelectric/ferroelectric layers can potentially result in induced ferroelectricity in the dielectric or paraelectric layers, resulting in the entire stack being switched, an exciting avenue for new discoveries. This thawing of "frozen ferroelectrics", paraelectrics and potentially dielectrics, promises a large class of new ferroelectrics with exciting prospects for previously unrealizable domain-patterned optoelectronic and memory technologies., Comment: 42 pages including 7 figures and 4 Appendices. To be submitted to Physical Review

Ferroelectric materials underlie key optical technologies in optical communications, integrated optics and quantum computing. Yet, there is a lack of a consistent thermodynamic framework to predict the optical properties of ferroelectrics and the mutual connections among ferroelectric polarization, optical properties, and optical dispersion. For example, there is no existing thermodynamic model for establishing the relationship between the ferroelectric polarization and the optical properties in the visible spectrum. Here we present a thermodynamic theory of the linear optical and electro-optic properties of ferroelectrics by separating the lattice and electronic contributions to the total polarization. We introduce a biquadratic coupling between the lattice and electronic contributions validated by both first-principles calculations and experimental measurements. As an example, we derive the temperature and wavelength-dependent anisotropic optical properties of BaTiO3, including the full linear optical dielectric tensor and the linear electro-optic (Pockels) effect through multiple ferroelectric phase transitions, which are in excellent agreement with existing experimental data and first principles calculations. This general framework incorporates essentially all optical properties of materials, including coupling between the ionic and electronic order parameters, as well as their dispersion and temperature dependence, and thus offers a powerful theoretical tool for analyzing light-matter interactions in ferroelectrics-based optical devices., Comment: 29 pages, 8 figures

The single variant spin cycloid and associated antiferromagnetic order in multiferroic BiFeO3 can provide a direct and predictable magnetoelectric coupling to ferroelectric order for deterministic switching, and also a key to fundamental understanding of spin transport and magnon-based applications in the system. (111) oriented BiFeO3 supplies an easy magnetic plane for the spin cycloid, but despite previous efforts, achieving deterministic switching of the single spin cycloid over multiple cycles remains challenging due to the presence of multiple spin cycloid domains in the (111) plane. Here we show that anisotropic in-plane strain engineering can stabilize a single antiferromagnetic domain and provide robust, deterministic switching. We grow BiFeO3 on orthorhombic NdGaO3 (011)o [(111)pc] substrates, breaking the spin cycloid degeneracy and by imposing a uniaxial strain in the (111) plane. This stabilization is confirmed through direct imaging with scanning NV microscopy and non-resonant X-ray diffraction. Remarkably, we achieved deterministic and non-volatile 180{\deg} switching of ferroelectric and associated antiferromagnetic domains over 1,000 cycles, significantly outperforming existing approaches. Our findings underscore that anisotropic strain engineering opens up exciting possibilities for (111)pc monodomain BiFeO3 in potential magnetoelectric and emerging magnonic applications., Comment: 17 pages, 4 figures

In CALPHAD-type thermodynamic databases, nonstoichiometric compounds are typically described by sublattice models where the sublattice site fractions represent the occupation probability of different atomic, ionic or defect species on different sublattices. Here, we develop a general procedure and corresponding linear algebra tools for converting the sublattice site fractions to a combination of independent component compositions and internal process order parameters describing the extent of internal atomic exchange, electronic redox and defect generation reactions. We apply them to a number of nonstoichiometric phases in thermodynamic databases and literature. The general procedure can be applied to constructing thermodynamic databases in terms of internal process order parameters for nonstoichiometric phases in multicomponent systems such as high-entropy oxides and alloys, which can be utilized to model their kinetics of nonequilibrium processes and microstructure evolution.

The incorporation of rare-earth elements in wurtzite nitride semiconductors, e.g., scandium alloyed aluminum nitride (ScAlN), promises dramatically enhanced piezoelectric responses, critical to a broad range of acoustic, electronic, photonic, and quantum devices and applications. Experimentally, however, the measured piezoelectric responses of nitride semiconductors are far below what theory has predicted. Here, we show that the use of a simple, scalable, post-growth thermal annealing process can dramatically boost the piezoelectric response of ScAlN thin films. We achieve a remarkable 3.5-fold increase in the piezoelectric modulus, d33 for 30% Sc content ScAlN, from 12.3 pC/N in the as-grown state to 45.5 pC/N, which is eight times larger than that of AlN. The enhancement in piezoelectricity has been unambiguously confirmed by three separate measurement techniques. Such a dramatic enhancement of d33 has been shown to impact the effective electromechanical coupling coefficient kt2 : increasing it from 13.8% to 76.2%, which matches the highest reported values in millimeter thick lithium niobate films but is achieved in a 100 nm ScAlN with a 10,000 fold reduction in thickness, thus promising extreme frequency scaling opportunities for bulk acoustic wave resonators for beyond 5G applications. By utilizing a range of material characterization techniques, we have elucidated the underlying mechanisms for the dramatically enhanced piezoelectric responses, including improved structural quality at the macroscopic scale, more homogeneous and ordered distribution of domain structures at the mesoscopic scale, and the reduction of lattice parameter ratio (c/a) for the wurtzite crystal structure at the atomic scale. Overall, the findings present a simple yet highly effective pathway that can be extended to other material families to further enhance their piezo responses.

Solid-state homogeneous precipitation of nano-sized precipitates is one of the most effective processes to strengthen metal alloys, where the final density and size distribution of precipitates are largely controlled by the precipitation kinetics. Here, we report a strategy to inoculate the homogeneous precipitation of coherent precipitates to enhance the precipitation strengthening. Using the technologically important dilute Al-Zr alloys as an example, we demonstrate that an addition of a trace level of economical and readily available, non-L1$_{2}$ phase forming impurity atoms, X (X= Sn, Sb, Bi or Cd) and Si, can significantly enhance the diffusivity of Zr atoms and overturn the precipitation of L1$_{2}$-structured Al$_{3}$Zr nanoparticles from the classical homogeneous nucleation and growth pathway into a nonclassical nucleation pathway: Al$_{3}$Zr forms through the spontaneous formation of nano-scale local concentration fluctuations of Zr atoms on Zr-X(-Si)-vacancy clusters followed by a continuous increase of the concentration and chemical short-range ordering (CSRO). Such an impurity atoms induced heterogeneous nucleation based on a "spinodal decomposition like" mechanism dramatically accelerates the precipitation kinetics, leading to an order of magnitude higher number density of precipitates and a record high hardening efficiency of solute Zr atoms. By formulating the generalized selection principles for inoculating impurity elements, this inoculation strategy should be extendable to a broader range of materials to further explore the precipitation strengthening potentials.

We unravel the distinct roles each cation plays in phase evolution, stability, and properties within Mg1/5Co1/5Ni1/5Cu1/5Zn1/5O high-entropy oxide (HEO) by integrating experimental findings, thermodynamic analyses, and first-principles predictions. Our approach is through sequentially removing one cation at a time from the five-component high-entropy oxide to create five four-component derivatives. Bulk synthesis experiments indicate that Mg, Ni, and Co act as rock salt phase stabilizers whereas only Mg and Ni enthalpically enhance single-phase rock salt stability in thin film growth; synthesis conditions dictate whether Co is a rock salt phase stabilizer or destabilizer. By examining the competing phases and oxidation state preferences using pseudo-binary phase diagrams and first-principles calculations, we resolve the stability differences between bulk and thin film for all compositions. We systematically explore HEO macroscopic property sensitivity to cation selection employing both predicted and measured optical spectra. This study establishes a framework for understanding high-entropy oxide synthesizability and properties on a per-cation basis that is broadly applicable to tailoring functional property design in other high-entropy materials., Comment: Saeed S. I. Almishal and Jacob T. Sivak contributed equally to this work

An analytical theory is developed for predicting the nonlinear susceptibility of ionic polarization to continuous electromagnetic waves in both bulk and strained thin film ferroelectrics. Using a perturbation method for solving the nonlinear equation of motion for ionic polarization within the framework of Landau-Ginzburg-Devonshire theory, the full second-order nonlinear susceptibility tensor is derived as a function of frequency, temperature, and strain. The theory predicts the coexistence of a significantly enhanced second-order dielectric susceptibility and a relatively low dielectric loss in BaTiO3 films with a strain-stabilized monoclinic ferroelectric phase and in a strained SrTiO3 film near its temperature-driven second-order ferroelectric-to-paraelectric phase transition. This work establishes a theoretical framework for predicting and exploiting nonlinear interactions between THz waves and ferroelectric materials, and more generally, suggests exciting opportunities to strain-engineer nonlinear dynamical properties of ferroelectrics beyond the static and quasi-static limits.

Interest in high-entropy inorganic compounds originates from their ability to stabilize cations and anions in local environments that rarely occur at standard temperature and pressure. This leads to new crystalline phases in many-cation formulations with structures and properties that depart from conventional trends. The highest-entropy homogeneous and random solid-solution is a parent structure from which a continuum of lower-entropy offspring can originate by adopting chemical and/or structural order. This report demonstrates how synthesis conditions, thermal history, and elastic and chemical boundary conditions conspire to regulate this process in Mg0.2Co0.2Ni0.2Cu0.2Zn0.2O, during which coherent CuO nano-tweeds and spinel nano-cuboids evolve. We do so by combining structured synthesis routes, atomic-resolution microscopy and spectroscopy, density functional theory, and a phase field modeling framework that accurately predicts the emergent structure and local chemistry. This establishes a framework to appreciate, understand, and predict the macrostate spectrum available to a high-entropy system that is critical to rationalize property engineering opportunities.

We theoretically investigate the steady-state transmission of continuous terahertz (THz) wave across a freestanding ferroelectric slab. Based on the Landau-Ginzburg-Devonshire theory of ferroelectrics and the coupled equations of motion for polarization and electromagnetic (EM) waves, we derive the analytical expressions of the frequency- and thickness-dependent dielectric susceptibility and transmission coefficient at the thin slab limit in the harmonic excitation regime. When the slab thickness is much smaller than the THz wavelength in the ferroelectric, the analytical predictions agree well with the numerical simulations from a dynamical phase-field model that incorporates the coupled dynamics of strain, polarization, and EM wave in multiphase systems. At larger thicknesses, the transmission is mainly determined by the frequency-dependent attenuation of THz waves in the ferroelectric and the formation of a standing polarization/THz wave. Our results advance the understanding of the interaction between THz wave and ferroelectrics and suggest the potential of exploiting ferroelectrics to achieve low-heat-dissipation, nonvolatile voltage modulation of THz transmission for high-data-rate wireless communication.

Optical second harmonic generation (SHG) is a nonlinear optical effect widely used for nonlinear optical microscopy and laser frequency conversion. Closed-form analytical solution of the nonlinear optical responses is essential for evaluating the optical responses of new materials whose optical properties are unknown a priori. A recent open-source code, SHAARP(si), can provide such closed form solutions for crystals with arbitrary symmetries, orientations, and anisotropic properties at a single interface. However, optical components are often in the form of slabs, thin films on substrates, and multilayer heterostructures with multiple reflections of both the fundamental and up to ten different SHG waves at each interface, adding significant complexity. Many approximations have therefore been employed in the existing analytical approaches, such as slowly varying approximation, weak reflection of the nonlinear polarization, transparent medium, high crystallographic symmetry, Kleinman symmetry, easy crystal orientation along a high-symmetry direction, phase matching conditions and negligible interference among nonlinear waves, which may lead to large errors in the reported material properties. To avoid these approximations, we have developed an open-source package named Second Harmonic Analysis of Anisotropic Rotational Polarimetry in Multilayers (SHAARP(ml)). The reliability and accuracy are established by experimentally benchmarking with both the SHG polarimetry and Maker fringes predicted from the package using standard materials.

Multifunctionality as a paradigm requires materials exhibiting multiple superior properties. Integrating second-order optical nonlinearity and large bandgap with piezoelectricity could, for example, enable broadband, strain-tunable photonics. Though very different phenomena at distinct frequencies, both second-order optical nonlinearity and piezoelectricity are third-rank polar tensors present only in acentric crystal structures. However, simultaneously enhancing both phenomena is highly challenging since it involves competing effects with tradeoffs. Recently, a large switchable ferroelectric polarization of ~ 80 uC cm-2 was reported in Zn1-xMgxO films. Here, ferroelectric Zn1-xMgxO is demonstrated to be a platform that hosts simultaneously a 30% increase in the electronic bandgap, a 50% enhancement in the second harmonic generation coefficients, and a near 200% improvement in the piezoelectric coefficients over pure ZnO. These enhancements are shown to be due to a 400% increase in the electronic anharmonicity and a ~200% decrease in the ionic anharmonicity with Mg substitution. Precisely controllable periodic ferroelectric domain gratings are demonstrated down to 800 nm domain width, enabling ultraviolet quasi-phase-matched optical harmonic generation as well as domain-engineered piezoelectric devices.

Using the Landau-Ginzburg-Devonshire approach, we study light-induced phase transitions, evolution of polar state and domain morphology in photo-ferroelectric nanoparticles (NPs). Light exposure increases the free carrier density near the NP surface and may in turn induce phase transitions from the nonpolar paraelectric to the polar ferroelectric phase. Using the uniaxial photo-ferroelectric Sn2P2S6 as an example, we show that visible light exposure induces the appearance and vanishing of striped, labyrinthine or curled domains and changes in the polarization switching hysteresis loop shape from paraelectric curves to double, pinched and single loops, as well as the shifting in the position of the tricritical point. Furthermore, we demonstrate that an ensemble of non-interacting photo-ferroelectric NPs may exhibit superparaelectric-like features at the tricritical point, such as strongly frequency-dependent giant piezoelectric and dielectric responses, which can potentially be exploited for piezoelectric applications., Comment: 42 pages, 7 figures, including 14 pages Supplement with 6 figures

Modelling elastocaloric effect (eCE) is crucial for the design of environmentally friendly and energy-efficient eCE based solid-state cooling devices. Here, a thermodynamically consistent non-isothermal phase-field model (PFM) coupling martensitic transformation with mechanics and heat transfer is developed and applied for simulating eCE. The model is derived from a thermodynamic framework which invokes the microforce theory and Coleman--Noll procedure. To avoid the numerical issue related to the non-differentiable energy barrier function across the transition point, the austenite-martensite transition energy barrier in PFM is constructed as a smooth function of temperature. Both the indirect method using isothermal PFM with Maxwell relations and the direct method using non-isothermal PFM are applied to calculate the elastocaloric properties. The former is capable of calculating both isothermal entropy change and adiabatic temperature change ($\Delta T_{\text{ad}}$), but induces high computation cost. The latter is computationally efficient, but only yields $\Delta T_{\text{ad}}$. In a model Mn-22Cu alloy, the maximum $\Delta T_{\text{ad}}$ ($\Delta T_{\text{ad}}^{\text{max}}$) under a compressive stress of 100 MPa is calculated as 9.5 and 8.5 K in single crystal (3.5 and 3.8 K in polycrystal) from the indirect and direct method, respectively. It is found that the discrepancy of $\Delta T_{\text{ad}}^{\text{max}}$ by indirect and direct method is within 10% at stress less than 150 MPa, confirming the feasibility of both methods in evaluating eCE at low stress. The results demonstrate the developed PFM herein, combined with both indirect and direct method for eCE calculations, as a practicable toolkit for the computational design of elastocaloric devices., Comment: 26 pages, 10 figures

We propose an analytical Landau-Ginzburg theory of the charge density waves coupled with lattice and electronic long-range order parameters. Examples of long-range order include electronic wave function of superconducting Cooper pairs, structural distortions, electric polarization, and magnetization. We formulate the Landau-Ginzburg free energy density as power expansion with respect to the charge density and other long-range order parameters, as well as their spatial gradients, and biquadratic coupling terms. We introduced a biquadratic coupling between the charge density gradient and long-range order parameters, as well as nonlinear higher gradients of the long-range order parameters. The biquadratic gradient coupling is critical to the appearance of different spatially-modulated phases in charge-ordered ferroics and high-temperature superconductors. We derived the thermodynamic conditions for the stability of the spatially-modulated phases, which are the intertwined spatial waves of charge density and lattice/electronic long-range order. The analytical expressions for the energies of different phases, corresponding order parameters, charge density waves amplitudes and modulation periods, obtained in this work, can be employed to guide the comprehensive physical explanation, deconvolution and Bayesian analysis of experimental data on quantum materials ranging from charge-ordered ferroics to high-temperature superconductors., Comment: 40 pages, 5 figures, Appendix with 7 figures

It has been a long challenge to analytically construct the quantitative temperature-dependent multi-well free-energy landscape over the space of order parameters describing phase transitions and associated critical phenomena. Here we propose a simple analytical model for the free energy landscape based on a priori concept of Boltzmann thermal mixing among multiple parabolic potentials representing different energetically degenerated ground states in contrast to the popular Landau theory using the high-temperature disordered state as a reference. The model recovers both the Weiss molecular field theory and the temperature dependent behaviors of the second order Landau coefficient. It is rather remarkable that such a simple analytical expression can describe a wide variety of properties across the ferroelectric phase transition in PbTiO3, including the temperature dependences of the spontaneous polarization, the heat capacities, the permittivity tensor, and the lattice parameters. The approach also allows parametrization of the free energy function across phase transitions based directly on the 0 K thermodynamics data that can be obtained from DFT calculations. We anticipate that the approach can be equally applied to other types of critical phenomena, e.g., the superconducting phase transition, the metal-insulator transition, and the magnetic transitions., Comment: 18 pages, 5 figures. arXiv admin note: text overlap with arXiv:2204.07067

Exploiting the emerging nanoscale periodicities in epitaxial, single-crystal thin films is an exciting direction in quantum materials science: confinement and periodic distortions induce novel properties. The structural motifs of interest are ferroelastic, ferroelectric, multiferroic, and, more recently, topologically protected magnetization and polarization textures. A critical step towards heterostructure engineering is understanding their nanoscale structure, best achieved through real-space imaging. X-ray Bragg coherent diffractive imaging visualizes sub-picometer crystalline displacements with tens of nanometers spatial resolution. Yet, it is limited to objects spatially confined in all three dimensions and requires highly coherent, laser-like x-rays. Here we lift the confinement restriction by developing real-space imaging of periodic lattice distortions: we combine an iterative phase retrieval algorithm with unsupervised machine learning to invert the diffuse scattering in conventional x-ray reciprocal-space mapping into real-space images of polar and elastic textures in thin epitaxial films. We first demonstrate our imaging in PbTiO3/SrTiO3 superlattices to be consistent with published phase-field model calculations. We then visualize strain-induced ferroelastic domains emerging during the metal-insulator transition in Ca2RuO4 thin films. Instead of homogeneously transforming into a low-temperature structure (like in bulk), the strained Mott insulator splits into nanodomains with alternating lattice constants, as confirmed by cryogenic scanning transmission electron microscopy. Our study reveals the type, size, orientation, and crystal displacement field of the nano-textures. The non-destructive imaging of textures promises to improve models for their dynamics and enable advances in quantum materials and microelectronics.

This work demonstrates a first-principles-based approach to obtaining finite temperature thermal and electronic transport properties which can be employed to model and understand mesoscale structural evolution during electronic, magnetic, and structural phase transitions. A computationally tractable model was introduced to estimate the temperature dependence of the electron relaxation time. The model is applied to Ca3Ru2O7 with a focus on understanding its electrical resistivity across the electronic phase transition at 48 K. A quasiharmonic phonon approach to the lattice vibrations was employed to account for thermal expansion while the Boltzmann transport theory including spin-orbit coupling was used to calculate the electron-transport properties, including the temperature dependence of electrical conductivity., Comment: 4 figures

Optical second harmonic generation is a second-order nonlinear process that combines two photons of a given frequency into a third photon at twice the frequency. Due to the symmetry constraints, it is widely used as a sensitive probe to detect broken inversion symmetry and local polar order. Analytical modeling of the electric-dipole SHG response is essential to extract fundamental properties of materials from experiments. However, complexity builds up dramatically in the analytical model when the probed crystal is of a low bulk crystal symmetry, with a low-symmetry surface orientation, exhibits absorption and dispersion, and consists of multiple interfaces. As a result, there is a largely uneven landscape in the literature on the SHG modeling of new materials, involving numerous approximations and a wide range of (in)accuracies, leading to a rather scattered dataset of reported SHG nonlinear susceptibility. Towards streamlining the reliability and accuracy of this process, we have developed an open-source package called the Second Harmonic Analysis of Anisotropic Rotational Polarimetry (SHAARP) which derives analytical solutions and performs numerical simulations of reflection SHG from a single interface for homogeneous crystals. Five key generalizations in SHG modeling are implemented, including all crystal symmetries down to triclinic, any crystal orientation, complex dielectric tensor (refractive indices) with frequency dispersion, and general polarization states of the light. SHAARP enables accurate anisotropic modeling of SHG response for a broad range of materials systems. The method is extendible to multiple interfaces. The code is free to download from https://github.com/Rui-Zu/SHAARP

The research on topological phenomena in ferroelectric materials has revolutionized the way we understand polar order. Intriguing examples are polar skyrmions, vortex/anti-vortex structures and ferroelectric incommensurabilties, which promote emergent physical properties ranging from electric-field-controllable chirality to negative capacitance effects. Here, we study the impact of topologically protected vortices on the domain formation in improper ferroelectric ErMnO3 polycrystals, demonstrating inverted domain scaling behavior compared to classical ferroelectrics. We observe that as the grain size increases, smaller domains are formed, which we relate to the interaction of the topological vortices with local strain fields. The inversion of the domain scaling behavior has far-reaching implications, providing fundamentally new opportunities for topology-based domain engineering and the tuning of the electromechanical and dielectric performance of ferroelectrics in general.

Magnetism in topological materials creates phases exhibiting quantized transport phenomena with applications in spintronics and quantum information. The emergence of such phases relies on strong interaction between localized spins and itinerant states comprising the topological bands, and the subsequent formation of an exchange gap. However, this interaction has never been measured in any intrinsic magnetic topological material. Using a multimodal approach, this exchange interaction is measured in MnBi2Te4, the first realized intrinsic magnetic topological insulator. Interrogating nonequilibrium spin dynamics, itinerant bands are found to exhibit a strong exchange coupling to localized Mn spins. Momentum-resolved ultrafast electron scattering and magneto-optic measurements reveal that itinerant spins disorder via electron-phonon scattering at picosecond timescales. Localized Mn spins, probed by resonant X-ray scattering, disorder concurrently with itinerant spins, despite being energetically decoupled from the initial excitation. Modeling the results using atomistic simulations, the exchange coupling between localized and itinerant spins is estimated to be >100 times larger than superexchange interactions. This implies an exchange gap of >25 meV should occur in the topological surface states. By directly quantifying local-itinerant exchange coupling, this work validates the materials-by-design strategy of utilizing localized magnetic order to create and manipulate magnetic topological phases, from static to ultrafast timescales.

The coexistence of different ferroelectric phases enables the tunability of the macroscopic properties and extensive applications from piezoelectric transducers to nonvolatile memories. Here we develop a thermodynamic model to predict the stability and low-energy orientations of boundaries between different phases in ferroelectrics. Taking lead zirconate titanate and bismuth ferrite as two examples, we demonstrate that the low-energy orientations of interphase boundaries are largely determined by minimizing the electrostatic and elastic energies. Phase-field simulations are employed to analyze the competition between the interfacial energy and the electrostatic and elastic energies. Our simulation results demonstrate that the lowering of crystal symmetry could occur due to the electrical and mechanical incompatibilities between the two phases, which can be used to explain the experimentally observed low-symmetry phases near morphotropic phase boundaries. Our work provides theoretical foundations for understanding and controlling the interphase boundaries in ferroelectric materials for multifunctional applications.

In condensed-matter systems, competition between ground states at phase boundaries can lead to significant changes in material properties under external stimuli, particularly when these ground states have different crystal symmetries. A key scientific and technological challenge is to stabilize and control coexistence of symmetry-distinct phases with external stimuli. Using BiFeO3 (BFO) layers confined between layers of the dielectric TbScO3 as a model system, we stabilize the mixed-phase coexistence of centrosymmetric and non-centrosymmetric BFO phases with antipolar, insulating and polar, semiconducting behavior, respectively at room temperature. Application of in-plane electric (polar) fields can both remove and introduce centrosymmetry from the system resulting in reversible, nonvolatile interconversion between the two phases. This interconversion between the centrosymmetric insulating and non-centrosymmetric semiconducting phases coincides with simultaneous changes in the non-linear optical response of over three orders of magnitude, a change in resistivity of over five orders of magnitude, and a change in the polar order. Our work establishes a materials platform allowing for novel cross-functional devices which take advantage of changes in optical, electrical, and ferroic responses., Comment: L Caretta and Y-T Shao contributed equally. 64 pages, 4 figures, 25 supplemental figures, 1 supplemental table

Recently various topological polar structures have been discovered in oxide thin films. Despite the increasing evidence of their switchability under electrical and/or mechanical fields, the dynamic property of isolated ones, which is usually required for applications such as data storage, is still absent. Here, we show the controlled nucleation and motion of isolated three-fold vertices under an applied electric field. At the PbTiO3/SrRuO3 interface, a two-unit-cell thick SrTiO3 layer provides electrical boundary conditions for the formation of three-fold vertices. Utilizing the SrTiO3 layer and in situ electrical testing system, we find that isolated three-fold vertices can move in a controllable and reversible manner with a velocity up to ~629 nm s-1. Microstructural evolution of the nucleation and propagation of isolated three-fold vertices is further revealed by phase-field simulations. This work demonstrates the ability to electrically manipulate isolated three-fold vertices, shedding light on the dynamic property of isolated topological polar structures.

Order-disorder transitions are widely explored in various vortex structures in condensed matter physics, i.e., in the type-II superconductors and Bose-Einstein condensates. In this study, we have investigated the ordering of the polar vortex phase in the (PZT)n/(STO)n superlattice systems through phase-field simulations. An antiorder state is discovered for short periodicity superlattice on an SSO substrate, owing to the huge interfacial coupling between PZT and STO as well as the giant in-plane polarization in STO layers due to the large tensile strain. Increasing the periodicity leads to the anti-order to disorder transition, resulting from the loss of interfacial coupling and disappearance of the polarization in STO layers. On the other hand, for short periodicity superlattices, order-disorder-antiorder transition can be engineered by mediating the substrate strain, due to the delicate competition between the depoling effect, interfacial coupling, and strain effect. We envision this study to spur further interest towards the understanding of order-disorder transition in ferroelectric topological structures., Comment: 20 Pages 5 figures

In this work, we present a software package in Python for high-throughput first-principles calculations of thermodynamic properties at finite temperatures, which we refer to as DFTTK (Density Functional Theory Tool Kit). DFTTK is based on the atomate package and integrates our experiences in the last decades on the development of theoretical methods and computational software. It includes task submissions on all major operating systems and task execution on high-performance computing environments. The distribution of the DFTTK package comes with examples of calculations of phonon density of states, heat capacity, entropy, enthalpy, and free energy under the quasi-harmonic phonon scheme for the stoichiometric phases of Al, Ni, Al3Ni, AlNi, AlNi3, Al3Ni4, and Al3Ni5, and the fcc solution phases treated using the special quasirandom structures at the compositions of Al3Ni, AlNi, and AlNi3., Comment: 50 pages, 18 figures

How topological defects, unavoidable at symmetry-breaking phase transitions in a wide range of systems, evolve through consecutive phase transitions with different broken symmetries remains unexplored. Nd2SrFe2O7, a bilayer ferrite, exhibits two intriguing structural phase transitions and dense networks of the so-called type-II Z8 structural vortices at room temperature, so it is an ideal system to explore the topological defect evolution. From our extensive experimental investigation, we demonstrate that the cooling rate at the second-order transition (1290oC) plays a decisive role in determining the vortex density at room temperature, following the universal Kibble-Zurek mechanism. In addition, we discovered a transformation between topologically-distinct vortices (Z8 to Z4 vortices) at the first-order transition (550oC), which conserves the number of vortex cores. Remarkably, the Z4 vortices consist of two phases with an identical symmetry but two distinct magnitudes of an order parameter. Furthermore, when lattice distortion is enhanced by chemical doping, a new type of topological defects emerges: loop domain walls with orthorhombic distortions in the tetragonal background, resulting in unique pseudo-orthorhombic twins. Our findings open a new avenue to explore the evolution of topological defects through multiple phase transitions., Comment: 21 pages, 8 figures

Topological phases such as polar skyrmions have been a fertile playground for ferroelectric oxide superlattices, with exotic physical phenomena such as negative capacitance. Herein, using phase-field simulations, we demonstrate the local control of the skyrmion phase with electric potential applied through a top electrode. Under a relatively small electric potential, the skyrmions underneath the electrode can be erased and recovered reversibly. A topologically protected transition from the symmetric to asymmetric skyrmion bubbles is observed at the edge of the electrode. While a topological transition to a labyrinthine domain requires a high applied potential, it can switch back to the skyrmion state with a relatively small electric potential. The topological transition from +1 to 0 occurs before the full destruction of the bubble state. It is shown that the shrinking and bursting of the skyrmions leads to a large reduction in the dielectric permittivity, the magnitude of which depends on the size of the electrode., Comment: 17 Pages 4 figures

Strongly correlated materials that exhibit an insulator-metal transition are key candidates in the search for new computing platforms. Understanding the pathways and timescales underlying the electrically-driven insulator-metal transition is crucial for uncovering the fundamental limits of device operation. Using stroboscopic electron diffraction, we perform synchronized time-resolved measurements of atomic motions and electronic transport in operating vanadium dioxide switches. We discover an electrically-triggered, isostructural state that forms transiently on microsecond timescales, stabilized by local heterogeneities and interfacial interactions between the equilibrium phases. This metastable phase bears striking similarity to that formed under photoexcitation within picoseconds, suggesting a universal transformation pathway across eight orders of magnitude of timescale. Our results establish a new route for uncovering non-equilibrium and metastable phases in correlated materials, and open avenues for engineering novel dynamical behavior in nanoelectronics.

The collective dynamics of topological structures have been of great interest from both fundamental and applied perspectives. For example, the studies of dynamical properties of magnetic vortices and skyrmions not only deepened the understanding of many-body physics but also led to potential applications in data processing and storage. Topological structures constructed from electrical polarization rather than spin have recently been realized in ferroelectric superlattices, promising for ultrafast electric-field control of topological orders. However, little is known about the dynamics of such complex extended nanostructures which in turn underlies their functionalities. Using terahertz-field excitation and femtosecond x-ray diffraction measurements, we observe ultrafast collective polarization dynamics that are unique to polar vortices, with orders of magnitude higher frequencies and smaller lateral size than those of experimentally realized magnetic vortices. A previously unseen soft mode, hereafter referred to as a vortexon, emerges as transient arrays of nanoscale circular patterns of atomic displacements, which reverse their vorticity on picosecond time scales. Its frequency is significantly reduced at a critical strain, indicating a condensation of structural dynamics. First-principles-based atomistic calculations and phase-field modeling reveal the microscopic atomic arrangements and frequencies of the vortex modes. The discovery of subterahertz collective dynamics in polar vortices opens up opportunities for applications of electric-field driven data processing in topological structures with ultrahigh speed and density.

Polar skyrmions are predicted to emerge from the interplay of elastic, electrostatic and gradient energies, in contrast to the key role of the anti-symmetric Dzyalozhinskii-Moriya interaction in magnetic skyrmions. With the discovery of topologically-stable polar skyrmions, it is of both fundamental and practical interest to understand the microscopic nature and the possibility of temperature- and strain-driven phase transitions in ensembles of such polar skyrmions. Here, we explore the reversible transition from a skyrmion state (topological charge of -1) to a two-dimensional, tetratic lattice of merons (with topological charge of -1/2) upon varying the temperature and elastic boundary conditions in [(PbTiO3)16/(SrTiO3)16]8 lifted-off membranes. This topological phase transition is accompanied by a change in chirality, from zero-net chirality (in meronic phase) to net-handedness (in skyrmionic phase). To map these changes microscopically required developing new imaging methods. We show how scanning convergent beam electron diffraction provides a robust measure of the local polarization simultaneously with the strain state at sub-nm resolution, while also directly mapping the chirality of each skyrmion. Using this, we demonstrate strain as a crucial order parameter to drive isotropic-to-anisotropic structural transitions of chiral polar skyrmions to non-chiral merons, validated with X-ray reciprocal space mapping and theoretical phase-field simulations. These results revealed by our new measurement methods provide the first illustration of systematic control of rich variety of topological dipole textures by altering the mechanical boundary conditions, which may offer a promising way to control their functionalities in ferroelectric nanodevices using the local and spatial distribution of chirality and order., Comment: 39 pages, 5 figures, 14 supporting figures

A comprehensive thermochemical database is constructed based on high-throughput first-principles phonon calculations of over 3000 atomic structures in Ni, Fe, and Co alloys involving a total of 26 elements including Al, B, C, Cr, Cu, Hf, La, Mn, Mo, N, Nb, O, P, Re, Ru, S, Si, Ta, Ti, V, W, Y, and Zr, providing thermochemical data largely unavailable from existing experiments. The database can be employed to predict the equilibrium phase compositions and fractions at a given temperature and an overall chemical composition directly from first-principles by minimizing the chemical potential. It is applied to the additively manufactured nickel-based IN718 superalloy to analyze the phase evolution with temperature. In particular, we successfully predicted the formation of L1$_0$-FeNi, $\gamma'$-Ni$_3$(Fe,Al), $\alpha$-Cr, $\gamma$-Ni$_3$(Nb,Mo), $\gamma''$-Ni$_3$Nb , and $\eta$-Ni$_3$Ti at low temperatures, $\gamma'$-Ni$_3$Al, $\delta$-Ni$_3$Nb, $\gamma''$-Ni$_3$Nb, $\alpha$-Cr, and $\gamma$-Ni(Fe,Cr,Mo) at intermediate temperatures, and $\delta$-Ni$_3$Nb and $\gamma$-Ni(Fe,Cr,Mo) at high temperatures in IN718. These predictions are validated by EDS mapping of compositional distributions and corresponding identifications of phase distributions. The database is expected to be a valuable source for future thermodynamic analysis and microstructure prediction of alloys involving the 26 elements.

The Lorenz number (L) contained in the Wiedemann-Franz law represents the ratio of two kinetic parameters of electronic charge carriers: the electronic contribution to the thermal conductivity (K_el) and the electrical conductivity (sigma), , and can be expressed as LT=K_el/sigma where T is temperature. We demonstrate that the Lorenz number simply equals to the ratio of two thermodynamic quantities: the electronic heat capacity (c_el) and the electrochemical capacitance (c_N) through LT=c_el/c_N , a purely thermodynamic quantity, and thus it can be calculated solely based on the electron density of states of a material. It is shown that our thermodynamic formulation for the Lorenz number leads to: i) the well-known Sommerfeld value L=pi^2/3(k_B/e)^2 at the low temperature limit, ii) the Drude value L=3/2(k_B/e)^2 at the high temperature limit with the free electron gas model, and iii) possible higher values than the Sommerfeld limit for semiconductors. It is also demonstrated that the purely electronic contribution to the thermoelectric figure-of-merit can be directly computed using high-throughput DFT calculations without resorting to the computationally more expensive Boltzmann transport theory to the electronic thermal conductivity and electrical conductivity.

Ferroelectric topological objects (e.g. vortices, skyrmions) provide a fertile ground for exploring emerging physical properties that could potentially be utilized in future configurable nanoelectronic devices. Here, we demonstrate quasi-one-dimensional metallic high conduction channels along two types of exotic topological defects, i.e. the topological cores of (i) a quadrant vortex domain structure and (ii) a center domain (monopole-like) structure confined in high quality BiFeO3 nanoisland array, abbreviated as the vortex core and the center core. We unveil via phase-field simulations that the superfine (< 3 nm) metallic conduction channels along center cores arise from the screening charge carriers confined at the core whereas the high conductance of vortex cores results from a field-induced twisted state. These conducting channels can be repeatedly and reversibly created and deleted by manipulating the two topological states via an electric field, leading to an apparent electroresistance effect with an on/off ratio higher than 103. These results open up the possibility of utilizing these functional one-dimensional topological objects in high-density nanoelectronic devices such as ultrahigh density nonvolatile memory., Comment: manuscript (21 pages, 4 figures) plus supplementary information

Metal-ion doping can effectively regulate the metal-insulator transition temperature in $\mathrm{VO}_2$. Experiments found that the pentavalent and hexavalent ion doping dramatically reduces the transition temperature while the trivalent ion doping increases the transition temperature and induces intermediate phases. Based on the phase-field model of the metal-insulator transition in $\mathrm{VO}_2$ we developed previously, we formulate a Landau potential of the metal-ion-doped $\mathrm{VO}_2$ taking account of the effects of doping on the electron correlation and lattice structure. The effect of metal-ion doping on the lattice structure is accounted for in a phenomenological way. Using the Landau potential, we calculate the temperature-dopant-concentration phase diagrams of $\mathrm{VO}_2$ doped with various metal ions consistent with the experiments and provide explanation to the different behaviors of different metal-ion doping. The phenomenological theory can provide estimations of phase diagrams of $\mathrm{VO}_2$ doped with other metal ions., Comment: 6 pages, 2 figures

Nanoelectronic devices based on ferroelectric domain walls (DWs), such as memories, transistors, and rectifiers, have been demonstrated in recent years. Practical high-speed electronics, on the other hand, usually demand operation frequencies in the giga-Hertz (GHz) regime, where the effect of dipolar oscillation is important. In this work, an unexpected giant GHz conductivity on the order of 103 S/m is observed in certain BiFeO3 DWs, which is about 100,000 times greater than the carrier-induced dc conductivity of the same walls. Surprisingly, the nominal configuration of the DWs precludes the ac conduction under an excitation electric field perpendicular to the surface. Theoretical analysis shows that the inclined DWs are stressed asymmetrically near the film surface, whereas the vertical walls in a control sample are not. The resultant imbalanced polarization profile can then couple to the out-of-plane microwave fields and induce power dissipation, which is confirmed by the phase-field modeling. Since the contributions from mobile-carrier conduction and bound-charge oscillation to the ac conductivity are equivalent in a microwave circuit, the research on local structural dynamics may open a new avenue to implement DW nano-devices for RF applications.

Cycloidal spin orders are common in multiferroics. One of the prototypical examples is BiFeO3 (BFO) which shows a large polarization and a cycloidal antiferromagnetic order at room temperature. Here we employ Landau theory and phase-field simulations to analyze the coupled switching dynamics of polarization and cycloidal antiferromagnetic orders in BFO. We are able to identify 14 types of transitional spin structures between two cycloids and 9 electric-field-induced spin switching paths. We demonstrate the electric-field-induced rotation of wave vectors of the cycloidal spins and discover 2 types of cycloidal spin switching dynamics: fast local spin flips and slow rotation of wave vectors. Also, we construct road maps to achieve the switching between any two spin cycloids through multi-step applications of electric fields. The work provides a theoretical framework for the phenomenological description of spin cycloids and a fundamental understanding of the switching mechanisms to achieve electrical control of magnetic orders., Comment: 23 pages, 8 figures, with supplementary materials

Flexoelectricity is a universal electromechanical coupling effect whereby all dielectric materials polarize in response to strain gradients. In particular, nanoscale flexoelectricity promises exotic phenomena and functions, but reliable characterization methods are required to unlock its potential. Here, we report anomalous mechanical control of quantum tunnelling that allows for characterizing nanoscale flexoelectricity. By applying strain gradients with an atomic force microscope tip, we systematically polarize an ultrathin film of otherwise nonpolar SrTiO$_{3}$, and simultaneously measure tunnel current across it. The measured tunnel current exhibits critical behaviour as a function of strain gradients, which manifests large modification of tunnel barrier profiles via flexoelectricity. Further analysis of this critical behaviour reveals significantly enhanced flexocoupling strength in ultrathin SrTiO$_{3}$, compared to that in bulk, rendering flexoelectricity more potent at the nanoscale. Our study not only suggests possible applications exploiting dynamic mechanical control of quantum effect, but also paves the way to characterize nanoscale flexoelectricity., Comment: Main and Supplementary combined, 40 pages

Very sensitive responses to external forces are found near phase transitions. However, phase transition dynamics and pre-equilibrium phenomena are difficult to detect and control. We have directly observed that the equilibrium domain structure following a phase transition in BaTiO3, a ferroelectric and ferroelastic material, is attained by halving of the domain periodicity, sequentially and multiple times. The process is reversible, displaying periodicity doubling as temperature is increased. This observation is backed theoretically and can explain the fingerprints of domain period multiplicity observed in other systems, strongly suggesting this as a general model for pattern formation during phase transitions in ferroelastic materials.

Electric current has been experimentally demonstrated to be able to drive the insulator-to-metal transition (IMT) in VO$_2$. The main mechanisms involved are believed to be the Joule heating effect and the strong electron-correlation effect. These effects are often entangled with each other in experiments, which complicates the understanding of the essential nature of the observations. We formulate a phase-field model to investigate theoretically in mesoscale the pure correlation effect brought by the current on the IMT in VO$_2$, i.e., the isothermal process under the current. We find that a current with a large density ($\sim 10^1$ nA/nm$^2$) induces a few-nanosecond ultrafast switch in VO$_2$, in agreement with the experiment. The temperature-current phase diagram is further calculated, which reveals that the current may induce the M2 phase at low temperatures. The current is also shown capable of driving domain walls to move. Our work may assist related experiments and provide guidance to the engineering of VO$_2$-based electric switching devices., Comment: 7 pages, 4 figures, 1 table

Two dimensional (2D) materials continue to hold great promise for future electronics, due to their atomic-scale thicknesses and wide range of tunable properties. However, commercial efforts in this field are relatively recent, and much progress is required to fully realize 2D materials for commercial success. Here, we present a roadmap for the realization of electronic-grade 2D materials. We discuss technology drivers, along with key aspects of synthesis and materials engineering required for development of these materials. Additionally, we highlight several fundamental milestones required for realization of electronic-grade 2D materials, and intend this article to serve as a guide for researchers in the field.

A thermodynamically consistent phase-field model is developed to study the non-isothermal grain coalescence during the sintering process, with a potential application to the simulation in unconventional sintering techniques, e.g. spark plasma sintering, field-assisted sintering, and selective laser sintering, where non-equilibrium and high temperature gradient exist. In the model, order parameters are adopted to represent the bulk and atmosphere/pore region, as well as the crystallographic orientations. Based on the entropy analysis, the temperature-dependent free energy density is developed, which includes contributions from the internal energy (induced by the change of temperature and order parameters) and the order parameter related configurational entropy. The temperature-dependent model parameters are determined by using the experimental data of surface and grain boundary energies and interface width. From laws of thermodynamics, the kinetics for the order parameters and the order-parameter-coupled heat transfer are derived. The model is numerically implemented by the finite element method. Grain coalescence from two identical particles shows that non-isothermal condition leads to the unsymmetric morphology and curved grain boundary due to the gradients of on-site surface and grain-boundary energies induced by the local temperature inhomogeneity. More simulations on the non-isothermal grain coalescence from two non-identical and multiple particles present the temporal evolution of grain shrinkage/growth, neck growth, and porosity, demonstrating the capability and versatility of the model. It is anticipated that the work could provide a contribution to the research community of unconventional sintering techniques that can be used to model the non-isothermal related microstructural features.

Physical systems evolve from one state to another along paths of least energy barrier. Without a priori knowledge of the energy landscape, multidimensional search methods aim to find such minimum energy pathways between the initial and final states of a kinetic process. However in many cases, the user has to repeatedly provide initial guess paths, thus ensuring that the reliability of the final result is heavily user-dependent. Recently, the idea of "distortion symmetry groups" as a complete description of the symmetry of a path has been introduced. Through this, a new framework is enabled that provides a powerful means of classifying the infinite collection of possible pathways into a finite number of symmetry equivalent subsets, and then exploring each of these subsets systematically using rigorous group theoretical methods. The method is shown to lead to the discovery of new, previously hidden pathways for the case studies of bulk ferroelectric switching and domain wall motion in proper and improper ferroelectrics, as well as in multiferroic switching. This approach is applicable to a wide variety of physical phenomena ranging from structural, electronic and magnetic distortions, diffusion, and phase transitions in materials.

Tuning the lattice degree of freedom in nanoscale functional crystals is critical to exploit the emerging functionalities such as piezoelectricity, shape-memory effect, or piezomagnetism, which are attributed to the intrinsic lattice-polar or lattice-spin coupling. Here it is reported that a mechanical probe can be a dynamic tool to switch the ferroic orders at the nanoscale multiferroic phase boundaries in BiFeO 3 with a phase mixture, where the material can be reversibly transformed between the "soft" tetragonal-like and the "hard" rhombohedrallike structures. The microscopic origin of the nonvolatile mechanical switching of the multiferroic phase boundaries, coupled with a reversible 180{\deg} rotation of the in-plane ferroelectric polarization, is the nanoscale pressure-induced elastic deformation and reconstruction of the spontaneous strain gradient across the multiferroic phase boundaries. The reversible control of the room-temperature multiple ferroic orders using a pure mechanical stimulus may bring us a new pathway to achieve the potential energy conversion and sensing applications.

Oxygen vacancies, especially their distribution, are directly coupled to the electromagnetic properties of oxides and related emergent functionalities that have implication in device applications. Here using a homoepitaxial strontium titanate thin film, we demonstrate a controlled manipulation of the oxygen vacancy distribution using the mechanical force from a scanning probe microscope tip. By combining Kelvin probe force microscopy imaging and phase-field simulations, we show that oxygen vacancies can move under a stress-gradient-induced depolarisation field. When tailored, this nanoscale flexoelectric effect enables a controlled spatial modulation. In motion, the scanning probe tip thereby deterministically reconfigures the spatial distribution of vacancies. The ability to locally manipulate oxygen vacancies on-demand provides a tool for the exploration of mesoscale quantum phenomena, and engineering multifunctional oxide devices., Comment: 35 pages, Main text and the supplementary information combined

Theoretical descriptions of the Seebeck coefficient in terms of the differential electrical conductivity given by Cutler and Mott is the foundation of later works in terms of transmission function from the thermoelectric transport theory. On the other hand, recent studies in the literature have shown the relation between the Seebeck coefficient and chemical potential of electrons. In this work, this relation is rigorously derived from fundamental thermodynamics, and an formalism for the parameter-free calculation of the Seebeck coefficient based on the electronic density-of-states from first-principles calculations is presented. Numerical results are given using the n-type La3-xTe4 thermoelectric material as the prototype. With the rigid band approximation, the calculated temperature dependences of the Seebeck coefficients of La3-xTe4 as a function of carrier concentration show excellent agreements with experimental data.

Thermoelectric effects, measured by the Seebeck coefficients, refer to the phenomena in which a temperature difference or gradient imposed across a thermoelectric material induces an electrical potential difference or gradient, and vice versa, enabling the direct conversion of thermal and electric energies. All existing understanding and first-principles calculations of Seebeck coefficients have been based on the Boltzmann kinetic transport theory. Here we demonstrate that the Seebeck coefficient is a well-defined thermodynamic quantity that can be determined from the change in the chemical potential of electrons induced by the temperature change and thus can be efficiently computed solely based on the electronic density of states through first-principles calculations at different temperatures. The proposed approach is demonstrated using the prototype PbTe and SnSe thermoelectric materials. The proposed thermodynamic approach dramatically simplifies the calculations of Seebeck coefficients, making it possible to search for high performance thermoelectric materials using high-throughput first-principles calculations.

An incommensurate phase refers to a solid state in which the period of a superstructure is incommensurable with the primitive unit cell. Recently the incommensurate phase is induced by applying an in-plane strain to hexagonal manganites, which demonstrates single chiral modulation of six domain variants. Here we employ Landau theory in combination with the phase-field method to investigate the incommensurate phase in hexagonal manganites. It is shown that the equilibrium wave length of the incommensurate phase is determined by temperature and the magnitude of the applied strain, and a temperature-strain phase diagram is constructed for the stability of the incommensurate phase. Temporal evolution of domain structures reveals that the applied strain not only produces the force pulling the vortices and anti-vortices in opposite directions, but also results in the creation and annihilation of vortex-antivortex pairs., Comment: 23 pages, 7 figures

Hexagonal manganites REMnO3 (RE, rare earths) have attracted significant attention due to their potential applications as multiferroic materials and the intriguing physics associated with the topological defects. The two-dimensional (2D) and 3D domain and vortex structure evolution of REMnO3 is predicted using the phase-field method based on a thermodynamic potential constructed from first-principles calculations. In 3D spaces, vortex lines show three types of topological changes, i.e. shrinking, coalescence, and splitting, with the latter two caused by the interaction and exchange of vortex loops. Compared to the coarsening rate of the isotropic XY model, the six-fold degeneracy gives rise to negligible differences with the vortex-antivortex annihilation controlling the scaling dynamics, whereas the anisotropy of interfacial energy results in a deviation. The temporal evolution of domain and vortex structures serves as a platform to fully explore the mesoscale mechanisms for the 0-D and 1-D topological defects., Comment: 19 pages, 5 figures

Phase decomposition is a well-known process leading to the formation of two-phase mixtures. Here we show that a strain imposed on a ferroelastic crystal promotes the formation of mixed phases and domains, i.e., domain and phase de-strain with local strains determined by a common tangent construction on the free energy versus strain curves. It is demonstrated that a domain structure can be understood using the concepts of domain/phase rule, lever rule, coherent and incoherent de-strain within the de-strain model description, in a complete analogy to phase decomposition. The proposed de-strain model is validated using phase-field simulations and experimental observations of PbTiO_3 and BiFeO_3 thin films as examples. The de-strain model provides a simple tool to guide and design domain structures of ferroelastic systems., Comment: 18 pages, 4 figures, submitted to PRL on 22Jan16, and reported on the 2016 MRS spring meeting