Base Properties of H2CO in the Excited 1n→π* State

DFT + F12 QFFs for Cost-Effective Rovibrational Spectral Data Predictions of Ground and Excited Electronic States

The quest for faster computation of anharmonic vibrational frequencies of both ground and excited electronic states has led to combining coupled cluster theory harmonic force constants with density functional theory cubic and quartic force constants for defining a quartic force field (QFF) utilized in conjunction with vibrational perturbation theory at second order (VPT2). This work shows that explicitly correlated coupled cluster theory at the singles, doubles, and perturbative triples levels [CCSD(T)-F12] provides accurate anharmonic vibrational frequencies and rotational constants when conjoined with any of B3LYP, CAM-B3LYP, BHandHLYP, PBE0, and ωB97XD for roughly one-quarter of the computational time of the CCSD(T)-F12 QFF alone for our test set. As the number of atoms in the molecule increases, however, the anharmonic terms become a greater portion of the QFF, and the cost comparison improves with HOCO+ and formic acid, requiring less than 15 and 10% of the time, respectively. In electronically excited states, PBE0 produces more consistently accurate results. Additionally, as the size of the molecule and, in turn, QFF increase, the cost savings for utilizing such a hybrid approach for both ground- and excited-state computations grows. As such, these methods are promising for predicting accurate rovibrational spectral properties for electronically excited states. In cases where well-behaved potentials for a small selection of targeted excited states are needed, such an approach should reduce the computational cost compared to that of methods requiring semiglobal potential surfaces or variational treatments of the rovibronic Hamiltonian. Such applications include spectral characterization of comets, exoplanets, or any situation in which gas phase molecules are being excited by UV-vis radiation.

Blue-Shifted A−H Stretching Modes and Cooperative Hydrogen Bonding. 1. Complexes of Substituted Formaldehyde with Cyclic Hydrogen Fluoride and Water Clusters

The structures and vibrational spectra of the intermolecular complexes formed by insertion of substituted formaldehyde molecules HRCO (R = H, Li, F, Cl) into cyclic hydrogen fluoride and water clusters are studied at the MP2/aug-cc-pVTZ computational level. Depending on the nature of the substituent R, the cluster type, and its size, the C-H stretching modes of HRCO undergo large blue and partly red shifts, whereas all the F-H and O-H stretching modes of the conventional hydrogen bonds are strongly red-shifted. It is shown that (i) the mechanism of blue shifting can be explained within the concept of the negative intramolecular coupling between C-H and C=O bonds that is inherent to the HRCO monomers, (ii) the blue shifts also occur even if no hydrogen bond is formed, and (iii) variation of the acceptor X or the strength of the C-H...X hydrogen bond may either amplify the blue shift or cause a transition from blue shift to red shift. These findings are illustrated by means of intra- and intermolecular scans of the potential energy surfaces. The performance of the negative intramolecular coupling between C-H and C=O bonds of H(2)CO is interpreted in terms of the NBO analysis of the isolated H(2)CO molecule and H(2)CO interacting with (H2O)n and (HF)n clusters.

C2H+H2CO: a new route for formaldehyde removal

The title unknown reaction is theoretically studied at various levels to probe the interaction mechanism between the ethynyl radical (HC triple bond C) and formaldehyde (H(2)C double bond O). The most feasible pathway is a barrier-free direct H-abstraction process leading to acetylene and formyl radical (C(2)H(2)+HCO) via a weakly bound complex, and then the product can take secondary dissociation to the final product C(2)H(2)+CO+H. The C-addition channel leading to propynal plus H-atom (HCCCHO+H) has the barrier of only 3.6, 2.9, and 2.1 kcal/mol at the CCSD(T)/6-311+G(3df,2p)MP2//6-311G(d,p)+ZPVE, CCSD(T)/6-311+G(3df,2p)//QCISD/6-311G(d,p)+ZPVE, and G3//MP2 levels, respectively [CCSD(T)--coupled cluster with single, double, and triple excitations; ZPVE--zero-point vibrational energy; QCISD--quadratic configuration interaction with single and double excitations; G3//MP2-Gaussian-3 based on Moller-Plesset geometry]. The O addition also leading to propynal plus H atom needs to overcome a higher barrier of 5.3, 8.7, and 3.0 kcalmol at the three corresponding levels. The title no-barrier reaction presents a new efficient route to remove the pollutant H(2)CO, and should be included in the combustion models of hydrocarbons. It may also represent the fastest radical-H(2)CO reaction among the available theoretical data. Moreover, it could play an important role in the interstellar chemistry where the zero- or minute-barrier reactions are generally favored. Discussions are also made on the possible formation of the intriguing propynal in space via the title reaction on ice surface.

