SEARCH · Engineering Papers

Results for “Linear response 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Correlation effects on coupled electronic and structural properties of doped rare-earth trihydrides

Rare-earth trihydride (𝑅⁢H 3 ) compounds exhibit intriguing coupled electronic and structural properties as a function of doping, hydrogen vacancies, and thermodynamic conditions. Theoretical studies of these materials typically rely on density functional theory (DFT), including the use of small supercells that may underestimate strong correlation effects and structural distortions which in turn may influence their metallicity. Here, we elucidate the roles of lattice distortions and correlation effects on the electronic properties of pristine and doped 𝑅⁢H 3 compounds by adopting DFT +U and quantum Monte Carlo (QMC) methods. Linear-response constrained DFT (LR-cDFT) methods find Hubbard U ≈ 2 eV for 𝑅 𝑑 orbitals and U ≈ 6 eV for H 𝑠 ⁡/N 𝑝 /O 𝑝 orbitals. The small U on Lu 𝑑 orbitals is consistent with QMC calculations on LuH 3 and LuH 2.875 ⁢N 0.125 . In pure face-centered-cubic (FCC) 𝑅⁢H 3 (𝑅 = Lu,Y) compounds, neither DFT nor DFT +U with the self-consistently determined U is enough to create a band gap, however a supercell with hydrogen distortions creates a small gap whose magnitude increases when performing DFT +U with self-consistently determined U values. Furthermore, correlation effects, in turn, have a moderate influence on the coupled structural and electronic properties of doped 𝑅⁢H 3 compounds and may be important when considering the competition between structural distortions and superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Predictable third harmonic generation in GaAs metasurfaces through group theory inverse design of meta-atoms

Here we report the group theory-based inverse design of meta-atoms for a dielectric metasurface in GaAs with predictable optical linear response and third harmonic generation (THG). Six sharp Fano resonances have been observed with a corresponding polarization dependence as predicted by group theory and the meta-atom’s symmetry in the D2h point group. THG has been observed for two modes under x-polarization excitation and one mode for y-polarization, in agreement with theoretical symmetry predictions. The polarization-dependent THG aspect ratio was observed to reach as high as 108. Through strategic structural or symmetry-preserving perturbations, it was shown that the THG can be enhanced or reduced by a factor of 7. The highest THG conversion efficiency was estimated to be 3.1×10-7 at the pump intensity of 1.51 MW/cm2. This high THG conversion efficiency indicates that our group theory approach to modal engineering opens a new path towards optical nonlinearity tuning in dielectric metasurfaces.

Aoueille, Andrew↗

DOE OSTI2025

On the bulk compaction of brittle granular materials, Part II: Representing interparticle comminution by K COMP

An applied and empirical compaction fracture resistance parameter, K COMP , was developed to represent the compression-induced interparticle-comminution response of granular brittle materials (e.g., ceramic or glass particles). The development of K COMP is an outcome from Part I of this three-paper series. The K COMP represents a continuum response of compaction fracture resistance and was adapted from established Griffith linear elastic fracture mechanics theory. From that, the K COMP relates macroscopically applied compressive stress to a corresponding inverse-square-root representative particle size, where the latter is estimated from an entire particle size distribution or specific surface area measurement. In conclusion, the K COMP value was not constant over a wide range of stresses for all examined brittle granular materials, and this is indicative of a change in the dominant mode of permanent deformation at higher compaction stresses.

brittle↗

DOE OSTI2025

Spin qubit properties of the boron-vacancy/carbon defect in the two-dimensional hexagonal boron nitride

Spin qubit defects in two-dimensional materials have a number of advantages over those in three-dimensional hosts including simpler technologies for defect creation and control, as well as qubit accessibility. In this work, we select the V B C B defect in the hexagonal boron nitride (hBN) as a possible optically controllable spin qubit and explain its triplet ground state and neutrality. In this defect a boron vacancy is combined with a carbon dopant substituting the closest boron atom to the vacancy. Our density-functional-theory calculations confirmed that the system has dynamically stable spin triplet and singlet ground states. As revealed from our linear response GW calculations, the spin-sensitive electronic states are localized around the three undercoordinated N atoms and make local peaks in the density of electronic states within the bandgap. Using the triplet and singlet ground state energies, as well as the energies of the optically excited states, obtained from solution to the Bethe–Salpeter equation, we construct the spin-polarization cycle, which is found to be favorable for the spin qubit initialization. The calculated zero-field splitting parameters ensure that the splitting energy between the spin projections in the triplet ground state is comparable to that of the known spin qubits. We thus propose the V B C B defect in hBN as a promising spin qubit.

