Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Landau theory”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Analytical model and dynamical phase-field simulation of terahertz transmission across ferroelectrics

Here, 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.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Predicting phase transitions in PbTi⁢O 3 using zentropy through quasiharmonic phonon calculations

According to x-ray diffraction (XRD) measurements, PbTi⁢O 3 undergoes a phase transition from a tetragonal ferroelectric (FE) phase to a cubic paraelectric phase at 763 K. However, x-ray absorption fine-structure (XAFS) measurements indicate that PbTi⁢O 3 is locally tetragonal even after the phase transition. The difference in these results is because XAFS measurements can probe local features of a structure, while XRD averages over such local features. For both measurements to be consistent, PbTi⁢O 3 is macroscopically cubic but locally tetragonal after the phase transition. Despite this, most models, such as the Landau-Ginsburg-Devonshire theory and effective Hamiltonians, are still unable to explain this phenomenon. Moreover, these methods involve model parameters fitted to experimental or theoretical data and do not consider other tetragonal configurations, such as domain walls, to predict the phase transition. In our previous study, we used our zentropy approach to predict the phase transition by considering the tetragonal FE ground-state configuration and the tetragonal 90° and 180° domain wall configurations with their total energies at 0 K. Here, in this paper, the Helmholtz energies of the three configurations are obtained from density functional theory calculations through energy-volume curves and phonon calculations. The predicted phase transition temperature using the meta-GGA 𝑟 2⁢ SCAN and revised multiplicities of configurations is 716 K, showing good agreement with the experimental value of 763 K.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Landau-Forbidden Quantum Criticality in Rydberg Quantum Simulators

The Landau-Ginzburg-Wilson theory of phase transitions precludes a continuous transition between two phases that spontaneously break distinct symmetries. However, quantum mechanical effects can intertwine the symmetries, giving rise to an exotic phenomenon called deconfined quantum criticality (DQC). Here, in this Letter, we study the ground state phase diagram of a one-dimensional array of individually trapped neutral atoms interacting strongly via Rydberg states, and demonstrate through extensive numerical simulations that it hosts a variety of symmetry-breaking phases and their transitions including DQC. We show how an enlarged, emergent continuous symmetry arises at the DQCs, which can be experimentally observed in the joint distribution of two distinct order parameters, obtained within measurement snapshots in the standard computational basis. Our findings highlight quantum simulators of Rydberg atoms not only as promising platforms to experimentally realize such exotic phenomena, but also as unique ones allowing access to physical properties not obtainable in traditional experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Langevin Dynamics modeling of gas-phase ion-ion recombination (Final Technical Report)

A self-consistent trajectory simulation approach to model MN reactions (Fig. 1) which incorporates the probability of electron transfer as a Monte Carlo operator (Fig. 2) was developed and published as Liu et al. J. Chem. Phys. 159, 114111 (2023). The electron transfer probability p ET estimated using the two-state Landau-Zener (LZ) theory was incorporated into classical trajectory simulations to elicit predictions of MN reaction cross-section σ (vacuum) or rate constant β (finite pressure). Electronic structure calculations with multireference configuration interaction (MRCI) and large correlation consistent basis sets were used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the incorporation of electron transfer probability and the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of neutral gas on ions. For H + – H - and Li + – H(D) - pairs, our approach quantitatively agrees with measured speed-dependent cross-sections for up to ~10 5 m/s. For the ion pair Ne + – Cl - , our predictions of the MN rate constant at ~1 torr are a factor of ~2 – 3 higher than the experimentally measured value. Similarly, for Xe + – F - in the pressure range of ~20000 – 80000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. The paradigm of using trajectory simulations to self-consistently model MN reactions is the basis for inclusion of additional non-classical, and static magnetic and electric field effects. Subsequent work, published as Roy et al. focused on modeling recombination rate constant for three ion pairs (rare gas Ar + cation and halide anions): Ar + – Cl - , Ar + – Br - , Ar + – I - , 2) considering spin-orbit couplings in the electronic structure calculations to obtain high-fidelity estimates of the electron transfer probability and incorporated within the classical trajectory simulations to elicit predictions. In addition to calculations of ion-ion recombination rate constants, a classical trajectory simulation technique (published as Roy et al. J. Chem. Phys. 162(9), 094104 (2023)) that uses quaternions to represent orientation of non-spherical particles (ions or aerosol particles) was developed to simulate the recombination of diatomic or more generally, polyatomic molecules. Finally, several other ion pairs such as Ne + – Cl - , Kr + – Cl - , were explored using the developed semi-classical trajectory simulations to understand various challenges in tackling electronic structure calculations. Using empirical approaches to parameterize the electron transfer radius, trajectory simulations were also used to probe the effect of ion number density on MN rate constant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Landau discriminants

