SEARCH · Engineering Papers

Results for “dynamic insulation”

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

SourceResearch connectionFrom yearThrough year

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

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Ultrafast laser-driven many-body dynamics and Kondo coherence collapse

Ultrafast laser pulse has provided a systematic way to inspect the dynamics of electrons in condensed matter systems. In this paper, by means of time-dependent density matrix renormalization group, we study an ultrafast laser-driven Kondo lattice model, in which conduction electrons are strongly coupled with magnetically local moments. The single-particle spectral function due to strong correlation effects and photon emission in the nonequilibrium states under laser driving are calculated. We find laser field excited collective doublon-hole pairs and an associated transient melting of Kondo coherence phase, signifying the collapse of Kondo energy gap. Moreover, we show that the photon emission, induced by a strong laser field, exhibits a different intensity characteristics than in the equilibrium Kondo insulator, which could be explained by the Kondo collapse and related suppression of both intraband and interband contribution in Kondo melting liquid. This theoretical insight is accessible with time- and angle-resolved photoemission spectroscopy and high-harmonic generation spectroscopy and will stimulate the investigation of nonequilibrium dynamics and nonlinear phenomenon in heavy fermion systems.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Hydrodynamic behavior of wrapped coil water heater tank

In this work, a numerical investigation using computational fluid dynamics (CFD) has been carried out to analyze the thermal behavior of a water tank integrated to a heat pump water heating system (HPWH). The heat pump condenser is a long coil wrapped around the water tank and insulated from outside from which heat is transferred from the heat pump cycle into the water stored inside the tank. The simulated water tank has a capacity of about 170 L with the charging inlet located at the bottom segment of the tank and the discharging outlet located at the top segment of the tank like commercially available water heaters. The tank was simulated during charging/discharging process with isothermal boundary condition at the tank wall to represent the majority of the process for refrigerant condensation inside the heat pump condenser (excluding superheated and sub-cooled section). Three flow rates were selected to represent low, medium and high flow pattern. Both mean and turbulent flow properties were calculated and analyzed. Heat transfer mechanisms inside the tank are quantified by calculating Richardson and Péclet numbers to show the relative contribution of driving forces. The calculations show that advection is dominant over diffusion. Mixing factor is calculated and presented. It is found that the rate of increase of mixing is nonlinear with the increase in flow rate. Mixing increased by a factor of four when the flow rate increased by a factor of two.

stratification↗

DOE OSTI2022

Ultrafast signatures of spin and orbital order in antiferromagnetic α-Sr2CrO4

Abstract The antiferromagnetic Mott insulator α -Sr 2 CrO 4 possesses multiple spin and orbital ordered phases, but their unique interplay is still relatively unexplored. Here, we used femtosecond optical spectroscopy to study ultrafast spin and orbital ordering dynamics in α -Sr 2 CrO 4 through their non-equilibrium response to photoexcitation. By varying the pump photon energy, we selectively drove inter-site spin hopping between neighboring Cr t 2 g orbitals and charge transfer-type transitions between oxygen 2 p and Cr e g orbitals. The resulting transient reflectivity dynamics revealed temperature-dependent anomalies across the Néel temperature for spin ordering as well as the transition temperatures linked to different types of orbital order. Our results reveal distinct relaxation timescales for spin and orbital orders in α -Sr 2 CrO 4 and provide experimental evidence for the phase transition at T O , possibly related to antiferro-type orbital ordering.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Tensor train continuous time solver for quantum impurity models

The simulation of strongly correlated quantum impurity models is a significant challenge in modern condensed matter physics that has multiple important applications. Thus far, the most successful methods for approaching this challenge involve Monte Carlo techniques that accurately and reliably sample perturbative expansions to any order. However, the cost of obtaining high precision through these methods is high. Recently, tensor train decomposition techniques have been developed as an alternative to Monte Carlo integration. In this study, we apply these techniques to the single-impurity Anderson model at equilibrium by calculating the systematic expansion in power of the hybridization of the impurity with the bath. Furthermore, we demonstrate the performance of the method in a paradigmatic application, examining the first-order phase transition on the infinite-dimensional Bethe lattice, which can be mapped to an impurity model through dynamical mean field theory. Our results indicate that using tensor train decomposition schemes allows the calculation of finite-temperature Green's functions and thermodynamic observables with unprecedented accuracy. The methodology holds promise for future applications to frustrated multiorbital systems, using a combination of partially summed series with other techniques pioneered in diagrammatic and continuous time quantum Monte Carlo.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Origin of metal-insulator transitions in the parent compounds of ruthenium-pnictide superconductors

