The loss of nitrous oxide from the atmosphere is very simple, involving photolysis in the stratospheric UV windows (90%);
In 1997 Yung and Miller presented the "zero point energy" model,^{1} which revised our understanding of the role of photolysis in fractionating isotopomers. In this model the position of the absorption band of an isotopically substituted species is taken to be that of the main species shifted by the change in vibrational zero point energy. The ZPE model incorporates the invariance of the potential energy surfaces to isotopic substitution (Born-Oppenheimer approximation), and concludes that since the excited state is repulsive, increasing the vertical excitation energy by the ZPE change can approximate the isotope effect. Such a small shift of the band position would result in a significantly slower photolysis rate for heavy isotopomers, removing the imbalance in the isotopic budget and the need for exotic in situ sources. On the other hand, subsequent experiments showed isotopic enrichments about twice those predicted by the ZPE model, although qualitative agreement was observed.^{2}
In the present work, the ultraviolet absorption cross sections of six nitrous oxide isotopomers (^{14}N^{14}N^{16}O, ^{14}N^{14}N^{17}O, ^{14}N^{14}N^{18}O, ^{14}N^{15}N^{16}O, ^{15}N^{14}N^{16}O, and ^{15}N^{15}N^{16}O) were computed using the wave packet propagation technique to explore the influence of excited state dynamics, transition dipole surface and initial vibrational state. The wave packet propagation technique used here was the "real wave packet" method .^{3} The wave packet is represented using evenly spaced Fourier grids for R and r coordinates, and a associated Legendre polynominals based discrete variable representation (DVR) for the angle θ in Jacobi coordinates. The absorption potential used in the real wave packet propagation is taken as
where x =R , x_{abs}= 7.0 a_{0}, and C_{abs}= 0.25, and where x =r , x_{abs}= 6.0 a_{0}, and C_{abs}= 0.25. To simulate the dissociation process induced by the photoexcitation, one must determine the initial wave packet correctly. In the time-dependent formulation, the product of the transition dipole moment function with the initial vibrational function of the electronic ground state gives the initial wave function.
where μ(R ,r,θ) is the transition dipole moment function for the corresponding electronic transition, and Ψ_{time-independent} (R ,r,θ) is the initial vibrational wave function of the electronic ground state. The resulting wave function Φ_{init} (R ,r,θ; t = 0) is propagated on the potential energy surface (PES) of the electronically excited state. Therefore, we have calculated the vibrational wave functions of the ground state using the Hamiltonian matrix similar to that of the wave packet propagation; the PES was switched with that of the excited state and the implicitly restarted Lanczos method in ARPACK^{4} was used to diagonalize the matrix. After determining the initial wave packet, it was evolved with time on the PES of the electronic excited state, and finally the autocorrelation function was computed by numerically integrating the product of the initial wave packet and the evolving wave packet at each time step. Total photo-dissociation cross section as a function of the excitation energy was computed as the Fourier transform of the autocorrelation function as a function of time.
where C_{k} is the order-dependent autocorrelation function in the Chebyshev order domain. In the present work, photodissociation cross sections were computed assuming the total angular momentum J = 0 and 1.
Numerous wave packet simulations were carried out and used to calculate the temperature-dependent photodissociation cross sections of the six isotopomers. The photolytic isotopic fractionation constants determined using the calculated cross sections are in good agreement with recent experiments (see Figure 1).