Search NASA⌕ Search

SEARCH · Search NASA

Results for “Non-equilibrium”

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 145 records · Page 8

Generative Physics-Informed Neural Network Solving Multi-Scale and Multi-Phase Plasma Chemical Flow Field

Low-temperature plasmas (LTPs) are non-equilibrium systems with near-room-temperature gas and highly energetic electrons. This makes them ideal for delicate applications in biomedicine and semiconductor manufacturing, enabling processes like wound healing, sterilization, etching, and plasma-enhanced chemical vapor deposition without thermal damage. However, LTPs involve complex chemistries, with hundreds of species and thousands of reactions, complicating their diagnosis, prediction, and control. Conventional diagnostics, such as Fourier-transform infrared spectroscopy (FTIR), laser-induced fluorescence (LIF), and optical emission spectroscopy (OES), offer limited species detection, while mass spectrometry (MS) struggles with low-sensitivity species. Additionally, LTP simulations face multi-scale challenges, as macroscopic fluid dynamics and microscopic particle collisions operate on vastly different timescales. To address these issues, we developed an artificial intelligence (AI) based diagnostic system: a generative physics-informed neural network (PINN-Gen) that can predict spatially resolved species concentrations and temperatures in LTPs by integrating experimental data from planar LIF with microscopic plasma chemical kinetics and macroscopic fluid mechanics, including plasma-liquid interactions at the interface between two phases. PINN-Gen solves no equations but checks the errors of physical laws by substituting the output from neural network, and the comparison with the experimental results. Thus, it naturally avoids the multi-scale difficulty of numerical simulations and predicts the results of conventionally unsolvable multi-scale and multi-phase problems. The real-time prediction will be robust due to the physical information used in the training of such a neural network, and only very limited input of condition required due to its generative feature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Applying deep learning methods to develop new models of molecular charge transfer, nonadiabatic dynamics, and nonlinear spectroscopy in the condensed phase

Photon- and field-induced charge transfer has central importance in the generation and storage of electricity, the novel properties of materials, photo-induced catalysis, and electro-optic activity (e.g., photovoltaic cells, fuel cells, and organic chromophores for use in optical fibers and light-emission diodes). These non-equilibrium electronic and chemical transformations are probed by ultrafast, nonlinear spectroscopies. Accurate simulations play a crucial role in our ability to understand, optimize, and control these transformations. This project applies modern deep learning and machine learning (ML) methods to dramatically improve models of electronic dynamics, electronic-nuclear dynamics, and spectroscopic measurements for improved simulations of chemistry in complex environments, far from equilibrium phenomena, and processes in extreme environments, such as materials exposed to strong or resonant fields. This project develops accurate neural net models that go beyond predictive capability to also provide new insight into the fundamental physics underlying electron and nuclear dynamics. To achieve its objectives, this project explores and develops customized versions of high-capacity deep learning algorithms/models. These techniques are developed with an emphasis on fundamental chemical insight, not just predictive accuracy, to assist the development of the next generation of quantum simulation methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Enabling Low-Temperature (LTP) Ignition Technologies for Multi-Mode Engines through the Development of a Validated High-Fidelity LTP Model for Predicative Simulations Tools

The goal of multi-mode engine architectures is to extend current lean-burn dilution limits with renewable fuels, which requires spark plugs to deposit high energies (hundreds of mJ) in order to initiate ignition and complete combustion. At elevated energy deposition rates, spark plugs experience increased electrode erosion and thermal losses, which ultimately shortens the spark-plug lifetime and lowers ignition efficiency. As such, in order to safeguard the efficiency gains of multi-mode concepts, new and improved ignition technologies are required. Recently, non-equilibrium low-temperature plasmas (LTP) have been shown to promote energy-efficient ignition via quenching and transport of electronically excited atoms and molecules, selective radical production and fast heating of hydrocarbon/air mixtures [1-2]. Thus, LTP is seen as a technology that can potentially improve the energy extraction efficiency of fuels, while enabling kinetically controlled combustion modes towards fuel leaner conditions to realize current DOE VTO goals of improving the sustainability of future mobility [3]. Although many previous studies have demonstrated the efficacy of plasma-assisted ignition to enhance combustion, the detailed enhancement mechanisms remain largely unknown, especially for oxygenated fuels and at elevated pressures that are most relevant to practical engine conditions. These barriers hinder the development of accurate and comprehensive numerical models that seek to describe LTP-based ignition in existing engine design software tools and methods. Current state-of-the-art simulation capabilities for LTP ignition systems are in need of improvements since they deliver qualitative results only due to important limitations of existing approaches. Firstly, validated kinetic models with elementary steps for plasma discharges in oxygenated fuel/air mixtures of relevance to the transportation sector are required. Such kinetic models do not exist at present and will be developed and validated within this project. Secondly, plasma discharges and reactive mixture ignition are multi-scale, unsteady processes requiring high-performance numerical methods and software that execute efficiently on DOE supercomputers. Such software does not exist at present and will be developed and applied to practical LTP ignition scenarios as part of this project. Thirdly, experimental databases that are tailored to serve as benchmark in support of the development of predictive computational models of LTP ignition do not exist and will be part of this project.

