SEARCH · Engineering Papers

Results for “master equation”

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 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Group-kinetic theory and modeling of atmospheric turbulence

A group kinetic method is developed for analyzing eddy transport properties and relaxation to equilibrium. The purpose is to derive the spectral structure of turbulence in incompressible and compressible media. Of particular interest are: direct and inverse cascade, boundary layer turbulence, Rossby wave turbulence, two phase turbulence; compressible turbulence, and soliton turbulence. Soliton turbulence can be found in large scale turbulence, turbulence connected with surface gravity waves and nonlinear propagation of acoustical and optical waves. By letting the pressure gradient represent the elementary interaction among fluid elements and by raising the Navier-Stokes equation to higher dimensionality, the master equation was obtained for the description of the microdynamical state of turbulence.

Tchen, C. M.↗

NASA NTRS1994

Rotational relaxation of molecular hydrogen at moderate temperatures

Using a coupled rotation-vibration-dissociation model the rotational relaxation times for molecular hydrogen as a function of final temperature (500-5000 K), in a hypothetical scenario of sudden compression, are computed. The theoretical model is based on a master equation solver. The bound-bound and bound-free transition rates have been computed using a quasiclassical trajectory method. A review of the available experimental data on the rotational relaxation of hydrogen is presented, with a critical overview of the method of measurements and data reduction, including the sources of errors. These experimental data are then compared with the computed results.

COUPLED ROTATION VIBRATION DIS↗

DOE OSTI2025

Current-based metrology with two-terminal mesoscopic conductors

The traditional approach to quantum parameter estimation focuses on the quantum state, deriving fundamental bounds on precision through the quantum Fisher information. In most experimental settings, however, performing arbitrary quantum measurements is highly unfeasible. In open quantum systems, an alternative approach to metrology involves the measurement of stochastic currents flowing from the system to its environment. However, the present understanding of current-based metrology is mostly limited to Markovian master equations. Considering a parameter estimation problem in a two-terminal mesoscopic conductor, we identify the key elements that determine estimation precision within the Landauer-Büttiker formalism. Crucially, this approach allows us to address arbitrary coupling and temperature regimes. Furthermore, we obtain analytical results for the precision in linear-response and zero-temperature regimes. For the specific parameter estimation task that we consider, we demonstrate that the boxcar transmission function is optimal for current-based metrology in all parameter regimes.

Landauer formula↗

DOE OSTI2024

Unraveling Electronic and Vibrational Coherences Following a Charge Transfer Process in a Photosystem II Reaction Center

A reaction center is a unique biological system that performs the initial charge separation within a Photosystem II (PSII) multiunit enzyme, which eventually drives the catalytic water-splitting in plants and algae. The possible role of quantum coherences coinciding with the energy and charge transfer processes in PSII reaction center is one of the active areas of research. Here, we study these quantum coherences by using a numerically exact method on an excitonic dimer model, including linear vibronic coupling and employing optimal parameters from experimental two-dimensional coherent spectroscopic measurements. This enables us to precisely capture the excitonic interaction between pigments and the dissipation of the energy from electronic and charge-transfer (CT) states to the protein environment. We employ the time nonlocal (TNL) quantum master equation to calculate the population dynamics, which yields numerically reliable results. The calculated results show that, due to the strong dissipation, the lifetime of electronic coherence is too short to have direct participation in the charge transfer processes. However, there are long-lived vibrational coherences present in the system at frequencies close to the excitionic energy gap. These are strongly coupled with the electronic coherences, which makes the detection of the electronic coherences with conventional techniques very challenging. Additionally, we unravel the strong excitonic interaction of radical pair (PD1 and PD2) in the reaction center, which results in a long-lived electronic coherence of >100 fs, even at room temperature. Our work provide important physical insight to the charge separation process in PSII reaction center, which may be helpful for better understanding of photophysical processes in other natural and artificial light-harvesting systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2026

