Search

A density functional theory investigation of the chemical reactivity of the S=C=S...C1 complex toward 2,3-dimethylbutane (DMB) has been reported and compared to the related reaction of a C1 atom with DMB. The DFT calculations indicate that the reaction of C1 atom with DMB would have little tertiary selectivity toward chlorination consistent with experimental results for photochlorination of alkanes in invert solvents.

The complete active space SCF method has been used to fully optimize the potential energy surfaces controlling hydrazoic acid dissociation to HN + N2 and H + N3 in the ground, excited triplet and singlet states. This research gives an insight into the mechanism of the ultraviolet protodissociation of HN3 at a wavelength range from 355 nm to 248 nm.

Research was conducted to examine the characterize the minimum-energy structure of the 7-hydroxyquinoline (7HQ)-(H2O)(sub n) (n = 1-3) complexes and the potential energy surfaces governing triple proton transfer in 7HQ-(H2O)2 and quadruple proton transfer in 7HQ-(H2O)3 using different ab initio quantum mechanical methods. The calculated results and experimental data were compared where available.

Machine learning is capable of effectively predicting the potential energies of molecules in the presence of high-quality data sets. Its application in the construction of ground- and excited-state potential energy surfaces is attractive to accelerate nonadiabatic molecular dynamics simulations of photochemical reactions. Because of the huge computational cost of excited-state electronic structure calculations, the construction of a high-quality data set becomes a bottleneck. In the present work, we first built two data sets. One was obtained from surface hopping dynamics simulations at the semiempirical OM2/MRCI level. Another was extracted from the dynamics trajectories at the CASSCF level, which was reported previously. The ground- and excited-state potential energy surfaces of ethylene-bridged azobenzene at the CASSCF computational level were constructed based on the former low-level data set. Although non-neural network machine learning methods can achieve good or modest performance during the training process, only neural network models provide reliable predictions on the latter external test data set. The BPNN and SchNet combined with the Δ-ML scheme and the force term in the loss functions are recommended for dynamics simulations. Then, we performed excited-state dynamics simulations of the photoisomerization of ethylene-bridged azobenzene on machine learning potential energy surfaces. Compared with the lifetimes of the first excited state (S 1 ) estimated at different computational levels, our results on the E isomer are in good agreement with the high-level estimation. However, the overestimation of the Z isomer is unimproved. It suggests that smaller errors during the training process do not necessarily translate to more accurate predictions on high-level potential energies or better performance on nonadiabatic dynamics simulations, at least in the present case.

Herein, we employed linear-response time-dependent functional theory nonadiabatic dynamic simulations to explore the photoinduced exciton dynamics of a chiral single-walled carbon nanotube CNT(6,5) covalently doped with a 4-nitrobenzyl group (CNT65-NO2). The results indicate that the introduction of a sp 3 defect leads to the splitting of the degenerate VBM/VBM-1 and CBM/CBM+1 states. Both the VBM upshift and the CBM downshift are responsible for the experimentally observed redshifted E 11 * trapping state. The simulations reveal that the photoinduced exciton relaxation dynamics completes within 500 fs, which is consistent with the experimental work. On the other hand, we also conducted the nonadiabatic carrier (electron and hole) dynamic simulations, which completely ignore the excitonic effects. The comparison demonstrates that excitonic effects are indispensable. Deep analyses show that such effects induce several dark states, which play an important role in regulating the photoinduced dynamics of CNT65-NO2. The present work demonstrates the importance of including excitonic effects in simulating photoinduced processes of carbon nanotubes. In addition, it not only rationalizes previous experiments but also provides valuable insights that will help in the future rational design of novel covalently doped carbon nanotubes with superior photoluminescent properties.

