SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Structure characterization and classification is frequently based on local environment information of all or selected atomic sites in the crystal structure.
Several families of layered perovskite oxide ferroelectrics exhibit a coupling between polarization and structural order parameters, such as octahedral rotation distortions. This coupling provides opportunities for novel electric field-based manipulation of material properties, and it also stabilizes complex domain patterns and domain-wall vortices. Among layered perovskites with such coupled orders, the Aurivillius-phase oxides SrBi 2 B 2 O 9 (B = Ta,Nb) are well-known for their excellent room-temperature ferroelectric performance. This work combines group theoretic analysis with density functional theory calculations to examine the ferroelectric switching processes of SrBi2B2O9. Low-energy two-step ferroelectric switching paths are identified, with polarization reversal facilitated by structural order-parameter rotations. Analysis of the domain structure reveals how the relative energetics of the coupled order parameters translates into a network of several distinct domain-wall types linked by domain-wall vortex structures. Comparisons are made between the ferroelectric switching and domain structure of SrBi 2 B 2 O 9 and those of the layered n = 2 Ruddlesden-Popper hybrid improper ferroelectrics. The results provide insight into how ferroelectric properties may be optimized by engineering the complex crystal structures of Aurivillius-phase oxides.
Understanding the driving forces behind the nucleation of different polymorphs is of great importance for material sciences and the pharmaceutical industry. This includes understanding the reaction coordinate that governs the nucleation process and correctly calculating the relative free energies of different polymorphs. Here, we demonstrate, for the prototypical case of urea nucleation from the melt, how one can learn such a one-dimensional reaction coordinate as a function of prespecified order parameters and use it to perform efficient biased all-atom molecular dynamics simulations. The reaction coordinate is learnt as a function of the generic thermodynamic and structural order parameters using the “spectral gap optimization of order parameters (SGOOP)” approach [Tiwary, P. and Berne, B. J. Proc. Natl. Acad. Sci. U.S.A. (2016)] and is biased using well-tempered metadynamics simulations. The reaction coordinate gives insights into the role played by different structural and thermodynamics order parameters, and the biased simulations obtain accurate relative free energies for different polymorphs. This includes an accurate prediction of the approximate pressure at which urea undergoes a phase transition and one of the metastable polymorphs becomes the most stable conformation. Here, we believe the ideas demonstrated in this work will facilitate efficient sampling of nucleation in complex, generic systems.
Alloying and site ordering play complementary roles in dictating a material’s properties. However, deconvolving the impacts of these separate phenomena can be challenging. In this work, we simulate structures of Zn(Sn,Ge)N 2 with varied Sn content and site ordering to determine the impacts of order and composition on structural and electronic properties. We assess the formation enthalpy, lattice parameters, elastic constants, spontaneous polarization, and piezoelectric coefficients. In mostly disordered structures (order parameters ranging from 0.2 to 0.4), the formation enthalpy exhibits local extrema as a function of the order parameter, deviating from the more linear trends seen in both fully disordered and fully ordered systems. This anomalous deviation from the otherwise linear trend in formation enthalpy with order manifests in each of the other properties calculated. This range of order parameters of interest may be caused by a transition in the ordering of the quaternary material similar to phase changes seen in ternary compounds but stretched over a region inclduing 20% of the order parameter range. Most parameters calculated are more sensitive to order than to composition in the limited composition range tested; however, the lattice parameter c, piezoelectric coefficient e 33 , and elastic moduli C 12 , C 13 , and C 23 are more sensitive to composition. Of the properties compared, the piezoelectric coefficients are influenced most significantly by changes in both the composition and order parameter. Lattice parameters undergo the smallest changes with order and composition, but these small differences appear to impart large trends in the other properties. Better understanding the effects of disorder and group IV alloying in Zn(Sn,Ge)N 2 allows for more accurate modeling of characteristics of this material system for solid state lighting and other applications.
The excitonic insulator is an electronically driven phase of matter that emerges upon the spontaneous formation and Bose condensation of excitons. Detecting this exotic order in candidate materials is a subject of paramount importance, as the size of the excitonic gap in the band structure establishes the potential of this collective state for superfluid energy transport. However, the identification of this phase in real solids is hindered by the coexistence of a structural order parameter with the same symmetry as the excitonic order. Only a few materials are currently believed to host a dominant excitonic phase, Ta 2 NiSe 5 being the most promising. Here, we test this scenario by using an ultrashort laser pulse to quench the broken-symmetry phase of this transition metal chalcogenide. Tracking the dynamics of the material’s electronic and crystal structure after light excitation reveals spectroscopic fingerprints that are compatible only with a primary order parameter of phononic nature. We rationalize our findings through state-of-the-art calculations, confirming that the structural order accounts for most of the gap opening. Our results suggest that the spontaneous symmetry breaking in Ta 2 NiSe 5 is mostly of structural character, hampering the possibility to realize quasi-dissipationless energy transport.
A comprehensive framework of characterizing complex self-assembled structures with a set of orientational order parameters is presented. It is especially relevant in the context of using anisotropic building blocks with various symmetries. Two classes of tensor order parameters are associated with polyhedral nematic and bond orientational order. For the latter, a variation of classical bond order parameters (BOPs) is introduced, which takes advantage of the symmetry of constituent particles and/or expected crystalline phases. These symmetrized BOPs can be averaged over an entire system or assigned locally to an individual bond. By combining that with the bond percolation procedure, one is able to identify coherent domains within a self-assembled structure. As a demonstration of the proposed framework, we apply it to a simulated hybrid system that combines isotropic and patchy particles with octahedral symmetry. In conclusion, not only does the methodology allow one to identify individual crystalline domains but also it detects coherent clusters of a peculiar compact amorphous structure that is not space-filling and lacks any long-range order.
In rare earth nickelates, the metal-to-insulator transition observed as a function of temperature can be described using an electronic and a structural order parameter. The electronic one is linked to the electronic disproportionation observed below the transition temperature and the structural one characterizes the breathing mode that develops in the low temperature phase. Here SmNiO 3 /NdNiO 3 superlattices are grown with insulating LaAlO 3 spacer layers to study the unusual coupling of metal-to- insulator transitions observed at SmNiO 3 /NdNiO 3 interfaces and determine the role of the two order parameters. Temperature-dependent transport measurements reveal that a single unit cell of LaAlO 3 inserted between SmNiO 3 and NdNiO 3 leads to a complete decoupling of the metal-to-insulator transitions suggesting that the order parameter controlling the coupling is the electronic one.
Spin density waves, based on modulated local moments, are usually associated with metallic materials, but have recently been reported in insulators which display coupled magnetic and structural order parameters. Here, we discuss one such example, the multiferroic Cu 3 Nb 2 O 8 , which is reported to undergo two magnetic phase transitions, first to an unknown antiferromagnetic phase at $T_N$ ≈ 26.5 K , and then to a helicoidal structure coupled to an electric polarization below $T_2$ ≈ 24 K which breaks the crystallographic inversion symmetry. By analogy with other complex oxides, one might naturally expect this intermediate phase to be a spin density wave phase. We apply spherical polarimetry to confirm the low-temperature magnetic structure, yet only observe a single magnetic phase transition to helicoidal order. We argue that the reported unknown phase actually supports an imitation spin density wave which originates from a decoupling of the components of the magnetic order parameter, as allowed by symmetry and driven by thermal fluctuations. This provides a mechanism for the magnetic, but not nuclear, structure to break inversion symmetry thereby creating an intermediate phase in the proximity of $T_N$ which imitates a spin density wave. As the temperature is reduced, this intermediate structure destabilizes the crystal such that a structural chirality is induced, as reflected by the emergence of the electric polarization, and the imitation spin density wave relaxes into a generic helicoid. This scenario in which critical fluctuations allow the magnetic structure to break inversion symmetry while the crystal structure remains centrosymmetric might be relevant to other complex multiferroics.
Even though nucleation is ubiquitous in different science and engineering problems, investigating nucleation is extremely difficult due to the complicated ranges of time and length scales involved. In this work, we simulate NaCl nucleation in both molten and aqueous environments using enhanced sampling of all-atom molecular dynamics with deep-learning-based estimation of reaction coordinates. By incorporating various structural order parameters and learning the reaction coordinate as a function thereof, we achieve significantly improved sampling relative to traditional ad hoc descriptions of what drives nucleation, particularly in an aqueous medium. Our results reveal a one-step nucleation mechanism in both environments, with reaction coordinate analysis highlighting the importance of local ion density in distinguishing solid and liquid states. However, although fluctuations in the local ion density are necessary to drive nucleation, they are not sufficient. Our analysis shows that near the transition states, descriptors such as enthalpy and local structure become crucial. Furthermore, our protocol proposed here enables robust nucleation analysis and phase sampling and could offer insights into nucleation mechanisms for generic small molecules in different environments.
LuFe 4 Ge 2 crystallizes in the ZrFe 4 Si 2 -type structure, hosting chains of Fe tetrahedra giving rise to geometric frustration and low dimensionality. The compound orders antiferromagnetically at around 36 K accompanied by a simultaneous structural transition from a tetragonal phase to an orthorhombic phase. The hydrostatic pressure dependence of the magnetic and structural transitions is investigated using electrical transport, AC magnetic susceptibility, AC calorimetry, Mössbauer, muon-spin relaxation (μSR), and x-ray-diffraction measurements. External pressure suppresses the first-order transition to the antiferromagnetic phase (AFM1) around 1.8 GPa. The structural transition is largely unaffected by pressure and remains between 30 to 35 K for pressures up to 2 GPa. A second antiferromagnetic phase (AFM2) is observed at higher pressures. The transition from the paramagnetic to the AFM2 phase is of second-order nature and appears to be connected to the structural transition. The magnetic volume fraction obtained from μSR and Mössbauer measurements reveal that the entire sample undergoes magnetic ordering in both magnetic phases. In addition, similar low-temperature muon-precession frequencies in AFM1 and AFM2 phases point at similar ordered moments and magnetic structures in both phases. Our results further indicate enhanced magnetic fluctuations in the pressure-induced AFM2 phase. The experimental observations together with density functional theory calculations suggest that the magnetic- and structural-order parameters in LuFe 4 Ge 2 are linked by magnetic frustration, causing the simultaneous magnetostructural transition.
Results of switching behavior of the improper ferroelectric LuFeO 3 are presented. Here, using a model set of films prepared under controlled chemical and growth-rate conditions, it is shown that defects can reduce the quasi-static switching voltage by up to 40% in qualitative agreement with first-principles calculations. Switching studies show that the coercive field has a stronger frequency dispersion for the improper ferroelectrics compared to a proper ferroelectric such as PbTiO 3 . It is concluded that the primary structural order parameter controls the switching dynamics of such improper ferroelectrics.
The interaction of many-body systems with intense light pulses may lead to novel emergent phenomena far from equilibrium. Recent discoveries, such as the optical enhancement of the critical temperature in certain superconductors and the photo-stabilization of hidden phases, have turned this field into an important research frontier. Here, we demonstrate nonthermal charge-density-wave (CDW) order at electronic temperatures far greater than the thermodynamic transition temperature. Using time- and angle-resolved photoemission spectroscopy and time-resolved X-ray diffraction, we investigate the electronic and structural order parameters of an ultrafast photoinduced CDW-to-metal transition. Tracking the dynamical CDW recovery as a function of electronic temperature reveals a behaviour markedly different from equilibrium, which we attribute to the suppression of lattice fluctuations in the transient nonthermal phonon distribution. A complete description of the system’s coherent and incoherent order-parameter dynamics is given by a time-dependent Ginzburg-Landau framework, providing access to the transient potential energy surfaces.
We show that unsupervised machine learning can be used to learn chemical transformation pathways from observational Scanning Transmission Electron Microscopy (STEM) data. To enable this analysis, we assumed the existence of atoms, a discreteness of atomic classes, and the presence of an explicit relationship between the observed STEM contrast and the presence of atomic units. With only these postulates, we developed a machine learning method leveraging a rotationally invariant variational autoencoder (VAE) that can identify the existing molecular fragments observed within a material. The approach encodes the information contained in STEM image sequences using a small number of latent variables, allowing the exploration of chemical transformation pathways by tracing the evolution of atoms in the latent space of the system. The results suggest that atomically resolved STEM data can be used to derive fundamental physical and chemical mechanisms involved, by providing encodings of the observed structures that act as bottom-up equivalents of structural order parameters. The approach also demonstrates the potential of variational (i.e., Bayesian) methods in the physical sciences and will stimulate the development of more sophisticated ways to encode physical constraints in the encoder–decoder architectures and generative physical laws and causal relationships in the latent space of VAEs.
Abstract A graph-based order parameter, based on the topology of the graph itself, is introduced for the characterization of atomistic structures. The order parameter is universal to any material/chemical system and is transferable to all structural geometries. Four sets of data are used to validate both the generalizability and accuracy of the algorithm: (1) liquid lithium configurations spanning up to 300 GPa, (2) condensed phases of carbon along with nanotubes and buckyballs at ambient and high temperature, (3) a diverse set of aluminum configurations including surfaces, compressed and expanded lattices, point defects, grain boundaries, liquids, nanoparticles, all at nonzero temperatures, and (4) eleven niobium oxide crystal phases generated with ab initio molecular dynamics. We compare our proposed method to existing, state-of-the-art methods for the cases of aluminum and niobium oxide. Our order parameter uniquely classifies every configuration and outperforms all studied existing methods, opening the door for its use in a multitude of complex application spaces that can require fine structure-level characterization of atomistic graphs.
Density-functional theory (DFT) is used to identify phase-equilibria in multi-principal-element and high-entropy alloys (MPEAs/HEAs), including duplex-phase and eutectic microstructures. Here, a combination of composition-dependent formation energy and electronic-structure-based ordering parameters were used to identify a transition from FCC to BCC favoring mixtures, and these predictions experimentally validated in the Al-Co-Cr-Cu-Fe-Ni system. A sharp crossover in lattice structure and dual-phase stability as a function of composition were predicted via DFT and validated experimentally. The impact of solidification kinetics and thermodynamic stability was explored experimentally using a range of techniques, from slow (castings) to rapid (laser remelting), which showed a decoupling of phase fraction from thermal history, i.e., phase fraction was found to be solidification rate-independent, enabling tuning of a multi-modal cell and grain size ranging from nanoscale through macroscale. Strength and ductility tradeoffs for select processing parameters were investigated via uniaxial tension and small-punch testing on specimens manufactured via powder-based additive manufacturing (directed-energy deposition). This work establishes a pathway for design and optimization of next-generation multiphase superalloys via tailoring of structural and chemical ordering in concentrated solid solutions.
A physical understanding of the nature of spatially modulated phases (SMPs) in rare-earth-doped antiferrodistortive (AFD) multiferroics and how they behave close to surfaces and interfaces is lacking. Here the emergence of the antiferroelectric (AFE), ferroelectric (FE), or ferrielectric (AFE-FE) spatial modulation in the vicinity of the morphotropic phase transition in La x Bi 1 – x Fe O 3 ( x ~ 0.2 ) is explored on the atomic level using high-resolution scanning transmission electron microscopy (HRSTEM). The suppression, or “melting,” of the AFE-type SMP in the vicinity of the AFD twin wall/surface junction is revealed by HRSTEM in La 0.22 Bi 0.78 Fe O 3 films and explained by the hybrid approach combining Landau-Ginzburg-Devonshire (LGD) phenomenology and the semimicroscopic four-sublattice model (FSM). The LGD-FSM approach reduces the problem of AFE (or AFE-FE) SMP emergence and stability to the thermodynamic analysis of the free-energy functional with AFE, FE, and AFD long-range order parameters and two master parameters: the FE-AFE coupling strength between four neighboring A sites and the nonstoichiometry factor, which are proportional to the variations of La concentration in La x Bi 1 – x Fe O 3 films. Finally, we establish that the surface-induced melting of SMPs and the associated broadening of AFE AFD domain walls minimize the film free energy under certain conditions imposed on the master parameters and gradient energy below the critical value. The observed behavior provides insight into the origin of SMPs in AFD multiferroics.
We introduce the group-equivariant autoencoder (GE autoencoder), a deep neural network (DNN) method that locates phase boundaries by determining which symmetries of the Hamiltonian have spontaneously broken at each temperature. We use group theory to deduce which symmetries of the system remain intact in all phases, and then use this information to constrain the parameters of the GE autoencoder such that the encoder learns an order parameter invariant to these “never-broken” symmetries. This procedure produces a dramatic reduction in the number of free parameters such that the GE-autoencoder size is independent of the system size. We include symmetry regularization terms in the loss function of the GE autoencoder so that the learned order parameter is also equivariant to the remaining symmetries of the system. By examining the group representation by which the learned order parameter transforms, we are then able to extract information about the associated spontaneous symmetry breaking. We test the GE autoencoder on the 2D classical ferromagnetic and antiferromagnetic Ising models, finding that the GE autoencoder (1) accurately determines which symmetries have spontaneously broken at each temperature; (2) estimates the critical temperature in the thermodynamic limit with greater accuracy, robustness, and time efficiency than a symmetry-agnostic baseline autoencoder; and (3) detects the presence of an external symmetry-breaking magnetic field with greater sensitivity than the baseline method. Lastly, we describe various key implementation details, including a quadratic-programming-based method for extracting the critical temperature estimate from trained autoencoders and calculations of the DNN initialization and learning rate settings required for fair model comparisons.