Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “box relaxation”

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.

14 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

One-dimensional phase retrieval: Regularization, box relaxation and uniqueness

Recovering a signal from its Fourier magnitude is referred to as phase retrieval, which occurs in different fields of engineering and applied physics. This paper gives a new characterization of the phase retrieval problem. Particularly useful is the analysis revealing that the common gradient-based regularization does not restrict the set of solutions to a smaller set. Specifically focusing on binary signals, we show that a box relaxation is equivalent to the binary constraint for Fourier-types of phase retrieval. Here, we further prove that binary signals can be recovered uniquely up to trivial ambiguities under certain conditions. Finally, we use the characterization theorem to develop an efficient denoising algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Quantum dynamics of a nucleon in the Fermi accelerator

Highlights: • In nuclear physics, the Fermi accelerator is a coherent open quantum system. • There is not a true equilibrium for a nucleon‘s mean-energy in the Fermi accelerator. • Occupation of higher energy eigen-states remains even after the mean-energy plateaus. • This work paves the way for work on quantum tunnelling in the Fermi accelerator. The quantum dynamics of a particle in a one-dimensional box with an oscillating wall (the Fermi accelerator) is investigated. The model is applied to the motion of a single nucleon in the mean-field potential of a heavy atomic nucleus whose surface vibrates. By directly solving the time-dependent Schrödinger equation, both the state of the particle and its mean-energy are studied. The effects of the frequency of the wall oscillation on the nucleon’s energy are addressed. Its energy oscillates in phase with the moving wall for all frequencies, showing no chaotic behaviour. There is a large initial peak of the nucleon’s energy as the particle adjusts to the sudden change in the size of the box and a varying relaxation time as it plateaus towards lower energy and a partial equilibrium. Small oscillations in energy continue, since there cannot be a true equilibrium while the wall is moving. The quantum coherence between the different parts of the nucleon’s wave-function in real space is very much preserved. This research lays the foundation for future investigations into quantum tunnelling in the Fermi accelerator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

ReaxFF reactive force field model enables accurate prediction of physiochemical and mechanical properties of crystalline and amorphous shape‐memory polyurethane

Abstract We investigate reactive force field (ReaxFF) prediction of physiochemical and mechanical properties, and stimuli response behavior of a shape‐memory polyurethane (SMPU) rigid segment. We used SMPUs as a platform to link the rigid/hard segment domain crystallinity with properties that can also be experimentally measured. Specifically, we developed models for crystalline and amorphous forms of a common rigid segment of SMPUs, for example, 4,4′‐diphenylmethane diisocyanate (MDI) with n‐butanediol (BDO), denoted as MDI‐BDO. ReaxFF molecular dynamics simulations with equilibrium dynamics are performed by controlling mass, temperature, and pressure or volume to predict the structural and physiochemical properties of crystalline and amorphous rigid segments of SMPU systems, followed by simulating their XRD and FTIR patterns, respectively. Non‐equilibrium ReaxFF simulations are implemented by using uniaxial box deformation technique to study response‐behavior of crystalline SMPU under tensile loading and during stress‐relaxation. The model and simulation scheme presented in this study demonstrates a robust theoretical pathway toward application‐dictated, on‐demand design of high‐performance SMPUs in applications such as soft robotics.

Afroz, Mohammad M.↗
DOE OSTI2024

Kinetic relaxation and nucleation of Bose stars in self-interacting wave dark matter

We revisit kinetic relaxation and soliton/boson star nucleation in fuzzy scalar dark matter featuring short-ranged self-interactions H int = − λ | ψ | 4 / 2 m 2 , alongside gravitational self-interactions. We map out the full curve of nucleation timescale for both repulsive ( λ < 0 ) and attractive ( λ > 0 ) short-ranged self-interaction strength and in doing so reveal two new points. Firstly, besides the two usual terms, ∝ G 2 and ∝ λ 2 , in the total relaxation rate Γ relax , there is an additional cross term ∝ G λ arising due to interference between gravitational and short-ranged self-interaction scattering amplitudes. This yields a critical repulsive interaction strength λ cr ≃ − 2 π G m 2 / v 0 2 , at which the relaxation rate is smallest and serves as the transition point between typical net attractive self-interaction ( λ ≳ λ cr ) and net repulsive self-interaction ( − λ ≳ − λ cr ). Secondly, while in the net attractive regime, nucleation timescale is similar to inverse relaxation timescale τ nuc ∼ Γ relax − 1 ; in the net repulsive regime, nucleation occurs at a delayed time τ nuc ∼ ( λ / λ cr ) Γ relax − 1 . We confirm our analytical understanding by performing 3D field simulations with varying average mass density ρ ¯ , box size L and grid size N . Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2025