Scattering amplitudes in quantum field theories have intricate analytic properties as functions of the energies and momenta of the scattered particles. In perturbation theory, their singularities are governed by a set of nonlinear polynomial equations, known as Landau equations, for each individual Feynman diagram. The singularity locus of the associated Feynman integral is made precise with the notion of the Landau discriminant, which characterizes when the Landau equations admit a solution. In order to compute this discriminant, we present approaches from classical elimination theory, as well as a numerical algorithm based on homotopy continuation. These methods allow us to compute Landau discriminants of various Feynman diagrams up to 3 loops, which were previously out of reach. For instance, the Landau discriminant of the envelope diagram is a reducible surface of degree 45 in the three-dimensional space of kinematic invariants. We investigate geometric properties of the Landau discriminant, such as irreducibility, dimension and degree. In particular, we find simple examples in which the Landau discriminant has codimension greater than one. Furthermore, we describe a numerical procedure for determining which parts of the Landau discriminant lie in the physical regions. In order to study degenerate limits of Landau equations and bounds on the degree of the Landau discriminant, we introduce Landau polytopes and study their facet structure. Finally, we provide an efficient numerical algorithm for the computation of the number of master integrals based on the connection to algebraic statistics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Topological Charge Pumping in Excitonic Insulators

Here, we show that in excitonic insulators with s-wave electron-hole pairing, an applied electric field (either pulsed or static) can induce a p-wave component to the order parameter, and further drive it to rotate in the s + ip plane, realizing a Thouless charge pump. In one dimension, each cycle of rotation pumps exactly two electrons across the sample. Higher dimensional systems can be viewed as a stack of one-dimensional chains in momentum space in which each chain crossing the Fermi surface contributes a channel of charge pumping. Physics beyond the adiabatic limit, including in particular dissipative effects is discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Vortex dynamics induced by scanning SQUID susceptometry

In this work, we measured the local magnetic response of a niobium thin film by applying a millitesla-scale AC magnetic field using a micron-scale field coil and detecting the response with a micron-scale pickup loop in a scanning superconducting quantum interference device (SQUID) susceptometry measurement. Near the film's critical temperature, we observed a steplike nonlinear and dissipative magnetic response due to the dynamics of a small number of vortex-antivortex pairs induced in the film by the local applied AC field. We modeled the dynamics of the measurement using a combined two-dimensional London-Maxwell and time-dependent Ginzburg-Landau approach, allowing us to construct a detailed real-space picture of the vortex motion causing the observed dissipative response. This work pushes scanning SQUID susceptometry of two-dimensional superconductors beyond the regime of linear response and lays the foundation for microscopic studies of vortex dynamics and pinning in superconducting devices and more exotic materials systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number (n) and frequency (ω) in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field (E||), ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfv\'en wave (KAW), ion sound wave (ISW), low-n kink, high-n ion temperature gradient mode (ITG) and kinetic ballooning mode (KBM). Moreover, MAS is successfully applied to model the Alfv\'en eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of reversed shear Alfv\'en eigenmode (RSAE), the tunneling damping of toroidal Alfv\'en eigenmode (TAE), as well as the polarization characteristics of kinetic beta-induced Alfv\'en eigenmode (KBAE) and beta-induced Alfv\'en-acoustic eigenmode (BAAE) being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1981

Two-component Fermi-liquid theory - Equilibrium properties of liquid metallic hydrogen

It is reported that the transition of condensed hydrogen from an insulating molecular crystal phase to a metallic liquid phase, at zero temperature and high pressure, appears possible. Liquid metallic hydrogen (LMH), comprising interpenetrating proton and electron fluids, would constitute a two-component Fermi liquid with both a very high component-mass ratio and long-range, species-dependent bare interactions. The low-temperature equilibrium properties of LMH are examined by means of a generalization to the case of two components of the phenomenological Landau Fermi-liquid theory, and the low-temperature specific heat, compressibility, thermal expansion coefficient and spin susceptibility are given. It is found that the specific heat and the thermal expansion coefficient are vastly greater in the liquid than in the corresponding solid, due to the presence of proton quasiparticle excitations in the liquid.

Oliva, J.↗
DOE OSTI2021

Non-perturbative BRST symmetry and the spectral structure of the ghost propagator

