Skip to main content

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home
  • Log in
  • My Cart

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
Physical Sciences

Carbon dioxide photolysis from 150 to 210 nm: Singlet and triplet channel dynamics, UV-spectrum, and isotope effects

Johan A. Schmidt, Matthew S. Johnson, and Reinhard Schinke
  1. aDepartment of Chemistry, University of Copenhagen, DK-2100 Copenhagen Ø, Denmark; and
  2. bMax-Planck-Institut für Dynamik und Selbstorganisation, D-37077 Göttingen, Germany

See allHide authors and affiliations

PNAS October 29, 2013 110 (44) 17691-17696; https://doi.org/10.1073/pnas.1213083110
Johan A. Schmidt
aDepartment of Chemistry, University of Copenhagen, DK-2100 Copenhagen Ø, Denmark; and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: johanalbrechtschmidt@gmail.com msj@kiku.dk
Matthew S. Johnson
aDepartment of Chemistry, University of Copenhagen, DK-2100 Copenhagen Ø, Denmark; and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: johanalbrechtschmidt@gmail.com msj@kiku.dk
Reinhard Schinke
bMax-Planck-Institut für Dynamik und Selbstorganisation, D-37077 Göttingen, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  1. Edited by Mark H. Thiemens, University of California, San Diego, La Jolla, CA, and approved May 23, 2013 (received for review September 21, 2012)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Abstract

We present a first principles study of the carbon dioxide (CO2) photodissociation process in the 150- to 210-nm wavelength range, with emphasis on photolysis below the carbon monoxide + Graphic singlet channel threshold at ∼167 nm. The calculations reproduce experimental absorption cross-sections at a resolution of ∼0.5 nm without scaling the intensity. The observed structure in the 150- to 210-nm range is caused by excitation of bending motion supported by the deep wells at bent geometries in the Graphic and Graphic potential energy surfaces. Predissociation below the singlet channel threshold occurs via spin-orbit coupling to nearby repulsive triplet states. Carbon monoxide vibrational and rotational state distributions in the singlet channel as well as the triplet channel for excitation at 157 nm satisfactorily reproduce experimental data. The cross-sections of individual CO2 isotopologues (12C16O2, 12C17O16O, 12C18O16O, 13C16O2, and 13C18O16O) are calculated, demonstrating that strong isotopic fractionation will occur as a function of wavelength. The calculations provide accurate, detailed insight into CO2 photoabsorption and dissociation dynamics, and greatly extend knowledge of the temperature dependence of the cross-section to cover the range from 0 to 400 K that is useful for calculations of propagation of stellar light in planetary atmospheres. The model is also relevant for the interpretation of laboratory experiments on mass-independent isotopic fractionation. Finally, the model shows that the mass-independent fractionation observed in a series of Hg lamp experiments is not a result of hyperfine interactions making predissociation of 17O containing CO2 more efficient.

  • mass-independent fractionation
  • photodissociation dynamics
  • fine interaction
  • magnetic isotope effect
  • Mars

Carbon dioxide (CO2) is the main component of the atmospheres of Mars, Venus, and the Hadean Earth (1). Its photoabsorption screens solar UV light, determining altitude-dependent photolysis rates, and its concentrations and infrared absorptions make it a powerful greenhouse gas. Graphic photodissociation is the basis of these atmospheres’ photochemistry and is the primary source of carbon monoxide (CO) and O2. Although the initially high concentration of CO2 during Earth’s Hadean era decreased as carbonate rocks accumulated, CO2 continued to be a prominent atmospheric gas, enhancing surface temperature and attenuating UV light, with a partial pressure Graphic mbar in the Archean (2) and Graphic mbar in the Proterozoic (3). Its influence on the radiative properties and chemical composition of Earth’s atmosphere has continued to the present day; one example is that CO2 photolysis is the main source of mesospheric CO (4). Variations in the abundances of naturally occurring stable isotopes, including reaction mechanisms exhibiting mass-independent fractionation, are central to efforts to interpret environmental records ranging from sedimentary rocks to oceanic carbonates to glacial ice (5). Thus, Graphic photolysis is both directly and indirectly linked to isotopic variations found in the environment.

The UV absorption band of CO2 Graphic consists of two broad, strongly overlapping bands peaking around 145 nm and 133 nm, respectively (6). Both exhibit pronounced vibrational structures; those overlaid on the first band are complex, whereas those on the second are more regular. The triplet, Graphic, and singlet, Graphic, dissociation channels open at 227.5 nm and 167.2 nm, respectively. Because of the very small actinic flux below Graphic nm, atmospheric photolysis of CO2 occurs almost exclusively via the triplet channel except at very high altitudes.

The electronic structure of CO2 has been the subject of several theoretical studies (7⇓–9). The first quantum dynamics analysis of the UV photoabsorption of Graphic was provided only recently by Grebenshchikov (10). This work determined 3D potential energy surfaces (PESs) of six singlet states and calculated the absorption cross-section starting near the singlet channel threshold. In contrast, in the present study, we focus on the wavelengths longer than 167 nm, which are important for most atmospheric applications.

