Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Methods in superconductivity”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

On-the-fly machine learned force fields for the study of warm dense matter: Application to diffusion and viscosity of CH

We develop a framework for on-the-fly machine learned force field (MLFF) molecular dynamics (MD) simulations of warm dense matter (WDM). In particular, we employ an MLFF scheme based on the kernel method and Bayesian linear regression, with the training data generated from the Kohn–Sham density functional theory (DFT) using the Gauss spectral quadrature method, within which we calculate energies, atomic forces, and stresses. We verify the accuracy of the formalism by comparing the predicted properties of warm dense carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that ab initio MD simulations of WDM can be accelerated by up to three orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the diffusion coefficients and shear viscosity of CH at a density of 1 g/cm3 and temperatures in the range of 75 000–750 000 K. We find that the self- and inter-diffusion coefficients and the viscosity obey a power law with temperature, and that the diffusion coefficient results suggest a weak coupling between C and H in CH. In addition, we find agreement within standard deviation with previous results for C and CH but disagreement for H, demonstrating the need for ab initio calculations as presented here.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Techniques for Polytypic Transformations in Silicon Carbide

Abstract—: Two main polytype transformations in silicon carbide, namely, 2H → 6H and 3C → 6H, have been studied by ab initio methods. It has been shown that the intermediate phases with trigonal symmetry P3m1 and monoclinic symmetry Cm make it much easier to move the close-packed layers in such transitions by breaking them up into separate stages. It has been found that these two polytype transformations proceed in completely different ways. The links being relocated noticeably tilted compared to their initial position at the transition 2H → 6H, which allows the compression of the SiC links in the plane (11 2-bar 0). The transition 3C → 6H is carried out through the formation of Si–Si and C–C auxiliary links, living for a short time and helping densely packed layers to swap places. As a result, the activation barrier of the transformation 2H → 6H (1.7 eV/atom) is significantly less than the activation barrier of the transformation 3C → 6H (3.6 eV/atom), which means that the second transition should occur at the temperatures by 750–800°C higher than the first one. The energy profiles of this polytypic transformations, as well as the geometry of all intermediate and transition phases have been calculated. It has been shown that all the transition states have monoclinic symmetry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Cumulant methods for electron-phonon problems. I. Perturbative expansions

Here in this work, we investigate the ability of the cumulant expansion (CE) to capture one-particle spectral information in electron-phonon coupled systems at both zero and finite temperatures. In particular, we present a comprehensive study of the second- and fourth-order CEs for the one-dimensional Holstein model as compared with numerically exact methods. We investigate both finite sized systems as well as the approach to the thermodynamic limit, drawing distinctions, and connections between the behavior of systems in and away from the thermodynamic limit that enable a greater understanding of the ability of the CE to capture real-frequency information across the full range of wave vectors. We find that for zero electronic momentum, the spectral function is well described by the second-order CE at low and high temperatures. However, for nonzero electronic momenta, the CE is only accurate at high temperature. We analyze the fourth-order cumulant and find that while it improves the description of the short-time dynamics encoded in the one-particle Green's function, it can introduce divergences in the time domain as well as unphysical negative spectral weight in the spectral function. When well-behaved, the fourth-order CE does provide notable accurate corrections to the second-order CE. Finally, we use our results to comment on the use of the CE as a tool for calculating transport behavior in the realistic ab initio modeling of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Electronic Structure, Lattice Dynamics, and Superconducting Properties of Mercury-Alkaline Earth Metal Compounds: a First-Principles Study

The electronic structure, lattice dynamics, and superconducting properties of SrHg and BaHg intermetallic compounds are investigated using plane wave pseudo-potential and linear response methods based on density functional theory. The ground state properties such as lattice parameters are in good agreement with previously reported results. The positive phonon frequencies reveal the stability of both compounds in B{sub 2} (CsCl-type) structure. The metallic nature of SrHg and BaHg is confirmed by the electronic band structure as well as phonon dispersion curves. The superconducting behavior described in terms of Eliashberg spectral function (α{sup 2}F(ω)) and electron-phonon coupling constant (λ). The calculated values of superconducting transition temperatures (T{sub c}) show good agreement with available experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

The cosmic DANCe of Perseus