Interface orientation between zinc phthalocyanine (ZnPc) and fullerene (C 60 ) affects their interfacial charge separation dynamics; however, the underlying physical origin is still elusive. In this work, we have employed the time-dependent density functional theory (TDDFT) method to explore excited-state properties of ZnPc and C 60 heterojunctions with both face-on and edge-on configurations. Spectroscopically bright absorption is from locally excited (LE) singlet excitons within ZnPc. In the face-on configuration, LE excitons are much higher in energy than charge-transfer (CT) excitons, thereby making charge separation process favorable. However, in the edge-on configuration, LE excitons are the lowest ones and CT ones are higher in energy; thus, charge separation is not efficient. Subsequently, we have carried out TDDFT-based nonadiabatic dynamics method to simulate photoinduced exciton and charge separation dynamics of ZnPc and C 60 heterojunctions with both edge-on and face-on configurations. In the former, there are no exciton transfer and charge separation processes observed within 300 fs simulation time; while, in the latter, fragment-based electronic transition density matrix analysis reveals that only LE excitons |C 60 ZnPC*⟩ and CT excitons |C 60 - ZnPC + ⟩ are involved. The exciton transfer from |C 60 ZnPC*⟩ to |C 60 - ZnPC + ⟩ is completed within about 100 fs in which charge separation takes place with electron-hole distances increasing from 1.0 to 4.5 Å. This exciton transfer process is essentially in company with electron transfer from ZnPc to C 60 but almost not involving hole transfer. These gained insights not only rationalize experiments but also enrich our knowledge to design donor-acceptor orientations to optimize organic photovoltaic performance.

Recently, we have developed a multilayer energy-based fragment (MLEBF) method to describe excited states of large systems in which photochemically active and inert regions are separately treated with multiconfigurational and single-reference electronic structure method and their mutual polarization effects are naturally described within the many-body expansion framework. This MLEBF method has been demonstrated to provide highly accurate energies and gradients. In this work, we have further derived the MLEBF method with which highly accurate excited-state Hessian matrices of large systems are efficiently constructed. Moreover, in combination with recently proposed embedded atom neural network (EANN) model we have developed a machine learning (ML) accelerated MLEBF method (i.e., ML-MLEBF) in which photochemically inert region is entirely replaced with trained ML models. ML-MLEBF is found to improve computational efficiency of Hessian matrices in particular for large systems. Furthermore, both MLEBF and ML-MLEBF methods are highly parallel and exhibit low-scaling computational cost with multiple CPUs. The present developments could motivate combining various ML techniques with fragment-based electronic structure methods to explore Hessian-matrix-based excited-state properties of large systems.

Two-dimensional (2D) heterostructures of black phosphorus (BP)/bismuth vanadate (BVO) have attracted much attention due to their potential uses in photocatalytic water splitting. However, the interfacial photoinduced electron- and hole-transfer dynamics are not explored computationally. Herein, we have used density functional theory (DFT) calculations and DFT-based fewest-switches surface-hopping dynamics simulations to investigate the light-driven electron and hole dynamics taking place at the interface of BP and the BVO(010) surface. Our results show that the BP monolayer is adsorbed on BVO(010) via van der Waals interaction. Upon irradiation, the electron transfer takes place from BP to BVO(010) within 500 fs but with two distinct processes. In the first phase, the electron transfer proceeds adiabatically and is mainly driven by atomic motions. In the second phase, the electron transfer decays very slowly. The hole-transfer dynamics from BVO(010) to BP exhibits a similar ultrafast decay in the first stage followed by a slow decay; however, there is a comparable amount of hole trapped in a BP state due to a large energy gap from its higher state. These insights may be useful for the design of novel photocatalytic water-splitting materials.

