G. Mulas

Electronic and optical properties of families of polycyclic aromatic hydrocarbons: A systematic (time-dependent) density functional theory study

Abstract Homologous classes of polycyclic aromatic hydrocarbons (PAHs) in their crystalline state are among the most promising materials for organic opto-electronics. Following previous works on oligoacenes we present a systematic comparative study of the electronic, optical, and transport properties of oligoacenes, phenacenes, circumacenes, and oligorylenes. Using density functional theory (DFT) and time-dependent DFT we computed: (i) electron affinities and first ionization energies; (ii) quasiparticle correction to the highest occupied molecular orbital (HOMO)-lowest unoccupied molecular orbital (LUMO) gap; (iii) molecular reorganization energies and (iv) electronic absorption spectra of neutral and ±1 charged systems. The excitonic effects are estimated by comparing the optical gap and the quasiparticle corrected HOMO–LUMO energy gap. For each molecular property computed, general trends as a function of molecular size and charge state are discussed. Overall, we find that circumacenes have the best transport properties, displaying a steeper decrease of the molecular reorganization energy at increasing sizes, while oligorylenes are much more efficient in absorbing low-energy photons in comparison to the other classes.

Electronic absorption spectra of PAHs up to vacuum UV Towards a detailed model of interstellar PAH photophysics

We computed the absolute photo-absorption cross-sections up to the vacuum ultaviolet (VUV) of a total of 20 Polycyclic Aromatic Hydrocarbons (PAHs) and their respective cations, ranging in size from naphthalene (C10H8 )t o di- coronylene (C48H20). We used an implementation in real time and real space of the Time-Dependent Density Functional Theory (TD-DFT), an approach which was proven to yield accurate results for conjugated molecules such as benzene. Concerning the low-lying excited states of π ∗ ← π character occurring in the near-IR, visible and near-UV spectral range, the computed spectra are in good agreement with the available experimental data, predicting vertical excitation energies precise to within a few tenths of eV, and the corresponding oscillator strengths to within experimental errors, which are indeed the typical accuracies currently achievable by TD-DFT. We find that PAH cations, like their parent molecules, display strong π ∗ ← π electronic transitions in the UV, a piece of information which is particularly useful since a limited amount of laboratory data is available concerning the absorption properties of PAH ions in this wavelength range. Moreover, a detailed discussion of the UV-VUV properties of both neutral and cation species is presented. Concerning neutrals, the agreement with existing laboratory data is very good, the specific TD-DFT implementation used in this work apparently being able to reproduce the overall far-UV behaviour, including the broad absorption peak dominated by σ ∗ ← σ transitions, which matches well both in position and width. The implications of these results are discussed in conjunction with the contribution PAH-like molecules are expected to give to the interstellar extinction curve.

Theoretical electron affinities of PAHs and electronic absorption spectra of their mono-anions

We present theoretical electron affinities, calculated as total energy differences, for a large sample of polycyclic aromatic hydrocarbons (PAHs), ranging in size from azulene (C10H8) to dicoronylene (C48H20). For 20 out of 22 molecules under study we obtained electron affinity values in the range 0.4-2.0 eV, showing them to be able to accept an additional electron in their LUMO π � orbital. For the mono-anions we computed the absolute photo-absorption cross-sections up to the vacuum ultraviolet (VUV) using an implementation in real time and real space of the Time-Dependent Density Functional Theory (TD-DFT), an approach which has already been proven to yield accurate results for neutral and cationic PAHs. Comparison with available experimental data hints that this is the case for mono-anions as well. We find that PAH anions, like their parent molecules and the corresponding cations, display strong π ∗ ← π electronic transitions in the UV. The present results provide a quantitative foundation to estimate the fraction of specific PAHs which can be singly negatively charged in various interstellar environments, to simulate their photophysics in detail and to evaluate their contribution to the interstellar extinction curve.

The role of the charge state of PAHs in ultraviolet extinction

