Search NASA⌕ Search

SEARCH · Search NASA

Results for “Structural error”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,117 records · Page 62

Equation of state of boron nitride combining computation, modeling, and experiment

The equation of state (EOS) of materials at warm dense conditions poses significant challenges to both theory and experiment. We report a combined computational, modeling, and experimental investigation leveraging new theoretical and experimental capabilities to investigate warm-dense boron nitride (BN). The simulation methodologies include path integral Monte Carlo (PIMC), several density functional theory (DFT) molecular dynamics methods [plane-wave pseudopotential, Fermi operator expansion (FOE), and spectral quadrature (SQ)], activity expansion (actex), and all-electron Green's function Korringa-Kohn-Rostoker (mecca), and compute the pressure and internal energy of BN over a broad range of densities and temperatures. Our experiments were conducted at the Omega laser facility and the Hugoniot response of BN to unprecedented pressures (1200–2650 GPa). The EOSs computed using different methods cross validate one another in the warm-dense matter regime, and the experimental Hugoniot data are in good agreement with our theoretical predictions. By comparing the EOS results from different methods, we assess that the largest discrepancies between theoretical predictions are ≲ 4% in pressure and ≲ 3% in energy and occur at 10 6 K, slightly below the peak compression that corresponds to the K -shell ionization regime. At these conditions, we find remarkable consistency between the EOS from DFT calculations performed on different platforms and using different exchange-correlation functionals and those from PIMC using free-particle nodes. This provides strong evidence for the accuracy of both PIMC and DFT in the high-pressure, high-temperature regime. Moreover, the recently developed SQ and FOE methods produce EOS data that have significantly smaller statistical error bars than PIMC, and so represent significant advances for efficient computation at high temperatures. The shock Hugoniot predicted by PIMC, actex, and mecca shows a maximum compression ratio of 4.55 ± 0.05 for an initial density of 2.26 g/cm 3 , higher than the Thomas-Fermi predictions by about 5%. In addition, we construct tabular EOS models that are consistent with the first-principles simulations and the experimental data. Our findings clarify the ionic and electronic structure of BN over a broad range of temperatures and densities and quantify their roles in the EOS and properties of this material. The tabular models may be utilized for future simulations of laser-driven experiments that include BN as a candidate ablator material.

36 MATERIALS SCIENCE↗

Performance of the strongly constrained and appropriately normed density functional for solid-state materials

Constructed to satisfy 17 known exact constraints for a semilocal density functional, the strongly constrained and appropriately normed (SCAN) meta-generalized-gradient-approximation functional has shown early promise for accurately describing the electronic structure of molecules and solids. One open question is how well SCAN predicts the formation energy, a key quantity for describing the thermodynamic stability of solid-state compounds. To answer this question, we perform an extensive benchmark of SCAN by computing the formation energies for a diverse group of nearly 1000 crystalline compounds for which experimental values are known. Due to an enhanced exchange interaction in the covalent bonding regime, SCAN substantially decreases the formation energy errors for strongly bound compounds, by approximately 50% to 110 meV/atom, as compared to the generalized gradient approximation of Perdew, Burke, and Ernzerhof (PBE). However, for intermetallic compounds, SCAN performs moderately worse than PBE with an increase in formation energy error of approximately 20%, stemming from SCAN's distinct behavior in the weak bonding regime. The formation energy errors can be further reduced via elemental chemical potential fitting. We find that SCAN leads to significantly more accurate predicted crystal volumes, moderately enhanced magnetism, and mildly improved band gaps as compared to PBE. Altogether, SCAN represents a significant improvement in accurately describing the thermodynamics of strongly bound compounds.

36 MATERIALS SCIENCE↗

An Adaptive Virtual Oscillator Control Structure for Grid-Forming Inverters

The electrical grid is facing unprecedented challenges due to the increasing penetration of inverter-based resources. Grid forming inverters (GFMIs) are a promising technology to address these challenges. Recently, the virtual oscillator based GFMI control is attracting more attention due to its superior performance over other control strategies. In this article, an adaptive control strategy is proposed to provide flexible operation and transition between grid-connected and islanded modes. The controller adapts the virtual oscillator's parameter values depending on the operation mode. It also provides a feedback signal using a measured frequency to account for any steady-state errors and to allow a seamless transition from grid-connected to islanded mode. Finally, to show the feasibility of the proposed controller, this article discusses the simulation results from the implementation of the controller on a single inverter system and on a group of inverters on a large practical utility feeder, the IEEE 13 node feeder, using the DIgSILENT simulation environment.

