Chemistry and Distribution of Daughter Species in the Circumstellar Envelopes of O-Rich AGB Stars

Context. Thanks to the advent of Herschel and ALMA, new high-quality observations of molecules present in the circumstellar envelopes of asymptotic giant branch (AGB) stars are being reported that reveal large differences from the existing chemical models. New molecular data and more comprehensive models of the chemistry in circumstellar envelopes are now available.
Aims: The aims are to determine and study the important formation and destruction pathways in the envelopes of O-rich AGB stars and to provide more reliable predictions of abundances, column densities, and radial distributions for potentially detectable species with physical conditions applicable to the envelope surrounding IK Tau.
Methods: We use a large gas-phase chemical model of an AGB envelope including the effects of CO and N2 self-shielding in a spherical geometry and a newly compiled list of inner-circumstellar envelope parent species derived from detailed modeling and observations. We trace the dominant chemistry in the expanding envelope and investigate the chemistry as a probe for the physics of the AGB phase by studying variations of abundances with mass-loss rates and expansion velocities.
Results: We find a pattern of daughter molecules forming from the photodissociation products of parent species with contributions from ion-neutral abstraction and dissociative recombination. The chemistry in the outer zones differs from that in traditional PDRs in that photoionization of daughter species plays a significant role. With the proper treatment of self-shielding, the N → N2 and C+→ CO transitions are shifted outward by factors of 7 and 2, respectively, compared with earlier models. An upper limit on the abundance of CH4 as a parent species of (≲2.5 × 10-6 with respect to H2) is found for IK Tau, and several potentially observable molecules with relatively simple chemical links to other parent species are determined. The assumed stellar mass-loss rate, in particular, has an impact on the calculated abundances of cations and the peak-abundance radius of both cations and neutrals: as the mass-loss rate increases, the peak abundance of cations generally decreases and the peak-abundance radius of all species moves outwards. The effects of varying the envelope expansion velocity and cosmic-ray ionization rate are not as significant.

Modelling of Ultraviolet and Vissible Emissions in Comets

