Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “condensed phase systems”

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.

At least 55 records · Page 3

Understanding Adsorption and Reactions at Aqueous Oxide Interfaces with Neural Network Potential Molecular Dynamics

Chemical processes at metal oxide−water interfaces are of central importance in geochemistry, biology, and energy technologies. A better understanding of these processes would allow us to make a significant step toward optimizing and controlling them, which could in turn lead to broader impacts. Computational modeling is indispensable to accomplishing this task because complexity and disorder often make it difficult to extract atomistic information from experiments. Balancing computational cost and accuracy, simulation schemes based on efficient machine learning representations of the potential energy surface (PES) predicted by ab initio calculations have become increasingly popular over the past decade. In particular, several studies have demonstrated the ability of machine learning models to accurately reproduce the complex ab initio PESs of aqueous oxide interfaces, allowing simulations of systems and processes that are not accessible with ab initio methods. In this Account, we review our recent efforts to understand adsorption processes and reactions at aqueous oxide interfaces using deep potential molecular dynamics (DPMD), a simulation scheme employing deep neural networks (DNNs), which has proven to be quite successful in accurately describing many different systems in the condensed phase. After summarizing the DPMD methodology, we first review our work on the acid−base chemistry of oxide surfaces in contact with water, a fundamental characteristic that controls proton transfer and surface charge at the interface. We focus on the aqueous interface of rutile IrO 2 , an oxide material thus far considered the best catalyst for the oxygen evolution reaction (OER). We show that this interface is characterized by a large fraction of dissociated water and a strong Brønsted acidity of the surface sites, in good agreement with the experimentally measured value of the point of zero proton charge. In our second example, we investigate how the adsorption of organic species from ambient air or water affects the structure and wettability of the aqueous interfaces of TiO 2 , a prototypical photocatalytic material. This is a question that is relevant to understanding the UV-induced hydrophilicity of TiO 2 surfaces, a property at the basis of self-cleaning windows and related applications. Specifically focusing on formic and acetic acids, the two most common atmospheric organic acids, our simulations reveal that these acids control the wettability of TiO 2 largely through acid−base chemistry at the interface rather than chemisorption on the oxide surface, a finding that could help improve the design of self-cleaning surfaces and photocatalytic devices. Finally, we review our recent study of methanol at TiO 2 −water interfaces, a system whose interest is largely motivated by the role of methanol in enhancing photocatalytic hydrogen evolution on TiO 2 . Our simulations provide mechanistic insights into the coupled roles of the organic adsorbate and water at the TiO 2 interface, with implications for how methanol enhances the activity of H 2 evolution.

adsorption↗

eReaxFF force field development for BaZr 0.8 Y 0.2 O 3-δ solid oxide electrolysis cells applications

The use of solid-oxide materials in electrocatalysis applications, especially in hydrogen-evolution reactions, is promising. However, further improvements are warranted to overcome the fundamental bottlenecks to enhancing the performance of solid-oxide electrolysis cells (SOECs), which is directly linked to the more-refined fundamental understanding of complex physical and chemical phenomena and mass exchanges that take place at the surfaces and in the bulk of electrocatalysis materials. Here, we developed an eReaxFF force field for barium zirconate doped with 20 mol% of yttrium, BaZr 0.8 Y 0.2 O 3-δ (BZY20) to enable a systematic, large-length-scale, and longer-timescale atomistic simulation of solid-oxide electrocatalysis for hydrogen generation. All parameters for the eReaxFF were optimized to reproduce quantum-mechanical (QM) calculations on relevant condensed phase and cluster systems describing oxygen vacancies, vacancy migrations, electron localization, water adsorption, water splitting, and hydrogen generation on the surfaces of the BZY20 solid oxide. Using the developed force field, we performed both zero-voltage (excess electrons absent) and non-zero-voltage (excess electrons present) molecular dynamics simulations to observe water adsorption, water splitting, proton migration, oxygen-vacancy migrations, and eventual hydrogen-production reactions. Based on investigations offered in the present study, we conclude that the eReaxFF force field-based approach can enable computationally efficient simulations for electron conductivity, electron leakage, and other non-zero-voltage effects on the solid oxide materials using the explicit-electron concept. Moreover, we demonstrate how the eReaxFF force field-based atomistic-simulation approach can enhance our understanding of processes in SOEC applications and potentially other renewable-energy applications.

08 HYDROGEN↗

Thermally reconfigurable monoclinic nematic colloidal fluids