adaptive control↗

A Low-Rank Multigrid Method for the Stochastic Steady-State Diffusion Problem

We study a multigrid method for solving large linear systems of equations with tensor product structure here. Such systems are obtained from stochastic finite element discretization of stochastic partial differential equations such as the steady-state diffusion problem with random coefficients. When the variance in the problem is not too large, the solution can be well approximated by a low-rank object. In the proposed multigrid algorithm, the matrix iterates are truncated to low rank to reduce memory requirements and computational effort. The technique is proved convergent with an analytic error bound. Numerical experiments show its effectiveness in solving the Galerkin systems compared to the original multigrid solver, especially when the number of degrees of freedom associated with the spatial discretization is large.

97 MATHEMATICS AND COMPUTING↗

Preconditioning of Interplanetary Space Due to Transient CME Disturbances

Interplanetary space is characteristically structured mainly by high-speed solar wind streams emanating from coronal holes and transient disturbances such as coronal mass ejections (CMEs). While high-speed solar wind streams pose a continuous outflow, CMEs abruptly disrupt the rather steady structure, causing large deviations from the quiet solar wind conditions. For the first time, we give a quantification of the duration of disturbed conditions (preconditioning) for interplanetary space caused by CMEs. To this aim, we investigate the plasma speed component of the solar wind and the impact of in situ detected interplanetary CMEs (ICMEs), compared to different background solar wind models (ESWF, WSA, persistence model) for the time range 2011–2015. We quantify in terms of standard error measures the deviations between modeled background solar wind speed and observed solar wind speed. Using the mean absolute error, we obtain an average deviation for quiet solar activity within a range of 75.1–83.1 km s{sup −1}. Compared to this baseline level, periods within the ICME interval showed an increase of 18%–32% above the expected background, and the period of two days after the ICME displayed an increase of 9%–24%. We obtain a total duration of enhanced deviations over about three and up to six days after the ICME start, which is much longer than the average duration of an ICME disturbance itself (∼1.3 days), concluding that interplanetary space needs ∼2–5 days to recover from the impact of ICMEs. The obtained results have strong implications for studying CME propagation behavior and also for space weather forecasting.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantified Gamow shell model interaction for p s d -shell nuclei

Background: The structure of weakly bound and unbound nuclei close to particle drip lines is one of the major science drivers of nuclear physics. A comprehensive understanding of these systems goes beyond the traditional configuration interaction approach formulated in the Hilbert space of localized states (nuclear shell model) and requires an open quantum system description. The complex-energy Gamow shell model (GSM) provides such a framework as it is capable of describing resonant and nonresonant many-body states on equal footing. Purpose: To make reliable predictions, quality input is needed that allows for the full uncertainty quantification of theoretical results. In this study, we carry out the optimization of an effective GSM (one-body and two-body) interaction in the psdf-shell-model space. The resulting interaction is expected to describe nuclei with 5 ≤ A ≲ 12 at the p-sd-shell interface. Method: The one-body potential of the 4 He core is modeled by a Woods-Saxon + spin-orbit + Coulomb potential, and the finite-range nucleon-nucleon interaction between the valence nucleons consists of central, spin-orbit, tensor, and Coulomb terms. The GSM is used to compute key fit observables. The χ 2 optimization is performed using the Gauss-Newton algorithm augmented by the singular value decomposition technique. The resulting covariance matrix enables quantification of statistical errors within the linear regression approach. Results: The optimized one-body potential reproduces nucleon- 4 He scattering phase shifts up to an excitation energy of 20 MeV. The two-body interaction built on top of the optimized one-body field is adjusted to the bound and unbound ground-state binding energies and selected excited states of the helium, lithium, and beryllium isotopes up to A = 9 . A very good agreement with experimental results was obtained for binding energies. First applications of the optimized interaction include predictions for two-nucleon correlation densities and excitation spectra of light nuclei with quantified uncertainties. In conclusion: The new interaction will enable comprehensive and fully quantified studies of structure and reactions aspects of nuclei from the psd region of the nuclear chart.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Accelerating Embedding Potential Optimization by Reconstructing the Pseudo-Valence Electron Density