A Theory of Ultrafast Charge Transfer Relaxation with Non-Innocent Solvent Molecules

In this work we revisit the photodynamics of tetracyanoethylene-hexamethylbenzene (TCNE- HMB), the molecular complex studied by Hochstrasser et al. [J. Chem. Phys. 100, 4797–4810, 1994] that has long challenged the applicability of Marcus’s theory of elec- tron transfer for predicting photochemical reactions. Using a novel black-box electronic structure algorithm (time-dependent density functional theory with one double, TD- DFT-1D) to efficiently run molecular dynamics that can treat charge recombination, we run ab initio surface hopping molecular dynamics and confirm that, for a polar solvent, charge recombination rates can be incredibly fast (indeed faster than the sol- vent relaxation time); for non-polar solvents, the rate is much slower. We demonstrate that, although Marcus theory cannot be directly applied, these nonequilibrium (and sometimes incredibly fast) photoexcited dynamics can be effectively explained within a two-state model without any evidence of a transition through a conical intersection. Most importantly, for this paradigmatic model system, we are able to identify two nuclear coordinates of interest (rather than the single coordinate predicted by Marcus or a full set of internal quantum modes studied by Bixon and Jortner): the solvent relaxation in the first shell (that strongly modulates the energies of the charge trans- fer state and differentiates time scales for relaxation) and a nuclear displacement in the TCNE-HMB complex arising from a handful of vibrations that induces non-Born Oppenheimer motion and eventually facilitates an abrupt electronic transition to the ground state. Altogether, these findings suggest a tractable generalization of Marcus theory for future simulations of photochemistry with non-innocent solvent environ- ments in the spirit of a Hamiltonian suggested by Stuchebrukhov (J. Chem. Phys. 107, 3821, 1997).

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

Pressure Effects on the Relaxation of an Excited Ethane Molecule in High-Pressure Bath Gases

Here, we use molecular dynamics to calculate the rotational and vibrational energy relaxation of C 2 H 6 in Ar, Kr, and Xe bath gases over a pressure range of 10 to 400 atm and at temperatures of 300 K and 800 K. The C 2 H 6 is instantaneously excited by 80 kcal/mol randomly distributed into both vibrational and rotational modes. The computed relaxation rates show little sensitivity to the identity of the noble gas in the bath. Vibrational relaxation rates show a non-linear pressure dependence at 300 K. At 800 K the reduced range of bath gas densities covered by the range of pressures do not yet show any non-linearity in the pressure dependence. Rotational relaxation is characterized with two relaxation rates. The slower rate is comparable to the vibrational relaxation rate. The faster rate has a linear pressure dependence at 300 K but an irregular, nonlinear pressure dependence at 800 K. To understand this, a model was developed based on approximating the periodic box used in the molecular dynamics simulations by an equal-volume collection of cubes where each cube is sized to allow only single occupancy by the noble gas or the molecule. Combinatorial statistics then leads to a pressure and temperature dependent analytic distribution of the bath gas species the molecule encounters in a collision. This distribution, the dissociation energy of molecule/bath gas complexes and bath gas clusters, and the computed energy release per collision combine to show that only at 300 K is the energy release sufficient to dissociate likely complexes and clusters. This suggests that persistent and pressure-dependent clusters and complexes at 800 K may be responsible for the non-linear pressure dependence of rotational relaxation.

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

Scalable semidefinite programming approach to variational embedding for quantum many-body problems

In quantum embedding theories, a quantum many-body system is divided into localized clusters of sites which are treated with an accurate ‘high-level’ theory and glued together self-consistently by a less accurate ‘low-level’ theory at the global scale. The recently introduced variational embedding approach for quantum many-body problems combines the insights of semidefinite relaxation and quantum embedding theory to provide a lower bound on the ground-state energy that improves as the cluster size is increased. The variational embedding method is formulated as a semidefinite program (SDP), which can suffer from poor computational scaling when treated with black-box solvers. Here, we exploit the interpretation of this SDP as an embedding method to develop an algorithm which alternates parallelizable local updates of the high-level quantities with updates that enforce the low-level global constraints. Moreover, we show how translation invariance in lattice systems can be exploited to reduce the complexity of projecting a key matrix to the positive semidefinite cone.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Glutamine-rich regions of the disordered CREB transactivation domain mediate dynamic intra- and intermolecular interactions

