Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TIME DEPENDENCE”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Localization dynamics in a centrally coupled system

In systems in which interactions couple a central degree of freedom and a bath, one would expect signatures of the bath's phase to be reflected in the dynamics of the central degree of freedom. This has been recently explored in connection with many-body localized baths coupled with a central qubit or a single-cavity mode - systems with growing experimental relevance in various platforms. Such models also have an interesting connection with Floquet many-body localization via quantizing the external drive, although this has been relatively unexplored. Here we adapt the multilayer multiconfigurational time-dependent Hartree (ML-MCTDH) method, a well-known tree tensor network algorithm, to numerically simulate the dynamics of a central degree of freedom, represented by a d-level system (qudit), coupled to a disordered interacting one-dimensional spin bath. ML-MCTDH allows us to reach ≈10 2 lattice sites, a far larger system size than what is feasible with exact diagonalization or kernel polynomial methods. From the intermediate time dynamics, we find a well-defined thermodynamic limit for the qudit dynamics upon appropriate rescaling of the system-bath coupling. The spin system shows similar scaling collapse in the Edward-Anderson spin-glass order parameter or entanglement entropy at relatively short times. At longer timescales, we see slow growth of the entanglement, which may arise from dephasing mechanisms in the localized system or long-range interactions mediated by the central degree of freedom. Similar signs of localization are shown to appear as well with unscaled system-bath coupling.

1-dimensional spin chains↗
DOE OSTI2023

Photoionization of C 60 at high energies

Calculations of the photoionization cross section of C60 including all valence and core subshells have been carried out for photon energies up to 1 keV using density functional theory (DFT) and time-dependent DFT (TDDFT) within the framework of a fully molecular model. The high-energy valence subshell cross sections behave rather differently than the results from model potentials where the cross section falls far too rapidly with energy. This unphysical behavior of the model calculations is traced to smearing out of the carbon nuclei in the model potentials which makes it difficult for momentum to be conserved at the higher energies. Comparison with 60 times carbon atomic cross section sheds light on similarities and differences. The high energy behavior of individual valence cross sections reflects the amount of C 2s/2p mixing in the molecular orbitals.

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

Valence-to-core X-ray emission spectroscopy of vanadium oxide and lithiated vanadyl phosphate materials

In this work, we report valence-to-core (VTC) X-ray emission spectroscopy (XES) measurements and theoretical calculations of the electrochemical sequence ε-VOPO 4 , ε-LiVOPO 4 , ε-Li 2 VOPO 4 and the reference oxides V 2 O 3 , VO 2 , and V 2 O 5 . In our analysis of these results, we establish a framework for interrogating chemical bonding that is generally applicable to a wide range of systems, including complex, extended inorganic compounds. While this latter regime has garnered less focused application than, e.g. , metalloenzymes in many excellent catalysis studies, we show that the technique provides high utility in materials-focused energy storage research. Here, sensitivities to the local atomic structure and hybridization schemes are discussed in detail. Similarly, the effect of lithiation on oxidation, delocalization, and shifts in ligand valence energy levels are all readily apparent in the analyzed results. Finally, the TDDFT projections clearly reveal the directional dependencies of the valence band at each of the vanadium sites. Our results demonstrate laboratory-based X-ray spectroscopy instrumentation is a viable route for attaining well-resolved VTC-XES features for inorganic compounds of 3d transition metals, even for samples of limited quantity or suffering from sensitivity to the atmosphere. The experimental results are in good agreement with results produced by real-space Green's function and time-dependent density functional theory (TDDFT) methods, respectively. Hence, we propose that VTC-XES, when equipped with appropriate theoretical support, can be a valuable complement to X-ray absorption pre-edge features for more detailed characterization of a compound's electronic structure. We expect similar analyses will find application in a broad range of materials chemistry research and provide both fundamental and applied insights.

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

A Feynman-Kac based numerical method for the exit time probability of a class of transport problems

