Vibronic coupling in the ground and excited states of the imidazole radical cation.

Vibronic interactions in the ground and two excited states of the imidazole radical cation, X2A″ (π-1), A2A' (nσ-1), and B2A″ (π-1), and the associated nuclear dynamics were studied theoretically. The results were used to interpret the recent photoelectron measurements [M. Patanen et al., J. Chem. Phys. 155, 054304 (2021)]. The present high-level electronic structure calculations employing, in particular, the single, double, and triple excitations and equation-of-motion coupled-cluster method accounting for single and double excitation approaches and complete basis set extrapolation technique for the evaluation of the vertical ionization energies of imidazole indicate that the A 2A' and B 2A″ states are very close in energy and subject to non-adiabatic effects. Our modeling confirms the existence of pronounced vibronic coupling of the A 2A' and B 2A″ states. Moreover, despite the large energy gap of nearly 1.3 eV, the ground state X 2A″ is efficiently coupled to the A 2A' state. The modeling was performed within the framework of the three-state linear vibronic coupling problem employing Hamiltonians expressed in a basis of diabatic electronic states and parameters derived from ab initio calculations. The ionization spectrum was computed using the multi-configuration time-dependent Hartree method. The calculated spectrum is in good agreement with the experimental data, allowing for some interpretation of the observed features to be proposed.

Vibronic coupling in the ground and excited states of the pyridine radical cation.

Vibronic interactions in the pyridine radical cation ground state, 2A1, and its lowest excited states, 2A2 and 2B1, are studied theoretically. These states originate from the ionization out of the highest occupied orbitals of pyridine, 7a1 (nσ), 1a2 (π), and 2b1 (π), respectively, and give rise to the lowest two photoelectron maxima. According to our previous high-level ab initio calculations [Trofimov et al., J. Chem. Phys. 146, 244307 (2017)], the 2A2 (π-1) excited state is very close in energy to the 2A1 (nσ-1) ground state, which suggests that these states could be vibronically coupled. Our present calculations confirm that this is indeed the case. Moreover, the next higher excited state, 2B1 (π-1), is also involved in the vibronic interaction with the 2A1 (nσ-1) and 2A2 (π-1) states. The three-state vibronic coupling problem was treated within the framework of a linear vibronic coupling model employing parameters derived from the ionization energies of pyridine computed using the linear response coupled-cluster method accounting for single, double, and triple excitations (CC3). The potential energy surfaces of the 2A1 and 2A2 states intersect in the vicinity of the adiabatic minimum of the 2A2 state, while the surfaces of the 2A2 and 2B1 states intersect near the 2B1 state minimum. The spectrum computed using the multi-configuration time-dependent Hartree (MCTDH) method accounting for 24 normal modes is in good qualitative agreement with the experimental spectrum of pyridine obtained using high-resolution He I photoelectron spectroscopy and allows for some assignment of the observed features.

An experimental and theoretical study of the photoelectron spectra of cis-dichloroethene: Valence shell vertical ionization and vibronic coupling in the low-lying cationic states.

The valence shell photoelectron spectrum of cis-dichloroethene has been studied both experimentally and theoretically. Photoelectron spectra have been recorded with horizontally and vertically plane polarized synchrotron radiation, thereby allowing the anisotropy parameters, characterizing the angular distributions, to be determined. The third-order algebraic-diagrammatic construction approximation scheme for the one-particle Green's function has been employed to compute the complete valence shell ionization spectrum. In addition, the vertical ionization energies have been calculated using the outer valence Green's function (OVGF) method and the equation-of-motion coupled-cluster, with single and double substitutions for calculating ionization potentials (EOM-IP-CCSD) model. The theoretical results have enabled assignments to be proposed for most of the structure observed in the experimental spectra, including the inner-valence regions dominated by satellite states. The linear vibronic coupling model has been employed to study the vibrational structure of the lowest photoelectron bands, using parameters obtained from ab initio calculations. The ground state optimized geometries and vibrational frequencies have been computed at the level of the second-order Møller-Plesset perturbation theory, and the dependence of the ionization energies on the nuclear configuration has been evaluated using the OVGF method. While the adiabatic approximation holds for the X̃2B1 state photoelectron band, the Ã2B2, B̃2A1, and C̃2A2 states interact vibronically and form a complex photoelectron band system with four distinct maxima. The D̃2B1 and Ẽ2B2 states also interact vibronically with each other. The potential energy surface of the D̃2B1 state is predicted to have a double-minimum shape with respect to the out-of-plane a2 deformations of the molecular structure. The single photoelectron band resulting from this interaction is characterized by a highly irregular structure, reflecting the non-adiabatic nuclear dynamics occurring on the two coupled potential energy surfaces forming a conical intersection close to the minimum of the Ẽ2B2 state.