33 ADVANCED PROPULSION SYSTEMS↗

Understanding Mechanisms of Selective Crystal Nucleation and Growth Using In Situ Electron Microscopy

Separation processes are essential for purifying and isolating critical materials, but are challenging due to the natural tendency of systems to mix. This project explores far-from-equilibrium separation methods using external electric and magnetic fields to drive selective nucleation and crystal growth. Using in situ electrochemical transmission electron microscopy (TEM), we observed ion transport and crystallization behavior under applied electric fields. We found that ions migrate toward the working electrode, where nucleation preferentially occurs. These findings advance our understanding of non-equilibrium separation mechanisms and support the development of efficient, sustainable methods for extracting critical materials from unconventional domestic sources.

36 MATERIALS SCIENCE↗

A spline-based method to obtain spatially dependent viscosity in confined flows

Coupling chemical physics to continuum theories is a critical step to understanding multi-scale phenomena. This paper will connect non-equilibrium molecular dynamics simulations to a continuum-based Navier-Stokes equation that has relaxed the assumption of spatial uniformity in viscosity. Using a form for viscosity based on spline interpolation, viscosity as a function of position is obtained from the least squares fit of the velocity profile measured from molecular simulations of flow in a nanochannel. Viscosity can vary widely, particularly near the channel boundaries, indicating that uniform viscosity is no longer appropriate. Variations of the viscosity near the channel surfaces imply that considering solution and surface chemistry could be necessary to rigorously understand molecular-scale flows in nanochannels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Symposium MT02: Statistical Mechanics-Based Computational Tools for the Study of Phase Transformation in Complex Materials (Final Report)

Symposium MT02 brought together a diverse and interdisciplinary community of scientists specializing in Statistical Mechanics-based computational modeling to investigate phase transformations in materials exhibiting complex disordered structures. As the demand for materials with extreme performance metrics grows—from aerospace components to next-generation optical fibers—the ability to predict microstructural evolution under non-equilibrium conditions has become paramount. The primary goal of this symposium was to identify, evaluate, and discuss advanced computational tools capable of designing precise manufacturing conditions to tailor material properties efficiently. By fostering a dialogue between computational theorists and experimentalists, the symposium sought to establish new protocols for predicting how processing history—such as cooling rates or strain paths—dictates the final microstructure.

36 MATERIALS SCIENCE↗

Demonstrating magnetic field robustness and reducing temporal T1 noise in transmon qubits through magnetic field engineering

The coherence of superconducting transmon qubits is often disrupted by fluctuations in the energy relaxation time (T1), limiting their performance for quantum computing. While background magnetic fields can be harmful to superconducting devices, we demonstrate that both trapped magnetic flux and externally applied static magnetic fields can suppress temporal fluctuations in T1 without significantly degrading its average value or qubit frequency. Using a three-axis Helmholtz coil system, we applied calibrated magnetic fields perpendicular to the qubit plane during cooldown and operation. Remarkably, transmon qubits based on tantalum-capped niobium (Nb/Ta) capacitive pads and aluminum-based Josephson junctions (JJs) maintained T1 lifetimes near 300 μs even when cooled in fields as high as 600 mG. Both trapped flux up to 600 mG and applied fields up to 400 mG reduced T1 fluctuations by more than a factor of two, while higher field strengths caused rapid coherence degradation. We attribute this stabilization to the polarization of paramagnetic impurities, the role of trapped flux as a sink for non-equilibrium quasiparticles (QPs), and partial saturation of fluctuating two-level systems (TLSs). These findings challenge the conventional view that magnetic fields are inherently detrimental and introduce a strategy for mitigating noise in superconducting qubits, offering a practical path toward more stable and scalable quantum systems.

Abdisatarov, Bektur [Fermilab; Jefferson Lab]↗

