Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “discrete quantum measurement”

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.

52 records · Page 3

Conditioned quantum-assisted deep generative surrogate for particle-calorimeter interactions

Particle collisions at accelerators like the Large Hadron Collider (LHC), recorded by experiments such as ATLAS and CMS, enable precise standard model measurements and searches for new phenomena. Simulating these collisions significantly influences experiment design and analysis but incurs immense computational costs, projected at millions of CPU-years annually during the high luminosity LHC (HL-LHC) phase. Currently, simulating a single event with Geant4 consumes around 1000 CPU seconds, with calorimeter simulations especially demanding. To address this, we propose a conditioned quantum-assisted generative model, integrating a conditioned variational autoencoder (VAE) and a conditioned restricted Boltzmann machine (RBM). Our RBM architecture is tailored for D-Wave’s Pegasus-structured advantage quantum annealer for sampling, leveraging the flux bias for conditioning. This approach combines classical RBMs as universal approximators for discrete distributions with quantum annealing’s speed and scalability. We also introduce an adaptive method for efficiently estimating effective inverse temperature, and validate our framework on Dataset 2 of CaloChallenge.

97 MATHEMATICS AND COMPUTING↗

A Discrete Component Low-Noise Preamplifier Readout for a Linear (1x16) SiC Photodiode Array

A compact, low-noise and inexpensive preamplifier circuit has been designed and fabricated to optimally readout a common cathode (1x16) channel 4H-SiC Schottky photodiode array for use in ultraviolet experiments. The readout uses an operational amplifier with 10 pF capacitor in the feedback loop in parallel with a low leakage switch for each of the channels. This circuit configuration allows for reiterative sample, integrate and reset. A sampling technique is given to remove Johnson noise, enabling a femtoampere level readout noise performance. Commercial-off-the-shelf acquisition electronics are used to digitize the preamplifier analogue signals. The data logging acquisition electronics has a different integration circuit, which allows the bandwidth and gain to be independently adjusted. Using this readout, photoresponse measurements across the array between spectral wavelengths 200 nm and 370 nm are made to establish the array pixels external quantum efficiency, current responsivity and noise equivalent power.

low-noise preamplifier↗

Interference Control of Fano Resonances and Dynamical Imaging of an Electron Wave Packet

The Fano line shape, arising from the interference of pathways for the excitation of discrete and continuum states, plays a fundamental role in many branches of physics, chemistry, and materials science. Exciting the resonance with a high harmonic provides naturally a phase delay between the pathways leading to a complex asymmetry parameter. We demonstrate that its amplitude and phase can be controlled on the femtosecond and attosecond time scales, respectively. With our high-energy-resolution (10-meV) experiment, we dynamically image a resonance-enhanced electron wave packet during its temporal evolution, extracting both the amplitude and the phase. Calculations reproduce our experimental results. Our approach constitutes a method for measuring the photoionization delays of a resonance and enables the reconstruction of the electron wave packet in the time domain. This concept of an interference-controlled Fano line shape is a step toward attosecond quantum optics with potential ramifications into nanoscience and next-generation optical materials.

