Search NASA⌕ Search

SEARCH · Search NASA

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

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

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Generalized quantum master equations can improve the accuracy of semiclassical predictions of multitime correlation functions

Multitime quantum correlation functions are central objects in physical science, offering a direct link between the experimental observables and the dynamics of an underlying model. While experiments such as 2D spectroscopy and quantum control can now measure such quantities, the accurate simulation of such responses remains computationally expensive and sometimes impossible, depending on the system’s complexity. A natural tool to employ is the generalized quantum master equation (GQME), which can offer computational savings by extending reference dynamics at a comparatively trivial cost. However, dynamical methods that can tackle chemical systems with atomistic resolution, such as those in the semiclassical hierarchy, often suffer from poor accuracy, limiting the credence one might lend to their results. By combining work on the accuracy-boosting formulation of semiclassical memory kernels with recent work on the multitime GQME, here we show for the first time that one can exploit a multitime semiclassical GQME to dramatically improve both the accuracy of coarse mean-field Ehrenfest dynamics and obtain orders of magnitude efficiency gains.

Chemistry↗
DOE OSTI2025

On the discretization error of the discrete generalized quantum master equation

The transfer tensor method (TTM) [Cerrillo and Cao, Phys. Rev. Lett. 112 , 110401 (2014)] can be considered a discrete-time formulation of the Nakajima–Zwanzig quantum master equation (NZ-QME) for modeling non-Markovian quantum dynamics. A recent paper [Makri, J. Chem. Theory Comput. 21 , 5037 (2025)] raised concerns regarding the consistency of the TTM discretization, particularly a spurious term at the initial time t = 0. Here, this work presents a detailed analysis of the discretization structure of the TTM, clarifying the origin of the initial-time correction and establishing a consistent relationship between the TTM discrete-time memory kernel K N and the continuous-time NZ-QME kernel $\mathscr{K}$( N Δ t ). This relationship is validated numerically using the spin-boson model, demonstrating convergence of reconstructed memory kernels and accurate dynamical evolution as Δ t → 0. While the TTM provides a consistent discretization, we note that alternative schemes are also viable, such as the midpoint derivative/midpoint integral scheme proposed in Makri’s work. The relative performance of various schemes for either computing accurate $\mathscr{K}$( N Δ t ) from exact dynamics or obtaining accurate dynamics from exact $\mathscr{K}$( N Δ t ) warrants further investigation.

Density-matrix↗
DOE OSTI2024

Efficient formulation of multitime generalized quantum master equations: Taming the cost of simulating 2D spectra

Modern 4-wave mixing spectroscopies are expensive to obtain experimentally and computationally. In certain cases, the unfavorable scaling of quantum dynamics problems can be improved using a generalized quantum master equation (GQME) approach. However, the inclusion of multiple (light–matter) interactions complicates the equation of motion and leads to seemingly unavoidable cubic scaling in time. In this paper, we present a formulation that greatly simplifies and reduces the computational cost of previous work that extended the GQME framework to treat arbitrary numbers of quantum measurements. Specifically, we remove the time derivatives of quantum correlation functions from the modified Mori–Nakajima–Zwanzig framework by switching to a discrete-convolution implementation inspired by the transfer tensor approach. We then demonstrate the method’s capabilities by simulating 2D electronic spectra for the excitation-energy-transfer dimer model. In our method, the resolution of data can be arbitrarily coarsened, especially along the t2 axis, which mirrors how the data are obtained experimentally. Even in a modest case, this demands O(103) fewer data points. We are further able to decompose the spectra into one-, two-, and three-time correlations, showing how and when the system enters a Markovian regime where further measurements are unnecessary to predict future spectra and the scaling becomes quadratic. This offers the ability to generate long-time spectra using only short-time data, enabling access to timescales previously beyond the reach of standard methodologies.

Chemistry↗
DOE OSTI2025

QEpsilon v0.1.0

QEpsilon is a Python package designed to minimize the effort required to build a data-driven quantum master equation of an open quantum system and to perform time evolution of the master equation. Applications of QEpsilon span from quantum computing to condensed matter systems.