Understanding of photoinduced interfacial carrier dynamics in organic-transition metal dichalcogenides heterostructures is very important for the enhancement of their potential photoelectronic conversion efficiencies. In this work we have used density functional theory (DFT) calculations and DFT-based fewest-switches surface-hopping dynamics simulations to explore the photoinduced hole transfer and subsequent nonadiabatic electron-hole recombination dynamics taking place at the interface of pentacene and MoS 2 in pentacene@MoS 2 . Upon photoexcitation the electronic transition mainly occurs on the MoS 2 monolayer, which corresponds to moving an electron to the MoS 2 conduction band. As a result, a hole is left in the valence band. This hole state is energetically lower than certain occupied states of the pentacene molecule; thus, the interfacial hole transfer from MoS 2 to pentacene is favorable in energy. In terms of nonadiabatic dynamics simulations, the hole transfer time to the HOMO-1 state of the pentacene is estimated to be about 600 fs; however, the following hole relaxation process from HOMO-1 to HOMO takes much longer time of ca. 15 ps due to the large energy gap between HOMO-1 and HOMO. Moreover, our results also show that the subsequent radiationless recombination process between the hole transferred to the pentacene molecule and the remaining electron on the MoS 2 CBM needs about 10.2 ns. The computational results shed important mechanistic insights on the interfacial carrier dynamics of mixed-dimensional pentacene@MoS 2 . These insights could help to design excellent interfaces for organic-TMDs heterostructures.

The mixed quantum-classical dynamical approaches have been widely used to study nonadiabatic phenomena in photochemistry and photobiology, in which the time evolutions of the electronic and nuclear subsystems are treated based on quantum and classical mechanics, respectively. The key issue is how to deal with coherence and decoherence during the propagation of the two subsystems, which has been the subject of numerous investigations for a few decades. A brief description on Ehrenfest mean-field and surface-hopping (SH) methods is first provided, and then different algorithms for treatment of quantum decoherence are reviewed in the present paper. More attentions were paid to quantum trajectory mean-field (QTMF) method under the picture of quantum measurements, which is able to overcome the overcoherence problem. Furthermore, the combined QTMF and SH algorithm is proposed in the present work, which takes advantages of the QTMF and SH methods. The potential to extend the applicability of the QTMF method was briefly discussed, such as the generalization to other type of nonadiabatic transitions, the combination with multiscale computational models, and possible improvements on its accuracy and efficiency by using machine-learning techniques.

Herein, we have employed OM2/MRCI-based full-dimensional nonadiabatic dynamics simulations to explore the photoisomerization and subsequent excited-state decay of a macrocyclic cyclobiazobenzene molecule. Two S 1 /S 0 conical intersection structures are found to be responsible for the excited-state decay. Related to these two conical intersections, we found two stereoselective photoisomerization and excited-state decay pathways, which correspond to the clockwise and counterclockwise rotation motions with respect to the N═N bond of the azo group. In both pathways, the excited-state isomerization is ultrafast and finishes within ca. 69 fs, but the clockwise isomerization channel is much more favorable than the counterclockwise one with a ratio of 74% versus 26%. Importantly, the present work demonstrates that stereoselective pathways exist not only in the photoisomerization of isolated azobenzene (AB)-like systems but also in macrocyclic systems with multiple ABs. This finding could provide useful insights for understanding and controlling the photodynamics of macrocyclic nanostructures with AB units as the main building units.

Photoinduced chemical reactions of organic compounds at the marine boundary layer have recently attracted significant experimental attention because this kind of photoreactions has been proposed to have substantial impact on local new particle formation and their photoproducts could be a source of secondary organic aerosols. In this work, we have employed first-principles density functional theory method combined with cluster models to systematically explore photochemical reaction pathways of nonanoic acids (NAs) to form volatile saturated and unsaturated C 9 and C 8 aldehydes at air-water interfaces. On the basis of the results, we have found that the formation of C 9 aldehydes is not initiated by intermolecular Norrish type II reaction between two NAs but by intramolecular T 1 C-O bond fission of NA generating acyl and hydroxyl radicals. Subsequently, saturated C 9 aldehydes are formed through hydrogenation reaction of acyl radical by another intact NA. Following two dehydrogenation reactions, unsaturated C 9 aldehydes are generated. In parallel, the pathway to C 8 aldehydes is initiated by T 1 C-C bond fission of NA, which generates octyl and carboxyl radicals; then, an octanol is formed through recombination reaction of octyl with hydroxyl radical. In the following, two dehydrogenation reactions result into an enol intermediate from which saturated C 8 aldehydes are produced via NA-assisted intermolecular hydrogen transfer. Finally, two dehydrogenation reactions generate unsaturated C 8 aldehydes. In these reactions, water and NA molecules are found to play important roles. They significantly reduce relevant reaction barriers. Our work has also explored oxygenation reactions of NA with molecular oxygen and radical-radical dimerization reactions.