Han, Meng [ETH Zürich (Switzerland); Kansas State ↗

Investigation of boosted proton energies through proton radiography of target normal sheath acceleration fields in the multi-ps regime

Multi-kilojoule, multi-picosecond short-pulse lasers, such as the National Ignition Facility-Advanced Radiographic Capability laser and the OMEGA-Extended Performance laser, which have been constructed over the last two decades, enable exciting opportunities to produce high-brightness, high-energy laser-driven proton sources for applications in high-energy-density science like proton fast ignition for inertial fusion energy, particle radiography, and materials science studies. Results on these platforms have demonstrated enhanced accelerated proton energies and electron temperatures when compared to established scaling laws. Recent work has developed a new scaling for proton TNSA in the multi-ps regime. However, this new physics in the multi-ps regime motivates the need to understand the origin of the enhancement in proton energies. Toward this goal, here this work presents the first measurements of the TNSA accelerating sheath field in the multi-ps regime for pulse durations of 0.6, 5, and 10 ps. This measurement was achieved by using a separate TNSA proton source to radiograph the spatiotemporal profile of the accelerating sheath that is responsible for proton acceleration. The use of stacked radiochromic film detectors allows for a discrete time profile of the radiographs, thus enabling the measurement of the temporal and spatial evolution of the accelerating field. In performing this measurement, we extract quantities such as the sheath strength as a function of time and pulse duration, which shows that longer pulse durations sustain a stronger electric field for a longer duration when compared to sub-ps laser pulses, which may enable the observed boosted proton energies and proton conversion efficiencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Digitization and subduction of S U ( N ) gauge theories

The simulation of lattice gauge theories on quantum computers necessitates digitizing gauge fields. One approach involves substituting the continuous gauge group with a discrete subgroup, but the implications of this approximation still need to be clarified. To gain insights, we investigate the subduction of S U ( 2 ) and S U ( 3 ) to discrete crystal-like subgroups. Using classical lattice calculations, we show that subduction offers valuable information based on subduced direct sums, helping us identify additional terms to incorporate into the lattice action that can mitigate the effects of digitization. Furthermore, we compute the static potentials of all irreducible representations of Σ ( 360 × 3 ) at a fixed lattice spacing. Our results reveal a percent-level agreement with the Casimir scaling of S U ( 3 ) for irreducible representations that subduce to a single Σ ( 360 × 3 ) irreducible representation. This provides a diagnostic measure of approximation quality, as some irreducible representations closely match the expected results while others exhibit significant deviations. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Comparison of Integrated Radiation Transport Models with TEPC Measurements for the Average Quality Factors in Spaceflights

The purpose of this work is to test our theoretical model for the interpretation of radiation data measured in space. During the space missions astronauts are exposed to the complex field of radiation type and kinetic energies from galactic cosmic rays (GCR), trapped protons, and sometimes solar particle events (SPEs). The tissue equivalent proportional counter (TEPC) is a simple time-dependent approach for radiation monitoring for astronauts on board the International Space Station. Another and a newer approach to Microdosimetry is the use of silicon-on-insulator (SOI) technology launched on the MidSTAR-1 mission in low Earth orbit (LEO). In the radiation protection practice, the average quality factor of a radiation field is defined as a function of linear energy transfer (LET), Q(sub ave)(LET). However, TEPC measures the average quality factor as a function of the lineal energy y, Q(sub ave)(y), defined as the average energy deposition in a volume divided by the average chord length of the volume. Lineal energy, y, deviates from LET due to energy straggling, delta-ray escape or entry, and nuclear fragments produced in the detector volume. Monte Carlo track structure simulation was employed to obtain the response of a TEPC irradiated with charged particle for an equivalent site diameter of 1 micron of wall-less counter. The calculated data of the energy absorption in the wall-less counter were compiled for various y values for several ion types at various discrete projectile energy levels. For the simulation of TEPC response from the mixed radiation environments inside a spacecraft, such as, Space Shuttle and International Space Station, the complete microdosimetric TEPC response, f( y, E, Z), were calculated with the Monte Carlo theoretical results by using the first order Lagrangian interpolation for a monovariate function at a given y value (y = 0.1 keV/micron 5000 keV/micron) at any projectile energy level (E = 0.01 MeV/u to 50,000 MeV/u) of each specific radiation type (Z = 1 to 28). Because the anomalous response has been observed at large event sizes in the experiment due to the escape of energy out of sensitive volume by delta-rays and the entry of delta-rays from the high-density wall into the low-density gas-volume cavity, Monte Carlo simulation was also made for the response of a walled-TEPC with wall thickness 2 mm and density 1 g/cm(exp 3). The radius of cavity was set to 6.35 mm and a gas density 7.874 x 10(exp -5) g/cm(exp 3). The response of the walled- and the wall-less counters were compared. The average quality factor Q(sub ave)(y) for trapped protons on STS-89 demonstrated the good agreement between the model calculations and flight TEPC data as shown. Using an integrated space radiation model (this includes the transport codes HZETRN and BRYNTRN, the quantum nuclear interaction model QMSFRG) and the resultant response distribution functions of walled-TEPC from Monte-Carlo track simulations, we compared model calculations with walled-TEPC measurements from NASA missions in LEO and made predictions for the lunar and the Mars missions. The Q(sub ave)(y) values for the trapped or the solar protons ranged from 1.9-2.5. This over-estimates the Qave(LET) values which ranged from 1.4-1.6. Both quantities increase with shield thickness due to nuclear fragmentation. The Q(sub ave)(LET) for the complete GCR spectra was found to be 3.5-4.5, while flight TEPCs measured 2.9-3.4 for Q(sub ave)(y). The GCR values are decreasing with the shield thickness. Our analysis for a proper interpretation of data supports the use of TEPCs for monitoring space radiation environment.

Kim, Myung-Hee Y.↗

Development of Position-Sensitive Magnetic Calorimeters for X-Ray Astronomy

Metallic magnetic calorimeters (MMC) are one of the most promising devices to provide very high energy resolution needed for future astronomical x-ray spectroscopy. MMC detectors can be built to large detector arrays having thousands of pixels. Position-sensitive magnetic (PoSM) microcalorimeters consist of multiple absorbers thermally coupled to one magnetic micro calorimeter. Each absorber element has a different thermal coupling to the MMC, resulting in a distribution of different pulse shapes and enabling position discrimination between the absorber elements. PoSMs therefore achieve the large focal plane area with fewer number of readout channels without compromising spatial sampling. Excellent performance of PoSMs was achieved by optimizing the designs of key parameters such as the thermal conductance among the absorbers, magnetic sensor, and heat sink, as well as the absorber heat capacities. Micro fab ri - cation techniques were developed to construct four-absorber PoSMs, in which each absorber consists of a two-layer composite of bismuth and gold. The energy resolution (FWHM full width at half maximum) was measured to be better than 5 eV at 6 keV x-rays for all four absorbers. Position determination was demonstrated with pulse-shape discrimination, as well as with pulse rise time. X-ray microcalorimeters are usually designed to thermalize as quickly as possible to avoid degradation in energy resolution from position dependence to the pulse shapes. Each pixel consists of an absorber and a temperature sensor, both decoupled from the cold bath through a weak thermal link. Each pixel requires a separate readout channel; for instance, with a SQUID (superconducting quantum interference device). For future astronomy missions where thousands to millions of resolution elements are required, having an individual SQUID readout channel for each pixel becomes difficult. One route to attaining these goals is a position-sensitive detector in which a large continuous or pixilated array of x-ray absorbers shares fewer numbers of temperature sensors. A means of discriminating the signals from different absorber positions, however, needs to be built into the device for each sensor. The design concept for the device is such that the shape of the temperature pulse with time depends on the location of the absorber. This inherent position sensitivity of the signal is then analyzed to determine the location of the event precisely, effectively yielding one device with many sub-pixels. With such devices, the total number of electronic channels required to read out a given number of pixels is significantly reduced. PoSMs were developed that consist of four discrete absorbers connected to a single magnetic sensor. The design concept can be extended to more than four absorbers per sensor. The thermal conductance between the sensor and each absorber is different by design and consequently, the pulse shapes are different depending upon which absorber the xrays are received, allowing position discrimination. A magnetic sensor was used in which a paramagnetic Au:Er temperature-sensitive material is located in a weak magnetic field. Deposition of energy from an x-ray photon causes an increase in temperature, which leads to a change of magnetization of the paramagnetic sensor, which is subsequently read out using a low noise dc-SQUID. The PoSM microcalorimeters are fully microfabricated: the Au:Er sensor is located above the meander, with a thin insulation gap in between. For this position-sensitive device, four electroplated absorbers are thermally linked to the sensor via heat links of different thermal conductance. One pixel is identical to that of a single-pixel design, consisting of an overhanging absorber fabricated directly on top of the sensor. It is therefore very strongly thermally coupled to it. The three other absorbers are supported directly on a silicon-nitride membrane. These absorbers are thermally coupled to the sensor via Ti (5 nm)/Au250 nm) metal links. The strength of the links is parameterized by the number of gold squares making up the link. For detector performance, experimentally different pulse-shapes were demonstrated with 6 keV x-rays, which clearly show different rise times for different absorber positions. For energy resolution measurement, the PoSM was operated at 32 mK with an applied field that was generated using a persistent current of 50 mA. Over the four pixels, energy resolution ranges from 4.4 to 4.7 eV were demonstrated.