Density functional embedding theory (DFET) enables use of electronic structure methods with higher accuracy than density functional theory in a local region, with applications thus far ranging from (photo/electro)catalysis to reactions in solution. DFET partitions a large collection of atoms into smaller groups that interact via a shared embedding (interaction) potential V emb , determined via functional optimization. The optimized effective potential (OEP) process used to optimize V emb is time-consuming and becomes a computational bottleneck due to sharp, oscillating features of V emb near nuclei. Here, similar to pseudopotential theory, by reconstructing electron densities used in the OEP process from smoother pseudo-valence-only (PVO) electron densities as proxies for total densities of the full system and subsystems, we can retain accuracy in the embedded electronic structure calculations while potentially reducing the overhead of V emb construction, within the projector augmented-wave (PAW) formalism. We explore three different chemical reactions as exemplars to test PVO–DFET, namely, H 2 dissociative adsorption on a Cu(111) surface, H 2 O adsorption on a Pt(111) surface, and aqueous [Ca 2+ –SO 4 2– ] ion-pair formation. The PVO approximation works well for all three systems with minimal loss of accuracy (∼10–70 meV error relative to the original exact-derivative (ED) approach) while accelerating V emb generation for the Cu and Pt systems respectively by 20× and 5×. Given proper numerical convergence parameters, the spatial distributions of differences between PVO- and ED-based V emb outside the core regions are small, explaining the exceptional agreement between the two approaches. Finally, we anticipate that this more efficient PVO–DFET approximation will be useful whenever computation of V emb is much more expensive than subsequent embedded high-level electron correlation calculations.

approximation↗

Demonstration of an AI-driven workflow for dynamic x-ray spectroscopy

X-ray absorption near edge structure (XANES) spectroscopy is a powerful technique for characterizing the chemical state and symmetry of individual elements within materials, but requires collecting data at many energy points which can be time-consuming. While adaptive sampling methods exist for efficiently collecting spectroscopic data, they often lack domain-specific knowledge about the structure of XANES spectra. Here we demonstrate a knowledge-injected Bayesian optimization approach for adaptive XANES data collection that incorporates understanding of spectral features like absorption edges and pre-edge peaks. We show this method accurately reconstructs the absorption edge of XANES spectra using only 15–20% of the measurement points typically needed for conventional sampling, while maintaining the ability to determine the x-ray energy of the sharp peak after the absorption edge with errors less than 0.03 eV, the absorption edge with errors less than 0.1 eV; and overall root-mean-square errors less than 0.005 compared to traditionally sampled spectra. Our experiments on battery materials and catalysts demonstrate the method’s effectiveness for both static and dynamic XANES measurements, improving data collection efficiency and enabling better time resolution for tracking chemical changes. This approach advances the degree of automation in XANES experiments, reducing the common errors of under- or over-sampling points near the absorption edge and enabling dynamic experiments that require high temporal resolution or limited measurement time.

Bayesian optimization↗

BEYONDPLANCK IV. On end-to-end simulations in CMB analysis — Bayesian versus frequentist statistics

End-to-end simulations play a key role in the analysis of any high-sensitivity cosmic microwave background (CMB) experiment, providing high-fidelity systematic error propagation capabilities that are unmatched by any other means. In this paper, we address an important issue regarding such simulations, namely, how to define the inputs in terms of sky model and instrument parameters. These may either be taken as a constrained realization derived from the data or as a random realization independent from the data. We refer to these as posterior and prior simulations, respectively. We show that the two options lead to significantly different correlation structures, as prior simulations (contrary to posterior simulations) effectively include cosmic variance, but they exclude realization-specific correlations from non-linear degeneracies. Consequently, they quantify fundamentally different types of uncertainties. We argue that as a result, they also have different and complementary scientific uses, even if this dichotomy is not absolute. In particular, posterior simulations are in general more convenient for parameter estimation studies, while prior simulations are generally more convenient for model testing. Before BEYONDPLANCK, most pipelines used a mix of constrained and random inputs and applied the same hybrid simulations for all applications, even though the statistical justification for this is not always evident. BEYONDPLANCK represents the first end-to-end CMB simulation framework that is able to generate both types of simulations and these new capabilities have brought this topic to the forefront. The BEYONDPLANCK posterior simulations and their uses are described extensively in a suite of companion papers. In this work, we consider one important applications of the corresponding prior simulations, namely, code validation. Specifically, we generated a set of one-year LFI 30 GHz prior simulations with known inputs and we used these to validate the core low-level BEYONDPLANCK algorithms dealing with gain estimation, correlated noise estimation, and mapmaking.

79 ASTRONOMY AND ASTROPHYSICS↗