CO2 is a linear molecule with 16 valence electrons, properties common to two other key terrestrial trace gases: nitrous oxide (N2O) and carbonyl sulfide (OCS)(6). We have recently completed first principle computations of the UV absorption spectra of these two molecules using accurate PESs, the transition dipole moment (TDM) functions coupling them with the ground state, and the quantum mechanical wave packet methodology (11⇓⇓⇓⇓–16). Our calculations do an excellent job of reproducing available experimental data, including the OCS isotopologue absorption spectra (17, 18) and the temperature- and isotopologue-dependent N2O cross-sections (14, 19).

The wave packet methodology extracts a wealth of information from the potential energy and transition dipole surfaces, including vibrational frequencies; temperature- and isotopologue-dependent absorption spectra; and detailed descriptions of the photodissociation dynamics, including product quantum state, and velocity and angular distributions (20). In practice, the method is limited mainly by the quality of the surfaces. Using the method we have built and tested on N2O and OCS, we are able to present a similar analysis of the CO2 UV absorption and dissociation from the onset of the UV spectrum (Graphic nm) to ca. 150 nm. Our calculations reproduce a variety of existing experimental data; therefore, we can report reliable absorption cross-sections for a set of isotopologues and for a wide range of temperatures by changing the isotopic masses and/or the populations of vibrational states.

Low-Energy Absorption and Dissociation Mechanisms

Six excited electronic states are energetically accessible below Graphic eV (Graphic nm) in the Franck–Condon (FC) region (Fig. 1): the three singlet states Graphic (A), Graphic (B), and Graphic (C) and the three triplet states Graphic (a), Graphic (b), and Graphic (c) (letters in parentheses indicate abbreviations used below). Fig. 1A, which shows 1D cuts along one of the two C-O bonds, Graphic, provides an overview of these states and how they connect with the dissociation channels (Table S1); the potential of the ground electronic state Graphic (X) is also included. Fig. 2B shows potential cuts along the O-C-O bending angle α. The potential energies were calculated from first principles using the multiconfiguration reference internally contracted configuration interaction (MRCI) theory (21, 22) including the Davidson correction (MRCI+Q) (SI Text). The corresponding potential cuts for OCS and Graphic are qualitatively very similar (12, 13).

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

1D potential energy curves for the lowest four singlet (full curves) and lowest three triplet (broken curves) electronic states as a function of Graphic, one of the two C-O bond lengths (A), and as a function of the O-C-O bending angle α (B). The Graphic states (X, A, and a) are shown in black, and the Graphic states (B, b, C, and c) are shown in red. In this figure and throughout this article, energy is measured relative to the global minimum of the X state. deg., degree.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

2D representations of the A state PES as a function of the two C-O bond lengths Graphic and Graphic (A) and as a function of one C-O bond length and the O-C-O bond angle α (B). The red contour is at 7 eV, and the green and blue lines mark the seam of intersections with the triplet states b Graphic and c Graphic, respectively. The spacing between the contours is 0.5 eV.

The three singlet states A, B, and C have similar energies near linearity Graphic, and they are coupled to one another (10, 23). States A and C form a Renner–Teller pair and are coupled by interaction between nuclear and electronic rotation about the O-C-O axis. The two Graphic states B and C are coupled by nonadiabatic coupling matrix elements (NACMEs) (20). A rigorous theoretical treatment, including the three singlet states (not to mention the Graphic triplet states), is therefore very challenging. The Renner–Teller and nonadiabatic couplings are largest for linear O-C-O. For total energies below ca. 8 eV, however, absorption takes place primarily at angles smaller than Graphic (Fig. 1B). Therefore, the region close to linearity is not important, and a model in which the couplings between A, B, and C are ignored and the absorption cross-sections are calculated separately for each state is justified. In our previous studies of photodissociation in Graphic and OCS, which, near linearity, have electronic structures identical to that of Graphic, we used the same model and were able to reproduce the absorption cross-sections very well both at low energies and throughout the first UV band (12, 13).

We have calculated PESs for the ground state X, the three excited singlet states A–C, and the three triplet states a–c using the MRCI+Q method on large 3D grids (SI Text, Figs. S1–S3, and Table S1). The calculated energies for the three fundamental vibrations in Graphic deviate from the experimental data by 1% (Table S2), which gives an indication of the quality of the potentials. Our earlier studies of Graphic (13) and OCS (11, 12) found equally good agreement between experimental and theoretical vibrational excitation energies. As an example of an excited state PES, Fig. 2 shows 2D representations of the A state PES, Graphic, which is the most important one; the corresponding contour plots for the B state are very similar. Graphic has a deep well at bent geometries, whose equilibrium energy is 5.532 eV [i.e., 1.961 eV below the (classical) singlet dissociation channel]. The well leads to substantial trapping of the O and CO products during dissociation, even for energies above the singlet channel threshold, and therefore has a strong effect on the absorption cross-section and on the product state distributions. Although the A and B PESs of OCS and Graphic have similar potential wells, they are not as deep and produce only minor structures in the respective absorption spectra (12, 13).

