Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion length”

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 901 records · Page 50

Particle-in-Cell Simulation of Collisionless Driven Reconnection with Open Boundaries

First results are discussed from an ongoing study of driven collisionless reconnection using a 2 1/2-dimensional electromagnetic particle-in-cell simulation model with open inflow and outflow boundaries. An extended electron diffusion region (EEDR) is defined as that region surrounding a reconnecting neutral line in which the out-of-plane nonideal electric field is positive. It is shown that the boundaries of this region in the directions of the outflow jets are at the positions where the electrons make the transition from unfrozen meandering motion in the current sheet to outward drifting with the magnetic field in the outflow jets; a turning length scale is defined to mark these positions, The initial width of the EEDR in the inflow directions is comparable to the electron bounce width. Later. as shoulders develop to form a two-scale structure. thc EEDR width expands to the ion bounce width scale. The inner portion of the EEDR or the electron diffusion region proper remains at the electron bounce width. Two methods are introduced for predicting the reconnection electric field using the dimensions of the EEDR. These results are interpreted as further evidence that the EEDR is the region that is relevant to understanding the electron role in the neutral line vicinity.

Kimas, Alex↗

Generative modeling enables molecular structure retrieval from Coulomb explosion imaging

Capturing the structural changes that molecules undergo during chemical reactions in real space and time is a long-standing dream and an essential prerequisite for understanding and ultimately controlling femtochemistry. A key approach to tackle this challenging task is Coulomb explosion imaging, which benefited decisively from recently emerging high-repetition-rate X-ray free-electron laser sources. With this technique, information on the molecular structure is inferred from the momentum distributions of the ions produced by the rapid Coulomb explosion of molecules. Retrieving molecular structures from these distributions poses a highly non-linear inverse problem that remains unsolved for molecules consisting of more than a few atoms. Here, we address this challenge using a diffusion-based Transformer neural network. We show that the network reconstructs unknown molecular geometries from ion-momentum distributions with a mean absolute error below one Bohr radius, which is half the length of a typical chemical bond.

Artificial Intelligence (cs.AI)↗

Convective shells in the interior of Cepheid variable stars: Overshooting models based on hydrodynamic simulations

Context. Because Cepheid variable stars have long been used as a cosmic benchmark for scaling distances in our Galaxy and beyond, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. However, our understanding of the dynamics in the interiors of these physically complex stars is limited. Aims. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we studied the structure of intermediate-mass stars that cross the instability strip. Then, we performed two-dimensional hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations did not model the radial pulsations but focused on the interior structure of this family of stars. We developed and applied a new statistical analysis to examine convection and convective boundary mixing in the interior of these stellar simulations. Results. Based on a grid of MESA models, we demonstrated that a common structure for intermediate mass Cepheids includes an interior convective shell as well as a thin outer convective envelope. Using the extreme value theory approach to analyze our MUSIC simulation data, we found that overshooting above the convective shell fills the space between these convectively unstable layers. We developed a new statistical analysis that provides a clearer picture of how overshooting fills this layer; it also allowed us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak overshooting layer and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of gamma distributions. Using our mixture model, we showed that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We proposed a new form for the diffusion coefficient that addresses the need for overlapping overshooting layers between convective shells. We introduced the idea of a “super-mixing layer” where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.

79 ASTRONOMY AND ASTROPHYSICS↗

Weak localization and universal conductance fluctuations in large-area twisted bilayer graphene

We study diffusive magnetotransport in highly 𝑝-doped large-area twisted bilayer graphene in 1∘, 7∘, 9∘, and 20∘ samples. All samples exhibit weak localization, from which we extract the phase coherence length and intervalley scattering lengths, and from that determine that dephasing is caused by electron-electron scattering and intervalley scattering is caused by point defects. We observe signatures of universal conductance fluctuations in the 9∘ sample, which has high mobility and is near the van Hove singularity. Further improvements in sample quality and applications to large-area moiré materials will open new avenues to observe quantum interference effects.

Talkington, Spenser [University of Pennsylvania]↗

Heat transfer and forces on concave surfaces in free molecule flow.

A Monte Carlo modeling technique is described for mathematically simulating free molecular flows over a concave spherical surface and a concave cylindrical surface of finite length. The half-angle of the surfaces may vary from 0 to 90 degrees, and the incident flow may have an arbitrary speed ratio and an arbitrary angle of attack. Partial diffuse reflection and imperfect energy accommodation for molecules colliding with the surfaces are also considered. Results of heat transfer, drag and lift coefficients are presented for a variety of flow conditions. The present Monte Carlo results are shown to be in very good agreement with certain available theoretical solutions.