2D BN↗

DOE OSTI2025

Colossal Cryogenic Electro‐Optic Response Through Metastability in Strained BaTiO 3 Thin Films

The search for thin film electro-optic materials that can retain superior performance under cryogenic conditions has become critical for quantum computing. Barium titanate thin films show large linear electro-optic coefficients in the tetragonal phase at room temperature, which is severely degraded down to ≈200 pm V −1 in the rhombohedral phase at cryogenic temperatures. There is immense interest in manipulating these phase transformations and retaining superior electro-optic properties down to liquid helium temperature. Utilizing the thermodynamic theory of optical properties, a large low-temperature electro-optic response is designed by engineering the energetic competition between different ferroelectric phases, leading to a low-symmetry monoclinic phase with a massive electro-optic response. The existence of this phase is demonstrated in a strain-tuned BaTiO 3 thin film that exhibits a linear electro-optic coefficient of 2516 ± 100 pm V −1 at 5 K, which is an order of magnitude higher than the best reported performance thus far. Importantly, the electro-optic coefficient increases by 100 × during cooling, unlike the conventional films, where it degrades. Further, at the lowest temperature, significant higher order electro-optic responses also emerge. These results represent a new framework for designing materials with property enhancements by stabilizing highly tunable metastable phases with strain.

barium titanate↗

DOE OSTI2025

Galactic echoes

Gaia has revealed a variety of substructures in the phase space of stars in the Solar neighbourhood, including the vertical ‘Snail’ $(z, v_z)$ in space. Such substructures are often interpreted as the incompletely phase-mixed response of the disc stars to a single perturbation, such as an impulsive encounter with a satellite galaxy. In this paper, we consider the possibility that such structures contain manifestations of phase-space echoes. First established in plasma physics in the 1960s, echoes arise when a collisionless system is perturbed twice: the macroscopic responses to both perturbations mix to small scales in phase space, whereupon they couple non-linearly, producing a third macroscopic ‘echo’ response without the need for a third perturbation. We derive the galactic analogue of the plasma echo theory using angle-action variables and apply it to a one-dimensional model of vertical motion in the Milky Way. We verify the predicted echo behaviour using idealized test particle simulations, both with and without the inclusion of diffusion through orbital scattering off molecular clouds. While we conclude that the Gaia Snail itself is unlikely a (pure) echo effect, the basic physics we uncover is sufficiently generic that we expect phase-space echoes to be common in disc galaxies.

Galaxy: evolution↗

DOE OSTI2025

Nonlinear Topological Photonics: Capturing Nonlinear Dynamics and Optical Thermodynamics

Combining multiple optical resonators or engineering dispersion of complex media has provided an effective method for demonstrating topological physics controlling photons in unprecedented ways such as unidirectional light propagation and spatially localized modes between an interface or on a corner. Further, adding nonlinear responses to those topological photonic systems has enabled achieving diverse phases of photons in both space and time, allowing for more functionalities in photonic devices that provide a new playground for studying dynamic features of nonlinear topological systems. However, most methods for describing nonlinear topological photonic systems rely on linear topological theories, making it challenging to accurately characterize the topology of nonlinear systems. Thus, substantial efforts have focused on rigorously describing nonlinear topological phases and developing effective tools to analyze nonlinear topological effects. Meanwhile, coupled multimode optical waveguides with nonlinear dynamic responses provide an excellent platform for the statistical description of photons, opening a new paradigm called “optical thermodynamics”. This review will introduce the basic concepts of nonlinear topological photonics and the recent development of theoretical approaches focusing on data-driven approaches for creating phase diagrams as well as the spectral localizer framework and the pseudospectrum method for understanding optical nonlinearities in topological systems. In addition, the new concept of optical thermodynamics will be introduced with some recent theoretical works.

