SEARCH · Engineering Papers

Results for “First principles”

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 1,225 records · Page 68

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Unraveling the transformation pathway of the 𝛽 to 𝛾 phase transition in Ga 2 ⁢O 3 from atomistic simulations

Defect spinel 𝛾−Ga 2 ⁢O 3 is the least stable polymorph of Ga 2 ⁢O 3 , so its frequent appearance as a structural defect within or on the surface of monoclinic 𝛽−Ga 2 ⁢O 3 remains a mystery. Through first-principles calculations, we explore potential pathways for the phase transition from 𝛽−Ga 2⁢ O 3 to 𝛾−Ga 2 ⁢O 3 , and examine two key driving forces: tensile strain and Ga deficiency. When configurational entropy contributions to phase energies are included, the 𝛾 phase becomes energetically competitive with the 𝛽 phase, with the free energy difference between these phases diminishing even further under Ga-deficient conditions. Notably, a stability crossover occurs at room temperature at high vacancy concentrations ([V$^{3−}_{Ga}$]>3%) . A simple model 𝛽 → 𝛾 transformation pathway is identified, comprising two primary reactions, that enables the formation of the 𝛾 phase via simultaneous migration of Ga atoms from tetrahedral lattice sites to octahedral interstitial positions. The transformation barriers are prohibitively large in pristine Ga 2 ⁢O 3 , but can be substantially reduced by: (1) the presence of Ga vacancies, (2) elongational strains along the crystallographic 𝑎-axis, and (3) when volumetric relaxations are possible during transformation. These results elucidate prior experimental observations, where 𝛾−Ga 2⁢ O 3 is seen on damaged surfaces or in highly 𝑛-type 𝛽−Ga 2⁢ O 3 environments, which support Ga deficiency and mechanical strain. The insights into the driving forces and mechanisms of 𝛾−Ga 2⁢ O 3 formation enhance understanding of how localized strain and nonequilibrium defect concentrations may facilitate its formation from the 𝛽 phase.

Defects↗

DOE OSTI2022

Enhanced mobility of cations and anions in the redox state: The polaronium mechanism

Hugely enhanced slow-ion diffusivity has been widely observed under extreme redox conditions and for unclear reasons. Aided by first-principles calculations on model systems of ZrO 2 , CeO 2 , BaTiO 3 and Li 4/3 Mn 2/3 O 2 , here we successfully explained the intriguing phenomenon by a polaronium mechanism. We found a polaronium, defined as a transitory complex of a polaronic electron or hole and a migrating counterion, becomes highly mobile when the counterion comes from a d 0 or f 0 cation (e.g., Zr 4+ , Ti 4+ and Ce 4 ) or a p 6 anion (e.g., O 2- ) in the host compound. Upon a redox reaction, the complex attains a d 1 /f 1 or p 5 configuration, which spontaneously forms because it is favored by an electron-phonon interaction (manifest as the Jahn-Teller effect in high symmetry systems) that enables local relaxation and lowers the system energy. Here our calculations found such interaction reaching its peak at the saddle point where the local environment is softest, so soft that it allows a reorientation of the anisotropic d/f/p orbital to minimize the electron repulsion locally. Since the complex may dissolve after a successful ion-migration event, the redox electron/hole can be recycled to form another free-radical-like polaronium elsewhere, thereby enhancing ion migration repeatedly. The proposed polaronium mechanism, which also operates in ceramics doped with mixed-valence cations, is most relevant under dynamic and extreme thermal/field/irradiation conditions where extra electrons/holes are abundantly generated by non-equilibrium redox reactions. For such operations, some with emerging applications, our diffusion-enhancing mechanism may provide new theoretical insight to help understand their material/microstructure stability and performance.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Diffusion of acceptor dopants in monoclinic 𝛽−Ga 2⁢ O 3

