Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “forward modeling”

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,279 records · Page 71

High-bandwidth image-based predictive laser stabilization via optimized Fourier filters

Controlling the delivery of kHz-class pulsed lasers is of interest in a variety of industrial and scientific applications, from next-generation laser-plasma acceleration to laser-based x-ray emission and high-precision manufacturing. The transverse position of the laser pulse train on the application target is often subject to fluctuations by external drivers (e.g., room cooling and heating systems, motorized optics stages and mounts, vacuum systems, chillers, and/or ground vibrations). For typical situations where the disturbance spectrum exhibits discrete peaks on top of a broad-bandwidth lower-frequency background, traditional PID (proportional-integral-derivative) controllers may struggle, since as a general rule PID controllers can be used to suppress vibrations up to only about 5%–10% of the sampling frequency. Here, a predictive feed-forward algorithm is presented that significantly enhances the stabilization bandwidth in such laser systems (up to the Nyquist limit at half the sampling frequency) by online identification and filtering of one or a few discrete frequencies using optimized Fourier filters. Furthermore, the system architecture demonstrated here uses off-the-shelf CMOS cameras and piezo-electric actuated mirrors connected to a standard PC to process the alignment images and implement the algorithm. To avoid high-end, high-cost components, a machine-learning-based model of the piezo mirror’s dynamics was integrated into the system, which enables high-precision positioning by compensating for hysteresis and other hardware-induced effects. A successful demonstration of the method was performed on a 1 kHz laser pulse train, where externally-induced vibrations of up to 400 Hz were attenuated by a factor of five, far exceeding what could be done with a standard PID scheme.

Natal, Joseph↗

Development and Application of a Data-Driven Methodology for Validation of Risk-Informed Safety Margin Characterization Models

The document is the Final Technical Progress Report for the Nuclear Energy University Program’s Integrated Research Project (IRP) on “Development and Application of a Data-Driven Methodology for Validation of Risk-Informed Safety Margin Characterization Models” that supports the LWR Sustainability Program’s RISMC R&D Pathway. The project goal is to develop and demonstrate a data-driven methodology for validation of advanced computer models used in nuclear power plant safety analysis. Specifically, the advanced computer models are those in the toolkit developed to support risk-informed safety margin characterization (RISMC), an integrated deterministic/probabilistic safety analysis methodology developed in the Department of Energy’s Light Water Reactor Sustainability (LWR-S) program. The report reflects the progress made towards the project’s stated goal by contributions by researchers and graduate students from universities, with support from researchers from national laboratories and industry companies. The project organization, effort coordination and technical implementation are summarized, followed by discussion of main findings, issues, and path forward. Selected chapters provide a more detailed description of tasks, approaches and respective findings and recommendations. Noteworthy are contributions that serve as guidelines for methodology development. It is also noted that this report is complemented by other milestone reports (as stand-alone deliverables) that provide detailed discussion of the technical developments. The project results have been documented in a number (12) dissertations and these, 50+ peer-reviewed publications in technical journals and conference proceedings.

42 ENGINEERING↗

2022 Standard Scenarios Report: A U.S. Electricity Sector Outlook

This report documents the eighth edition of the annual Standard Scenarios. It summarizes 70 forward-looking scenarios of the U.S. electricity sector that have been designed to capture a wide range of possible futures. In August 2022, the United States Congress passed the Inflation Reduction Act (IRA), a law aimed at accelerating U.S. decarbonization, clean energy manufacturing, and deployment of new power and end-use technologies. This year’s scenarios include representations of the main electricity-sector provisions from IRA and the potential impact on electricity demand. The Standard Scenarios are simulated using the Regional Energy Deployment System (ReEDS) model, which projects utility-scale electricity sector evolution for the contiguous United States using a system-wide, least-cost approach subject to policy and operational constraints. A subset of the scenarios are simulated in the PLEXOS production cost model to obtain a broader suite of metrics at the hourly resolution, which are made available through the National Renewable Energy Laboratory’s (NREL’s) annual Cambium data sets. The scenarios can be viewed and downloaded from NREL’s Scenario Viewer. Annual results are available for the full suite of scenarios in the Standard Scenarios projects in the viewer, whereas the Cambium projects contain hourly data for a subset of scenarios. The Standard Scenarios includes a scenario called the Mid-case, which has central or median values for core inputs such as technology costs and fuel prices, moderately paced demand growth averaging 1.3% per year, and electricity sector policies as they existed in September 2022 (including IRA). The remaining 69 scenarios are created by varying inputs such as technology and fuel prices, resource availability, demand growth, whether nascent generation technologies are allowed, and by introducing national decarbonization constraints.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