Here we study the interplay of the structural phase transition, flat electronic band dispersion, and metal-to-insulator transition (MIT) in the parent compounds of the Ru-pnictide superconductors by using first-principles calculations. Our electron and phonon calculations reveal that Ru(P,As) undergo MIT accompanied by orthorhombic to monoclinic distortion at low temperature, but RuSb stays orthorhombic and metallic in agreement with the experimental findings. It is demonstrated that electronic correlation, as treated in DFT + U, DFT + Gutzwiller, and dynamical mean-field theory (DMFT), cannot induce MIT in the undistorted crystal structure. We find that, although small monoclinic distortion can remove the van Hove singularity at the Fermi level, it does not immediately gap out the Fermi surface and a large value of monoclinic distortion is necessary for a clear MIT, suggesting the possibility of an intermediate pseudogapped monoclinic metallic phase. Furthermore, we predict a light-induced two-step insulator-to-metal and structural transitions in the monoclinic phases of RuP and RuAs, which can be tested in future ultrafast pump-probe experiments as an alternative ideal playground to VO 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Unraveling the multi-step crystallization mechanism of polytetrafluoroethylene, modified polytetrafluoroethylene, and their nanocomposites with boron nitride nanobarbs: Experimental insights and theoretical analysis

The non-isothermal crystallization behavior and kinetics of polytetrafluoroethylene (PTFE) composites with boron nitride nanobarb (BNNB), a new generation nanostructure with unique surface morphology and mechanical “barbs” have been analyzed, understanding these properties is essential for their high-end applications as thermal interface materials (TIM) for microwave, 5G and microelectronic devices. The analysis of the crystallization parameters includes crystallization onset, peak and end temperatures, crystallization half-life and overall crystallinity of PTFE, modified PTFE and their BNNB composites. The results were further analyzed using theoretical models such as the combined Avrami-Ozawa model. It was found that BNNB supports crystallization in the modified PTFE but shows minimal effect on the crystallization of PTFE. Due to the limitation of the classical theoretical models used above in fully characterizing the multi-step crystallization process of PTFE, an in-depth analysis using the model-free advanced isoconversional computation was used to characterize the PTFE crystallization based on the evolution of activation energy with fractional crystallinity and for the first time with temperature. Three kinetic regions were identified in the crystallization mechanism. Here, this study investigated the molecular organization and microstructural evolution of PTFE, modified PTFE and their composites during non-isothermal crystallization using advanced X-ray scattering measurements. An insight into the changes undergone by the material's microstructural units including crystallite size and morphology, lamellar thickness and lamellar interfacial layer thickness, and crystallographic phase dynamics during non-isothermal cooling from the melt, was provided here in this work. The effect of copolymer modification of PTFE and the inclusion of pristine and functionalized BNNB (a thermally conductive and electrically insulating ceramic) are both new investigations that provide valuable knowledge for the development of materials with strong matrix-nanofiller interaction and guidance for optimizing sintering and cooling cycles, two key steps in PTFE processing that largely affect the material microstructural features. Overall, the result of the three-part investigation demonstrates that BNNB supports crystallization in the modified PTFE up to 20 wt% concentration and at low and high cooling rates typically used in the industrial processing of PTFE.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Discerning element and site-specific fluctuations of the charge-orbital order in Fe 3 O 4 below the Verwey transition

Despite countless experimental probes into magnetite's electronic structure across the Verwey transition Fe 3 O 4 , the exact origin of this archetypical metal-insulator transition remains a puzzle. Advanced x-ray diffraction techniques have mostly resolved the monoclinic structure of the insulating phase, including interatomic bond lengths, but the complexity of the charge-orbitally ordered state is difficult to disentangle. Here we combined resonant elastic x-ray scattering and x-ray photon correlation spectroscopy to probe charge-orbital fluctuations in the insulating state of magnetite. By accessing the Bragg forbidden $(00\frac{1}{2})_{c}$ peak at the oxygen K-edge, we complement our previous study on the iron L 3 – edge to reveal the dynamics of the iron 3d and oxygen 2p orbital domains. Our new results reveal a decoupling of the orbital correlation lengths between the oxygen 2p states and site-specific iron 3d states, and we further show charge-orbital domain fluctuations at the iron t 2g orbital sites of trimeron chains. These results also demonstrate an experimental method capable of distinguishing electronic dynamics between the oxygen ligands and the transition metal that underpins emergent behaviors in complex oxides.