Few-femtosecond passage of conical intersections in the benzene cation.

Observing the crucial first few femtoseconds of photochemical reactions requires tools typically not available in the femtochemistry toolkit. Such dynamics are now within reach with the instruments provided by attosecond science. Here, we apply experimental and theoretical methods to assess the ultrafast nonadiabatic vibronic processes in a prototypical complex system-the excited benzene cation. We use few-femtosecond duration extreme ultraviolet and visible/near-infrared laser pulses to prepare and probe excited cationic states and observe two relaxation timescales of 11 ± 3 fs and 110 ± 20 fs. These are interpreted in terms of population transfer via two sequential conical intersections. The experimental results are quantitatively compared with state-of-the-art multi-configuration time-dependent Hartree calculations showing convincing agreement in the timescales. By characterising one of the fastest internal conversion processes studied to date, we enter an extreme regime of ultrafast molecular dynamics, paving the way to tracking and controlling purely electronic dynamics in complex molecules.

Ab Initio Benchmark Study of Nonadiabatic S1-S2 Photodynamics of cis- and trans-Hexatriene.

The dynamics of the nonadiabatically coupled lowest singlet excited states of cis- and trans-hexatriene are studied theoretically, in a comprehensive electronic structure and quantum dynamical investigation. At the ground state equilibrium geometry the relevant S2 and S1 states carry the A1 (Ag) and B2 (Bu) symmetry labels, for the cis (trans) isomer. Various high-level electronic structure methods are used, including the recently reparametrized DFT/MRCI method, and the results are critically compared. Key parameters of interest are the vertical energy gap and the strength of vibronic coupling between the interacting electronic states. To estimate their influence, suitable comparison calculations are performed. The results are used as the basis for quantum dynamical calculations on the UV absorption spectrum and electronic population transfer involving the S1 and S2 states. Up to nine nonseparable degrees of freedom are included in the calculations. The experimental UV absorption spectrum in the 5-5.2 eV energy range can be very well reproduced. The time-dependent wavepacket propagations reveal a population transfer on the order of 30-50 fs, which becomes increasingly complete with more degrees of freedom included in the calculation. The results are briefly compared with analogous data for the s-trans-butadiene system treated by some of us recently.

Carnosine treatment in combination with ACE inhibition in diabetic rats.

In humans, we reported an association of a certain allele of carnosinase gene with reduced carnosinase activity and absence of nephropathy in diabetic patients. CN1 degrades histidine dipeptides such as carnosine and anserine. Further, we and others showed that treatment with carnosine improves renal function and wound healing in diabetic mice and rats. We now investigated the effects of carnosine treatment alone and in combination with ACE inhibition, a clinically established nephroprotective drug in diabetic nephropathy. Male Sprague-Dawley rats were injected i.v. with streptozotocin (STZ) to induce diabetes. After 4 weeks, rats were unilaterally nephrectomized and randomized for 24 weeks of treatment with carnosine, lisinopril or both. Renal CN1 protein concentrations were increased under diabetic conditions which correlated with decreased anserine levels. Carnosine treatment normalized CN1 abundance and reduced glucosuria, blood concentrations of glycosylated hemoglobin (HbA1c), carboxyl-methyl lysine (CML), N-acetylglucosamine (GlcNac; all p<0.05 vs. non-treated STZ rats), reduced cataract formation (p<0.05) and urinary albumin excretion (p<0.05), preserved podocyte number (p<0.05) and normalized the increased renal tissue CN1 protein concentration. Treatment with lisinopril had no effect on HbA1C, glucosuria, cataract formation and CN1 concentration, but reduced albumin excretion rate more effectively than carnosine treatment (p<0.05). Treatment with both carnosine and lisinopril combined the effects of single treatment, albeit without additive effect on podocyte number or albuminuria. Increased CN1 amount resulted in decreased anserine levels in the kidney. Both carnosine and lisinopril exert distinct beneficial effects in a standard model of diabetic nephropathy. Both drugs administered together combine the respective effects of single treatment, albeit without exerting additive nephroprotection.

Reaction surface approach to multimode vibronic coupling problems: general framework and application to furan.