The exit time probability, which gives the likelihood that an initial condition leaves a prescribed region of the phase space of a dynamical system at, or before, a given time, is arguably one of the most natural and important transport problems. In this work, we present an accurate and efficient numerical method for computing this probability for systems described by non-autonomous (time-dependent) stochastic differential equations (SDEs) or their equivalent Fokker-Planck partial differential equations. The method is based on the direct approximation of the Feynman-Kac formula that establishes a link between the adjoint Fokker-Planck equation and the forward SDE. The Feynman-Kac formula is approximated using the Gauss-Hermite quadrature rules and piecewise cubic Hermite interpolating polynomials, and a GPU accelerated matrix representation is used to compute the entire time evolution of the exit time probability using a single pass of the algorithm. The method is unconditionally stable, exhibits second order convergence in space, first order convergence in time, and it is straightforward to parallelize. Applications are presented to the advection diffusion of a passive tracer in a fluid flow exhibiting chaotic advection, and to the runaway acceleration of electrons in a plasma in the presence of an electric field, collisions, and radiation damping. Benchmarks against analytical solutions as well as comparisons with explicit and implicit finite difference standard methods for the adjoint Fokker-Planck equation are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Measurement of C P asymmetries in B 0 → K S 0 K S 0 K S 0 decays at Belle II

We report a measurement of decay-time-dependent charge-parity ( C P ) asymmetries in B 0 → K S 0 K S 0 K S 0 decays. We use 387 × 10 6 B B ¯ pairs collected at the ϒ ( 4 S ) resonance with the Belle II detector at the SuperKEKB asymmetric-energy electron-positron collider. We reconstruct 220 signal events and extract the C P -violating parameters S and C from a fit to the distribution of the decay-time difference between the two B mesons. The resulting confidence region is consistent with previous measurements in B 0 → K S 0 K S 0 K S 0 and B 0 → ( c c ¯ ) K 0 decays and with predictions based on the standard model. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2021

Ambiguities in the hadro-chemical freeze-out of Au+Au collisions at SIS18 energies and how to resolve them

The thermal fit to preliminary HADES data of Au+Au collisions at $\sqrt{s_{NN}}$=2.4 GeV shows two degenerate solutions at T≈50 MeV and T≈70 MeV. The analysis of the same particle yields in a transport simulation of the UrQMD model yields the same features, i.e. two distinct temperatures for the chemical freeze-out. While both solutions yield the same number of hadrons after resonance decays, the feeddown contribution is very different for both cases. This highlights that two systems with different chemical composition can yield the same multiplicities after resonance decays. The nature of these two minima is further investigated by studying the time-dependent particle yields and extracted thermodynamic properties of the UrQMD model. It is confirmed, that the evolution of the high temperature solution resembles cooling and expansion of a hot and dense fireball. The low temperature solution displays an unphysical evolution: heating and compression of matter with a decrease of entropy. These results imply that the thermal model analysis of systems produced in low energy nuclear collisions is ambiguous but can be interpreted by taking also the time evolution and resonance contributions into account.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

The Dynamics of Debris Disk Creation in Neutron Star Mergers

Abstract The detection of GW170817/AT2017gfo inaugurated an era of multimessenger astrophysics, in which gravitational-wave and multiwavelength photon observations complement one another to provide unique insight into astrophysical systems. A broad theoretical consensus exists, in which the photon phenomenology of neutron star mergers largely rests upon the evolution of the small amount of matter left on bound orbits around the black hole or massive neutron star remaining after the merger. Because this accretion disk is far from inflow equilibrium, its subsequent evolution depends very strongly on its initial state, yet very little is known about how this state is determined. Using both snapshot and tracer particle data from a numerical relativity/MHD simulation of an equal-mass neutron star merger that collapses to a black hole, we show how gravitational forces arising in a nonaxisymmetric, dynamical spacetime supplement hydrodynamical effects in shaping the initial structure of the bound debris disk. The work done by hydrodynamical forces is ∼10 times greater than that due to time-dependent gravity. Although gravitational torques prior to remnant relaxation are an order of magnitude larger than hydrodynamical torques, their intrinsic sign symmetry leads to strong cancellation; as a result, hydrodynamical and gravitational torques have a comparable effect. We also show that the debris disk’s initial specific angular momentum distribution is sharply peaked at roughly the specific angular momentum of the merged neutron star’s outer layers, a fewr g c, and identify the regulating mechanism.

Astronomy & Astrophysics↗
DOE OSTI2023

