SEARCH · Search NASA

Results for “quantum state estimation”

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 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Real classical shadows

Efficiently learning expectation values of a quantum state using classical shadow tomography has become a fundamental task in quantum information theory. In a classical shadows protocol, one measures a state in a chosen basis $\mathcal{W}$ after it has evolved under a unitary transformation randomly sampled from a chosen distribution $\mathcal{U}$. In this work we study the case where $\mathcal{U}$ corresponds to either local or global orthogonal Clifford gates, and $\mathcal{W}$ consists of real-valued vectors. Our results show that for various situations of interest, this ‘real’ classical shadow protocol improves the sample complexity over the standard scheme based on general Clifford unitaries. For example, when one is interested in estimating the expectation values of arbitrary real-valued observables, global orthogonal Cliffords typically decrease the required number of samples by a factor of two. More dramatically, for k-local observables composed only of real-valued Pauli operators, sampling local orthogonal Cliffords leads to a reduction by an exponential-in-k factor in the sample complexity over local unitary Cliffords. Finally, we show that by measuring in a basis containing complex-valued vectors, orthogonal shadows can, in the limit of large system size, exactly reproduce the original unitary shadows protocol.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Assessing the difficulty of capturing the distribution function of neutrinos in neutron star merger simulations

The collision of two neutron stars is a rich source of information about nuclear physics. In particular, the kilonova signal following a merger can help us elucidate the role of neutron stars in nucleosynthesis, and informs us about the properties of matter above nuclear saturation. Approximate modeling of neutrinos remains an important limitation to our ability to make predictions for these observables. Part of the problem is the fermionic nature of neutrinos. By the exclusion principle, the expected value 𝑓 𝜈 for the number of neutrinos in a quantum state is at most 1. Any process producing neutrinos is suppressed by a blocking factor (1 −𝑓 𝜈 ). Recent simulations focused on neutrino physics mostly use a gray two-moment scheme to evolve neutrinos. This evolves integrals of 𝑓 𝜈 over momentum space, preventing direct calculations of blocking factors. Monte Carlo methods may be an attractive alternative, providing access to the full distribution of neutrinos. Their current implementation is, however, inadequate to estimate 𝑓 𝜈 : in our most recent simulations, a single Monte Carlo packet causes, in the worst cases, estimates of 𝑓 𝜈 to jump from 𝑓 𝜈 =0 to 𝑓 𝜈 ∼10 5 . While this is concerning, this brazen violation of the fermionic nature of neutrinos has been largely inconsequential, as the interactions used in simulations avoid direct calculations of 𝑓 𝜈 . We are, however, reaching a level of modeling at which this problem can no longer be ignored. Here, we discuss the relatively simple origin of this issue. We then show that very rough estimates of 𝑓 𝜈 can in theory be obtained in merger simulations, but that they will require a combination of unintuitive weighting schemes for Monte Carlo packets and smoothing of the neutrino distribution at coarser resolution than what the merger simulation uses.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Nanoscale Defects as Probes of Time-Reversal Symmetry Breaking

Nanoscale defects such as nitrogen-vacancy (NV) centers can serve as sensitive and noninvasive probes of electromagnetic fields and fluctuations from materials, which in turn can be used to characterize these systems. Here we specifically discuss how NV centers can probe time-reversal symmetry breaking (TRSB) phenomena in low-dimensional electronic systems. We argue that the difference in relaxation rates Γ$_{±\hat{z}}$ of NV centers starting from 𝑚 = ±1 spin states to the ground state with 𝑚 = 0 directly probes TRSB. The effect arises from the difference in the fluctuation spectrum of left- and right-polarized electromagnetic fields emanating from such materials. In the quantum Hall setting, the NV center experiences (nearly zero) large additional contribution to its relaxation due to the presence of the material when its magnetic dipole (anti)aligns with the external field. More generally, the difference in the relaxation rates is sensitive to the imaginary part of the wave-vector-dependent Hall conductivity. We argue that this can be used to determine the Hall viscosity, which can potentially distinguish candidate fractional quantum Hall states and be used to infer the pairing angular momentum in TRSB superconductors. We compute the average relaxation rate [Γ$_{+\hat{z}}$ + Γ$_{-\hat{z}}$] near thin film superconductors and find that it exhibits a Hebel-Slichter-like enhancement below 𝑇 𝑐 . The difference Γ$_{+\hat{z}}$ −Γ$_{-\hat{z}}$ also inherits this peak but is only nonzero for 𝑇 < 𝑇 𝑐 and only if the superconductivity is chiral. We provide concrete estimates for observing this effect in stacked twisted bismuth strontium calcium copper oxide flakes.