In this article, structures and energies of cyclopropenone in the low-lying electronic states have been determined by the CASSCF and MS-CASPT2 calculations with different basis sets. Two minimum-energy conical intersections (CI-1 and CI-2) between S0 and S1 were obtained and their topographic characters were characterized by the SA4-CAS(10,9) calculated energy gradients and nonadiabatic coupling vectors. The AIMS method was used to carry out nonadiabatic dynamics simulation with ab initio calculation performed at the SA4-CAS(10,9) level. On the basis of time evolution of wave functions simulated here, the S1 lifetime is fitted to be 125 fs with a pure exponential decay for the S1 electronic population. The CI-1 intersection is mainly responsible for ultrafast S1→S0 nonadiabatic transition and the photoinduced decarbonylation is a sequential process, where the first C-C bond is broken in the S1 state and fission of the second C-C bond occurs in the S0 state as a result of the S1→S0 internal conversion via the CI-1 region. As a minor channel through the CI-2 region, the decarbonylation proceeds in an asynchronous concerted way. Effects of the S1 excess energies and the S1-S0 energy gap on the nonadiabatic dynamics were examined, which reveals that the S1→S0 nonadiabatic transition occurs within a small energy gap and high-energy conical intersection regions can play an important role. The present study provides new insights into mechanistic photochemistry of cyclopropenones and reveals that the AIMS dynamics simulation at a high-accuracy ab initio level is a powerful tool for exploring a mechanism of an ultrafast photochemical reaction.

Intramolecularly hydrogen-bonded aromatic azo compound 1-cyclopropyldiazo-2-naphthol (CPDNO) exhibits complicated excited-state behaviors, e.g., wavelength-dependent photoinduced proton transfer and photoproducts. Its photochemistry differs from that of common aromatic azo compounds in which cis-trans photoisomerization is dominant. To rationalize the intriguing photochemistry of CPDNO at the atomic level, we have in this work employed the complete active space self-consistent field (CASSCF) and its second-order perturbation (CASPT2) methods to explore the S0, S1, and S2 potential-energy profiles relevant to enol-keto proton transfer and isomerization reactions. It is found that the proton transfer along the bright diabatic (1)ππ* potential-energy profile is almost barrierless, quickly forming the fluorescent (1)ππ* keto minimum. In this process, the dark (1)nπ* state is populated via a (1)ππ*/(1)nπ* crossing point, but the proton transfer on this dark state is suppressed heavily as a result of a large barrier. In addition, two deactivation paths that decay the S1 enol and keto minima to the S0 state, respectively, were uncovered. For the former, it is exoenergetic and thereby thermodynamically favorable; for the latter, it is a little endothermic (ca. 5 kcal/mol). Both are energetically allowable concerning the available total energy. Finally, on the basis of the present results, the experimentally observed wavelength-dependent photoproducts were explained very well.

This work uses the photodissociation of acetyl chloride to assess the utility of a recently proposed impulsive model when the dissociation occurs on an excited electronic state that is not repulsive in the Franck-Condon region. The impulsive model explicitly includes an average over the vibrational quantum states of acetyl chloride when it calculates an impact parameter for fission of the C-Cl bond, as well as the distribution of thermal energy in the photolytic precursor. The experimentally determined stability of the resulting acetyl radical to subsequent dissociation is the key observable that allows us to test the model's ability to predict the partitioning of energy between rotation and vibration of the radical. We compare the model's predictions for three different assumed geometries at which the impulsive force might act, generating the relative kinetic energy and the concomitant rotational energy in the acetyl radical. Assuming that the impulsive force acts at the transition state for C-Cl fission on the S(1) excited state gives a poor prediction; the model predicts far more energy in rotation of the acetyl radical than is consistent with the measured velocity map imaging spectrum of the stable radicals. The best prediction results from using a geometry derived from the classical trajectory calculations on the excited state potential energy surface. We discuss the insight gained into the excited state dissociation dynamics of acetyl chloride and, more generally, the utility of using the impulsive model in conjunction with excited state trajectory calculations to predict the partitioning of internal energy between rotation and vibration for radicals produced from the photolysis of halogenated precursors.