Fundamental relationships are believed to exist between the symmetries of building blocks and the condensed matter phases that they form. For example, constituent molecular and colloidal rods and disks impart their uniaxial symmetry onto nematic liquid crystals, such as those used in displays. Low-symmetry organizations could form in mixtures of rods and disks, but entropy tends to phase-separate them at the molecular and colloidal scales, whereas strong elasticity-mediated interactions drive the formation of chains and crystals in nematic colloids. To have a structure with few or no symmetry operations apart from trivial ones has so far been demonstrated to be a property of solids alone, but not of their fully fluid condensed matter counterparts, even though such symmetries have been considered theoretically and observed in magnetic colloids. In this work, we show that dispersing highly anisotropic charged colloidal disks in a nematic host composed of molecular rods provides a platform for observing many low-symmetry phases. Depending on the temperature, concentration and surface charge of the disks, we find nematic, smectic and columnar organizations with symmetries ranging from uniaxial to orthorhombic and monoclinic. With increasing temperature, we observe unusual transitions from less- to more-ordered states and re-entrant phases. Most importantly, we demonstrate the presence of reconfigurable monoclinic colloidal nematic order, as well as the possibility of thermal and magnetic control of low-symmetry self-assembly. Our experimental findings are supported by theoretical modelling of the colloidal interactions between disks in the nematic host and may provide a route towards realizing many low-symmetry condensed matter phases in systems with building blocks of dissimilar shapes and sizes, as well as their technological applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Lagrangian formulation of nuclear–electronic orbital Ehrenfest dynamics with real-time TDDFT for extended periodic systems

Here, we present a Lagrangian-based implementation of Ehrenfest dynamics with nuclear–electronic orbital (NEO) theory and real-time time-dependent density functional theory for extended periodic systems. In addition to a quantum dynamical treatment of electrons and selected protons, this approach allows for the classical movement of all other nuclei to be taken into account in simulations of condensed matter systems. Furthermore, we introduce a Lagrangian formulation for the traveling proton basis approach and propose new schemes to enhance its application for extended periodic systems. Validation and proof-of-principle applications are performed on electronically excited proton transfer in the o-hydroxybenzaldehyde molecule with explicit solvating water molecules. These simulations demonstrate the importance of solvation dynamics and a quantum treatment of transferring protons. This work broadens the applicability of the NEO Ehrenfest dynamics approach for studying complex heterogeneous systems in the condensed phase.

Calculus of variations↗

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car–Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson’s equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. Furthermore, this is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500–1000 atoms (e.g., (H₂O)₂₅₆) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantum phase transitions in a model Hamiltonian exhibiting entangled simultaneous fermion-pair and exciton condensations

Quantum states of a novel Bose-Einstein condensate, in which both fermion-pair and exciton condensations are simultaneously present, have recently been realized theoretically in a model Hamiltonian system. Here, in this study, we identify quantum phase transitions in that model between fermion-pair and exciton condensations based on a geometric analysis of the convex set of ground-state two-particle reduced density matrices (2-RDMs). The 2-RDM set provides a finite representation of the infinite parameter space of Hamiltonians that readily reveals a fermion-pair condensate phase and two distinct exciton condensate phases, as well as the emergence of first- and second-order phase transitions as the particle number of the system is increased. The set, furthermore, shows that the fermion-exciton condensate (FEC) lies along the second-order phase transition between the exciton and fermion-pair condensate phases. The detailed information about the exciton and fermion-pair phases, the forces behind these phases, as well as their associated transitions provides additional insight into the formation of the FEC condensate, which we anticipate will prove useful in its experimental realization.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Inferring the Isotropic-Nematic Phase Transition with Generative Machine Learning

Generative machine learning models are capable of learning the phase behavior in condensed matter systems such as the Ising model. We utilize a score-based modeling procedure called thermodynamic maps to describe the isotropic-nematic phase transition in a melt of Gay-Berne ellipsoids. When trained on samples from a single temperature on either side of the phase transition, this generative machine learning approach infers effectively the nematic order parameter at intermediate temperatures. Furthermore, these results demonstrate score-based models’ ability to learn the physics of a nontrivial liquid crystal phase transition.

Critical exponents↗

A Solar Fuels Nexus: Molecules and Materials for Light-Driven Catalysis