Xie, Pinchen [Lawrence Berkeley National Laborator↗
DOE OSTI2025

Tree tensor network hierarchical equations of motion based on time-dependent variational principle for efficient open quantum dynamics in structured thermal environments

In this work, we introduce an efficient method, TTN-HEOM, for exactly calculating the open quantum dynamics for driven quantum systems interacting with highly structured bosonic baths by combining the tree tensor network (TTN) decomposition scheme with the bexcitonic generalization of the numerically exact hierarchical equations of motion (HEOM). The method yields a series of quantum master equations for all core tensors in the TTN that efficiently and accurately capture the open quantum dynamics for non-Markovian environments to all orders in the system–bath interaction. These master equations are constructed based on the time-dependent Dirac–Frenkel variational principle, which isolates the optimal dynamics for the core tensors given the TTN ansatz. The dynamics converges to the HEOM when increasing the rank of the core tensors, a limit in which the TTN ansatz becomes exact. We introduce TENSO, tensor equations for non-Markovian structured open systems, as a general-purpose Python code to propagate the TTN-HEOM dynamics. We implement three general propagators for the coupled master equations: two fixed-rank methods that require a constant memory footprint during the dynamics and one adaptive-rank method with a variable memory footprint controlled by the target level of computational error. We exemplify the utility of these methods by simulating a two-level system coupled to a structured bath containing one Drude–Lorentz component and eight Brownian oscillators, which is beyond what can presently be computed using the standard HEOM. Our results show that the TTN-HEOM is capable of simulating both dephasing and relaxation dynamics of driven quantum systems interacting with structured baths, even those of chemical complexity, with an affordable computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2023

Thermalization at low temperatures via weakly damped multisite baths

We study the thermalization properties of one-dimensional open quantum systems coupled to baths at their boundary. The baths are driven to their thermal states via Lindblad operators, while the system undergoes Hamiltonian dynamics. We specifically consider multisite baths and investigate the extent to which the late-time steady state resembles a Gibbs state at some controllable temperature set by the baths. We study three models: a noninteracting fermion model accessible via free-fermion technology, and two interacting models, the XZ model and the chiral clock model, which are accessible via tensor network methods. We show that, by tuning towards the weak coupling and slow relaxation limits, one can engineer low temperatures in the bulk of the system provided the bath size is big enough. Here, we use this capability to study energy transport in the XZ model at lower temperatures than previously reported. Our work paves the way for future studies of interacting open quantum systems at low temperatures.

1-dimensional spin chains↗
DOE OSTI2025

Quantum algorithm to simulate Lindblad master equations

We present a quantum algorithm for simulating a family of Markovian master equations that can be realized through a probabilistic application of unitary channels and state preparation. Our approach employs a second-order product formula for the Lindblad master equation, achieved by decomposing the dynamics into dissipative and Hamiltonian components and replacing the dissipative segments with randomly compiled, easily implementable elements. The sampling approach eliminates the need for ancillary qubits to simulate the dissipation process and reduces the gate complexity in terms of the number of jump operators. We provide a rigorous performance analysis of the algorithm. We also extend the algorithm to time-dependent Lindblad equations, generalize the family of Markovian master equations it can be applied to, and explore applications beyond the Markovian noise model. A new error bound, in terms of the diamond norm, for second-order product formulas for time-dependent Liouvillians is provided that might be of independent interest. Published by the American Physical Society 2025

Borras, Evan (ORCID:000900017709037X)↗
DOE OSTI2024

Polariton creation in coupled cavity arrays with spectrally disordered emitters

Abstract Integrated photonics has been a promising platform for analog quantum simulation of condensed matter phenomena in strongly correlated systems. To that end, we explore the implementation of all-photonic quantum simulators in coupled cavity arrays with integrated ensembles of spectrally disordered emitters. Our model is reflective of color center ensembles integrated into photonic crystal cavity arrays. Using the Quantum Master equation and the Effective Hamiltonian approaches, we study energy band formation and wavefunction properties in the open quantum Tavis–Cummings–Hubbard framework. We find conditions for polariton creation and (de)localization under experimentally relevant values of disorder in emitter frequencies, cavity resonance frequencies, and emitter-cavity coupling rates. To quantify these properties, we introduce two metrics, the polaritonic and nodal participation ratios, that characterize the light-matter hybridization and the node delocalization of the wavefunction, respectively. These new metrics combined with the Effective Hamiltonian approach prove to be a powerful toolbox for cavity quantum electrodynamical engineering of solid-state systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS2016

Computational Role of Tunneling in a Programmable Quantum Annealer

Quantum tunneling is a phenomenon in which a quantum state tunnels through energy barriers above the energy of the state itself. Tunneling has been hypothesized as an advantageous physical resource for optimization. Here we present the first experimental evidence of a computational role of multiqubit quantum tunneling in the evolution of a programmable quantum annealer. We developed a theoretical model based on a NIBA Quantum Master Equation to describe the multi-qubit dissipative cotunneling effects under the complex noise characteristics of such quantum devices.We start by considering a computational primitive, the simplest non-convex optimization problem consisting of just one global and one local minimum. The quantum evolutions enable tunneling to the global minimum while the corresponding classical paths are trapped in a false minimum. In our study the non-convex potentials are realized by frustrated networks of qubit clusters with strong intra-cluster coupling. We show that the collective effect of the quantum environment is suppressed in the critical phase during the evolution where quantum tunneling decides the right path to solution. In a later stage dissipation facilitates the multiqubit cotunneling leading to the solution state. The predictions of the model accurately describe the experimental data from the D-WaveII quantum annealer at NASA Ames. In our computational primitive the temperature dependence of the probability of success in the quantum model is opposite to that of the classical paths with thermal hopping. Specially, we provide an analysis of an optimization problem with sixteen qubits,demonstrating eight qubit cotunneling that increases success probabilities. Furthermore, we report results for larger problems with up to 200 qubits that contain the primitive as subproblems.

hard problems↗
DOE OSTI2025

Simulating quantum-classical interfaces via the Lindblad master equation

In hybrid quantum systems, the interface between quantum and classical domains is essential for the generation, control, and measurement of quantum states. Quantum-classical interfaces (QCIs) are ubiquitous in devices such as optical modulators, quantum sensors, and signal processors, where classical signals influence quantum dynamics. In this paper, we employ the Lindblad master equation to simulate the evolution of a quantum system interacting with a classical control system. Our model captures both linear and nonlinear interactions by incorporating first- and second-order susceptibilities, and it quantifies the influence of externally applied control parameters on decoherence and state evolution. As an illustrative example, we analyze an optical modulator and demonstrate how variations in material response and drive conditions affect photon statistics, coherence, and phase-space distributions. In conclusion, the findings offer a path to an all-encompassing model for understanding and optimizing QCIs, with wide-ranging implications for the performance, design, and robustness of next-generation quantum devices.

Quantum engineering↗
DOE OSTI2025

Molecular axis distribution moments in ultrafast transient absorption spectroscopy: A path toward ultrafast quantum state tomography

In ultrafast time-resolved experiments with gas phase molecules, the alignment of the molecular axis relative to the polarization of the interacting laser pulses plays a crucial role in determining the dynamics following this light–matter interaction. The molecular axis distribution is influenced by the interacting pulses and is intrinsically linked to the electronic coherences of the excited molecules. However, in typical theoretical calculations of such interactions, the signal is either calculated for a single molecule in the molecular frame or averaged over all possible molecular orientations to compare with the experiment. Such averaging removes information about anisotropy in the molecular-axis distribution, even though anisotropic contributions can play a significant role in the measured experimental signal. Here, we calculate the laboratory frame transient electronic first-order polarization [P (1) ] spectra in terms of separated molecular frame and laboratory frame quantities. The laboratory frame polarizations are compared with orientation-averaged quantum master equation calculations, demonstrating that orientation-averaging captures only the isotropic contributions. We show that our formalism also allows us to evaluate the anisotropic contributions to the spectrum. Lastly, we discuss the application of this approach to achieve ultrafast quantum state tomography using transient absorption spectroscopy and field observables in nonlinear spectroscopy.

74 ATOMIC AND MOLECULAR PHYSICS↗
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 OSTI2025

Universal time scalings of sensitivity in Markovian quantum metrology

Assuming Markovian time evolution of a quantum sensing system, we study the general characterization of the optimal sensitivity scalings with time, under most general quantum control protocols. We allow the estimated parameter to influence both the Hamiltonian as well as the dissipative part of the quantum master equation and focus on the asymptotic-time along with the short-time sensitivity scalings. We find that via simple algebraic conditions (in terms of the Hamiltonian, the jump operators as well as their parameter derivatives), one can characterize the four classes of metrological models that represent: quadratic-linear, quadratic-quadratic, linear-linear, and linear-quadratic time scalings. We also investigate the relevant time scales on which the transition between the two regimes appears. Additionally, we provide universal numerical methods to obtain quantitative bounds on sensitivity that are the tightest that exist in the literature. Simplicity and universality of our results make it suitable for diverse applications in quantum metrology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

On the Mechanism of Polaritonic Rate Suppression from Quantum Transition Paths

Polariton chemistry holds promise for facilitating mode-selective chemical reactions, but the underlying mechanism behind the rate modifications observed under strong vibrational coupling is not well understood. Using the recently developed quantum transition path theory, we have uncovered a mechanism of resonant suppression of a thermal reaction rate in a simple model polaritonic system consisting of a reactive mode in a bath confined to a lossless microcavity with a single photon mode. Here we observed the formation of a polariton during rate-limiting transitions on reactive pathways and identified the concomitant rate suppression as being due to hybridization between the reactive mode and the cavity mode, which inhibits bath-mediated tunneling. The transition probabilities that define the quantum master equation can be directly translated into a visualization of the corresponding polariton energy landscape. This landscape exhibits a double funnel structure with a large barrier between the initial and final states.

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

Interplay between Ultrafast Electronic and Librational Dynamics in Liquid Nitrobenzene Probed with Two-Color Four-Wave Mixing

Here, we present an experimental and theoretical study of the interplay between ultrafast electron dynamics and librational dynamics in liquid nitrobenzene. A femtosecond ultraviolet pulse and two femtosecond near-infrared pulses interact with nitrobenzene molecules, generating a four-wave mixing nonlinear signal measured in the Optical Kerr Effect geometry. The signal is measured to be nonzero only at negative time delays, corresponding to the near-infrared pulses arriving before the ultraviolet pulse. We perform time-dependent Quantum Master Equation calculations with classical libration to simulate the experiment. The simulations support the conclusion that the near-infrared pulses launch librational motion while creating electronic coherences resulting in a libration-modulated electronic nonlinear response. The analysis of the phase-matched four-wave mixing signals suggests a nonparametric process leaving the molecules in an excited electronic state, providing new insight into ultrafast nonlinear optical interactions in liquids and advancing toward probing ultrafast electronic coherences in complex molecular liquids.

Shivaram, Niranjan [Lawrence Berkeley National Lab↗
DOE OSTI2025

Engineering long-lived entanglement through dissipation in quantum hybrid solid-state platforms

Spin squeezing, a form of many-body entanglement, is a crucial resource in quantum metrology and information processing. While experimentally viable protocols for generating stable spin squeezing have been proposed in quantum optics setups, there is growing interest in quantum hybrid solid-state systems as alternative platforms for both engineering and exploring many-body quantum phenomena. In this work, we propose a scheme to generate long-lived spin squeezing in an ensemble of solid-state qubits interacting with electromagnetic noise emitted by a squeezed solid-state bath. We identify the conditions under which quantum correlations within the bath can be transferred to the qubit array, driving it into an entangled state independently of its initial configuration. To assess the experimental feasibility of our approach, we analyze the dynamics of an array of solid-state spin defects coupled to a common ferromagnetic bath, which is driven into a non-equilibrium squeezed state through its interaction with a surface acoustic wave mode. Our results demonstrate that the ensemble can exhibit steady-state spin squeezing under suitable conditions, opening new pathways for the generation of robust many-body entanglement in solid-state spin ensembles.

NV centers↗
DOE OSTI2026

Quantum speed limit for the out-of-time-ordered correlator from an open-system perspective

Scrambling, the delocalization of initially localized quantum information, is commonly characterized by the out-of-time-ordered correlator (OTOC). Employing the OTOC–Renyi-2 entropy theorem, we derive a quantum speed limit for the OTOC, which sets a lower bound for the rate with which information can be scrambled. This bound becomes particularly tractable by describing the scrambling of information in a closed quantum system as an effective decoherence process of an open system interacting with an environment. We prove that decay of the OTOC can be bounded by the strength of the system-environment coupling and two-point environmental correlation functions. We validate our analytic bound numerically using the nonintegrable transverse field Ising model. Furthermore, our results provide a universal and model-agnostic quantitative framework for understanding the dynamical limits of information spreading across quantum many-body physics, condensed matter systems, and engineered quantum platforms.

Fermions↗