Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 379 records · Page 21

Thermodynamics of Liquid Uranium from Atomistic and Ab Initio Modeling

We present thermodynamic properties for liquid uranium obtained from classical molecular dynamics (MD) simulations and the first-principles theory. The coexisting phases method incorporated within MD modeling defines the melting temperature of uranium in good agreement with the experiment. The calculated melting enthalpy is in agreement with the experimental range. Classical MD simulations show that ionic contribution to the total specific heat of uranium does not depend on temperature. The density of states at the Fermi level, which is a crucial parameter in the determination of the electronic contribution to the total specific heat of liquid uranium, is calculated by ab initio all electron density functional theory (DFT) formalism applied to the atomic configurations generated by classical MD. The calculated specific heat of liquid uranium is compared with the previously calculated specific heat of solid γ-uranium at high temperatures. The liquid uranium cannot be supercooled below T sc ≈ 800 K or approximately about 645 K below the calculated melting point, although, the self-diffusion coefficient approaches zero at T D ≈ 700 K. Uranium metal can be supercooled about 1.5 times more than it can be overheated. The features of the temperature hysteresis are discussed.

36 MATERIALS SCIENCE↗

A Panoramic View of MXenes via an Atomic Coordination‐Based Design Strategy

Two‐dimensional (2D) transition metal carbides and nitrides, known as MXenes, possess unique physical and chemical properties, enabling diverse applications in fields ranging from energy storage to communication, catalysis, sensing, healthcare, and beyond. Despite extensive research and notable advancements, a fundamental understanding of MXenes’ phase diversity and its connection to their hierarchical precursors, including the intermediate MAX phases and the ancestral bulk phases, remains limited. Here, in this study, it is hypothesized that the atomic coordination environments adopted by transition metal and nonmetallic atoms in their three‐dimensional (3D) bulk precursors may persist in 2D MXenes to govern their phase diversity. Using high‐throughput modeling based on first‐principles density functional theory, a wide range of MXene phases is unveiled and comprehensively evaluate their relative stabilities across a large chemical space. The key to the approach lies in considering various atomic coordination environments drawn from four types of ancestral bulk phases. Through this comprehensive structural library of MXenes, general guiding principles are uncovered, such as a close alignment between the phase stability of MXenes and that of their 3D precursors. These findings introduce a new design strategy in which the atomic coordination environments in bulk phases can serve as reliable predictors for accessing the diverse structural landscape of MXenes.

MXenes↗

Tuning the Magnetic Properties of CrI 3 Using Ni Thin Film Deposition for Applications in Spintronic Devices

Chromium triiodide (CrI 3 ), a van der Waals magnet, has recently been shown to host Ising ferromagnetism down to the monolayer limit. It is a potentially important material in 2D magnet-based applications, such as magnetic sensors and spintronic devices. Prior studies have revealed the coexistence of two different types of interlayer magnetic coupling, the antiferromagnetic (AFM) coupling near the surface and ferromagnetic (FM) coupling in the deep bulk layers, in pristine CrI 3 crystals below the Curie temperature. In this study, we used Ni thin film deposition to tune the surface magnetic states in bulk CrI 3 . A nanometer thickness (4 nm thick) of Ni was deposited on the surface of CrI 3 using electron-beam evaporation to form a Ni/CrI 3 heterostructure. The deposited Ni thin layer forms nanoclusters that completely cover the CrI 3 surface. Magnetic states of CrI 3 before and after Ni deposition are probed by ultralow-frequency magneto-Raman spectroscopy. Instead of seeing three spin wave branches in Raman scattering below 2 T, as in pristine bulk CrI 3 , in Ni/CrI 3 , we observe only a single spin wave branch, which is softened compared to that in pristine CrI 3 and displays a Zeeman shift under an out-of-plane magnetic field up to 7 T. This observation reveals that the AFM layers on pristine CrI 3 crystal surfaces transit into FM layers, so the entire CrI 3 crystal is in the FM state after Ni deposition. First-principles calculations show that Ni atoms tend to diffuse into the CrI 3 lattice, and the FM interlayer coupling has a much lower energy than AFM coupling in the presence of Ni atoms. As a result, our studies show that the magnetic state of CrI 3 can be modified through deposition of a thin metal layer on the surface, offering a route for controlling the spin degree of freedom in van der Waals magnets.

36 MATERIALS SCIENCE↗

Competing magnetic phases in Li-Fe-Ge kagome systems

