The Effect of ß-Hydrogens on the Tropospheric Photochemistry of Aldehydes: Norrish Type 1, Triple Fragmentation, and Methylketene Formation from Propanal.

Wavelength and pressure dependent quantum yields (ϕ, QYs) of propanal photolysis have been measured for photolysis wavelengths, λ = 300-330 nm, and buffer gases of 3-10 Torr propanal and 0-757 Torr N2. Following laser photolysis, three photochemical pathways were established, using Fourier transform infrared spectroscopy of the stable end-products. Photolysis is dominated by the Norrish Type 1 reaction, which has been reported previously, but with inconsistent quantum yields. The propanal α-hydrogen leads to a 4-center elimination of H2, as observed in CH3CHO, here leading to methylketene. The presence of hydrogen attached to the ß-carbon allows a new photochemical pathway: concerted triple fragmentation into CO + H2 + C2H4 via a 5-center transition state. Neither of these channels has been reported previously. No evidence for the previously reported C2H6 + CO, C2H4 + H2CO or CH3 + CH2CHO channels, nor for phototautomerization to 1-propenol (CH3CH═CHOH) was found. Modeling of the wavelength, pressure and collision partner dependence of the QYs allows us to reconcile the previous NT1a results and make recommendations for the quantum yields of all three channels under tropospheric conditions. The general impact of ß-hydrogen atoms in the photochemistry of aldehydes is to open up new pathways from cyclic transition states and to reduce the importance of other photolysis or isomerization channels.

Photophysical oxidation of HCHO produces HO2 radicals.

Formaldehyde, HCHO, is the highest-volume carbonyl in the atmosphere. It absorbs sunlight at wavelengths shorter than 330 nm and photolyses to form H and HCO radicals, which then react with O2 to form HO2. Here we show HCHO has an additional HO2 formation pathway. At photolysis energies below the energetic threshold for radical formation we directly detect HO2 at low pressures by cavity ring-down spectroscopy and indirectly detect HO2 at 1 bar by Fourier-transform infrared spectroscopy end-product analysis. Supported by electronic structure theory and master equation simulations, we attribute this HO2 to photophysical oxidation (PPO): photoexcited HCHO relaxes non-radiatively to the ground electronic state where the far-from-equilibrium, vibrationally activated HCHO molecules react with thermal O2. PPO is likely to be a general mechanism in tropospheric chemistry and, unlike photolysis, PPO will increase with increasing O2 pressure.

Rotational resonances in the H2CO roaming reaction are revealed by detailed correlations.

Since its discovery 16 years ago, roaming has become a ubiquitous mechanism in molecular photochemistry. Its general features are now understood, but little detail is known about how the potential energy surface (PES) determines reaction outcomes. We performed detailed experiments on formaldehyde (H2CO) photodissociation and determined fully correlated quantum state distributions of the molecular hydrogen and carbon monoxide products. These experiments reveal previously undetected bimodal carbon monoxide rotational distributions. Insights from classical trajectory calculations demonstrate that these features arise from resonances as the PES directs the reaction into cis and trans O-C-H···H critical geometries, which produce rebound and stripping mechanisms, respectively. These subtle and pervasive effects demonstrate additional complexity in this prototypical roaming reaction, which we expect to be general. They also provide detailed benchmarks for predictive theories of roaming.

Structural Effects on the Norrish Type I α-Bond Cleavage of Tropospherically Important Carbonyls.

Norrish Type I (NTI) α-bond cleavage is the dominant photolysis mechanism in small carbonyls and is an important source of radicals in the troposphere. In nonsymmetric species two cleavages are possible, NTIa and NTIb, forming larger and smaller alkyl radicals, respectively. For a data set of 20 small, atmospherically relevant carbonyls we predict NTIa and NTIb thresholds on the S0, S1, and T1 electronic states. The calculated NTIa T1 thresholds give a mean absolute deviation (MAD) of 5.8 kJ/mol with respect to the available experimental thresholds of five carbonyls. In addition, the intrinsic barrier heights to dissociation on the S0, S1, and T1 electronic states are predicted. We find RI-B2GP-PLYP/def2-TZVP calculations on S0 and unrestricted RI-B2GP-PLYP/def2-TZVP calculations on T1 give MADs of 6.1 kJ/mol for S0 asymptotic energies and 6.3 kJ/mol for S0 → T1 0-0 excitation energies, with respect to available experimental data. A composite method is used to determine S1 thresholds, with bt-STEOM-CCSD/cc-pVQZ calculation of vertical excitation energies and TD-RI-B3LYP/def2-TZVP calculations on S1, which achieves a MAD of 7.2 kJ/mol, with respect to experimental 0-0 excitation energies. Our calculations suggest, with the exception of bifunctional carbonyls and enones, NTI reactions on S1 are unlikely to be important at tropospherically relevant photolysis energies (<400 kJ/mol). In contrast, at these energies almost all possible NTI channels on T1 are open, and all barrierless S0 NTI dissociations are accessible. Our calculations allow a number of structural effects on both 0-0 excitation energies and intrinsic reaction barriers, on a given electronic state, to be elucidated and rationalized.

