First principles simulation of the UV absorption spectrum of ethylene using the vertical Franck-Condon approach

Environmental effects in computational spectroscopy: accuracy and interpretation

Spectroscopic techniques are valuable tools for understanding the structure and dynamics of complex systems, such as biomolecules or nanomaterials. Most of the current research is devoted to the development of new experimental techniques for improving the intrinsic resolution of different spectra. However, the subtle interplay of several different effects acting at different length and time scales still makes the interpretation and analysis of such spectra a very difficult task. In this respect, computational spectroscopy is becoming a needful and versatile tool for the assignment and interpretation of experimental spectra. It is in fact possible nowadays to model with relatively high accuracy the physical-chemical properties of complex molecules in different environments, and to link spectroscopic evidence directly to the structural and dynamical properties of optically or magnetically active solvated probes. In this Review, significant steps toward the simulation of entire spectra in condensed phases are presented together with some basic aspects of computational spectroscopy, which highlight how intramolecular and intermolecular degrees of freedom influence several spectroscopic parameters.

Combining the nuclear-electronic orbital approach with vibronic coupling theory: calculation of the tunneling splitting for malonaldehyde

The nuclear-electronic orbital (NEO) method is combined with vibronic coupling theory to calculate hydrogen tunneling splittings in polyatomic molecules. In this NEO-vibronic coupling approach, the transferring proton and all electrons are treated quantum mechanically at the NEO level, and the other nuclei are treated quantum mechanically using vibronic coupling theory. The dynamics of the molecule are described by a vibronic Hamiltonian in a diabatic basis of two localized nuclear-electronic states for the electrons and transferring proton. This ab initio approach is computationally practical and efficient for relatively large molecules, and the accuracy can be improved systematically. The NEO-vibronic coupling approach is used to calculate the hydrogen tunneling splitting for malonaldehyde. The calculated tunneling splitting of 24.5 cm(-1) is in excellent agreement with the experimental value of 21.6 cm(-1). This approach also enables the identification of the dominant modes coupled to the transferring hydrogen motion and provides insight into their roles in the hydrogen tunneling process.

Analysis and prediction of absorption band shapes, fluorescence band shapes, resonance Raman intensities, and excitation profiles using the time-dependent theory of electronic spectroscopy

A general method for the simulation of absorption (ABS) and fluorescence band shapes, resonance-Raman (rR) spectra, and excitation profiles based on the time-dependent theory of Heller is discussed. The following improvements to Heller's theory have been made: (a) derivation of new recurrence relations for the time-dependent wave packet overlap in the case of frequency changes between the ground and electronically excited states, (b) a new series expansion that gives insight into the nature of Savin's preresonance approximation, (c) incorporation of inhomogeneous broadening effects into the formalism at no additional computational cost, and (d) derivation of a new and simple short-time dynamics based equation for the Stokes shift that remains valid in the case of partially resolved vibrational structure. Our implementation of the time-dependent theory for the fitting of experimental spectra and the simulation of model spectra as well as the quantum mechanical calculation of the model parameters is discussed. The implementation covers all electronic structure approaches which are able to deliver ground- and excited-state energies and transition dipole moments. The technique becomes highly efficient if analytic gradients for the excited-state surface are available. In this case, the computational cost for the simultaneous prediction of ABS, fluorescence, and rR spectra is equal to that of a single excited-state geometry optimization step while the limitations of the short-time dynamics approximation are completely avoided. As a test case we discuss the well-known case of the strongly allowed 1 (1)A(g) --> 1 (1)B(u) transition in 1,3,5 trans-hexatriene in detail using method ranging from simple single-reference treatments to elaborate multireference electronic structure approaches. At the highest computational level, the computed spectra show the best agreement that has so far been obtained with quantum chemical methods for this problem.

A momentum-conserving Franck-Condon approximation: theory and application to the photodissociation of Li2+ in an intense laser field

A new, computationally efficient approximation to calculating matrix elements between vibrational-electronic states that does not assume stationary nuclei is introduced. The approach emphasizes the importance of the conservation of nuclear momenta. The calculated quantities can be used wherever equivalent quantities from standard Franck-Condon treatments can be used. The new method is illustrated by a time-dependent perturbation theory description of the photodissociation of the Li(2) (+) molecular cation in an intense laser field (I=0.9x10(12) W/cm(2)), which was previously identified in a detailed study of the electronic structure [Khait et al., J. Chem. Phys. 122, 094111 (2005)] as likely to have unusual sensitivity to the initial vibrational state in dynamics. The current study confirms this speculation and shows this to be a specific instance of a situation in which nuclear dynamics during an electronic optical transition cannot be ignored.