Bandler, SImon↗

Comparing three generations of D-Wave quantum annealers for minor embedded combinatorial optimization problems

Abstract Quantum annealing (QA) is a novel type of analog computation that aims to use quantum mechanical fluctuations to search for optimal solutions of Ising problems. QA in the transverse Ising model, implemented on D-Wave quantum processing units, are available as cloud computing resources. In this study we report concise benchmarks across three generations of D-Wave quantum annealers, consisting of four different devices, for the NP-hard discrete combinatorial optimization problems unweighted maximum clique and unweighted maximum cut on random graphs. The Ising, or equivalently quadratic unconstrained binary optimization, formulation of these problems do not require auxiliary variables for order reduction, and their overall structure and weights are not highly variable, which makes these problems simple test cases to understand the sampling capability of current D-Wave quantum annealers. All-to-all minor embeddings of size 52, with relatively uniform chain lengths, are used for a direct comparison across the Chimera, Pegasus, and Zephyr device topologies. A grid-search over annealing times and the minor embedding chain strengths is performed in order to determine the level of reasonable performance for each device and problem type. Experiment metrics that are reported are approximation ratios for non-broken chain samples, chain break proportions, and time-to-solution for the maximum clique problem instances. How fairly the quantum annealers sample optimal maximum cliques, for instances which contain multiple maximum cliques, is quantified using entropy of the measured ground state distributions. The newest generation of quantum annealing hardware, which has a Zephyr hardware connectivity, performed the best overall with respect to approximation ratios and chain break frequencies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Anomalous Loss Reduction Below Two‐Level System Saturation in Aluminum Superconducting Resonators