Optically controlling the competition between spin flips and intersite spin transfer in a Heusler half-metal on sub–100-fs time scales

The direct manipulation of spins via light may provide a path toward ultrafast energy-efficient devices. However, distinguishing the microscopic processes that can occur during ultrafast laser excitation in magnetic alloys is challenging. Here, we study the Heusler compound Co 2 MnGa, a material that exhibits very strong light-induced spin transfers across the entire M-edge. By combining the element specificity of extreme ultraviolet high-harmonic probes with time-dependent density functional theory, we disentangle the competition between three ultrafast light-induced processes that occur in Co 2 MnGa: same-site Co-Co spin transfer, intersite Co-Mn spin transfer, and ultrafast spin flips mediated by spin-orbit coupling. By measuring the dynamic magnetic asymmetry across the entire M-edges of the two magnetic sublattices involved, we uncover the relative dominance of these processes at different probe energy regions and times during the laser pulse. Our combined approach enables a comprehensive microscopic interpretation of laser-induced magnetization dynamics on time scales shorter than 100 femtoseconds.

36 MATERIALS SCIENCE↗
DOE OSTI2024

SUNDIALS time integrators for exascale applications with many independent systems of ordinary differential equations

Many complex systems can be accurately modeled as a set of coupled time-dependent partial differential equations (PDEs). However, solving such equations can be prohibitively expensive, easily taxing the world’s largest supercomputers. One pragmatic strategy for attacking such problems is to split the PDEs into components that can more easily be solved in isolation. This operator splitting approach is used ubiquitously across scientific domains, and in many cases leads to a set of ordinary differential equations (ODEs) that need to be solved as part of a larger “outer-loop” time-stepping approach. The SUNDIALS library provides a plethora of robust time integration algorithms for solving ODEs, and the U.S. Department of Energy Exascale Computing Project (ECP) has supported its extension to applications on exascale-capable computing hardware. In this paper, we highlight some SUNDIALS capabilities and its deployment in combustion and cosmology application codes (Pele and Nyx, respectively) where operator splitting gives rise to numerous, small ODE systems that must be solved concurrently.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Dual-Color Optical Recording of Bioelectric Potentials by Polymer Electrochromism

Optical recording based on voltage-sensitive fluorescent reporters allows for spatial flexibility of measuring from desired cells, but photobleaching and phototoxicity of the fluorescent labels often limit their sensitivity and recording duration. Voltage-dependent optical absorption, rather than fluorescence, of electrochromic materials, would overcome these limitations to achieve long-term optical recording of bioelectrical signals. Electrochromic materials such as PEDOT:PSS possess the property that an applied voltage can either increase or decrease the light absorption depending on the wavelength. In this work, we harness this anticorrelated light absorption at two different wavelengths to significantly improve the signal detection. With dual-color detection, electrical activity from cells produces signals of opposite polarity, while artifacts, mechanical motions, and technical noises are uncorrelated or positively correlated. Using this technique, we are able to optically record cardiac action potentials with a high signal-to-noise ratio, 10 kHz sampling rate, >15 min recording duration, and no time-dependent degradation of the signal. Furthermore, we can reliably perform multiple recording sessions from the same culture for over 25 days.

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

Non-equilibrium rate theory for polariton relaxation dynamics

We derive an analytic expression of the non-equilibrium Fermi’s golden rule (NE-FGR) expression for a Holstein–Tavis–Cumming Hamiltonian, a universal model for many molecules collectively coupled to the optical cavity. These NE-FGR expressions capture the full-time-dependent behavior of the rate constant for transitions from polariton states to dark states. The rate is shown to be reduced to the well-known frequency domain-based equilibrium Fermi’s golden rule (E-FGR) expression in the equilibrium and collective limit and is shown to retain the same scaling with the number of sites in non-equilibrium and non-collective cases. We use these NE-FGR to perform population dynamics with a time-non-local and time-local quantum master equation and obtain accurate population dynamics from the initially occupied upper or lower polariton states. Furthermore, NE-FGR significantly improves the accuracy of the population dynamics when starting from the lower polariton compared to the E-FGR theory, highlighting the importance of the non-Markovian behavior and the short-time transient behavior in the transition rate constant.

Chemical dynamics↗
DOE OSTI2023