A theoretical study of the chiroptical properties of molecules with isotopically engendered chirality

There has been a considerable interest in the chiroptical properties of molecules whose chirality is exclusively due to an isotopic substitution and numerous examples for the electronic circular dichroism (CD) spectra of isotopically chiral systems have been reported in literature. Four different explanations have been proposed for the mechanism as to how the isotopic substitution induces a chiral perturbation of the otherwise achiral electronic wave function; however, up to now no conclusive answer has been given about the dominating effect responsible for the experimental observations. In this study we will present, for the first time, fully quantum-mechanical calculations of the CD spectra of three different molecular systems with isotopically engendered chirality. As examples, we consider the spectra of organic molecules with ketone and alpha-diketone carbonyl and diene chromophores. The effect of vibronic couplings for the reorientation of the electric and magnetic transition dipole moments is taken into account within the Herzberg-Teller approximation. The ground and excited state geometries and vibrational normal modes are obtained with (time-dependent) density functional theory [(TD)DFT], while the vibronic coupling effects are calculated at the TDDFT and density functional theory/multireference configuration interaction (DFT/MRCI) levels of theory. Generally, the band shapes of the experimental CD spectra are reproduced very well, and also the absolute CD intensities from the simulations are of the right order of magnitude. The sign and the intensity of the CD band are determined by a delicate balance of the contributions of a large number of individual vibronic transitions, and it is found that the vibrational normal modes with a large displacement are dominant. The separation of the calculated CD spectrum into the different contributions due to the overlap of the in-plane and out-of-plane components (regarding the symmetry plane of the unsubstituted molecule) of the electric and magnetic transition dipole moments yields information about the influence of the vibronic coupling effects for the reorientation of the corresponding transition dipole moments. In conclusion, the calculations clearly show that vibronic effects are responsible or at least dominant for the chiroptical properties of isotopically chiral organic molecules.

Theoretical study of the electronic structure of MCH2+(M=Fe,Co,Ni)

State of the art coupled cluster (CC) methods are applied to accurately characterize the ground state electronic structure and photoelectron spectra of transition metal carbene ions MCH(2) (+) (M=Fe, Co, and Ni). The geometries and energies of the lowest energy quartet, triplet, and doublet electronic states as well as several low-lying vertical excitation energies of FeCH(2) (+), CoCH(2) (+), and NiCH(2) (+) are reported. The excitation energies are computed using the equation-of-motion CC and for states of different symmetries, by the energy differences of single reference ground and excited states (Delta-CC). The latter employ several reference states; the unrestricted Hartree-Fock, restricted open shell Hartree-Fock, and unrestricted Kohn-Sham. We conclude that the (2)A(1) electronic ground state of NiCH(2) (+) is separated by about 30.0 kJ/mol from the next highest state, and the lowest (4)B(1) and (4)B(2) states of FeCH(2) (+) as well as the (3)A(2) and (3)A(1) states of CoCH(2) (+) are nearly degenerate. The presence of metal-pi*(MCH(2)) charge transfer states with significant oscillator strengths in the visible/near-UV energy domain of the theoretical spectra of FeCH(2) (+) and CoCH(2) (+) are at the origin of the photofragmentation of these compounds observed after irradiation between 310 and 360 nm.

Vibronic coupling in the excited cationic states of ethylene: simulation of the photoelectron spectrum between 12 and 18 eV

The effect of vibronic coupling on structure and spectroscopy is investigated in the excited cationic states of ethylene. It is found from equation of motion coupled cluster singles and doubles method for ionization potential electronic structure calculations in a triple-zeta plus double polarization basis set that ethylene in its third (B (2)A(g)) and fourth (C (2)B(2u)) ionized states does not have a stable minimum-energy geometry. The potential-energy surfaces of these states are energetically distinct and well separated at the ground-state geometry of ethylene, but in a geometry optimization as the structure of the ion relaxes, these surfaces end up in conical intersections and finally in the stable equilibrium geometry of the second ionized state (A (2)B(3g)). The topology of the potential-energy surfaces can be clearly understood using a vibronic model Hamiltonian. Furthermore, by diagonalizing this model Hamiltonian, the photoelectron spectrum of ethylene corresponding to the second, third, and fourth ionized states (12-18 eV) is simulated. Spectra from vibronic simulations including up to quartic coupling constants and using various normal-mode basis sets are compared to those from vertical Franck-Condon simulations to understand the importance of vibronic coupling and nonadiabatic effects and to examine the influence of individual normal modes on the spectrum.