Dynamics and quantum yields of H2 + CH2CO as a primary photolysis channel in CH3CHO.

The first experimental observation of the primary photochemical channel of acetaldehyde leading to the formation of ketene (CH2CO) and hydrogen (H2) molecular products is reported. Acetaldehyde (CH3CHO) was photolysed in a molecular beam at 305.6 nm and the resulting H2 product characterized using velocity-map ion (VMI) imaging. Resonance-enhanced multiphoton ionization (REMPI), via two-photon excitation to the double-well EF 1Σ state, was used to state-selectively ionize the H2 and determine angular momentum distributions for H2 (ν = 0) and H2 (ν = 1). Velocity-map ion images were obtained for H2 (ν = 0 and 1, J = 5), allowing the total translational energy release of the photodissociation process to be determined. Following photolysis of CH3CHO in a gas cell, the CH2CO co-fragment was identified, using Fourier transform infrared spectroscopy, by its characteristic infrared absorption at 2150 cm-1. The measured quantum yield of the CH2CO + H2 product channel at 305.0 nm is φ = 0.0075 ± 0.0025 for both 15 Torr of neat CH3CHO and a mixture with 745 Torr of N2. Although small, this result has implications for the atmospheric photochemistry of carbonyls and this reaction represents a new tropospheric source of H2. Quasi-classical trajectory (QCT) simulations on a zero-point energy corrected reaction-path potential are also performed. The experimental REMPI and VMI image distributions are not consistent with the QCT simulations, indicating a non reaction-path mechanism should be considered.

Zero-point energy conservation in classical trajectory simulations: Application to H2CO.

A new approach for preventing zero-point energy (ZPE) violation in quasi-classical trajectory (QCT) simulations is presented and applied to H2CO "roaming" reactions. Zero-point energy may be problematic in roaming reactions because they occur at or near bond dissociation thresholds and these channels may be incorrectly open or closed depending on if, or how, ZPE has been treated. Here we run QCT simulations on a "ZPE-corrected" potential energy surface defined as the sum of the molecular potential energy surface (PES) and the global harmonic ZPE surface. Five different harmonic ZPE estimates are examined with four, on average, giving values within 4 kJ/mol-chemical accuracy-for H2CO. The local harmonic ZPE, at arbitrary molecular configurations, is subsequently defined in terms of "projected" Cartesian coordinates and a global ZPE "surface" is constructed using Shepard interpolation. This, combined with a second-order modified Shepard interpolated PES, V, allows us to construct a proof-of-concept ZPE-corrected PES for H2CO, V eff, at no additional computational cost to the PES itself. Both V and V eff are used to model product state distributions from the H + HCO → H2 + CO abstraction reaction, which are shown to reproduce the literature roaming product state distributions. Our ZPE-corrected PES allows all trajectories to be analysed, whereas, in previous simulations, a significant proportion was discarded because of ZPE violation. We find ZPE has little effect on product rotational distributions, validating previous QCT simulations. Running trajectories on V, however, shifts the product kinetic energy release to higher energy than on V eff and classical simulations of kinetic energy release should therefore be viewed with caution.

Photo-tautomerization of acetaldehyde as a photochemical source of formic acid in the troposphere.

Organic acids play a key role in the troposphere, contributing to atmospheric aqueous-phase chemistry, aerosol formation, and precipitation acidity. Atmospheric models currently account for less than half the observed, globally averaged formic acid loading. Here we report that acetaldehyde photo-tautomerizes to vinyl alcohol under atmospherically relevant pressures of nitrogen, in the actinic wavelength range, λ = 300-330 nm, with measured quantum yields of 2-25%. Recent theoretical kinetics studies show hydroxyl-initiated oxidation of vinyl alcohol produces formic acid. Adding these pathways to an atmospheric chemistry box model (Master Chemical Mechanism) demonstrates increased formic acid concentrations by a factor of ~1.7 in the polluted troposphere and a factor of ~3 under pristine conditions. Incorporating this mechanism into the GEOS-Chem 3D global chemical transport model reveals an estimated 7% contribution to worldwide formic acid production, with up to 60% of the total modeled formic acid production over oceans arising from photo-tautomerization.