The Lowest Singlet (n,π*) and (π,π*) Excited States of the Hydrogen-Bonded Complex between Water and Pyrazine

The hydrogen bonding between water and pyrazine in its ground, lowest (n,pi*), and lowest (pi,pi*) states is investigated using density-functional theory (DFT), time-dependent density function theory (TD-DFT), coupled-cluster singles and doubles (CCSD) theory and equation-of-motion coupled cluster (EOM-CCSD) theory. For all states, the minimum-energy configuration is found to be an orthodox linear hydrogen-bonded species, with the bond strength increasing by 0.4 kcal mol-1 upon formation of the (pi,pi*) state and decreasing by 1.0 kcal mol-1 upon formation of the (n,pi*) state. The calculated solvent shifts for the complexes match experimental data and provide a basis for the understanding of the aqueous solvation of pyrazine, and the excited-state complexes are predicted to be only short-lived, explaining the failure of molecular beam experiments to observe them. Quite a different scenario for hydrogen bonding to the (n,pi*) excited state is found compared to those of H2O:pyridine and H2O:pyrimidine: for pyridine linear hydrogen bonds are unstable and hydrogen bonds to the electron-enriched pi cloud are strong, whereas for pyrimidine the excitation localizes on the nonbonded nitrogen leaving the hydrogen-bonding unaffected. For H2O:pyrazine, the (n,pi*) excitation remains largely delocalized, providing a distinct intermediary scenario.

Strongly Blue-Shifted C−H Stretches: Interaction of Formaldehyde with Hydrogen Fluoride Clusters

The equilibrium structures, binding energies, and vibrational spectra of the cyclic, hydrogen-bonded complexes formed between formaldehyde, H(2)CO, and hydrogen fluoride clusters, (HF)(1< or =n < or =4), are investigated by means of large-scale second-order Møller-Plesset calculations with extended basis sets. All studied complexes exhibit marked blue shifts of the C-H stretching frequencies, exceeding 100 cm(-1) for n = 2-4. It is shown that these blue shifts are, however, only to a minor part caused by blue-shifting hydrogen bonding via C-H...F contacts. The major part arises due to the structural relaxation of the H(2)CO molecule under the formation of a strong C=O...H-F hydrogen bond which strengthens as n increases. The close correlation between the different structural parameters in the studied series of complexes is demonstrated, and the consequences for the frequency shifts in the complexes are pointed out, corroborating thus the suggestion of the primary role of the C=O...H-F hydrogen bonding for the C-H stretching frequency shifts. This particular behavior, that the appearance of an increasingly stronger blue shift of the C-H stretching frequencies is mainly induced by the formation of a progressively stronger C=O...H-F hydrogen bond in the series of H(2)CO...(HF)(1< or =n < or =4), complexes and only to a lesser degree by the formation of the so-called blue-shifting C-H...F hydrogen bond, is rationalized with the aid of selected sections of the intramolecular H(2)CO potential energy surface and by performing a variety of structural optimizations of the H(2)CO molecule embedded in external, differently oriented dipole electric fields, and also by invoking a simple analytical force-field model.

First Singlet (n,π*) Excited State of Hydrogen-Bonded Complexes between Water and Pyrimidine

Hydrogen bonds from water to excited-state formaldehyde and from water to excited-state pyridine have been shown to display novel motifs to traditional hydrogen bonds involving ground states, with, in particular for H2O:pyridine, strong interactions involving the electron-rich pi cloud dominating the (n,pi) excited state. We investigate H2O:pyrimidine and various dihydrated species and reveal another motif, one in which the hydrogen bonding can dramatically alter the electronic structure of the excited state. Such effects are rare for ground-state interactions for which hydrogen bonding usually acts to merely perturb the electronic structure of the participating molecules. It arises as the (n,pi*) excitation of isolated pyrimidine is delocalized over both nitrogens but asymmetric hydrogen bonding causes it to localize on just the noninteracting atom. As a result, the excited-state hydrogen bond in H2O:pyrimidine is suprisingly very similar to the ground-state structure. These results lead to an improved understanding of the spectroscopy of pyrimidine in liquid water, and to the prediction that stable excited-state hydrogen bonds in H2O:pyrimidine should be observable, despite failure of experiments to actually do so. They also provide a simple model for the intricate control over primary charge separation in photosynthesis exerted by hydrogen bonding, and for solvent-induced electron localization in symmetric mixed-valence complexes. All conclusions are based on strong parallels found between the results of calculations performed using density-functional theory (DFT) and time-dependent DFT (TDDFT), complete-active-space self-consistent-field (CASSCF) with second-order perturbation-theory correction (CASPT2) theory, and equation-of-motion coupled cluster (EOM-CCSD) theory, calculations that are verified through detailed comparison of computed properties with experimental data for both the isolated molecules and the ground-state hydrogen bond.