NV centers↗

DOE OSTI2025

Quantum chromodynamics resolution of the ATOMKI anomaly in 4 He nuclear transitions

Observations of anomalous angular correlations in electron-positron pairs produced from excited states of 4 He, 8 Be, and 12 C nuclei have been suggested as being due to the creation and subsequent decay of a new light particle of mass ~17MeV. In this work, we investigate the possibility that the source of the observed signals is a set of new excitation channels created by the 12-quark hidden-color Fock state within the 4 He nuclear wave function dubbed the “hexadiquark.” We calculate the invariant e + ⁢ e - mass spectrum for the electromagnetic transition from a new excitation of 4 He, estimating its differential and total decay width. We find that we can fit the shape of the anomalous signal with the QCD Fock state at excitation energy E* = 17.9 ± 1MeV and a Gaussian form factor for the electromagnetic decay. We address the physical issues with the fit parameters using properties of the hexadiquark state, in particular the three weakly repulsive 6 C interactions of SU⁡(3) C between diquark pairs. Experimental tests of our model are described in detail. In light of this work, we emphasize the need for independent experimental confirmation or refutation of the ATOMKI results as well as dedicated experiments to search for the proposed new excitations of 4 He and other a-cluster nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Characterizing non-Markovian and coherent errors in quantum simulation

Quantum simulation of many-body systems, particularly using ultracold atoms and trapped ions, presents a unique form of quantum control—it is a direct implementation of a multi-qubit gate generated by the Hamiltonian. As a consequence, it also faces a unique challenge in terms of benchmarking, because the well-established gate benchmarking techniques are unsuitable for this form of quantum control. Here we show that the symmetries of the target many-body Hamiltonian can be used not only to benchmark but to characterize experimental errors in the quantum simulation. We use our results to develop protocols to characterize these errors, which can be implemented using state-of-the-art technology. We consider two forms of errors: (i) unitary errors arising out of systematic errors in the applied Hamiltonian and (ii) canonical non-Markovian errors arising out of random shot-to-shot fluctuations in the applied Hamiltonian. We show that the dynamics of the expectation value of the target Hamiltonian itself, which is ideally constant in time, can be used to characterize these errors. In the presence of errors, the expectation value of the target Hamiltonian shows a characteristic thermalization dynamics, when it satisfies the operator thermalization hypothesis (OTH). That is, an oscillation in the short time followed by relaxation to a steady-state value in the long time limit. We show that while the steady-state value can be used to characterize the coherent errors, the amplitude of the oscillations can be used to estimate the non-Markovian errors. We prove a sandwich theorem to establish a linear relation between the amplitude of the oscillations and the magnitude of the non-Markovian errors. Moreover, by varying the initial state, we show that the steady state values can be used to completely construct the generator of the coherent errors. Using these results, we develop two experimental protocols to characterize the unitary errors based on these results, one of which requires single-qubit addressing and the other one doesn't. We also develop a protocol to characterize non-Markovian errors. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1987

Absolute I(asterisk) quantum yields for the ICN A state by diode laser gain-vs-absorption spectroscopy

Absolute I(asterisk) quantum yields have been measured as a function of wavelength for room temperature photodissociation of the ICN A state continuum. The yields are obtained by the technique of time-resolved diode laser gain-vs-absorption spectroscopy. Quantum yields are evaluated at seven wavelengths from 248 to 284 nm. The yield at 266 nm is 66.0 + or - 2 percent and it falls off to 53.4 + or - 2 percent and 44.0 + or - 4 percent at 284 and 248 nm, respectively. The latter values are significantly higher than those obtained by previous workers using infrared fluorescence. Estimates of I(asterisk) quantum yields obtained from analysis of CN photofragment rotational distributions, as discussed by other workers, are in good agreement with the I(asterisk) yields reported here. The results are considered in conjunction with recent theoretical and experimental work on the CN rotational distributions and with previous I(asterisk) quantum yield results.

Hess, Wayne P.↗

NASA NTRS1991

The energy separation between the classical and nonclassical isomers of protonated acetylene - An extensive study in one- and n-particle space saturation

