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.

576 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Metal-insulator transition and magnetism of SU(3) fermions in the square lattice

We study the SU(3) symmetric Fermi-Hubbard model (FHM) in the square lattice at 1/3-filling using numerically exact determinant quantum Monte Carlo and numerical linked-cluster expansion techniques. We present the different regimes of the model in the T–U plane, which are characterized by local and short-range correlations, and capture signatures of the metal-insulator transition and magnetic crossovers. These signatures are detected as the temperature scales characterizing the rise of the compressibility, and an interaction-dependent change in the sign of the diagonal spin-spin correlation function. The analysis of the compressibility estimates the location of the metal-insulator quantum critical point at U c /t ~ 6, and provides a temperature scale for observing Mott physics at finite T. Furthermore, from the analysis of the spin-spin correlation function we observe that for U/t ≳ 6 and T ~ J = 4⁢t 2 /U there is a development of a short-range two-sublattice (2SL) antiferromagnetic structure, as well as an emerging three-sublattice (3SL) antiferromagnetic structure as the temperature is lowered below T/J ≲ 0.57. This crossover from 2SL to 3SL magnetic ordering agrees with Heisenberg limit predictions, and has observable effects on the density of on-site pairs. Finally, we describe how the features of the regimes in the T–U plane can be explored with alkaline-earth-like atoms in optical lattices with currently achieved experimental techniques and temperatures. Furthermore, the results discussed in this paper provide a starting point for the exploration of the SU(3) FHM upon doping.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Uncovering extreme nonlinear dynamics in solids through time-domain field analysis

Time-domain analysis of harmonic fields with sub-cycle resolution is now experimentally viable due to the emergence of sensitive, on-chip techniques for petahertz-scale optical-field sampling. We demonstrate how such a time-domain, field-resolved analysis uncovers the extreme nonlinear electron dynamics responsible for high-harmonic generation within solids. Time-dependent density functional theory was used to simulate harmonic generation from a solid-state band-gap system driven by near- to mid-infrared waveforms. Particular attention was paid to regimes where both intraband and interband emission mechanisms play a critical role in shaping the nonlinear response. We show that a time-domain analysis of the harmonic radiation fields identifies the interplay between intra- and interband dynamical processes underlying the nonlinear light generation. With further analysis, we show that changes to the dominant emission regime can occur after only slight changes to the peak driving intensity and central driving wavelength. Furthermore, time-domain analysis of harmonic fields also reveals, for the first time, the possibility of rapid changes in the dominant emission mechanism within the temporal window of the driving pulse envelope. Finally, we examine the experimental viability of performing time-domain analysis of harmonic fields with sub-cycle resolution using realistic parameters.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Time-dependent phenomena in correlated materials

Understanding time-dependent processes and light-matter interaction in strongly correlated materials, and the interplay between electronic, orbital, vibrational, and spin degrees of freedom, is a cornerstone of condensed matter. These mechanisms can be proven by measuring the response of the systems to time-dependent perturbations. The corresponding time scales are dictated by the way light couples to the different excitations, and how these excitations exchange energy and momentum. Our research advances our understanding of these processes, and the interpretation of different equilibrium and time-resolved spectroscopies. Our project encompasses two main themes: (i) developing and refining computational techniques to study non-equilibrium spectroscopies including non-perturbative effects and (ii) applications to non-equilibrium phenomena. We have developed a new computational approach that works directly in the time domain: by including all the degrees of freedom involved in the scattering process (e.g. electrons, photons, neutrons), we solve the time dependent problem: a faithful numerical simulation of the experiment. By measuring the energy and momentum of the outgoing particles, we can extract information about the energy and momentum absorbed by the system. Prior to our work, people attempting to model and calculate non-equilibrium spectral functions relied on a description of the scattering cross section based on a formulation in the frequency domain, a treatment that is extremely cumbersome and complex. Our technique works in and out of equilibrium and can reproduce spectra by several spectroscopic techniques, such as time-resolved photoemission, neutron scattering, Raman, X-ray spectroscopies (RIXS, Auger, XAS, XMCD), and, by not relying on analytical approximations, yields results that reveal novel overlooked transient mechanisms. These tools provide sorely needed intuition for understanding the phenomenology of strongly correlated materials and will help experimentalists in identifying signatures of relevant excitations in pump-probe experiments, such as those conducted in DOE supported facilities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Toward Complete Leading-Order Predictions for Neutrinoless Double β Decay