Superconducting resonators are widely used in many applications such as qubit readout for quantum computing, and kinetic inductance detectors. These resonators are susceptible to numerous loss and noise mechanisms, especially the dissipation due to two-level systems (TLS) which become the dominant source of loss in the few-photon and low temperature regime. In this study, capacitively-coupled aluminum half-wavelength coplanar waveguide resonators are investigated. Surprisingly, the loss of the resonators is observed to decrease with a lowering temperature at low excitation powers and temperatures below the TLS saturation. This behavior is attributed to the reduction of the TLS resonant response bandwidth with decreasing temperature and power to below the detuning between the TLS and the resonant photon frequency in a discrete ensemble of TLS. When response bandwidths of TLS are smaller than their detunings from the resonance, the resonant response and thus the loss is reduced. At higher excitation powers, the loss follows a logarithmic power dependence, consistent with predictions from the generalized tunneling model (GTM). A model combining the discrete TLS ensemble with the GTM is proposed and matches the temperature and power dependence of the measured internal loss of the resonator with reasonable parameters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Joint Acceptance Attenuation Factor of Integrated Pressure with Unsteady Pressure-Sensitive Paint Measurements

The Unsteady Pressure-Sensitive Paint (uPSP) is widely used to measure the surface pressure of the scaled model in the wind tunnel tests. Compared to the conventional pressure transducers, uPSP has the advantage of high spatial resolution. With multiple high-speed cameras, the uPSP data collected with the camera pixels are mapped to the surface grid of the scaled model and converted to pressure. The shot noise, which is associated with the quantum processes both in the generation of the photons by the luminescent uPSP and in the conversion of the photons into electrons within the camera, is the dominant component of the noise in the uPSP measurement. The integrated pressure is usually computed as an area-weighted average of the uPSP measured pressures on the grid nodes of a user-defined patch. The effect of the shot noise is reduced in the integrated pressure; however, the measurement of the aerodynamic pressure is also attenuated. This paper discusses the Joint Acceptance Attenuation Factor (JAAF) of the integrated pressure with uPSP measurements. The JAAF is defined as the ratio of the auto or cross Power Spectral Density (PSD) of the integrated aerodynamic pressure on the patch to the average auto PSD of the aerodynamic pressure on the grid nodes of the patch. Under the assumption that the cross spectrum of the aerodynamic pressure can be described with the Corcos model, the closed-form formula of the JAAF is derived for the integrated pressure on rectangular patches via the summation of uPSP measurements on discrete grid nodes. The closed-form formula is verified via the integration of the pressure field of Corcos model over continuous areas. The properties of the JAAF with different parameters are presented and the applications in the comparison of the spectrum of the pressure transducers and uPSP measurements, which were collected in the NASA Space Launch System (SLS) tests with the 11-by-11-foot transonic wind tunnel at NASA Ames Research Center in December 2017 and September 2019 respectively, are also discussed. The work described in this paper is a part of NASA’s development of a new state-of-the-art uPSP capability in production wind tunnels. Funding was provided by the NASA Aerosciences Evaluation and Test Capabilities(AETC) Portfolio Office.

