Cyanopolyynes in hot cores: modelling G305.2+0.2

We present results from a time-dependent gas-phase chemical model of a hot core based on the physical conditions of G305.2+0.2. While the cyanopolyyne HC3N has been observed in hot cores, the longer chained species, HC5N, HC7N and HC9N, have not been considered as the typical hot-core species. We present results which show that these species can be formed under hot core conditions. We discuss the important chemical reactions in this process and, in particular, show that their abundances are linked to the parent species acetylene which is evaporated from icy grain mantles. The cyanopolyynes show promise as ‘chemical clocks’ which may aid future observations in determining the age of hot core sources. The abundance of the larger cyanopolyynes increases and decreases over relatively short time-scales, ~10^2.5 yr. We present results from a non-local thermodynamic equilibrium statistical equilibrium excitation model as a series of density, temperature and column density dependent contour plots which show both the line intensities and several line ratios. These aid in the interpretation of spectral-line data, even when there is limited line information available. In particular, non-detections of HC5N and HC7N in Walsh et al. are analysed and discussed.

The infrared fundamental intensities of some cyanopolyynes

Some cyanopolyynes, HCnN (n = 1, 3, ... , 17), are investigated by means of calculations at the MP2/cc-pVTZ and CCSD/cc-pVDZ levels. Although the MP2/cc-pVTZ results for geometries and molecular dipole moments are encouraging, the CCSD/cc-pVDZ level was superior for the study of infrared fundamental intensities. The main bands are also analyzed with a charge-charge flux-dipole flux (CCFDF) partition model based on quantities given by the Quantum Theory of Atoms in Molecules (QTAIM). The intensity of vibrations corresponding to the stretching of CH bonds (3471-3473 cm(-1)) increases in line with the number of carbon atoms (from 61 to 146 km mol(-1) between HCN and HC13N). This increase is due to the charge flux contribution while the other contributions remain roughly unaltered except for HCN. Moreover, the hydrogen atom loses an almost constant amount of electronic charge during the CH bond enlargement and a small fraction of this charge spreads to atoms farther and farther away from hydrogen as the molecule size increases. The band associated with the doubly degenerate CH bending vibrations (643-732 cm(-1)) presents approximately the same intensity in all the studied cyanopolyynes (from 67 to 76 km mol(-1)). The CCFDF/QTAIM contributions are also nearly the same for these bending modes in HC5N and larger systems. The intensity of the mode mostly identified as CN stretching (around 2378-2399 cm(-1) except for HCN) increases from HCN up to HC7N (from 0.3 to 83 km mol(-1)) and nearly stabilizes around 80-90 km mol(-1) for larger systems. The CCFDF/QTAIM contributions for this mode also change significantly up to HC7N and remain almost constant in larger systems. We also observed the appearing of a very relevant band between 2283 and 2342 cm(-1). This mode is mainly associated with the symmetric stretching of CC triple bonds near the molecule center and exhibits large charge fluxes while the other contributions are almost negligible in the largest cyanopolyynes. The two vibrational bands associated with the smallest frequencies are also studied and extrapolation equations are suggested to predict their positions in larger cyanopolyynes. (C) 2012 Elsevier B.V. All rights reserved.

High-temperature kinetics of the reaction between CN and hydrocarbons using a novel high-enthalpy flow tube

More than 70 molecules of varied nature have been identified in the envelopes of carbon-rich stars through their spectral fingerprints in the microwave or far infrared regions. Many of them are carbon chain molecules and radicals, and a significant number are unique to the circumstellar medium. The determination of relevant laboratory kinetics data is critical to keep up with the development of the high spectral and spatial resolution observations and of the refinement of chemical models. Neutralneutral reactions of the CN radical with unsaturated hydrocarbons could be a dominant route in the formation of cyanopolyynes, even at low temperatures and deserve a detailed laboratory investigation. The approach we have developed aims to bridge the temperature gap between resistively heated flow tubes and shock tubes. The present kinetic measurements are obtained using a new reactor combining a high-enthalpy source with a flow tube and a pulsed laser photolysislaser-induced fluorescence system to probe the undergoing chemical reactions. The high-enthalpy flow tube has been used to measure the rate constant of the reaction of the CN radical with propane (C3H8), propene (C3H6), allene (C3H4), 1,3-butadiene (1,3-C4H6), and 1-butyne (C4H6) over a temperature range extending from 300 to 1200 K. All studied reactions of CN with unsaturated hydrocarbons are rapid, with rate coefficients greater than 10-10 cm3 center dot molecule-1 center dot s-1 and exhibit slight negative temperature dependence above room temperature. (c) 2012 Wiley Periodicals, Inc. Int J Chem Kinet 44: 753766, 2012

Density-Enhanced Gas and Dust Shells in a New Chemical Model for IRC+10216

A new chemical model is presented for the carbon-rich circumstellar envelope (CSE) of the asymptotic giant branch star IRC+10216. The model includes shells of matter with densities that are enhanced relative to the surrounding circumstellar medium. The chemical model uses an updated reaction network including reactions from the RATE06 database and a more detailed anion chemistry. In particular, new mechanisms are considered for the formation of CN-, C3N-, and C2H-, and for the reactions of hydrocarbon anions with atomic nitrogen and with the most abundant cations in the CSE. New reactions involving H- are included which result in the production of significant amounts of C2H- and CN- in the inner envelope. The calculated radial molecular abundance profiles for the hydrocarbons C2H, C4H, and C6H and the cyanopolyynes HC3N and HC5N show narrow peaks which are in better agreement with observations than previous models. Thus, the narrow rings observed in molecular microwave emission surrounding IRC+10216 are interpreted as arising in regions of the envelope where the gas and dust densities are greater than the surrounding circumstellar medium. Our models show that CN- and C2H- may be detectable in IRC+10216 despite the very low theorized radiative electron attachment rates of their parent neutral species. We also show that magnesium isocyanide (MgNC) can be formed in the outer envelope through radiative association involving Mg+ and the cyanopolyyne species.