The density functional theory (DFT) method was used to explore the light-induced O(2) formation from H(2)O promoted by Ru(II) PNN complex in the present work. The elimination of H(2)O(2) was found to be highly endothermic, which is not in competition with the H(2)O elimination and hydrogen transfer. The calculated results reported here do not support the mechanism proposed in a recent experiment, where H(2)O(2) was suggested as an important intermediate for formation of O(2). We proposed a new mechanism for formation of the triplet O(2) molecule, which contains the two steps of the concerted hydrogen transfer and dehydration. The light-induced O(2) evolution from water promoted by the Ru(II) complex was found to be a nonadiabatic process. The O-O bond is formed along the T(1) pathway as a result of the efficient S(1) → T(1) intersystem crossing. All experimental findings on the light-induced O(2) evolution can be explained by the mechanism proposed in the present work.

Heme is a key cofactor of hemoproteins in which porphyrin is often found to be preferentially metalated by the iron cation. In our previous work [Feng, X. T.; Yu, J. G.; Lei, M.; Fang, W. H.; Liu, S. B. J. Phys. Chem. B 2009, 113, 13381], conceptual density functional theory (CDFT) descriptors have been applied to understand the metal-binding specificity of porphyrin. We found that the iron-porphyrin complex significantly differs in many aspects from porphyrin complexes with other metal cations except Ru, for which similar behaviors for the reactivity descriptors were discovered. In this study, we employ the spin-polarized version of CDFT to investigate the reactivity for a series of (pyridine)(n)-M(ll)-porphyrin complexes-where M = Mg, Ca, Cr, Mn, Co, Ni, Cu, Zn, Ru, and Cd, and n = 0, 1, and 2-to further appreciate the metal-binding specificity of porphyrin. Both global and local descriptors were examined within this framework. We found that, within the spin resolution, not only chemical reactivity descriptors from CDFT of the iron complex are markedly different from that of other metal complexes, but we also discovered substantial differences in reactivity descriptors between Fe and Ru complexes. These results confirm that spin properties play a highly important role in physiological functions of hemoproteins. Quantitative reactivity relationships have been revealed between global and local spin-polarized reactivity descriptors. These results contribute to our better understanding of the metal binding specificity and reactivity for heme-containing enzymes and other metalloproteins alike.

In the present work, density functional theory and canonical nonadiabatic Monte Carlo transition state theory have been used to investigate the histidine dissociation process from hexacoordinate heme in Ngb protein. The potential energy surfaces (PES) of the lowest singlet, triplet, and quintet states are calculated by stepwise optimization along with the histidine dissociation pathway. Based on the calculated two-dimensional PES, the histidine dissociation rates for the spin-forbidden processes via singlet to triplet and singlet to quintet transitions have been calculated by the nonadiabatic Monte Carlo transition state theory in canonical ensemble. The present study provides a quantitative description on spin-forbidden histidine dissociation processes.

Resonance Raman spectra were obtained for 2-hydroxybenzaldehyde (OHBA) in cyclohexane solution with excitation wavelengths in resonance with the first charge-transfer/proton-transfer (CT/PT) band absorption. These spectra indicate that the Franck-Condon region photodissociation dynamics have multidimensional character with motion predominantly along the nominal C=CH in-plane bend+ring deformation modes (nu9, nu10, nu14, nu16, nu18, nu19, nu20, nu26, nu30, nu31, and nu35) accompanied by a smaller amount of motion along the nominal C=O stretch mode (nu7), the nominal C=C-C(=O) in-plane bend modes (nu33 and nu37), and the nominal ring C-O-H in-plane bend modes (nu9 and nu14). A preliminary resonance Raman intensity analysis was done, and these results for the OHBA molecule were compared to results previously reported for the 2-hydroxyacetophenone (OHAP) molecule. Several proton-transfer tautomers in the ground and excited states were predicted from the results of B3LYP/cc-PVTZ, UB3LYP/cc-PVTZ, and CASSCF/cc-PVDZ level of theory computations. The differences and similarities between the CT/PT band resonance Raman spectra and the vibrational reorganizational energies for the OHBA molecule relative to those for the OHAP molecule are briefly discussed.