Path integrals with higher order actions: Application to realistic chemical systems.

Quantum thermodynamic parameters can be determined using path integral Monte Carlo (PIMC) simulations. These simulations, however, become computationally demanding as the quantum nature of the system increases, although their efficiency can be improved by using higher order approximations to the thermal density matrix, specifically the action. Here we compare the standard, primitive approximation to the action (PA) and three higher order approximations, the Takahashi-Imada action (TIA), the Suzuki-Chin action (SCA) and the Chin action (CA). The resulting PIMC methods are applied to two realistic potential energy surfaces, for H2O and HCN-HNC, both of which are spectroscopically accurate and contain three-body interactions. We further numerically optimise, for each potential, the SCA parameter and the two free parameters in the CA, obtaining more significant improvements in efficiency than seen previously in the literature. For both H2O and HCN-HNC, accounting for all required potential and force evaluations, the optimised CA formalism is approximately twice as efficient as the TIA formalism and approximately an order of magnitude more efficient than the PA. The optimised SCA formalism shows similar efficiency gains to the CA for HCN-HNC but has similar efficiency to the TIA for H2O at low temperature. In H2O and HCN-HNC systems, the optimal value of the a1 CA parameter is approximately 13, corresponding to an equal weighting of all force terms in the thermal density matrix, and similar to previous studies, the optimal α parameter in the SCA was ∼0.31. Importantly, poor choice of parameter significantly degrades the performance of the SCA and CA methods. In particular, for the CA, setting a1 = 0 is not efficient: the reduction in convergence efficiency is not offset by the lower number of force evaluations. We also find that the harmonic approximation to the CA parameters, whilst providing a fourth order approximation to the action, is not optimal for these realistic potentials: numerical optimisation leads to better approximate cancellation of the fifth order terms, with deviation between the harmonic and numerically optimised parameters more marked in the more quantum H2O system. This suggests that numerically optimising the CA or SCA parameters, which can be done at high temperature, will be important in fully realising the efficiency gains of these formalisms for realistic potentials.

The energy dependence of CO(v,J) produced from H2CO via the transition state, roaming, and triple fragmentation channels.

The dynamics of CO production from photolysis of H2CO have been explored over a 8000 cm-1 energy range (345 nm-266 nm). Two-dimensional ion imaging, which simultaneously measures the speed and angular momentum distribution of a photofragment, was used to characterise the distribution of rotational and translational energy and to quantify the branching fraction of roaming, transition state (TS), and triple fragmentation (3F) pathways. The rotational distribution for the TS channel broadens significantly with increasing energy, while the distribution is relatively constant for the roaming channel. The branching fraction from roaming is also relatively constant at 20% of the observed CO. Above the 3F threshold, roaming decreases in favour of triple fragmentation. Combining the present data with our previous study on the H-atom branching fractions and published quantum yields for radical and molecular channels, absolute quantum yields were determined for all five dissociation channels for the entire S1←S0 absorption band, covering almost 8000 cm-1 of excitation energy. The S0 radical and TS molecular channels are the most important over this energy range. The absolute quantum yield of roaming is fairly constant ∼5% at all energies. The T1 radical channel is important (20%-40%) between 1500 and 4000 cm-1 above the H + HCO threshold, but becomes unimportant at higher energy. Triple fragmentation increases rapidly above its threshold reaching a maximum of 5% of the total product yield at the highest energy.

Infrared Spectra of Gas-Phase 1- and 2-Propenol Isomers.

Fourier transform infrared spectra of isolated 1-propenol and 2-propenol in the gas-phase have been collected in the range of 900-3800 cm-1, and the absolute infrared absorption cross sections reported for the first time. Both cis and trans isomers of 1-propenol were observed with the trans isomer in greater abundance. Syn and anti conformers of both 1- and 2-propenol were also observed, with abundance consistent with thermal population. The FTIR spectrum of the smaller ethenol (vinyl alcohol) was used as a benchmark for our computational results. As a consequence, its spectrum has been partially reassigned resulting in the first report of the anti-ethenol conformer. Electronic structure calculations were used to support our experimental results and assign vibrational modes for the most abundant isomers, syn-trans-1-propenol and syn-2-propenol.