𝛽−Ga 2 ⁢O 3 is a promising material for next-generation power electronics because of its ultrawide band gap and high critical breakdown voltage. However, realizing its full potential requires precise control over dopant incorporation and stability. In this work, we use first-principles calculations to systematically assess the diffusion behavior of eight potential deep-level substitutional acceptors (Au, Ca, Co, Cu, Fe, Mg, Mn, and Ni) in 𝛽−Ga 2 ⁢O 3 . We consider two key diffusion mechanisms: (i) interstitial diffusion under nonequilibrium conditions relevant to ion implantation, and (ii) trap-limited diffusion (TLD) under near-equilibrium thermal annealing conditions. Our results reveal a strong diffusion anisotropy along the 𝑏 and 𝑐 axes, with dopant behavior governed by competition between diffusion and incorporation (or dissociation) activation energies. Under interstitial diffusion, Ca$^{2+}_{i}$ and Mg$^{2+}_{i}$ show the most favorable combination of low migration and incorporation barriers, making them promising candidates for efficient doping along the 𝑏 and 𝑐 axes, respectively. In contrast, Au$^{+}_{i}$ diffuses readily, but exhibits an incorporation barrier that exceeds 5 eV, rendering it ineffective as a dopant. From a thermal stability perspective, Co$^{2+}_{i}$ shows poor activation but high diffusion barriers, which may suppress undesirable migration at elevated temperatures. Under trap-limited diffusion, the dissociation of dopant-host complexes controls mobility. Mg$^{2+}_{i}$ again emerges as a leading candidate, exhibiting the lowest dissociation barriers along both axes, whereas Co$^{2+}_{i}$ and Fe$^{2+}_{i}$ display the highest barriers, suggesting improved dopant retention under thermal stress. In conclusion, our findings guide dopant selection by balancing activation and thermal stability, essential for robust semi-insulating substrates.

Defects↗

DOE OSTI2025

Coverage-dependent structures and thermodynamic stability of intercalated Gd layers beneath buffer-layer graphene on SiC(0001)

Electronic properties of two-dimensional (2D) materials are strongly influenced by their atomic arrangements, making the theoretically-aided characterization of experimentally-synthesized 2D structures crucial. Using first-principles density functional theory, we analyze nearly 200 configurations of intercalated Gd layers beneath buffer-layer graphene on SiC(0001) over a Gd coverage range of 0.01 < θ < 1.2. By fully relaxing selectively-constructed configurations at each coverage within a large, low-strain supercell, we determine the coverage dependence of the chemical potential for intercalated Gd structures. Thermodynamically-preferred configurations below θ ≈ 0.8 form single-atom-thick monolayers, while 3D-like or multilayer structures emerge beyond θ ≈ 0.9. Most structures are amorphous-like, including the configuration at the chemical potential minimum around θ ≈ 0.4. In contrast, a strongly stretched Gd(0001)-like monolayer at θ = 1/3 and a nearly perfect Gd(0001) monolayer at θ = 1 are significantly less favorable with 0.16 eV and 0.82 eV higher chemical potentials above the minimum, respectively. Furthermore, the graphene layer decoupled by intercalated Gd near the chemical potential minimum is significantly flatter compared to its morphology above intercalated 3D structures at higher coverages and nearly isolated Gd atoms in the lowest coverage region. In conclusion, these findings align with our experimental results and underscore the need for further research on this unique intercalated system, which holds significant potential for diverse applications.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Unveiling the nature of Ga-based chalcogenides for electrical switching selectors

Three-dimensional phase-change memory with stackable crossbar architecture is a promising technology to meet the urgent demands for high-density storage and rapid information processing in the era of explosive data growth. The performance depends strongly on the properties of ovonic threshold switching (OTS) selectors, which control the on/off states of memory units. Amorphous GaS serves as an outstanding OTS material, distinguished by its sizable mobility gap and high crystallization temperature, while the underlying mechanism continues to be inadequately comprehended. Here, in this work, we systematically studied the structural and electronic properties of amorphous Ga-X (X = S/Se/Te) using first-principles calculations. The results show that Ga atoms adopt tetrahedral motifs, while S/Se/Te atoms predominantly exhibit the structure of a distorted triangular pyramid. This structural arrangement is ascribed to the substantial dative bonds formed by the lone-pair electrons of the anions and the vacant sp3 orbitals around Ga atoms. Large mobility gaps (e.g., GaS: 2.43 eV, GaSe: 1.76 eV, GaTe: 1.26 eV) and distinct mid-gap states (e.g., ∼0.66 eV above valence band tail) ensure that these three chalcogenide glasses can be switched on under an external electric field while effectively suppressing leakage current without a bias, and the defect electronic states originate from short, robust Ga-Ga bonds due to the formation of distorted chain-like local structures. Our research elucidates the mechanisms of amorphous Ga-X as OTS materials, enriching the spectrum of electrical switching selectors by incorporating III-VI chalcogenides. This inclusion offers novel opportunities for the refinement and optimization of high-density integrated memory systems.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Highly Anisotropic Quasi‐Direct Organic Metal Halide Hybrids: A Platform for Polarization‐Sensitive Optoelectronics