Context. Star-forming regions are excellent benchmarks for testing and validating theories of star formation and stellar evolution. The Perseus star-forming region, being one of the youngest (< 10 Myr), closest (280-320 pc), and most studied in the literature, is a fundamental benchmark. Aims. We aim to study the membership, phase-space structure, mass, and energy (kinetic plus potential) distribution of the Perseus star-forming region using public catalogues (Gaia, APOGEE, 2MASS, and Pan-STARRS). Methods. We used Bayesian methodologies that account for extinction to identify the Perseus physical groups in the phase-space, retrieve their candidate members, derive their properties (age, mass, 3D positions, 3D velocities, and energy), and attempt to reconstruct their origin. Results. We identify 1052 candidate members in seven physical groups (one of them new) with ages between 3 and 10 Myr, dynamical super-virial states, and large fractions of energetically unbounded stars. Their mass distributions are broadly compatible with that of Chabrier for masses ≳0.1 M ⊙ and do not show hints of over-abundance of low-mass stars in NGC 1333 with respect to IC 348. These groups’ ages, spatial structure, and kinematics are compatible with at least three generations of stars. Future work is still needed to clarify if the formation of the youngest was triggered by the oldest. Conclusions. The exquisite Gaia data complemented with public archives and mined with comprehensive Bayesian methodologies allow us to identify 31% more members than previous studies, discover a new physical group (Gorgophone: 7 Myr, 191 members, and 145 M ⊙ ), and confirm that the spatial, kinematic, and energy distributions of these groups support the hierarchical star formation scenario.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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↗
DOE OSTI2020

A Virtual Element Method for the Full MHD system

In this manuscript we present a novel discretization for the incompressible MHD system. Our approach follows the framework of the Virtual Element Method and offers two main advantages. The method can be implemented in unstructured meshes making it highly versatile and capable of handling a wide array of problems involving interfaces, free-boundaries or adaptive refinements on the mesh. The second advantage involves the divergence of the magnetic field, our approach guarantees that it remains solenoidal. We include a theoretical proof of the condition on the magnetic field as well as energy estimates and a well-posedness study. The latter sheds light as to the stability properties of the method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Phase diagram of twisted bilayer graphene at filling factor $v = ±3$

Here, we study the correlated insulating phases of twisted bilayer graphene (TBG) in the absence of lattice strain at integer filling $v = ±3$. Using the self-consistent Hartree-Fock method on a particle-hole symmetric model and allowing translation symmetry breaking terms, we obtain the phase diagram with respect to the ratio of AA interlayer hopping (w 0 ) and AB interlayer hopping (w 1 ). When the interlayer hopping ratio is close to the chiral limit (w 0 /w 1 ≲ 0.5), a quantum anomalous Hall state with Chern number v c = ±1 can be observed consistent with previous studies. Around the realistic value w 0 /w 1 ≈ 0.8, we find a spin and valley polarized, translation symmetry breaking, state with C 2⁢z T symmetry, a charge gap and a doubling of the moiré unit cell, dubbed theC 2⁢z T stripe phase. The real-space total charge distribution of this C 2⁢z T stripe phase in the flat band limit does not have modulation between different moiré unit cells, although the charge density in each layer is modulated, and the translation symmetry is strongly broken. Other symmetries, including C 2⁢z , C 2⁢x , particle-hole symmetry P, and the topology of the C 2⁢z T stripe phase, are also discussed in detail. We observed braiding and annihilation of the Dirac nodes by continuously turning on the order parameter to its fully self-consistent value, and provide a detailed explanation of the mechanism for the charge gap opening despite preserving C 2⁢z T and valley U⁡(1) symmetries. In the transition region between the quantum anomalous Hall phase and the C 2⁢z T stripe phase, we find an additional competing state with comparable energy corresponding to a phase with a tripling of the moiré unit cell.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Real time time-dependent density functional theory using higher order finite-element methods