Path integral Monte Carlo simulations of H2 adsorbed to lithium-doped benzene: A model for hydrogen storage materials.

Finite temperature quantum and anharmonic effects are studied in H2-Li(+)-benzene, a model hydrogen storage material, using path integral Monte Carlo (PIMC) simulations on an interpolated potential energy surface refined over the eight intermolecular degrees of freedom based upon M05-2X/6-311+G(2df,p) density functional theory calculations. Rigid-body PIMC simulations are performed at temperatures ranging from 77 K to 150 K, producing both quantum and classical probability density histograms describing the adsorbed H2. Quantum effects broaden the histograms with respect to their classical analogues and increase the expectation values of the radial and angular polar coordinates describing the location of the center-of-mass of the H2 molecule. The rigid-body PIMC simulations also provide estimates of the change in internal energy, ΔUads, and enthalpy, ΔHads, for H2 adsorption onto Li(+)-benzene, as a function of temperature. These estimates indicate that quantum effects are important even at room temperature and classical results should be interpreted with caution. Our results also show that anharmonicity is more important in the calculation of U and H than coupling-coupling between the intermolecular degrees of freedom becomes less important as temperature increases whereas anharmonicity becomes more important. The most anharmonic motions in H2-Li(+)-benzene are the "helicopter" and "ferris wheel" H2 rotations. Treating these motions as one-dimensional free and hindered rotors, respectively, provides simple corrections to standard harmonic oscillator, rigid rotor thermochemical expressions for internal energy and enthalpy that encapsulate the majority of the anharmonicity. At 150 K, our best rigid-body PIMC estimates for ΔUads and ΔHads are -13.3 ± 0.1 and -14.5 ± 0.1 kJ mol(-1), respectively.

H2 Adsorption in a Porous Crystal: Accurate First-Principles Quantum Simulation.

A general method is presented for constructing, from ab initio quantum chemistry calculations, the potential energy surface (PES) for H2 absorbed in a porous crystalline material. The method is illustrated for the metal-organic framework material MOF-5. Rigid body quantum diffusion Monte Carlo simulations are used in the construction of the PES and to evaluate the quantum ground state of H2 in MOF-5, the zero-point energy, and the enthalpy of adsorption at 0 K.

"Plug-and-Play" potentials: Investigating quantum effects in (H2)2-Li(+)-benzene.

Quantum and anharmonic effects are investigated in (H2)2-Li(+)-benzene, a model for hydrogen adsorption in metal-organic frameworks and carbon-based materials, using rigid-body diffusion Monte Carlo (RBDMC) simulations. The potential-energy surface (PES) is calculated as a modified Shepard interpolation of M05-2X/6-311+G(2df,p) electronic structure data. The RBDMC simulations yield zero-point energies (ZPE) and probability density histograms that describe the ground-state nuclear wavefunction. Binding a second H2 molecule to the H2-Li(+)-benzene complex increases the ZPE of the system by 5.6 kJ mol(-1) to 17.6 kJ mol(-1). This ZPE is 42% of the total electronic binding energy of (H2)2-Li(+)-benzene and cannot be neglected. Our best estimate of the 0 K binding enthalpy of the second H2 to H2-Li(+)-benzene is 7.7 kJ mol(-1), compared to 12.4 kJ mol(-1) for the first H2 molecule. Anharmonicity is found to be even more important when a second (and subsequent) H2 molecule is adsorbed; use of harmonic ZPEs results in significant error in the 0 K binding enthalpy. Probability density histograms reveal that the two H2 molecules are found at larger distance from the Li(+) ion and are more confined in the θ coordinate than in H2-Li(+)-benzene. They also show that both H2 molecules are delocalized in the azimuthal coordinate, ϕ. That is, adding a second H2 molecule is insufficient to localize the wavefunction in ϕ. Two fragment-based (H2)2-Li(+)-benzene PESs are developed. These use a modified Shepard interpolation for the Li(+)-benzene and H2-Li(+)-benzene fragments, and either modified Shepard interpolation or a cubic spline to model the H2-H2 interaction. Because of the neglect of three-body H2, H2, Li(+) terms, both fragment PESs lead to overbinding of the second H2 molecule by 1.5 kJ mol(-1). Probability density histograms, however, indicate that the wavefunctions for the two H2 molecules are effectively identical on the "full" and fragment PESs. This suggests that the 1.5 kJ mol(-1) error is systematic over the regions of configuration space explored by our simulations. Notwithstanding this, modified Shepard interpolation of the weak H2-H2 interaction is problematic and we obtain more accurate results, at considerably lower computational cost, using a cubic spline interpolation. Indeed, the ZPE of the fragment-with-spline PES is identical, within error, to the ZPE of the full PES. This fragmentation scheme therefore provides an accurate and inexpensive method to study higher hydrogen loading in this and similar systems.

