Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DENSITY OF STATES”

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 235 records · Page 13

Analytical derivatives of the individual state energies in ensemble density functional theory. II. Implementation on graphical processing units (GPUs)

Conical intersections control excited state reactivity, and thus, elucidating and predicting their geometric and energetic characteristics are crucial for understanding photochemistry. Locating these intersections requires accurate and efficient electronic structure methods. Unfortunately, the most accurate methods (e.g., multireference perturbation theories such as XMS-CASPT2) are computationally challenging for large molecules. The state-interaction state-averaged restricted ensemble referenced Kohn–Sham (SI-SA-REKS) method is a computationally efficient alternative. The application of SI-SA-REKS to photochemistry was previously hampered by a lack of analytical nuclear gradients and nonadiabatic coupling matrix elements. We have recently derived analytical energy derivatives for the SI-SA-REKS method and implemented the method effectively on graphical processing units. We demonstrate that our implementation gives the correct conical intersection topography and energetics for several examples. Furthermore, our implementation of SI-SA-REKS is computationally efficient, with observed sub-quadratic scaling as a function of molecular size. This demonstrates the promise of SI-SA-REKS for excited state dynamics of large molecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Topological and magnetic properties of a noncollinear spin state on a honeycomb lattice in a magnetic field

Here, this paper studies the topological and magnetic properties of a noncollinear spin state on a honeycomb lattice that evolves from coplanar to ferromagnetic with a magnetic field applied along the z axis. The coplanar state is stabilized by nearest-neighbor ferromagnetic interactions, single-ion anisotropy along z, and DzyaloshinskiiMoriya interactions between next-nearest-neighbor sites. Below the critical field H$_c$ that aligns the spins, the magnetic unit cell contains six sites and the spin dynamics contains six magnon modes. Although the classical energy is degenerate with respect to the twist angle φ between nearest-neighbor spins, the dependence of the free energy on φ at low temperatures is dominated by the magnon zero-point energy, which contains extremum at φ = πl/3 for integer l. The only unique ground states GS(φ) have l = 0 or 1. For H < H$'_c$, the zero-point energy has minima at even l and the ground state is GS(0); for H$'_c$ < H < H$'_c$, the zero-point energy has minima at odd l and the ground state is GS(π/3). In GS(0), the magnon density of states exhibits five distinct topological phases with increasing field associated with the opening and closing of energy gaps between two or three magnonic bands. While the Berry curvature vanishes for the coplanar φ = 0 phase in zero field, the Berry curvature and Chern numbers exhibit signatures of the five topological phases below H$'_c$. Whereas the Berry curvature and Chern number are sensitive to changes in the magnon density of states within GS(π/3), the inelastic spectrum S(k,ω) is sensitive to changes in the intensity of the magnon modes in the different magnetic phases GS(0) and GS(π/3) rather than the five topological phases within GS(π/3).