36 MATERIALS SCIENCE↗

DOE OSTI2022

On the initiation and evolution of dielectric breakdown in auto-magnetizing liner experiments

Auto-magnetizing (AutoMag) liners are cylindrical tubes composed of discrete metallic helices encapsulated in insulating material; when driven with a ~2 MA, ~100-ns prepulse on the 20 MA, 100-ns rise time Z accelerator, AutoMag targets produced >150 T internal axial magnetic fields. Once the current rise rate of the pulsed power driver reaches sufficient magnitude, the induced electric fields in the liner cause dielectric breakdown of the insulator material and, with sufficient current, the cylindrical target radially implodes. The dielectric breakdown process of the insulating material in AutoMag liners has been studied in experiments on the 500–900 kA, ~100-ns rise time Mykonos accelerator. Multi-frame gated imaging enabled the first time-resolved observations of photoemission from dynamically evolving plasma distributions during the breakdown process in AutoMag targets. Using magnetohydrodynamic simulations, we calculate the induced electric field distribution and provide a detailed comparison to the experimental data. We find that breakdown in AutoMag targets does not primarily depend on the induced electric field in the gaps between conductive helices as previously thought. Finally, to better control the dielectric breakdown time, a 12–32 mJ, 170 ps ultraviolet (λ = 266 nm) laser was implemented to irradiate the outer surface of AutoMag targets to promote breakdown in a controlled manner at a lower internal axial field. Here, the laser had an observable effect on the time of breakdown and subsequent plasma evolution, indicating that pulsed UV lasers can be used to control breakdown timing in AutoMag.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Second-Order Josephson Effect in Excitonic Insulators

We show that in electron-hole bilayers with excitonic orders arising from conduction and valence bands formed by atomic orbitals that have different parities, nonzero interlayer tunneling leads to a second-order Josephson effect. This means the interlayer electrical current is related to the phase of the excitonic order parameter as J = J c sin 2θ instead of J = J c sin θ and that the system has two degenerate ground states at θ = 0, π that can be switched by an interlayer voltage pulse. When generalized to a three dimensional stack of alternating electron-hole planes or a two dimensional stack of chains, the ac Josephson effect implies that electric field pulses perpendicular to the layers and chains can steer the order parameter phase between the two degenerate ground states, making these devices ultrafast memories. Here, the order parameter steering also applies to the excitonic insulator candidate Ta 2 NiSe 5 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Development of the INL Thermal Energy Distribution System (TEDS) in the Modelica Eco-System for Validation and Verification

This report provides an overview on the development of the Thermal Energy Distribution System (TEDS) in the Modelica process model eco-system at Idaho National Laboratory as part of the integrated energy system initiative. Model development has led to the creation of a dynamic process model of the experimental TEDS facility housed within the Engineering Science Laboratory (ESL) at INL. The model was then used during the pre-construction phase of the experimental effort to inform experimental design (insulation requirements, bypass line placement, expected performance of components) and to test innovative control schemes prior to the initial operation. The TEDS model developed in Modelica includes the primary components of the TEDS experimental unit including: a 200kW Chromalox heater, a single-tank packed bed thermal energy storage system filled with 0.125 inch Alumina (Al 2 O 3 ), an ethylene-glycol to Therminol-66 heat exchanger, system piping, six control valves, and all associated temperature, pressure, and mass flow sensors. The model does not include nitrogen fill gas tanks, or associated overfill tanks, as these are not part of standard system control. Two simulation sets were run. The first is a five-hour test that operates as a shakedown test for the facility. This simulation puts the facility through all five potential operating modes and showcases the ability of valving, control sensors, and component controllers to meet the system demands. The second case imposes a typical summer day demand on the system from a region with mixed commercial and residential electrical needs where the generator alone cannot meet peak demand but instead requires the thermal storage unit to act as a peaking unit. Over the course of the fifteen-hour discharge cycle heat losses account for a heat loss of approximately 15°C in outlet temperature. Such heat losses would lead to a decrease in overall system efficiency and may lead to additional design changes for such designs in the future that may include a topping heater. Through commencement of this work, a systems level model of TEDS with associated control systems, sensors, piping diameters, and component capabilities has been created. This model has been utilized in the pre-experimental phase to inform system design, insulation thicknesses, and potential control schemes to operate the system effectively and safely. Once TEDS begins operation, this model will be refined, tuned, and used for validation and verification purposes.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2023