A multi-agent quantum Monte Carlo model for charge transport: Application to organic field-effect transistors.

We have developed a multi-agent quantum Monte Carlo model to describe the spatial dynamics of multiple majority charge carriers during conduction of electric current in the channel of organic field-effect transistors. The charge carriers are treated by a neglect of diatomic differential overlap Hamiltonian using a lattice of hydrogen-like basis functions. The local ionization energy and local electron affinity defined previously map the bulk structure of the transistor channel to external potentials for the simulations of electron- and hole-conduction, respectively. The model is designed without a specific charge-transport mechanism like hopping- or band-transport in mind and does not arbitrarily localize charge. An electrode model allows dynamic injection and depletion of charge carriers according to source-drain voltage. The field-effect is modeled by using the source-gate voltage in a Metropolis-like acceptance criterion. Although the current cannot be calculated because the simulations have no time axis, using the number of Monte Carlo moves as pseudo-time gives results that resemble experimental I/V curves.

Generating accurate dipole moment surfaces using modified Shepard interpolation.

We outline an approach for building molecular dipole moment surfaces using modified Shepard interpolation. Our approach is highly automated, requires minimal parameterization, and is iteratively improvable. Using the water molecule as a test case, we investigate how different aspects of the interpolation scheme affect the rate of convergence of calculated IR spectral line intensities. It is found that the interpolation scheme is sensitive to coordinate singularities present at linear geometries. Due to the generally monotonic nature of the dipole moment surface, the one-part weight function is found to be more effective than the more complicated two-part variant, with first-order interpolation also giving better-than-expected results. Almost all sensible schemes for choosing interpolation reference data points are found to exhibit acceptable convergence behavior.

Quantum effects and anharmonicity in the H2-Li(+)-benzene complex: a model for hydrogen storage materials.

Quantum and anharmonic effects are investigated in H2-Li(+)-benzene, a model for hydrogen adsorption in metal-organic frameworks and carbon-based materials. Three- and 8-dimensional quantum diffusion Monte Carlo (QDMC) and rigid-body diffusion Monte Carlo (RBDMC) simulations are performed on potential energy surfaces interpolated from electronic structure calculations at the M05-2X/6-31+G(d,p) and M05-2X/6-311+G(2df,p) levels of theory using a three-dimensional spline or a modified Shepard interpolation. These calculations investigate the intermolecular interactions in this system, with three- and 8-dimensional 0 K H2 binding enthalpy estimates, ΔH(bind) (0 K), being 16.5 kJ mol(-1) and 12.4 kJ mol(-1), respectively: 0.1 and 0.6 kJ mol(-1) higher than harmonic values. Zero-point energy effects are 35% of the value of ΔH(bind) (0 K) at M05-2X/6-311+G(2df,p) and cannot be neglected; uncorrected electronic binding energies overestimate ΔHbind (0 K) by at least 6 kJ mol(-1). Harmonic intermolecular binding enthalpies can be corrected by treating the H2 "helicopter" and "ferris wheel" rotations as free and hindered rotations, respectively. These simple corrections yield results within 2% of the 8-dimensional anharmonic calculations. Nuclear ground state probability density histograms obtained from the QDMC and RBDMC simulations indicate the H2 molecule is delocalized above the Li(+)-benzene system at 0 K.

Experimental and theoretical investigation of triple fragmentation in the photodissociation dynamics of H2CO.