A new general framework for treating the dynamics on intersecting multidimensional potential energy surfaces is presented. It rests on a sub-division of the nuclear coordinates into different classes, one of primary importance with large-amplitude displacements during the process of interest and another one with smaller displacements, thus permitting a more approximate description. The latter are treated within the well-known linear + quadratic vibronic coupling scheme, where, however, the expansion "coefficients" are general functions of the "primary" coordinates. This may be augmented by an effective-mode approach for further degrees of freedom acting as an environment for the dynamics of the original modes. Following the general considerations, the approach is applied to the nonadiabatic photodynamics of furan and is shown to allow for an eight-dimensional quantum treatment, of higher dimension than was possible so far. The influence of the various degrees of freedom on the dynamics and lifetime of furan due to nonadiabatic ring-opening is discussed.

Ab initio quantum study of nonadiabatic S(1)-S(2) photodynamics of s-trans-butadiene.

The nonadiabatic photodynamics of s-trans-butadiene in its lowest singlet excited states is studied theoretically, using a fully quantal approach. The coupled 1Bu and 2Ag states are considered in the calculation, representing the lowest dipole-allowed electronic transition, and the dipole-forbidden state with substantial double-excitation character, respectively. Up to six nuclear degrees of freedom, including out-of-plane dihedral angles, are included. The calculation of the underlying potential energy surfaces relies on the CASPT2 method, where widely different CAS spaces have been compared. The ultrafast electronic population decay is confirmed, proceeding on a time scale of 30-40 fs. Pronounced out-of-plane distortions are obtained for the first time from fully quantal calculations. The complexity of the electronic absorption spectrum increases substantially upon including additional vibrational modes in the calculation. Further computations were performed to facilitate inclusion of the coupling to the ground state in subsequent work.

Quantum dynamics through conical intersections: combining effective modes and quadratic couplings.

We present a detailed study for the short-time dynamics through conical intersections in molecular systems related to the quadratic vibronic coupling (QVC) Hamiltonian [Müller, H.; Köppel, H.; Cederbaum, L. S. New J. Chem. 1993, 17, 7-29] and the effective-mode formalism [Cederbaum, L. S.; Gindensperger, E.; Burghardt, I. Phys. Rev. Lett. 2005, 94, 113003]. Our approach is based on splitting the nuclear degrees of freedom of the whole system into system modes and environment modes. It was found that only three-effective environmental modes together with the system's modes are needed to describe the short-time dynamics of the complex system correctly. In addition, a detailed mathematical proof is given in the appendix to demonstrate that the exact cumulants are recovered up to the second order within the cumulant expansion of the autocorrelation function. The butatriene molecule is studied as an explicit showcase example to stress the viability of our proposed scheme and to compare with other systems.

Ab initio quantum dynamical study of photoinduced ring opening in furan.