The remarkable inefficiency of stratocumulus

Marine stratocumulus clouds play a central role in Earth's climate system by reflecting incoming solar radiation and exerting a strong cooling effect. Their organization into open and closed mesoscale cellular morphologies can be thought of as an example of bistable dynamics driven by aerosol–cloud interactions and mesoscale processes. From the perspective of non-equilibrium thermodynamics, these structures are an example of a far-from-equilibrium open system that continuously produces and exports entropy. While entropy production has been studied in idealized deep convective systems, it has not yet been quantified for shallow clouds. Here, we compute and decompose the internal entropy production of open- and closed-cell stratocumulus using an ensemble of large-eddy simulations. We show that the overall entropy production of stratocumulus is low, reflecting the limited vertical extent and corresponding reduced ability to utilize the energy fluxes at the system's boundaries. Moist processes dominate the overall irreversibility, which, combined with their low entropy production, leads to a mechanical efficiency about an order of magnitude smaller than in deep convective systems. Although the dominant irreversible processes differ between open- and closed-cell regimes, the distributions of total entropy production largely overlap across the ensemble, limiting the ability to distinguish the dynamics of individual cases based solely on total entropy production.

54 ENVIRONMENTAL SCIENCES↗

Technical note: A modified formulation of dynamic energy budget theory for faster computation of biological growth

Abstract. The mass conservation equation in the presence of boundary fluxes and chemical reactions from non-equilibrium thermodynamics is used to derive a modified dynamic energy budget (mDEB) model. Compared to the standard dynamic energy budget (sDEB) model (Kooijman, 2009), this modified formulation does not place the dilution effect in the mobilization kinetics of reserve biomass, and it maintains the partition principle for reserve mobilization dynamics for both linear and non-linear kinetics. Overall, the mDEB model shares most features with the sDEB model. However, for biological growth that requires multiple nutrients, the mDEB model is computationally much more efficient by not requiring numerical iterations for obtaining the specific growth rate. In an example of modeling the growth of Thalassiosira weissflogii in a nitrogen-limiting chemostat, the mDEB model was found to have almost the same accuracy as the sDEB model while requiring almost half of the computing time of the sDEB model. Since the sDEB model has been successfully applied in numerous studies, we believe that the mDEB model can help improve the modeling of biological growth and the associated ecosystem processes in various contexts.

Tang, Jinyun↗

RTE: A computer code for Rocket Thermal Evaluation

The numerical model for a rocket thermal analysis code (RTE) is discussed. RTE is a comprehensive thermal analysis code for thermal analysis of regeneratively cooled rocket engines. The input to the code consists of the composition of fuel/oxidant mixture and flow rates, chamber pressure, coolant temperature and pressure. dimensions of the engine, materials and the number of nodes in different parts of the engine. The code allows for temperature variation in axial, radial and circumferential directions. By implementing an iterative scheme, it provides nodal temperature distribution, rates of heat transfer, hot gas and coolant thermal and transport properties. The fuel/oxidant mixture ratio can be varied along the thrust chamber. This feature allows the user to incorporate a non-equilibrium model or an energy release model for the hot-gas-side. The user has the option of bypassing the hot-gas-side calculations and directly inputting the gas-side fluxes. This feature is used to link RTE to a boundary layer module for the hot-gas-side heat flux calculations.

Naraghi, Mohammad H. N.↗

A Technique for Determining the Nozzle-Flow Properties of Air in an Equilibrium, Nonequilibrium, or Frozen State

One of the problems associated with the development of high-enthalpy hypersonic test facilities hos been the determination of nozzle flow properties. This note presents a technique for determining the flow properties of air in an equilibrium, non-equilibrium, or frozen state using two test section measurements. A knowledge of two stagnation properties is required before the technique can be applied. The stagnation pressure can be measured, and reference 1 offers a method for determining stagnation enthalpy. Results of the calculations using the method of reference 1 are presented in reference 2 in a chart (chart 21) which can be readily used to determine the stagnation enthalpy from measurements of stagnation pressure, mass flow of air through the tunnel, and throat diameter.

Hayman, Lovick O.↗

Trajectory Engineering with Modular Patched Conics for Entry Systems and TPS (TEMPEST)