Low-dimensional organic–inorganic metal halide hybrids (OMHHs) exhibit remarkable optical properties and enhanced environmental stability. We investigate a 1D OMHH with formula C 4 N 2 H 14 PbBr 4 , consisting of Pb–Br chains separated by organic cations, which shows a large Stokes shift (0.83 eV) and broadband emission. Through first-principles calculations and polarized Raman spectroscopy, we characterize the material's vibrational properties and identify the specific phonon modes that drive exciton self-trapping. Our novel GW/Bethe-Salpeter equation force formalism reveals that low-frequency phonons (∼ 100 cm −1 , primarily involving Pb–Br motions) couple strongly with excitons, with a remarkably high Huang-Rhys factor of 137 ± 4, and gives a pathway for ultrafast structural analysis during the absorption process. This phonon-exciton coupling mechanism explains the material's broadband emission and provides a pathway for controlling optical properties through vibrations and for tuning vibrations through optical excitations. The material also exhibits highly anisotropic optical properties and electronic transport, with bands that are dispersive along the Pb–Br chains but nearly flat in perpendicular directions, resulting in direction-dependent electrical conductivity that is calculated to be an order of magnitude higher along the chain direction and consistent with measurements. These combined properties make this system an excellent platform for polarization-sensitive optoelectronic devices.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Lattice anharmonicity effects in fluorite oxide single crystals and anomalous increase in phonon lifetime in ceria at elevated temperature

Here, we investigate the temperature dependence of the frequency and linewidth of the triply degenerate T 2g zone-centered optical phonon in flux-grown ceria and hydrothermally synthesized thoria single crystals from room temperature to 1273 K using Raman spectroscopy. Both crystals exhibit an expected increase in the phonon linewidth with temperature due to enhanced phonon–phonon scattering. However, ceria displays an anomalous linewidth reduction in the temperature range of 1023–1123 K. First-principles phonon linewidth calculations considering cubic and quartic phonon interactions within temperature-independent phonon dispersion fail to describe this anomaly. A parameterization of the temperature-dependent second-order interatomic force constants based on previously reported phonon dispersion measured at room and high temperatures predicts a deviation from the monotonic linewidth increase, albeit at temperatures lower than those observed experimentally for ceria. The qualitative agreement in the trend of temperature-dependent linewidth suggests that lattice anharmonicity-induced phonon renormalization plays a role in phonon lifetime. Specifically, a change in the overlap between softened acoustic and optical branches in the dispersion curve reduces the available phonon scattering phase space of the Raman-active mode at the zone center, leading to an increased phonon lifetime within a narrow temperature interval. These findings provide insights into higher-order anharmonic interactions in ceria and thoria, motivating further investigations into the role of anharmonicity-induced phonon renormalization on phonon lifetimes at high temperatures.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Engineering 2D Square Lattice Hubbard Models in 90° Twisted GeX/SnX (X =S, Se) Moiré Superlattices

Because of the large-period superlattices emerging in moiré two-dimensional (2D) materials, electronic states in such systems exhibit low energy flat bands that can be used to simulate strongly correlated physics in a highly tunable setup. While many investigations have thus far focused on moiré flat bands and emergent correlated electron physics in triangular, honeycomb, and quasi-one-dimensional lattices, tunable moiré realizations of square lattices subject to strong correlations remain elusive. Here, in this work, we propose a feasible scheme to construct moiré square lattice systems by twisting two or more layers of 2D materials in a rectangular lattice by 90°. We demonstrate the concept with twisted GeX/SnX (X =S, Se) moiré superlattices and calculate their electronic structures from first principles. We show that the lowest conduction flat band in these systems can be described by a square lattice Hubbard model with parameters which can be controlled by varying the choice of host materials, number of layers, and external electric fields. In particular, twisted double bilayer GeSe realizes a square lattice Hubbard model with strong frustration due to the next-nearest-neighbor hopping that could host unconventional superconductivity, in close analogy to the Hubbard model for copper-oxygen planes of cuprate high-temperature superconductors. The presented scheme uses 90° twisted 2D materials with rectangular unit cells as a promising platform for realizing the physical phenomena of square lattice Hubbard models, establishing a new route for studying its rich phase diagram of magnetism, charge order, and unconventional superconductivity in a highly tunable setting.