Comets are the spectacular objects in the night sky since the dawn of mankind. Due to their giant apparitions and enigmatic behavior, followed by coincidental calamities, they were termed as notorious and called as `bad omens'. With a systematic study of these objects modern scienti c community understood that these objects are part of our solar system. Comets are believed to be remnant bodies of at the end of evolution of solar system and possess the material of solar nebula. Hence, these are considered as most pristine objects which can provide the information about the conditions of solar nebula. These are small bodies of our solar system, with a typical size of about a kilometer to a few tens of kilometers orbiting the Sun in highly elliptical orbits. The solid body of a comet is nucleus which is a conglomerated mixture of water ice, dust and some other gases. When the cometary nucleus advances towards the Sun in its orbit the ices sublimates and produces the gaseous envelope around the nucleus which is called coma. The gravity of cometary nucleus is very small and hence can not in uence the motion of gases in the cometary coma. Though the cometary nucleus is a few kilometers in size they can produce a transient, extensive, and expanding atmosphere with size several orders of magnitude larger in space. By ejecting gas and dust into space comets became the most active members of the solar system. The solar radiation and the solar wind in uences the motion of dust and ions and produces dust and ion tails, respectively. Comets have been observed in di erent spectral regions from rocket, ground and space borne optical instruments. The observed emission intensities are used to quantify the chemical abundances of di erent species in the comets. The study of various physical and chemical processes that govern these emissions is essential before estimating chemical abundances in the coma. Cameron band emission of CO molecule has been used to derive CO2 abundance in the comets based on the assumption that photodissociation of CO2 mainly produces these emissions. Similarly, the atomic oxygen visible emissions have been used to probe H2O in the cometary coma. The observed green ([OI] 5577 A) to red-doublet emission ([OI] 6300 and 6364 A) ratio has been used to con rm H2O as the parent species of these emissions. In this thesis a model is developed to understand the photochemistry of these emissions and applied to several comets. The model calculated emission intensities are compared with the observations done by space borne instruments like International Ultraviolet Explorer (IUE) and Hubble Space Telescope (HST) and also by various ground based telescopes.

Electron-Nuclear Correlation in Laser-Driven Diatomic Molecules

The interaction of short intense laser pulses with atoms/molecules produces a multitude of highly nonlinear processes requiring a non-perturbative treatment. Detailed study of these highly nonlinear processes by numerically solving the time-dependent Schrodinger equation becomes a daunting task when the number of degrees of freedom is large. Also the coupling between the electronic and nuclear degrees of freedom further aggravates the computational problems. In the present work we show that the time-dependent Hartree (TDH) approximation, which neglects the correlation effects, gives unreliable description of the system dynamics both in the absence and presence of an external field. A theoretical framework is required that treats the electrons and nuclei on equal footing and fully quantum mechanically. To address this issue we discuss two approaches, namely the multicomponent density functional theory (MCDFT) and the multiconfiguration time-dependent Hartree (MCTDH) method, that go beyond the TDH approximation and describe the correlated electron-nuclear dynamics accurately. In the MCDFT framework, where the time-dependent electronic and nuclear densities are the basic variables, we discuss an algorithm to calculate the exact Kohn-Sham (KS) potentials for small model systems. By simulating the photodissociation process in a model hydrogen molecular ion, we show that the exact KS potentials contain all the many-body effects and give an insight into the system dynamics. In the MCTDH approach, the wave function is expanded as a sum of products of single-particle functions (SPFs). The MCTDH method is able to describe the electron-nuclear correlation effects as the SPFs and the expansion coefficients evolve in time and give an accurate description of the system dynamics. We show that the MCTDH method is suitable to study a variety of processes such as the fragmentation of molecules, high-order harmonic generation, the two-center interference effect, and the lochfrass effect. We discuss these phenomena in a model hydrogen molecular ion and a model hydrogen molecule. Inclusion of absorbing boundaries in the mean-field approximation and its consequences are discussed using the model hydrogen molecular ion. To this end, two types of calculations are considered: (i) a variational approach with a complex absorbing potential included in the full many-particle Hamiltonian and (ii) an approach in the spirit of time-dependent density functional theory (TDDFT), including complex absorbing potentials in the single-particle equations. It is elucidated that for small grids the TDDFT approach is superior to the variational approach.

Time-resolved gas-phase kinetic and quantum chemical studies of reactions of silylene with chlorine-containing species. 2. CH3Cl

Time-resolved kinetic studies of the reaction of silylene, SiH2, generated by laser flash photolysis of both silacyclopent-3-ene and phenylsilane, have been carried out to obtain second-order rate constants for its reaction with CH3Cl. The reaction was studied in the gas phase at six temperatures in the range 294-606 K. The second-order rate constants gave a curved Arrhenius plot with a minimum value at T approximate to 370 K. The reaction showed no pressure dependence in the presence of up to 100 Torr SF6. The rate constants, however, showed a weak dependence on laser pulse energy. This suggests an interpretation requiring more than one contributing reaction pathway to SiH2 removal. Apart from a direct reaction of SiH2 with CH3Cl, reaction of SiH2 with CH3 (formed by photodissociation of CH3Cl) seems probable, with contributions of up to 30% to the rates. Ab initio calculations (G3 level) show that the initial step of reaction of SiH2 with CH3Cl is formation of a zwitterionic complex (ylid), but a high-energy barrier rules out the subsequent insertion step. On the other hand, the Cl-abstraction reaction leading to CH3 + ClSiH2 has a low barrier, and therefore, this seems the most likely candidate for the main reaction pathway of SiH2 with CH3Cl. RRKM calculations on the abstraction pathway show that this process alone cannot account for the observed temperature dependence of the rate constants. The data are discussed in light of studies of other silylene reactions with haloalkanes.

Spectroscopy of highly excited vibrational states of HCN in its ground electronic state

An experimental technique based on a scheme of vibrationally mediated photodissociation has been developed and applied to the spectroscopic study of highly excited vibrational states in HCN, with energies between 29 000 and 30 000 cm(-1). The technique consists of four sequential steps: in the first one, a high power laser is used to vibrationally excite the sample to an intermediate state, typically (0,0,4), the nu(3) mode being approximately equivalent to the C-H stretching vibration. Then a second laser is used to search for transitions between this intermediate state and highly vibrationally excited states. When one of these transitions is found, HCN molecules are transferred to a highly excited vibrational state. Third, a ultraviolet laser photodissociates the highly excited molecules to produce H and CN radicals in its A (2)Pi electronic state. Finally, a fourth laser (probe) detects the presence of the CN(A) photofragments by means of an A-->B-->X laser induced fluorescence scheme. The spectra obtained with this technique, consisting of several rotationally resolved vibrational bands, have been analyzed. The positions and rotational parameters of the states observed are presented and compared with the results of a state-of-the-art variational calculation. (C) 2004 American Institute of Physics.

Highly excited vibrational states of HCN around 30 000 cm(-1)

This article describes the analysis and interpretation of rovibrational spectra involving highly excited vibrational states in the molecule of HCN. The spectra were obtained by means of the vibrationally mediated photodissociation technique. Analysis of the spectra revealed four bands with Sigma-Sigma structures that, once fitted, provided the energies and rotational constants of four new, highly excited vibrational states in the region of the potential energy surface near and above 30 000 cm(-1). All the states could be identified with the help of a state-of-the-art variational calculation. Together with the states already observed in previous works, eight highly excited states have so far been identified in this region. (c) 2006 American Institute of Physics.

CO migration in native and mutant myoglobin: atomistic simulations for the understanding of protein function

Molecular dynamics simulations of the events after the photodissociation of CO in the myoglobin mutant L29F in which leucine is replaced by phenylalanine are reported. Using both classical and mixed quantum-classical molecular dynamics calculations, we observed the rapid motion of CO away from the distal heme pocket to other regions of the protein, in agreement with recent experimental results. The experimentally observed and calculated infrared spectra of CO after dissociation are also in good agreement. We compared the results with data from simulations of WT myoglobin. As the time resolution of experimental techniques is increased, theoretical methods and models can be validated at the atomic scale by direct comparison with experiment.

Theoretical investigation of infrared spectra and pocket dynamics of photodissociated carbonmonoxy myoglobin

Molecular dynamics simulations of the photodissociated state of carbonmonoxy myoglobin (MbCO) are presented using a fluctuating charge model for CO. A new three-point charge model is fitted to high-level ab initio calculations of the dipole and quadrupole moment functions taken from the literature. The infrared spectrum of the CO molecule in the heme pocket is calculated using the dipole moment time autocorrelation function and shows good agreement with experiment. In particular, the new model reproduces the experimentally observed splitting of the CO absorption spectrum. The splitting of 3–7 cm−1 (compared to the experimental value of 10 cm−1) can be directly attributed to the two possible orientations of CO within the docking site at the edge of the distal heme pocket (the B states), as previously suggested on the basis of experimental femtosecond time-resolved infrared studies. Further information on the time evolution of the position and orientation of the CO molecule is obtained and analyzed. The calculated difference in the free energy between the two possible orientations (Fe···CO and Fe···OC) is 0.3 kcal mol−1 and agrees well with the experimentally estimated value of 0.29 kcal mol−1. A comparison of the new fluctuating charge model with an established fixed charge model reveals some differences that may be critical for the correct prediction of the infrared spectrum and energy barriers. The photodissociation of CO from the myoglobin mutant L29F using the new model shows rapid escape of CO from the distal heme pocket, in good agreement with recent experimental data. The effect of the protein environment on the multipole moments of the CO ligand is investigated and taken into account in a refined model. Molecular dynamics simulations with this refined model are in agreement with the calculations based on the gas-phase model. However, it is demonstrated that even small changes in the electrostatics of CO alter the details of the dynamics.

The reaction between silylene and ammonia: some gas-phase kinetic and quantum chemical studies

Time-resolved kinetic studies of the reaction of silylene, SiH2, generated by 193 nm laser flash photolysis of silacyclopent-3-ene, have been carried out in the presence of ammonia, NH3. Second order kinetics were observed. The reaction was studied in the gas phase at 10 Torr total pressure in SF6 bath gas at each of the three temperatures, 299, 340 and 400 K. The second order rate constants (laser pulse energy of 60 mJ/pulse) fitted the Arrhenius equation: log(k/cm3 molecule-1 s-1) = (-10.37 ± 0.17) + (0.36 ± 1.12 kJ mol-1)/RTln10 Experiments at other pressures showed that these rate constants were unaffected by pressure in the range 10-100 Torr, but showed small decreases in value at 3 and 1 Torr. There was also a weak intensity dependence, with rate constants decreasing at laser pulse energies of 30 mJ/pulse. Ab initio calculations at the G3 level of theory, show that SiH2 + NH3 should form an initial adduct (donor-acceptor complex), but that energy barriers are too great for further reaction of the adduct. This implies that SiH2 + NH3 should be a pressure dependent association reaction. The experimental data are inconsistent with this and we conclude that SiH2 decays are better explained by reaction of SiH2 with the amino radical, NH2, formed by photodissociation of NH3 at 193 nm. The mechanism of this previously unstudied reaction is discussed.

Modelling the warm H(2) infrared emission of the Helix nebula cometary knots

Molecular hydrogen emission is commonly observed in planetary nebulae. Images taken in infrared H(2) emission lines show that at least part of the molecular emission is produced inside the ionized region. In the best studied case, the Helix nebula, the H(2) emission is produced inside cometary knots (CKs), comet-shaped structures believed to be clumps of dense neutral gas embedded within the ionized gas. Most of the H(2) emission of the CKs seems to be produced in a thin layer between the ionized diffuse gas and the neutral material of the knot, in a mini-photodissociation region (mini-PDR). However, PDR models published so far cannot fully explain all the characteristics of the H(2) emission of the CKs. In this work, we use the photoionization code AANGABA to study the H(2) emission of the CKs, particularly that produced in the interface H(+)/H(0) of the knot, where a significant fraction of the H(2) 1-0 S(1) emission seems to be produced. Our results show that the production of molecular hydrogen in such a region may explain several characteristics of the observed emission, particularly the high excitation temperature of the H(2) infrared lines. We find that the temperature derived from H(2) observations, even of a single knot, will depend very strongly on the observed transitions, with much higher temperatures derived from excited levels. We also proposed that the separation between the H alpha and [N II] peak emission observed in the images of CKs may be an effect of the distance of the knot from the star, since for knots farther from the central star the [N II] line is produced closer to the border of the CK than H alpha.

Entropy Drives an Attached Water Molecule from the C- to N-Terminus on Protonated Proline

Results from infrared photodissociation (IRPD) spectroscopy and kinetics of singly hydrated, protonated proline indicate that the water molecule hydrogen bonds preferentially to the formally neutral carboxylic acid at low temperatures and at higher temperatures to the protonated N-terminus, which bears the formal charge. Hydration isomer populations obtained from IRPD kinetic data as a function of temperature are used to generate a van`t Hoff plot that reveals that C-terminal binding is enthalpically favored by 4.2-6.4 kJ/mol, whereas N-terminal binding is entropically favored by 31-43 J/(mol K), consistent with a higher calculated barrier for water molecule rotation at the C-terminus.