Brief Presenter Biography (35 word limit): Bohdan Wesely is an Aerospace Engineer in the Entry Systems and Technology Division at Ames. He has worked on a variety of projects for NASA including integrated TPS (thermal protection system) flight hardware deliveries and testing services for commercial partners. Introduction: TEMPEST is a new trajectory analysis framework that is designed to fill the gap between dedicated flight mechanics tools and aerothermal and TPS sizing tools. The project started as an SJSU master’s thesis and has since evolved into a general conceptual design tool capable of studying a wide variety of entry problems. Development is ongoing in the Entry Systems and Technology Division at NASA ARC. Why TEMPEST: Space missions involving entry into a planetary atmosphere involve a series of unique requirements across multiple disciplines. Whether it is traditional entry descent and landing (EDL), or aerocapture, the vehicle must navigate to its target landing location or orbit state, and the TPS must protect the payload during entry. The design process typically involves iterative handoffs between various flight mechanics, flow solver, and material response level tools. During the early conceptual phase, a wide variety of feasible trajectories are simulated in a Monte Carlo scenario which broadly satisfy the mission or landing requirements. Next, computational fluid dynamics (CFD), direct simulation Monte Carlo (DSMC), and other flow solver analyses are performed at various key trajectory points to generate an aero-database, heating and TPS design requirements also emerge at this stage. At this point, with updated aerodynamics from the various flow solvers, trajectories can be re-run, this in turn can change the required freestream conditions for the CFD tools, and as a project progresses, these analyses converge, and uncertainty is reduced. However, there is always a “hand-off” occurring between two inherently coupled phenomena. Analysis Description: One of the goals with TEMPEST is to use a variety of first principles estimation methods coupled with an atmosphere model to predict vehicle aerothermodynamics across the entire flight regime while propagating a 3 or 6 degree of freedom (DoF) trajectory. Aerodynamics methods include modified Newtonian, Maxwell and Cercignani- Lampis-Lord (CLL) for continuum, transitional, and free molecular flow regimes. Aerothermodynamics include boundary layer and reference enthalpy methods, and Mutation++ for non-equilibrium chemistry modeling. TEMPEST is also capable of stitching multiple trajectory segments together to study mission scenarios like multi-pass aerocapture and aero-gravity assists. Most of the program is implemented in MATLAB using modern system objects, it relies on several C++ shared libraries for supporting tools like Gmsh, the Global Reference Atmospheric Model (GRAM), and Mutation++. The various first principles aerothermal estimation methods are discretized across either a structured axisymmetric panel mesh or an unstructured tri-mesh generated from an open-source tool such as Gmsh, this allows solutions on the same mesh to be compared across tools such as CB-Aero. CFD Coupling. A physics-aware, gaussian process CFD anchoring scheme is proposed to adjust the various first principles methods as a CFD database is populated. One goal for this anchoring module is to inform the project where CFD should be run. Full knowledge of the entire trajectory, atmosphere, and aerothermodynamics allows for easier identification of high sensitivity areas and uncertainty quantification. While the first principles effects are well known and proven accurate in existing tools such as CB- Aero and Cart3D, a physics aware CFD anchoring scheme increases tool credibility across a project lifecycle. Material Response Modeling. Correct TPS sizing is critical for optimizing mass for science payloads and ensuring mission success. The process typically involves a thermal analysis along the trajectory with surface heating environments as a boundary condition. Several design constraints are maximum bondline temperature and maximum recession with various margining techniques. The material response tool FIAT, developed out of NASA Ames, is currently being integrated into the TEMPEST environment. TPS recession, shape change, mass loss, and mass property alteration are all factors that can perturb an entry trajectory. For missions like Mars 2020, recession was minimal and was safely handled separately as a post process. For missions such as Jupiter Galileo with a high TPS mass fraction or asteroid entries, recession plays a major role. The proposed fully coupled scheme is to use an epoch-based approach where the trajectory integration is halted after a recession threshold, the energy balance and FIAT are solved at each panel, the mesh, aerodynamics, and mass properties are updated, and the trajectory continues. Several computational tradeoffs have been made during the development of TEMPEST to limit the cost of a single trajectory and preserve its utility as a conceptual, rapid iteration tool. Conclusion: Development of TEMPEST is ongoing and the project is still in its infancy. This talk aims to showcase its unique capabilities to support future NASA entry systems missions.

Bohdan O Wesely↗

Thermochemically-Closed Sonic-Flow Inversion for Enthalpy and Temperature in Multispecies Arc-Jet Flows