Tuning anisotropic bonding via chemistry and pressure in layered pnictides and chalcogenides (Final Report)

The overarching goal of this research program, as originally delineated in the proposal “Tuning anisotropic bonding via chemistry and pressure in layered pnictides and chalcogenides” was to develop a predictive, chemistry-driven understanding of the impact of the phonon behavior on thermal properties of bulk layered materials. At finite temperatures, atomic vibrations (phonons) strongly impact the thermodynamics, thermal and electrical transport, and phase-switching properties of functional materials. In particular, soft phonon modes and strongly anharmonic potentials can have spectacular consequences, including structural phase transitions (for example in ferroelectrics and phase-change memory materials), metal-insulator transitions, and extreme thermal resistance preventing heat propagation. Bulk materials with highly-anisotropic bonding may provide unique strategies to induce soft-phonon modes and lattice instabilities. Recently, increasingly detailed investigations of the lattice dynamics in layered materials have been made possible by the advent of first-principles phonon calculations and advanced characterization techniques based on neutron and X-ray scattering. However, due to the lack of studies in which composition and bonding character are systematically varied, there are still fundamental questions regarding the impacts of anisotropic bonding and anharmonicity on lattice stability and thermal transport. One of the major goals of this research program is therefore to address this gap by coherently tuning bonding anisotropy and anharmonicity across families of related compounds. Such approaches have revealed new strategies for exploiting structural anisotropy in quasi-1D and 2D bulk materials to obtain tailored functional properties. This project systematically explored the lattice dynamics, phase stability, and transport properties in bulk layered materials by using both composition and applied pressure to tune the degree of bonding anisotropy and anharmonicity. To accomplish this work, we combined i) single-crystal growth of key material systems with tunable anisotropy, ii) in-situ high-temperature/high-pressure characterization of structure and phonons to probe bonding anisotropy and anharmonicity, including state-of-the-art inelastic X-ray scattering (IXS) and inelastic neutron scattering (INS), and iii) first-principles simulations leveraging large-scale computing to identify the fundamental origins of the observed effects, by relating atomic structure and dynamics to electronic orbital interactions. Finally, we modeled and verified the impact of the phonon behavior on thermal transport to identify new strategies for a-priori design of thermal conductivity. Our integrated collaborative approach helped to systematically unravel the effects of anisotropy and bonding anharmonicity on phonon transport, thermodynamics, and thermal properties of complex anisotropic materials.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2022

Tuning anisotropic bonding via chemistry and pressure in layered pnictides and chalcogenides. Final Report

36 MATERIALS SCIENCE↗

DOE OSTI2025

Applications of Nickelate perovskites for neuromorphic computing from electronic structure and Machine Learning

While the limit of Moore's law is presently being reached with current microelectronic technologies, we need to develop new paradigms that overcome this limitation. In that respect, neuromorphic computing is a concept that emulates the neural behavior and response of the human brain, and it has been recognized as a promising alternative approach. In this research project, we will perform multi-fidelity scale bridging to explore the potential use of materials with metal to insulator transition for neuromorphic applications. In particular, rare earth nickelates are promising for such purposes, as the transition in these materials is quite sensitive to a broad set of different external stimuli. Our multi-fidelity approach will bridge the high-fidelity electronic structure calculations with classical potentials. We will bridge dynamical mean field theory with a classical atomistic representation via a deep learning force field. The neural network is trained with energies, charges, and forces obtained by accurate electronic structure theories based on Dynamical Mean Field Theory. The configurational space is generated from known crystal phases, ab initio molecular dynamics with exchange-correlation functionals corrected with the Hubbard model, disordered phases with different concentrations of oxygen vacancies, and nonsymmetrical positions and induced strain by grain interfaces or contact with a substrate. Strategies to train the model with a reduced number of training examples are obtained from active learning methods, and new structures for improving the learning process are generated by using machine learning autoencoders. This classical potential will be validated through a diversity of electronic structure methods and represents an important step to combine the flexibility and accuracy of first-principles with the speed of classical potentials. The generated multi-fidelity surrogate model will be used to understand the role of strain, oxygen vacancies, proton doping, the variation of the crystal phase, substrate effects, vibrational effects as the octahedral rotation, grain boundaries and defect effects on the response of a Metal to Insulator Transition (MIT) in correlated materials. Long time and large-scale simulations will help understand the role of different stimuli to control the hysteresis of the MIT, as it has been experimentally suggested. Selected configurations will be analyzed with higher-level theories to provide an accurate electronic description and to study how the orbitals and charges are rearranged under different conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Imaging of Electrothermal Filament Formation in a Mott Insulator