2025 Advances in NekRS: Supporting improved performance for nuclear applications

This report presents several 2025 advancements in NekRS, a high-fidelity spectral element CFD code developed at Argonne National Laboratory to support the NEAMS thermal-hydraulics program. The forthcoming v25 release consolidates several of these advances, adding new features for portability across heterogeneous GPU architectures, real-time in situ visualization, improved turbulence modeling, and conjugate heat transfer coupling. Over the past year, NekRS has demonstrated strong scalability and performance on DOE’s leading exascale platforms, including Aurora and Frontier, confirming its readiness for some of the largest and most complex simulations attempted to date. These achievements provide a powerful new platform for high-fidelity data generation, which in turn supports the development and validation of advanced closure models critical for reactor safety and design. Significant algorithmic innovations have also been introduced. A new global runtime h-refinement capability simplifies workflows by reducing mesh preparation burdens and enabling coarse-to-fine restarts. Building on this, a novel multigrid strategy was implemented to accelerate pressure and transport solves at scale, addressing long-standing bottlenecks in exascale CFD. Together, these developments improve both the efficiency and accessibility of high-fidelity simulations for reactor-relevant problems. Collectively, these enhancements represent a major step forward in simulation technology, positioning NekRS as a cornerstone of NEAMS efforts to enable accurate, efficient, and scalable high-fidelity analysis of advanced nuclear systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The NuMI Flux Prediction at ICARUS

The Deep Underground Neutrino Experiment (DUNE) is a next-generation long-baseline neutrino oscillation experiment seeking to probe fundamental symmetries within the structure of the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix, and perform precision measurements its parameters including the neutrino mass ordering via the sign of $\Delta m^2_{31}$, and the charge-parity violating phase, $\delta_{CP}$. To make these measurements with high precision, DUNE will require external $\nu$-Ar scattering cross section data as a crucial input to the oscillation fit. Imaging Cosmic And Rare Underground Signals (ICARUS) is a 476 t liquid argon neutrino detector located at Fermi National Accelerator Laboratory (FNAL) where it is serving as the far detector for the Short-Baseline Neutrino (SBN) program along the Booster Neutrino Beam (BNB) axis. ICARUS additionally lies 795 m downstream and 100.1 mrad off-axis of the Neutrinos at the Main Injector (NuMI) neutrino beam. From this position, ICARUS is exposed to a large flux of NuMI (anti-)electron and (anti-)muon neutrinos, and poses a unique opportunity to provide high-statistics measurements of quasi-elastic and single pion-production cross sections for four neutrino flavors ($\nu_{\mu}$, $\nu_{e}$, $\bar{\nu}_{\mu}$, $\bar{\nu}_{e}$). This dissertation is centered around accurately characterizing the models and estimating their precision for use in making these measurements. This includes identifying major sources of uncertainty in the models such that they can be properly propagated to the cross section measurements. Specifically, this work focused on the model of the NuMI beamline and its impact on the neutrino fluxes, but also delved into the detector response model and its impact on reconstructed observables in the detector. Significant efforts were made to improve the characterization to enhance precision and thus reduce the level of propagated uncertainty. In particular, the NuMI flux was determined to be composed of 57% $\nu_{\mu}$, 38% $\bar{\nu}_{\mu}$, 3% $\nu_{e}$, and 2% $\bar{\nu}_{e}$ while the horns are operating in the positive-particle focusing configuration. The total uncertainty on the $\nu_{\mu} + \bar{\nu}_{\mu}$ ($\nu_{e} + \bar{\nu}_{e}$) flux while operating in the forward horn operating mode was determined to be 10.84% (9.04%). Compared to the on-axis flux, mesons that eventually decay to neutrinos more frequently reinteract within the NuMI structure, resulting in elevated uncertainty as these processes are not well-constrained by existing hadron interaction cross section measurements. Covariance matrices were calculated to propagate the flux uncertainty characterization to NuMI analyses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Quantum complexity of time evolution with chaotic Hamiltonians