A thermochemically-closed sonic-flow inversion framework (TSIF) is developed to infer bulk enthalpy and total temperature upstream of a choked nozzle in arc-jet flows. The formulation recasts a pressure-rise total enthalpy quantification technique as an inverse problem in characteristic-velocity c * space using measured mass flow rate, upstream total pressure, gas composition, and nozzle throat geometry as inputs. Unlike calorimetric energy-balance approaches or optical diagnostics, the method relies primarily on routinely measured facility quantities combined with explicit thermochemical closure. Thermochemical states are obtained using NASA’s open-source Chemical Equilibrium with Applications (CEA) code, enabling construction of a chemistry-consistent relation between characteristic velocity, total enthalpy, and total temperature under equilibrium or frozen assumptions. A discharge coefficient is self-calibrated using cold-flow (arc-off) operation data and applied to hot-flow (arc-on) measurements, enabling upstream losses to be accounted for without empirical correlations. The framework is applied to air, N 2 , and CO 2 –N 2 arc-jet flows and demonstrates expected trends for the inferred thermochemical states as function of arc power, specific energy input, mass-flow, heater configuration, and test gas. In the air limit, under equilibrium assumptions, the method recovers the classical high-enthalpy asymptotic correlation of Winovich with a mean residual of 4.4%, demonstrating compatibility with established sonic-flow scaling, while extending applicability to arbitrary multi-species mixtures and non-equilibrium chemistry. The framework provides a mixture-flexible methodology for determining bulk thermochemical states in modern arc-jet environments using routine facility pressure, mass-flow, gas-composition, and nozzle-geometry information together with a cold-flow calibration.

stagnation heat flux↗

A Machine Learning Framework for Error Compensation in Radiative Transfer Calculations

Radiative heat transfer influences the amount of heat flux transferred to the surface of the hypersonic vehicle, which is essential to evaluate the performance of thermal protection systems. The radiative heat flux is found to be computationally prohibitive while accounting for the variation in spatial, angular, and spectral domains. A new methodology has been recently developed to alleviate the cost of computation in the spectral domain by constructing flow-agnostic reduced-order models (ROMs). The developed spectral ROM databases provide grouping strategies that account for non-equilibrium absorption and emission as well as interaction between disparate species due to spectral overlap in associated radiative processes. However, the developed ROMs need to be optimized for a specific combination of interacting gas species and would need to re-calibrated in case individual species are added/omitted. In this work, we use various machine learning (ML) techniques to approximate the radiative intensities determined by a ROM optimized for a specific gas mixture. The ML model relies on the ROM databases developed for a single species which ignores any spectral overlap. Thus, radiation evaluation starts with a simple summation of radiative intensities predicted using these non-calibrated ROMs for the contributing species. The ML framework then provides a correction to account for the interplay in the frequency, i.e., emission of photons by one species and absorption by another, and yields mixture-specific radiation fields. Once trained on the individual ROM databases, the ML framework offers instantaneous corrections that serves as a time/cost effective alternative to the optimization of ROMs for a specific gas mixture. The ML framework is trained on both the high fidelity and ROM evaluated line of sight (LOS) data from Orion, Stardust, and FIRE II cases to obtain a general purpose correction model for earth re-entry scenarios when radiation contributions from both atomic nitrogen and atomic oxygen are considered. A geometric length scale parameter is used in the training process to account for errors introduced in the ROM databases as a consequence of high optical thickness. The efficacy of the ML framework is underscored through extensive analysis of train and test errors with respect to all the re-entry scenarios. The applicability of such an ML framework was further corroborated by embedding it in a state-of-the-art US3D - NERO system for determining the radiative heat flux transferred to the hypersonic vehicle surface.

Radiation↗

Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo

Accurate prediction of aerothermodynamic loads in thermal non-equilibrium flows requires precise modeling of internal energy exchange. While previous direct simulation Monte Carlo frameworks have successfully implemented discrete rotational energy models for diatomic species, the treatment of polyatomic molecules has traditionally relied on continuous energy assumptions that break down at low temperatures and neglect critical high-temperature corrections. This study extends the discrete rotational energy models of Boyd and Gimelshein to fully encompass polyatomic molecules. The proposed framework implements quantized rotational energy level sampling for linear, spherical, and symmetric/asymmetric top rotors. Crucially, the model incorporates centrifugal distortion to address the limitations of the rigid-rotor assumption at hypersonic temperatures, and accounts for nuclear spin parity, which dictates the permissible rotational states and macroscopic specific heats at low temperatures. The model is verified through equilibrium sampling procedures, demonstrating agreement with theoretical quantum Boltzmann distributions and accurately reproducing thermophysical properties across a wide range of temperatures.

DSMC↗