2-dimensional systems↗

DOE OSTI2025

Atomic Structure, Dynamics, Changes in Chemical Bonding and Semiconductor-Metal Transition in Sb 2 Se 3 : A Remarkable Material for Quantum Networks and Energy Applications

Antimony sesquiselenide has become an outstanding functional material for photovoltaics, energy storage and transformation, memory and photonic applications. Sb 2 Se 3 is one of the most successful emerging solar light absorbers and has also been identified as a highly promising ultralow-loss phase-change material (PCM) for next-generation coherent nanophotonic processors, photonic tensor cores, quantum and neuromorphic networks. Unlike benchmark telluride PCMs, Sb 2 Se 3 features a quasi-one-dimensional (1D) crystalline structure consisting of (Sb 4 Se 6 ) ∞ ribbons, lacks the typical PCM chemical bonding, and undergoes an extended semiconductor-metal transition above the melting point. Consequently, the origin of high optical contrast between crystalline (SET) and amorphous (RESET) logic states remains elusive and presents a significant challenge. Using high-energy X-ray diffraction and Raman spectroscopy over a wide temperature range, supported by first-principles simulations and complemented by thermal, optical and electrical measurements, as well as by 121 Sb-Mossbauer spectroscopy, the quasi-1D network of orthorhombic antimony sesquiselenide was found to undergo significant evolution in amorphous and supercooled Sb 2 Se 3 , leading to lower coordination, shorter interatomic distances and a higher p-electron density on antimony, indicating changes in chemical bonding. The observed novel Sb 2 Se 3 nanocrystalline polymorph, characterized by trigonal antimony coordination and more isolated Sb-Se ribbons, could help reduce multiple trapping defect states in the bandgap, which are typical of orthorhombic Sb 2 Se 3 , thereby enhancing the power-conversion efficiency of photovoltaic devices. Semimetallic and metallic liquid Sb 2 Se 3 exhibit a gradual transformation into a denser 2D and/or 3D network with higher antimony coordination. Localized electron states in the pseudogap are becoming extended, leading to an increase in electronic conductivity σ following the relationship σ ∝ N(E F ) 2 . Liquid Sb 2 Se 3 also appears to be strongly fragile, with a nonmonotonic change in viscosity and higher atomic mobility in the metallic liquid. Furthermore, these results explain extraordinary functionalities of Sb 2 Se 3 for photonic and energy applications.

antimony↗

DOE OSTI2026

Development of an ab initio learned model of electron deposition range in deuterium-tritium plasmas through time-dependent density functional theory calculations and machine learning

Accurate hydrodynamic modeling for laser-direct-drive (LDD) inertial-confinement-fusion (ICF) relies on precise calculations of the electron thermal conduction in all target materials. The nonlocal stopping range of electrons in ICF plasmas directly influences thermal conduction; yet, no first principles model exists for the electron mean free path in the conduction-zone regime. This work utilized time-dependent stochastic density-functional theory (TD-sDFT) to calculate the electron stopping power in deuterium-tritium (DT) plasmas at (ρ, T) conditions relevant to the conduction zone and the compressed shell in ICF. Using a combination of our TD-sDFT data and already established analytical models, we developed and trained an artificial neural network to create a global model for the nonlocal electron deposition range, λ E . We compared our machine-learning (ML) based model for λ E to the currently-used modified-Lee-More model in LDD radiation-hydrodynamic codes, such as lilac, and saw an overall decrease in the deposition range. To understand the effects of λ E on LDD ICF implosion dynamics, we implemented the ML-based model into lilac; specifically, we looked at designs consistent with a current experiment on the OMEGA laser and for a newly designed LDD-ICF target for the future OMEGA-Next facility. In both cases, we saw an overall drop in predicted ablation pressure, peak areal density, and neutron yield due to the reduced thermal conduction (smaller λ E ) in DT plasmas. Comparisons with the experiment on OMEGA are also made.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Rise and fall of Mott insulating gaps in YNiO 3 paramagnets as a reflection of symmetry breaking and remaking