Aims. We explore the relation between the charge state of polycyclic aromatic hydrocarbons (PAHs) and the extinction curve morphology. Methods. We fit extinction curves with a dust model including core-mantle spherical particles of mixed chemical composition (silicate core, sp 2 and sp 3 carbonaceous layers), and an additional molecular component. We use exact methods to calculate the extinction due to classical particles and accurate computed absorption spectra of PAHs in different charge states, for the contribution due to the molecular component, along five different lines of sight. Results. A combination of classical dust particles and mixtures of real PAHs satisfactorily matches the observed interstellar extinction curves. Variations of the spectral properties of PAHs in different charge states produce changes consistent with the varying relative strengths of the bump and non-linear far-UV rise.

Theoretical evaluation of PAH dication properties

Aims. We present a systematic theoretical study on 40 polycyclic aromatic hydrocarbons dications (PAHs++) containing up to 66 carbon atoms. Methods. We performed our calculations using well established quantum-chemical techniques in the framework of the density functional theory (DFT) to obtain the electronic ground-state properties, and of the time-dependent DFT (TD-DFT) to evaluate the excited\textendash state properties. Results. For each PAH++ considered we computed the absolute visible-UV photo-absorption cross-section up to about 30 eV. We also evaluated their vibrational properties and compared them to those of the corresponding neutral and singly-ionised species. We estimated the adiabatic and vertical second ionisation energy DeltaI through total energy differences. Conclusions. The DeltaI values obtained fall in the energy range 8-13 eV, confirming that PAHs could reach the doubly-ionised state in HI regions. The total integrated IR absorption cross-sections show a marked increase upon ionization, being on the average about two and five times larger for PAHs++ than for PAHs+ and PAHs, respectively. The visible-UV photo-absorption cross-sections for the 0, +1 and +2 charge-states show comparable features but PAHs++ are found to absorb slightly less than their parent neutrals and singly ionized species between ~7 and ~12 eV. Combining these pieces of information we found that PAHs++ should actually be more stable against photodissociation than PAHs and PAHs+,_if_ dissociation tresholds are nearly unchanged by ionization.

Dehydrogenated polycyclic aromatic hydrocarbons and UV bump (Research Note)

Context. Recent calculations have shown that the UV bump at about 217.5 nm in the extinction curve can be explained by a complex mixture of polycyclic aromatic hydrocarbons (PAHs) in several ionisation states. Other studies proposed that the carriers are a restricted population made of neutral and singly-ionised dehydrogenated coronene molecules (C24Hn, n ≤ 3), in line with models of the hydrogenation state of interstellar PAHs predicting that medium-sized species are highly dehydrogenated. Aims. To assess the observational consequences of the latter hypothesis we have undertaken a systematic theoretical study of the electronic spectra of dehydrogenated PAHs. We use our first results to see whether such spectra show strong general trends upon dehydrogenation. Methods. We performed calculations using state-of-the-art techniques in the framework of the density functional theory (DFT) to obtain the electronic ground-state geometries, and of the time-dependent DFT to evaluate the electronic excited-state properties. Results. We computed the absorption cross-section of the species C24Hn (n = 12, 10, 8, 6, 4, 2, 0) in their neutral and cationic charge-states. Similar calculations were performed for other PAHs and their fully dehydrogenated counterparts. Conclusions. π-electron energies are always found to be strongly affected by dehydrogenation. In all cases we examined, progressive dehydrogenation translates into a correspondingly progressive blue shift of the main electronic transitions. In particular, the π → π ∗ collective resonance becomes broader and bluer with dehydrogenation. Its calculated energy position is therefore predicted to fall in the gap between the UV bump and the far-UV rise of the extinction curve. Since this effect appears to be systematic, it poses a tight observational limit on the column density of strongly dehydrogenated medium-sized PAHs.

Estimated IR and phosphorescence emission fluxes for specific polycyclic aromatic hydrocarbons in the Red Rectangle