The vibrational progressions of the N→V electronic transition of ethylene: A test case for the computation of Franck-Condon factors of highly flexible photoexcited molecules

The vibrational progressions of the N-->V electronic transition of ethylene--a test case for the computation of Franck-Condon factors between electronic states exhibiting very different equilibrium geometries--have been calculated by using both the Cartesian and the curvilinear internal coordinate representations of the normal modes of vibration. The comparison of the theoretical spectra with the experimental one shows that the Cartesian representation yields vibrational progressions which are not observed in the experimental spectrum, whereas the curvilinear one gives a very satisfying agreement, even in harmonic approximation.

Nonadiabatic effects in the photoelectron spectrum of the pyrazolide-d3 anion: Three-state interactions in the pyrazolyl-d3 radical

The 351.1 nm photoelectron spectrum of the 1-pyrazolide-d(3) anion has been measured. The photoelectron angular distributions indicate the presence of nearly degenerate electronic states of the 1-pyrazolyl-d(3) radical. Equation-of-motion ionization potential coupled-cluster singles and doubles (EOMIP-CCSD) calculations have been performed to study the low-lying electronic states. The calculations strongly suggest that three electronic states, energetically close to each other, are accessed in the photodetachment process. Strong interactions of the pseudo-Jahn-Teller type in each pair of the three states are evident in the calculations for the radical at the anion geometry. Model diabatic potentials of the three states have been constructed around the anion geometry in terms of the anion reduced normal coordinates up to the second order. An analytic method to parametrize the quadratic vibronic coupling (QVC) model potentials has been introduced. Parameters of the QVC model potentials have been determined from the EOMIP-CCSD and CCSD(T) calculations. Simulations of the 1-pyrazolide-d(3) spectrum have been performed with the model Hamiltonian, treating all vibronic interactions amongst the three states simultaneously. The simulation reproduces the fine structure of the observed spectrum very well, revealing complicated nonadiabatic effects in the low-lying states of the radical. The ground state of the 1-pyrazolyl-d(3) radical is (2)A(2) and the electron affinity is 2.935+/-0.006 eV. The first excited state is (2)B(1) with a term energy of 32+/-1 meV. While the high-symmetry (C(2v)) stationary points of the X (2)A(2) and A (2)B(1) states are minima, that of the state is a saddle point as a result of the pseudo-Jahn-Teller interactions with the other two states. The topology of the adiabatic potential energy surfaces is discussed.

Franck-Condon factors based on anharmonic vibrational wave functions of polyatomic molecules

Franck-Condon (FC) integrals of polyatomic molecules are computed on the basis of vibrational self-consistent-field (VSCF) or configuration-interaction (VCI) calculations capable of including vibrational anharmonicity to any desired extent (within certain molecular size limits). The anharmonic vibrational wave functions of the initial and final states are expanded unambiguously by harmonic oscillator basis functions of normal coordinates of the respective electronic states. The anharmonic FC integrals are then obtained as linear combinations of harmonic counterparts, which can, in turn, be evaluated by established techniques taking account of the Duschinsky rotations, geometry displacements, and frequency changes. Alternatively, anharmonic wave functions of both states are expanded by basis functions of just one electronic state, permitting the FC integral to be evaluated directly by the Gauss-Hermite quadrature used in the VSCF and VCI steps [Bowman et al., Mol. Phys. 104, 33 (2006)]. These methods in conjunction with the VCI and coupled-cluster with singles, doubles, and perturbative triples [CCSD(T)] method have predicted the peak positions and intensities of the vibrational manifold in the X 2B1 photoelectron band of H2O with quantitative accuracy. It has revealed that two weakly visible peaks are the result of intensity borrowing from nearby states through anharmonic couplings, an effect explained qualitatively by VSCF and quantitatively by VCI, but not by the harmonic approximation. The X 2B2 photoelectron band of H2CO is less accurately reproduced by this method, likely because of the inability of CCSD(T)/cc-pVTZ to describe the potential energy surface of open-shell H2CO+ with the same high accuracy as in H2O+.