The YNiO 3 nickelate is a paradigm d-electron oxide that manifests the intriguing temperature-mediated sequence of three phases transitions from (i) magnetically ordered insulator to (ii) paramagnetic (PM) insulator and then to (iii) PM metal. Such phenomena raised the question of the nature of the association of magnetism and structural symmetry breaking with the appearance in (i) and (ii) and disappearance in (iii) of insulating band gaps. It is demonstrated here that first-principles mean-field–like density-functional theory (DFT), driven by molecular dynamics temperature evolution, can describe not only the origin of the magnetically long-range ordered insulating phase (i), but also the creation of an insulating paramagnet (ii) that lacks spin- long-range order, and of a metallic paramagnet (iii) as temperature rises. This approach provides the patterns of structural and magnetic symmetry breaking at different temperatures, in parallel with band gaps obtained when the evolving geometries are used as input to DFT electronic band-structure calculations. This disentangles the complex interplay among spin, charge, and orbital degrees of freedom. Analysis shows that the success in describing the rise and fall of the insulating band gaps along the phase transition sequence is enabled by allowing sufficient flexibility in describing diverse local structural and magnetic motifs as input to DFT. This entails the use of sufficiently large supercells that allow expressing structural disproportionation of octahedra, as well as a description of PM phases as a distribution of local magnetic moments (rather than using a single averaged moment). It appears that the historic dismissal of mean-field–like DFT as being unable to describe such Mott-like transitions was premature, as it was based on consideration of averaged crystallographic unit cells, a description that washes out local symmetry-breaking motifs. The magnetically ordered insulating YNiO 3 phase (i) and the PM insulating phase (ii) result in DFT from allowing symmetry breaking, evident already by considering the athermal internal energy. In contrast, the PM metallic phase (iii) is formed thermally by smearing out thus weakening symmetry breaking. Analysis of snapshots of the different forms of structural vs magnetic symmetry breaking shows that only the loss of the polymorphous distribution of magnetic moments existing in (ii) causes the fall of the band gap, resulting in the metallic state in (iii). Finally, the interesting conclusion is that such a description of the rise [in phases (i) and (ii)] and fall [in phase (iii)] of the insulating gap does not rely on the traditional Mott-like strong correlation understanding, but on breaking and remaking of magnetic and structural symmetries reflected in energy lowering.

36 MATERIALS SCIENCE↗

NASA NTRS1985

Modeling of pressure-induced far-infrared absorption spectra Molecular hydrogen pairs

Meyer et al. (1985) have calculated the accurate induced dipole moment function of H2-H2 from first principles, using highly correlated wave functions for the first time in such work. The present paper is concerned with the collision-induced translational-rotational absorption coefficient for molecular hydrogen pairs, taking into account computations on the basis of the fundamental theory considered by Meyer et al. Data have been obtained for temperatures in the range from 40 to 300 K. Criteria are developed for choosing among various model line shapes. It is found that certain models are capable of approximating the quantum profiles closely, with rms errors of only a few percent.

Borysow, J.↗

DOE OSTI2020

An attempt to simulate laser-induced all-optical spin switching in a crystalline ferrimagnet

Interest in all-optical spin switching (AOS) is growing rapidly. The recent discovery of AOS in Mn 2 RuGa provides a much needed clean case of crystalline ferrimagnets for theoretical simulations. Here, we attempt to simulate it using the state-of-the-art first-principles method combined with the Heisenberg exchange model. We first compute the spin moments at two inequivalent manganese sites and then feed them into our model Hamiltonian. We employ an ultrafast laser pulse to switch the spins. We find that there is a similar optimal laser field amplitude to switch spins. However, we find that the exchange interaction has a significant effect on the system switchability. Weakening the exchange interaction could make the system unswitchable. This provides a crucial insight into the switching mechanism in ferrimagnets.

36 MATERIALS SCIENCE↗

NASA NTRS2018

ALMA Spectral Line and Imaging Survey of a Low and a High Mass-Loss Rate AGB Star Between 335 and 362 GHz