Quantum tracking control of the orientation of symmetric-top molecules

The goal of quantum tracking control is to identify shaped fields to steer observable expectation values along designated time-dependent tracks. The fields are determined via an iteration-free procedure, which is based on inverting the underlying dynamical equations governing the controlled observables. In this paper, we generalize the ideas in [Phys. Rev. A 98, 043429 (2018)] to the task of orienting symmetric top molecules in three dimensions. To this end, we derive equations for the control fields capable of directly tracking the expected value of the three-dimensional dipole orientation vector along a desired path in time. In conclusion, we show this framework can be utilized for tracking the orientation of linear molecules as well, and present numerical illustrations of these principles for symmetric-top tracking control problems.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Elucidation of Ultrafast Decay, Vibrational Beating, and Slow Decay Processes of Excited Azulene

Azulene’s nonradiative decay dynamics and kinetics from its singlet excited-states were studied using mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) combined with trajectory surface-hopping nonadiabatic molecular dynamics (NAMD) and Fermi’s golden rule (FGR) rate theory. The NAMD dynamics reproduce experimental observations that the S 1 → S 0 decay is accelerated and exhibits a crossover from mono- to biexponential kinetics with increasing excess vibrational energy. Minimum-energy-path analyses reveal a continuous S 1 /S 0 conical-intersection seam slightly above the S 1 minimum, providing readily accessible funnels. Time-resolved normal-mode projections reveal selective energy funneling into C–C stretching modes at 1272, 1528, and 1692 cm –1 . Constructive combinations of the latter two modes appear to promote rapid access to the conical-intersection seam, whereas their beating at ∼ 160 fs imposes a natural limit on the decrease of the decay time upon a further increase of excess energies, suggesting that their interference delays the decay. Here, the FGR rate calculation data for the S 1 → S 0 transition reaffirm that the nonradiative decay proceeds mainly near or through the conical intersection, rather than via simple nonadiabatic derivative coupling around the minimum of S 1 . On the other hand, the FGR rates indicate that S 2 decays predominantly to S 0 , as the S 2 –S 1 vibronic coupling is exceptionally weak, which serves as a primary cause for azulene’s characteristic anti-Kasha emission.

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

DAmodel: hierarchical Bayesian modelling of DA white dwarfs for spectrophotometric calibration

We use hierarchical Bayesian modelling to calibrate a network of 32 all-sky faint DA white dwarf (DA WD) spectrophotometric standards (⁠16.5 < V , 19.5⁠) alongside three CALSPEC standards, from 912 Å to 32 μm. The framework is the first of its kind to jointly infer photometric zero points and WD parameters (surface gravity log g⁠, effective temperature T eff ⁠, extinction A V ⁠, dust relation parameter R V ) by simultaneously modelling both photometric and spectroscopic data. We model panchromatic Hubble Space Telescope Wide Field Camera 3 (HST/WFC3) UVIS and IR photometry, HST/STIS UV spectroscopy, and ground-based optical spectroscopy to sub-per cent precision. Photometric residuals for the sample are the lowest yet yielding < 0.004 mag RMS on average from the UV to the NIR, achieved by jointly inferring time-dependent changes in system sensitivity and WFC3/IR count-rate nonlinearity. Our GPU-accelerated implementation enables efficient sampling via Hamiltonian Monte Carlo, critical for exploring the high-dimensional posterior space. The hierarchical nature of the model enables population analysis of intrinsic WD and dust parameters. Inferred spectral energy distributions from this model will be essential for calibrating the James Webb Space Telescope as well as next-generation surveys, including Vera Rubin Observatory’s Legacy Survey of Space and Time and the Nancy Grace Roman Space Telescope.

methods: statistical↗
DOE OSTI2021

COVID-19 dynamics across the US: A deep learning study of human mobility and social behavior

