Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “strongly interacting quantum 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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

NESSi : The N on- E quilibrium S ystems S imulation package

The nonequilibrium dynamics of correlated many-particle systems is of interest in connection with pump–probe experiments on molecular systems and solids, as well as theoretical investigations of transport properties and relaxation processes. Nonequilibrium Green’s functions are a powerful tool to study interaction effects in quantum many-particle systems out of equilibrium, and to extract physically relevant information for the interpretation of experiments. Here, we present the open-source software package NESSi (The Non-Equilibrium Systems Simulation package) which allows to perform many-body dynamics simulations based on Green’s functions on the L-shaped Kadanoff–Baym contour. NESSi contains the library libcntr which implements tools for basic operations on these nonequilibrium Green’s functions, for constructing Feynman diagrams, and for the solution of integral and integro-differential equations involving contour Green’s functions. The library employs a discretization of the Kadanoff–Baym contour into time points and a high-order implementation of integration routines. The total integrated error scales up to $\mathcal{O}(N^{-7})$, which is important since the numerical effort increases at least cubically with the simulation time. A distributed-memory parallelization over reciprocal space allows large-scale simulations of lattice systems. We provide a collection of example programs ranging from dynamics in simple two-level systems to problems relevant in contemporary condensed matter physics, including Hubbard clusters and Hubbard or Holstein lattice models. The libcntr library is the basis of a follow-up software package for nonequilibrium dynamical mean-field theory calculations based on strong-coupling perturbative impurity solvers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Anisotropic Two-Dimensional Disordered Wigner Solid

The interplay between the Fermi sea anisotropy, electron-electron interaction, and localization phenomena can give rise to exotic many-body phases. An exciting example is an anisotropic two-dimensional (2D) Wigner solid (WS), where electrons form an ordered array with an anisotropic lattice structure. Such a state has eluded experiments up to now as its realization is extremely demanding: First, a WS entails very low densities where the Coulomb interaction dominates over the kinetic (Fermi) energy. Attaining such low densities while keeping the disorder low is very challenging. Second, the low-density requirement has to be fulfilled in a material that hosts an anisotropic Fermi sea. Here, we report transport measurements in a clean (low-disorder) 2D electron system with anisotropic effective mass and Fermi sea. The data reveal that at extremely low electron densities, when the r s parameter, the ratio of the Coulomb to the Fermi energy, exceeds ≃ 38, the current-voltage characteristics become strongly nonlinear at small dc biases. Several key features of the nonlinear characteristics, including their anisotropic voltage thresholds, are consistent with the formation of a disordered, anisotropic WS pinned by the ubiquitous disorder potential.

42 ENGINEERING↗
DOE OSTI2025

Unconventional Superconductivity Mediated by Exciton Density Wave Fluctuations

Synthetic platforms afford an unparalleled degree of controllability in realizing strongly correlated phases of matter. In this Letter, we study the possibility of electrically tunable exciton-mediated superconductivity arising in charge-imbalanced bilayer semiconductors. Focusing on the case of a bilayer semiconductor heterostructure, we identify the gating conditions required to achieve exciton density wave order within a self-consistent Hartree-Fock approximation. We analyze the role of the coupling of excitonic fluctuations to the fermionic charge carriers to find that the Goldstone mode of the density wave order can mediate attractive interactions leading to superconductivity. Furthermore, when the system is close to the density wave ordering, the interactions mediated by low-energy exciton modes can support an interlayer pair-density wave superconductor of anisotropic character. We discuss experimental signatures associated with these phenomena.

excitons↗
DOE OSTI2022

Gate-defined wires in twisted bilayer graphene: From electrical detection of intervalley coherence to internally engineered Majorana modes