Topological photonics↗

DOE OSTI2025

Continuum model of strong light-matter coupling for molecular polaritons

Strong coupling between light and matter generates hybrid polariton modes. We present a continuum formalism that expresses the polariton modes in terms of light and matter densities of states (DOS). We derive exact expressions for the light and matter DOS for a planar cavity containing a strongly dispersive dielectric.We show that these DOS depend exclusively on the linear response of the cavity components, i.e., the reflectance of the mirrors and the susceptibility of the dielectric. Here, we further show that, within the strong coupling regime, the light and matter DOS are well approximated by coupled mode theory. Altogether, our continuum formalism offers a unified treatment of polaritons that connects the framework of light-matter coupling with that of linear dispersion.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Non-resonant Raman optical activity from phase-space electronic structure theory

In order to model experimental non-resonant Raman optical activity, chemists must compute a host of second-order response tensors (e.g., the electric-dipole–magnetic-dipole polarizability) and their nuclear derivatives along a set of vibrational modes. While these response functions are almost always computed within a Born–Oppenheimer (BO) framework, here we provide a natural interpretation of the electric-dipole–magnetic-dipole polarizability within phase space electronic structure theory, a beyond-BO model whereby the electronic structure depends on nuclear momentum (P) in addition to nuclear position (R). By coupling to nuclear momentum, phase space electronic structure theory is able to capture the asymmetric response of the electronic properties to an external field, in so far as for a vibrating (non-stationary) molecule, $\frac{∂μ}{∂B}$≠$\frac{∂m}{∂F}$, where μ and m are the electrical linear and magnetic dipoles, and F and B are electric and magnetic fields. As an example, for a prototypical methyloxirane molecule, we show that phase space electronic structure theory is able to deliver a reasonably good match with experimental results in a manner that is formally invariant to gauge origin G 0 —provided that one uses a complete basis or, alternatively, gauge invariant atomic orbitals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Experimental and computational study of phase space dynamics in strongly coupled plasmas with steep density gradients

Understanding how plasmas thermalize when density gradients are steep remains a fundamental challenge in plasma physics, with direct implications for fusion experiments and astrophysical phenomena. Standard hydrodynamic models break down in these regimes, and kinetic theories make predictions that have never been directly tested. Here, we present the first detailed phase-space measurements of a strongly coupled plasma as it evolves from sharp density gradients to thermal equilibrium. Using laser-induced fluorescence imaging of an ultracold calcium plasma, we track the complete ion distribution function f(x,v,t). We discover that commonly used kinetic models (Bhatnagar–Gross–Krook and Lenard–Bernstein) overpredict thermalization rates, even while correctly capturing the initial counterstreaming plasma formation. Our measurements reveal that the initial ion acceleration response scales linearly with electron temperature, and that the simulations underpredict the initial ion response. In our geometry we demonstrate the formation of well-controlled counterpropagating plasma beams. This experimental platform enables precision tests of kinetic theories and opens new possibilities for studying plasma stopping power and flow-induced instabilities in strongly coupled systems.

Bergeson, Scott (ORCID:0000000231249226)↗

DOE OSTI2024

Prediction of photodynamics of 200 nm excited cyclobutanone with linear response electronic structure and ab initio multiple spawning

Simulations of photochemical reaction dynamics have been a challenge to the theoretical chemistry community for some time. In an effort to determine the predictive character of current approaches, we predict the results of an upcoming ultrafast diffraction experiment on the photodynamics of cyclobutanone after excitation to the lowest lying Rydberg state (S 2 ). A picosecond of nonadiabatic dynamics is described with ab initio multiple spawning. Herein we use both time dependent density functional theory (TDDFT) and equation-of-motion coupled cluster singles and doubles (EOM-CCSD) theory for the underlying electronic structure theory. We find that the lifetime of the S 2 state is more than a picosecond (with both TDDFT and EOM-CCSD). The predicted ultrafast electron diffraction spectrum exhibits numerous structural features, but weak time dependence over the course of the simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Random Phase Approximation Correlation Energy Using Real-Space Density Functional Perturbation Theory