The potential energy surfaces of dissociation and elimination reactions for CH(3)COCl in the ground (S0) and first excited singlet (S1) states have been mapped with the different ab inito calculations. Mechanistic photodissociation of CH(3)COCl has been characterized through the intrinsic reaction coordinate and ab initio molecular dynamics calculations. The alpha-C-C bond cleavage along the S1 pathway leads to the fragments of COCl((2)A' ') and CH(3) ((2)A') in an excited electronic state and a high barrier exists on the pathway. This channel is inaccessible in energy upon photoexcitation of the CH(3)COCl molecules at 236 nm. The S1 alpha-C-Cl bond cleavage yields the Cl((2)P) and CH(3)CO(X(2)A') fragments in the ground state and there is very small or no barrier on the pathway. The S1 alpha-C-Cl bond cleavage proceeds in a time scale of picosecond in the gas phase, followed by CH(3)CO decomposition to CH(3) and CO. The barrier to the C-Cl bond cleavage on the S1 surface is significantly increased by effects of the argon matrix. The S1 alpha-C-Cl bond cleavage in the argon matrix becomes inaccessible in energy upon photoexcitation of CH(3)COCl at 266 nm. In this case, the excited CH(3)COCl(S1) molecules cannot undergo the C-Cl bond cleavage in a short period. The internal conversion from S1 to S0 becomes the dominant process for the CH(3)COCl(S1) molecules in the condensed phase. As a result, the direct HCl elimination in the ground state becomes the exclusive channel upon 266 nm photodissociation of CH(3)COCl in the argon matrix at 11 K.

The s-cis and s-trans isomer radical cations of hexafluoro-1,3-butadiene (s-cis-HFBD+ and s-trans-HFBD+) were generated by a gamma-irradiated solid solution of the neutral HFBD molecule in solid matrix at 77 K and observed by means of electron spin resonance (ESR) and electronic spectroscopies. In comparing the experimental isotropic and anisotropic 19F hyperfine splittings with the computational ones by the DFT B3LYP and MP2 methods, the generated s-cis-HFBD+ and s-trans-HFBD+ radical cations were confirmed to be 2A2 and 2Bg electronic ground states in C2v and C2h symmetries, respectively. The present spectroscopic study revealed that the relative abundance of s-cis-HFBD+ to s-trans-HFBD+ was 4.0 immediately after being formed by gamma-irradiation, and subsequently most s-cis-HFBD+ was isomerized to s-trans-HFBD+ by visible-light illumination with 500-600 nm wavelength. The process of nonplanar HFBD ionizing to form stable planar s-cis- and s-trans-HFBD+ and the reaction mechanism of the cis-to-trans photoisomerization were discussed by (MS-)CASPT2//CASSCF calculated vertical excitation energies (Tv) and torsional potential energy curves (TPECs) of HFBD and HFBD+.

A computational study on the transamination reaction of molecular complexes that consist of NH2CH2COOH + CH2O + nH2O, where n = 0, 1, 2, is presented. This work has allowed the description of the geometries of all the intermediates and transition states of the reactions, which can be described by five steps: carbinolamine formation, dehydration, 1,3 proton transfer, hydrolysis, and carbinolamine elimination. Among the five steps of the reaction, hydrolysis and elimination occur with the existence of general acid catalysis related to the carboxylic group. The water molecules can be involved in the reaction by performing as a proton-transfer carrier and a stabilizing zwitterion. It can be predicted from our calculations that in the transamination between alpha-amino acids and alpha-keto acids, the carbinolamine is formed with small barrier or even barrierless while the dehydration occurs easily at room temperature. However, without heating the 1,3 proton transfer could not occur as the barrier is 26.7 kcal/mol relative to the reactant complex when including two water molecules. Our results are in good agreement with experimental conclusions.