Twisted bilayer graphene (TBG) realizes a highly tunable, strongly interacting system featuring superconductivity and various correlated insulating states. We establish gate-defined wires in TBG with proximity-induced spin-orbit coupling as (i) a tool for revealing the nature of correlated insulators and (ii) a platform for Majorana-based topological qubits. We show that the band structure of a gate-defined wire immersed in an intervalley coherent correlated insulator inherits electrically detectable fingerprints of symmetry breaking native to the latter. Surrounding the wire by a superconducting TBG region on one side and an intervalley coherent correlated insulator on the other further enables the formation of Majorana zero modes—possibly even at zero magnetic field depending on the precise symmetry-breaking order present. Our proposal not only introduces a highly gate-tunable topological qubit medium relying on internally generated proximity effects but can also shed light on the Cooper-pairing mechanism in TBG.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Thermal nature of confining strings

We investigate the quantum statistical properties of the confining string connecting a static fermion-antifermion pair in the massive Schwinger model. By analyzing the reduced density matrix of the subsystem located in between the fermion and antifermion, we demonstrate that as the interfermion separation approaches the string-breaking distance, the overlap between the microscopic density matrix and an effective thermal density matrix exhibits a pronounced, narrow peak, approaching unity at the onset of string breaking. This behavior reveals that the confining flux tube evolves toward a genuinely thermal state as the separation between the charges grows, even in the absence of an external heat bath. In other words, one cannot tell whether a reduced state of the subsystem arises from a surrounding heat bath or from entanglement with the rest of the system. The entanglement spectrum near the critical string-breaking distance exhibits a rapid transition from the dominance of a single state describing the confining electric string towards a strongly entangled state containing virtual fermion-antifermion pairs. Our findings establish a quantitative link between confinement, entanglement, and emergent thermality, and suggest that string breaking corresponds to a microscopic thermalization transition within the flux tube.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Solvent Dynamics Are Critical to Understanding Carbon Dioxide Dissolution and Hydration in Water

Simulations of carbon dioxide (CO 2 ) in water may aid in understanding the impact of its accumulation in aquatic environments and help advance technologies for carbon capture and utilization (via, e.g., mineralization). Quantum mechanical (QM) simulations based on static molecular models with polarizable continuum solvation poorly reproduce the energetics of CO 2 hydration to form carbonic acid in water, independent of the level of QM theory employed. Only with density-functional-theory-based molecular dynamics and rare-event sampling, followed by energy corrections based on embedded correlated wavefunction theory (in conjunction with density functional embedding theory), can a close agreement between theory and experiment be achieved. Such multilevel simulations can serve as benchmarks for simpler, less costly models, giving insight into potential errors of the latter. The strong influence of sampling/averaging over dynamical solvent configurations on the energetics stems from the difference in polarity of both the transition state and product (both polar) versus the reactant (nonpolar). When a solute undergoes a change in polarity during reaction, affecting its interaction with the solvent, careful assessment of the energetic contribution of the solvent response to this change is critical. Here, we show that static models (without structural sampling) that incorporate three explicit water molecules can yield far superior results than models with more explicit water molecules because fewer water molecules yield less configurational artifacts. Static models intelligently incorporating both explicit (molecules directly participating in the reaction) and implicit solvation, along with a proper QM theory, e.g., CCSD(T) for closed-shell systems, can close the accuracy gap between static and dynamic models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Interaction effects in graphene in a weak magnetic field

A weak perpendicular magnetic field, B, breaks the chiral symmetry of each valley in the electron spectrum of graphene, preserving the overall chiral symmetry in the Brillouin zone. Here we explore the consequences of this symmetry breaking for the interaction effects in graphene. In particular, we demonstrate that the electron-electron interaction lifetime acquires an anomalous B dependence. Also, the ballistic zero-bias anomaly, δv(ω), where ω is the energy measured from the Fermi level, emerges at a weak B and has the form δv(B) ~ B 2 /ω 2 . Temperature dependence of the magnetic-field corrections to the thermodynamic characteristics of graphene is also anomalous. We discuss experimental manifestations of the effects predicted. The microscopic origin of the B-field sensitivity is an extra phase acquired by the electron wave function resulting from the chirality-induced pseudospin precession.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Decoding the 2D IR spectrum of the aqueous proton with high-level VSCF/VCI calculations