Resistive switching—the current- and voltage-induced change of electrical resistance—is at the core of memristive devices, which play an essential role in the emerging field of neuromorphic computing. This study is about resistive switching in a Mott insulator, which undergoes a thermally driven metal-to-insulator transition. Two distinct switching mechanisms are reported for such a system: electric-field-driven resistive switching and electrothermal resistive switching. The latter results from an instability caused by Joule heating. Here, we present the visualization of the reversible resistive switching in a planar V 2 O 3 thin-film device using high-resolution wide-field microscopy in combination with electric transport measurements. We investigate the interaction of the electrothermal instability with the strain-induced spontaneous phase separation in the V 2 O 3 thin film at the Mott transition. The photomicrographs show the formation of a narrow metallic filament with a minimum width ≲500 nm. Although the filament formation and the overall shape of the current-voltage characteristics (IVCs) are typical of an electrothermal breakdown, we also observe atypical effects such as oblique filaments, filament splitting, and hysteretic IVCs with sawtoothlike jumps at high currents in the low-resistance regime. We are able to reproduce the experimental results in a numerical model based on a two-dimensional resistor network. This model demonstrates that resistive switching in this case is indeed electrothermal and that the intrinsic heterogeneity is responsible for the atypical effects. This heterogeneity is strongly influenced by strain, thereby establishing a link between switching dynamics and structural properties.

42 ENGINEERING↗

DOE OSTI2026

Terahertz‐Nanoscale Visualization of the Microscopic Spin‐Charge Architecture of Colossal Magnetoresistive Switching

Resolving sub-10 nm spin switching and the associated terahertz (THz) electrodynamics during the colossal magnetoresistance (CMR) transition is a definitive frontier in reaching the fundamental spatial, temporal, and energy-dissipation limits of spin-electronics. Yet, simultaneous control of high magnetic field, cryogenic environment, and nanometer resolution has remained an elusive benchmark for THz nanoscopy, leaving the local THz dynamics of these transitions largely unexplored. Here, we overcome these limitations by utilizing a custom-built cryogenic magneto-THz scattering-type scanning near-field optical microscopy (cm-THz-sSNOM) to resolve the near-field THz spectroscopic evolution of the magnetic field-driven CMR transition in a manganite single crystal. Our measurements provide a nanoscale visualization of the THz conductivity, capturing the moment that magnetic-field-induced spin switching triggers the transition from an antiferromagnetic insulator to a ferromagnetic metal. An ellipsoidal near-field model reveals a multi-scale transition initiated by 1–2 nm isolated spin-flip sites at low magnetic fields, which coalesce into ∼15 nm conducting regions as the threshold field is approached. These results provide an nano-THz view of CMR switching, establishing an analysis framework for mapping spin–charge–lattice–orbit–coupled dynamics at spatial scales that transcend the nominal sSNOM resolution.