The potential energy surfaces of isomerization and dissociation reactions for CH2CHCOCl in the S0, T1, T2, and S1 states have been mapped with DFT, CASSCF, MP2, and MR-CI calculations. Rate constants for adiabatic and nonadiabatic processes have been calculated with the RRKM rate theory, in conjugation with the vibronic interaction method. Mechanistic photochemistry of CH2CHCOCl at 230-310 nm has been characterized through the computed potential energy surfaces and rate constants. Upon photoexcitation of CH2CHCOCl at 310 nm, the S1-->T1 intersystem crossing is the dominant primary process, which is followed by the 1,3-Cl migration along the T1 pathway. Meanwhile, the S1-->S0 internal conversion occurs with considerable probability and the subsequent trans-cis isomerization proceeds in the ground state. The C-Cl bond cleavage is an exclusive primary channel upon photoexcitation of gaseous CH2CHCOCl at 230 nm. The direct C-Cl bond cleavage is partially blocked by effects of the matrix, and the internal conversion from S1 to S0 becomes an important process for the excited molecule to deactivate in the condensed phase. The present calculations not only provide a reasonable explanation of the experimental findings, but also give new insight into the mechanistic photochemistry of CH2CHCOCl.

Photodissociation of aqueous formic acid has been investigated with the CASSCF, DFT, and MR-CI methods. Solvent effects are considered as a combination of the hydrogen-bonding interaction from explicit H2O molecules and the effects from the bulk surrounding H2O molecules using the polarizable continuum model. It is found that the hydrogen-bonding effect from the explicit water in the complex is the major factor to influence properties of aqueous formic acid, while the bulk surrounding H2O molecules has a noticeable influence on the structures of the complex. The direct C-O bond fission along the S1 pathway is predicted to be an important channel upon photolysis of aqueous formic acid at 200 nm, which is consistent with experimental observation that aqueous formic acid dissociates predominantly into fragments of HCO and OH. The existence of a dark channel upon photolysis of aqueous formic acid at 200 nm is assigned as fast relaxation from the S1 Franck-Condon geometry to the T1/S1 intersection and subsequent S1-->T1 intersystem crossing process. S1-->S0 internal conversion followed by molecular elimination to CO+H2O is the most probable primary process for formation of carbon monoxide, which was observed with considerable yield upon photolysis of aqueous formic acid at 253.7 nm.

All possible unimolecular processes upon photolysis of ClC(O)SCl in the UV-visible region have been characterized in the present paper through the optimized stationary structures and computed potential-energy profiles of the S0, S1, T2, and S2 states with the MP2, B3LYP, CASSCF, and MR-CI methods in conjugation with the cc-pVDZ basis set. Upon photoexcitation in the range of 300-400 nm, the ClC(O)SCl molecules are excited to the S1 state. From this state, the dissociation into ClC(O)S + Cl takes place immediately and subsequently Cl2 and SCO are formed. The C-Cl and C-S bond fissions that start from the S2 state are the dominant channels upon photodecomposition of ClC(O)SCl in the gas and condensed phases in the wavelength range of 200-248 nm. The formed Cl, C(O)SCl, ClCO, and SCl radicals are very reactive, and the Cl2, SCO, CO, and SCl2 molecules are subsequently produced as stable products in the condensed phase.

The potential energy surfaces of the C-O cleavage, rotational isomerization, keto-enolic tautomerization, and dehydration reactions of acetylacetone in the lowest triplet and ground states have been determined using the complete active space self-consistent field and density functional theory methods. The main photochemical mechanism obtained indicates that the acetylacetone molecule in the S(2)((1)pipi*) state can relax to the T(1)((3)pipi*) state via the S(2)-S(1) vibronic interaction and an S(1)/T(1)/T(2) intersection. The C-O fission pathway is the predominant dissociation process in the T(1)((3)pipi) state. Rotational isomerization reactions proceed difficultly in the ground state but very easily in the T(1)((3)pipi*) state. Keto-enolic tautomerization takes place with little probability for acetylacetone in the gas phase.