Fishman, Randy S. [Oak Ridge National Laboratory (↗

Echo state network for coarsening dynamics of charge density waves

An echo state network (ESN) is a type of reservoir computer that uses a recurrent neural network with a sparsely connected hidden layer. Compared with other recurrent neural networks, one great advantage of ESN is the simplicity of its training process. Yet, despite the seemingly restricted learnable parameters, ESN has been shown to successfully capture the spatial-temporal dynamics of complex patterns. Here we build an ESN to model the coarsening dynamics of charge-density waves (CDWs) in a semiclassical Holstein model, which exhibits a checkerboard electron density modulation at half-filling stabilized by a commensurate lattice distortion. The inputs to the ESN are local CDW order parameters in a finite neighborhood centered around a given site, while the output is the predicted CDW order of the center site at the next time step. Special care is taken in the design of couplings between hidden layer and input nodes to ensure lattice symmetries are properly incorporated into the ESN model. Since the model predictions depend only on CDW configurations of a finite domain, the ESN is scalable and transferrable in the sense that a model trained on dataset from a small system can be directly applied to dynamical simulations on larger lattices. Furthermore, our work opens avenues for efficient dynamical modeling of pattern formations in functional electron materials.

2-dimensional systems↗

Multiband superconductivity in strongly hybridized 1 T ' - WTe 2 / NbSe 2 heterostructures

The interplay of topology and superconductivity has become a subject of intense research in condensed-matter physics for the pursuit of topologically nontrivial forms of superconducting pairing. An intrinsically normal-conducting material can inherit superconductivity via electrical contact to a parent superconductor via the proximity effect, usually understood as Andreev reflection at the interface between the distinct electronic structures of two separate conductors. However, at high interface transparency, strong coupling inevitably leads to changes in the band structure, locally, owing to hybridization of electronic states. Here, we investigate such strongly proximity-coupled heterostructures of monolayer $1T$' - $\mathrm{WTe_2}$ grown on NbSe 2 by van der Waals epitaxy. The superconducting local density of states, resolved in scanning tunneling spectroscopy down to 500 mK, reflects a hybrid electronic structure well described by a multiband framework based on the McMillan equations which captures the multiband superconductivity inherent to the NbSe 2 substrate and that is induced by proximity to $\mathrm{WTe_2}$ self-consistently. Our material-specific tight-binding model captures the hybridized heterostructure quantitatively and confirms that strong interlayer hopping gives rise to a semimetallic density of states in the two-dimensional $\mathrm{WTe_2}$ bulk, even for nominally band-insulating crystals. The model further accurately predicts the measured order parameter Δ ≅ 0.6 meV induced in the $\mathrm{WTe_2}$ monolayer bulk, stable beyond a 2 T magnetic field. So we believe that our detailed multiband analysis of the hybrid electronic structure provides a useful tool for sensitive spatial mapping of induced order parameters in proximitized atomically thin topological materials.

36 MATERIALS SCIENCE↗

Neutron thermalization in nuclear graphite: A modern story of a classic moderator

Graphite has been of great interest since the dawn of the nuclear era. In this work, we will focus on the thermal neutron scattering of two types of nuclear graphite that have different microstructures, namely, the historic PGA which has a very anisotropic medium-grained needle coke, and the new G347A which is a near-isotropic fine-grained graphite. We report the results of room temperature inelastic neutron scattering measurements for the scattering function, double differential scattering cross-section, and phonon density of states. The obtained results are discussed in conjunction with results extracted from first-principles lattice dynamics simulations of crystalline graphite, using the incoherent approximation and by also including the coherent-one phonon contribution. Excellent agreements have been found between the measured and calculated scattering functions and phonon densities of states. The inclusion of the coherent-one phonon scattering shows a significant improvement of the calculated scattering functions, double differentials, and inelastic scattering cross-sections. We also highlight the physical concerns about the current thermal neutron scattering ENDF/B-VIII of nuclear graphite library and its corresponding phonon density of states. This work aims to provide a new presentation of neutron thermalization by graphite and offers a new set of high quality experimental and computational data for consideration.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Strain engineering of doped hydrogen passivated silicon quantum dots

Silicon quantum dots are nanomaterials that are attractive candidates for photovoltaic applications. Doping of these materials creates p-n junctions and is important for solar cells. In this work, we present a first-principles study of the coupled influence of doping and strain on the stability, energy gap, Fermi level, electronic density, and density of states of hydrogen-passivated silicon quantum dots. We find that the cohesive energy and the energy gap decrease with increasing quantum dot size and are strongly influenced by strain. Furthermore, the response to strain also depends on the size of the quantum dot and dopant type. We present expressions of cohesive energy and energy gap as power-law of size and polynomial dependence on strain. We also show that the Fermi energy increases with size for pristine and p-type doping but decreases with size for n-type doping. We also discuss the influence of strain and dopant type on the density of states and electron density of the quantum dots.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Thermal Neutron Scattering Law for Beryllium Hydride and Critical Mass Calculations

The thermal neutron scattering law (TSL) for crystalline beryllium hydride (BeH 2 ) is developed from first-principles ab initio lattice dynamics calculations and the impact of neutron thermalization in this material on critical mass is estimated. BeH 2 has a body-centered orthorhombic crystal structure with 12 molecules per unit cell and a theoretical density of 0.755 g/cm 3 . The vibrational (phonon) densities of states for H and Be bound in BeH 2 are determined using VASP density functional theory and PHONON lattice dynamics calculations. The TSLs for H bound in BeH 2 , H(BeH 2 ), and Be bound in BeH 2 , Be(BeH 2 ), are then evaluated in the incoherent approximation from the calculated H and Be partial phonon density of states using FLASSH. Finally, critical mass as a function of 235 U loading density for bare and reflected BeH 2 moderated spheres is predicted from MC21 Monte Carlo neutron transport calculations using ENDF/B-VIII.0 cross sections and the H(BeH 2 ) and Be(BeH 2 ) TSL evaluations. Comparisons are made to water (H 2 O), polyethylene (CH 2 ), and beryllium oxide (BeO) as moderators. These critical mass predictions are a refinement upon the prior work by Rao and Srinivasan that neglected thermal neutron scattering effects. The minimum critical mass of a BeH 2 moderated assembly is estimated to be 0.207 kg 235 U for a 0.20 m thick BeO reflected sphere and 0.178 kg 235 U for a 0.40 m thick BeO reflected sphere.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Marcus–Hush–Chidsey kinetics at electrode–electrolyte interfaces

Electrochemical kinetics at electrode–electrolyte interfaces limit the performance of devices including fuel cells and batteries. While the importance of moving beyond Butler–Volmer kinetics and incorporating the effect of electronic density of states of the electrode has been recognized, a unified framework that incorporates these aspects directly into electrochemical performance models is still lacking. In this work, we explicitly account for the density functional theory-calculated density of states numerically in calculating electrochemical reaction rates for a variety of electrode–electrolyte interfaces. We first show the utility of this for two cases related to Li metal electrodeposition and stripping on a Li surface and a Cu surface (anode-free configuration). The deviation in reaction rates is minor for cases with flat densities of states such as Li, but is significant for Cu due to nondispersive d-bands creating large variation. Finally, we consider a semiconducting case of a solid-electrolyte interphase consisting of LiF and Li2CO3 and note the importance of the Fermi level at the interface pinned by the redox reaction occurring there. We identify the asymmetry in reaction rates as a function of discharge/charge naturally within this approach.

Kurchin, Rachel (ORCID:0000000221474809)↗

Spatial BCS-BEC crossover in superconducting $\textit{p–n}$ junctions

Here, we present a theory of superconducting $\textit{p–n}$ junctions. To this end, we consider a two band model of doped bulk semiconductors with attractive interactions between the charge carriers and derive the superconducting order parameter, the quasiparticle density of states, and the chemical potential as a function of the semiconductor gap $Δ_0$ and the doping level ε. We verify previous results for the quantum phase diagram for a system with constant density of states in the conduction and valence band, which show BCS-superconductor to Bose-Einstein-condensation (BEC) and BEC-to-insulator transitions as a function of doping level and the size of the band gap. Then, we extend this formalism to a density of states which is more realistic for 3D systems and derive the corresponding quantum phase diagram, where we find that a BEC phase can only exist for small band gaps $Δ_0 < Δ ^{\ast}_0$. For larger band gaps, we find rather a direct transition from an insulator to a BCS phase. Next, we apply this theory to study the properties of superconducting $\textit{p–n}$ junctions. We derive the spatial variation of the superconducting order parameter along the $\textit{p–n}$ junction. As the potential difference across the junction leads to energy band bending, we find a spatial crossover between a BCS and BEC condensate, as the density of charge carriers changes across the $\textit{p–n}$ junction. For the two-dimensional system, we find two possible regimes, when the bulk is in a BCS phase, a BCS-BEC-BCS junction with a single BEC layer in the space charge region, and a BCS-BEC-I-BEC-BCS junction with two layers of BEC condensates separated by an insulating layer. In three dimensions we find that there can also be a conventional BCS-I-BCS junction for semiconductors with band gaps exceeding $Δ^{\ast}_0$. Thus, we find that there can be BEC layers in the well controlled setting of doped semiconductors, where the doping level can be varied to change and control the thickness of BEC and insulator layers, making Bose-Einstein condensates thereby possibly accessible to experimental transport and optical studies in solid-state materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Microscopic calculations of nuclear level densities with the Lanczos method

A new method for computing the density of states in nuclei making use of an extrapolated form of the tridiagonal matrix obtained from the Lanczos method is presented. It will be shown that the global, average properties of the entire Lanczos matrix can be predicted from just four Lanczos iterations. The extrapolated Lanczos matrix (ELM) approach provides for an accurate computation of the density of states described within the configuration space, which, in some cases, is sufficient to accurately calculate the density of states at, or near, the neutron separation energy. Comparisons between theory and experiment are shown for $^{57}Fe$, $^{74}Ge$, and $^6Ge$. In addition, we show results for the J -dependence of moments and the level density for these three nuclei.

Shell model↗

Absence of bulk charge density wave order in the normal state of UTe 2

A spatially modulated superconducting state, known as pair density wave (PDW), is a tantalizing state of matter with unique properties. Recent scanning tunneling microscopy (STM) studies revealed that spin-triplet superconductor UTe 2 hosts an unprecedented spin-triplet, multi-component PDW whose three wavevectors are indistinguishable from a preceding charge-density wave (CDW) order that survives to temperatures well above the superconducting critical temperature, T c . Whether the PDW is the mother or a subordinate order remains unsettled. Here, based on a systematic search for bulk charge order above T c using resonant elastic X-ray scattering (REXS), we show that the structure factor of charge order previously identified by STM is absent in the bulk within the sensitivity of REXS. Our results invite two scenarios: either the density-wave orders condense simultaneously at T c in the bulk, in which case PDW order is likely the mother phase, or the charge modulations are restricted to the surface.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Scaling of pulsed nanosecond capillary plasmas at different specific energy deposition

Nano-second, capillary discharges (nCDs) are unique plasma sources in their ability to sustain high specific energy deposition ω dep approaching 10 eV/molecule in molecular gases. This high energy loading on short timescales produces both high plasma densities and high densities of molecular exited states. These high densities of electrons and excited states interact with each other during the early afterglow through electron collision quenching and associative ionization. In this paper we discuss results from a two-dimensional computational investigation of a nCD sustained in air at a pressure of 28.5 mbar and with a voltage amplitude 20 kV. Discharges were investigated for two circuit configurations—a floating low voltage electrode and with the low voltage electrode connected to ground through a ballast resistor. The first configuration produced a single ionization wave from the high to low voltage electrode. The second produced converging ionization waves beginning at both electrodes. With a decrease of the tube radius, the velocity of the ionization fronts decreased while the shape of the ionization wave changed from the electron density being distributed smoothly in the radial direction, to being hollow shaped where there is a higher electron density near the tube wall. For sufficiently small tubes, the near-wall maxima merge to have the higher density on the axis of the capillary tube. In the early afterglow, the temporal and radial behavior of the N 2 (C 3 Π u ) density is a sensitive function of ω dep due to electron collision quenching. Finally, these trends indicate that starting from ω dep $\geqslant$ 0.3 eV/molecule, it is necessary to take into account interactions of electrons with electronically excited species during the discharge and early afterglow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Atomistic nature of amorphous graphite

This paper focuses on the structural, electronic, and vibrational features of amorphous graphite [R. Thapa et al, Phys. Rev. Lett., 2022, 128, 236402]. The structural order in amorphous graphite is discussed and compared with graphite and amorphous carbon. The electronic density of states and localization in these phases were analyzed. Spatial projection of charge densities in the π bands showed a high charge concentration on participating atoms in connecting hexagons. A vibrational density of states was computed and is potentially an experimentally testable fingerprint of the material. An analysis of the vibrational modes was carried out using the phase quotient, and the mode stretching character. The average thermal conductivity calculated for aG was 0·85 and 0·96 W/cmK at room temperature and 1000 K, respectively.

Chemistry↗

Atomic-Site-Specific Surface Valence-Band Structure from X-Ray Standing-Wave Excited Photoemission

X-ray standing-wave (XSW) excited photoelectron emission was used to measure the site-specific valence band (VB) for ½ monolayer (ML) Pt grown on a SrTiO 3 (001) surface. The XSW induced modulations in the core level (CL), and VB photoemission from the surface and substrate atoms were monitored for three ℎ⁡kl substrate Bragg reflections. The XSW CL analysis shows the Pt to have a face-centered-cubic-like cube-on-cube epitaxy with the substrate. The XSW VB information compares well to a density functional theory calculated projected density of states from the surface and substrate atoms. Altogether, this Letter represents a novel method for determining the contribution to the density of states by valence electrons from specific atomic surface sites.

74 ATOMIC AND MOLECULAR PHYSICS↗

Chiral condensate and the equation of state at nonzero baryon density from the hadron resonance gas model with a repulsive mean field

We study the QCD equation of state and the chiral condensate using the hadron resonance gas model with repulsive mean-field interactions. We find that the repulsive interactions improve the agreement with the lattice results on the derivatives of the pressure with respect to the baryon chemical potential up to eighth order. From the temperature dependence of the chiral condensate we estimate the crossover temperature as a function of baryon chemical potential, T p c ( μ B ) . We find that the chiral crossover line starts to deviate significantly from the chemical freeze-out line already for μ B > 400 MeV . Furthermore, we find that the chiral pseudocritical line can be parametrized as T p c ( μ B ) / T p c ( 0 ) = 1 − κ 2 [ μ B / T p c ( 0 ) ] 2 − κ 4 [ μ B / T p c ( 0 ) ] 4 with κ 2 = 0.0150 ( 2 ) and κ 4 = 3.1 ( 6 ) × 10 − 5 , which are in agreement with lattice QCD results for small values of μ B . For the first time we find a tiny but nonzero value of κ 4 in our study. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Phase-transition-induced thermal hysteresis in Type-II weyl semimetals M o T e 2 and M o 1 − x W x T e 2

The resistivity versus temperature measurement is commonly used for identifying temperature-induced phase-change and the resulting hysteresis loop. However, while the resistance is influenced by both the density of states and the carrier lifetimes, the Seebeck coefficient is influenced predominantly by the density of states, and hence is a better probe of the phase of the material. Here, 1T'-T d temperature-induced phase transition in MoTe 2 is studied using temperature-dependent X-ray diffraction, resistivity, and Seebeck coefficient measurements. A more distinct hysteresis is observed when measuring the Seebeck coefficient which is consistent with direct measurements of the crystallographic angle using the temperature-dependent X-ray diffraction. The Seebeck and electrical resistivity measurements indicate a competing contribution of the electrons and holes. The contribution of electron pockets becomes more dominant when molybdenum atoms are replaced by tungsten. In MoTe 2 , a topologically induced enhancement of the Nernst coefficient is observed at low temperatures and a relatively large phase-transition induced Thomson coefficient of 111 μV∙K -1 is measured at 254 K.

36 MATERIALS SCIENCE↗

Electron scattering with ethane adsorbed on rare gas multilayers: Hole transfer, coulomb decay, and ion dissociation

Positive ion desorption following electron impact dissociative ionization of ethane adsorbed on Ar, Kr, and Xe multilayers has been studied as a function of incident electron energy from threshold to 100 eV. Based on the dependence of ion yields on the identity of the rare gas, it is likely that the majority of ethane molecules undergo indirect ionization following hole transfer from the ionized underlying rare gas. Further, this has also been corroborated by density of states calculations showing the energetic alignment of the outer valence states of ethane and the condensed rare gas ionization energies. Due to the near-resonant nature of charge transfer for single-hole states, the ethane molecular ion is excited to different final ionic states on different rare gases, which leads to differences in ion desorption yields and branching ratios. The quantitative yields increase with increasing ionization energy gap between the rare gas and ethane, in the order Ar > Kr >Xe. The large increase in yields from 25 eV onwards for all rare gases is likely due to the formation and decay of two-hole states on neighboring rare gas and ethane molecules due to interatomic and intermolecular Coulomb decay (ICD) and not electron transfer mediated decay (ETMD). The ICD and ETMD pathways become accessible when the incoming electron has sufficient energy to excite the inner valence ns level of the rare gas to a Rydberg state or ionize it. The experimental findings are supported by calculations of thresholds, density of states for the final configurations of these processes, and coupling strengths for hole transfer between ethane and rare gases. The fragment ion branching ratios vary with energy from threshold to about 35 eV, showing the fragmentation pattern changes with the mode of hole transfer and availability of excess energy. Sigma C–C bonds are more likely to break than C–H bonds in the mid-20 eV range, and this effect is most pronounced for Xe, followed by Kr, and then Ar.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Techniques for studying materials under extreme states of high energy density compression

The properties of materials under extreme conditions of pressure and density are of key interest to a number of fields, including planetary geophysics, materials science, and inertial confinement fusion. In geophysics, the equations of state of planetary materials, such as hydrogen and iron, under ultrahigh pressure and density provide a better understanding of their formation and interior structure [Celliers et al., “Insulator-metal transition in dense fluid deuterium,” Science 361, 677–682 (2018) and Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)]. The processes of interest in these fields occur under conditions of high pressure (100 GPa–100 TPa), high temperature (>3000 K), and sometimes at high strain rates (>103 s−1) depending on the process. With the advent of high energy density (HED) facilities, such as the National Ignition Facility (NIF), Linear Coherent Light Source, Omega Laser Facility, and Z, these conditions are reachable and numerous experimental platforms have been developed. To measure compression under ultrahigh pressure, stepped targets are ramp-compressed and the sound velocity, measured by the velocity interferometer system for any reflector diagnostic technique, from which the stress-density of relevant materials is deduced at pulsed power [M. D. Knudson and M. P. Desjarlais, “High-precision shock wave measurements of deuterium: Evaluation of exchange-correlation functionals at the molecular-to-atomic transition,” Phys. Rev. Lett. 118, 035501 (2017)] and laser [Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)] facilities. To measure strength under high pressure and strain rates, experimenters measure the growth of Rayleigh–Taylor instabilities using face-on radiography [Park et al., “Grain-size-independent plastic flow at ultrahigh pressures and strain rates,” Phys. Rev. Lett. 114, 065502 (2015)]. The crystal structure of materials under high compression is measured by dynamic x-ray diffraction [Rygg et al., “X-ray diffraction at the national ignition facility,” Rev. Sci. Instrum. 91, 043902 (2020) and McBride et al., “Phase transition lowering in dynamically compressed silicon,” Nat. Phys. 15, 89–94 (2019)]. Medium range material temperatures (a few thousand degrees) can be measured by extended x-ray absorption fine structure techniques, Yaakobi et al., “Extended x-ray absorption fine structure measurements of laser-shocked V and Ti and crystal phase transformation in Ti,” Phys. Rev. Lett. 92, 095504 (2004) and Ping et al., “Solid iron compressed up to 560 GPa,” Phys. Rev. Lett. 111, 065501 (2013), whereas more extreme temperatures are measured using x-ray Thomson scattering or pyrometry. This manuscript will review the scientific motivations, experimental techniques, and the regimes that can be probed for the study of materials under extreme HED conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