The American Chemical Society (ACS) selects two groups of graduate students each year to plan and host a one-day symposium at each national meeting (both fall and spring).This year our Graduate Student Symposium Planning Committee (GSSPC), composed of seven students from four universities, proposal entitled “A Solar Fuels Nexus: Molecules and Materials for Light-Driven Catalysis” was selected for the “Crossroads in Chemistry” ACS Meeting that will take place March 23-26, 2023 in Indianapolis, IN. All members of our GSSPC are affiliated with the DOE Fuels from Sunlight Energy Innovation Hub, with two from the Liquid Sunlight Alliance (LiSA) and five from the Center for Hybrid Approaches in Solar Energy to Liquid Fuels (CHASE). Here we request funds to support this symposium. This symposium will highlight research progress and perspectives in the solar fuels generation field and seeks to advance the four priority research objectives (PROs) established by the Department of Energy’s Office of Basic Energy Sciences (DOE-BES) Roundtable Report that are also central to many research goals within LiSA and CHASE. The symposium will consist of research presentations from 10 invited senior researcher speakers on topics such as molecular catalyst design, computational modeling of electron transfer systems, microenvironmental effects on CO2 reduction and H2O oxidation catalysis, and intelligent design of semiconductor interfaces with ample time for discussions. These research topics fit very well with the Solar Photochemistry supported research areas of “light-driven electron and energy transfer in condensed phase and interfacial molecular systems,” “electrocatalysis and photocatalysis of solar fuels reactions,” and “semiconductor photoelectrochemistry.” More broadly, this symposium seeks to advance the DOE-BES’s mission to: “support fundamental research to understand, predict, and ultimately control matter and energy at the level of electrons, atoms, and molecules” by providing a diverse atmosphere where such research will be disseminated, discussed, and debated. There will be a strong focus on Diversity, Equity, and Inclusivity (DEI) in our symposium. Of our 10 speakers, 7 will be from underrepresented demographics in STEM, including 5 who identify as women. Furthermore, we have representatives from academia accompanied by one national lab scientist and one officer from the Office of Fossil Energy and Carbon Management at the DOE. All speakers will be holding a short DEI moment ahead of their talks. In order to support the career development of attending early career scientists, we will also be hosting a luncheon specifically for graduate students and postdocs to provide them opportunities to network with the distinguished speakers and other attendees. DOE funds for this symposium will be used to support the attendance and participation of 15 graduate students from US institutions by defraying travel and registration costs. These funds will promote engagement and conversation between early career scientists in the solar fuels field, while disseminating solar fuels research funded by and relevant to the DOE.

30 DIRECT ENERGY CONVERSION↗

The Gross-Neveu-Yukawa archipelago

We perform a bootstrap analysis of a mixed system of four-point functions of bosonic and fermionic operators in parity-preserving 3d CFTs with O(N) global symmetry. Our results provide rigorous bounds on the scaling dimensions of the O(N)-symmetric Gross-Neveu-Yukawa (GNY) fixed points, constraining these theories to live in isolated islands in the space of CFT data. We focus on the cases N = 1, 2, 4, 8, which have applications to phase transitions in condensed matter systems, and compare our bounds to previous analytical and numerical results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Quantum criticality and topology in non-equilibrium systems

A key goal of condensed matter physics research is to identify new phases of matter, and to understand the universal features of the phase transitions between them. In the past decade, physicists have uncovered a wealth of new phases with interesting surface properties, exemplified by the theoretical prediction and subsequent experimental discovery of topological insulators and superconductors. Traditional condensed matter systems are usually in a thermal equilibrium state and typically at very low temperature. Very recently, experimental advances have sparked interest in the non-equilibrium setting. Non-equilibrium systems can host new phases and phenomena with no equilibrium counterpart, and could also enable robust ways to build quantum memory devices to store and manipulate quantum information in a coherent manner. These phases and phenomena are inherently “dynamical”: they are described not by changes in the arrangement or structure of the constituent particles, but instead marked by sharp distinctions in how the particles move and exchange energy or quantum information. The discovery of robust non-equilibrium phases raises many fundamental questions: Can we develop a systematic theory of states of matter and of dynamical transitions between such states? How can such states be realized and probed experimentally? The main goal of this project was to explore the emergence of topological phases and quantum criticality (two cornerstones of modern condensed matter physics in equilibrium) in such non-equilibrium quantum systems. Specific goals included (1) using tensor networks to efficiently represent non-equilibrium states of matter and their phase transitions; (2) studying and designing new probes for periodically driven systems; and (3) developing analytic and numerical tools to analyze non-equilibrium topological phase transitions. Advances in these directions were achieved using novel techniques appropriate to study the non-equilibrium dynamics in many-body quantum systems combining strong interactions and randomness. Taken together, these results provide a new conceptual framework for understanding the emergence of quantum critical and topological properties in quantum systems far from thermal equilibrium.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A high-dimensional neural network potential for molecular dynamics simulations of condensed phase nickel and phase transitions