The cyclic AMP response element (CRE) binding protein (CREB) is a transcription factor that contains a 280-residue N-terminal transactivation domain and a basic leucine zipper that mediates interaction with DNA. The transactivation domain comprises three subdomains, the glutamine-rich domains Q1 and Q2 and the kinase inducible activation domain (KID). NMR chemical shifts show that the isolated subdomains are intrinsically disordered but have a propensity to populate local elements of secondary structure. The Q1 and Q2 domains exhibit a propensity for formation of short β-hairpin motifs that function as binding sites for glutamine-rich sequences. These motifs mediate intramolecular interactions between the CREB Q1 and Q2 domains as well as intermolecular interactions with the glutamine-rich Q1 domain of the TATA-box binding protein associated factor 4 (TAF4) subunit of transcription factor IID (TFIID). Using small-angle X-ray scattering, NMR, and single-molecule Förster resonance energy transfer, we show that the Q1, Q2, and KID regions remain dynamically disordered in a full-length CREB transactivation domain (CREB TAD ) construct. The CREB TAD polypeptide chain is largely extended although some compaction is evident in the KID and Q2 domains. Paramagnetic relaxation enhancement reveals transient long-range contacts both within and between the Q1 and Q2 domains while the intervening KID domain is largely devoid of intramolecular interactions. Phosphorylation results in expansion of the KID domain, presumably making it more accessible for binding the CBP/p300 transcriptional coactivators. In conclusion, our study reveals the complex nature of the interactions within the intrinsically disordered transactivation domain of CREB and provides molecular-level insights into dynamic and transient interactions mediated by the glutamine-rich domains.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Scattering phase shift in quantum mechanics on quantum computers

Here, we investigate the feasibility of extracting infinite volume scattering phase shift on quantum computers in a simple one-dimensional quantum mechanical model, using the formalism established in the work by Guo and Gasparian [Phys. Rev. D 108, 074504 (2023)] that relates the integrated correlation functions for a trapped system to the infinite volume scattering phase shifts through a weighted integral. The system is first discretized in a finite box with periodic boundary conditions, and the formalism in real time is verified by employing a contact interaction potential with exact solutions. Quantum circuits are then designed and constructed to implement the formalism on current quantum computing architectures. To overcome the fast oscillatory behavior of the integrated correlation functions in real-time simulation, different methods of postdata analysis are proposed and discussed. Test results on IBM hardware show that good agreement can be achieved with two qubits, but complete failure ensues with three qubits due to two-qubit gate operation errors and thermal relaxation errors.