We study the quantum complexity of time evolution in large-N chaotic systems, with the SYK model as our main example. This complexity is expected to increase linearly for exponential time prior to saturating at its maximum value, and is related to the length of minimal geodesics on the manifold of unitary operators that act on Hilbert space. Using the Euler-Arnold formalism, we demonstrate that there is always a geodesic between the identity and the time evolution operator e -iHt whose length grows linearly with time. This geodesic is minimal until there is an obstruction to its minimality, after which it can fail to be a minimum either locally or globally. We identify a criterion — the Eigenstate Complexity Hypothesis (ECH) — which bounds the overlap between offdiagonal energy eigenstate projectors and the k-local operators of the theory, and use it to argue that the linear geodesic will at least be a local minimum for exponential time. We show numerically that the large-N SYK model (which is chaotic) satisfies ECH and thus has no local obstructions to linear growth of complexity for exponential time, as expected from holographic duality. In contrast, we also study the case with N = 2 fermions (which is integrable) and find short-time linear complexity growth followed by oscillations. Our analysis relates complexity to familiar properties of physical theories like their spectra and the structure of energy eigenstates and has implications for the hypothesized computational complexity class separations PSPACE * BQP/poly and PSPACE * BQSUBEXP/subexp, and the “fast-forwarding” of quantum Hamiltonians.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Measurement of nuclear effects on ψ(2S) production in p-Pb collisions at $ \sqrt{s_{\mathrm{NN}}}$ = 8.16 TeV

Inclusive ψ (2S) production is measured in p-Pb collisions at the centre-of-mass energy per nucleon-nucleon pair \( \sqrt{s_{\mathrm{NN}}} \) = 8.16 TeV, using the ALICE detector at the CERN LHC. The production of ψ (2S) is studied at forward (203 < y cms < 3 . 53) and backward (-4.46 < y cms < -2.96) centre-of-mass rapidity and for transverse momentum p T < 12 GeV/ c via the decay to muon pairs. In this paper, we report the integrated as well as the y cms - and p T -differential inclusive production cross sections. Nuclear effects on ψ (2S) production are studied via the determination of the nuclear modification factor that shows a strong suppression at both forward and backward centre-of-mass rapidities. Comparisons with corresponding results for inclusive J/ ψ show a similar suppression for the two states at forward rapidity (p-going direction), but a stronger suppression for ψ (2S) at backward rapidity (Pb-going direction). As a function of p T , no clear dependence of the nuclear modification factor is found. The relative size of nuclear effects on ψ (2S) production compared to J/ ψ is also studied via the double ratio of production cross sections [ σ ψ (2S) /σ J /ψ ]pPb/[ σ ψ (2S) /σ J /ψ ] pp between p-Pb and pp collisions. The results are compared with theoretical models that include various effects related to the initial and final state of the collision system and also with previous measurements at \( \sqrt{s_{\mathrm{NN}}} \) = 5.02 TeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An assessment of the harvesting and fuel performance of a single-pass cut-and-chip harvester in commercial-scale short-rotation poplar crops as influenced by crop and weather conditions

There is a need for data on commercial-scale harvesting operations in hybrid poplar short rotation crops to better understand costs and logistics, especially for modeling large scale biorefineries. An investigation was conducted on the in-field performance of a single-pass cut-and-chip harvester operating in commercial fields for over 370 individual wagon loads representing a range of crop and field conditions. Mean material capacity peaked at 70 Mg h -1 in dry conditions in lower biomass (<60 Mg ha -1 ), but was only 28 Mg h -1 during wet conditions and often higher standing biomass. Wet ground conditions require the harvester to divert additional power to maintaining forward movement, which results in decreased material capacity and increased fuel consumption. Crops with higher standing biomass had taller trees that do not always feed as smoothly into the harvester, slowing forward progress and lowering material capacity. Mean crop specific fuel consumption, L of fuel per Mg of biomass processed, generally decreased as standing biomass increased. When standing biomass was above 40 Mg ha -1 mean crop specific fuel consumption (FCC) was 1.69 L Mg -1 in dry conditions and 3.98 L Mg -1 in wet conditions, but when standing biomass was below 40 Mg ha-1 FCC increased drastically (as high as 5 L Mg -1 ) because the harvester is putting more effort into forward speed instead of material processing. Developing relationships between material capacity and fuel consumption based on standing biomass and ground conditions at representative scales are essential for conducting environmental and economic analyses of these systems.