The TDMs of A, B, and C with the ground electronic state (e.g., Graphic) have been calculated at the MRCI level on small coordinate grids around the FC region (SI Text). The transition to the A state is dominant because Graphic is the largest (Fig. S4). It lies in the plane of the molecule and has two components, Graphic and Graphic. The B state TDM, Graphic, has a single component perpendicular to the molecular plane. The Graphic increases rapidly as the molecule bends, whereas Graphic and Graphic increase along the asymmetrical stretch coordinate and have similar values. The C state TDM is very small compared with Graphic and Graphic; therefore, excitation of the C state has been neglected in the calculations.

There is an additional complication. The X and A states belong to different irreducible representations in the Graphic configuration; therefore, their PESs are allowed to cross along the Graphic symmetry line [conical intersection (CI)]. This happens around Graphic, as seen in Fig. 1B and Fig. S1. When Graphic, the two states belong to the same irreducible representation Graphic and the potentials form an avoided crossing. In the vicinity of the CI, states A and X are coupled through a NACME and the manifolds of vibrational states belonging to either A or X are mixed. The low-energy part of the measured absorption spectrum reflects this mixing in the form of irregular subnanometer structures.

Below the singlet channel threshold (Graphic eV), dissociation from X, A, or B is only possible via coupling to the repulsive triplet state a, b, or c. Spin-orbit (SO) coupling splits each (zero-order) triplet state into three components with Graphic. The SO coupling elements between the singlet and triplet states have been calculated as described for Graphic (16) (SI Text), and the elements involving A and B are shown in Figs. S5 and S6. Most important is the coupling between A on one hand and b0 and c0 on the other; the corresponding SO elements exceed 0.01 eV in the vicinity of the intersection seams, which are indicated in Fig. 2B.

Absorption Cross-Section

Temperature-dependent absorption cross-sections were calculated for states A and B separately, using wave packet propagation (Materials and Methods). The total cross-section is the sum of Graphic and Graphic. The cross-section for state C was neglected in the analysis because of the very small TDM. Fig. 3 shows the calculated absorption cross-section, shifted by 100 cm−1 to lower energies, at different temperatures in comparison to experimental data. Similar shifts were applied in our studies of Graphic (13) and OCS (12), and are necessary due to small inaccuracies in the calculated electronic excitation energies. We emphasize that the intensity of the calculated cross-section was not scaled, as was done for other molecules [e.g., OCS (12) and Graphic (13)]. Fig. 3A shows an overview from 140 to 170 nm. The underestimation below 150 nm is due to the neglect of higher electronic states [i.e., the more intense 133-nm band (10)]. Above 150 nm, the agreement with the measured cross-sections is very good even at the longest wavelengths, where the cross-section has fallen by about six orders of magnitude. The temperature dependence is also reproduced. Karaiskou et al.’s data (26), shown in Fig. 3D, have orders of magnitude higher resolution, so they can only be used to show trends.

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Comparison of calculated (black lines) and measured (red and blue lines and symbols) absorption cross-sections in different wavelength regions and for several temperatures. The experimental data in A–C are from Yoshino and colleagues (24, 25) (www.cfa.harvard.edu/amp/ampdata/co296/co2.html), and the single-wavelength results in D are from Karaiskou et al. (26). T, temperature.

At very long wavelengths, excitation of the B state is dominant, whereas at shorter wavelengths Graphic predominates. The ratio Graphic is about 9 at 150 nm; the crossover occurs at around 200 nm. Excitation of vibrationally excited states of X plays a significant role in the studied region. At 170 K, for example, the Graphic ratio is Graphic at 166 nm and Graphic at 180 nm; at 300 K, these ratios increase to Graphic and Graphic, respectively. This behavior is driven by the slope of the transition dipole surfaces (Fig. S4).

The theoretical and experimental cross-sections show undulations whose positions and amplitudes are largely reproduced by the calculation. The spacing is about 650 cm−1 at 200 nm and 520 cm−1 at 150 nm. The undulations correspond to a clear recurrence of the (0,0,0) autocorrelation function with a period of ≈60 fs. As for OCS (12), the structure mainly arises from bending motion in the deep wells of Graphic and Graphic at bent geometries, above and below the singlet channel threshold.

Superimposed on the undulations, the calculated spectrum, and, to a lesser extent, the measured spectrum exhibit additional structure (Fig. 3 B–D) that arises from specific highly excited stretching/bending states of A and B. Moreover, the experimental spectra show details finer than the calculated spectra (Fig. 3C). These structures probably represent vibrational states in the dense quasicontinuum of the X state, which is coupled to the A state as described above. Because the X state is not included in the model, the calculated spectrum does not show these fine irregular spectral features.

Full reproduction of the experimental spectrum in the long-wavelength region would require prohibitively high accuracy for all PESs and coupling surfaces. In this work, we focus on reproducing the low-resolution undulations, the general increase with photon energy, and the magnitude of the absorption cross-section. This allows us to calculate cross-sections and fractionation spectra for all isotopologues and for temperatures from 0 to 400 K, with sufficient accuracy to derive reliable loss rates and isotopic fractionation for broadband photolysis in planetary atmospheres (SI Text, Datasets S1–S6). However, because of the resolution limit, the results cannot be directly applied to the question of self-shielding.

Product State Distributions for Excitation at 157 nm