We present a real-space method for computing the random phase approximation (RPA) correlation energy within Kohn–Sham density functional theory, leveraging the low-rank nature of the frequency-dependent density response operator. In particular, we employ a cubic-scaling formalism based on density functional perturbation theory that circumvents the calculation of the response function matrix, instead relying on the ability to compute its product with a vector through the solution of the associated Sternheimer linear systems. We develop a large-scale parallel implementation of this formalism using the subspace iteration method in conjunction with the spectral quadrature method while employing the Kronecker product-based method for the application of the Coulomb operator and the conjugate orthogonal conjugate gradient method for the solution of the linear systems. We demonstrate convergence with respect to key parameters and verify the method’s accuracy by comparing with plane-wave results. We show that the framework achieves good strong scaling to many thousands of processors, reducing the time to solution for a lithium hydride system with 128 electrons to around 150 s on 4608 processors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Simplified spin dependence in dark matter direct detection

The interactions of dark matter with Standard Model particles can be systematically studied in the language of effective field theories. We investigate dark matter interactions with Standard Model particles, including spin-dependent interactions, for direct detection experiments and demonstrate that, although the scattering rate generally depends on multiple types of material response functions, certain linear combinations of these material response functions vanish if the initial and final electronic states share the same Hamiltonian. We also find that several other response functions vanish in parity-symmetric materials, making these systems as simple as isotropic detectors in some respects. Finally, we present the scattering rate for an anisotropic, possibly chiral detector, for generic dark matter-electron spin interactions. These relations reduce the number of independent response functions needed, thereby simplifying the computational complexity for a broad class of dark matter models. Our results provide a complete and efficient toolkit for analyzing electron recoil signals in diverse detector materials.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Looking at extremal black holes from very far away

Near-extremal black holes are subject to large quantum effects, which modify their low-temperature thermodynamic behavior. Hitherto, these quantum effects were analyzed by separating the geometry into the near-horizon region and its exterior. It is desirable to understand and reproduce such corrections from the full higher-dimensional asymptotically flat or AdS geometry’s perspective. We address this question in this article and fill this gap. Specifically, we find off-shell eigenmodes of the quadratic fluctuation operator of the Euclidean gravitational dynamics, with eigenvalues that vanish linearly with temperature. We illustrate this for BTZ and neutral black holes with hyperbolic horizons in AdS in Einstein-Hilbert theory, and for the charged black holes in Einstein-Maxwell theory. The linear scaling with Matsubara frequency, which is a distinctive feature of the modes, together with the fact that their wavefunctions localize close to the horizon as we approach extremality, identifies them as responsible for the aforementioned quantum effects. We provide a contour prescription to deal with the sign indefiniteness of the Euclidean Einstein-Maxwell action, which we derive to aid our analysis. We also resolve a technical puzzle regarding modes associated with rotational isometries in stationary black hole spacetimes.

AdS-CFT Correspondence↗

DOE OSTI2024

The nonlinear behavior of generic tail fins for small wind turbines

This paper describes analysis and measurements of the yaw response of tail fins for small wind turbines. It is based on an extension of unsteady slender body theory (USBT) to cover non-slender fins and high angles of incidence, both of which make the theory nonlinear. We provide three main additions to the substantial literature on linearized USBT for tail fins. First, USBT is extended to high angles by modeling the nonlinear vortex dynamics. Second, the restriction to slender bodies is removed by modeling the chordwise load variation. Third, we consider the effect of time-varying wind speed. Further, the extended theory is compared to wind tunnel measurements of the yaw behavior of delta, elliptical, and rectangular tail fins without a rotor and nacelle. The fins were released from initial yaw angles of -40° and -80°; the latter is of sufficient magnitude to show the importance of the nonlinear yaw dynamics. Generally good agreement was found between the theory and measurements, and the theory was shown to be more accurate than a “polar” or quasi-steady model which uses only the lift and drag of a delta planform. Of the three planforms, the rectangular one showed the lowest accuracy in terms of frequency but the damping was accurately predicted. Overall, the results demonstrate the importance of nonlinearity in the response of a yawing tail fin, particularly for the higher aspect ratio fins at large yaw angles.

17 WIND ENERGY↗

DOE OSTI2025

Predicted thresholds for RMP ELM suppression access in double-null configurations