Unsteady Pressure-Sensitive Paint↗

Long-Wavelength Infrared (LWIR) Quantum Dot Infrared Photodetector (QDIP) Focal Plane Array

We have exploited the artificial atomlike properties of epitaxially self-assembled quantum dots for the development of high operating temperature long wavelength infrared (LWIR) focal plane arrays. Quantum dots are nanometer-scale islands that form spontaneously on a semiconductor substrate due to lattice mismatch. QDIPs are expected to outperform quantum well infrared detectors (QWIPs) and are expected to offer significant advantages over II-VI material based focal plane arrays. QDIPs are fabricated using robust wide bandgap III-V materials which are well suited to the production of highly uniform LWIR arrays. We have used molecular beam epitaxy (MBE) technology to grow multi-layer LWIR quantum dot structures based on the InAs/InGaAs/GaAs material system. JPL is building on its significant QWIP experience and is basically building a Dot-in-the-Well (DWELL) device design by embedding InAs quantum dots in a QWIP structure. This hybrid quantum dot/quantum well device offers additional control in wavelength tuning via control of dot-size and/or quantum well sizes. In addition the quantum wells can trap electrons and aide in ground state refilling. Recent measurements have shown a 10 times higher photoconductive gain than the typical QWIP device, which indirectly confirms the lower relaxation rate of excited electrons (photon bottleneck) in QDPs. Subsequent material and device improvements have demonstrated an absorption quantum efficiency (QE) of approx. 3%. Dot-in-the-well (DWELL) QDIPs were also experimentally shown to absorb both 45 deg. and normally incident light. Thus we have employed a reflection grating structure to further enhance the quantum efficiency. JPL has demonstrated wavelength control by progressively growing material and fabricating devices structures that have continuously increased in LWIR response. The most recent devices exhibit peak responsivity out to 8.1 microns. Peak detectivity of the 8.1 micrometer devices has reached approx. 1 x 10(exp 10) Jones at 77 K. Furthermore, we have fabricated the first long-wavelength 640x512 pixels QDP focal plane array. This QDIP focal plane may has produced excellent infrared imagery with noise equivalent temperature difference of 40 mK at 60K operating temperature. In addition, we have managed to increase the quantum efficiency of these devices from 0.1% (according to the data published in literature) to 20% in discrete devices. This is a factor of 200 increase in quantum efficiency. With these excellent results, for the first time QDIP performance has surpassed the QWIP performance. Our goal is to operate these long-wavelength detectors at much higher operating temperature than 77K which can be passively achieved in space. This will be a huge leap in high performance infrared detectors specifically applicable to space science instruments.

infrared imaging↗

Non-Markovian relaxation spectroscopy of fluxonium qubits