Haeuser, Samuel [Iowa State University, Ames, IA (↗

DOE OSTI2020

Compositional phase stability of correlated electron materials within DFT+DMFT

Predicting the compositional phase stability of strongly correlated electron materials is an outstanding challenge in condensed matter physics, requiring precise computations of total energies. In this work, we employ the density functional theory plus dynamical mean-field theory (DFT+DMFT) formalism to address local correlations due to transition metal $\textit{d}$ electrons on compositional phase stability in the prototype rechargeable battery cathode material Li x CoO 2 , and detailed comparisons are made with the simpler DFT+$\textit{U}$ approach (i.e., the Hartree-Fock solution of the DMFT impurity problem). Local interactions are found to strongly impact the energetics of the band insulator LiCoO 2 , most significantly via the $E_g$ orbitals, which are partially occupied via hybridization with $\textit{O p}$ states. We find CoO 2 and Li 1/2 CoO 2 to be moderately correlated Fermi liquids with quasiparticle weights of 0.6–0.8 for the $T_{2g}$ states, which are most impacted by the interactions. Further, as compared to DFT+$\textit{U}$, DFT+DMFT considerably dampens the increase in total energy as $\textit{U}$ is increased, which indicates that dynamical correlations are important to describe this class of materials despite the relatively modest quasiparticle weights. Unlike DFT+$\textit{U}$, which can incorrectly drive Li x CoO 2 toward spurious phase separating or charge-ordered states, DFT+DMFT correctly captures the system's phase stability and does not exhibit a strong charge-ordering tendency. Most importantly, the error within DFT+$\textit{U}$ varies strongly as the composition changes, challenging the common practice of artificially tuning $\textit{U}$ within DFT+$\textit{U}$ to compensate the errors of Hartree-Fock. DFT+DMFT predicts the average intercalation voltage decreases relative to DFT, opposite to the result of DFT+$\textit{U}$, which would yield favorable agreement with experiment in conjunction with the overprediction of the voltage by the strongly constrained and appropriately normed DFT functional.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Nanosecond Structural Dynamics during Electrical Melting of Charge Density Waves in 1⁢T-TaS 2

Electrical control of charge density waves has been of immense interest, as the strong underlying electron-lattice interactions potentially open new, efficient pathways for manipulating their ordering and, consequently, their electronic properties. However, the transition mechanisms are often unclear as electric field, current, carrier injection, heat, and strain can all contribute and play varying roles across length scales and timescales. Here, we provide insight on how electrical stimulation melts the room temperature charge density wave order in 1T-TaS 2 by visualizing the atomic and mesoscopic structural dynamics from quasistatic to nanosecond pulsed melting. Using a newly developed ultrafast electron microscope setup with electrical stimulation, we reveal the order and strain dynamics during voltage pulses as short as 20 ns. The order parameter dynamics across a range of pulse amplitudes and durations support a thermally driven mechanism even for fields as high as 19 kV cm -1 . In addition, time-resolved imaging reveals a heterogeneous, mesoscopic strain response across the flake, including MHz-scale acoustic resonances that emerge during sufficiently short pulsed excitation which may modulate the order. In conclusion, these results suggest that metallic charge density wave phases like studied here may be more robust to electronic switching pathways than insulating ones, motivating further investigations at higher fields and currents in this and other related systems.

36 MATERIALS SCIENCE↗

DOE OSTI2020

A density-functional theory study of the Al/AlOx/Al tunnel junction

The aluminum oxide tunnel junction is a key component of the majority of superconducting quantum devices. For high-quality, reproducible, and scalably manufacturable qubits, the ability to fabricate Josephson junctions (JJs) with a targeted critical current and high uniformity is essential. In this study, we use first-principles modeling to assess fundamental aspects of the atomic structure of both amorphous and crystalline aluminum oxide tunnel junctions and relate the structure to predicted performance metrics. We use modified ab initio molecular dynamics to develop realistic models of the tunnel junction, from which interface roughness and local thickness fluctuations are analyzed in an unbiased manner by training a neural network to identify the boundary between metal and oxide. We show that the effective thickness of the insulating part of the junction can be different from the apparent physical thickness. We calculate the rate of Cooper pair tunneling for the atomically resolved electrostatic potential using direct numerical solution in 3D, which shows a channeling effect that impacts the junction critical current. The predicted critical current is a useful JJ design parameter that can be accessed from the ab initio calculations without fitting parameters. To assess the limits of uniformity and fabrication choices (e.g., oxidation vs epitaxy), we compare the amorphous junctions to crystalline models, which show order of magnitude more efficient tunneling compared to the amorphous case, underlining the connection between atomistic structure and Cooper pair tunneling efficiency. Further, this work provides a foundation for ab initio materials design and evaluation to help accelerate future development of improved tunnel junctions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