Thermochemically-Closed Sonic-Flow Inversion for Enthalpy and Temperature in Multispecies Arc-Jet Flows

A thermochemically-closed sonic-flow inversion framework (TSIF) is developed to infer bulk enthalpy and total temperature upstream of a choked nozzle in arc-jet flows. The formulation recasts a pressure-rise total enthalpy quantification technique as an inverse problem in characteristic-velocity c * space using measured mass flow rate, upstream total pressure, gas composition, and nozzle throat geometry as inputs. Unlike calorimetric energy-balance approaches or optical diagnostics, the method relies primarily on routinely measured facility quantities combined with explicit thermochemical closure. Thermochemical states are obtained using NASA’s open-source Chemical Equilibrium with Applications (CEA) code, enabling construction of a chemistry-consistent relation between characteristic velocity, total enthalpy, and total temperature under equilibrium or frozen assumptions. A discharge coefficient is self-calibrated using cold-flow (arc-off) operation data and applied to hot-flow (arc-on) measurements, enabling upstream losses to be accounted for without empirical correlations. The framework is applied to air, N 2 , and CO 2 –N 2 arc-jet flows and demonstrates expected trends for the inferred thermochemical states as function of arc power, specific energy input, mass-flow, heater configuration, and test gas. In the air limit, under equilibrium assumptions, the method recovers the classical high-enthalpy asymptotic correlation of Winovich with a mean residual of 4.4%, demonstrating compatibility with established sonic-flow scaling, while extending applicability to arbitrary multi-species mixtures and non-equilibrium chemistry. The framework provides a mixture-flexible methodology for determining bulk thermochemical states in modern arc-jet environments using routine facility pressure, mass-flow, gas-composition, and nozzle-geometry information together with a cold-flow calibration.

inviscid theory↗

Mesoscale Modeling of Dislocation Cell Structure Evolution and Radiation-Induced Segregation in Additively Manufactured Austenitic Stainless Steel

Structural alloys under irradiation develop radiation-induced segregation (RIS) at point defect sinks, which lead to undesired changes in the alloy's properties. Additively manufactured austenitic stainless steels are expected to yield a distinct response to irradiation damage owing to their unique as-printed dislocation cell structure. In this talk, we present the development of a mesoscale model for dislocation cell structure evolution coupled with RIS using the MOOSE framework. The evolution of network dislocations within spatially distinct cells and cell walls are modeled using climb-mediated edge annihilation and generation processes. Starting with initial microstructures of as-printed non-equilibrium segregations at dislocation cell walls, our 2D and 3D simulations explore the effects of irradiation temperature, dose rate, dislocation sink bias, and specimen thickness on microstructure evolution. The results will be compared against experimental characterizations of in-situ and ex-situ ion-irradiated samples, and the implications of evolving cell structure on irradiation damage response will be discussed.

additive manufacturing↗

Vortex-Controlled Quasiparticle Multiplication and Self-Growth Dynamics in Superconducting Resonators

Even in the quantum limit, non-equilibrium quasiparticle (QP) populations induce QP poisoning that irreversibly relaxes the quantum state and significantly degrades the coherence of transmon qubits. A particularly detrimental yet previously unexplored mechanism arises from QP multiplication facilitated by vortex trapping in superconducting quantum circuits, where a high-energy QP relaxes by breaking additional Cooper pairs and amplifying the QP population due to the locally reduced excitation gap and enhanced quantum confinement within the vortex core. Here we directly resolve this elusive QP multiplication process by revealing vortex-controlled QP self-generation in a highly nonequilibrium regime preceding the phonon bottleneck of QP relaxation. At sufficiently low fluence, femtosecond-resolved magneto-reflection spectroscopy directly reveals a continuously increasing QP population that is strongly dependent on magnetic-field-tuned vortex density and absent at higher excitation fluences. Quantitative analysis of the emergent QP pre-bottleneck dynamics further reveals that, although the phonon population saturates within $\simeq$10~ps, both free and trapped QPs continue to grow in a self-sustained manner--hallmarks of the long-anticipated QP-vortex interactions in nonequilibrium superconductivity. We estimate a substantial increase of $\sim$34% in QP density at vortex densities of $\sim$ 100 magnetic flux quanta per $\mathrm{μm^{2}}$. Our findings establish a powerful spectroscopic tool for uncovering QP multiplication and reveal vortex-assisted QP relaxation as a critical materials bottleneck whose mitigation will be essential for resolving QP poisoning and enhancing coherence in superconducting qubits.

Park, Joong M. [Ames Lab]↗