Importance of vibronic effects on the circular dichroism spectrum of dimethyloxirane

We present a theoretical study on the vibrational structure of a circular dichroism (CD) spectrum using time-dependent density-functional theory in combination with a Franck-Condon-type approach. This method is applied to analyze the complex CD spectrum of dimethyloxirane, which involves delicate cancellations of positive and negative CD bands. Our approach reveals that these cancellations are strongly affected by the shapes of the CD bands, and that it is vital for an accurate simulation of the spectrum to take the different envelopes of these bands into account. One crucial point in some former theoretical studies on this compound, which were restricted to vertical excitations, was the appearance of a strong negative CD band in the energy range of 7.0-7.5 eV, which is not present in the experimental spectrum. We can explain the disappearance of this 2B band by a strong vibrational progression along normal modes with C-O stretching character, so that the band extends over an energy range of almost 1.1 eV. Thus, it overlaps with many other (mostly positive) CD bands, leading to a cancellation of its intensity. The dominant vibrational features in the experimental spectrum can be assigned to the 1B, 3B, and 5B bands, which show several clear vibrational peaks and a total bandwidth of only 0.3-0.5 eV. In order to obtain close agreement between the simulated and the experimental spectrum we have to apply small shifts to the vertical excitation energies that enter the calculation. These shifts account both for possible errors in the time-dependent density-functional theory calculations and for the neglect of differential zero-point energy between ground and excited states in our gradient-based vertical Franck-Condon approach.

Molecular applications of the intermediate Hamiltonian Fock-space coupled-cluster method for calculation of excitation energies

The intermediate Hamiltonian Fock-space coupled-cluster (FS-CC) method with singles and doubles is applied to calculate vertical excitation energies (EEs) for some molecular systems. The calculations are performed for several small molecules, such as H2O, N2, and CO, and for larger systems, such as C2H4, C4H6, and C6H6. Due to the intermediate Hamiltonian formulation, which provides a robust computational scheme for solving the FS-CC equations, and the efficient factorization strategy, relatively large basis sets and model spaces are employed permitting a comparison of the calculated vertical EEs with the experimental data.

Calculation of Franck–Condon factors including anharmonicity: Simulation of the C2H4+X̃B3u2←C2H4X̃Ag1 band in the photoelectron spectrum of ethylene

Our new simple method for calculating accurate Franck-Condon factors including nondiagonal (i.e., mode-mode) anharmonic coupling is used to simulate the C2H4+X2B3u<--C2H4X1A(g) band in the photoelectron spectrum. An improved vibrational basis set truncation algorithm, which permits very efficient computations, is employed. Because the torsional mode is highly anharmonic it is separated from the other modes and treated exactly. All other modes are treated through the second-order perturbation theory. The perturbation-theory corrections are significant and lead to a good agreement with experiment, although the separability assumption for torsion causes the C2D4 results to be not as good as those for C2H4. A variational formulation to overcome this circumstance, and deal with large anharmonicities in general, is suggested.

A Case Study of State-Specific and State-Averaged Brueckner Equation-of-Motion Coupled-Cluster Theory: The Ionic-Covalent Avoided Crossing in Lithium Fluoride

State-specific Brueckner equation-of-motion coupled-cluster theory (SS-B-EOMCC) is summarized, which can be considered an internally contracted version of a state-selective multireference coupled-cluster theory, which, however, is not entirely size-consistent. The method is applicable to general multireference problems, adheres to the space and spin symmetries of the molecular system, is straightforwardly extended to a state-averaged version, and has an associated perturbative variant which yields results close to the full coupled-cluster treatment. A key strength is that Brueckner orbitals are used, such that orbitals are optimized in the presence of dynamic correlation. A number of variations on the theme of SS-EOMCC is applied to study the ionic-covalent avoided crossing in LiF in a 6-311++G(3df,3pd) basis set. While reasonable results are obtained at the state-averaged level, the iterative solution process does not consistently converge for SS-EOMCC, due to the non-Hermiticity of the transformed Hamiltonian which may yield complex eigenvalues upon truncated diagonalization. This leads to an irrevocable breakdown of the state-specific EOMCC approach. We indicate some future directions that can resolve some of the problems with the SS-EOMCC methodology, as revealed by the demanding test case of the LiF potential energy curves.