The potential energy surfaces of isomerization, dissociation, and elimination reactions for CH3CH2COCl in the S0 and S1 states have been mapped with the different ab initio calculations. Mechanistic photodissociation of CH3CH2COCl at 266 nm has been characterized through the computed potential energy surfaces, the optimized surface crossing structure, intrinsic reaction coordinate, and ab initio molecular dynamics calculations. Photoexcitation at 266 nm leads to the CH3CH2COCl molecules in the S1 state. From this state, the C-Cl bond cleavage proceeds in a time scale of picosecond in the gas phase. The barrier to the C-Cl bond cleavage on the S1 surface is significantly increased by effects of the matrix and the internal conversion to the ground state prevails in the condensed phase. The HCl eliminations as a result of internal conversion to the ground state become the dominant channel upon photodissociation of CH3CH2COCl in the argon matrix at 10 K.

The CASSCF and DFT methods have been used to determine geometric and electronic structures of the benzoic acid monomer in the S(0), S(1), S(2), T(1), and T(2) electronic states. The S(1)/T(2)/T(1) three-surface intersection was found by the state-averaged CASSCF calculations, which, in combination with features of the five lowest electronic states, provides new insights into photophysical processes of the benzoic acid monomer. The potential energy profiles of the alpha C-C and C-O bond fissions as well as decarboxylation reaction in different electronic states have been determined for the benzoic acid monomer. The alpha C-O bond cleavage starts from the T(2) state and leads to the fragments of C(6)H(5)CO(X(2)A') and OH(X(2)Pi) in the ground state, which is predicted to be the most possible channel upon photoexcitation of the benzoic acid monomer at 270 nm or shorter wavelengths.

The potential energy profiles for proton-transfer reactions of 2-hydroxypyridine and its complexes with water were determined by MP2, CASSCF and MR-CI calculations with the 6-31G** basis set. The tautomerization reaction between 2-hydroxypyridine (2HP) and 2-pyridone (2PY) does not take place at room temperature because of a barrier of approximately 35 kcal/mol for the ground-state pathway. The water-catalyzed enol-keto tautomerization reactions in the ground state proceed easily through the concerted proton transfer, especially for the two-water complex. The S1 tautomerization between the 2HP and 2PY monomers has a barrier of 18.4 kcal/mol, which is reduced to 5.6 kcal/mol for the one-water complex and 6.4 kcal/mol for the two-water complex. The results reported here predict that the photoinduced tautomerization reaction between the enol and keto forms involves a cyclic transition state having one or two water molecules as a bridge.

The potential energy surfaces for Cl(2)CS dissociation into ClCS + Cl in the five lowest electronic states have been determined with the combined complete active space self-consistent field (CASSCF) and MR-CI method. The wavelength-dependent photodissociation dynamics of Cl(2)CS have been characterized through computed potential energy surfaces, surface crossing points, and CASSCF molecular dynamics calculations. Irradiation of the Cl(2)CS molecules at 360-450 nm does not provide sufficient internal energy to overcome the barrier on S(1) dissociation, and the S(1)/T(2) intersection region is energetically inaccessible at this wavelength region; therefore, S(1) --> T(1) intersystem crossing is the dominant process, which is the main reason S(1)-S(0) fluorescence breaks off at excess energies of 3484-9284 cm(-1). Also, the S(1) --> T(2) intersystem crossing process can take place via the S(1)-T(2) vibronic interaction in this range of excess energies, which is mainly responsible for the quantum beats observed in the S(1) emission. Both S(2) direct dissociation and S(2) --> S(3) internal conversion are responsible for the abrupt breakoff of S(2)-S(0) fluorescence at higher excess energies. S(2) direct dissociation leads to the formation of the fragments of Cl(X(2)P) + ClCS(A(2)A' ') in excited electronic states, while S(2) --> S(3) internal conversion followed by direct internal conversion to the ground electronic state results in the fragments produced in the ground state.