09 BIOMASS FUELS↗

Understanding structure-processing relationships in metal additive manufacturing via featurization of microstructural images

Understanding and predicting accurate property-structure-processing relationships for additively manufactured components is important for both forward and inverse design of robust, reliable parts and assemblies. While direct mapping of process parameters to properties is sometimes plausible, it is often rendered difficult due to poor microstructural control. Exploring the direct relationship between processing conditions and microstructural features can thus provide significant physical insights and aid the overall design process. Here, in this study, we develop an automated high-throughput framework to simulate an uncertainty-aware additive manufacturing (AM) process, characterize microstructural images, and extract meaningful features/descriptors. A kinetic Monte Carlo (KMC) based model of the AM process is used to simulate microstructural evolution for a diverse set of experimentally relevant processing conditions. We perform a parametric study to explore the relationship between microstructural features and processing conditions. Our results indicate that a many-to-one mapping can exist between processing conditions and typical descriptors; therefore, multiple descriptors are thus necessary to unambiguously represent microstructural images. Our work provides crucial quantitative and qualitative in-formation that would aid in the selection of features for microstructural images. Featurized microstructures could then be utilized to build data-driven models for predictive control of microstructures and thereby properties of additively manufactured components.

36 MATERIALS SCIENCE↗

First measurement of the forward rapidity gap distribution in pPb collisions at $\sqrt{s_{NN}}$ = 8.16 TeV

For the first time at LHC energies, the forward rapidity gap spectra from proton-lead collisions for both proton and lead dissociation processes are presented. The analysis is performed over 10.4 units of pseudorapidity at a center-of-mass energy per nucleon pair of $\sqrt{s_{NN}}$ = 8.16 TeV , almost 300 times higher than in previous measurements of diffractive production in proton-nucleus collisions. For lead dissociation processes, which correspond to the pomeron-lead event topology, the EPOS-LHC generator predictions are a factor of 2 below the data, but the model gives a reasonable description of the rapidity gap spectrum shape. For the pomeron-proton topology, the EPOS-LHC , QGSJET II , and HIJING predictions are all at least a factor of 5 lower than the data. The latter effect might be explained by a significant contribution of ultraperipheral photoproduction events mimicking the signature of diffractive processes. These data may be of significant help in understanding the high energy limit of quantum chromodynamics and for modeling cosmic ray air showers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The iron oxidation state of Ryugu samples

The Hayabusa2 mission sampled Ryugu, an asteroid that did not suffer extensive thermal metamorphism, and returned rocks to the Earth with no significant air exposure. It therefore offers a unique opportunity to study the redox state of carbonaceous Cb-type asteroids and evaluate the overall redox state of the most primitive rocks of the solar system. An analytical framework was developed to investigate the iron mineralogy and valence state in extraterrestrial material at the micron scale by combining x-ray diffraction, conventional Mössbauer (MS), and nuclear forward scattering (NFS) spectroscopies. An array of standard minerals was analyzed and cross-calibrated between MS and NFS. Then, MS and NFS spectra on three Ryugu grains were collected at the bulk and the micron scales. In Ryugu samples, iron is essentially accommodated in magnetite, clay minerals (serpentine–smectite), and sulfides. Only a single set of Mössbauer parameters was necessary to account for the entire variability observed in MS and NFS spectra, at all spatial scales investigated. These parameters therefore make up a fully consistent iron mineralogical model for the Ryugu samples. As far as MS and NFS spectroscopies are concerned, Ryugu grains are overall similar to each other and share most of their mineralogical features with CI-type chondrites. In detail however, no ferrihydrite is found in Ryugu particles even at the very sensitive scale of Mössbauer spectroscopy. The typical Fe 3+ /Fe tot of clay minerals is much lower than typical redox ratios measured in CI chondrites (Fe 3+ /Fe tot = 85%–90%). Furthermore, magnetite from Ryugu is stoichiometric with no significant maghemite component, whereas up to 12% of maghemite was previously identified in the Orgueil's so-called magnetite. These differences suggest that most CI meteorites suffered terrestrial alteration and that the preterrestrial composition of these carbon-rich samples was less oxidized than previously measured. However, it is not clear yet whether or not the parent bodies of CI chondrites were as reduced as Ryugu. Finally, the high spatial resolution of NFS allows to disentangle the redox state and the crystal chemistry of iron accommodated in serpentine and smectite. The most likely polytype of serpentine is lizardite, containing <35% of Fe 3+ , a fraction of which being tetrahedrally coordinated. Smectite is more oxidized (Fe 3+ /Fe tot > 65%) and mainly contains octahedral ferric iron. In conclusion, this finding implies that these clays formed from highly alkaline fluids and the spatial variability highlighted here may suggest a temporal evolution or a spatial variability of the nature of this fluid.