Recent studies have shown that parasitic two-level systems (TLS) in superconducting qubits, which are a leading source of decoherence, can have relaxation times longer than the qubits themselves. However, the standard techniques used to characterize qubit relaxation is only valid for measuring T 1 under the Born-Markov approximation and could mask environmental memory effects in practice. Here, we introduce two-timescale relaxometry, a technique to probe the qubit and environment relaxation simultaneously and efficiently. We apply it to high-coherence fluxonium qubits over a frequency range of 0.1-0.4 GHz, and reveal a discrete spectrum of TLS with millisecond lifetimes. Our analysis of the spectrum is consistent with a random distribution of TLS in the aluminum oxide tunnel barrier of the Josephson junction chain of the fluxonium, with a spectral and volumetric density and average electric dipole similar to previous TLS studies at much higher frequencies. Our study suggests that investigating and mitigating TLS in the junction chain is crucial to the development of various types of noise-protected qubits in circuit QED.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Staggered scheme for the compressible fluctuating hydrodynamics of multispecies fluid mixtures

Here, we present a numerical formulation for the solution of nonisothermal, compressible Navier-Stokes equations with thermal fluctuations to describe mesoscale transport phenomena in multispecies fluid mixtures. The novelty of our numerical method is the use of staggered grid momenta along with a finite volume discretization of the thermodynamic variables to solve the resulting stochastic partial differential equations. The key advantages of the numerical scheme are that it significantly simplifies the discretization of diffusive and stochastic momentum fluxes into a more compact form, and it provides an unambiguous prescription of boundary conditions involving pressure. The staggered grid scheme more accurately reproduces the equilibrium static structure factor of hydrodynamic fluctuations in gas mixtures compared to a collocated scheme described previously by Balakrishnan et al. [Phys. Rev. E 89, 013017 (2014)1539-375510.1103/PhysRevE.89.013017]. The numerical method is tested for ideal noble gases mixtures under various nonequilibrium conditions, such as applied thermal and concentration gradients, to assess the role of cross-diffusion effects, such as Soret and Dufour, on the long-ranged correlations of hydrodynamic fluctuations, which are also more accurately reproduced compared to the collocated scheme. We numerically study giant nonequilibrium fluctuations driven by concentration gradients and fluctuation-driven Rayleigh-Taylor instability in gas mixtures. Wherever applicable, excellent agreement is observed with theory and measurements from the direct simulation Monte Carlo method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Accurate numerical simulations of open quantum systems using spectral tensor trains

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

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ambiguities in the partial-wave analysis of the photoproduction of pairs of pseudoscalar mesons

Applying the technique of partial-wave analysis, there are cases where more than one set of underlying complex-valued amplitudes can describe the measured observables. These ambiguities can sometimes be resolved using additional information, but assumptions are often required. It is known that the partial-wave analysis of two-pseudoscalar meson systems produced in photoproduction with a linearly polarized photon beam is free from discrete ambiguities stemming from the Barrelet zeros when the nucleon spin is ignored. In this article, we show that continuous ambiguities are possible for certain wave sets, even though the discrete ambiguities do not appear. We also explore ways to resolve these ambiguities and determine the maximal amount of information that can be obtained from analyses that suffer from these continuous ambiguities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Perturbation theory under the truncated Wigner approximation: How system-environment entanglement formation drives quantum decoherence

In discrete quantum systems, quantum decoherence is the disappearance of simple phase relations as a result of interactions with an environment. For many applications, the question is not necessarily how to avoid (inevitable) system-environment interactions, but rather how to design environments that optimally preserve a system's phase relations in spite of such interactions. The formation of system-environment entanglement is a major driving mechanism for decoherence, and a detailed understanding of this process could inform strategies for conserving coherence optimally. This requires scalable, flexible, and systematically improvable quantum dynamical methods that retain detailed information about the entanglement properties of the environment, yet very few current methods offer this combination of features. Here, in this work, we address this need by introducing a theoretical framework wherein we combine the truncated Wigner approximation with standard time-dependent perturbation theory, allowing for computing expectation values of operators in the combined system-environment Hilbert space. We demonstrate the utility of this framework by applying it to the spin-boson model, representative of qubits and simple donor-acceptor systems. For this model, our framework provides an analytical description of perturbative contributions to expectation values. We monitor how quantum decoherence at zero temperature is accompanied by entanglement formation with individual environmental degrees of freedom. Based on this entanglement behavior, we find that the selective suppression of low-frequency environmental modes is particularly effective for mitigating quantum decoherence.

74 ATOMIC AND MOLECULAR PHYSICS↗