Context. Low and intermediate mass stars are known to power strong stellar winds when evolving through the asymptotic giant branch (AGB) phase. Initial mass, luminosity, temperature, and composition determine the pulsation characteristics of the star and the dust species formed in the pulsating photospheric layers. Radiation pressure on these grains triggers the onset of a stellar wind. However, as of today, we still cannot predict the wind mass-loss rates and wind velocities from first principles neither do we know which species are the first to condense in the upper atmospheric regions. Aims. We aim to characterise the dominant physical, dynamical, and chemical processes in the inner wind region of two archetypical oxygen-rich (C/O<1) AGB stars, that is, the low mass-loss rate AGB star R Dor (M~ 1x10(exp 7) M☉/yr) and the high mass-loss rate AGB star IK Tau (M~5x10(exp 6) M☉/yr). The purpose of this study is to observe the key molecular species contributing to the formation of dust grains and to cross-link the observed line brightnesses of several species to the global and local properties of the star and its wind. Methods. A spectral line and imaging survey of IK Tau and R Dor was made with ALMA between 335 and 362 GHz (band 7) at a spatial resolution of ~150 mas, which corresponds to the locus of the main dust formation region of both targets. Results. Some two hundred spectral features from 15 molecules (and their isotopologues) were observed, including rotational lines in both the ground and vibrationally excited states (up to v=5 for SiO). Detected species include the gaseous precursors of dust grains such as SiO, AlO, AlOH, TiO, and TiO2.We present a spectral atlas for both stars and the parameters of all detected spectral features. A clear dichotomy for the sulphur chemistry is seen: while CS, SiS, SO, and SO2 are abundantly present in IK Tau, only SO and SO2 are detected in R Dor. Also other species such as NaCl, NS, AlO, and AlOH display a completely different behaviour. From some selected species, the minor isotopologues can be used to assess the isotopic ratios. The channel maps of many species prove that both large and small-scale inhomogeneities persist in the inner wind of both stars in the form of blobs, arcs, and/or a disk. The high sensitivity of ALMA allows us to spot the impact of these correlated density structures in the spectral line profiles. The spectral lines often display a half width at zero intensity much larger than expected from the terminal velocity, v1, previously derived for both objects (36 km/s versus v(∞)~17.7 km/s for IK Tau and 23 km/s versus v(∞)~5.5 km/s for R Dor). Both a more complex 3D morphology and a more forceful wind acceleration of the (underlying) isotropic wind can explain this trend. The formation of fractal grains in the region beyond _400 mas can potentially account for the latter scenario. From the continuum map, we deduce a dust mass of ~3.7 x 10(exp 7) M☉ and ~2 x 10(ex0 8) M☉ for IK Tau and R Dor, respectively. Conclusions. The observations presented here provide important constraints on the properties of these two oxygen-dominated AGB stellar winds. In particular, the ALMA data prove that both the dynamical and chemical properties are vastly different for this high mass-loss rate (IK Tau) and low mass-loss rate (R Dor) star.

Stars: AGB and post-AGB, Stars: mass loss, Stars: ↗

NASA NTRS2018

ALMA Spectral Line and Imaging Survey of a Low and a High Mass-Loss Rate AGB Star Between 335 and 362 GHz

Context. Low and intermediate mass stars are known to power strong stellar winds when evolving through the asymptotic giant branch (AGB) phase. Initial mass, luminosity, temperature, and composition determine the pulsation characteristics of the star and the dust species formed in the pulsating photospheric layers. Radiation pressure on these grains triggers the onset of a stellar wind. However, as of today, we still cannot predict the wind mass-loss rates and wind velocities from first principles neither do we know which species are the first to condense in the upper atmospheric regions. Aims. We aim to characterise the dominant physical, dynamical, and chemical processes in the inner wind region of two archetypical oxygen-rich (C/O < 1) AGB stars, that is, the low mass-loss rate AGB star R Dor (M’ ∼ 1 × 10(exp −7) M/yr) and the high mass-loss rate AGB star IK Tau (M’ ∼ 5 × 10(exp −6) M/yr). The purpose of this study is to observe the key molecular species contributing to the formation of dust grains and to cross-link the observed line brightnesses of several species to the global and local properties of the star and its wind. Methods. A spectral line and imaging survey of IK Tau and R Dor was made with ALMA between 335 and 362 GHz (band 7) at a spatial resolution of ∼150 mas, which corresponds to the locus of the main dust formation region of both targets. Results. Some two hundred spectral features from 15 molecules (and their isotopologues) were observed, including rotational lines in both the ground and vibrationally excited states (up to v = 5 for SiO). Detected species include the gaseous precursors of dust grains such as SiO, AlO, AlOH, TiO, and TiO2. We present a spectral atlas for both stars and the parameters of all detected spectral features. A clear dichotomy for the sulphur chemistry is seen: while CS, SiS, SO, and SO2 are abundantly present in IK Tau, only SO and SO2 are detected in R Dor. Also other species such as NaCl, NS, AlO, and AlOH display a completely different behaviour. From some selected species, the minor isotopologues can be used to assess the isotopic ratios. The channel maps of many species prove that both large and small-scale inhomogeneities persist in the inner wind of both stars in the form of blobs, arcs, and/or a disk. The high sensitivity of ALMA allows us to spot the impact of these correlated density structures in the spectral line profiles. The spectral lines often display a half width at zero intensity much larger than expected from the terminal velocity, v∞, previously derived for both objects (36 km/s versus v(∞) ∼ 17.7 km/s for IK Tau and 23 km/s versus v(∞) ∼ 5.5 km/s for R Dor). Both a more complex 3D morphology and a more forceful wind acceleration of the (underlying) isotropic wind can explain this trend. The formation of fractal grains in the region beyond ∼400 mas can potentially account for the latter scenario. From the continuum map, we deduce a dust mass of ∼3.7 × 10(exp −7) Mand ∼2 × 10(exp −8) M for IK Tau and R Dor, respectively. Conclusions. The observations presented here provide important constraints on the properties of these two oxygen-dominated AGB stellar winds. In particular, the ALMA data prove that both the dynamical and chemical properties are vastly different for this high mass-loss rate (IK Tau) and low mass-loss rate (R Dor) star.

