Search NASA⌕ Search

SEARCH · Search NASA

Results for “Relaxation method”

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

Accurate numerical simulations of open quantum systems using spectral tensor trains

Decoherence between qubits is a major bottleneck in quantum computations. Decoherence results from intrinsic quantum and thermal fluctuations as well as noise in the external fields that perform the measurement and preparation processes. With prescribed colored noise spectra for intrinsic and extrinsic noise, we present a numerical method, Quantum Accelerated Stochastic Propagator Evaluation (Q-ASPEN), to solve the time-dependent noise-averaged reduced density matrix in the presence of intrinsic and extrinsic noise. Q-ASPEN is arbitrarily accurate and can be applied to provide estimates for the resources needed to error-correct quantum computations. We employ spectral tensor trains, which combine the advantages of tensor networks and pseudospectral methods, as a variational ansatz to the quantum relaxation problem and optimize the ansatz using methods typically used to train neural networks. Here, the spectral tensor trains in Q-ASPEN make accurate calculations with tens of quantum levels feasible. We present benchmarks for Q-ASPEN on the spin-boson model in the presence of intrinsic noise and on a quantum chain of up to 32 sites in the presence of extrinsic noise. In our benchmark, the memory cost of Q-ASPEN scales as a low-order polynomial in the size of the system once the number of system states surpasses the number of basis functions used in the spectral expansion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Data-driven equation-free dynamics applied to many-protein complexes: The microtubule tip relaxation

Microtubules (MTs) constitute the largest components of the eukaryotic cytoskeleton and play crucial roles in various cellular processes, including mitosis and intracellular transport. The property allowing MTs to cater to such diverse roles is attributed to dynamic instability, which is coupled to the hydrolysis of GTP (guanosine-5'-triphosphate) to GDP (guanosine-5'-diphosphate) within the β-tubulin monomers. Understanding the equilibrium dynamics and the structural features of both GDP- and GTP-complexed MT tips, especially at an all-atom level, remains challenging for both experimental and computational methods because of their dynamic nature and the prohibitive computational demands of simulating large, many-protein systems. This study employs the “equation-free” multiscale computational method to accelerate the relaxation of all-atom simulations of MT tips toward their putative equilibrium conformation. Using large MT lattice systems (14 protofilaments × 8 heterodimers) comprising ~21-38 million atoms, we applied this multiscale approach to leapfrog through time and nearly double the computational efficiency in realizing relaxed all-atom conformations of GDP- and GTP-complexed MT tips. Commencing from an initial 4 μs unbiased all-atom simulation, we interleave coarse projective “equation-free” jumps with short bursts of all-atom molecular dynamics simulation to realize an additional effective simulation time of 1.875 μs. Our 5.875 μs of effective simulation trajectories for each system expose the subtle yet essential differences in the structures of MT tips as a function of whether β-tubulin monomer is complexed with GDP or GTP, as well as the lateral interactions within the MT tip, offering a refined understanding of features underlying MT dynamic instability. Furthermore, the approach presents a robust and generalizable framework for future explorations of large biomolecular systems at atomic resolution.