Guo, Peng [Dakota State Univ., Madison, SD (United↗
DOE OSTI2025

The hierarchical growth of bright central galaxies and intracluster light as traced by the magnitude gap

Using a sample of 2800 galaxy clusters identified in the Dark Energy Survey across the redshift range 0.20 < z < 0.60, we characterize the hierarchical assembly of bright central galaxies (BCGs) and the surrounding intracluster light (ICL). To quantify hierarchical formation we use the stellar mass–halo mass (SMHM) relation, comparing the halo mass, estimated via the mass–richness relation, to the stellar mass within the BCG + ICL system. Moreover, we incorporate the magnitude gap (M14), the difference in brightness between the BCG (measured within 30 kpc) and fourth brightest cluster member galaxy within 0.5 $R_{200,c}$, as a third parameter in this linear relation. The inclusion of M14, which traces BCG hierarchical growth, increases the slope and decreases the intrinsic scatter, highlighting that it is a latent variable within the BCG + ICL SMHM relation. Moreover, the correlation with M14 decreases at large radii. However, the stellar light within the BCG + ICL transition region (30 –80 kpc) most strongly correlates with halo mass and has a statistically significant correlation with M14. Since the transition region and M14 are independent measurements, the transition region may grow due to the BCG’s hierarchical formation. Additionally, as M14 and ICL result from hierarchical growth, we use a stacked sample and find that clusters with large M14 values are characterized by larger ICL and BCG + ICL fractions, which illustrates that the merger processes that build the BCG stellar mass also grow the ICL. Furthermore, this may suggest that M14 combined with the ICL fraction can identify dynamically relaxed clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Constraining the phase shift of relativistic species in DESI BAOs

In the early Universe, neutrinos decouple quickly from the primordial plasma and propagate without further interactions. The impact of free-streaming neutrinos is to create a temporal shift in the gravitational potential that impacts the acoustic waves known as baryon acoustic oscillations (BAOs), resulting in a non-linear spatial shift in the Fourier-space BAO signal. In this work, we make use of and extend upon an existing methodology to measure the phase shift amplitude $\beta _{\phi }$ and apply it to the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) BAOs with an anisotropic BAO fitting pipeline. We validate the fitting methodology by testing the pipeline with two publicly available fitting codes applied to highly precise cubic box simulations and realistic simulations representative of the DESI DR1 data. We find further study towards the methods used in fitting the BAO signal will be necessary to ensure accurate constraints on $\beta _{\phi }$ in future DESI data releases. Using DESI DR1, we present individual measurements of the anisotropic BAO distortion parameters and the $\beta _{\phi }$ for the different tracers, and additionally a combined fit to $\beta _{\phi }$ resulting in $\beta _{\phi } = 2.7 \pm 1.7$. After including a prior on the distortion parameters from constraints using Planck we find $\beta _{\phi } = 2.7^{+0.60}_{-0.67}$ suggesting $\beta _{\phi } > 0$ at 4.3$\sigma$ significance. This result may hint at a phase shift that is not purely sourced from the standard model expectation for $N_{\rm {eff}}$ or could be a upwards statistical fluctuation in the measured $\beta _{\phi }$; this result relaxes in models with additional freedom beyond Lambda-cold dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Gradient-informed Hamiltonian Monte Carlo for multicomponent CALPHAD model optimization and uncertainty quantification

CALPHAD model parameter optimization is inherently challenging due to non-smooth objective functions, high-dimensional parameter spaces, and the need for uncertainty quantification (UQ). Traditional weighted nonlinear least squares approaches are computationally efficient but local, whereas black-box global optimizers and ensemble Markov Chain Monte Carlo (MCMC) methods provide broader exploration at substantial computational cost. The objective of this work is to combine the global exploration capability of gradient-informed Hamiltonian Monte Carlo – specifically the No-U-Turn Sampler (NUTS) – with local deterministic refinement using BFGS to efficiently optimize multicomponent CALPHAD models with minimal manual intervention. Analytic gradients are computed via the Jansson derivative framework. The methodology is demonstrated on the Cr—Fe binary system and extended to the Cr—Fe—Ni ternary system with 32 degrees of freedom. For Cr—Fe, NUTS achieves comparable or superior optimality relative to ensemble MCMC while requiring over an order-of-magnitude fewer likelihood evaluations. Parameter uncertainties are quantified through NUTS sampling and propagated to thermodynamic observables using local expansion, demonstrating a novel modular approach that combines binary and ternary parameter subsets without requiring global relaxation. These results establish gradient-informed exploration as a scalable strategy for multicomponent CALPHAD optimization and provide a practical route towards efficient higher-order database development with quantified uncertainty.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Systematic Construction of Time-Dependent Hamiltonians for Microwave-Driven Josephson Circuits

Time-dependent electromagnetic drives are fundamental for controlling complex quantum systems, including superconducting Josephson circuits. In these devices, accurate time-dependent Hamiltonian models are imperative for predicting their dynamics and designing high-fidelity quantum operations. Existing numerical methods, such as black-box quantization (BBQ) and energy-participation ratio (EPR), excel at modeling the static Hamiltonians of Josephson circuits. However, these techniques do not fully capture the behavior of driven circuits stimulated by external microwave drives, nor do they include a generalized approach to account for the inevitable noise and dissipation that enter through microwave ports. Here, we introduce numerical techniques that leverage classical microwave simulations, efficiently executable in finite-element solvers, to obtain the time-dependent Hamiltonian of microwave-driven superconducting circuits with arbitrary geometries under charge, flux, or mixed electromagnetic modulation. Importantly, our techniques do not rely on a lumped-element description of the superconducting circuit, in contrast to previous approaches to tackling this problem. We demonstrate the versatility of our approach by characterizing the driven properties of realistic circuit devices in complex electromagnetic environments, including coherent dynamics due to charge and flux modulation, as well as drive-induced relaxation and dephasing. Our techniques offer a powerful toolbox for optimizing circuit designs and advancing practical applications in superconducting quantum computing.

Lu, Yao [Yale U.; Yale U. (main); Fermilab] (ORCID↗
DOE OSTI2025

Constraining gas motion and non-thermal pressure beyond the core of the Abell 2029 galaxy cluster with XRISM

We report on a detailed spectroscopic study of the gas dynamics and hydrostatic mass bias of the galaxy cluster Abell 2029, utilizing high-resolution observations from XRISM Resolve. Abell 2029, known for its cool core and relaxed X-ray morphology, provides an excellent opportunity to investigate the influence of gas motions beyond the central region. Expanding upon prior studies that revealed low turbulence and bulk motions within the core, our analysis covers regions out to the scale radius $R_{2500}$ (670 kpc) based on three radial pointings extending from the cluster center toward the northern side. We obtain accurate measurements of bulk and turbulent velocities along the line of sight. The results indicate that non-thermal pressure accounts for no more than 2% of the total pressure at all radii, with a gradual decrease outward. The observed radial trend differs from many numerical simulations, which often predict an increase in non-thermal pressure fraction at larger radii. These findings suggest that deviations from hydrostatic equilibrium are small, leading to a hydrostatic mass bias of around 2% across the observed area.

X-rays: galaxies: clusters↗