The amplitude for the neutrinoless double β (0νββ) decay of the two-neutron system nn → ppe – e – constitutes a key building block for nuclear-structure calculations of heavy nuclei employed in large-scale 0νββ searches. Assuming that the 0νββ process is mediated by a light-Majorana-neutrino exchange, a systematic analysis in chiral effective field theory shows that already at leading order a contact operator is required to ensure renormalizability. Here, we develop a method to estimate the numerical value of its coefficient (in analogy to the Cottingham formula for electromagnetic contributions to hadron masses) and validate the result by reproducing the charge-independence-breaking contribution to the nucleon-nucleon scattering lengths. Our central result, while derived in dimensional regularization, is given in terms of the renormalized amplitude Aν(|p|, |p'|) , matching to which will allow one to determine the contact-term contribution in regularization schemes employed in nuclear-structure calculations. Our results thus greatly reduce a crucial uncertainty in the interpretation of searches for 0νββ decay.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Roadmap for warm dense matter physics

This roadmap presents the state-of-the-art, current challenges and near future developments anticipated in the thriving field of warm dense matter (WDM) physics. Originating from strongly coupled plasma physics, high pressure physics and high energy density science, the WDM physics community has recently taken a giant leap forward. This is due to spectacular developments in laser technology, diagnostic capabilities, and computer simulation techniques. Only in the last decade has it become possible to perform accurate enough simulations & experiments to truly verify theoretical results as well as to reliably design experiments based on predictions. Consequently, this roadmap discusses recent developments of and contemporary challenges for theoretical methods and experimental techniques needed to describe, create and diagnose WDM. A large part of this roadmap is dedicated to specific WDM systems and applications in astrophysics, inertial confinement fusion and novel material synthesis.

dense astrophysical objects↗
DOE OSTI2025

Ambiguities in the partial-wave analysis of the photoproduction of pairs of pseudoscalar mesons

Applying the technique of partial-wave analysis, there are cases where more than one set of underlying complex-valued amplitudes can describe the measured observables. These ambiguities can sometimes be resolved using additional information, but assumptions are often required. It is known that the partial-wave analysis of two-pseudoscalar meson systems produced in photoproduction with a linearly polarized photon beam is free from discrete ambiguities stemming from the Barrelet zeros when the nucleon spin is ignored. In this article, we show that continuous ambiguities are possible for certain wave sets, even though the discrete ambiguities do not appear. We also explore ways to resolve these ambiguities and determine the maximal amount of information that can be obtained from analyses that suffer from these continuous ambiguities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Collective Phonon–Polaritonic Modes in Silicon Carbide Subarrays

Localized surface phonon polaritons (LSPhPs) can be implemented to engineer light–matter interactions through nanoscale patterning for a range of midinfrared application spaces. However, the polar material systems studied to date have mainly focused on simple designs featuring a single element in the periodic unit cell. Increasing the complexity of the unit cell can serve to modify the resonant near-fields and intra- and inter-unit-cell coupling as well as to dictate spectral tuning in the far-field. In this work, we exploit more complicated unit-cell structures to realize LSPhP modes with additional degrees of design freedom, which are largely unexplored. Collectively excited LSPhP modes with distinctly symmetric and antisymmetric near-fields are supported in these subarray designs, which are based on nanopillars that are scaled by the number of subarray elements to ensure a constant unit-cell size. Moreover, we observe an anomalous mode-matching of the collective symmetric mode in our fabricated subarrays that is robust to changing numbers of pillars within the subarrays as well as to defects intentionally introduced in the form of missing pillars. This work therefore illustrates the hierarchical design of tailored LSPhP resonances and modal near-field profiles simultaneously for a variety of IR applications such as surface-enhanced spectroscopies and biochemical sensing.

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

Gauging nexus between topological and fracton phases

Coupled layer constructions are a valuable tool for capturing the universal properties of certain interacting quantum phases of matter in terms of the simpler data that characterizes the underlying layers. In the study of fracton phases, the X-Cube model in 3+1D can be realized via such a construction by starting with a stack of 2+1D Toric Codes and turning on a coupling which condenses a composite "particle-string" object. In a recent work [Phys. Rev. B 112, 125124 (2025)], we have demonstrated that in fact, the particle-string can be viewed as a symmetry defect of a topological 1-form symmetry. In this paper, we study the result of gauging this symmetry in depth. We unveil a rich gauging web relating the X-Cube model to symmetry protected topological (SPT) phases protected by a mix of subsystem and higher-form symmetries, subsystem symmetry fractionalization in the 3+1D Toric Code, and non-trivial extensions of topological symmetries by subsystem symmetries. Here, our work emphasizes the importance of topological symmetries in non-topological, geometric phases of matter.

Anyons↗
DOE OSTI2020

Top Loading Helium Cryostat Integrated with High-Pressure Cell with Fast Remote Pressure Control (Final Report for SBIR Phase 1)

A large number of physical phenomena, such as superconductivity and quantum critical phenomena, often appear only at very low temperatures below 5 K. There is an immense interest to investigate these phenomena at high pressure as a means of tuning interatomic distances, and thus the interaction parameters controlling these phenomena, in a continuous and controlled fashion. The current P-T condition for neutron scattering experiments are limited to either relatively low pressures of about 2 GPa at temperatures below 5 K, or to relatively high temperatures at pressures or tens of GPa. The purpose of this research is to develop neutron sample environment instrumentation for reaching 10-20 GPa at 1-2K with rapid and reliable online pressure and temperature control (i.e. without having to interrupt the experiment). The ultimate goal of the project is to design an integrated fast-cooling low-temperature sample environment cryogenic system compatible with state of the art neutron diamond anvil cells suitable for single-crystal neutron scattering experiments for temperatures down to 2K and pressures of several tens of GPa. The integrated system will consist of top-loading Helium flow cryostat with in-situ sample alignment mechanisms, large-volume diamond anvil cells (DAC) made from novel superalloy Pascalloy and optimized for fast cooling and heating, and a compact remote pressure control mechanism for the DAC based on a novel concept of inflatable bellows integrated with a lever-arm based force amplifier. In Phase I, we have designed, manufactured, and tested prototypes of the novel compact force-amplified pneumatic pressure control mechanism for Neutron Diamond Anvil Cell (nDAC), which allows to use pneumatic bellows system for smooth remote pressure control in the nDAC inside a top-loading cryostat with bore size of 70 mm or larger. This allows significant minimum temperature decrease in remotely controlled nDAC from 5-10 K down to 2 K. We also studied mechanical properties of a novel non-magnetic superalloy Pascalloy with different heat-treatment conditions and preliminary results indicate that up to date this is probably the strongest and most suitable material for making more compact and lightweight cryogenic nDACs for neuron scattering experiments at extreme conditions, allowing much faster cooling and heating. Potential applications: The new developments will allow to create a wide range of compact Diamond Anvil Cells for neutron diffraction which can reach several tens of GPA pressure at 2-4 K while preserving accurate remote pressure control capabilities. The new concepts can be used for developing sample environment instrumentation outside the neutron scattering field. The new development have very strong potential for expanding experimental capabilities in materials sciences and have high commercialization potential. The proposed new cryogenic high-pressure system or any of its individual components will be in high demand not only in neutron scattering facilities, but also at synchrotron beamlines and other high-pressure research facilities around the world.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Polariton relaxation under vibrational strong coupling: Comparing cavity molecular dynamics simulations against Fermi’s golden rule rate

Under vibrational strong coupling (VSC), the formation of molecular polaritons may significantly modify the photo-induced or thermal properties of molecules. In an effort to understand these intriguing modifications, both experimental and theoretical studies have focused on the ultrafast dynamics of vibrational polaritons. Here, following our recent work [Li et al., J. Chem. Phys. 154, 094124 (2021)], we systematically study the mechanism of polariton relaxation for liquid CO 2 under a weak external pumping. Classical cavity molecular dynamics (CavMD) simulations confirm that polariton relaxation results from the combined effects of (i) cavity loss through the photonic component and (ii) dephasing of the bright-mode component to vibrational dark modes as mediated by intermolecular interactions. The latter polaritonic dephasing rate is proportional to the product of the weight of the bright mode in the polariton wave function and the spectral overlap between the polariton and dark modes. Both these factors are sensitive to parameters such as the Rabi splitting and cavity mode detuning. Compared to a Fermi’s golden rule calculation based on a tight-binding harmonic model, CavMD yields a similar parameter dependence for the upper polariton relaxation lifetime but sometimes a modest disagreement for the lower polariton. We suggest that this disagreement results from polariton-enhanced molecular nonlinear absorption due to molecular anharmonicity, which is not included in our analytical model. Finally, we also summarize recent progress on probing nonreactive VSC dynamics with CavMD.

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

Beyond fragmented dopant dynamics in quantum spin lattices: Robust localization and non-Gaussian diffusion

The motion of dopants in magnetic spin lattices has received tremendous attention for at least four decades due to its connection to high-temperature superconductivity. Despite these efforts, we lack a complete understanding of their behavior, especially out of the equilibrium and at nonzero temperatures. In this paper, we take a significant step towards a much deeper understanding based on state-of-the-art matrix-product-state calculations. In particular, we investigate the nonequilibrium dynamics of a dopant in two-leg 𝑡−𝐽 ladders with antiferromagnetic XXZ spin interactions. In the Ising limit, we find that the dopant is localized for all investigated nonzero temperatures due to an emergent disordered potential, with a localization length controlled by the underlying correlation length of the spin lattice, which increases exponentially with decreasing temperature. The dopant, hereby, only delocalizes asymptotically in the zero temperature limit. This greatly generalizes the localization effect discovered recently in Hilbert space fragmented models [Phys. Rev. Res. 6, 023325 (2024), SciPost Phys. Core 7, 054 (2024)]. In the presence of spin-exchange processes at rate 𝛼, the dopant diffuses with a diffusion coefficient, 𝐷 ℎ , depending nonmonotonically on 𝛼. It initially increases linearly as 𝐷 ℎ ∝ 𝛼 for 𝛼 ≪ 1 before dropping off as 𝛼 −1 for 𝛼 > 1. Moreover, we show that the underlying spin dynamics at infinite temperature behaves qualitatively the same, albeit with important quantitative differences. We substantiate these findings by showing that the dynamics features self-similar scaling behavior, which strongly deviates from the Gaussian behavior of regular diffusion, especially for weak spin exchange. Finally, we show that the diffusion coefficient 𝐷 ℎ follows an Arrhenius relation at high temperatures, whereby it is exponentially suppressed upon cooling.

Anomalous diffusion↗
DOE OSTI2023

Free energy and specific heat near a quantum critical point of a metal

We analyze the free energy and the specific heat for fermions interacting with a gapless boson at a quantum-critical point (QCP) in a metal. We use the Luttinger-Ward-Eliashberg formula for the free energy in the normal state, which includes contributions from bosons, fermions, and their interaction, all expressed via fully dressed fermionic and bosonic propagators. The sum of the last two contributions is the free energy F γ of an effective low-energy model of fermions with boson-mediated dynamical 4-fermion interaction V(Ω m ) ∝ 1/|Ω m | γ (the γ model). This purely electronic model has been used to analyze the interplay between non-Fermi liquid behavior and pairing near a QCP, which are both independent of the upper energy cutoff Λ. However, the specific heat C γ (T), obtained from F γ , does depend on Λ. We argue that this dependence is spurious and cancels out, once we include the contribution from bosons. We further argue that the full C(T) is the sum of the contribution from free fermions and the one from a critical boson, with the fully dressed propagator, other terms cancel out. We compare the full C(T) with the C γ (T), obtained using recently proposed regularization of F γ . We argue that for γ < 1, the full C(T) and the regularized C γ (T) differ by a γ-dependent prefactor, while for γ > 1, the full C(T) and C γ (T) differ by the positive contribution from free massless fermions (a positive constant for the electron-phonon case γ = 2). For these γ, C γ (T) is negative, but the full C(T) is positive. We argue that only the full C(T) matters as the positive and the negative contributions originate from the term in C(T) which contains the fully dressed bosonic propagator. Here, we then argue that the normal state remains stable until the pairing instability develops.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Interaction-enhanced nesting in spin-fermion and Fermi-Hubbard models

The spin-fermion (SF) model postulates that the dominant coupling between low-energy fermions in near critical metals is mediated by collective spin fluctuations (paramagnons) peaked at the Néel wave vector, Q N , connecting hot spots on opposite sides of the Fermi surface. It has been argued that strong correlations at hot spots lead to a Fermi surface deformation (FSD) featuring flat regions and increased nesting. This conjecture was confirmed in the perturbative self-consistent calculations when the paramagnon propagator dependence on momentum deviation from Q N is given by χ − 1 ∝ | Δ q | . Using diagrammatic Monte Carlo (diagMC) technique we show that such a dependence holds only at temperatures orders of magnitude smaller than any other energy scale in the problem, indicating that a different mechanism may be at play. Instead, we find that a χ − 1 ∝ | Δ q | 2 dependence yields a robust finite- T scenario for achieving FSD. To link phenomenological and microscopic descriptions, we applied the connected determinant diagMC method to the ( t − t ′ ) Hubbard model and found that at large U / t > 5.5 before the formation of electron and hole pockets (i) the FSD defined as a maximum of the spectral function is not very pronounced; instead, it is the lines of zeros of the renormalized dispersion relation that deforms toward nesting, and (ii) the static spin susceptibility is well described by χ − 1 ∝ | Δ q | 2 . Flat FS regions yield a nontrivial scenario for realizing a non-Fermi liquid state. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Single-Crystal Diffuse Neutron Scattering Study of the Dipole-Octupole Quantum Spin-Ice Candidate Ce 2⁢ Zr 2 ⁢O 7 : No Apparent Octupolar Correlations Above 𝑇 = 0.05 K

The insulating magnetic pyrochlore Ce 2 ⁢Zr 2 ⁢O 7 has gained attention as a quantum spin-ice candidate with dipole-octupole character that arises from the crystal-electric-field ground-state doublet for the Ce 3+ Kramers ion. This dipole-octupole character permits both spin-ice phases based on magnetic dipoles and those based on more-exotic octupoles. This work reports low-temperature neutron diffraction measurements on single-crystal Ce 2 ⁢Zr 2⁢ O 7 with 𝑄 coverage both at low 𝑄, where the magnetic form factor for dipoles is near maximal, and at high 𝑄, covering the region where the magnetic form factor for Ce 3+ octupoles is near maximal. This study was motivated by recent powder neutron diffraction studies of other Ce-based dipole-octupole pyrochlores, Ce 2 ⁢Sn 2 ⁢O 7 and Ce 2 ⁢Hf 2 ⁢O 7 , which each showed temperature-dependent diffuse diffraction at high 𝑄, interpreted as arising from octupolar correlations. Our measurements use an optimized single-crystal diffuse scattering instrument that allows us to screen against strong Bragg scattering from Ce 2 ⁢Zr 2 ⁢O 7 . The temperature-difference neutron diffraction reveals a low-𝑄 peak consistent with dipolar spin-ice correlations reported in previous work, and an alternation between positive and negative net intensity at higher 𝑄. These features are consistent with our numerical-linked-cluster calculations using pseudospin interaction parameters previously reported for Ce 2 ⁢Zr 2⁢ O 7 , Ce 2 ⁢Sn 2 ⁢O 7 , and Ce 2 ⁢Hf 2 ⁢O 7 . Importantly, neither the measured data nor any of the NLC calculations show evidence for increased scattering at high 𝑄 resulting from octupolar correlations. We conclude that at the lowest attainable temperature for our measurements (𝑇 = 0.05 K), scattering from octupolar correlations in Ce 2 ⁢Zr 2 ⁢O 7 is not present in the neutron diffraction signal on the level of our observation threshold of around 0.1% of the low-𝑄 dipole scattering. We compare these results to those obtained earlier on powder Ce 2 ⁢Sn 2 ⁢O 7 and Ce 2⁢ Hf 2⁢ O 7 , and to low-energy inelastic neutron scattering from single-crystal Ce 2 ⁢Zr 2 ⁢O 7 .

36 MATERIALS SCIENCE↗
DOE OSTI2024

Symmetry Breaking and Ascending in the Magnetic Kagome Metal FeGe

Spontaneous symmetry breaking—the phenomenon in which an infinitesimal perturbation can cause the system to break the underlying symmetry—is a cornerstone concept in the understanding of interacting solid-state systems. In a typical series of temperature-driven phase transitions, higher-temperature phases are more symmetric due to the stabilizing effect of entropy that becomes dominant as the temperature is increased. However, the opposite is rare but possible when there are multiple degrees of freedom in the system. Here, we present such an example of a symmetry-ascending phenomenon upon cooling in a magnetic kagome metal FeGe by utilizing neutron Larmor diffraction and Raman spectroscopy. FeGe has a kagome lattice structure with simple A-type antiferromagnetic order below Néel temperature T N ≈ 400 K and a charge density wave (CDW) transition at T CDW ≈ 110 K , followed by a spin-canting transition at around 60 K. In the paramagnetic state at 460 K, we confirm that the crystal structure is indeed a hexagonal kagome lattice. On cooling to around T N , the crystal structure changes from hexagonal to monoclinic with in-plane lattice distortions on the order of 10 − 4 and the associated splitting of the double-degenerate phonon mode of the pristine kagome lattice. Upon further cooling to T CDW , the kagome lattice shows a small negative thermal expansion, and the crystal structure gradually becomes more symmetric upon further cooling. A tendency of increasing the crystalline symmetry upon cooling is unusual; it originates from an extremely weak structural instability that coexists and competes with the CDW and magnetic orders. These observations are against the expectations for a simple model with a single order parameter and hence can only be explained by a Landau free energy expansion that takes into account multiple lattice, charge, and spin degrees of freedom. Thus, the determination of the crystalline lattice symmetry as well as the unusual spin-lattice coupling is a first step towards understanding the rich electronic and magnetic properties of the system, and it sheds new light on intertwined orders where the lattice degree of freedom is no longer dominant. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Hund’s flat band in a frustrated spinel oxide

Electronic flat bands associated with quenched kinetic energy and heavy electron mass have attracted great interest for promoting strong electronic correlations and emergent phenomena such as high-temperature charge fractionalization and superconductivity. Intense experimental and theoretical research has been devoted to establishing the rich nontrivial metallic and heavy fermion phases intertwined with such localized electronic states. Here, in this study, we investigate the transition metal oxide spinel LiV 2 O 4 , an enigmatic heavy fermion compound lacking localized f orbital states. We use angle-resolved photoemission spectroscopy and dynamical mean-field theory to reveal a kind of correlation-induced flat band with suppressed interatomic electron hopping arising from intra-atomic Hund’s coupling. The appearance of heavy quasiparticles is ascribed to a proximate orbital-selective Mott state characterized by fluctuating local moments as evidenced by complementary magnetotransport measurements. The spectroscopic fingerprints of long-lived quasiparticles and their disappearance with increasing temperature further support the emergence of a high-temperature “bad” metal state observed in transport data. This work resolves a long-standing puzzle on the origin of heavy fermion behavior and unconventional transport in LiV 2 O 4 . Simultaneously, it opens a path to achieving flat bands through electronic interactions in d-orbital systems with geometrical frustration, potentially enabling the realization of exotic phases of matter such as the fractionalized Fermi liquids.

36 MATERIALS SCIENCE↗
DOE OSTI2025

COLUMBUS─An Efficient and General Program Package for Ground and Excited State Computations Including Spin–Orbit Couplings and Dynamics

The COLUMBUS program system provides the tools for performing high-level multireference (MR) computations, including the multireference configuration interaction (MRCI) method and its multireference averaged quadratic coupled cluster (MR-AQCC) extension, allowing computations on a wide range of fascinating atomic and molecular systems, including the treatment of open-shells and complicated excited state phenomena. The inclusion of spin−orbit coupling (SOC) directly within the MRCI step enables the description of systems containing heavy elements, such as lanthanides and actinides, whose properties are strongly influenced by SOC. Analytic energy gradients and nonadiabatic couplings at the correlated MRCI level provide the foundation for a variety of dynamics studies, giving insight into ultrafast photochemistry. New and ongoing method developments in COLUMBUS include the computation of spin densities, improved descriptions of ionic states, enhancements to the AQCC method, and the porting of COLUMBUS to graphical processing units (GPUs). New external interfaces enable an enhanced description of electronic resonances and molecules in strong laser fields. This work highlights these new developments while providing a detailed account of the diverse applications of COLUMBUS in recent years.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