Following tentative identification of the blue luminescence in the Red Rectangle by Vijh et al., we computed absolute fluxes for the vibrational IR emission and phosphorescence bands of three small polycyclic aromatic hydrocarbons. The calculated IR spectra were compared with available ISO observations. A subset of the emission bands are predicted to be observable using presently available facilities, and can be used for an immediate, independent, discriminating test of their alleged presence in this well-known astronomical object.

Search for corannulene (C20H10) in the Red Rectangle

Polycyclic Aromatic Hydrocarbons (PAHs) are widely accepted as the carriers of the Aromatic Infrared Bands (AIBs), but an unambiguous identification of any specific interstellar PAH is still missing. For polar PAHs, pure rotational transitions can be used as spectral fingerprints for identification. Combining dedicated experiments, detailed simulations and observations, we explore d the mm wavelength domain to search for specific rotational transitions of corannulene (C20H10). We performed high-resolution spectroscopic measurements and a simulation of the emission spectrum of ultraviolet-excited C20H10 in the environment of the Red Rectangle (RR), calculating its synthetic rotational spectrum. Based on these results, we conducted a first observational campaign at the IRAM 30-m telescope towards this source to search for several high-J rotational transitions of C20H10. The laboratory detection of the J = 112 134 transition at 137.615 GHz. Comparing the noise level with the synthetic spectrum, we are able to estimate an upper limit to the fraction of carbon locked in corannulene of about 1.0 x 10(-5) relative to the total abundance of carbon in PAHs. The sensitivity achieved in this work shows that radio spectroscopy can be a powerful tool to search for polar PAHs. We compare this upper limit with models for the PAH size distribution, emphasizing that small PAHs are much less abundant than predicted. We show that this cannot be explained by destruction but is more likely related to the chemistry of their formation in the environment of the RR.

Efficient calculation of van der Waals dispersion coefficients with time-dependent density functional theory in real time: Application to polycyclic aromatic hydrocarbons

The van der Waals dispersion coefficients of a set of polycyclic aromatic hydrocarbons, ranging in size from the single-cycle benzene to circumovalene (C(66)H(20)), are calculated with a real-time propagation approach to time-dependent density functional theory (TDDFT). In the nonretarded regime, the Casimir-Polder integral is employed to obtain C(6), once the dynamic polarizabilities have been computed at imaginary frequencies with TDDFT. On the other hand, the numerical coefficient that characterizes the fully retarded regime is obtained from the static polarizabilities. This ab initio strategy has favorable scaling with the size of the system--as demonstrated by the size of the reported molecules--and can be easily extended to obtain higher order van der Waals coefficients.

Mass motions and chromospheres of RGB stars in the globular cluster NGC 2808

We present the results of the first observations, taken with FLAMES during Science Verification, of red giant branch (RGB) stars in the globular cluster NGC 2808. A total of 137 stars was observed, of which 20 at high resolution (R = 47 000) with UVES and the others at lower resolution (R = 19 000-29 000) with GIRAFFE in MEDUSA mode, monitoring ∼3 mag down from the RGB tip. Spectra were taken of the Ha, Na I D and Ca II H and K lines. This is by far the largest and most complete collection of such data in globular cluster giants, both for the number of stars observed within one cluster, and for monitoring all the most important optical diagnostics of chromospheric activity/mass motions. Evidence of mass motions in the atmospheres was searched from asymmetry in the profiles and coreshifts of the Ha, Na I D and Ca II K lines, as well as from Ha emission wings. We have set the detection thresholds for the onset of Ha emission, negative Na D 2 coreshifts and negative K 3 coreshifts at log L/L ○. ∼ 2.5, 2.9 and 2.8, respectively. These limits are significantly fainter than the results found by nearly all previous studies. Also the fraction of stars where these features have been detected has been increased significantly with respect to the previous studies. Our observations confirm the widespread presence of chromospheres among globular cluster giants, as it was found among Population I red giants. Some of the above diagnostics suggest clearly the presence of outward mass motions in the atmosphere of several stars.