Large molecules in the envelope surrounding IRC+10216

A new chemical model of the circumstellar envelope surrounding the carbon-rich star IRC+10216 is developed that includes carbon-containing molecules with up to 23 carbon atoms. The model consists of 3851 reactions involving 407 gas-phase species. Sizeable abundances of a variety of large molecules - including carbon clusters, unsaturated hydrocarbons and cyanopolyynes - have been calculated. Negative molecular ions of chemical formulae C-n(-) and CnH- (7 less than or equal to n less than or equal to 23) exist in considerable abundance, with peak concentrations at distances from the central star somewhat greater than their neutral counterparts. The negative ions might be detected in radio emission, or even in the optical absorption of background field stars. The calculated radial distributions of the carbon-chain CnH radicals are looked at carefully and compared with interferometric observations.

Discovery of Interstellar Anions in Cepheus and Auriga

We report the detection of microwave emission lines from the hydrocarbon anion C6H- and its parent neutral C6H in the star-forming region L1251A (in Cepheus), and the pre-stellar core L1512 (in Auriga). The carbon-chain-bearing species C4H, HC3N, HC5N, HC7N and C3S are also detected in large abundances. The observations of L1251A constitute the first detections of anions and long- chain polyynes and cyanopolyynes (with more than 5 carbon atoms) in the Cepheus Flare star- forming region, and the first detection of anions in the vicinity of a protostar outside of the Taurus molecular cloud complex, highlighting a wider importance for anions in the chemistry of star formation. Rotational excitation temperatures have been derived from the HC3N hyperfine structure lines, and are found to be 6.2 K for L1251A and 8.7 K for L1512. The anion-to-neutral ratios are 3.6% and 4.1%, respectively, which are within the range of values previously observed in the interstellar medium, and suggest a relative uniformity in the processes governing anion abundances in different dense interstellar clouds. This research contributes towards the growing body of evidence that carbon chain anions are relatively abundant in interstellar clouds throughout the Galaxy, but especially in the regions of relatively high density and high depletion surrounding pre-stellar cores and young, embedded protostars.

New detections of HC5N towards hot cores associated with 6.7 GHz methanol masers

We present new detections of cyanodiacetylene (HC5N) toward hot molecular cores, observed with the Tidbinbilla 34 m radio telescope (DSS–34). In a sample of 79 hot molecular cores, HC5N was detected towards 35. These results are counter to the expectation that long chain cyanopolyynes, such as HC5N, are not typically found in hot molecular cores, unlike their shorter chain counterpart HC3N. However it is consistent with recent models which suggest HC5N may exist for a limited period during the evolution of hot molecular cores.

Large-scale coupled-cluster calculations

Coupled-cluster theory provides one of the most successful concepts in electronic-structure theory. This work covers the parallelization of coupled-cluster energies, gradients, and second derivatives and its application to selected large-scale chemical problems, beside the more practical aspects such as the publication and support of the quantum-chemistry package ACES II MAB and the design and development of a computational environment optimized for coupled-cluster calculations. The main objective of this thesis was to extend the range of applicability of coupled-cluster models to larger molecular systems and their properties and therefore to bring large-scale coupled-cluster calculations into day-to-day routine of computational chemistry. A straightforward strategy for the parallelization of CCSD and CCSD(T) energies, gradients, and second derivatives has been outlined and implemented for closed-shell and open-shell references. Starting from the highly efficient serial implementation of the ACES II MAB computer code an adaptation for affordable workstation clusters has been obtained by parallelizing the most time-consuming steps of the algorithms. Benchmark calculations for systems with up to 1300 basis functions and the presented applications show that the resulting algorithm for energies, gradients and second derivatives at the CCSD and CCSD(T) level of theory exhibits good scaling with the number of processors and substantially extends the range of applicability. Within the framework of the ’High accuracy Extrapolated Ab initio Thermochemistry’ (HEAT) protocols effects of increased basis-set size and higher excitations in the coupled- cluster expansion were investigated. The HEAT scheme was generalized for molecules containing second-row atoms in the case of vinyl chloride. This allowed the different experimental reported values to be discriminated. In the case of the benzene molecule it was shown that even for molecules of this size chemical accuracy can be achieved. Near-quantitative agreement with experiment (about 2 ppm deviation) for the prediction of fluorine-19 nuclear magnetic shielding constants can be achieved by employing the CCSD(T) model together with large basis sets at accurate equilibrium geometries if vibrational averaging and temperature corrections via second-order vibrational perturbation theory are considered. Applying a very similar level of theory for the calculation of the carbon-13 NMR chemical shifts of benzene resulted in quantitative agreement with experimental gas-phase data. The NMR chemical shift study for the bridgehead 1-adamantyl cation at the CCSD(T) level resolved earlier discrepancies of lower-level theoretical treatment. The equilibrium structure of diacetylene has been determined based on the combination of experimental rotational constants of thirteen isotopic species and zero-point vibrational corrections calculated at various quantum-chemical levels. These empirical equilibrium structures agree to within 0.1 pm irrespective of the theoretical level employed. High-level quantum-chemical calculations on the hyperfine structure parameters of the cyanopolyynes were found to be in excellent agreement with experiment. Finally, the theoretically most accurate determination of the molecular equilibrium structure of ferrocene to date is presented.