The photodissociation dynamics of H2CO molecules at energies bracketing the triple fragmentation threshold were investigated using velocity map ion imaging of the H-atom fragments. An algorithm was developed to model the experimental results as a two-step process: initially barrierless C-H bond fission on the S0 potential energy surface to form H + HCO, followed by secondary fragmentation of those HCO radicals with sufficient internal energy to overcome the small exit channel barrier on the HCO surface to form H + CO. Our model treats the first step using phase space theory (PST) and the second using a combined PST-impulsive model, with a tunneling correction. Experimentally, triple fragmentation reaches 25% of the radical (H + HCO) channel photochemical yield at energies about 1500 cm(-1) above the barrier for breaking the second bond. In addition, the triplet (T1) channel appears to reduce in importance after the barrier on the T1 surface is exceeded, slowly decreasing to <10% of the total radical yield at higher energy. The double PST-impulsive model provides a good fit to the experimental H-atom speed and energy distributions for H2CO dissociation on S0 spanning >7000 cm(-1) of available energy.

A phase space theory for roaming reactions.

We describe a new, simple theory for predicting the branching fraction of products in roaming reactions, compared to the analogous barrierless bond dissociation products. The theory uses a phase space theory (PST) formalism to divide reactive states in the bond dissociation channel into states with enough translational energy to dissociate and states that may roam. Two parameters are required, ΔEroam, the energy difference between the bond dissociation threshold and the roaming threshold, and the roaming probability, Proam, the probability that states that may roam do roam rather than recombine to form reactants. The PST-roaming theory is tested against experimental and theoretical data on the dissociation dynamics of H2CO, NO3, and CH3CHO. The theory accurately models the relative roaming to bond dissociation branching fraction over the experimental or theoretical energy range available in the literature for each species. For H2CO, fixing ΔEroam = 146 cm(-1), the midpoint of the experimental bounds for the roaming threshold, we obtain Proam = 1. The best-fit value, ΔEroam = 161 cm(-1), is also consistent with the experimental bounds. Using this value, the relative roaming to dissociation branching ratios are predicted to be similar in D2CO and H2CO, consistent with experimental observation. For NO3, we fix ΔEroam = 258.6 cm(-1), the experimental threshold for NO + O2 production, and we model low-temperature experimental branching fractions using the experimental rotational and vibrational temperatures of Trot = 0 K and Tvib = 300 K. The best fit to the experimental data is obtained for Proam = 0.0075, with this very small Proam being consistent with the known geometric constraints to formation of NO + O2. Using Proam = 0.0075, our PST-roaming theory also accurately predicts the low-temperature NO yield spectrum and quantum yield data for room-temperature NO3 photolysis. For CH3CHO, we fix ΔEroam = 385 cm(-1), based on theoretical calculations, and obtain a best-fit value of Proam = 0.21, fitting to reduced dimensional trajectory calculations. These values of ΔEroam and Proam yield PST-roaming theory results that are also consistent with two experimental room-temperature data points. The combination of other kinetic theories and the PST-roaming theory will provide rate coefficients for roaming reactions.

Modeling molecular response in uniform and non-uniform electric fields.

The response of a molecule to an electric field E, often a model of environment, can be expressed in terms of a sum of power series expansions. We investigate the accuracy and limits of applicability of this expression using one-, two-, and three-dimensional models of the hydrogen-bonded complex, ClH:NH(3). Energetic, structural, and vibrational spectroscopic characteristics are determined at first- and second-order in E and [nabla]E and compared with ab initio values for a range of uniform and non-uniform electric fields chosen to simulate molecular environments. It is found that even at field strengths large enough to cause dramatic structural change in the complex, energetic, structural, and vibrational spectroscopic characteristics are accurately calculated using only terms linear in E and [nabla]E. These results suggest that knowledge of the zero-field molecular potential energy, dipole, and quadrupole moment surfaces may be sufficient to accurately model the interaction of a molecule with a wide range of chemical environments.

Photo-tautomerization of acetaldehyde to vinyl alcohol: a potential route to tropospheric acids.

Current atmospheric models underestimate the production of organic acids in the troposphere. We report a detailed kinetic model of the photochemistry of acetaldehyde (ethanal) under tropospheric conditions. The rate constants are benchmarked to collision-free experiments, where extensive photo-isomerization is observed upon irradiation with actinic ultraviolet radiation (310 to 330 nanometers). The model quantitatively reproduces the experiments and shows unequivocally that keto-enol photo-tautomerization, forming vinyl alcohol (ethenol), is the crucial first step. When collisions at atmospheric pressure are included, the model quantitatively reproduces previously reported quantum yields for photodissociation at all pressures and wavelengths. The model also predicts that 21 ± 4% of the initially excited acetaldehyde forms stable vinyl alcohol, a known precursor to organic acid formation, which may help to account for the production of organic acids in the troposphere.