The aqueous proton is a common and long-studied species in chemistry, yet there is currently intense interest devoted to understanding its hydration structure and transport dynamics. Typically described in terms of two limiting structures observed in gas-phase clusters, the Zundel H 5 O 2 + and Eigen H 9 O 4 + ions, the aqueous structure is less clear due to the heterogeneity of hydrogen bonding environments and room-temperature structural fluctuations in water. The linear infrared (IR) spectrum, which reports on structural configurations, is challenging to interpret because it appears as a continuum of absorption, and the underlying vibrational modes are strongly anharmonically coupled to each other. Recent two-dimensional IR (2D IR) experiments presented strong evidence for asymmetric Zundel-like motifs in solution, but true structure–spectrum correlations are missing and complicated by the anharmonicity of the system. In this study, we employ high-level vibrational self-consistent field/virtual state configuration interaction calculations to demonstrate that the 2D IR spectrum reports on a broad distribution of geometric configurations of the aqueous proton. We find that the diagonal 2D IR spectrum around 1200 cm –1 is dominated by the proton stretch vibrations of Zundel-like and intermediate geometries, broadened by the heterogeneity of aqueous configurations. There is a wide distribution of multidimensional potential shapes for the proton stretching vibration with varying degrees of potential asymmetry and confinement. In conclusion, we find specific cross peak patterns due to aqueous Zundel-like species. These studies provide clarity on highly debated spectral assignments and stringent spectroscopic benchmarks for future simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Symmetrizing relativistic three-body partial wave amplitudes

S matrix principles and symmetries impose constraints on three-particle scattering amplitudes, which can be formulated as a class of integral equations for their partial wave projections. However, these amplitudes are typically expressed in an asymmetric basis, where one of the initial and final state particles is singled out, and all quantum numbers are defined relative to this spectator. In this work, we show how to construct symmetric partial wave amplitudes, which have been symmetrized over all possible spectator combinations, using their asymmetric counterparts and sets of recoupling coefficients. We derive these recoupling coefficients for arbitrary angular momentum and isospin for arbitrary systems of spinless particles with SU(2) flavor symmetry. We propose a simple intensity observable suitable for visualizing the structure of three-body dynamics in Dalitz distributions. Finally, we provide some numerical examples of Dalitz distributions relevant for future lattice QCD calculations of 3⁢𝜋 systems and provide numerical evidence that the symmetrization procedure is consistent with expected symmetries of Dalitz plots.

Lattice QCD↗
DOE OSTI2023

One-dimensionality signature in optical conductivity of heavy-fermion CeIr 3 B 2

In low dimensions, the combined effects of interactions and quantum fluctuations can lead to dramatically new physics distinct from that existing in higher dimensions. Here, we investigate the electronic and optical properties of CeIr 3 ⁢B 2 , a quasione-dimensional (1D) Kondo lattice system, using ab initio calculations. The Ce atoms in the hexagonal crystal structure form 1D chains along the c axis, with extremely short Ce-Ce distances. The quasi-1D nature of the crystal structure is well reflected in its electronic structure. Extremely flat bands emerge within the ab plane of the Brillouin zone, yielding sharp optical transitions in the corresponding optical conductivity. Further, our calculations indicate that these prominent peaks in the optical conductivity provide a clear signature of quasi-1D heavy fermion systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Electrically tuned hyperfine spectrum in neutral Tb( II )(Cp iPr5 ) 2 single-molecule magnet

