No Arabic abstract
Advances in merged-beams instruments have allowed experimental studies of the mutual neutralisation (MN) processes in collisions of both Li$^+$ and Na$^+$ ions with D$^-$ at energies below 1 eV. These experimental results place constraints on theoretical predictions of MN processes of Li$^+$ and Na$^+$ with H$^-$, important for non-LTE modelling of Li and Na spectra in late-type stars. We compare experimental results with calculations for methods typically used to calculate MN processes, namely the full quantum (FQ) approach, and asymptotic model approaches based on the linear combination of atomic orbitals (LCAO) and semi-empirical (SE) methods for deriving couplings. It is found that FQ calculations compare best overall with the experiments, followed by the LCAO, and the SE approaches. The experimental results together with the theoretical calculations, allow us to investigate the effects on modelled spectra and derived abundances and their uncertainties arising from uncertainties in the MN rates. Numerical experiments in a large grid of 1D model atmospheres, and a smaller set of 3D models, indicate that neglect of MN can lead to abundance errors of up to 0.1 dex (26%) for Li at low metallicity, and 0.2 dex (58%) for Na at high metallicity, while the uncertainties in the relevant MN rates as constrained by experiments correspond to uncertainties in abundances of much less than 0.01~dex (2%). This agreement for simple atoms gives confidence in the FQ, LCAO and SE model approaches to be able to predict MN with the accuracy required for non-LTE modelling in stellar atmospheres.
Deuterated molecules are important chemical tracers of prestellar and protostellar cores. Up to now, the titular reaction has been assumed to contribute to the generation of these deuterated molecules. We have measured the merged-beams rate coefficient for this reaction as function of the relative collision energy in the range of about 10 meV to 10 eV. By varying the internal temperature of the reacting H$_3^+$ molecules, we found indications for the existence of a reaction barrier. We have performed detailed theoretical calculations for the zero-point-corrected energy profile of the reaction and determined a new value for the barrier height of $approx$ 68 meV. Furthermore, we have calculated the tunneling probability through the barrier. Our experimental and theoretical results show that the reaction is essentially closed at astrochemically relevant temperatures. We derive a thermal rate coefficient of $<1times 10^{-12}$ cm$^3$ s$^{-1}$ for temperatures below 75 K with tunneling effects included and below 155 K without tunneling.
The ratio of deuterium to hydrogen (D/H ratio) of Solar System bodies is an important clue to their formation histories. Here we fit a Neptunian atmospheric model to Gemini Near Infrared Spectrograph (GNIRS) high spectral resolution observations and determine the D/H ratio in methane absorption in the infrared H-band ($sim$ 1.6 {mu}m). The model was derived using our radiative transfer software VSTAR (Versatile Software for the Transfer of Atmospheric Radiation) and atmospheric fitting software ATMOF (ATMOspheric Fitting). The methane line list used for this work has only become available in the last few years, enabling a refinement of earlier estimates. We identify a bright region on the planetary disc and find it to correspond to an optically thick lower cloud. Our preliminary determination of CH$_{rm 3}$D/CH$_{rm 4}$ is 3.0$times10^{-4}$, which is in line with the recent determination of Irwin et al. (2014) of 3.0$^{+1.0}_{-0.9}simtimes10^{-4}$, made using the same model parameters and line list but different observational data. This supports evidence that the proto-solar ice D/H ratio of Neptune is much less than that of the comets, and suggests Neptune formed inside its present orbit.
The influence of the uncertainties in the rate coefficient data for electron-impact excitation and ionization on non-LTE Li line formation in cool stellar atmospheres is investigated. We examine the electron collision data used in previous non-LTE calculations and compare them to recent calculations that use convergent close-coupling (CCC) techniques and to our own calculations using the R-matrix with pseudostates (RMPS) method. We find excellent agreement between rate coefficients from the CCC and RMPS calculations, and reasonable agreement between these data and the semi-empirical data used in non-LTE calculations up to now. The results of non-LTE calculations using the old and new data sets are compared and only small differences found: about 0.01 dex (~ 2%) or less in the abundance corrections. We therefore conclude that the influence on non-LTE calculations of uncertainties in the electron collision data is negligible. Indeed, together with the collision data for the charge exchange process Li(3s) + H <-> Li^+ + H^- now available, and barring the existence of an unknown important collisional process, the collisional data in general is not a source of significant uncertainty in non-LTE Li line formation calculations.
Cosmological models can be constrained by determining primordial abundances. Accurate predictions of the He I spectrum are needed to determine the primordial helium abundance to a precision of $< 1$% in order to constrain Big Bang Nucleosynthesis models. Theoretical line emissivities at least this accurate are needed if this precision is to be achieved. In the first paper of this series, which focused on H I, we showed that differences in $l$-changing collisional rate coefficients predicted by three different theories can translate into 10% changes in predictions for H I spectra. Here we consider the more complicated case of He atoms, where low-$l$ subshells are not energy degenerate. A criterion for deciding when the energy separation between $l$ subshells is small enough to apply energy-degenerate collisional theories is given. Moreover, for certain conditions, the Bethe approximation originally proposed by Pengelly & Seaton (1964) is not sufficiently accurate. We introduce a simple modification of this theory which leads to rate coefficients which agree well with those obtained from pure quantal calculations using the approach of Vrinceanu et al. (2012). We show that the $l$-changing rate coefficients from the different theoretical approaches lead to differences of $sim 10$% in He I emissivities in simulations of H II regions using spectral code Cloudy.
We performed 3-D simulations of proton-rich material entrainment into czw-rich He-shell flash convection and the subsequent H-ingestion flash that took place in the post-AGB star Sakurais object. Observations of the transient nature and anomalous abundance features are available to validate our method and assumptions, with the aim to apply them to very low metallicity stars in the future. We include nuclear energy feedback from H burning and cover the full $4pi$ geometry of the shell. Runs on $768^3$ and $1536^3$ grids agree well with each other and have been followed for $1500mathrm{min}$ and $1200mathrm{min}$. After a $850mathrm{min}$ long quiescent entrainment phase the simulations enter into a global non-spherical oscillation that is launched and sustained by individual ignition events of H-rich fluid pockets. Fast circumferential flows collide at the antipode and cause the formation and localized ignition of the next H-overabundant pocket. The cycle repeats for more than a dozen times while its amplitude decreases. During the global oscillation the entrainment rate increases temporarily by a factor $approx 100$. Entrained entropy quenches convective motions in the upper layer until the burning of entrained H establishes a separate convection zone. The lower-resolution run hints at the possibility that another global oscillation, perhaps even more violent will follow. The location of the H-burning convection zone agrees with a 1-D model in which the mixing efficiency is calibrated to reproduce the light curve. The simulations have been performed at the NSF Blue Waters supercomputer at NCSA.