Modeling of DIII-D plasmas spanning shapes from single to double-null (DN) reveals new insights into the nature of resonant magnetic perturbation (RMP) conditions necessary for edge-localized mode (ELM) suppression. The suppression of ELMs with RMPs has proven difficult in DN configurations, where no device has thus far reported any hints of suppression. Modeling using the GPEC code finds a reduced high-field side response closer to DN shaping. The resulting synthetic diagnostic measurements are consistent with what has been observed in experiments on DIII-D, validating the plasma response model in this regime. While common metrics for suppression do not illustrate a clear distinction between single-null (SN) and DN cases, the pedestal top resonant field does show a ∼20% decrease at DN shaping in modeling. Field penetration is assessed using linear tearing theory, which demonstrates a lack of sufficient pedestal top resonant flux in the DN shape, requiring at least 1.5 × greater RMP coil currents than what was used in experiment. Analysis from drift kinetic simulations further indicate up to 2 × larger critical island widths are required at the pedestal top for tearing mode growth compared to SN cases. Effective island widths inferred from 3D ideal MHD are also analyzed, where maximum widths in lab coordinates indicate a threshold of ∼18–24× the ion gyroradius for sufficient profile flattening for ELM suppression. These results suggest that ELM suppression may be possible in DN with sufficiently large RMP coil amplitude. Future prospects of achieving RMP ELM suppression in the DN configuration may involve going to lower triangularity, which is also highlighted in this work.

RMP ELM suppression↗

DOE OSTI2025

Hydrodynamic characterization of the coastal pioneer array ocean observing system

Ocean observation buoys require relatively small amounts of power, yet traditionally necessitate costly resupply trips for battery replacement. With the offshore location of the buoys and small power requirements, wave energy may be an effective solution for providing consistent and reliable power to support the buoy instrumentation. The US National Science Foundation Ocean Observatories Initiative (OOI) includes arrays of point absorber-like buoy systems used for ocean observation that have been deployed at multiple locations including the Southern Mid-Atlantic Bight. A study is currently underway to design a pitch resonator wave energy converter to supplement existing renewable energy generation for powering observation instrumentation. This paper details field measurements from surface moorings of the OOI Coastal Pioneer Array, which informs the subsequent development of a numerical model for the moored observation system. The model is developed in Wave Energy Converter Simulator (WEC-Sim), which leverages the Simscape multibody solver within the MATLAB/Simulink framework and linear potential flow theory to simulate the hydrodynamic interactions and multibody dynamics in 6 degrees of freedom. Multiple tuning variables are considered to produce a model for the system that matches well with empirical data (about 8% error). In conclusion, the WEC-Sim model will serve as a platform for integrating the pitch resonator wave energy converter concept and deployment preparation (detailed design including power take-off and control systems, response evaluation, etc.).

hydrodynamic modeling↗

DOE OSTI2026

Perturbative second-order optical susceptibility of bulk materials: a symmetry-enforced return to non-orthogonal localized basis sets

The second-order optical susceptibility of semiconductors $\chi^{(2)}_{ijk}(-2\omega;\omega,\omega)$ finds application in metrology, spectroscopy, telecommunications, material characterization, and quantum information. Pioneering calculations of $\chi^{(2)}_{ijk}(-2\omega;\omega,\omega)$ utilized non-orthogonal Gaussian orbitals centered at atoms. That formulation transitioned into plane-wave-based algorithms as time went by. As of late, nevertheless, multiple tools for calculating optical susceptibilities have recast the problem using Wannier (i.e. localized) orbitals, making a comeback onto frameworks based on localized basis sets. Here, in this work, we present an approach for calculating $\chi^{(2)}_{ijk}(-2\omega;\omega,\omega)$ reliant on numerical pseudo-atomic orbitals (PAOs) within perturbation theory in the velocity gauge. Its salient feature is a calculation of ‘Slater–Koster-like’ two-center integrals of the momentum operator in between PAOs identified by symmetry. The approach was successfully tested on paradigmatic cubic silicon carbide (3C-SiC) and gallium arsenide, for which linear responses are contributed as well.

Huamán, Angiolo [Univ. of Arkansas, Fayetteville, ↗