Molecular spin qubits with long spin coherence time as well as non-invasive operation methods on such qubits are in high demand. It was shown that both molecular electronic and nuclear spin levels can be used as qubits. In solid state systems with dopants, an electric field was shown to effectively change the spacing between the nuclear spin qubit levels when the electron spin density is high at the nucleus of the dopant. Inspired by such solid-state systems, we propose that divalent lanthanide (Ln) complexes with an unusual electronic configuration of Ln 2+ have a strong interaction between the Ln nuclear spin and the electronic degrees of freedom, which renders electrical tuning of the interaction. Furthermore, as an example, we study electronic structure and hyperfine interaction of the 159 Tb nucleus in a neutral Tb(II)(Cp iPr5 ) 2 single-molecule magnet (SMM), which exhibits unusually long magnetization relaxation time, using the complete active space self-consistent field (CASSCF) method with spin–orbit interaction included within the restricted active space state interaction (RASSI). Our calculations show that the low-energy states arise from 4f 8 (6s,5d z 2 ) 1 , 4f 8 (5d x 2 -y 2 ) 1 , and 4f 8 (5d xy ) 1 configurations. We compute the hyperfine interaction parameters and the electronic–nuclear spectrum within our multiconfigurational approach. We find that the hyperfine interaction is about one order of magnitude greater than that for Tb(III)Pc 2 SMMs. This stems from the strong Fermi contact interaction between the Tb nuclear spin and the electron spin density at the nucleus that originates from the occupation of the (6s,5d) orbitals. We also uncover that the response of the Fermi contact term to electric field results in electrical tuning of the electronic–nuclear level separations. This hyperfine Stark effect may be useful for applications of molecular nuclear spins for quantum computing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Ground state energy and magnetization curve of a frustrated magnetic system from real-time evolution on a digital quantum processor

Models of interacting many-body quantum systems that may realize new exotic phases of matter, notably quantum spin liquids, are challenging to study using even state-of-the-art classical methods such as tensor network simulations. Quantum computing provides a promising route for overcoming these difficulties to find ground states, dynamics, and more. In this paper, we argue that recently developed hybrid quantum-classical algorithms based on real-time evolution are promising methods for solving a particularly important model in the search for spin liquids, the antiferromagnetic Heisenberg model on the two-dimensional kagome lattice. We show how to construct efficient quantum circuits to implement time evolution for the model and to evaluate key observables on the quantum computer, and we argue that the method has favorable scaling with increasing system size. We then restrict to a 12-spin star plaquette from the kagome lattice and a related 8-spin system, and we give an empirical demonstration on these small systems that the hybrid algorithms can efficiently find the ground state energy and the magnetization curve. For these demonstrations, we use four levels of approximation: exact state vectors, exact state vectors with statistical noise from sampling, noisy classical emulators, and (for the 8-spin system only) real quantum hardware, specifically the Quantinuum H1-1 processor; for the noisy simulations and hardware demonstration, we also employ error mitigation strategies based on the symmetries of the Hamiltonian. Our results strongly suggest that these hybrid algorithms present a promising direction for studying quantum spin liquids and more generally for resolving important unsolved problems in condensed matter theory and beyond.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Systematic Improvement of Quantum Monte Carlo Calculations in Transition Metal Oxides: sCI-Driven Wavefunction Optimization for Reliable Band Gap Prediction

Accurate determination of the electronic properties of correlated oxides remains a significant challenge for computational theory. Traditional Hubbard-corrected density functional theory (DFT+U) frequently encounters limitations in precisely capturing electron correlation, particularly in predicting band gaps. We introduce a systematic methodology to enhance the accuracy of diffusion Monte Carlo (DMC) simulations for both ground and excited states, focusing on LiCoO 2 as a case study. By employing a selected configuration interaction (sCI) approach, we demonstrate the capability to optimize wavefunctions beyond the constraints of single-reference DFT+U trial wavefunctions. Here, we show that the sCI framework enables accurate prediction of band gaps in LiCoO 2 , closely aligning with experimental values and substantially improving traditional computational methods. The study uncovers a nuanced mixed state of t 2g and e g orbitals at the band edges that is not captured by conventional single-reference methods, further elucidating the limitations of PBE+U in describing d-d excitations. Our findings advocate for the adoption of beyond-DFT methodologies, such as sCI, to capture the essential physics of excited-state wavefunctions in strongly correlated materials. The improved accuracy in band gap predictions and the ability to generate more reliable trial wavefunctions for DMC calculations underscore the potential of this approach for broader applications in the study of correlated oxides. This work not only provides a pathway for more accurate simulations of electronic structures in complex materials but also suggests a framework for future investigations of the excited states of other challenging systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Quenched disorder in the triangular lattice antiferromagnet YbZn 2⁢ GaO 5