A separate model was constructed explicitly considering the SO coupling of the A state to the b0 and c0 triplet states (Materials and Methods). This model was used to calculate the Graphic branching ratio, as well as the rotational (j) and vibrational (v) distributions of product CO in the singlet and triplet channels for photodissociation at 157 nm. The calculated Graphic yield of 12% agrees reasonably well with the experimental estimate of Graphic (27). The singlet channel rotational distribution (Fig. S7) for Graphic shows irregular oscillations from the onset at Graphic up to Graphic, where the distribution is essentially zero in good agreement with experiments (28). The triplet channel rotational distributions (Fig. S7) are qualitatively similar but extend to much higher values of j. The vibrational distribution of the singlet channel is narrow and peaks at Graphic. The Graphic ratio is 4.15, in excellent agreement with the measured value of Graphic (28). The triplet vibrational distribution also peaks at Graphic but falls off much more slowly with v, extending to Graphic. The total kinetic energy distribution constructed from the vibrational-rotational distribution agrees qualitatively with the measurement (29). A more detailed discussion will be published elsewhere.

Isotope Effects

Wavelength-dependent isotopic fractionation constants Graphic (Materials and Methods) for three isotopologues are shown in Fig. 4 B and C. In Fig. 4B, the ε constants are averaged over 0.25 nm, which is sufficiently narrow to show the range of fractionation generated by the undulations. Isotopic substitution alters the vibrational wave functions and energy levels of both the ground and excited states. The isotope shifting brings the undulations in the cross-sections for the heavier isotopologues out of phase with the undulations of the light 12C16O2 (Fig. 4A), leading to the periodic oscillations. They vary rapidly with wavelength and can range from Graphic to 500‰ over an interval of just a few nanometers.

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

(A) Calculated cross-sections for 12C16O2 and 13C16O2. (B and C) Calculated fractionation constants Graphic as a function of wavelength for different isotopologues of CO2. Results were averaged over a Gaussian window with FWHMs of 0.25 nm (B) and 2.5 nm (C). The broken red line in C shows Graphic calculated using the ZPE method (30). per mil, parts per thousand.

Fig. 4C shows the fractionation constants averaged over a 2.5-nm window, chosen to illustrate broadband fractionation effects. Photolytic fractionation is increasingly negative with increasing wavelength. At very low photon energies (long wavelengths), absorption takes place mainly into the lower vibrational states, which are localized around the bent minimum of the A and B state PESs. Because the X-state wave functions, localized near linearity, are narrower for the heavier isotopologues, the overlap with the bent vibrational states of A and B becomes smaller with isotopic substitution. This general effect is amplified by the fact that the TDM decreases with the bending angle (i.e., when Graphic becomes more linear).

We have included a zero point energy (ZPE) method-based prediction of the 12C18O16O fractionation in Fig. 4C. The magnitude of the fractionation is approximately threefold smaller for the ZPE method compared with the wave packet results (Fig. 4C). A similar underprediction by the ZPE method, arising from the slope of the transition dipole surface and changes of the vibrational wave functions with isotopic substitution, is observed for N2O (14) and OCS (17, 18). The ZPE method is a semiempirical approach, which uses the cross-section for one isotopologue (typically the most abundant) and the ZPEs of the different isotopologues to estimate the cross sections of the remaining (rare) isotopologues. The ZPE method is convenient because the required data (cross-sections and ZPEs) are often readily available. Unfortunately, the method has numerous shortcomings; for example, it does not consider changes to the width of the cross-sections, changes in the vibrational energy levels of the excited state, or changes in the overall magnitude of the cross-sections due to the TDM. The ZPE method was introduced by Yung and Miller (30) to investigate the isotopic fractionation in N2O photolysis. The ZPE method was later applied to O-isotope fractionation in CO2 photolysis (31). Our results demonstrate that the ZPE method is not well suited to describe fractionation in CO2 photolysis.

Fig. 5 shows results from a simulated photolysis experiment using a Gaussian lamp spectral function with the center varied from 160 to 190 nm and a width of 2.5 nm (Materials and Methods). The simulation continues to 50% photolysis of the initial 12C16O2. There is a significant increase in Graphic and Graphic of the remaining CO2 with λ. The value of Graphic has both positive and negative excursions from zero. The average value of Graphic from 160 to 190 nm is 0.8‰, which is small considering the large value of Graphic (58.8‰) and the large extent of reaction, showing that the general trend in fractionation of oxygen isotopes is close to being mass-dependent. The Graphic value shows larger departures from zero. The sign of Graphic depends intricately on λ; the average value in the plotted interval is Graphic. Because actinic flux decreases sharply by ca. Graphic between 190 and 180 nm, there is a distinct possibility of producing mass-independent fractionation and a clumped isotope anomaly when CO2 is photolyzed in planetary atmospheres. These implications will be explored in detail in a forthcoming publication.

Fig. 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 5.

Simulated Graphic, Graphic, Graphic, and Graphic of the remaining CO2 for photolysis using a Gaussian lamp spectrum with a FWHM of 2.5 nm and a moving center at λ. The full black line shows an example of the lamp spectrum at Graphic nm. The photolytic yield was set to 50%.