Fan, C.↗

Warm neutral halos around molecular clouds. II - H I and CO (J = 1-0) observations

Observational data are presented relating to the boundary regions between dense, molecular clouds and the more diffuse atomic gas. The data set consists of H I 21 cm and CO (J = 1-0) strip maps in 14 molecular clouds which traverse 62 known CO cloud boundaries. The angular lengths of the strips range from 0.5 deg to several degrees of arc, which is generally adequate to eliminate confusion with background (and foreground) H I. Most of the crossed boundaries (49 of 62) show clear evidence for the presence of H I emission maxima, or halos, lying outside the molecular clouds' gas. There is also general evidence for limb brightening in the CO emission.

Wannier, Peter G.↗

The discovery of a radio galaxy of at least 5 Mpc

Context. Giant radio galaxies (GRGs, or colloquially ‘giants’) are the Universe’s largest structures generated by individual galaxies. They comprise synchrotron-radiating active galactic nucleus ejecta and attain cosmological (megaparsec-scale) lengths. However, the main mechanisms that drive their exceptional growth remain poorly understood. Aims. To deduce the main mechanisms that drive a phenomenon, it is usually instructive to study extreme examples. If there exist host galaxy characteristics that are an important cause for GRG growth, then the hosts of the largest GRGs are likely to possess them. Similarly, if there exist particular large-scale environments that are highly conducive to GRG growth, then the largest GRGs are likely to reside in them. For these reasons, we aim to perform a case study of the largest GRG available. Methods. We reprocessed the LOFAR Two-Metre Sky Survey DR2 by subtracting compact sources and performing multi-scale CLEAN de-convolutions at 60" and 90" resolution. The resulting images constitute the most sensitive survey yet for radio galaxy lobes, whose diffuse nature and steep synchrotron spectra have allowed them to evade previous detection attempts at higher resolution and shorter wavelengths. We visually searched these images for GRGs. Results. We have discovered Alcyoneus, a low-excitation radio galaxy with a projected proper length l p = 4.99 ± 0.04 Mpc. Both its jets and lobes are detected at very high significance, and the SDSS-based identification of the host, at spectroscopic redshift z spec = 0.24674 ± 6 × 10 -5 , is unambiguous. The total luminosity density at ν = 144 MHz is L ν = 8 ± 1 × 10 25 W Hz -1 , which is below average, though near median (percentile 45 ± 3%) for GRGs. The host is an elliptical galaxy with a stellar mass M* = 2.4 ± 0.4 × 10 11 M ⊙ and a super-massive black hole mass M• = 4 ± 2 × 10 8 M ⊙ , both of which tend towards the lower end of their respective GRG distributions (percentiles 25 ± 9% and 23 ± 11%). The host resides in a filament of the Cosmic Web. Through a new Bayesian model for radio galaxy lobes in three dimensions, we estimate the pressures in the megaparsec-cubed-scale northern and southern lobes to be P min,1 = 4.8 ± 0.3 × 10 -16 Pa and P min,2 = 4.9 ± 0.6 × 10 -16 Pa, respectively. The corresponding magnetic field strengths are B min,1 = 46 ± 1 pT and B min,2 = 46 ± 3 pT. Conclusions. We have discovered what is in projection the largest known structure made by a single galaxy – a GRG with a projected proper length l p = 4.99 ± 0.04 Mpc. The true proper length is at least l min = 5.04 ± 0.05 Mpc. Beyond geometry, Alcyoneus and its host are suspiciously ordinary: the total low-frequency luminosity density, stellar mass, and super-massive black hole mass are all lower than, though similar to, those of the medial GRG. Thus, very massive galaxies or central black holes are not necessary to grow large giants, and, if the observed state is representative of the source over its lifetime, neither is high radio power. A low-density environment remains a possible explanation. The source resides in a filament of the Cosmic Web, with which it might have significant thermodynamic interaction. The pressures in the lobes are the lowest hitherto found, and Alcyoneus therefore represents the most promising radio galaxy yet to probe the warm–hot inter-galactic medium.

79 ASTRONOMY AND ASTROPHYSICS↗

Nanoparticle Diffusion in Miscible Polymer Nanocomposite Melts