We investigate the crystal electric field (CEF) excitations of Yb 3+ ions in powder samples of the triangular-lattice rare-earth-based antiferromagnet YbZn 2⁢ GaO 5 using inelastic neutron scattering (INS). Three CEF excitations from the ground-state Kramers doublet were observed, each exhibiting significant broadening beyond instrumental resolution. Combining temperature-dependent INS and neutron powder diffraction, we identify a significant static contribution to this broadening and attribute it to heterogeneous coordination of Yb 3+ ions due to Ga 3+ /Zn 2+ site mixing. Rietveld refinement of neutron powder diffraction indicates that 35% of Ga occupies the Zn site and 60% of Zn occupies the Ga site. We show with a point charge model for the CEF Hamiltonian that heterogeneous coordination of Yb3+ ions leads to broadened CEF peaks. First-principles calculations demonstrate that the random Ga 3+ /Zn 2+ distribution can produce the distortions of the YbO 6 octahedra observed from neutron diffraction. Because the documented heterogeneity will extend to exchange interactions, our results suggest that disorder is a significant factor in the unusual magnetism previously reported in YbZn 2 ⁢GaO 5 , including broad low-energy magnetic excitations and the absence of magnetic ordering down to 0.3 K.

Crystal field excitations↗
DOE OSTI2025

Spontaneous generation of athermal phonon bursts within bulk silicon causing excess noise, low energy background events, and quasiparticle poisoning in superconducting sensors

Solid state phonon detectors used in the search for dark matter and coherent neutrino nucleus interactions (CE $v$ NS) require excellent energy resolution (eV-scale or below) and low backgrounds. An unknown source of phonon bursts, the low energy excess (LEE), dominates other above-threshold backgrounds and generates excess shot noise from subthreshold bursts. Here, in this paper, we measure these phonon bursts for 12 days after cooldown in two nearly identical 1 cm 2 silicon detectors that differ only in the thickness of their substrate (1 vs 4 mm thick). We find that both the channel-correlated shot noise and near-threshold shared LEE relax with time since cooldown. Additionally, both the correlated shot noise and LEE rates scale linearly with substrate thickness. When combined with previous measurements of other silicon phonon detectors with different substrate geometries and mechanical support strategies, these measurements strongly suggest that the dominant source of both above and below threshold LEE is the bulk substrate. By monitoring the relation between bias power and excess phonon shot noise, we estimate that the energy scale for subthreshold noise events is 0.68 ± 0.38 meV. In our final dataset, we report a world-leading energy resolution of 258.5 ± 0.4 meV in the 1 mm thick detector. Simple calculations suggest that these silicon substrate phonon bursts are likely a significant source of quasiparticle poisoning in superconducting qubits operated in well shielded and vibration free environments.