Wu, Jiangbo [University of Chicago, IL (United Sta↗
DOE OSTI2024

Hole Relaxation Bottlenecks in CdSe/CdTe/CdSe Lateral Heterostructures Lead to Bicolor Emission

Concentric lateral CdSe/CdTe/CdSe heterostructures show bicolor photoluminescence from both a red charge transfer band of the CdSe/CdTe interface and a green fluorescence from CdSe. Further, this work uses visible and near-infrared transient spectroscopy measurements to demonstrate that the deviation from Kasha's rule arises from a hole relaxation bottleneck from CdSe to CdTe. Hole transfer can take up to 1 ns, which permits radiative relaxation of excitons remaining in CdSe. Simulations indicate that the hole relaxation bottleneck arises due to the sparse density of states and poor spatial overlap of hole states at energies near the CdSe band edge. The divergent kinetics of transfer for band edge and hot holes is exploited to vary the ratio of green and red photoluminescence with excitation wavelength, providing another knob to control emission color. These findings support the use of lateral heterojunctions as a method for slowing carrier relaxation in two-dimensional materials.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

An improved guess for the variational calculation of charge-transfer excitations in large systems

Ab initio quantum-chemical methods that perform well for computing the electronic ground state are not straightforwardly transferable to electronically excited states, particularly in large molecular systems. Wave function theory offers high accuracy, but is often prohibitively expensive. Methods based on time-dependent density functional theory (TD-DFT) are crucially sensitive to the chosen exchange-correlation functional (XCF) parameterization, and system-specific tuning protocols were therefore proposed to address the method's robustness. Methods based on the variational relaxation of the excited-state electron density showcased promising results for the calculation of charge-transfer excitations, but the complex shape of the electronic hypersurface makes convergence to a specific excited state much more difficult than for the ground state when standard variational techniques are applied. We address the latter aspect by providing suitable initial guesses, which we obtain by two separate constrained algorithms. Combined with the squared-gradient minimization algorithm for all-electrons relaxation in a freeze-and-release scheme (FRZ-SGM), we demonstrate that orbital-optimized density functional theory (OO-DFT) calculations can reliably converge to the charge-transfer states of interest even for large molecular systems. We test the FRZ-SGM method on a phenothiazine-anthraquinone CT excitation in a supramolecular Pd(II) coordination cage complex as a function of the cage conformation. This compound has been studied experimentally prior to our work. We compare this freeze-and-release scheme to two XCF reparameterizations, which were recently proposed as low-cost TD-DFT-based alternatives to variational methods. Two dye-semiconductor complexes, which were previously investigated in the context of photovoltaic applications, serve as a second example to investigate the convergence and stability of the FRZ-SGM approach. Our results demonstrate that FRZ-SGM provides reliable convergence for charge-transfer excited states and avoids variational collapse to lower-lying electronic states, whereas time-dependent DFT calculations with an adequate tuning procedure for the range-separation parameter provide a computationally efficient initial estimate of the corresponding energies, with a computational cost comparable to that of configuration-interaction singles (CIS) calculations.

Bogo, Nicola↗
DOE OSTI2026

Two-Level System Spectroscopy from Correlated Multilevel Relaxation in Superconducting Qubits

Transmon qubits are a cornerstone of modern superconducting quantum computing platforms. Temporal fluctuations of energy relaxation in these qubits are widely attributed to microscopic two-level systems (TLSs) in device dielectrics and interfaces, yet isolating individual defects typically relies on tuning the qubit or the TLS into resonance. We demonstrate a novel spectroscopy method for fixed-frequency transmons based on multilevel relaxation: repeated preparation of the second excited state and simultaneous $T_1$ extraction of the first and second excited states reveals characteristic correlations in the decay rates of adjacent transitions. From these correlations we identify one or more dominant TLSs and reconstruct their frequency drift over time. Remarkably, we find that TLSs detuned by $\gtrsim 100\,\mathrm{MHz}$ from the qubit transition can still significantly influence relaxation. The proposed method provides a powerful tool for TLS spectroscopy without the need to tune the transmon frequency, either via a flux-tunable inductor or AC-Stark shifts.

Roy, Tanay [Fermilab] (ORCID:000000019442862X)↗
DOE OSTI2025

The Electron Thermal Conductivity of Pu and Zr Substituted $\mathcal{γ}$-U

Uranium alloys are attractive recycled nuclear fuels because of their high thermal conductivity (𝑘) and fissile density. Limited experimental studies of the 𝑘 of U-Pu-Zr alloys in the range of 15 to 20 wt% Pu and 6 to 15 wt% Zr indicate that increasing the content of either Zr or Pu tends to lower 𝑘. However, which element has the greater effect on 𝑘, and the associated mechanisms, remains unclear. Here, in this study, the electron thermal conductivity (𝑘 𝑒 ) of U-Pu-Zr compositions are calculated using density functional theory. The electronic structure is evaluated to understand the effects of plutonium (Pu) and zirconium (Zr) substitution on the 𝑘 𝑒 of 𝛾-U. Alloys of up to 37.5 at. % Pu and 37.5 at. % Zr are examined. Two methods are applied to calculate 𝑘 𝑒 ; we find that the accuracy of each method depends on the electronic and mass similarities between the solute and solvent atoms. Specifically, when the solute atom is similar in electronic structure and mass, the more accurate method is that which employs the electron relaxation time of 𝛾-U, while if the elements are dissimilar, a mixed method that mixes several parameters associated with JNW_S⁢3033426825100132 from each element in the alloy is best. The introduction of all alloying elements decreases 𝑘 𝑒 ; however, in binary compounds, Pu and Zr have different effects. Pu flattens the electronic bands but compensates for this deleterious effect by increasing electron density near the Fermi level. Zr flattens the electronic bands more severely without adding electron density near the Fermi level. Therefore, Zr decreases 𝑘 𝑒 more than Pu in binary compounds. In ternary compounds, the difference between Pu and Zr is minimal due to the phononic change from the large mass change of Zr substitution, even at 12.5 at. %. Thus, we predict that higher loadings of Pu, and potentially other actinides, can be added to U-Pu-Zr compositions for faster recycling of spent fuel without sacrificing 𝑘. We also note that these 𝑘 𝑒 calculation methods can be applied to non-fuel alloys that require 𝑘 𝑒 predictions, such as cladding, heat exchanger, and structural materials.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

A Decomposition-Based Learn-To-Optimize Approach with Feasibility Layer Assistance for Sub-Hourly Unit Commitment

Sub-hourly unit commitment (UC) with 15-min intervals is gaining significant attention as a way to respond rapidly to the fluctuations in electricity supply and demand introduced by renewable resources. However, the increased temporal resolution and complex inter-temporal dependencies pose substantial computational challenges for traditional optimization methods. To this end, this paper explores a decomposition-based learn-to-optimize approach. Building on recent advances in machine learning, our method revisits the long- overlooked Lagrangian relaxation framework, which is a classical decomposition technique that enables tractable subproblem solving. These smaller subproblems are inherently well-suited for machine learning, as their reduced dimensionality and structural regularity allow predictive models to efficiently learn and generalize solution patterns. We thus propose a generic predictive model, which embeds Gated Recurrent Units (GRUs) and Attention in the encoder-decoder structure, and integrate a rule-based feasibility layer to capture temporal dependencies, reduce training effort, and improve feasibility w.r.t. unit-level constraints. Our method has been validated on the IEEE 118-bus system, demonstrating promising performance in solving sub-hourly UC problems efficiently and feasibly.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Spectroscopic Signatures of Phonon Character in Molecular Electron Spin Relaxation

Spin–lattice relaxation constitutes a key challenge for the development of quantum technologies, as it destroys superpositions in molecular quantum bits (qubits) and magnetic memory in single molecule magnets (SMMs). Gaining mechanistic insight into the spin relaxation process has proven challenging owing to a lack of spectroscopic observables and contradictions among theoretical models. Here, we use pulse electron paramagnetic resonance (EPR) to profile changes in spin relaxation rates (T 1 ) as a function of both temperature and magnetic field orientation, forming a two-dimensional data matrix. For randomly oriented powder samples, spin relaxation anisotropy changes dramatically with temperature, delineating multiple regimes of relaxation processes for each Cu(II) molecule studied. We show that traditional T 1 fitting approaches cannot reliably extract this information. Single-crystal T 1 anisotropy experiments reveal a surprising change in spin relaxation symmetry between these two regimes. We interpret this switch through the concept of a spin relaxation tensor, enabling discrimination between delocalized lattice phonons and localized molecular vibrations in the two relaxation regimes. Variable-temperature T 1 anisotropy thus provides a unique spectroscopic method to interrogate the character of nuclear motions causing spin relaxation and the loss of quantum information.

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

Near-ideal relaxed MHD in slab geometry

We investigate the solutions of the relaxed magnetohydrodynamic (MHD) model (RxMHD) of R. Dewar and Z. Qu. This model generalizes Taylor relaxation by including the ideal Ohm's law constraint using an augmented Lagrangian method, providing a pathway to extend the multi-region relaxed MHD (MRxMHD) model. We present the first numerical solution of the RxMHD model by Dewar and Qu, demonstrating that it is mathematically well-defined and computationally feasible for constructing MHD equilibria in slab geometry. We also show that a cross-field flow can exist without enforcing an arbitrary constraint on the angular momentum, as is done in the case of MRxMHD with flow. Our results also demonstrate the self-organization of fully relaxed regions during the optimization, which was an important motivation behind developing this model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Robust 3D multi-material hydrodynamics using discontinuous Galerkin methods

A high-order discontinuous Galerkin (DG) method is presented for nonequilibrium multi-material (m ≥ 2) flow with sharp interfaces. Material interfaces are reconstructed using the algebraic THINC approach, resulting in a sharp interface resolution. The system assumes stiff velocity relaxation and pressure nonequilibrium. The presented DG method uses Dubiner's orthogonal basis functions on tetrahedral elements. This results in a unique combination of sharp multimaterial interfaces and high-order accurate solutions in smooth single-material regions. A novel shock indicator based on the interface conservation condition is introduced to mark regions with discontinuities. Slope limiting techniques are applied only in these regions so that nonphysical oscillations are eliminated while maintaining high-order accuracy in smooth regions. A local projection is applied on the limited solution to ensure discrete closure law preservation. The effectiveness of this novel limiting strategy is demonstrated for complex three-dimensional multi-material problems, where robustness of the method is critical. The presented numerical problems demonstrate that more accurate and efficient multi-material solutions can be obtained by the DG method, as compared to second-order finite volume methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

On the functional dependence of transition-potential coupled cluster

Orbital relaxation of the core region is a primary source of error in the computation of core ionization and core excitation energies. Recently, Transition-Potential Coupled Cluster (TP-CC) methods have been used to explicitly treat orbital relaxation using non-variational molecular orbitals determined by reoccupation of orbitals optimized for a fractional core occupation. The amount of fractional occupation is governed by parameter λ, and recommended values for accurate TP-CCSD and XTP-CCSD computations of carbon, nitrogen, oxygen, and fluorine K edges were previously determined. Herein, we explore the performance of several density functionals for generating the fractionally occupied orbitals used in TP-CCSD. These functionals include HF, BP86, BH&HLYP, B3LYP, M06-2X, and ωB97m-V. The fractionally occupied orbitals computed across the various functionals were subsequently employed as the initial orbitals for our TP-CCSD calculations of organic K-edge x-ray absorption and photoelectron spectra. Regardless of the functional used to generate the fractionally occupied orbitals, the TP-CCSD calculations yield accurate and comparable core ionization energies, core excitation energies, and oscillator strengths.

Coupled-cluster methods↗
DOE OSTI2025

𝐴𝑏 initio density-matrix approach to exciton coherence: Phonon scattering, Coulomb interactions, and radiative recombination

Relaxation processes following light excitation in semiconductors are key in materials-based quantum technology applications. These processes are broadly studied in atomically thin transition-metal dichalcogenides, quasi-two-dimensional excitonic semiconductors in which atomistic design allows for tunable excited-state properties, such as relaxation lifetimes and photoinduced coherence. In this work, we present a density-matrix-based approach to compute exciton relaxation within a many-body ab initio perspective. We expand our previously developed Lindblad density-matrix formalism to capture multichannel electron-hole pair relaxation processes, including phonon and Coulomb scattering as well as radiative recombination, and we study their effect on the time-resolved excited-state propagation. Using monolayer MoSe 2 as a prototypical example, we examine many-body effects on the time-dependent dynamics of photoactive excitations, exploring how the electron-hole pair interactions are reflected in variations of the excitation energy, spectral signature, and state coherence. In conclusion, our method supplies a detailed understanding of exciton relaxation mechanisms in realistic materials, offering a previously unexplored pathway to study excited-state dynamics in semiconductors from first principles.

Band structure methods↗
DOE OSTI2025

Role of Electron Correlation beyond the Active Space in Achieving Quantitative Predictions of Spin-Phonon Relaxation

Single-molecule magnets (SMMs) are promising candidates for molecular-scale data storage and processing due to their strong magnetic anisotropy and long spin relaxation times. However, as temperature rises, interactions between electronic states and lattice vibrations accelerate spin relaxation, significantly limiting their practical applications. Recently, ab initio simulations have made it possible to advance our understanding of phonon-induced magnetic relaxation, but significant deviations from experiments have often been observed. The description of molecules’ electronic structure has been mostly based on complete active space self-consistent field (CASSCF) calculations, and the impact of electron correlation beyond the active space remains largely unexplored. In this study, we provide the first systematic investigation of spin-phonon relaxation in SMMs with post-CASSCF multiconfigurational methods, specifically CAS followed by second-order perturbation theory and multiconfiguration pair-density functional theory. Taking Co(II)- and Dy(III)-based SMMs as case studies, we analyze how electron correlation influences spin-phonon relaxation rates across a range of temperatures, comparing theoretical predictions with experimental observations. Our findings demonstrate that post-CASSCF treatments make it possible to achieve quantitative predictions for Co(II)-based SMMs. For Dy(III)-based systems, however, accurate predictions require consideration of additional effects, underscoring the urgent necessity of further advancing the study of the effects of electronic correlation in these complex systems.

Energy↗
DOE OSTI2024

Towards dynamical low-rank approximation for neutrino kinetic equations. Part I: Analysis of an idealized relaxation model

Dynamical low-rank approximation (DLRA) is an emerging tool for reducing computational costs and provides memory savings when solving high-dimensional problems. Here, in this work, we propose and analyze a semi-implicit dynamical low-rank discontinuous Galerkin (DLR-DG) method for the space homogeneous kinetic equation with a relaxation operator, modeling the emission and absorption of particles by a background medium. Both DLRA and the discontinuous Galerkin (DG) scheme can be formulated as Galerkin equations. To ensure their consistency, a weighted DLRA is introduced so that the resulting DLR-DG solution is a solution to the fully discrete DG scheme in a subspace of the standard DG solution space. Similar to the standard DG method, we show that the proposed DLR-DG method is well-posed. We also identify conditions such that the DLR-DG solution converges to the equilibrium. Numerical results are presented to demonstrate the theoretical findings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

G EANT 4 atomic relaxation data for transfermium nuclei (Z = 101–104)

Advanced theoretical methods can accurately calculate various atomic observables and predict electronic structure. Still, systematic computations of the radiative and non-radiative transition probabilities and energies are missing for the actinides and all the transfermium elements. However, these compilations are needed for comprehensive Monte-Carlo simulations (such as GEANT4) of the radioactive decay of transfermium nuclei. These simulations can forma basis for data analysis of experiments, especially with complex detection setups. Investigation of the transfermium nuclei is crucial for understanding the nature of the nuclear force. In this study, simulations based on data from the Jena Atomic Calculator (JAC) and the data from the Evaluated Atomic Data Library (EADL) present in GEANT4 were found compatible for the three elements Ba(Z = 56), U(Z = 92), and Fm(Z = 100), thus, validating the JAC calculations. For Z> 100, we also found sound agreement between simulations that used data generated with JAC and experimental results involving No(Z = 102) and Rf(Z = 104) isotopes. In conclusion, these results demonstrate that JAC can produce reliable atomic data sets for transfermium elements, which will assist in analyzing nuclear-decay-spectroscopy experiments.

GEANT4↗
DOE OSTI2026

Experimental Quantification of Spin–Phonon Coupling in Molecular Qubits Using Inelastic Neutron Scattering

Electronic spin superposition states enable nanoscale sensing through their sensitivity to the local environment, yet their sensitivity to vibrational motion also limits their coherence times. In molecular spin systems, chemical tunability and atomicscale resolution are accompanied by a dense, thermally accessible phonon spectrum that introduces efficient spin relaxation pathways. Despite extensive theoretical work, there is little experimental consensus on which vibrational energies dominate spin relaxation or how molecular structure controls spin−phonon coupling (SPC). We present a fully experimental method to quantify SPC coefficients by combining temperature-dependent vibrational spectra from inelastic neutron scattering with spin relaxation rates measured by electron paramagnetic resonance. We apply this framework to two model S = 1/2 systems, copper(II) phthalocyanine (CuPc) and copper(II) octaethylporphyrin (CuOEP). Two distinct relaxation regimes emerge: below 40 K, weakly coupled lattice modes below 50 cm −1 dominate, whereas above 40 K, optical phonons above ∼185 cm −1 become thermally populated and drive relaxation with SPC coefficients nearly 3 orders of magnitude larger. Structural distortions in CuOEP that break planar symmetry soften the crystal lattice and enhance anharmonic scattering but also raise the energy of stretching modes at the molecular core where the spins reside. This redistributes vibrational energy toward the molecular periphery and out of plane, ultimately reducing SPC relative to CuPc and enabling room-temperature spin coherence in CuOEP. Although our method does not provide mode-specific SPC coefficients, it quantifies contributions from distinct spectral regions and establishes a broadly applicable, fully experimental link between crystal structure, lattice dynamics, and spin relaxation.

Lohaus, Stefan H. [California Institute of Technol↗
DOE OSTI2024

Phonon modal analysis of thermal transport in ThO 2 with point defects using equilibrium molecular dynamics

Defects can significantly degrade the thermal conductivity of ThO 2 , an advanced nuclear fuel material as well as a surrogate for other fluorite-structured materials. Here, we investigate how point defects in ThO 2 impact phonon mode-resolved thermal transport. By incorporating phonon modes from lattice dynamics, we decompose the trajectory and heat flux to phonon normal mode space and extract key phonon properties, including phonon relaxation times and their contributions to thermal conductivity. We implement two methods. The first method is based on the Green Kubo formalism to resolve the contribution of each phonon mode to thermal conductivity. The second resolves the lifetime of individual phonon modes and the thermal conductivity is calculated using the Boltzmann transport equation within relaxation time approximation. Notably, a lower contribution of acoustic modes is revealed compared to perturbative approaches considering only three-phonon scattering processes. The effects of four types of point defects are evaluated. The strongest impact on a reduction in thermal conductivity is from Th interstitials, followed by Th vacancies. O interstitials/vacancies have a similar impact, albeit smaller than defects on the thorium sublattice. These observations are consistent with previous studies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Comparison between explicit and implicit discretization strategies for a dissipative thermal environment

We investigate strategies for simulating open quantum systems coupled to dissipative baths by comparing explicit wave function-based discretization [via multi-layer multi-configuration time-dependent Hartree (ML-MCTDH)] and the implicit density matrix-based master equation method [via tree tensor network hierarchical equations of motion (TTN-HEOM)]. For dissipative baths characterized by exponentially decaying bath correlation functions, the implicit discretization approach of HEOM—rooted in bath correlation function decompositions—proves significantly more efficient than explicit discretization of the bath into discrete harmonic modes. Explicit methods, like ML-MCTDH, require extensive mode discretization to approximate continuum baths, leading to computational bottlenecks. Case studies for two-level systems and a Fenna–Matthews–Olson complex model highlight TTN-HEOM’s superiority in capturing dissipative dynamics with relaxations with a minimal number of auxiliary modes, while the explicit methods are as exact as the HEOM in pure dephasing regimes. This comparison is enabled by the TENSO package, which has both ML-MCTDH and TTN-HEOM implemented using the same computational structure and propagation strategy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