X-ray photon correlation spectroscopy measurements were used to quantify the dynamics of bare and bimodal grafted silica nanoparticles mixed with PEO melts of different molecular weights. In dilute polymer nanocomposite (PNC) samples, we find diffusive NP behavior as described by the Stokes–Einstein relationship so long as the adsorbed PEO polymer layer is taken into account in determining both the effective NP size and its role on composite viscosity. The size of this bound layer was found to be approximately 2R g , where R g is the chain radius of gyration. We also expanded our system to investigate how the dynamics of grafted NPs differ from bare NPs with an adsorbed layer. We showed that the dynamics again can be determined by an effective NP radius at a scale smaller than the effective interparticle spacing; however, at larger length scales, the morphology and grafting parameters play a major role in the system dynamics. These results allow us to quantify NP ordering driven by polymer crystallization. It has previously been speculated that behavior is controlled by the relative ratio of time scale of crystal growth and the diffusive time scale of the NPs, a Peclet number. When the former time scale is longer, then the NPs are expected to be segregated into the interlamellar amorphous zones, while the NPs are too slow to be reorganized in the opposite case. We show here that this conjecture is quantitatively correct and the demarcation in behavior occurs for Pe = 1. Thus, we provide a way to estimate a critical spherulite growth rate for any semicrystalline PNC, at which a given NP can be ordered. Together, the results of this study permit us to tunably design PNCs through directed dispersion of NPs in a polymer matrix.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Stereochemically‐Controlled Fluorinated Copolymers for Selectively Permeable Barrier Applications

Selective oxygen permeability coupled with low water vapor transmission is essential for biomedical and packaging applications requiring controlled oxygen flux under humid conditions. However, most high‐performance barrier polymers depend on perfluoroalkyl substances (PFAS), whose persistence and regulatory restrictions limit their long‐term applicability. We designed a series of stereocontrolled thiol‐yne‐based polyesters, including both fluorinated and non‐fluorinated variants, for selective oxygen permeability with considerable water barrier performance. Tailoring polymer crystallinity and morphology tuned both oxygen transport and mechanical properties. Fluorinated polymers demonstrated enhanced hydrophobicity and water resistance while maintaining oxygen diffusivity within a range relevant to oxygen‐sensing applications. Structure–property relationships were elucidated through small‐ and wide‐angle X‐ray scattering, revealing semi‐crystalline domains influenced by fluorine content and dithiol chain length. Barrier performance was rigorously evaluated via water vapor transmission rate and dynamic vapor sorption, showing reduced water uptake with increasing dithiol monomer length and crystallinity. In conclusion, this work introduces a PFAS‐free alternative to conventional barrier materials and establishes a tunable materials platform with potential relevance for biomedical devices and packaging systems requiring controlled oxygen permeability.

36 MATERIALS SCIENCE↗

Laser powder bed fusion of Inconel 718 on 316 stainless steel

The joining of Inconel 718 and 316 stainless steel has notable industrial importance. However, traditional arc-welding has many issues. In particular, the formation of the detrimental Nb-rich Laves phase in the heat affected zone can promote crack initiation and propagation, which leads to poor mechanical properties. The small heat input and the rapid solidification of laser powder bed fusion (PBF) has the potential to suppress this problem. Despite the infancy of the multi-material PBF, it is of interest to investigate the relationship among the process parameter, the microstructure at the materials joining point, and the mechanical performance. In this study, we used high-speed x-ray imaging to in-situ monitor the printing of Inconel 718 tracks on a 316H SS substrate via laser powder bed fusion, followed by post-build characterization and testing with SEM, XRD and nanoindentation. We observed that a higher linear energy density resulted in a more diffused interface and a smaller compositional gradient, but as a trade-off, generated more keyhole pores and thermal cracks. Furthermore, we also showed that the built layer exhibited chemical inhomogeneity at two length scales: a tens-of-mu m-scale compositional vortex due to the melt flow, and a sub-mu m-scale Nb-rich dendritic structure due to the non-equilibrium solidification that led to the formation of the Laves phase. It was inferred that the chemical inhomogeneity may benefit the mechanical properties at the interface by providing more interlocking between the two materials and also by limiting the effect of the brittle Laves phase to potentially short range.

36 MATERIALS SCIENCE↗

Measurement of Three Transport Coefficients and the Thermodynamic Factor in Block Copolymer Electrolytes with Different Morphologies