In a series of experiments, Bhattacharya and colleagues (32⇓–34) measured the isotopic fractionation caused by CO2 photolysis. The experiments used a mercury lamp with an emission at 184.9 nm. The widths of the emission were cited to be 10 nm (32) and 2 nm (33). This light is below the singlet channel dissociation threshold at 167 nm, and dissociation can only proceed via coupling to the triplet states. The experiments (32, 33) found that product CO and O2 are significantly enriched in 17O, although being neither enriched nor depleted in 18O compared with the initial CO2. This observed mass-independent fractionation was interpreted as evidence of a predissociative mechanism dominated by hyperfine interaction, particularly coupling between singlet and triplet electronic states dependent on both the electronic and nuclear spin (16O and 18O are spin 0, whereas 17O is spin Graphic). The hypothesis is that that CO2 containing one or more 17O nuclei predissociates more efficiently (i.e., faster) than other CO2 isotopologues, which therefore have a higher probability of being quenched radiatively or through collisions. Our results contradict this interpretation because it would require SO-induced predissociation onto the triplet surfaces to be inefficient; otherwise, SO coupling (fine interaction) would completely overshadow any predissociated induced effects sensitive to the nuclear spin (hyperfine interaction). However, in the case of CO2, the fine interactions are strong (e.g., Figs. S5 and S6) and give rise to fast predissociation on a time scale of several picoseconds. Our assertation, that hyperfine interactions are not causing the observed fractionation pattern, is further supported by the observed photolysis quantum yield of unity at 184.9 nm (35) for CO2, which leaves no room for a preferential predissociation of 17O-substituted CO2. The experimental results (32⇓–34) can be explained if the width of the Hg emission is narrower than indicated, with structure on a subnanometer scale. The experimental observations may arise from resonance between the narrowly structured emission lines of the lamp and the isotopologue-dependent fine structure of the Graphic cross-section. This would also explain the very large variation in the isotopic fractionation with temperature (34), because relatively small differences in temperature lead to significant differences in the population of the rotational states, which can have a large effect on the fine structure of the cross-sections. High-resolution measurements of the isotopologue-specific cross-sections of CO2 similar to studies performed for SO2 (36, 37) are needed to fully unravel isotopic fractionation fully on a subnanometer scale.

Materials and Methods

Potential energy energy surfaces for the Graphic, Graphic, Graphic, and Graphic electronic states were constructed using the MRCI+Q method (21, 22) and Dunning’s aug-cc-pVQZ (38) orbital basis set. Further details are given in SI Text.

Absorption cross-sections for different isotopologues of CO2 were calculated using time-dependent quantum mechanical methodology (20, 39). Wave packets were propagated on the uncoupled A and B excited state PESs, and the cross-sections were obtained by Fourier-transforming the autocorrelation functions. The calculations were performed in terms of the Jacobi coordinates R (distance from one O atom to the center-of-mass of the product CO), r (bond length of the product CO molecule), and γ (angle between R and r). The initial wave packets were defined as the product of a vibrational state of the X state and the modulus of the A or B state TDM. The wave packets were propagated in time using a Chebychev polynomial expansion of the time evolution operator (40). The action of the Hamiltonian on the wave packet was evaluated using the pseudospectral scheme of Le Quéré and Leforestier (41). The Fourier method (42) was used to evaluate the action of the radial kinetic energy operators, whereas the action of the angular kinetic energy operator was evaluated by transforming between a grid representation and a basis set representation consisting of associated Legendre polynomials, Graphic, with Graphic depending on the initial vibrational state. The lowest nine vibrational states were considered for the calculation of thermal cross-sections. The dynamics calculations were carried out using the Wavepacket program package (43).

All quantum dynamics calculations were performed with a total angular momentum Graphic. This is an approximation and does not satisfy the selection rules for dipole excitation. A rigorous treatment with Graphic is a formidable computational task, which, if the spectrum does not show rotational resolution, merely leads to some, usually less interesting, broadening. All our calculations for N2O and OCS used the Graphic approximation, and very good agreement with the measured spectra was achieved. In the present study, rotational effects were approximately accounted for by level-shifting the cross-sections using a linear rigid rotor-guided approach and the selection rule Graphic. This effectively led to a temperature-dependent broadening of the narrow lines. The bands are also broadened by coupling of all singlet states to the repulsive triplet states; however, this effect is overshadowed by the rotational broadening. Nonadiabatic coupling between the A and X states causes mixing of the vibrational states of A with the quasicontinuum of X, which leads to additional strong broadening. This effect was empirically accounted for by averaging the A state cross-sections with a Gaussian with a FWHM of 250 cm−1. Calculated absorption cross-sections for five isotopologues and temperatures from 120 to 395 K are available in SI Text (Datasets S1–S6).

When calculating the product state distributions and the singlet/triplet branching ratio, the SO coupling of the A state to the b0 and c0 triplet states was explicitly taken into account (i.e., three coupled wave packets were propagated simultaneously). These calculations were performed on a larger grid extending further along R than used in the calculations for the spectra. The methodology used was similar to that described by Schinke et al. (16) for N2O.

The wavelength-dependent fractionation constants were defined as Graphic, where σ, Graphic, Graphic, and Graphic are the cross-sections for 12C16O2, 12C17O16O, 12C18O16O, and 13C16O2, respectively. The simulation in Fig. 4C was carried out as follows: The rate of photolysis for all relevant isotopologues was calculated for different values of λ (the center of the lamp spectrum) asEmbedded Image

