We present dynamical studies of the dissociation of polycyclic aromatic hydrocarbon (PAH) radical cations in their ground electronic states with significant internal energy. Molecular dynamics ...simulations are performed, the electronic structure being described on-the-fly at the self-consistent-charge density functional-based tight binding (SCC-DFTB) level of theory. The SCC-DFTB approach is first benchmarked against DFT results. Extensive simulations are achieved for naphthalene , pyrene and coronene at several energies. Such studies enable one to derive significant trends on branching ratios, kinetics, structures and hints on the formation mechanism of the ejected neutral fragments. In particular, dependence of branching ratios on PAH size and energy were retrieved. The losses of H and C2H2 (recognized as the ethyne molecule) were identified as major dissociation channels. The H/C2H2 ratio was found to increase with PAH size and to decrease with energy. For , which is the most interesting PAH from the astrophysical point of view, the loss of H was found as the quasi-only channel for an internal energy of 30 eV. Overall, in line with experimental trends, decreasing the internal energy or increasing the PAH size will favour the hydrogen loss channels with respect to carbonaceous fragments.
This article is part of the themed issue ‘Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces’.
ABSTRACT
We present an analysis of photometric, spectroscopic, and spectropolarimetric data of the nearby, cool, magnetic DZ white dwarf PM J08186−3110. High-dispersion spectra show the presence of ...Zeeman splitted spectral lines due to the presence of a surface average magnetic field of 92 kG. The strong magnesium and calcium lines show extended wings shaped by interactions with neutral helium in a dense, cool helium-rich atmosphere. We found that the abundance of heavy elements varied between spectra taken 10 years apart but we could not establish a time-scale for these variations; such variations may be linked to surface abundance variations in the magnetized atmosphere. Finally, we show that volume-limited samples reveal that about 40 per cent of DZ white dwarfs with effective temperatures below 7000 K are magnetic.
Metal-polluted white dwarfs (DZ stars) are characterized by a helium-rich atmosphere contaminated by heavy elements traces originating from accreted rocky planetesimals. As a detailed spectroscopic ...analysis of those objects can reveal the composition of the accreted debris, there is a great interest in developing accurate DZ atmosphere models. However, the coolest DZ white dwarfs are challenging to model due to the fluidlike density of their atmospheres. Under such extreme conditions, spectral absorption lines are heavily broadened by interactions with neutral helium, and it is no longer justified to use the conventional Lorentzian profiles. In this work, we determine the theoretical profiles of the Ca i resonance line (the most prominent spectral line for the coolest DZ white dwarfs) in the dense atmospheres of cool DZ white dwarfs. To do so, we use a unified theory of collisional line profiles and accurate ab initio potential energies and transition dipole moments for the CaHe molecule. We present the resulting profiles for the full range of temperatures and helium densities relevant for the modeling of cool, metal-polluted white dwarfs (from 3000 to 6000 K and from 1021 to 1023 cm−3). We also implement these new profiles in our atmosphere models and show that they lead to improved fits to the Ca i resonance line of the coolest DZ white dwarfs.
Context. Line shapes of the magnesium resonance lines in white dwarf spectra are determined by the properties of magnesium atoms and the structure of the white dwarf atmosphere. Through their ...blanketing effect, these lines have a dominant influence on the model structure and thus on the determination from the spectra of other physical parameters that describe the stellar atmosphere and elemental abundances. Aims. In continuation of previous work on Mg+He lines in the UV, we present theoretical profiles of the resonance line of neutral Mg perturbed by He at the extreme density conditions found in the cool largely transparent atmosphere of DZ white dwarfs. Methods. We accurately determined the broadening of Mg by He in a unified theory of collisional line profiles using ab initio calculations of MgHe potential energies and transition matrix elements among the singlet electronic states that are involved for the observable spectral lines. Results. We computed the shapes and line parameters of the Mg lines and studied their dependence on helium densities and temperatures. We present results over the full range of temperatures from 4000 to 12 000 K needed for input to stellar spectra models. Atmosphere models were constructed for a range of effective temperatures and surface gravities typical for cool DZ white dwarfs. We present synthetic spectra tracing the behavior of the Mg resonance line profiles under the low temperatures and high gas pressures prevalent in these atmospheres. Conclusions. The determination of accurate opacity data of magnesium resonance lines together with an improved atmosphere model code lead to a good fit of cool DZ white dwarf stars. The broadening of spectral lines by helium needs to be understood to accurately determine the H/He and Mg/He abundance ratio in DZ white dwarf atmospheres. We emphasize that no free potential parameters or ad hoc adjustments were used to calculate the line profiles.
The influence of inelastic hydrogen atom collisions on non-LTE spectral line formation has been, and remains to be, a significant source of uncertainty for stellar abundance analyses, due to the ...difficulty in obtaining accurate data for low-energy atomic collisions either experimentally or theoretically. For lack of a better alternative, the classical “Drawin formula” is often used. Over recent decades, our understanding of these collisions has improved markedly, predominantly through a number of detailed quantum mechanical calculations. In this paper, the Drawin formula is compared with the quantum mechanical calculations both in terms of the underlying physics and the resulting rate coefficients. It is shown that the Drawin formula does not contain the essential physics behind direct excitation by H atom collisions, the important physical mechanism being quantum mechanical in character. Quantitatively, the Drawin formula compares poorly with the results of the available quantum mechanical calculations, usually significantly overestimating the collision rates by amounts that vary markedly between transitions.
The purpose of the present work is to make an exhaustive study of the line shape of the triplet 3p-4s Mg line (Mgb triplet), which is perturbed by He in the extreme physical conditions found in the ...cool atmosphere of DZ white dwarfs. This study is undertaken by inferring both a unified theory of spectral line broadening and ab initio potential energies. Cool white dwarfs require a specific treatment for line broadening owing to the high helium densities that are involved. Beyond the conventional symmetrical Lorentzian core at low density, we show that the line profiles are asymmetrical and have significant additional contributions on the short wavelength side. This blue asymmetry is a consequence of low maxima in the corresponding Mg-He potential energy difference curves at short and intermediate internuclear distances. The new profiles are shown to provide a good fit to an SDSS (Sloan Digital Sky Survey) observation.
Rate coefficients for inelastic Na+H collisions are calculated for all transitions between the ten levels up to and including the ionic state (ion-pair production), namely ...Na(3s,3p,4s,3d,4p,5s,4d,4f,5p)+H(1s) and Na++H-. The calculations are based on recent full quantum scattering cross-section calculations. The data are needed for non-LTE applications in cool astrophysical environments, especially cool stellar atmospheres, and are presented for a temperature range of 500-8000 K. From consideration of the sensitivity of the cross-sections to input quantum chemical data and the results of different methods for the scattering calculations, a measure of the possible uncertainties in the rate coefficients is estimated.
In the present work, all adiabatic potential energy curves, spectroscopic constants and dipole moments of CaH
+
molecular ion dissociating below the ionic limit Ca
2+
H
−
are presented. These curves ...are determined by an ab initio approach involving a non-empirical pseudo-potential for the Ca core, core-valence correlation accounted in operator form with two types of core polarization potentials (CPP) and full valence Configuration Interaction. The molecule is thus treated as a two-electron system. The potential energy curves and the spectroscopic constants are presented. In addition, the permanent and transition dipole moments are calculated for most of the states and reveal the underlying ionic states. A rather good agreement with the available theoretical works in the literature is obtained for the spectroscopic constants of the lowest states of the CaH
+
molecule.
Display omitted
•N2/He+ interaction comes into play when cold rare-gas plasmas leave the interior of their generators and get into contact with environment species.•Cold-plasmas exhibit useful ...properties in medicinal applications, surface treatment, food storage, and/or spacecraft propulsion.•The main goal of the present study is to perform a comprehensive series of test electronic-structure calculations using various ab initio approaches in order to (mainly) assess the reliability of employed methods and their practicability in dynamical calculations.
As a first step towards realistic modeling of transport properties of the N2+ ion in helium gas, detailed investigations of the electronic structure of the ion as well as the N2+/He collision complex have been performed with the main focus on computational efficiency and accuracy. A broad range of correlation consistent basis sets and representative orbital spaces have been considered for both computationally cheap multi-configuration self-consistent field (MCSCF) methods and benchmark multi-reference configuration interaction (MRCI) approaches. It has been found that the computationally advantageous MCSCF approach, even if combined with basis sets of moderate sizes, leads to an acceptable accuracy and is thus suitable for direct dynamics simulations of N2+/He collisions to be addressed in subsequent works.