The design and engineering of composite materials is one strategy to satisfy the materials needs of systems with multiple orthogonal property requirements. In the case of rechargeable batteries with lithium metal anodes, the system requires a separator with fast lithium ion transport and good mechanical strength. In this work, we focus on the system polystyrene-block-poly(ethylene oxide) (SEO) with bis(trifluoromethane)sulfonimide lithium salt (LiTFSI). Ion transport occurs in the salt-containing poly(ethylene oxide)-rich domains. Mechanical rigidity arises due to the glassy nature of polystyrene (PS). If we assume that the salt does not interact with the PS-rich domains, we can describe ion transport in the electrolyte by three transport parameters (ionic conductivity, κ, salt diffusion coefficient, D, and cation transference number, t + 0 ) and a thermodynamic factor, T f . By systematically varying the volume fraction of the conducting phase, Φ c between 0.29 and 1.0, and chain length, N between 80 and 8000, we elucidate the role of morphology on ion transport. We find that κ is the strongest function of morphology, varying by three full orders of magnitude, while D is a weaker function of morphology. To calculate t + 0 and T f , we measure the current fraction, ρ + , and the open circuit potential, U, of concentration cells. We find that ρ + and U follow universal trends as a function of salt concentration, regardless of chain length, morphology, or Φ c , allowing us to calculate t + 0 for any SEO/LiTFSI or PEO/LiTFSI mixture when κ and D are known. Finally, the framework developed in this paper enables predicting the performance of any block copolymer electrolyte in a rechargeable battery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Compact Diffusion Bonded Heat Exchanger Fatigue Life Simulations (Final Report)

Compact diffusion bonded heat exchangers are essential for high pressure heat exchange, but they are subject to thermal fatigue and ramp rate limitations. Simulation of these geometries is challenging with a large range of length and time scales from thousands of mm-sized microchannels inside a m-sized heat exchanger. Multi-physics simulations including thermal, fluid, and solid mechanics components are being used to predict stress within the heat exchangers under these conditions. These predictions can then be used to understand thermal ramp rate limitations while keeping maximum stresses low as well as fatigue life predictions from well-known empirical models.

42 ENGINEERING↗

Experimental investigation of an axisymmetric free jet with an initially uniform velocity profile

An experimental investigation was conducted to determine the flow characteristics of a circular free helium jet having an initially uniform velocity profile. Complete velocity profiles are presented at Reynolds numbers of 1027 and 4571 at 0, 3, 6, 10, 15, and 20 nozzle diameters (where possible) from the nozzle exit. Centerline velocity decay and potential core length were obtained over a range of Reynolds numbers from 155 to 5349 at distances up to and including 25 nozzle diameters from the nozzle exit. The angles of spread associated with the diffusion of the jet downstream of the nozzle are also given. Axial jet momentum flux and entrained mass flux, at various distances downstream of the nozzle, are presented as a function of the jet Reynolds number.

Labus, T. L.↗

Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows

The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.

Source record↗

The North Galactic Pole Rift and the Local Hot Bubble

The North Galactic Pole Rift (NGPR) is one of the few distinct neutral hydrogen clouds at high Galactic latitudes that have well-defined distances. It is located at the edge of the Local Cavity (LC) and provides an important test case for understanding the Local Hot Bubble (LHB), the presumed location for the hot diffuse plasma responsible for much of the observed 1/4 keV emission originating in the solar neighborhood. Using data from the ROSAT All- Sky Survey and the Planck reddening map, we find the path length within the LC (LHB plus Complex of Local Interstellar Clouds) to be 98 plus or minus 27 pc, in excellent agreement with the distance to the NGPR of 98 +/- 6 pc. In addition, we examine another 14 directions that are distributed over the sky where the LC wall is apparently optically thick at 1/4 keV. We find that the data in these directions are also consistent with the LHB model and a uniform emissivity plasma filling most of the LC.

ISM: bubbles aEuro" ISM: clouds aEuro" ISM: magnet↗

Bayesian Analysis of TRISO Fuel: Quantifying Model Inadequacy, Incorporating Lower-Length-Scale Effects, and Developing Parallel Active Learning Capabilities