In BRST-quantised Yang-Mills theory the existence of BRST symmetry imposes significant constraints on the analytic structure of the continuum theory. In particular, the presence of this symmetry in the non-perturbative regime implies that any on-shell state with vanishing norm must have an associated partner state with identical mass, but negative inner product. In this work we demonstrate that for quantum chromodynamics (QCD) this constraint gives rise to an interconnection between the ghost and gluon spectra, and in doing so provides a non-trivial test of whether BRST symmetry is realised non-perturbatively. By analysing infrared lattice data for the minimal Landau gauge ghost propagator in pure SU(3) Yang-Mills theory, and comparing this with previous results for the gluon propagator, we show that this interconnection is violated, and hence conclude that continuum and current lattice formulations of Yang-Mills theory in Landau gauge represent two distinct realisations of the theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Q-POP-Thermo: A general-purpose thermodynamics solver for ferroelectric materials

We report that Q-POP-Thermo is a program designed to compute thermodynamic monodomain equilibrium states and their properties for ferroelectric single crystals and thin films based on the Landau-Ginzburg-Devonshire (LGD) Theory. Utilizing symbolic manipulation with the SymPy Library, the governing equations along with appropriate boundary conditions are solved for speedy minimization of the free energy of a crystal. Utilizing the popular Differential Evolution algorithm, with appropriate hybridization, multiple phase diagrams, such as the pressure-temperature phase diagram for bulk single crystals and the common strain-temperature phase diagram for monodomain thin-film systems can be readily generated. Furthermore, a variety of material properties of stable ferroelectric phases, including dielectric, piezoelectric, and electrocaloric properties, can simultaneously be calculated. Validation studies are presented for both thin-film and single crystal systems to test the effectiveness and capability of the open-source program.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Electron– K -phonon interaction in twisted bilayer graphene

We develop an analytic theory to describe the interaction between electrons and K phonons and study its influence on superconductivity in the bare bands of twisted bilayer graphene (TBG). We find that, due to symmetry and the two-center approximation, only one optical K phonon (160meV) of graphene is responsible for the intervalley electron-phonon interaction. Here, this phonon has recently been found in angular-resolved photoemission spectroscopy to be responsible for replicas of the TBG flat bands. By projecting the interaction to the TBG flat bands, we perform the full symmetry analysis of the phonon-mediated attractive interaction and pairing channels in the Chern basis, and show that several channels are guaranteed to have gapless order parameters. From the linearized gap equations, we find that the highest T c pairing induced by this phonon is a singlet gapped s-wave inter-Chern-band order parameter, followed closely by a gapless nematic d-wave intra-Chern-band order parameter. We justify these results analytically, using the topological heavy-fermion mapping of TBG which has allowed us to obtain an analytic form of a phonon-mediated attractive interaction and to analytically solve the linearized and T = 0 gap equations. For the intra-Chern-band channel, the nematic state with nodes is shown to be stabilized in the chiral flat-band limit. While the flat-band Coulomb interaction can be screened sufficiently enough—around the Van Hove singularities—to allow for electron-phonon based superconductivity, it is unlikely that this effect can be maintained in the lower density of states excitation bands around the correlated insulator states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Atomic Alignment in PbS Nanocrystal Superlattices with Compact Inorganic Ligands via Reversible Oriented Attachment of Nanocrystals

Nanocrystals (NCs) serve as versatile building blocks for the creation of functional materials with NC self-assembly offering opportunities to enable novel material properties. Here, in this work, we demonstrate that PbS NCs functionalized with strongly negatively charged metal chalcogenide complex (MCC) ligands, such as Sn 2 S 6 4– and AsS 4 3– , can self-assemble into all-inorganic superlattices with both long-range superlattice translational and atomic-lattice orientational order. Structural characterizations reveal that the NCs adopt an unexpected edge-to-edge alignment, and numerical simulation clarifies that orientational order is thermodynamically stabilized by many-body ion correlations originating from the dense electrolyte. Furthermore, we show that the superlattices of Sn 2 S 6 4– -functionalized PbS NCs can be fully disassembled back into the colloidal state, which is highly unusual for orientationally attached superlattices with atomic-lattice alignment. The reversible oriented attachment of NCs, enabling their dynamic assembly and disassembly into effectively single-crystalline superstructures, offers a pathway toward designing reconfigurable materials with adaptive and controllable electronic and optoelectronic properties.

Landau-Ginzburg theory↗
DOE OSTI2020

Scaling behavior of low-temperature orthorhombic domains in the prototypical high-temperature superconductor La 1.875 Ba 0.125 CuO 4