The BOSS bispectrum analysis at one loop from the Effective Field Theory of Large-Scale Structure

Abstract We analyze the BOSS power spectrum monopole and quadrupole, and the bispectrum monopole and quadrupole data, using the predictions from the Effective Field Theory of Large-Scale Structure (EFTofLSS). Specifically, we use the one loop prediction for the power spectrum and the bispectrum monopole, and the tree level for the bispectrum quadrupole. After validating our pipeline against numerical simulations as well as checking for several internal consistencies, we apply it to the observational data. We find that analyzing the bispectrum monopole to higher wavenumbers thanks to the one-loop prediction, as well as the addition of the tree-level quadrupole, significantly reduces the error bars with respect to our original analysis of the power spectrum at one loop and bispectrum monopole at tree level. After fixing the spectral tilt to Planck preferred value and using a Big Bang Nucleosynthesis prior, we measureσ 8 = 0.794 ± 0.037,h = 0.692 ± 0.011, and Ω m = 0.311 ± 0.010 to about 4.7%, 1.6%, and 3.2%, at 68% CL, respectively.This represents an error bar reduction with respect to the power spectrum-only analysis of about 30%, 18%, and 13% respectively. Remarkably, the results are compatible with the ones obtained with a power-spectrum-only analysis, showing the power of the EFTofLSS in simultaneously predicting several observables. We find no tension with Planck.

Astronomy & Astrophysics↗

The Signatures of Self-interacting Dark Matter and Subhalo Disruption on Cluster Substructure

The abundance, distribution, and inner structure of satellites of galaxy clusters can be sensitive probes of the properties of dark matter. We run 30 cosmological zoom-in simulations with self-interacting dark matter (SIDM), with a velocity-dependent cross section, to study the properties of subhalos within cluster-mass hosts. We find that the abundance of subhalos that survive in the SIDM simulations are suppressed relative to their cold dark matter (CDM) counterparts. Once the population of disrupted subhalos—which may host orphan galaxies—is taken into account, satellite galaxy populations in CDM and SIDM models can be reconciled. However, even in this case, the inner structures of subhalos are significantly different in the two dark matter models. We study the feasibility of using the weak-lensing signal from the subhalo density profiles to distinguish between the cold and self-interacting dark matter while accounting for the potential contribution of orphan galaxies. We find that the effects of self-interactions on the density profile of subhalos can appear degenerate with subhalo disruption in CDM, when orphans are accounted for. With current error bars from the Subaru Hyper Suprime-Cam Strategic Program, we find that subhalos in the outskirts of clusters (where disruption is less prevalent) can be used to constrain dark matter physics. In the future, the Vera C. Rubin Observatory Legacy Survey of Space and Time will give precise measurements of the weak-lensing profile and can be used to constrain σ $_{T}$/m at the ∼1 cm$^{2}$ g$^{−1}$ level at v ∼ 2000 km s$^{−1}$.

79 ASTRONOMY AND ASTROPHYSICS↗

Image Characterization of Reactions Generated by an Aeroengine Micromixing Injector for Lean Direct Injection of Hydrogen and Hydrogen/Natural Gas Blend

An aeroengine micromixing injector, originally designed for lean direct injection of jet fuel, was adapted to work with hydrogen, natural gas, and any blend in between. The ultimate goal of the design was to achieve low NOx emissions when operated on pure hydrogen; to better characterize its performance, flame imaging diagnostics was implemented. Three types of cameras mounted parallel to the injector base and pointing towards the quartz combustor cylinder with the same angle were operated simultaneously: Nikon D90, Dynacolor FB-N9-U, and Phantom v7.1, to capture the visible spectrum, OH* chemiluminescence, and dynamic behavior of the flames, respectively. The current work also includes the overlapping of visible spectrum and OH* chemiluminescence, as an effort to qualitatively define the heat release over the flame area, i.e., UV over visible domain. Additionally, a z-type schlieren configuration was employed to reinforce the presence of some small-scale details occurring near the injector outlet ports. Using the airbox pressure drop, preheat temperature, fuel composition and flame temperature as the factors of study, a Box Behnken model was designed yielding into a 16-points matrix for the operability region, at atmospheric pressure. The main responses were extracted from the OH* images, because of its correlation with the flame heat release, these being the flame area and its average brightness, the heat release area and its center of gravity and leading edge. To assess the importance of the main factors and their interactions, an analysis of variance (ANOVA) was performed for each response, considering as significant each contributor with a p-value below 0.05. For these five responses it was found that all four main factors need to be included into the model (most of them because its own significance, others to ensure the model hierarchy), as well as some common interactions as the preheat and flame temperatures or pressure drop and fuel composition, and quadratic terms from fuel composition or flame temperature. Three of the five models presented R2 and coefficient of variance (C.V.) indicators around 0.97 and 5%, respectively; the leading edge ANOVA yielded into R2=0.78 and C.V.=18.88%; and the heat release area showed R2=0.91 and C.V.=24.13%. Lastly, the confinement ratio effect of the combustor over the flame structure was quantified, for three different conditions: the original 80 mm inner diameter and 200 mm length cylinder was compared against a shorter 150 mm tube (same inner diameter) and a narrower 47 mm inner diameter one (same 200 mm length). The length had a smaller effect than the inner diameter, even though both factors were not strong enough to make a significant deviation on most of the flame parameters: the difference was within the repeatability margin of error for the brightness, flame area, and center of gravity and leading edge of the heat release area.