A high-dimensional neural network interatomic potential was developed and used in molecular dynamics simulations of condensed phase Ni and Ni systems with liquid–solid phase coexistence. The reference data set was generated by sampling the potential energy surface over a broad temperature-pressure domain using ab initio MD simulations to train a unified potential. Excellent agreement was achieved between bulk face-centered cubic nickel thermal expansion simulations and relevant experimental data. The same potential also yields accurate structures and diffusivities in the liquid state. The phase transition between liquid and solid phases was simulated using the two-phase interface method. The predicted melting point temperature is within a few kelvins of the literature value. Here, the general methodology could be applied to describe crystals with much more complex phase behaviors.

74 ATOMIC AND MOLECULAR PHYSICS↗

Reconstruction and uncertainty quantification of lattice Hamiltonian model parameters from observations of microscopic degrees of freedom

The emergence of scanning probe and electron beam imaging techniques has allowed quantitative studies of atomic structure and minute details of electronic and vibrational structure on the level of individual atomic units. These microscopic descriptors, in turn, can be associated with local symmetry breaking phenomena, representing the stochastic manifestation of the underpinning generative physical model. In this work, we explore the reconstruction of exchange integrals in the Hamiltonian for a lattice model with two competing interactions from observations of microscopic degrees of freedom and establish the uncertainties and reliability of such analysis in a broad parameter-temperature space. In contrast to other approaches, we specifically specify a loss function inherent to thermodynamic systems and utilize it to estimate uncertainty in simulated realizations of different models. As an ancillary task, we develop a machine learning approach based on histogram clustering to predict phase diagrams efficiently using a reduced descriptor space. We further demonstrate that reconstruction is possible well above the phase transition and in the regions of parameter space when the macroscopic ground state of the system is poorly defined due to frustrated interactions. This suggests that this approach can be applied to the traditionally complex problems of condensed matter physics such as ferroelectric relaxors and morphotropic phase boundary systems, spin and cluster glasses, and quantum systems once the local descriptors linked to the relevant physical behaviors are known.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Geometric transformation and three-dimensional hopping of Hopf solitons

Arising in many branches of physics, Hopf solitons are three-dimensional particle-like field distortions with nontrivial topology described by the Hopf map. Despite their recent discovery in colloids and liquid crystals, the requirement of applied fields or confinement for stability impedes their utility in technological applications. Here we demonstrate stable Hopf solitons in a liquid crystal material without these requirements as a result of enhanced stability by tuning anisotropy of parameters that describe energetic costs of different gradient components in the molecular alignment field. Nevertheless, electric fields allow for inter-transformation of Hopf solitons between different geometric embodiments, as well as for their three-dimensional hopping-like dynamics in response to electric pulses. Numerical modelling reproduces both the equilibrium structure and topology-preserving out-of-equilibrium evolution of the soliton during switching and motions. Our findings may enable myriads of solitonic condensed matter phases and active matter systems, as well as their technological applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Direct observation of magnetic ordering induced via systematic lattice disorder in artificial rhombus spin ices

It is critical to understand the effect of lattice geometry on the order parameter of a condensed matter system, as it controls phase transitions in such systems. Artificial spin ices (ASIs) are two-dimensional lattices of Ising-like nanomagnets that provide an opportunity to explore such phenomena by lithographically controlling the lattice geometry to observe its influence on magnetic ordering and frustration effects. Here we report a systematic approach to studying the effects of disorder in rhombus ASIs generated from combinations of five vertex motifs. We investigate four geometries characterized by a geometric order parameter, with symmetries ranging from periodic to quasiperiodic to random. Lorentz transmission electron microscopy data indicates magnetic domain behavior depends on chains of strongly-coupled islands in the periodic and sixfold-twinned lattices, while the behavior of the disordered lattice is dominated by vertex motifs with large configurational degeneracy. Utilizing micromagnetic simulations, a quantitative analysis of the lattice energetics showed that the experimental rotationally-demagnetized state of the disordered ASI was closer in energy to the idealized ground state compared to other periodic and twinned ASIs. Finally, our work provides a unique pathway for using degeneracy, magnetic frustration, and order to control the magnetization behavior of designer disordered systems.

36 MATERIALS SCIENCE↗

Modeling the Electronic Absorption Spectra of the Indocarbocyanine Cy3