This paper presents a deep learning framework for epidemiology system identification from noisy and sparse observations with quantified uncertainty. The proposed approach employs an ensemble of deep neural networks to infer the time-dependent reproduction number of an infectious disease by formulating a tensor-based multi-step loss function that allows us to efficiently calibrate the model on multiple observed trajectories. The method is applied to a mobility and social behavior-based SEIR model of COVID-19 spread. The model is trained on Google and Unacast mobility data spanning a period of 66 days, and is able to yield accurate future forecasts of COVID-19 spread in 203 US counties within a time-window of 15 days. Interestingly, a sensitivity analysis that assesses the importance of different mobility and social behavior parameters reveals that attendance of close places, including workplaces, residential, and retail and recreational locations, has the largest impact on the effective reproduction number. Furthermore, the model enables us to rapidly probe and quantify the effects of government interventions, such as lock-down and re-opening strategies. Taken together, the proposed framework provides a robust workflow for data-driven epidemiology model discovery under uncertainty and produces probabilistic forecasts for the evolution of a pandemic that can judiciously provide information for policy and decision making. All codes and data accompanying this manuscript are available at https://github.com/PredictiveIntelligenceLab/DeepCOVID19.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

On the dynamics of the fluoroethylene carbonate generated solid electrolyte interphase on silicon anodes during calendar life aging

Here, the widespread use of silicon (Si)-rich anodes in lithium-ion batteries (LIBs) is impeded by an unstable solid electrolyte interphase (SEI) incurring insufficient cell life. Fluoroethylene carbonate (FEC) additive in the electrolyte significantly improves cycle life. However, the gains on calendar life remain unclear; the SEI structure still undergoes detrimental alterations at rest. Thus, elucidating the SEI dynamics during calendar aging is critical to mitigating time-dependent capacity degradation. ATR-FTIR, XPS, and ToF-SIMS are used herein to investigate the SEI structure before and after calendar aging. Si cycled without FEC exhibits no notable SEI chemistry changes Pre- and Post-aging, leaving poor passivation as the main failure pathway. Conversely, the FEC-SEI starts as short oligomeric species from FEC/EC electroreduction prior to aging; after calendar aging, polymerized carbonates become consistently more prominent. Unexpectedly, the deposition of self-polymerized FEC species results from time exposure to the delithiated Si specifically as opposed to the lithiated surface. This unexpected finding is supported by another recent Si calendar-aging research, which albeit not investigating FEC, finds global failure of the SEI upon delithiation resulting in ∼247 fold more reactive surface compared to the lithiated.

Batteries↗
DOE OSTI2022

Conjugation in Isomeric Cyclosilane Thioethers

Conjugation in thioether-functionalized cyclosilanes depends on both conformation and connectivity. Borane-catalyzed dehydrocoupling of thiophenols with partially hydrogenated cyclosilanes proceeded with high chemoselectivity for Si–H bonds in the presence of Si–Si bonds. Isomeric 1,3- and 1,4-functionalized cyclosilanes afforded isomeric cyclosilane thioethers that also differed in preference for chair or twist-boat conformations as demonstrated by both X-ray crystal structures and density functional theory (DFT) calculations. Here, the confounding effects of different conformations and connectivity could be distinguished with time-dependent DFT (TD-DFT) calculations and validated with experimentally obtained UV–vis spectra.

36 MATERIALS SCIENCE↗
DOE OSTI2025

An Optimization-Based Coupling of Reduced Order Models with an Efficient Reduced Adjoint Basis Generation Approach

Optimization-based coupling (OBC) is an attractive alternative to traditional Lagrange multiplier approaches in multiple modeling and simulation contexts. However, application of OBC to time-dependent problems has been hindered by the computational cost of finding the stationary points of the associated Lagrangian, which requires primal and adjoint solves. This issue can be mitigated by using OBC in conjunction with computationally efficient reduced order models (ROMs). To demonstrate the potential of this combination, in this paper, we develop an optimization-based ROM-ROM coupling for a transient advection-diffusion transmission problem. We pursue the “optimize-then-reduce” path toward solving the minimization problem at each time step and solve reduced space adjoint system of equations, where the main challenge in this formulation is the generation of adjoint snapshots and reduced bases for the adjoint systems required by the optimizer. One of the main contributions of the paper is a new technique for an efficient adjoint snapshot collection for gradient-based optimizers in the context of optimization-based ROM-ROM couplings. In conclusion, we present numerical studies demonstrating the accuracy of the approach along with comparison between various approaches for selecting a reduced order basis for the adjoint systems, including decay of snapshot energy, average iteration counts, and timings.

coupled problems↗