The U.S. Department of Energy (DOE)’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel-cycle systems. This program has been providing engineering-scale support for the continued development of BISON, a high-fidelity, high-resolution fuel performance tool. Fuel behavior in nuclear reactors is governed by a complex network of mechanisms that interact with various other physics aspects in the reactor system. Any model developed to represent fuel behavior will likely be idealized, resulting in uncertainties when comparing their predictions against the observed data. In Fiscal Year (FY)-23, we initiated the Uncertainty Quantification (UQ) work by using Bayesian methods to establish a level of model trustworthiness and further improve it, with a particular emphasis on TRI-Structural isOtropic (TRISO) nuclear fuel. This year, we further expanded on that UQ work by investigating an approach to quantifying model inadequacy and accounting for lower-length scale (LLS) effects in TRISO silver (Ag) release modeling. Furthermore, we are implementing parallel active learning capabilities to reduce the computational cost (i.e., required computational resources and elapsed time) of performing UQ. Specifically, we utilized The Kennedy O’Hagan framework for Bayesian uncertainty quantification (KOH) to account for model inadequacy in TRISO Ag release predictions made by BISON. The KOH framework represents an improvement over the standard Bayesian framework used in FY-23. Explicitly accounting for model inadequacy in the Bayesian framework helps establish the level of experimental noise uncertainty in the Advanced Gas Reactor (AGR) data. We compared the inverse UQ results obtained from both the standard Bayesian and KOH frameworks in light of the AGR-2/3/4 data, and also compared the predictive UQ results obtained from these two frameworks in light of the AGR-1 data. Next, we investigated the impact of considering LLS effects in the Ag release simulations. We developed an expanded database of LLS simulated effective diffusivities for Ag, covering a wide range of microstructures and temperatures. Using this database, we developed a framework for incorporating LLS effects into the engineering-scale Ag release UQ. We developed both parametric and non-parametric approaches for bridging the length scales. We then investigated the inverse UQ results in light of the AGR-2/3/4 data and the predictive UQ results in light of the AGR-1 data, and compared the LLS-informed approach and the Arrhenius equation, which does not include microstructure information. Finally, we discussed implementing parallel active learning capabilities in the Multiphysics Object Oriented Simulation Environment (MOOSE)/BISON to reduce the computational cost (i.e., computational resources and elapsed time) of Bayesian UQ. For verification purposes, we first tested these new capabil ities on a species interaction problem. We then demonstrated them on the TRISO Ag release application, showing that parallel active learning capabilities can enhance the accuracy of UQ while also substantially reducing the computational cost in comparison to the reference methods developed in FY-23.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Comprehensive Dynamics in a Polyelectrolyte Complex Coacervate

The linear viscoelastic response, LVR, of a hydrated polyelectrolyte complex coacervate, PEC, was evaluated over a range of frequencies, temperatures, and salt concentrations. The PEC was a nearly stoichiometric blend of a quaternary ammonium poly([3-(methacrylamido)propyl]trimethylammonium chloride), PMAPTAC, and poly(2-acrylamido-2-methyl-1-propanesulfonic acid sodium salt), PAMPS, an aliphatic sulfonate, selected because they remain fully charged over the conditions of use. Narrow molecular weight distribution polyelectrolytes were prepared using fractionation techniques. A partially deuterated version of PMAPTAC was incorporated to determine its coil radius of gyration, R g , within PECs using small-angle neutron scattering. Chain dimensions were determined to be Gaussian with a Kuhn length of 2.37 nm, which remained constant from 25 to 65 °C. The LVR for a series of matched molecular weight PECs, mostly above the entanglement threshold, exhibited crossovers of modulus versus frequency classically attributed to the reptation time, relaxation between entanglements, and the relaxation of a Kuhn length of units (the “monomer” time). The scaling for zero shear viscosity, η 0 , versus chain length, N, was η 0 ~ N 3.1 , in agreement with “sticky reptation” theory. The lifetime and activation energy, E p , of a pair between polyanion and polycation repeat units, Pol + Pol – , were determined from diffusion coefficients of salt ions within the PEC. The activation energy for LVR of salt-free PECs was 2E p , showing that the key mechanism limiting the dynamics of undoped PECs is pair exchange. An FTIR technique was used to distinguish whether SCN– acts as a counterion or a co-ion within PECs. Doping of PECs with NaSCN breaks Pol + Pol – pairing efficiently, which decreases effective cross-linking and decreases viscosity. Finally, an equation was derived that quantitatively predicts this effect.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impact of Annealing on Superconducting Properties of Tantalum Thin Films

The pursuit for higher-coherence superconducting qubits has led to the adoption of new materials and innovative fabrication techniques. Among the commonly used superconducting materials, tantalum-based or tantalum-capped superconducting qubits have demonstrated the longest lifetimes, outperforming niobium devices even in the presence of their native oxides. Given the losses introduced by the surface oxide, one method to further improve performance is to reduce or eliminate the native oxide on tantalum through UHV annealing. In this work, we examine the effect of dissolved oxygen as a function of annealing temperature on the superconducting properties of tantalum thin films. Dissolved oxygen is quantified using atom probe tomography (APT) and secondary ion mass spectrometry (TOF-SIMS) and correlated with electromagnetic transport measurements performed using a Physical Property Measurement System (PPMS). This study observes that annealing decomposes the native oxide and drives oxygen diffusion into the film matrix. Increasing oxygen content decreases the superconducting transition temperature (Tc) and the residual resistivity ratio (RRR). The critical field decreases with the incorporation of oxygen, while the film’s coherence length decreases relatively by 5%. These results, as well as findings from ongoing studies linking fabrication processes to the structural and chemical properties of the superconducting film and Josephson junction, provide insight into effective strategies for improving the performance of superconducting devices.

Shinde, Maithilee [IIT, Chicago; Fermilab]↗