The FWHM of the lamp spectrum was set to 2.5 nm. The photolysis yield was 50% [i.e., Graphic]. The initial isotopic composition of COGraphic was assumed to be stochastic, and the oxygen and carbon isotope ratios were set to those of Vienna Standard Mean Ocean Water and Vienna Pee Dee Belemnite (i.e., Graphic, Graphic, and Graphic). The δ and Graphic values were calculated asEmbedded Image

etc.,Embedded Image

andEmbedded Image

where Graphic is the mass 47-to-mass 44 ratio of CO2 and Graphic is the same ratio if the O and C isotopes are distributed stochastically (44). Only the dominant isotopologues, Graphic, Graphic, Graphic, and Graphic, were considered when calculating Graphic.

Acknowledgments

We thank the Initial Training Network in Mass Independent Fractionation sponsored by the European Community’s Seventh Framework Program (FP7/2007-2013) under Grant Agreement 237890 and the Gesellschaft für wissenschaftliche Datenverarbeitung mbH Göttingen.

Footnotes

  • ↵1To whom correspondence may be addressed. E-mail: johanalbrechtschmidt{at}gmail.com or msj{at}kiku.dk.
  • Author contributions: J.A.S., M.S.J., and R.S. designed research; J.A.S. performed research; J.A.S. and R.S. contributed new reagents/analytic tools; J.A.S., M.S.J., and R.S. analyzed data; and J.A.S., M.S.J., and R.S. wrote the paper.

  • The authors declare no conflict of interest.

  • This article is a PNAS Direct Submission.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1213083110/-/DCSupplemental.