We present a computationally efficient approach to solve the time-dependent Kohn-Sham equations in real time using higher order finite-element spatial discretization, applicable to both pseudopotential and all-electron calculations. To this end, we develop an a priori mesh adaption technique, based on the semidiscrete (discrete in space but continuous in time) error estimate on the time-dependent Kohn-Sham orbitals, to construct an efficient finite-element discretization. Subsequently, we obtain the full-discrete error estimate to guide our choice of the time step. We employ spectral finite elements along with special reduced order quadrature to render the overlap matrix diagonal, thereby simplifying the inversion of the overlap matrix that features in the evaluation of the discrete time-evolution operator. We use the second-order Magnus operator as the time-evolution operator, wherein the action of the discrete Magnus operator, expressed as exponential of a matrix, on the Kohn-Sham orbitals is obtained efficiently through an adaptive Lanczos iteration. We observe close to optimal rates of convergence of the dipole moment with respect to spatial and temporal discretization, for both pseudopotential and all-electron calculations. We demonstrate a staggering 100-fold reduction in the computational time afforded by higher order finite elements over linear finite elements, for both pseudopotential and all-electron calculations. Further, for similar level of accuracy, we obtain significant computational savings by our approach as compared to state-of-the-art finite-difference methods. Here, we also demonstrate the competence of higher order finite elements for all-electron benchmark systems. Lastly, we observe good parallel scalability of the proposed method on many hundreds of processors.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Characterizing the Distribution of Microwave Surface Resistance of HTS Film Based on Metal Ring

This paper presents a novel method for measuring the microwave surface resistance distribution of high-temperature superconducting thin films by a dielectric resonator with a metal ring. After calibration, the effect of metal ring loss on measurement accuracy is negligible. The measurement resolution is determined by the internal diameter of the metal ring, which makes the design of the device more flexible. Therefore, measurement with high resolution and convenience can be achieved simultaneously by using a metal ring with proper size and operating at a relatively low frequency. A test device working at 32 GHz with TE{sub 01(1 + δ)} mode was developed. Its resolution is a circular area with a diameter of 5 mm and an area of 19.6 mm{sup 2}. The unloaded quality factor of the test device is more than 140,000, when it is loaded by the calibration cavity. The test device has a relative uncertainty of less than 6%, which complies with the uncertainty (20%) specified by the international standard for the microwave surface resistance measurement of high-temperature superconducting (HTS) thin film. The surface resistance distribution of a 2-in. YBCO/LAO/YBCO superconducting film sample was measured by this test device.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Superstructures and Superconductivity Linked with Pd Intercalation in Nb 2 Pd x Se 5

Single crystals of Nb 2 Pd x Se 5 have been grown covering a broad range of palladium stoichiometries, from x = 0.68 to x = 0.95 (Pd2 occupancies 0.36 and 0.90 respectively). Miscibility gaps in the Pd stoichiometry range are related to superstructures resulting from staged filling and ordering of the partially filled Pd atom site (Pd2, Wyckoff site 2a). For Pd2 near 1/2 occupancy, a superstructure based on an incommensurate occupation wave with modulation along the monoclinic b-axis is identified. Here, a superconducting dome is found for Pd2 occupancies above 0.58, with maximum Tc for Pd2 occupancy of approximately 2/3. DFT calculations for the fully and 1/2 occupied Pd2 site show that Nb 2 Pd x Se 5 can be considered an intercalation compound for Pd2 occupancy above the minimum required for structural stability.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Emergence of low-energy spin waves in superconducting electron-doped cuprates

In order to fully utilize the technological potential of unconventional superconductors, an enhanced understanding of the superconducting mechanism is necessary. In the best performing superconductors, the cuprates, superconductivity is intimately linked with magnetism, although the details of this coupling remain elusive. Here, we address this gap by studying the electron-doped cuprate Nd 1.85 Ce 0.15 CuO 4−δ that has an antiferromagnetic ground state when synthesized and only becomes superconducting after a reductive annealing process. Using neutron spectroscopy, we show that the as-grown crystal exhibits a large spin pseudogap in the magnetic fluctuation spectrum. Annealing removes defects introduced by the commonly employed synthesis method and significantly reduces the spin pseudogap. While the spin pseudogap in the annealed sample likely arises from superconductivity, in the as-grown sample it results from the absence of long-wavelength spin waves. These results reveal a direct connection between defects, magnetism, and superconductivity, offering new insight into the mechanisms underlying high-temperature superconductivity and guiding the design of improved superconducting materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Thermal effects on the electronic properties of sodium electride under high pressures

Sodium is a simple metal at ambient conditions, while it transits to an electride phase at pressures above ~160 GPa along the room temperature isotherm. Here we explore the thermal effects on the electronic properties of the $hP4$ phase of sodium along the $ρ$ = 5.872 g/cm 3 isochore. We quantitatively classify this phase as an insulator based on the criterion of nearsightedness of the one-particle density matrix. Ab initio calculations suggest that the band gap of this insulator decreases with increasing temperature along the isochore, primarily because of ionic distortions; culminating in an insulator-to-metal transition upon melting at $T_m \approx 2100$ $\hbox{K}$. This transition is accompanied by residual electronic localization (in d orbitals) in the form of dynamic electron bubbles and a change in hybridization form $p–d$ to $s–p$ upon melting. This transition is explored by tracking the electronic and electro-optical properties along the isochore under consideration.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Lattice collective modes from a continuum model of magic-angle twisted bilayer graphene