Investigation of Iron (III) desferrioxamine B oxalate ligand exchange by UV-vis spectroscopy: the next step in developing a quantitative analytical method for measuring bioavailable iron in marine ecosystems

To date there are no analytical techniques designed to exclusively measure bioavailable iron in marine environments. The goal of this research is to develop such a technique by isolating the bioavailable iron using the terrestrial siderophore desferrioxamine B (DFB). This project contained many challenging aspects, but the specific goal of this study was to develop a robust analytical technique for quantification of Fe(III)-DFB complexes at nanomolar concentrations. Past work showed that oxalate (Ox) promotes photodissociation of Fe(III)-DFB to Fe(Il), and we are specifically interested in the mechanism of this process. A model was developed using known thermodynamic constants for Fe(III)-DFB and Fe(III) oxalato complexes and adjusting for ionic strength. The model was confirmed by monitoring the UV-VIS absorbance of the system at a variety of oxalate concentrations and pH. The model did not include ternary complexes. Next., the rate of Fe(1I) production during UV irradiation was examined. The results showed that the rate of Fe(II) production was based entirely on the [Fe(Ox)?]3- speciation, and that reoxidation of Fe(II) occurred via reactive oxygen intermediates. This reoxidation could be avoided by either decreasing the oxygen concentration or by adding a Fe(II) stabilizing reagent, such as ferrozine. Further studies need to be done to confirm that these results apply at sub nanomolar concentrations, and the issue of Fe(II) reoxidation at lower Fe concentrations needs to be addressed.