Competing interlayer magnetic interactions in kagome magnets can lead to diverse magnetic phases, which enable promising topological or quantum material properties. In this work, the electronic structure and magnetic properties are studied using first-principles calculations for the Li⁢Fe6⁢G⁡e6, Li⁢Fe6⁢G⁡e4, and Li⁢Fe6⁢G⁡e5 compounds sharing the kagome F⁡e3⁢Ge layer motif but with different interlayer arrangements. For Li⁢Fe6⁢G⁡e4 and Li⁢Fe6⁢G⁡e5, the predicted magnetic ground states are collinear antiferromagnetic (AFM) states involving a mix of ferromagnetic (FM) and AFM interlayer orientations. In contrast, for Li⁢Fe6⁢G⁡e6, an incommensurate cycloidal spin spiral ground state is stabilized, although its energy is close to a collinear A-type AFM state. The analysis of magnetic RKKY exchange coupling confirms the results of calculations. The values of atomic magnetic moments are in good agreement with existing experimental estimations. Our experiments on Li⁢Fe6⁢G⁡e6 single crystals reveals AFM ordering at ∼540 K and transition to another magnetic phase with a small FM component, possibly with spin canting, below ∼270 K.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Nanoscale modulation of flat bands via controllable charge density wave defects in 4⁢𝐻⁡𝑏−Ta⁢S 2

Electron correlation is a main driver of exotic quantum phases and their interplay. The 4H b-TaS 2 system, possessing an intrinsic heterostructure of 1T- and 1H -TaS 2 monolayers, offers a unique opportunity to control electron correlation by distorting the atomic lattice or tuning interlayer coupling. Here, we investigated intrinsically deformed charge-density waves (CDWs) in the 1T layer of 4H b-TaS 2 to elucidate and control their effects on flat bands using scanning tunneling microscopy and spectroscopy (STM/S) combined with first-principles calculations. We identified two types of CDW defects: Type 1 has structural distortion and locally suppressed flat bands, while Type 2 features an increased flat band filling factor of intact CDW structure. Density functional theory calculations indicate that a sulfur vacancy in the 1T layer distorts the CDW structure and gives rise to a Type 1, whereas a sulfur vacancy in the 1H layer reduces the interlayer charge transfer and leads to a Type 2. Furthermore, we demonstrated creating and erasing individual CDW defects via STM manipulation. Here, our findings provide a pathway to not only tune flat bands but also selectively manipulate the interaction between CDW, the atomic lattice, and interlayer coupling in strongly correlated systems with atomic precision.

Charge density waves↗

Doping Effects on the Ductility of a Lightweight Refractory High-Entropy Alloy: Grain Boundary and Bulk Lattice Aspects

Doping elements in small amounts often segregate to grain boundaries (GBs) in alloys and can significantly impact mechanical properties and performance. Refractory high-entropy alloys (RHEAs) are known for their poor ductility, especially at low temperatures. Promoting GB cohesion through segregation can be an effective approach to mitigate embrittlement. Here, in this study, first-principles density functional theory (DFT) calculations were performed to examine the effects of important interstitial dopants (O, B, C, and N) and substitutional dopants (Cr, Y, La, and Ce) on the Σ5(310) [001] tilt GB of a lightweight RHEA Nb 32.5 Ti 27.5 Mo 22.5 Ta 12.5 Hf 2.5 Zr 2.5 . The DFT calculations reveal that certain dopants, such as B, C, Cr, N, and Y, exhibit favorable GB strengthening effects by improving bonding interactions with the bulk alloy. The impact of doping on the ductility parameter of the bulk lattice, defined as the ratio of surface energy to unstable stacking fault energy for the {110} <111> slip system, was also studied; and the results show that doping reduces the intrinsic ductility of the alloy, decreasing the D-parameter from 2.90 to 2.55, depending on the specific dopant. The present findings provide a foundational understanding at atomic level of the effect of representative dopants on mechanical properties of RHEAs and can be used to guide future alloy design for improved mechanical properties.

deformation charge density↗

Band Structure-Driven Design of a α-CsPbI 3 Ammonia Sensor for Industrial Applications

We investigate the defect-dependent electronic structure and gas-sensing potential of cubic α-CsPbI 3 using first-principles density functional theory and nonadiabatic molecular dynamics. Among the intrinsic defects, interstitials, vacancies, antisites, and switches studied, the I Pb and Pb I antisite defects exhibit transition energy levels near the middle of the band gap, thus functioning as deep traps. Short-term adsorption of ammonia selectively modifies the electronic structure, coordinating with Pb at Pb I sites and Cs at I Pb sites, significantly altering recombination pathways. Detailed analysis reveals that NH 3 reduces anharmonicity at I Pb defects, enabling enhanced recombination at elevated temperatures, while trap-assisted recombination dominates at room temperature. Other analytes, including CH 3 NH 2 and NO 2 , show negligible impact on the band gap or recombination dynamics, highlighting the potential selectivity of NH 3 interactions. Ab initio nonadiabatic molecular dynamics simulations at 300 K and 600 K further demonstrate temperature-dependent modulation of carrier lifetimes, with NH 3 accelerating recombination at ambient conditions and suppressing certain pathways at higher temperatures. These findings suggest that α-CsPbI 3 can serve as a selective and sensitive ammonia sensor over a broad temperature range and offer insights for ammonia detection under industrially relevant conditions.

ammonia sensor↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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.↗

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: ↗