The nonadiabatic photoinduced ring opening occurring in the two lowest excited singlet states of furan is investigated theoretically, using wave-packet propagation techniques. The underlying multidimensional potential energy surfaces (PESs) are obtained from ab initio computations, using the equation-of-motion coupled cluster method restricted to single and double excitations (EOM-CCSD), reported in earlier recent work [E. V. Gromov, A. B. Trofimov, F. Gatti, and H. Köppel, J. Chem. Phys. 133, 164309 (2010)]. Up to five nuclear degrees of freedom are considered in the quantum dynamical treatment. Four of them represent in-plane motion for which the electronic states in question (correlating with the valence (1)B(2)(V) and Rydberg (1)A(2)(3s) states at the C(2v) ground-state molecular configuration) have different symmetries, A(') and A(''), respectively. The fifth mode, representing out-of-plane bending of the oxygen atom against the carbon-atom plane, leads to an interaction of these states, as is crucial for the photoreaction. The nonadiabatic coupling and conical intersection cause an electronic population transfer on the order of ∼10 fs. Its main features, and that of the wave-packet motion, are interpreted in terms of properties of the PES. The lifetime due to the ring-opening process has been estimated to be around 2 ps. The dependence of this estimate on the nuclear degrees of freedom retained in the computations is discussed.

Ab initio quantum dynamical study of the multi-state nonadiabatic photodissociation of pyrrole.

There has been a substantial amount of theoretical investigations on the photodynamics of pyrrole, often relying on surface hopping techniques or, if fully quantal, confining the study to the lowest two or three singlet states. In this study we extend ab initio based quantum dynamical investigations to cover simultaneously the lowest five singlet states, two π-σ∗ and two π-π∗ excited states. The underlying potential energy surfaces are obtained from large-scale MRCI ab initio computations. These are used to extract linear and quadratic vibronic coupling constants employing the corresponding coupling models. For the N-H stretching mode Q(24) an anharmonic treatment is necessary and also adopted. The results reveal a sub-picosecond internal conversion from the S(4) (π-π∗) state, corresponding to the strongly dipole-allowed transition, to the S(1) and S(2) (π-σ∗) states and, hence, to the ground state of pyrrole. The significance of the various vibrational modes and coupling terms is assessed. Results are also presented for the dissociation probabilities on the three lowest electronic states.

Theoretical study of photoinduced ring-opening in furan.

The potential energy surfaces (PESs) of the two lowest excited singlet states of furan [correlating with the Rydberg (1)A(2)(3s) and valence (1)B(2)(V) states at the C(2v) ground-state molecular configuration] have been studied in some detail with regard to the photoinduced ring-opening reaction. The surfaces have been characterized in terms of their stationary points and points of minimum energy conical intersections along the ring-opening pathway. The optimization of the geometrical parameters has been performed with the equation of motion coupled cluster singles and doubles method. The ab initio PESs have been modeled by energy grids and Taylor series. The resulting 11-dimensional PESs reproduce the ab initio results to a good accuracy and can be used in dynamical calculations.

Ab initio study of the VUV-induced multistate photodynamics of formaldehyde.

Although formaldehyde, H2CO, has been extensively studied there are still several issues not-well understood, specially regarding its dynamics in the VUV energy range, mainly due to the amount of nonadiabatic effects governing its dynamics. Most of the theoretical work on this molecule has focused on vertical excitation energies of Rydberg and valence states. In contrast to photodissociation processes involving the lowest-lying electronic states below 4.0 eV, there is little known about the photodynamics of the high-lying electronic states of formaldehyde (7-10 eV). One question of particular interest is why the (π, π*) electronic state is invisible experimentally even though it corresponds to a strongly dipole-allowed transition. In this work we present a coupled multisurface 2D photodynamics study of formaldehyde along the CO stretching and the symmetric HCH bending motion, using a quantum time-dependent approach. Potential energy curves along all the vibrational normal modes of formaldehyde have been computed using equation-of-motion coupled cluster including single and double excitations with a quadruply augmented basis set. In the case of the CO stretching coordinate, state-averaged complete active space self-consistent field followed by multireference configuration interaction was used for large values of this coordinate. 2D (for the CO stretching coordinate and the HCH angle) and 3D (including the out-of-plane distortion) potential energy surfaces have been computed for several Rydberg and valence states. Several conical intersections (crossings between potential energy surfaces of the same multiplicity) have been characterized and analyzed and a 2D 5 × 5 diabatic model Hamiltonian has been constructed. Based on this Hamiltonian, electronic absorption spectra, adiabatic and diabatic electronic populations and vibrational densities have been obtained and analyzed. The experimental VUV absorption spectrum in the 7-10 eV energy range is well reproduced, including the vibrational structure and the high irregularity in the regime of strong interaction between the (π, π*) electronic state and neighboring Rydberg states.

Ab initio potential energy surfaces and nonadiabatic interactions in the H+ +NO collision system.

Ab initio calculations on the H(+)+NO system have been carried out in Jacobi coordinates at the multireference configuration interaction level employing Dunning's correlation-consistent polarized valence triple zeta basis set to analyze the role of low-lying electronic excited states in influencing the collision dynamics relevant to the experimental collision energy range of 9.5-30 eV. The lowest two adiabatic potential energy surfaces, asymptotically correlating to H(+)+NO(X (2)Pi) and H((2)S)+NO(+)(X (1)Sigma(+)), have been obtained. Using ab initio procedures, the (radial) nonadiabatic couplings and the mixing angle between the lowest two electronic states (1 (2)A' and 2 (2)A') have been obtained to yield the corresponding quasidiabatic potential energy matrix. The strengths of the computed vibrational coupling matrix elements reflect a similar trend, as has been observed experimentally in the magnitudes of the state-to-state transition probability for the inelastic vibrational excitations [J. Krutein and F. Linder, J. Chem. Phys. 71, 559 (1979); F. A. Gianturco et al., J. Phys. B 14, 667 (1981)].

Multimode Jahn-Teller and pseudo-Jahn-Teller interactions in the cyclopropane radical cation: complex vibronic spectra and nonradiative decay dynamics.

The complex vibronic spectra and the nonradiative decay dynamics of the cyclopropane radical cation (CP+) are simulated theoretically with the aid of a time-dependent wave packet propagation approach using the multireference time-dependent Hartree scheme. The theoretical results are compared with the experimental photoelectron spectrum of cyclopropane. The ground and first excited electronic states of CP+ are of X2E' and A2E'' type, respectively. Each of these degenerate electronic states undergoes Jahn-Teller (JT) splitting when the radical cation is distorted along the degenerate vibrational modes of e' symmetry. The JT split components of these two electronic states can also undergo pseudo-Jahn-Teller (PJT)-type crossings via the vibrational modes of e'', a1'' and a2'' symmetries. These lead to the possibility of multiple multidimensional conical intersections and highly nonadiabatic nuclear motions in these coupled manifolds of electronic states. In a previous publication [J. Phys. Chem. A 2004, 108, 2256], we investigated the JT interactions alone in the X2E' ground electronic manifold of CP+. In the present work, the JT interactions in the A2E'' electronic manifold are treated, and our previous work is extended by considering the coupling between the X2E' and A2E'' electronic states of CP+. The nuclear dynamics in this coupled manifold of two JT split doubly degenerate electronic states is simulated by considering fourteen active and most relevant vibrational degrees of freedom. The vibronic level spectra and the ultrafast nonradiative decay of the excited cationic states are examined and are related to the highly complex entanglement of electronic and nuclear degrees of freedom in this prototypical molecular system.

The serotonin transporter gene polymorphism is not associated with smoking behavior.

Nicotine increases serotonin release in the brain and symptoms of nicotine withdrawal may be modulated by diminished serotonergic neurotransmission. The promoter region of the serotonin transporter gene, solute carrier family neurotransmitter transporter member 4 (SLC6A4), contains a functional tandem repeat polymorphism. The long (L) variant is more actively transcribed than the short (S) variant and is associated with a higher serotonin uptake. To investigate the potential role of this polymorphism for smoking behavior, SLC6A4 genotypes were determined in two different studies, the SMOKING GENES case-control study (470 current smokers and 419 subjects who had never smoked) and the cross-sectional Ludwigshafen risk and cardiovascular health (LURIC) study (777 current smokers and 1178 subjects who had never smoked). In the SMOKING GENES case-control study, SLC6A4 genotype frequencies were not statistically different between smokers (LL: 30.9%; LS: 46.8%; SS: 16.4%) and non-smokers (LL: 36.3%; LS: 41.8%; SS: 14.3%; P=0.13). Similar results were obtained in the cross-sectional LURIC study (smokers: LL, 36.5%, LS, 45.6%, SS, 17.9%; non-smokers: LL, 33.6%, LS, 48.9%, SS, 17.6%; P=0.33). SLC6A4 genotypes were furthermore not associated with Fagerstrom Tolerance Questionnaire score, packyears, number of cigarettes smoked per day or previous attempts to quit smoking. We conclude that the SLC6A4 promoter polymorphism is not a major determinant of smoking behavior in Caucasian.

On the (Emultiply sign in circlee)-Jahn-Teller conical intersections in the 3p(E') and 3d(E") Rydberg electronic states of triatomic hydrogen.

The static and dynamic aspects of the Jahn-Teller (JT) interactions in the 3p(E') and 3d(E") Rydberg electronic states of H3 are analyzed theoretically. The static aspects are discussed based on recent ab initio quantum chemistry results, and the dynamic aspects are examined in terms of the vibronic spectra and nonradiative decay behavior of these states. The adiabatic potential-energy surfaces of these degenerate electronic states are derived from extensive ab initio calculations. The calculated adiabatic potential-energy surfaces are diabatized following our earlier study on this system in its 2p(E') ground electronic state. The nuclear dynamics on the resulting conically intersecting manifold of electronic states is studied by a time-dependent wave-packet approach. Calculations are performed both for the uncoupled and coupled state situations in order to understand the importance of nonadiabatic interactions due to the JT conical intersections in these excited Rydberg electronic states.

Remarkable impact of intermode couplings on multimode vibronic dynamics: the photoelectron spectrum of CH3F.

Electronic and nuclear motions on intersecting potential energy surfaces are often intricately mixed and the spectrum can become very complex. Here we choose the strongly coupled Jahn-Teller system CH3F+ as a prototype example, and establish the importance of intermode coupling terms on multimode vibronic dynamics. The theoretical approach consists of a full second-order diabatic vibronic Hamiltonian, constructed from high-quality electronic structure calculations. Our results compare amazingly well with the experimental data. This highlights the success of the present theoretical approach in explaining the complex structure of vibronic spectra, ubiquitous in molecular systems.