58 GEOSCIENCES↗

Machine learning for domain transfer between simulated and experimental 2D X-ray diffraction patterns using generative adversarial networks

X-ray diffraction (XRD) is a well-established technique for analyzing materials at an atomic level. Dynamic compression experiments (DCE), in which materials are subject to extreme pressures, can provide fundamental understanding to pressure-induced phase transitions and compression of the crystal lattice. The analysis of XRD patterns from highly compressed samples is non-trivial given the sparsity of data, high experimental costs, and the fact that the data is often marred with X-ray background and other artifacts. While accurate computational frameworks exist, they solve the forward problem—from structures and orientations to XRD patterns. Solving the inverse problem for 2D experimental diffraction patterns is currently a complex manual process of matching and comparing experimentally observed patterns to computationally generated ones. Machine learning is a promising tool for automating the matching process but often requires data-intensive architectures. Here, in this study, we use a CycleGAN to translate the domain of limited experimental data to a domain in which there is readily available simulated data. This domain shift allows data-intensive machine learning models that have only been trained on simulated XRD patterns to be used in the analysis of experiments.

Brozak, Samantha Jean [Sandia National Laboratorie↗

An ultra-fast method for generating synthetic down-scattered neutron data for inertial confinement fusion implosions

In inertial confinement fusion experiments at the National Ignition Facility, asymmetries are probed by a variety of neutron diagnostics, including neutron imaging systems, real-time neutron activation diagnostics (RTNADs), and neutron spectrometers. It is often useful to generate synthetic data based on these diagnostics to validate and tune models. However, current methods of doing so using Monte Carlo particle tracing are time-consuming. In this paper, an ultra-fast method is presented for generating synthetic neutron images, RTNAD data, and spectrometry data using line integrals and 3D convolutions. While it does not contain as much physics as particle tracing codes, it is thousands of times faster and produces nearly identical data. This enables analysis techniques that depend on generating large amounts of synthetic data, which will prove very useful for the study of asymmetries going forward.

Deuterium↗

Third moments of nucleon unpolarized, polarized, and transversity parton distribution functions from physical-point lattice QCD

Using forward matrix elements of local leading-twist operators, we present a determination of the isovector third Mellin moments ⟨𝑥 2 ⟩ of nucleon unpolarized, polarized, and transversity parton distribution functions. Two lattice QCD ensembles at the physical pion mass are used, which were generated using a tree-level Symanzik-improved gauge action and 2+1 flavor tree-level improved Wilson Clover fermions coupling via 2-level HEX-smearing. Leveraging a wide set of operators, two extraction methods for the matrix elements, and the automatic inclusion of model uncertainties via bootstrapped model averages, we extract values of the third Mellin moments. Furthermore, this is the first direct calculation of these observables performed at the physical pion mass.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Commercial, industrial, and institutional discount rate estimation for efficiency standards analysis: Sector-level data 1998–2021

Underlying each of the Department of Energy’s (DOE’s) federal appliance and equipment energy conservation standards are a set of complex analyses of the projected costs and benefits of regulation. Any new or amended standard must be designed to achieve significant additional energy conservation, provided that it is technologically feasible and economically justified (42 U.S.C. 6295(o)(2)(A)). DOE determines economic justification based on whether the benefits exceed the burdens, considering a variety of factors, including the economic impact of the standard on consumers of the product and the savings in lifetime operating cost compared to any increase in price or maintenance expenses (42 U.S.C. 6295(o)(2)(B)). As part of this determination, DOE conducts a Life-Cycle Cost (LCC) analysis, which models the combined impact of appliance first cost and operating cost changes on a representative commercial building sample in order to identify the fraction of customers achieving LCC savings or incurring net cost at the considered efficiency levels. Thus, the commercial discount rate value(s) used to calculate the present value of energy cost savings within the LCC model implicitly plays a role in estimating the economic impact of potential standard levels. This report provides an in-depth discussion of the commercial discount rate estimation process. It is an update to previous reports on estimating commercial discount rates from firm-level financial data (Fujita, 2016). Major topics covered in this report include: Discount rate estimation methods and rationale; -Data sources used and data limitations; -Discount rate distributions for use in standards analysis; -Discount rate estimation methods and distributions specific to the small business subgroup analysis. Going forward, this report will be updated as data allow and analyses necessitate.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Gluon Sivers function from forward exclusive χ c 1 photoproduction on unpolarized protons

Exclusive production of a χ c 1 axial vector quarkonia in photon-proton scattering at high energies requires a C -odd t -channel exchange. In the limit of vanishing momentum transfer this occurs either via the exchange of a photon, the Primakoff process, where the spin of the proton does not change. For axial-vector meson production, as a consequence of the Landau-Yang theorem, the Primakoff cross section is finite as t → 0 . Alternatively, a C -odd spin dependent Odderon can be exchanged, which involves a spin flip of the proton. The resulting cross section is related to the square of the collinear trigluon correlator or the k ⊥ -moment of the gluon Sivers function. Using two models for the gluon Sivers function from the literature we compute the ratio of Sivers to Primakoff cross sections and the angular coefficient λ θ governing the angular distribution of the χ c 1 → J / ψ + γ decay as functions of x . We point out that these observables constrain the magnitude of the gluon Sivers function at small x which could be accessed in electron-proton scattering and ultraperipheral proton-proton and nucleus-proton collisions. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Tackling Disorder in γ‐Ga 2 O 3

Abstract Ga 2 O 3 and its polymorphs are attracting increasing attention. The rich structural space of polymorphic oxide systems such as Ga 2 O 3 offers potential for electronic structure engineering, which is of particular interest for a range of applications, such as power electronics. γ‐Ga 2 O 3 presents a particular challenge across synthesis, characterization, and theory due to its inherent disorder and resulting complex structure–electronic‐structure relationship. Here, density functional theory is used in combination with a machine‐learning approach to screen nearly one million potential structures, thereby developing a robust atomistic model of the γ‐phase. Theoretical results are compared with surface and bulk sensitive soft and hard X‐ray photoelectron spectroscopy, X‐ray absorption spectroscopy, spectroscopic ellipsometry, and photoluminescence excitation spectroscopy experiments representative of the occupied and unoccupied states of γ‐Ga 2 O 3 . The first onset of strong absorption at room temperature is found at 5.1 eV from spectroscopic ellipsometry, which agrees well with the excitation maximum at 5.17 eV obtained by photoluminescence excitation spectroscopy, where the latter shifts to 5.33 eV at 5 K. This work presents a leap forward in the treatment of complex, disordered oxides and is a crucial step toward exploring how their electronic structure can be understood in terms of local coordination and overall structure.

36 MATERIALS SCIENCE↗

Revealing the Origin and Nature of the Buried Metal‐Substrate Interface Layer in Ta/Sapphire Superconducting Films

Abstract Despite constituting a smaller fraction of the qubit's electromagnetic mode, surfaces and interfaces can exert significant influence as sources of high‐loss tangents, which brings forward the need to reveal properties of these extended defects and identify routes to their control. Here, we examine the structure and composition of the metal‐substrate interfacial layer that exists in Ta/sapphire‐based superconducting films. Synchrotron‐based X‐ray reflectivity measurements of Ta films, commonly used in these qubits, reveal an unexplored interface layer at the metal‐substrate interface. Scanning transmission electron microscopy and core‐level electron energy loss spectroscopy identified an intermixing layer (≈0.65 ± 0.05 nm) at the metal‐substrate interface containing Al, O, and Ta atoms. Density functional theory modeling reveals that the structure and properties of the Ta/sapphire heterojunctions are determined by the oxygen content on the sapphire surface prior to Ta deposition for two atomic terminations of sapphire. Using a multimodal approach, we gained deeper insights into the interface layer between the metal and substrate, which suggests that the orientation of deposited Ta films depend on the surface termination of sapphire. The observed elemental intermixing at the metal‐substrate interface influences the thermodynamic stability and electronic behavior of the film, which may also affect qubit performance.

36 MATERIALS SCIENCE↗