Mechanism, Thermochemistry, and Kinetics for the CH + N 2 Reaction Leading to Prompt NO Formation in Combustion

Here, the reaction of CH + N 2 forming H + NCN is a remarkable example of activation of the nitrogen triple bond and is an important source of prompt NO in combustion. The reaction pathway is complex and proceeds through two competing mechanisms: a cyclic addition channel initiated by c-HC(NN) and a chain-addition channel initiated by HCNN, both of which eventually form HNCN prior to dissociation to H + NCN. This work reinvestigates this reaction with composite coupled cluster protocols, including a novel spin-flip equation of motion coupled cluster scheme, combined with pragmatic two-dimensional master equation simulations of the resulting rate coefficients. These improved calculations predict the CH + N 2 rate coefficient between the two more recent previous theoretical results and reduce the uncertainties of the best theoretical models of this reaction to less than a factor of 1.3. Additionally, we provide a closer theoretical investigation of the simultaneous dependence of the CH + N 2 rate coefficient on pressure and temperature, and affirm that collisionally stabilized HNCN, another potential source of prompt NO, emerges as an appreciable product of this reaction under conditions relevant to automotive internal combustion engines and aircraft gas turbine engines.

Nguyen, Thanh Lam [Univ. of Florida, Gainesville, ↗

NASA NTRS1989

The thermal catastrophe model of substorms

Equations are presented for the physical mechanisms involved in the resonant absorption of Alfven waves in the plasma sheet boundary layer. It is shown that energy absorbed by the plasma-sheet particles is a function of the central plasma sheet temperature. The heating curve, when coupled with convective transport, yielded an equation of state for the steady state plasma sheet, whose solution has the form of a mathematical catastrophe. The master equation includes dynamic terms describing the transition across the thermal catastrophe, making it possible to evaluate the time scale for the catastrophe to occur.

Goertz, C. K.↗

NASA NTRS1984

Theory and modeling of atmospheric turbulence, part 2

Two dimensional geostrophic turbulence driven by a random force is investigated. Based on the Liouville equation, which simulates the primitive hydrodynamical equations, a group-kinetic theory of turbulence is developed and the kinetic equation of the scaled singlet distribution is derived. The kinetic equation is transformed into an equation of spectral balance in the equilibrium and non-equilibrium states. Comparison is made between the propagators and the Green's functions in the case of the non-asymptotic quasi-linear equation to prove the equivalence of both kinds of approximations used to describe perturbed trajectories of plasma turbulence. The microdynamical state of fluid turbulence is described by a hydrodynamical system and transformed into a master equation analogous to the Vlasov equation for plasma turbulence. The spectral balance for the velocity fluctuations of individual components shows that the scaled pressure strain correlation and the cascade transfer are two transport functions that play the most important roles.

Chen, C. M.↗

DOE OSTI2024

Classifying topology in photonic crystal slabs with radiative environments

Abstract In the recent years, photonic Chern materials have attracted substantial interest as they feature topological edge states that are robust against disorder, promising to realize defect-agnostic integrated photonic crystal slab devices. However, the out-of-plane radiative losses in those photonic Chern slabs has been previously neglected, yielding limited accuracy for predictions of these systems’ topological protection. Here, we develop a general framework for measuring the topological protection in photonic systems, such as in photonic crystal slabs, while accounting for in-plane and out-of-plane radiative losses. Our approach relies on the spectral localizer that combines the position and Hamiltonian matrices of the system to draw a real-picture of the system’s topology. This operator-based approach to topology allows us to use an effective Hamiltonian directly derived from the full-wave Maxwell equations after discretization via finite-elements method (FEM), resulting in the full account of all the system’s physical processes. As the spectral FEM-localizer is constructed solely from FEM discretization of the system’s master equation, the proposed framework is applicable to any physical system and is compatible with commonly used FEM software. Moving forward, we anticipate the generality of the method to aid in the topological classification of a broad range of complex physical systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

NASA NTRS2014

Calculation of State Specific Rate Coefficients for Non-Equilibrium Hypersonics Applications: from H(Psi) = E(Psi) to k(T) = A *exp(-E(sub a)/RT)

Development of High-Fidelity Physics-Based Models to describe hypersonic flight through the atmospheres of Earth and Mars is underway at NASA Ames Research Center. The goal is to construct chemistry models of the collisional and radiative processes that occur in the bow shock and boundary layers of spacecraft during atmospheric entry that are free of empiricism. In this talk I will discuss our philosophy and describe some of our progress. Topics to be covered include thermochemistry, internal energy relaxation, collisional dissociation and radiative emission and absorption. For this work we start by solving the Schrodinger equation to obtain accurate interaction potentials and radiative properties. Then we invoke classical mechanics to compute state-specific heavy particle collision cross sections and reaction rate coefficients. Finally, phenomenological rate coefficients and relaxation times are determined from master equation solutions.

aerothermodynamics↗

DOE OSTI2026

A fractional calculus framework for open quantum dynamics: From Liouville to Lindblad to memory kernels

Open quantum systems exhibit dynamics ranging from unitary evolution to irreversible dissipation. While the Gorini–Kossakowski–Sudarshan–Lindblad equation uniquely characterizes Markovian completely positive and trace-preserving (CPTP) evolution, many physical platforms display non-Markovian features such as algebraic relaxation and coherence backflow. Fractional calculus provides a natural way to model such long-memory behavior through power-law temporal kernels introduced by fractional time derivatives. Here, we develop a unified framework that embeds fractional master equations within the broader hierarchy of open-system formalisms. The fractional equation forms a structured subclass of memory-kernel models, reduces to the Lindblad form at unit order, and, through Bochner–Phillips subordination, admits a CPTP representation as an average over Lindblad semigroups. Its resolvent structure further connects fractional dynamics to established non-Markovian approaches, including Nakajima–Zwanzig kernels and hierarchical equations of motion, providing a compact surrogate for long-memory effects. This formulation positions fractional calculus as a rigorous and practical language for modeling non-Markovian quantum dynamics in chemical physics and physical chemistry, providing a CPTP-preserving, computationally efficient surrogate for structured condensed-phase environments where long-time memory and dissipation play a central role.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Quantum simulation of Lindbladian dynamics via repeated interactions

The Lindblad equation generalizes the Schrödinger equation to quantum systems that undergo dissipative dynamics. The quantum simulation of Lindbladian dynamics is therefore non-unitary, preventing a naive application of state-of-the-art quantum algorithms. Here, we make use of an approximate correspondence between Lindbladian dynamics and evolution based on repeated interaction (RI) CPTP maps to write down a Hamiltonian formulation of the Lindblad dynamics and derive a rigorous error bound on the master equation. Specifically, we show that the number of interactions needed to simulate the Liouvillian within error e scales in most physical scenarios as . This is significant because the error in the Lindbladian approximation to the dynamics is not explicitly bounded in existing quantum algorithms for open system simulations. We then provide quantum algorithms to simulate RI maps using an iterative qubitization approach and Trotter–Suzuki formulas, and specifically show that for iterative qubitization the number of operations needed to simulate the dynamics (for a fixed value of ?) scales as in the limit where a0 (the coefficient 1-norm for the system and bath Hamiltonians) asymptotically dominates over the corresponding factor for the interaction Hamiltonian, which is often the case in weak coupling. This scaling would appear to be optimal if the complexity of ? is not considered, which underscores the importance of considering the error in the Liouvillian that we reveal in this work.

Quantum Computing↗

NASA NTRS1984

Basic governing equations for the flight regimes of aeroassisted orbital transfer vehicles

The basic governing equations for the low-density, high-enthalpy flow regimes expected in the shock layers over the heat shields of the proposed aeroassisted orbital transfer vehicles are derived by combining and extending existing theories. The conservation equations are derived from gas kinetic principles for a four-component ionized gas consisting of neutral molecules, neutral atoms, singly ionized ions, and electrons, assuming a continuum flow. The differences among translational-rotational, vibrational, and electron temperatures are accounted for, as well as chemical nonequilibrium and electric-charge separation. Expressions for convective and viscous fluxes, transport properties, and the terms representing interactions among various energy modes are given explicitly. The expressions for the rate of electron-vibration energy transfer, which violates the Landau-Teller conditions, is derived by solving the system of master equations accounting for the multiple-level transitions.

Lee, J.-H.↗

NASA NTRS1985

Basic Governing Equations for the Flight Regimes of Aeroassisted Orbital Transfer Vehicles

The basic governing equations for the low-density, high-enthalpy flow regimes expected in the shock layers over the heat shields of the proposed aeroassisted orbital transfer vehicles are derived by combining and extending existing theories. The conservation equations are derived from gas kinetic principles for a four-component ionized gas consisting of neutral molecules, neutral atoms, singly ionized ions, and electrons, assuming a continuum flow. The differences among translational-rotational, vibrational, and electron temperatures are accounted for, as well as chemical nonequilibrium and electric-charge separation. Expressions for convective and viscous fluxes, transport properties, and the terms representing interactions among various energy modes are explicitly given. The expressions for the rate of electron-vibration energy transfer, which violates the Landau-Teller conditions, are derived by solving the system of master equations accounting for the multiple-level transitions.

Lee, Jong-Hun↗

NASA NTRS2016

Comparison of Quantum Mechanical and Empirical Potential Energy Surfaces and Computed Rate Coefficients for N2 Dissociation

Physics-based modeling of hypersonic flows is predicated on the availability of chemical reaction rate coefficients and cross sections for the collisional processes. This approach has been built around the use of quantum mechanical calculations to describe the interaction between the colliding particles. In this approach a potential energy surface (PES) is computed by solving the electronic Schrödinger equation and collision cross sections are determined for that PES using classical, semiclassical or quantum mechanical scattering methods. The rate coefficients are computed by integrating the thermally weighted cross sections. State-to-state rate coefficients are determined by only integrating over a thermal distribution of collisional energies. Finally, thermal rate coefficients are determined by summation of the state-to-state rate coefficients for reactions of molecules in all relevant ro-vibrational energy levels. If the flow is in thermal non-equilibrium, the translational, vibrational and rotational energy modes can be represented in different ways: three unique temperatures can be used to describe the distributions, the populations of individual ro-vibrational energy levels can be determined by solving the Master Equation, or through the use of direct simulation in particle-based Monte Carlo sampling. The PES-to-rate coefficient approach had been proposed and attempted in the early days of digital computing, but it is only in the last 15 years that computer hardware and software have been up to the task of calculating accurate interatomic and intermolecular potentials.

Jaffe, Richard L.↗

DOE OSTI2025

Efficient Simulation of Open Quantum Systems on NISQ Trapped‐Ion Hardware

Abstract Simulating open quantum systems, which interact with external environments, presents significant challenges on noisy intermediate‐scale quantum (NISQ) devices due to limited qubit resources and noise. In this study, an efficient framework is proposed for simulating open quantum systems on NISQ hardware by leveraging a time‐perturbative Kraus operator representation of the system's dynamics. This approach avoids the computationally expensive Trotterization method and exploits the Lindblad master equation to represent time evolution in a compact form, particularly for systems satisfying specific commutation relations. The efficiency of this method is demonstrated by simulating quantum channels, such as the continuous‐time Pauli channel and damped harmonic oscillators, on NISQ trapped‐ion hardware, including IonQ Harmony and Quantinuum H1‐1. Additionally, hardware‐agnostic error mitigation techniques are introduced, including Pauli channel fitting and quantum depolarizing channel inversion, to enhance the fidelity of quantum simulations. These results show strong agreement between the simulations on real quantum hardware and exact solutions, highlighting the potential of Kraus‐based methods for scalable and accurate simulation of open quantum systems on NISQ devices. This framework opens pathways for simulating more complex systems under realistic conditions in the near term.

Burdine, Colin [Department of Electrical and Compu↗

DOE OSTI2024

Ladder symmetries and Love numbers of Reissner-Nordström black holes

It is well known that asymptotically flat black holes in general relativity have vanishing tidal Love numbers. In the case of Schwarzschild and Kerr black holes, this property has been shown to be a consequence of a hidden structure of ladder symmetries for the perturbations. In this work, we extend the ladder symmetries to non-rotating charged black holes in general relativity. As opposed to previous works in this context, we adopt a more general definition of Love numbers, including quadratic operators that mix gravitational and electromagnetic perturbations in the point-particle effective field theory. We show that the calculation of a subset of those couplings in full general relativity is affected by an ambiguity in the split between source and response, which we resolve through an analytic continuation. As a result, we derive a novel master equation that unifies scalar, electromagnetic and gravitational perturbations around Reissner-Nordström black holes. The equation is hypergeometric and can be obtained from previous formulations via nontrivial field redefinitions, which allow to systematically remove some of the singularities and make the presence of the ladder symmetries more manifest.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Non -degenerate marginal-likelihood calibration with application to quantum characterization

Here, we propose a marginal likelihood strategy within the Kennedy-O’Hagan (KOH) Bayesian framework, where a Gaussian process (GP) models the discrepancy between a physical system and its simulator. Our approach introduces a novel marginalized likelihood by integrating out the degenerate eigenspace of the covariance matrix, rather than approximating the original likelihood. Unlike approximation methods that compromise accuracy for computational efficiency, our method defines an exact likelihood—distinct from the original but preserving all relevant information. This formulation achieves computational efficiency and stability, even for large datasets where the covariance matrix nears degeneracy. Applied to the characterization of a superconducting quantum device at Lawrence Livermore National Laboratory, the approach enhances the predictive accuracy of the Lindblad master equations for modeling Ramsey measurement data by effectively quantifying uncertainties consistent with the quantum data.

general physics↗

DOE OSTI2024

On the prediction of pressure effects for the combination kinetics of two alkyl radicals with the geometric mean rule

Here, alkyl combination reactions, along with their reverse alkane unimolecular decomposition reactions, are widely present in pyrolysis and combustion systems. These reactions occur mainly at temperatures at which their rate coefficients depend on pressure, and only a few studies have explored this dependence for large size alkyls. The accurate calculation of these kinetic data is limited by the size of the alkyl fragments involved, since the theoretical methods used are computationally expensive. Therefore, little reliable data is available for the modeling of real fuels alkane decomposition kinetics. In this paper, we propose an alternative approach that allows an accurate estimation of pressure-dependent kinetic rates for the combination of long-chain alkyl radicals, based on tabulated reaction models. The rates of combination of linear alkyls radicals forming ethane to pentane were calculated using the CASPT2-F12/cc-pVDZ-F12 level within VRC-TST theory combined with master-equation simulations. The computed high-pressure limit rate coefficients showed that n-propyl could be a minimum-size representative structure for larger n-alkyl chain, but this approach fails to capture the pressure dependence of the reaction that remains size-dependent for larger alkyls. The Geometric Mean Rule (GMR), a simple rule that links the self and cross-reactions of three alkyl radicals, is used to address this issue. Our theoretical calculations show for the first time that the GMR applies not only at the high-pressure limit, but also for every pressure explored. The GMR is shown to be accurate, with a maximum 22 % deviation in the falloff region for the estimation of rates for small alkyl combinations. Using this method, we propose a tabulated reaction rate rule for the combination of alkyl radicals up to C5+C5. The use of computed rate coefficients with the GMR approximation opens ways to develop pressure-dependent reaction rate rules for large n-alkyl combinations, based on accurate and computationally expensive theoretical calculations on smaller alkyl fragments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