Structural symmetry breaking and recovery in condensed-matter systems are closely related to exotic physical properties such as superconductivity (SC), magnetism, spin density waves, and charge density waves (CDWs). The interplay between different order parameters is intricate and often subject to intense debate, as in the case of CDW order and superconductivity. In La 1.875 Ba 0.125 CuO 4 (LBCO), the low-temperature structural domain walls are hypothesized as nanometer-scale pinning sites for the CDWs. Coherent x-ray diffraction techniques have been employed here to visualize the domain structures associated with these symmetry changes directly during phase transition. We have pushed Bragg coherent diffractive imaging (BCDI) into the cryogenic regime where most phase transitions in quantum materials reside. Utilizing BCDI, we image the structural evolution of LBCO microcrystal samples during the high-temperature tetragonal to low-temperature orthorhombic (LTO) phase transition. Our results show the formation of LTO domains close to the transition temperature and how the domain size decreases with temperature. In conclusion, the number of domains follows the secondary order parameter (or orthorhombic strain) measurement with a critical exponent that is consistent with the three-dimensional universality class.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Phenomenological description of bright domain walls in ferroelectric-antiferroelectric layered chalcogenides

Recently, a layered ferroelectric CuIn P 2 Se 6 was shown to exhibit domain walls with locally enhanced piezoresponse—a striking departure from the observations of nominally zero piezoresponse in most ferroelectrics. Although it was proposed that such "bright" domain walls are phase boundaries between ferri- and antiferroelectrically ordered regions of the materials, the physical mechanisms behind the existence and response of these boundaries remain to be understood. Here, using Landau-Ginzburg-Devonshire phenomenology combined with a four sublattices model, we describe quantitatively the bright-contrast and dark-contrast domain boundaries between the antiferroelectric, ferroelectric, or ferrielectric long-range-ordered phases in a layered ferroelectric-antiferroelectric ferroics, such as CuIn P 2 (S 1-y Se y ) 6 .

36 MATERIALS SCIENCE↗
DOE OSTI2020

Coexistence of antiferromagnetism and topological superconductivity on the honeycomb lattice Hubbard model

In this report, motivated by the recent numerical simulations for doped t-J model on the honeycomb lattice, we study superconductivity of singlet and triplet pairing on the honeycomb lattice Hubbard model. We show that a superconducting state with coexisting spin-singlet and spin-triplet pairings is induced by the antiferromagnetic order near half filling. The superconducting state we obtain has a topological phase transition that separates a topologically trivial state and a nontrivial state with Chern number two. Possible experimental realization of such a topological superconductivity is also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ultrafast formation of domain walls of a charge density wave in SmTe 3

We study ultrafast x-ray diffraction on the charge density wave (CDW) of SmTe 3 using an x-ray free-electron laser. The high momentum and time resolution afforded by the x-ray laser enabled capturing fine wave-vector and time-dependent features of the CDW that originate from fast (in time) and sharp (in real space) variations of the CDW lattice distortion, which we attribute to an inversion of the order parameter. These domain inversions occur near the surface and are caused by the short penetration depth of the near-infrared pump with the wavelength centered at 800 nm, resulting in CDW domain walls perpendicular to the sample surface. These domain walls break the CDW long-range order on the scale of the x-ray probe depth, controlled experimentally by the x-ray incidence angle and suppress the diffraction intensity of the CDW for times much longer than the ~1–ps recovery of the electronic gap observed in time and angle-resolved photoemission spectroscopy. We model the spatial and temporal dependences of the order parameter using a simple Ginzburg-Landau model with all the parameters obtained from the published literature. We find reasonable agreement between the calculated and the measured diffraction across the momentum, time, fluence, and incidence angle dependence without adjusting any parameters. In conclusion, we reconstruct the spatial and temporal dependences of the lattice order parameter and find that at long times, depending on the pump fluence, multiple domain walls remain at distances of a few nanometers from the surface.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Inhomogeneous time-reversal symmetry breaking in Sr 2 RuO 4

We show that the observed time-reversal symmetry breaking (TRSB) of the superconducting state in Sr 2 RuO 4 can be understood as originating from inhomogeneous strain fields near edge dislocations of the crystal. Specifically, we argue that, without strain inhomogeneities, Sr 2 RuO 4 is a single-component, time-reversal symmetric superconductor, likely with d x 2 -y 2 symmetry. However, due to the strong strain inhomogeneities generated by dislocations, a slowly decaying subleading pairing state contributes to the condensate in significant portions of the sample. As it phase winds around the dislocation, time-reversal symmetry is locally broken. Global phase locking and TRSB occur at a sharp Ising transition that is not accompanied by a change of the single-particle gap and yields a very small heat capacity anomaly. Our model thus explains the puzzling absence of a measurable heat capacity anomaly at the TRSB transition in strained samples and the dilute nature of the time-reversal symmetry broken state probed by muon spin rotation experiments. We propose that plastic deformations of the material may be used to manipulate the onset of broken time-reversal symmetry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