The energy separation between the classical and nonclassical forms of protonated acetylene has been reinvestigated in light of the recent experimentally deduced lower bound to this value of 6.0 kcal/mol. The objective of the present study is to use state-of-the-art ab initio quantum mechanical methods to establish this energy difference to within chemical accuracy (i.e., about 1 kcal/mol). The one-particle basis sets include up to g-type functions and the electron correlation methods include single and double excitation coupled-cluster (CCSD), the CCSD(T) extension, multireference configuration interaction, and the averaged coupled-pair functional methods. A correction for zero-point vibrational energies has also been included, yielding a best estimate for the energy difference between the classical and nonclassical forms of 3.7 + or - 1.3 kcal/mol.

Lindh, Roland↗

DOE OSTI2026

Quantum utility-scale error mitigation for quantum quench dynamics in Heisenberg spin chains

Here, we implement a quantum error mitigation method termed self-mitigation, which is comparable to zero-noise extrapolation, at large scales to achieve quantum utility on near-term, noisy quantum computers. We investigate the effectiveness of several quantum error mitigation strategies, including self-mitigation, by simulating quantum quench dynamics for Heisenberg spin chains with system sizes up to 104 qubits using IBM quantum processors. In particular, we discuss the limitations of zero-noise extrapolation and the advantages offered by self-mitigation at large scales. The self-mitigation method demonstrates stable accuracy with large systems of 104 qubits comprising more than 3,000 CNOT gates. Also, we combine the discussed quantum error mitigation methods with practical entanglement entropy measuring methods, and it shows a good agreement with the theoretical estimation. Our study illustrates the usefulness of near-term noisy quantum hardware in examining the quantum quench dynamics of many-body systems at large scales and lays the groundwork for surpassing classical simulations with quantum methods prior to the development of fault-tolerant quantum computers.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Surrogate optimization of variational quantum circuits

Variational quantum eigensolvers are touted as a near-term algorithm capable of impacting many applications. However, the potential has not yet been realized, with few claims of quantum advantage and high resource estimates, especially due to the need for optimization in the presence of noise. Finding algorithms and methods to improve convergence is important to accelerate the capabilities of near-term hardware for VQE or more broad applications of hybrid methods in which optimization is required. To this goal, we look to use modern approaches developed in circuit simulations and stochastic classical optimization, which can be combined to form a surrogate optimization approach to quantum circuits. Using an approximate (classical CPU/GPU) state vector simulator as a surrogate model, we efficiently calculate an approximate Hessian, passed as an input for a quantum processing unit or exact circuit simulator. This method will lend itself well to parallelization across quantum processing units. We demonstrate the capabilities of such an approach with and without sampling noise and a proof-of-principle demonstration on a quantum processing unit utilizing 40 qubits.