Here, we show that the insulating states of magic-angle twisted bilayer graphene support a series of collective modes corresponding to local particle-hole excitations on triangular lattice sites. Our theory is based on a continuum model of the magic angle flat bands. When the system is insulating at moiré band filling ν = -3, our calculations show that the ground state supports seven low-energy modes that lie well below the charge gap throughout the moiré Brillouin zone, one of which couples strongly to THz photons. The low-energy collective modes are faithfully described by a model with a local SU(8) degree of freedom in each moiré unit cell that we identify as the direct product of spin, valley, and an orbital pseudospin. Apart from spin and valley-wave modes, the collective mode spectrum includes a low-energy intraflavor exciton mode associated with transitions between flat valence and conduction band orbitals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Dynamical vertex approximation for many-electron systems with spontaneously broken SU(2) symmetry

We generalize the formalism of the dynamical vertex approximation (DΓA)—a diagrammatic extension of the dynamical mean-field theory (DMFT)—to treat magnetically ordered phases. To this aim, we start by concisely illustrating the many-electron formalism for performing ladder resummations of Feynman diagrams in systems with broken SU(2) symmetry associated to ferromagnetic (FM) or antiferromagnetic (AF) order. We then analyze the algorithmic simplifications introduced by taking the local approximation of the two-particle irreducible vertex functions in the Bethe-Salpeter equations, which defines the ladder implementation of DΓA for magnetic systems. The relation of this assumption with the DMFT limit of large coordination-number/high dimensions is explicitly discussed. As a last step, we derive the expression for the ladder DΓA self-energy in the FM- and AF-ordered phases of the Hubbard model. The physics emerging in the AF-ordered case is explicitly illustrated by means of approximated calculations based on a static mean-field input for DΓA equations. The results obtained capture fundamental aspects of both metallic and insulating ground states of two-dimensional antiferromagnets, providing a reliable compass for future, more extensive applications of our approach. Furthermore, possible routes to further develop diagrammatic-based treatments of magnetic phases in correlated electron systems are briefly outlined in the Conclusions.

2-dimensional systems↗
DOE OSTI2021

Competing magnetic states in transition metal dichalcogenide moiré materials

Small-twist-angle transition metal dichalcogenide (TMD) heterobilayers develop isolated flat moiré bands that are approximately described by triangular lattice generalized Hubbard models. In this paper we explore the metallic and insulating states that appear under different control conditions at a density of one electron per moiré period and the transitions between them. By combining fully self-consistent Hartree-Fock theory calculations with strong-coupling expansions around the atomic limit, we identify four different magnetic states and one nonmagnetic state near the model phase diagram’s metal-insulator phase transition line. Here, ferromagnetic insulating states, stabilized by nonlocal direct exchange interactions, are surprisingly prominent.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Fluctuations analysis of spin susceptibility: Néel ordering revisited in dynamical mean field theory

Here in this paper, we revisit the antiferromagnetic (AF) phase diagram of the single-band three-dimensional half-filled Hubbard model on a simple cubic lattice studied within the dynamical mean field theory (DMFT). Although this problem has been investigated extensively in the literature, a comprehensive understanding of the impact of the different one- and, in particular, two-particle local correlation functions of DMFT on the AF transition temperature is still missing. We have, hence, performed a fluctuation analysis of T N with respect to different local bosonic fluctuations (charge, spin, particle-particle) contained in the two-particle vertex of DMFT. Our results indicate that, beyond weak coupling, the screening of the DMFT vertex by local fluctuations leads to an enhancement of T N with respect to a random phase approximation (RPA) like calculation where this vertex is replaced by the bare interaction. The overall suppression of T N in DMFT with respect to RPA is then solely due to the incoherence introduced by the DMFT self-energy in the one-particle Green's functions. This illustrates the Janus-faced role of the local moment formation in the DMFT solution of the Hubbard model, which leads to completely opposite effects in the one- and two-particle correlation functions.

36 MATERIALS SCIENCE↗