Imaging diagnostics, flame characterization, OH* c↗

Quantum block encoding for one-pair semiseparable matrices

Quantum block encoding (QBE) is a crucial step in the development of most quantum algorithms, as it provides an embedding of a given matrix into a suitable larger unitary matrix. Historically, the development of efficient techniques for QBE has mostly focused on sparse matrices; less effort has been devoted to data-sparse (e.g., rank-structured) matrices. In this work we examine a particular case of rank structure, namely, one-pair semiseparable matrices. We present a new block encoding approach that relies on a suitable factorization of the given matrix as the product of triangular and diagonal factors. To encode the matrix, the algorithm needs $2\log(N)+7$ ancillary qubits. Assuming that the data input oracles can be implemented with polylogarithmic depth, or that a QRAM input model is available, our proposed method requires $\mathcal{O}({\rm polylog} (N))$ time and has an error of $\mathcal{O}(N^2)$, where $N$ is the matrix size.

Antonioli, Giacomo [Pisa U.; CERN] (ORCID:00090000↗

Accurate Effective Stress Measures: Predicting Creep Life for 3D Stresses Using 2D and 1D Creep Rupture Simulations and Data

Operating structural components experience complex loading conditions resulting in 3D stress states. Current design practice estimates multiaxial creep rupture life by mapping a general state of stress to a uniaxial creep rupture correlation using effective stress measures. The data supporting the development of effective stress measures are nearly always only uniaxial and biaxial, as 3D creep rupture tests are not widely available. This limitation means current effective stress measures must extrapolate from 2D to 3D stress states, potentially introducing extrapolation error. In this work, we use a physics-based, crystal plasticity finite element model to simulate uniaxial, biaxial, and triaxial creep rupture. Here, we use the virtual dataset to assess the accuracy of current and novel effective stress measures in extrapolating from 2D to 3D stresses and also explore how the predictive accuracy of the effective stress measures might change if experimental 3D rupture data was available. We confirm these conclusions, based on simulation data, against multiaxial creep rupture experimental data for several materials, drawn from the literature. The results of the virtual experiments show that calibrating effective stress measures using triaxial test data would significantly improve accuracy and that some effective stress measures are more accurate than others, particularly for highly triaxial stress states. Results obtained using experimental data confirm the numerical findings and suggest that a unified effective stress measure should include an explicit dependence on the first stress invariant, the maximum tensile principal stress, and the von Mises stress.

36 MATERIALS SCIENCE↗

Graphical Gaussian Process Regression Model for Aqueous Solvation Free Energy Prediction of Organic Molecules in Redox Flow Battery

The solvation free energy of organic molecules is a critical parameter in determining emergent properties such as solubility, liquid-phase equilibrium constants, and pKa and redox potentials in an organic redox flow battery. In this work, we present a machine learning (ML) model that can learn and predict the aqueous solvation free energy of an organic molecule using Gaussian process regression method based on a new molecular graph kernel. To investigate the performance of the ML model on electrostatic interaction, the nonpolar interaction contribution of solvent and the conformational entropy of solute in solvation free energy, three data sets with implicit or explicit water solvent models, and contribution of conformational entropy of solute are tested. We demonstrate that our ML model can predict the solvation free energy of molecules at chemical accuracy with a mean absolute error of less than 1 kcal/mol for subsets of the QM9 dataset and the Freesolv database. To solve the general data scarcity problem for a graph-based ML model, we propose a dimension reduction algorithm based on the distance between molecular graphs, which can be used to examine the diversity of the molecular data set. It provides a promising way to build a minimum training set to improve prediction for certain test sets where the space of molecular structures is predetermined.

25 ENERGY STORAGE↗

Probabilistic inversion of circular phase spectra: application to two-station phase-velocity dispersion estimation in western Canada

SUMMARY Periodic directional and temporal measurements are common in seismology, and necessitate specific statistical analyses that are appropriate for circular quantities. In this work, we explore the use of a von Mises distribution as a representation of errors on circular seismological observations. Specifically, we automate the estimation of surface-wave phase-velocity dispersion for the teleseismic two-station method, which generally suffers from a 2π phase ambiguity. The use of Bayesian inverse techniques, which aim to rigorously quantify model parameter uncertainty, have become widespread throughout seismology over the last decade. Here, we apply Bayesian inversion to measurements of surface-wave phase spectra in order to estimate 1-D, path-averaged Earth structure between station pairs. The dispersion curve and associated uncertainties are additional results of the inversion, which can then be used as input for subsequent analyses (e.g. tomography). We demonstrate this technique through application to surface-wave recordings from long-running seismic stations throughout western Canada. Our results for over 10 000 station pairs reveal first-order tectonic features consistent with previous studies, which provides confidence in our approach as well as an initial step towards resolving a full 3-D seismic velocity model for the region. This work also presents a foundation for the inversion of surface-wave phase spectra to estimate 3-D Earth structure directly. Finally, the ideas presented in this work are not limited to the inversion of surface-wave phase spectra, but can also be considered for Bayesian geophysical inversion of any circular quantities.

Gosselin, Jeremy M. (ORCID:0000000203754102)↗

Predicting initial dissolution rates using structural features from molecular dynamics simulations

Predicting chemical durability of glass materials is important for various applications from daily life such as drink glass and kitchen ware to advanced technologies such as nuclear waste disposal and biomedicine. In this work, we explored prediction of initial dissolution rate through structural features from molecular dynamics (MD) simulations for a wide range of glass compositions (total 28) including borosilicate and aluminosilicate glasses, ZrO 2 -containing and V 2 O 5 -containing boroaluminosilicate glasses. The initial dissolution rates (r 0 ) measured experimentally at 90 °C with varying solution conditions were correlated with structural features (e.g., polyhedral linkages and non-bridging oxygen species) obtained from MD simulations, either from this study or from literature. Since hydrolysis of the glass network through breaking of the network former linkages (e.g., Si-O-Si, Si-O-Al, etc.) is a critical step of network glass dissolution, the statistics of these linkages obtained from MD were correlated to r0 through linear regression, where the coefficient of determination (R 2 ) and root mean square error are found to be 0.949 and 0.681, respectively. This model was compared and discussed with existing models developed by various approaches including machine learning, the kinetic rate equation, topological constraint theory, and other descriptors from MD simulations. The discussion provides insights on future model improvements to predict glass dissolution. In addition, the impact of V 2 O 5 on the glass dissolution was examined in detail, implicating that the impact is not the same across all glass compositions and test conditions.

36 MATERIALS SCIENCE↗

Evaluating Radiation Impact on Superconducting Transmon Qubits in Above and Underground Facilities

Superconducting qubits can be sensitive to sudden energy deposits into the substrate from cosmic rays and ambient radioactivity. While previous studies have primarily focused on exploring correlated effects over time and distance due to cosmic rays, this study uniquely offers a direct comparison of a transmon qubit's response at Fermilab SQMS above-ground facilities and the deep-underground Gran Sasso Laboratory (INFN-LNGS, Italy). Despite the vastly different radiation environments, we observe a similar average qubit lifetime (T1) of approximately 80 microseconds at both sites. We then employ a fast decay detection protocol and investigate the time structure, sensitivity, and relative rates of triggered events due to radiation versus intrinsic noise. We compare the above and underground performance of several high-coherence qubits with similar designs. Our findings reveal no significant difference in events with radiation-like signatures between above-ground and underground conditions for these sapphire and niobium-based transmon qubits. This suggests that most such events are likely due to other noise sources, which predominantly contribute to single-qubit errors in contemporary transmon qubits. By exposing the qubits to gamma sources with varying activity levels, we further evaluate their response to radiation in a low-background environment. The results show that qubits can indeed detect particle impacts; however, this capability is only realized by substantially reducing other noise sources, which in turn affects detection efficiency [1].

Roy, Tanay↗