Gustafson, Erik J. [RIACS, Mtn. View] (ORCID:00000↗

DOE OSTI2024

Identifiability and characterization of transmon qutrits through Bayesian experimental design

Robust control of a quantum system is essential to utilize the current noisy quantum hardware to its full potential, such as quantum algorithms. To achieve such a goal, a systematic search for an optimal control for any given experiment is essential. The design of optimal control pulses requires accurate numerical models and, therefore, accurate characterization of the system parameters. We present an online Bayesian approach for quantum characterization of qutrit systems, which automatically and systematically identifies optimal experiments that provide maximum information on the system parameters, thereby greatly reducing the number of experiments that need to be performed on the quantum testbed. Unlike most characterization protocols that provide point-estimates of the parameters, the proposed approach is able to estimate their probability distribution. The applicability of the Bayesian experimental design technique was demonstrated on test problems, where each experiment was defined by a parameterized control pulse. In addition to this, we also present an approach for iterative pulse extension, which is robust under uncertainties in transition frequencies and coherence times, and shot noise, despite being initialized with wide uninformative priors. Furthermore, we provide a mathematical proof of the theoretical identifiability of the model parameters and present conditions on the quantum state under which the parameters are identifiable. The proof and conditions for identifiability are presented for both closed and open quantum systems using the Schrödinger equation and the Lindblad master equation, respectively.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2020

Remote Sensing of Diatom Bloom Succession

Marine diatoms are major biogeochemical and ecological “influencers” that contribute to a large fraction of the carbon export and supplying fisheries (Falkowski 2015). The fluxes of carbon transfer to the food web or to the deep ocean vary according to the stage of a diatom bloom (Du Toit 2018). Stages can be determined using inherent optical properties that reflect their physiological state, such as the chlorophyll fluorescence to particulate backscattering ratio (ChlF/b(sub bp), Cetinic et al. 2015). Identifying the bloom stage can potentially improve biogeochemical models of carbon export and fishery management. However, it is not yet possible to adequately determine the stage of phytoplankton blooms using satellites. Satellite-derived remote sensing reflectance R(sub rs)(𝜆) allow for remote identification of diatom blooms in the open ocean (Sathyendranath et al. 2004), and there are techniques to estimate the fluorescence quantum yield (𝜑) that, when high, can indicate the nutrient limitation that often takes place when blooms start to senesce (Behrenfeld et al. 2009). The goal of this study is to use the ratio between the normalized fluorescence line height from R(sub rs)(𝜆) (nFLH) and the particulate backscattering (b(sub bp)(443)) provided by satellites to identify exponentially growing and senescent diatom blooms from space.

Lange, Priscila K.↗

DOE OSTI2025

Lanczos algorithm for lattice QCD matrix elements

Recent work [M. L. Wagman, Lanczos, the transfer matrix, and the signal-to-noise problem, .] found that an analysis formalism based on the Lanczos algorithm allows energy levels to be extracted from Euclidean correlation functions with faster ground-state convergence than effective masses, convergent estimators for multiple states from a single correlator, and two-sided error bounds. After filtering out spurious eigenvalues and using outlier-robust estimators within a nested bootstrap framework, Lanczos estimators behave more like multistate fit results than effective masses—but without involving statistical fitting. We extend this formalism to the determination of matrix elements from three-point correlation functions and provide a physical picture of “spurious-state filtering” involving restriction to a Hermitian subspace. We demonstrate similar advantages for matrix elements as for spectroscopy through example applications to noiseless mock-data and (bare) forward matrix elements of the strange scalar current between both ground and excited states with the quantum numbers of the nucleon.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Two-Tower Quantum Matrix Chain Multiplication: Trading Qubits for Depth

Matrix chain multiplication -- computing $\mathcal{W} = M^{(0)}\cdots M^{(K-1)}$ where $M^{(k)} \in \mathbb{R}^{P_k \times P_{k+1}}$-- arises in scientific computing, machine learning, and graph analysis. Despite the importance of this problem, for chains of distinct matrices, the classical number of operations grows linearly with the chain length $K$ and polynomially in the matrix dimensions. We present \emph{Two-Tower Matrix Multiplication}, a quantum subroutine that encodes the product $\mathcal{W}$ of the $K$ matrices into a quantum state in circuit depth $\mathcal{O}(\max_{k} \mathrm{polylog} (P_k P_{k+1}))$, which is independent of~$K$ within the QRAM-based state-preparation model, whereas the qubit count is $\mathcal{O}\bigl(\sum_{k} \log P_k \bigr)$; the total gate count remains linear in $K$, so the gain is in the circuit depth. The construction interleaves state-preparation operators across two layers; within each layer, all operators act on disjoint registers and execute in parallel. This subroutine can be specialized for the chain-vector case, which computes the product of $K-1$ matrices applied to a vector. We prove the correctness of the subroutine for all $K$ and provide two implementations using the Qiskit and QCLAB frameworks. The subroutine is applicable to any downstream quantum algorithm that operates on a matrix encoded in the statevector, including norm estimation, graph-matrix powers, linear system solving, and quantum machine learning kernels.

Antonioli, Giacomo [Pisa U.] (ORCID:00090000668703↗

DOE OSTI2025

Studying QED 3 with radial quantization on the lattice: Free limit

To investigate the three-dimensional quantum electrodynamics in the radial quantization on the lattice, the lattice action is constructed and the free limit is studied on 𝑆 2 × ℝ. With the overlap fermion, it is numerically verified that the important symmetries of the theory can be realized on the lattice. The analytic correlators are derived and compared to the lattice results, which agree including the overall normalization. The 𝑂⁡(𝑎 2 ) scaling is confirmed toward the analytic value in the continuum limit, and the number of reproduced excited states is estimated heuristically for the first few refinement levels. Our study helps us identify the features of the theory that we can study on the icosahedral lattice without fine-tuning.

Boyle, Peter A. [Brookhaven National Laboratory (B↗

DOE OSTI2025

Exponentially Reduced Circuit Depths Using Trotter Error Mitigation

Product formulas are a popular class of digital quantum simulation algorithms due to their conceptual simplicity, low overhead, and performance, which often exceeds theoretical expectations. Recently, Richardson extrapolation and polynomial interpolation have been proposed to mitigate the Trotter error incurred by the use of these formulas. This work provides a rigorous, general analysis of these techniques for computing time-evolved observables, simplifying the interpolation algorithm in the process, and shows that extrapolation generically improves the performance of product formulas for this task. We demonstrate that, to achieve error 𝜖 in a simulation of time 𝑇 using a 𝑝 ⁢th-order product formula with extrapolation, circuit depths of 𝑂⁡(𝑇 1+1/𝑝 ⁢polylog (1/𝜖)) are sufficient—an exponential improvement in the precision over product formulas alone. Furthermore, we prove that these algorithms achieve commutator scaling, and improve the 𝑇 complexity for the interpolation algorithm. By relaxing the requirement of performing exact Chebyshev interpolation, our simplified algorithm eliminates the need for fractional implementations of Trotter steps, reducing computational overhead. Finally, we show these techniques can be combined with the classical shadows method to estimate many time-evolved local observables. Taken together, our findings provide the strongest evidence yet for the utility of Trotter error-mitigation techniques in algorithmic applications.

quantum algorithms & computation↗

NASA NTRS1998

Quantum Computation Using Optically Coupled Quantum Dot Arrays

A solid state model for quantum computation has potential advantages in terms of the ease of fabrication, characterization, and integration. The fundamental requirements for a quantum computer involve the realization of basic processing units (qubits), and a scheme for controlled switching and coupling among the qubits, which enables one to perform controlled operations on qubits. We propose a model for quantum computation based on optically coupled quantum dot arrays, which is computationally similar to the atomic model proposed by Cirac and Zoller. In this model, individual qubits are comprised of two coupled quantum dots, and an array of these basic units is placed in an optical cavity. Switching among the states of the individual units is done by controlled laser pulses via near field interaction using the NSOM technology. Controlled rotations involving two or more qubits are performed via common cavity mode photon. We have calculated critical times, including the spontaneous emission and switching times, and show that they are comparable to the best times projected for other proposed models of quantum computation. We have also shown the feasibility of accessing individual quantum dots using the NSOM technology by calculating the photon density at the tip, and estimating the power necessary to perform the basic controlled operations. We are currently in the process of estimating the decoherence times for this system; however, we have formulated initial arguments which seem to indicate that the decoherence times will be comparable, if not longer, than many other proposed models.

Pradhan, Prabhakar↗

DOE OSTI2026

Ab Initio Many Body Quantum Embedding and Local Correlation in Crystalline Materials using Interpolative Separable Density Fitting

We present an efficient implementation of ab initio many-body quantum embedding and local correlation methods for infinite periodic systems through translational symmetry adapted interpolative separable density fitting, an approach which reduces the scaling of the calculations to only linear with the number of k-points. Employing this methodology, we compute correlated ground-state coupled cluster energies within density matrix embedding and local natural orbital correlation frameworks for both weakly and strongly correlated solids, using up to 1000 k-points. By extrapolating the local correlation domains and k-point sampling we further obtain estimates of the full coupled cluster with singles, doubles, and perturbative triples ground-state energies in the thermodynamic limit.

Chemical Physics (physics.chem-ph)↗

NASA NTRS1988

Temperature dependence of the atmospheric photolysis rate coefficient for NO2

Accurate values for the photolysis rate coefficient of NO2 (j1) are required for studies related to the observed imbalance in the photostationary state of O3, NO, and NO2 in the troposphere. Direct measurements of the temperature dependence of j1 at temperatures from -70 to 30 C were made in sunlight for relatively cloudless summer days in Boulder, Colorado. The ratios of j1 (30 C)/j1 (T C) for T = -10 C and -70 C, respectively, were 1.046 + or - 0.040 and 1.070 + or - 0.031. The j1 ratios were independent of solar zenith angle. Theoretical estimates of the temperature-dependent j1 ratios based upon recently reported cross section (sigma) and quantum yield (phi) data are more consistent with these experimental measurements than those based upon the currently accepted sigma and phi data.

Shetter, Richard E.↗