References

  1. ↵
    1. Yung YL,
    2. DeMore WB
    (1999) Photochemistry of Planetary Atmospheres (Oxford Univ Press, New York).
  2. ↵
    1. Driese SG,
    2. et al.
    (2011) Neoarchean paleoweathering of tonalite and metabasalt: Implications for reconstructions of 2.69 Ga early terrestrial ecosystems and paleoatmospheric chemistry. Precambrian Res 189(1–2):1–17.
    OpenUrlCrossRef
  3. ↵
    1. Bao H,
    2. Lyons JR,
    3. Zhou C
    (2008) Triple oxygen isotope evidence for elevated CO2 levels after a Neoproterozoic glaciation. Nature 453(7194):504–506.
    OpenUrlCrossRefPubMed
  4. ↵
    1. Clerbaux C,
    2. et al.
    (2005) Carbon monoxide distribution from the ACE-FTS solar occultation measurements. Geophys Res Lett 32(16):L16S01.
    OpenUrlCrossRef
  5. ↵
    1. Thiemens MH
    (1999) Mass-independent isotope effects in planetary atmospheres and the early solar system. Science 283(5400):341–345.
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. McGlynn SP,
    2. Rabalais JW,
    3. McDonald JR,
    4. Scherr VM
    (1971) Electronic spectroscopy of isoelectronic molecules. II. Linear triatomic groupings containing sixteen valence electrons. Chem Rev 71(1):73–108.
    OpenUrlCrossRef
  7. ↵
    1. Knowles PJ,
    2. Rosmus P,
    3. Werner H-J
    (1988) On the assignment of the electronically excited singlet states in linear CO2. Chem Phys Lett 146(3–4):230–235.
    OpenUrlCrossRef
  8. ↵
    1. Spielfiedel A,
    2. et al.
    (1992) Bent valence excited states of CO2. J Chem Phys 97(11):8382–8388.
    OpenUrlCrossRef
  9. ↵
    1. Spielfiedel A,
    2. Feautrier N,
    3. Chambaud G,
    4. Rosmus P,
    5. Werner H-J
    (1993) The first dipole-allowed electronic transition Graphic of CO2. Chem Phys Lett 216(1–2):162–166.
    OpenUrlCrossRef
  10. ↵
    1. Grebenshchikov SY
    (2012) Communication: Multistate quantum dynamics of photodissociation of carbon dioxide between 120 nm and 160 nm. J Chem Phys 137(2):021101.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Schmidt JA,
    2. Johnson MS,
    3. McBane GC,
    4. Schinke R
    (2012) Communication: Multi-state analysis of the OCS ultraviolet absorption including vibrational structure. J Chem Phys 136(13):131101.
    OpenUrlCrossRefPubMed
  12. ↵
    1. Schmidt JA,
    2. Johnson MS,
    3. McBane GC,
    4. Schinke R
    (2012) The ultraviolet spectrum of OCS from first principles: Electronic transitions, vibrational structure and temperature dependence. J Chem Phys 137(5):054313.
    OpenUrlCrossRefPubMed
  13. ↵
    1. Schinke R
    (2011) Photodissociation of N2O: Potential energy surfaces and absorption spectrum. J Chem Phys 134(6):064313.
    OpenUrlCrossRefPubMed
  14. ↵
    1. Schmidt JA,
    2. Johnson MS,
    3. Schinke R
    (2011) Isotope effects in N2O photolysis from first principles. Atmos Chem Phys 11:8965–8975.
    OpenUrlCrossRef
  15. ↵
    1. Schmidt JA,
    2. Johnson MS,
    3. Lorenz U,
    4. McBane GC,
    5. Schinke R
    (2011) Photodissociation of N2O: Energy partitioning. J Chem Phys 135(2):024311.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Schinke R,
    2. Schmidt JA,
    3. Johnson MS
    (2011) Photodissociation of N2O: Triplet states and triplet channel. J Chem Phys 135(19):194303.
    OpenUrlCrossRefPubMed
  17. ↵
    1. Hattori S,
    2. et al.
    (2011) Ultraviolet absorption cross sections of carbonyl sulfide isotopologues OC32S, OC33S, OC34S, and O13CS: Isotopic fractionation in photolysis and atmospheric implications. Atmos Chem Phys 11:10293–10303.
    OpenUrlCrossRef
  18. ↵
    1. Schmidt JA,
    2. et al.
    (2013) OCS photolytic isotope effects from first principles: Sulfur and carbon isotopes, temperature dependence and implications for the stratosphere. Atmos Chem Phys 13:511–1520.
    OpenUrl
  19. ↵
    1. von Hessberg P,
    2. et al.
    (2004) Ultra-violet absorption cross sections of isotopically substituted nitrous oxide species: 14N14NO, 15N14NO, 14N15NO and 15N15NO. Atmos Chem Phys 4:1237–1253.
    OpenUrlCrossRef
  20. ↵
    1. Schinke R
    (1993) Photodissociation Dynamics, Cambridge Monographs on Atomic, Molecular, and Chemical Physics (Cambridge Univ Press, Cambridge, UK).
  21. ↵
    1. Werner H-J,
    2. Knowles PJ
    (1988) An efficient internally contracted multiconfiguration reference configuration interaction method. J Chem Phys 89(9):5803–5814.
    OpenUrlCrossRef
  22. ↵
    1. Knowles PJ,
    2. Werner H-J
    (1988) An efficient method for the evaluation of coupling coefficients in configuration interaction calculations. Chem Phys Lett 145(6):514–522.
    OpenUrlCrossRef
  23. ↵
    1. Grebenshchikov SY,
    2. Raffaele B
    (2012) Crossing electronic states in the Franck Condon zone of carbon dioxide: A five-fold closed seam of conical and glancing intersections. J Phys Chem Lett 3(21):3223–3227.
    OpenUrlCrossRef
  24. ↵
    1. Yoshino K,
    2. et al.
    (1996) Absorption cross section measurements of carbon dioxide in the wavelength region 118.7-175.5 nm and the temperature dependence. J Quant Spectrosc Radiat Transf 55(1):53–60.
    OpenUrlCrossRef
  25. ↵
    1. Parkinson WH,
    2. Rufus J,
    3. Yoshino K
    (2003) Absolute absorption cross section measurements of CO2 in the wavelength region 163-200 nm and the temperature dependence. Chem Phys 290(2–3):251–256.
    OpenUrlCrossRef
  26. ↵
    1. Karaiskou A,
    2. Vallance C,
    3. Papadakis V,
    4. Vardavas IM,
    5. Rakitzis TP
    (2004) Absolute absorption cross-section measurements of CO2 in the ultraviolet from 200 to 206 nm at 295 and 373 K. Chem Phys Lett 400(1–3):30–34.
    OpenUrlCrossRef
  27. ↵
    1. Stolow A,
    2. Lee YT
    (1993) Photodissociation dynamics of CO2 at 157.6 nm by photofragment-translational spectroscopy. J Chem Phys 98(3):2066–2076.
    OpenUrlCrossRef
  28. ↵
    1. Miller RL,
    2. Kable SH,
    3. Houston PL,
    4. Burak I
    (1992) Product distributions in the 157 nm photodissociation of CO2. J Chem Phys 96(1):332–338.
    OpenUrlCrossRef
  29. ↵
    1. Chen Z,
    2. Liu F,
    3. Jiang B,
    4. Yang X,
    5. Parker DH
    (2010) Imaging CO2 photodissociation at 157 nm: State-to-state correlations between Graphic and GraphicJ Phys Chem Lett 1(12):1861–1865.
    OpenUrlCrossRef
  30. ↵
    1. Yung YL,
    2. Miller CE
    (1997) Isotopic fractionation of stratospheric nitrous oxide. Science 278(5344):1778–1780.
    OpenUrlAbstract/FREE Full Text
  31. ↵
    1. Liang M-C,
    2. Blake GA,
    3. Yung YL
    (2008) Seasonal cycle of C16O16O, C16O17O, and C16O18O in the middle atmosphere: Implications for mesospheric dynamics and biogeochemical sources and sinks of CO2. J Geophys Res 113:D12305.
    OpenUrlCrossRef
  32. ↵
    1. Bhattacharya SK,
    2. Savarino J,
    3. Thiemens MH
    (2000) A new class of oxygen isotopic fractionation in photodissociation of carbon dioxide: Potential implications for atmospheres of Mars and Earth. Geophys Res Lett 27(10):1459–1462.
    OpenUrlCrossRef
  33. ↵
    1. Mahata S,
    2. Bhattacharya SK
    (2009) Anomalous enrichment of 17O and 13C in photodissociation products of CO2: Possible role of nuclear spin. J Chem Phys 130(23):234312.
    OpenUrlCrossRefPubMed
  34. ↵
    1. Mahata S,
    2. Bhattacharya SK
    (2009) Temperature dependence of isotopic fractionation in CO2 photolysis. Chem Phys Lett 477(1–3):52–56.
    OpenUrlCrossRef
  35. ↵
    1. DeMore WB,
    2. Mosesman M
    (1971) Photolysis of CO2 at 1849 Å. Journal of the Atmospheric Sciences 28(6):842–846.
    OpenUrlCrossRef
  36. ↵
    1. Danielache SO,
    2. et al.
    (2012) Photoabsorption cross-section measurements of 32S, 33S, 34S and 36S sulfur dioxide for the Graphic absorption band. J Geophys Res 117:D24301.
    OpenUrl
  37. ↵
    1. Hattori S,
    2. et al.
    (2013) SO2 photoexcitation mechanism links mass-independent sulfur isotopic fractionation in cryospheric sulfate to climate impacting volcanism. Proc Natl Acad Sci USA 110:17656–17661.
    OpenUrlAbstract/FREE Full Text
  38. ↵
    1. Dunning TH Jr.
    (1989) Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen. J Chem Phys 90(2):1007–1023.
    OpenUrlCrossRef
  39. ↵
    1. Balint-Kurti GG
    (2004) Wavepacket theory of photodissociation and reactive scattering. Advances in Chemical Physics 128:249–301.
    OpenUrl
  40. ↵
    1. Tal-Ezer H,
    2. Kosloff R
    (1984) An accurate and efficient scheme for propagating the time dependent Schrödinger equation. J Chem Phys 81(9):3967–3971.
    OpenUrlCrossRef
  41. ↵
    1. Le Quéré F,
    2. Leforestier C
    (1990) Quantum exact three-dimensional study of the photodissociation of the ozone molecule. J Chem Phys 92(1):247.
    OpenUrlCrossRef
  42. ↵
    1. Kosloff D,
    2. Kosloff R
    (1983) A Fourier method solution for the time dependent Schrödinger equation as a tool in molecular dynamics. J Comput Phys 52(1):35.
    OpenUrlCrossRef
  43. ↵
    Schmidt B, Lorenz U (2009) Wavepacket 4.6: A program package for quantum-mechanical wave packet propagation and time-dependent spectroscopy. Available at http://wavepacket.sourceforge.net. Accessed June 6, 2013.
  44. ↵
    1. Eiler JM
    (2007) Clumped-isotope geochemistry—The study of naturally occurring, multiply substituted isotopologues. Earth Planet Sci Lett 262(3–4):309–327.
    OpenUrlCrossRef
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Carbon dioxide photolysis from 150 to 210 nm: Singlet and triplet channel dynamics, UV-spectrum, and isotope effects
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
CO2 UV spectrum from first principles
Johan A. Schmidt, Matthew S. Johnson, Reinhard Schinke
Proceedings of the National Academy of Sciences Oct 2013, 110 (44) 17691-17696; DOI: 10.1073/pnas.1213083110

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
CO2 UV spectrum from first principles
Johan A. Schmidt, Matthew S. Johnson, Reinhard Schinke
Proceedings of the National Academy of Sciences Oct 2013, 110 (44) 17691-17696; DOI: 10.1073/pnas.1213083110
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley

Article Classifications

  • Physical Sciences
  • Earth, Atmospheric, and Planetary Sciences
Proceedings of the National Academy of Sciences: 110 (44)
Table of Contents

Submit

Sign up for Article Alerts

Jump to section

  • Article
    • Abstract
    • Low-Energy Absorption and Dissociation Mechanisms
    • Absorption Cross-Section
    • Product State Distributions for Excitation at 157 nm
    • Isotope Effects
    • Materials and Methods
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Water from a faucet fills a glass.
News Feature: How “forever chemicals” might impair the immune system
Researchers are exploring whether these ubiquitous fluorinated molecules might worsen infections or hamper vaccine effectiveness.
Image credit: Shutterstock/Dmitry Naumov.
Reflection of clouds in the still waters of Mono Lake in California.
Inner Workings: Making headway with the mysteries of life’s origins
Recent experiments and simulations are starting to answer some fundamental questions about how life came to be.
Image credit: Shutterstock/Radoslaw Lecyk.
Cave in coastal Kenya with tree growing in the middle.
Journal Club: Small, sharp blades mark shift from Middle to Later Stone Age in coastal Kenya
Archaeologists have long tried to define the transition between the two time periods.
Image credit: Ceri Shipton.
Illustration of groups of people chatting
Exploring the length of human conversations
Adam Mastroianni and Daniel Gilbert explore why conversations almost never end when people want them to.
Listen
Past PodcastsSubscribe
Panda bear hanging in a tree
How horse manure helps giant pandas tolerate cold
A study finds that giant pandas roll in horse manure to increase their cold tolerance.
Image credit: Fuwen Wei.

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Special Feature Articles – Most Recent
  • List of Issues

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Subscribers
  • Librarians
  • Press
  • Cozzarelli Prize
  • Site Map
  • PNAS Updates
  • FAQs
  • Accessibility Statement
  • Rights & Permissions
  • About
  • Contact

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490