Controle de dissociação caótica de moléculas diatômicas

In this work we investigate the dissociation of heteronuclear diatomic molecules subjected to laser pulses. This phenomenon can be modeled by the classical forced Morse oscillator. This system presents a chaotic dynamics associated with the anharmonicity of the internuclear potential and with the coupling of permanent dipole of molecule with the electric field of laser. We want to verify how the dissociation probability evolves while we change the intensity and frequency of laser. We study the phase space of molecules to have a better understanding of system dynamics. We make the calculations changing two parameters of laser (intensity and frequency) and checking how this parameters influences on molecule dissociation. We compare the results of HF molecule (Fluoride acid) and CO molecule (Carbon monoxide) to check how the dipole moment of each molecule can influence on laser interaction

Oxidierte Kohlenwasserstoffe in der oberen Troposphäre: Kinetische Untersuchungen zum Abbau und Einfluss auf HO x und Ozon-Chemie

Der erste Teil dieser Arbeit befasst sich mit der Kinetik der Reaktion des OH-Radikals mit Glykolaldehyd (HOCH2CHO). Die Geschwindigkeitskonstante k1 wurde für diese Reaktion temperaturabhängig bestimmt. Durch gepulste Photolyse wurden OH-Radikale erzeugt. Anschließend wurde die laserinduzierte Fluoreszenz der OH-Radikale bei 309 nm detektiert. Die ermittelte Geschwindigkeitskonstante k1 für die Reaktion von OH mit HOCH2CHO von (8,0 ± 0,8) x 10-12 cm3 Teilchen-1 s-1 erweist sich für den Temperaturbereich von 240 K < T < 362 K als temperaturunabhängig. Zwischen 60 und 250 Torr kann zudem keine Druckabhängigkeit für k1 beobachtet werden. Die unerwartet niedrigere Geschwindigkeitskonstante für die betrachtete Reaktion im Vergleich zur Reaktion von OH mit CH3CHO konnte anhand von Überlegungen zur Korrelation zwischen der C-H-Bindungsstärke und dem H-Abstraktionskanal erklärt werden. Im zweiten Teil dieser Arbeit wurde die Photochemie von Aceton (CH3C(O)CH3), Methylethylketon (C2H5C(O)CH3, MEK) und Acetylbromid (CH3C(O)Br) betrachtet. Für die Photolyse von Aceton (bei 248 nm und 266 nm), MEK (bei 248 nm) und Acetylbromid (bei 248 nm) wurden bei 298 ± 3 K druckabhängig zwischen 5 und 1600 Torr N2 Quantenausbeuten für die Methylbildung (Phi(CH3)) bestimmt. Nach gepulster Photolyse der betrachteten Moleküle wurden die transienten Absorptionssignale der Methylradikale bei 216,4 nm verfolgt. Die Quantenausbeuten wurden relativ zur Photolyse von Methyliodid (CH3I) unter gleichen Reaktionsbedingungen ermittelt. Die erhaltenen Quantenausbeuten für CH3-Radikale nehmen für die beiden Systeme Aceton / 248 nm (Phi(CH3, Aceton) = 1,42 – 0,99) und MEK / 248 nm (Phi(CH3, MEK) = 0,45 – 0,19) druckabhängig zu hohen Drücken ab. Die Druckabhängigkeit von Phi(CH3) wird auf die Konkurrenz zwischen Stoßrelaxation und Dissoziation der schwingungsangeregten Acetylradikale (CH3CO#) zurückgeführt. Für das System Aceton / 266 nm wird keine Druckabhängigkeit von Phi(CH3) = 0,93 ± 0,1 beobachtet. Dies wird damit erklärt, dass CH3CO# nicht genügend Energie besitzt, um die Barriere zur Dissoziation zu überschreiten. Bei der Photolyse von Acetylbromid bei 248 nm wird druckunabhängig Phi(CH3) = 0,92 ± 0,10 bestimmt. In diesem System dissoziieren die schwingungsangeregten Acetylradikale bei allen Drücken vollständig. Bei 266 nm wurde die Gesamtquantenausbeute für die Photodissoziation von Aceton (Phi(diss, 266nm)) bestimmt. Die nach Photolyse erhaltenen Methyl - und Acetylradikale wurden nach Titration mit Br2 durch die Resonanzfluoreszenz der Bromatome detektiert. Phi(diss, 266nm) wurde mit 0,92 ± 0,07 bestimmt.

Isoprene oxidation and its impacts on the atmospheric composition

Terrestrische Vegetation, vor allem tropischer Regenwald, emittiert grosse Mengen flüchtiger organischer Verbindungen (VOCs) in die rnAtmosphäre, die durch Oxidationsreaktionen und Deposition der Reaktionsprodukte wieder entfernt werden. Die Oxidation wird vor allem durch Hydroxyl-Radikale (OH) initiiert, die hauptsächlich durch Photodissoziation von Ozon gebildet werden. Zuvor ging man davon aus, dass biogene VOCs OH in unverschmutzter Luft abbauen und dadurch die atmosphärische Oxidationskapazität verringern. Umgekehrt, führt rndie Oxidation von VOCs in verschmutzter Luft durch die katalytische Wirkung von Stickstoffoxiden (NOx = NO + NO2) zu schädlicher Oxidationsmittelbildung. Flugzeugmessungen atmosphärischer Spurengase, die über dem unberührten Amazonas-Regenwald durchgeführt worden sind, haben jedoch unerwartet hohe OH-Konzentrationen aufgezeigt. Das VOC mit der höchsten Emission in dieser Region war Isopren, dessen Oxidation als stärkeste OH-Senke berechnet wurde. In dieser Arbeit wurde die Hypothese genauestens untersucht, dass die natürliche Isopren-Oxidation in niedrig-NOx Luft OH effizient erneuert. Es wurde ein sehr detaillierter Oxidationsmechanismus für Isopren entwickelt, in dem neueste experimentelle und theoretische Fortschritte umgesetzt worden sind. Die Haupt-OH-Rückgewinnungswege wurden angewendet wodurch gezeigt wurde, dass sie wesentlich zur Oxidation unter niedrig-NOx Bedingungen beitragen. Verstärkte OH-Konzentrationen blieben unter verminderten Lichtverhältnissen, wie sie unter dichten Vegetationsdächern typisch sind, dauerhaft erhalten. Im Vergleich zu Flugzeugmessungen, der neue Oxidationsmechanismus reproduziert die OH-Konzentrationen innerhalb des Unsicherheitsbereiches. Darüber hinaus zeigten Simulationen eine erhebliche Produktion eines Isopren-Dihydroxyepoxids, das ein potenziell wichtiger Vorläufer organischer Aerosole in der Atmosphäre sein könnte. Es wurde einen neuen vereinfachten Oxidationsmechanismus auf Basis des traditionellen Wissenstands entwickelt und seine Anwendung für globale atmosphärische Studien getestet. Die Eingliederung der neuen Oxidationswege in diesen Mechanismus ermöglicht es folgende Auswirkungen der verstärkten VOC-Oxidation zu studieren die Zusammensetzung der Atmosphäre, den Austausch zwischen Erdoberfläche und Atmosphäre, Aerosole und Klima.