Accurate modeling of optical spectra requires careful treatment of the molecular structures and vibronic, environmental, and thermal contributions. The accuracy of the computational methods used to simulate absorption spectra is limited by their ability to account for all the factors that affect the spectral shapes and energetics. The ensemble-based approaches are widely used to model the absorption spectra of molecules in the condensed-phase, and their performance is system dependent. The Franck–Condon approach is suitable for simulating high resolution spectra of rigid systems, and its accuracy is limited mainly by the harmonic approximation. In this work, the absorption spectrum of the widely used cyanine Cy3 is simulated using the ensemble approach via classical and quantum sampling, as well as, the Franck–Condon approach. The factors limiting the ensemble approaches, including the sampling and force field effects, are tested, while the vertical and adiabatic harmonic approximations of the Franck–Condon approach are also systematically examined. Our results show that all the vertical methods, including the ensemble approach, are not suitable to model the absorption spectrum of Cy3, and recommend the adiabatic methods as suitable approaches for the modeling of spectra with strong vibronic contributions. We find that the thermal effects, the low frequency modes, and the simultaneous vibrational excitations have prominent contributions to the Cy3 spectrum. The inclusion of the solvent stabilizes the energetics significantly, while its negligible effect on the spectral shapes aligns well with the experimental observations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multipolar condensates and multipolar Josephson effects

When single-particle dynamics are suppressed in certain strongly correlated systems, dipoles arise as elementary carriers of quantum kinetics. These dipoles can further condense, providing physicists with a rich realm to study fracton phases of matter. Whereas recent theoretical discoveries have shown that an unconventional lattice model may host a dipole condensate as the ground state, we show that dipole condensates prevail in bosonic systems due to a self-proximity effect. Our findings allow experimentalists to manipulate the phase of a dipole condensate and deliver dipolar Josephson effects, where supercurrents of dipoles arise in the absence of particle flows. The self-proximity effects can also be utilized to produce a generic multipolar condensate. The kinetics of the n -th order multipoles unavoidably creates a condensate of the ( n + 1)-th order multipoles, forming a hierarchy of multipolar condensates that will offer physicists a whole new class of macroscopic quantum phenomena.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

CHARMM-GUI Drude prepper for molecular dynamics simulation using the classical Drude polarizable force field

Explicit treatment of electronic polarizability in empirical force fields (FFs) represents an extension over a traditional additive or pairwise FF and provides a more realistic model of the variations in electronic structure in condensed phase, macromolecular simulations. To facilitate utilization of the polarizable FF based on the classical Drude oscillator model, Drude Prepper has been developed in CHARMM-GUI. Drude Prepper ingests additive CHARMM protein structures file (PSF) and pre-equilibrated coordinates in CHARMM, PDB, or NAMD format, from which the molecular components of the system are identified. These include all residues and patches connecting those residues along with water, ions, and other solute molecules. This information is then used to construct the Drude FF-based PSF using molecular generation capabilities in CHARMM, followed by minimization and equilibration. In addition, inputs are generated for molecular dynamics (MD) simulations using CHARMM, GROMACS, NAMD, and OpenMM. Validation of the Drude Prepper protocol and inputs is performed through conversion and MD simulations of various heterogeneous systems that include proteins, nucleic acids, lipids, polysaccharides, and atomic ions using the aforementioned simulation packages. Stable simulations are obtained in all studied systems, including 5 μs simulation of ubiquitin, verifying the integrity of the generated Drude PSFs. Additionally, the ability of the Drude FF to model variations in electronic structure is shown through dipole moment analysis in selected systems. Finally, the capabilities and availability of Drude Prepper in CHARMM-GUI is anticipated to greatly facilitate the application of the Drude FF to a range of condensed phase, macromolecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Asymmetry of the Ferroelectric Phase Transition in BaTiO 3

Symmetry changes during phase transformations fundamentally determine the behavior of thermodynamic systems, governing phenomena as diverse as water evaporation, fermion condensation, and epidemic spreading. Phase transitions are conventionally divided into two classes, and this classification is typically assumed to remain invariant upon reversing the transition. Here, an asymmetric phase transformation is uncovered in the ferroelectric–paraelectric transition of single-crystal BaTiO 3 , a model system for first-order transitions. Under slow temperature variation (≤0.1 °C min −1 ), thermodynamic, dielectric, and domain-structure measurements reveal that the ferroelectric-to-paraelectric transition exhibits latent heat, phase coexistence, and a discontinuous order parameter, while these signatures are absent upon cooling. Complementary phase-field simulations demonstrate similar behavior, attributing it to distinct elastic strain energy accumulation and release during heating and cooling. These findings reveal a first-order character upon heating but a second-order-like behavior upon cooling, challenging the conventional paradigm of symmetric phase-transition classification and suggesting new possibilities for ferroelectric-based energy storage.

Hershkovitz, Asaf [Technion-Israel Institute of Te↗