Comparison of CIS- and EOM-CCSD-Calculated Adiabatic Excited-State Structures. Changes in Charge Density on Going to Adiabatic Excited States

The CIS and EOM-CCSD adiabatic geometries for the first excited states of a set of small molecules (C2H4, C2H2, H2C=O, H2C=S, CS2, CO2, SO2, NO2) have been calculated using the 6-311++G** basis set to see if the former geometries can be good starting points for optimizations at the latter theoretical level. With most of the molecules, there is fairly good agreement between the results from the two methods, and EOM-CCSD gives good agreement with the available experimental data. A detailed discussion of the lowest-lying singlet excited states in CO2 and CS2 is presented, highlighting the pronounced differences in electronic character and equilibrium structure displayed by these isovalent species. The origins of the structural distortions that are frequently found for the adiabatic excited states are examined with the aid of deformation density plots and the electron localization function (ELF).

Hydrogen-Bonding-Induced Shifts of the Excitation Energies in Nucleic Acid Bases: An Interplay between Electrostatic and Electron Density Overlap Effects

The theoretically calculated dimerization-induced shifts of the lowest excitation energies in two model systems, adenine-thymine and guanine-cytosine base pairs, are analyzed. The applied formalism is based on first principles and allows one to study the influence of the microscopic environment of a given molecule on its ground- [Wesolowski, T. A.; Warshel, A. J. Phys. Chem. 1993, 97, 8050] and excited-state [Casida, M. E.; Wesolowski, T. A. Int. J. Quantum Chem. 2004, 96, 577] properties. The assessment of the relative importance of such effects as (a) Coulomb interactions, (b) orbital interactions, (c) electronic polarization of the environment, and (d) electron density overlap effects is straightforward in this formalism. In the applied formalism, electron density overlap effects can be further decomposed into the exchange-correlation component which provides a small attractive contribution and the repulsive kinetic energy-dependent component. It is shown that the shifts can be attributed to the electrostatic interactions and the repulsive overlap-dependent term in the embedding potential. The electronic polarization of the environment plays a significant role (up to 30% of the total shift) only in transitions involving the orbitals localized on hydrogen bond donor groups. For all analyzed shifts, the contribution of the intermolecular orbital interactions is negligible. The analysis of this work provides strong evidence supporting the use of the widely applied embedding-molecule strategy in computational studies of chromophores in a condensed phase even in such cases where only one end of the hydrogen bond is included in the quantum mechanical part.

Higher-order equation-of-motion coupled-cluster methods

The equation-of-motion coupled-cluster (EOM-CC) methods truncated after double, triple, or quadruple cluster and linear excitation operators (EOM-CCSD, EOM-CCSDT, and EOM-CCSDTQ) have been derived and implemented into parallel execution programs. They compute excitation energies, excited-state dipole moments, and transition moments of closed- and open-shell systems, taking advantage of spin, spatial (real Abelian), and permutation symmetries simultaneously and fully (within the spin-orbital formalisms). The related Lambda equation solvers for coupled-cluster (CC) methods through and up to connected quadruple excitation (CCSD, CCSDT, and CCSDTQ) have also been developed. These developments have been achieved, by virtue of the algebraic and symbolic manipulation program that automated the formula derivation and implementation altogether. The EOM-CC methods and CC Lambda equations introduce a class of second quantized ansatz with a de-excitation operator (Y), a number of excitation operators (X), and a physical (e.g., the Hamiltonian) operator (A), leading to the tensor contraction expressions that can be performed in the order of ((...((yx)x)...)x)a or ((...((ax)x)...)x)y at the minimal peak operation cost, where x, y, and a are basis-set representations (i.e., tensors) of the respective operators X, Y, and A. Any intermediate tensor resulting from either contraction order is shown to have at most six groups of permutable indices, and this knowledge is used to guide the computer-synthesized programs to fully exploit the permutation symmetry of any tensor to minimize the arithmetic and memory costs.