Chang, C. L. [Argonne National Laboratory 1 , 9700↗
DOE OSTI2024

Global quark spin correlations in relativistic heavy ion collisions

The observation of the vector meson’s global spin alignment by the STAR Collaboration reveals that strong spin correlations may exist for quarks and antiquarks in relativistic heavy-ion collisions in the normal direction of the reaction plane. We propose a systematic method to describe such correlations in the quark matter. We classify them as local and long range quark spin correlations in the system. We show in particular that the effective quark spin correlations contain the genuine spin correlations originated directly from the dynamical process as well as those induced by averaging over other degrees of freedom. We also show that such correlations can be studied by measuring the vector meson’s spin density matrix and hyperon-hyperon and hyperon-anti-hyperon spin correlations. We present the relationships between these measurable quantities and spin correlations of quarks and antiquarks. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Suppression of antiferromagnetic order and strong ferromagnetic spin fluctuations in Ca(Co 1-x Ni x ) 2-y As 2 single crystals

CaCo 2–y As 2 is a unique itinerant system having strong magnetic frustration. Here, we report the effect of electron doping on the physical properties resulting from Ni substitutions for Co. The single crystals of Ca(Co 1–x Ni x ) 2–y As 2 were characterized by single-crystal x-ray diffraction, energy-dispersive x-ray spectroscopy, magnetization M versus temperature T, magnetic field H, time t, and heat capacity C p (H, T) measurements. The A-type antiferromagnetic (AFM) transition temperature T N = 52 K for x = 0 decreases to 22 K with only 3% Ni substitution and is completely suppressed for x > 0.16. For 0.11 ≤ x ≤ 0.52 strong ferromagnetic (FM) fluctuations develop as revealed by magnetic susceptibility χ(T) = M(T) / H measurements. For x = 0.11 and 0.16 competing AFM and FM interactions result in a reentrant spin-glass behavior below T N , as evidenced by the observations of thermomagnetic hysteresis and magnetic relaxation. Enhanced FM fluctuations are also found for the x = 0.21 and 0.31 crystals, where χ c increases significantly at low T. A large χ anisotropy in these compositions where χ c is up to a factor of two larger than χ ab suggests that the FM spin fluctuations are quasi-1D in nature. Weak ferromagnetic contributions to the magnetization are found at T = 2 K for x = 0.11 –0.31. Heat-capacity Cp(T) measurements reveal the presence of FM quantum spin fluctuations for 0.11 ≤ x ≤ 0.52, where a logarithmic T dependence of Cp(T) / T is observed at low T. Here, the suppression of AFM order by the development of strong FM fluctuations in Ca(Co 1–x Ni x ) 2–y As 2 crystals suggests the presence of a FM quantum-critical point at x ≈ 0.20. Our density-functional theory (DFT) calculations confirm that FM fluctuations are enhanced by Ni substitutions for Co in CaCo 2–y As 2 . The Sommerfeld electronic heat-capacity coefficient is enhanced for x = 0, 0.21, and 0.42 by about a factor of two compared to DFT calculations of the density of states (DOS) at the Fermi energy, suggesting an enhancement of the DOS from electron-phonon and/or electron-electron interactions. The crystals with x > 0.52 do not exhibit FM spin fluctuations or magnetic order at T ≥ 1.8 K, which was found from the DFT calculations to arise from a Stoner transition. Superconductivity is not observed above 1.8 K for any of the compositions. Neutron-diffraction studies of crystals with x = 0.11 and 0.16 in the crossover regime (0.1 ≲ x ≲ 0.2) show no evidence of A-type ordering as observed in the parent compound with x = 0. Furthermore, no other common magnetic structures, such as ferromagnetic (FM), helical stacking of in-plane FM layers, or in-plane AFM structure, are found with an ordered moment greater than the uncertainty of 0.05 μ B per transition-metal atom.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

r 2 SCAN-D4: Dispersion corrected meta-generalized gradient approximation for general chemical applications

We combine a regularized variant of the strongly constrained and appropriately normed semilocal density functional with the latest generation semi-classical London dispersion correction. The resulting density functional approximation r 2 SCAN-D4 has the speed of generalized gradient approximations while approaching the accuracy of hybrid functionals for general chemical applications. We demonstrate its numerical robustness in real-life settings and benchmark molecular geometries, general main group and organo-metallic thermochemistry, and non-covalent interactions in supramolecular complexes and molecular crystals. Main group and transition metal bond lengths have errors of just 0.8%, which is competitive with hybrid functionals for main group molecules and outperforms them for transition metal complexes. The weighted mean absolute deviation (WTMAD2) on the large GMTKN55 database of chemical properties is exceptionally small at 7.5 kcal/mol. This also holds for metal organic reactions with an MAD of 3.3 kcal/mol. The versatile applicability to organic and metal–organic systems transfers to condensed systems, where lattice energies of molecular crystals are within the chemical accuracy (errors <1 kcal/mol).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