L. Decin↗

DOE OSTI2020

Multi-channel direct detection of light dark matter: theoretical framework

We present a unified theoretical framework for computing spin-independent direct detection rates via various channels relevant for sub-GeV dark matter — nuclear re- coils, electron transitions and single phonon excitations. Despite the very different physics involved, in each case the rate factorizes into the particle-level matrix element squared, and an integral over a target material- and channel-specific dynamic structure factor. We show how the dynamic structure factor can be derived in all three cases following the same procedure, and extend previous results in the literature in several aspects. For electron transitions, we incorporate directional dependence and point out anisotropic target materials with strong daily modulation in the scattering rate. For single phonon excitations, we present a new derivation of the rate formula from first principles for generic spin-independent couplings, and include the first calculation of phonon excitation through electron couplings. We also discuss the interplay between single phonon excitations and nuclear recoils, and clarify the role of Umklapp processes, which can dominate the single phonon production rate for dark matter heavier than an MeV. Our results highlight the complementarity between various search channels in probing different kinematic regimes of dark matter scattering, and provide a common reference to connect dark matter theories with ongoing and future direct detection experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Unveiling the strong interaction among hadrons at the LHC

One of the key challenges for nuclear physics today is to understand from first principles the effective interaction between hadrons with different quark content. First successes have been achieved using techniques that solve the dynamics of quarks and gluons on discrete space-time lattices. Experimentally, the dynamics of the strong interaction have been studied by scattering hadrons off each other. Such scattering experiments are difficult or impossible for unstable hadrons and so high-quality measurements exist only for hadrons containing up and down quarks. Here we demonstrate that measuring correlations in the momentum space between hadron pairs produced in ultrarelativistic proton–proton collisions at the CERN Large Hadron Collider (LHC) provides a precise method with which to obtain the missing information on the interaction dynamics between any pair of unstable hadrons. Specifically, we discuss the case of the interaction of baryons containing strange quarks (hyperons). We demonstrate how, using precision measurements of proton–omega baryon correlations, the effect of the strong interaction for this hadron–hadron pair can be studied with precision similar to, and compared with, predictions from lattice calculations. The large number of hyperons identified in proton–proton collisions at the LHC, together with accurate modelling of the small (approximately one femtometre) inter-particle distance and exact predictions for the correlation functions, enables a detailed determination of the short-range part of the nucleon-hyperon interaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1976

Study of the extra-ionic electron distributions in semi-metallic structures by nuclear quadrupole resonance techniques

A straightforward self-consistent method was developed to estimate solid state electrostatic potentials, fields and field gradients in ionic solids. The method is a direct practical application of basic electrostatics to solid state and also helps in the understanding of the principles of crystal structure. The necessary mathematical equations, derived from first principles, were presented and the systematic computational procedure developed to arrive at the solid state electrostatic field gradients values was given.

Murty, A. N.↗