Decoupling conditions for the vibronic equation in dimers

The analytical conditions are obtained and discussed under which the vibronic equations for dimers can be decoupled.     

An unusually large nonadiabatic error in the BNB molecule

The vibronic coupling model of Ko?uppel, Domcke, and Cederbaum in one dimension is introduced as a means to estimate the effects of electronic nonadiabaticity on the vibrational energy levels of molecules that exhibit vibronic coupling. For the BNB molecule, the nonadiabatic contribution to the nominal fundamental vibrational energy of the antisymmetric stretching mode is approximately -80?cm(-1). The surprisingly large effect for this mode, which corresponds to an adiabatic potential that is essentially flat near the minimum due to the vibronic interaction, is contrasted with another model system that also exhibits a flat potential (precisely, a vanishing quadratic force constant) but has a significantly larger gap between interacting electronic states. For the latter case, the nonadiabatic contribution to the level energies is about two orders of magnitude smaller even though the effect on the potential is qualitatively identical. A simple analysis shows that significant nonadiabatic corrections to energy levels should occur only when the affected vibrational frequency is large enough to be of comparable magnitude to the energy gap involved in the coupling. The results provide evidence that nonadiabatic corrections should be given as much weight as issues such as high-level electron correlation, relativistic corrections, etc., in quantum chemical calculations of energy levels for radicals with close-lying and strongly coupled electronic states even in cases where conical intersections are not obviously involved. The same can be said for high-accuracy thermochemical studies, as the zero-point vibrational energy of the BNB example contains a nonadiabatic contribution of approximately -70?cm(-1) (-0.9?kJ mol(-1)).     

Nonadiabatic scattering of NO off Au3 clusters: A simple and robust diabatic state manifold generation method for multiconfigurational wavefunctions

The presence of nonadiabatic effects during the interaction of small molecules with metals has been observed experimentally for the last decades. Specially remarkable are the effects found for NO/Au, where experiments have suggested the presence of very strong vibronic coupling during the molecular scattering. However, the accurate inclusion of the nonadiabatic effects in periodic boundary conditions (PBC) theoretical methods remain an unapproachable challenge. Here, aiming to give some theoretical insight to the strong vibronic coupling, we have adopted a pragmatic point of view, taking use of an auxiliary simplified system, NO/Au₃. We show the importance of nonadiabatic coupling, during the scattering of NO from a Au₃ cluster, using a diabatic representation of 12 electronic states of the system, including a few charge-transfer states. Our diabatic representation is obtained by rotating the orbital and configuration interaction (CI) vectors of a restricted active space (RAS) wavefunction. We present a strategy for extracting the best effective manifold of states relevant to the system, below some prescribed energy, directly from the RAS CI vectors. This scheme is able to disentangle a large dense manifold of adiabatic states with strong coupling and crossings. This approach is also shown to work for multireference configuration interaction (MRCI). By performing quantum propagations, we observed an increase in vibrational redistribution with increasing initial vibrational or translational energies. We suggest that these nonadiabatic effects should also be present at smaller energies in larger clusters. © 2018 Wiley Periodicals, Inc.     

Extremely narrow peaks in predissociation of sodium dimer due to rovibronic coupling

In sodium dimer the 2 (3)Pi(g), 3 (3)Pi(g), and 4 (3)Sigma(g) (+) electronic states are coupled; the coupling of the two (3)Pi(g) states is due to vibrational motion while the nonadiabatic interaction between the (3)Sigma(g) (+) and the (3)Pi(g) states-in particular, the 3 (3)Pi(g) state-is mediated by rotational interaction. The resulting vibronic problem is studied in some detail. The bound vibrational states of the 3 (3)Pi(g) and 4 (3)Pi(g) (+) states lie in the dissociation continuum of the 2 (3)Pi(g) state and become resonances due to the prevailing nonadiabatic coupling. The resonances are calculated using the complex scaling method and the available ab initio adiabatic potential energy curves. It is demonstrated that the resonances associated with rotational nonadiabatic coupling are narrower by several orders of magnitude than those that emerge from the vibrational nonadiabatic coupling. The predissociation cross section is computed and compared with experiment.     

Analytical determination of non-born—oppenheimer matrix elements between different molecular electronic states

Matrix elements for the nonadiabatic correction between different adiabatic Born—Oppenheimer (ABO) states are discussed in terms of a generating function (nonadiabatic coupling function, NAF) which was introduced earlier. Recursion relations are derived for the case where different vibronic states are coupled via a single non-totally symmetric mode. The formalism is applied to the estimation of non-Born—Oppenheimer coupling between vibronic states of pyrazine.     

Reduced dimensionality spin-orbit dynamics of CH3 + HCl ⇌ CH4 + Cl on ab initio surfaces

A reduced dimensionality quantum scattering method is extended to the study of spin-orbit nonadiabatic transitions in the CH(3) + HCl ? CH(4) + Cl((2)P(J)) reaction. Three two-dimensional potential energy surfaces are developed by fitting a 29 parameter double-Morse function to CCSD(T)/IB//MP2/cc-pV(T+d)Z-dk ab initio data; interaction between surfaces is described by geometry-dependent spin-orbit coupling functions fit to MCSCF/cc-pV(T+d)Z-dk ab initio data. Spectator modes are treated adiabatically via inclusion of curvilinear projected frequencies. The total scattering wave function is expanded in a vibronic basis set and close-coupled equations are solved via R-matrix propagation. Ground state thermal rate constants for forward and reverse reactions agree well with experiment. Multi-surface reaction probabilities, integral cross sections, and initial-state selected branching ratios all highlight the importance of vibrational energy in mediating nonadiabatic transition. Electronically excited state dynamics are seen to play a small but significant role as consistent with experimental conclusions.     

Nonadiabatic coupling and diabatic electronic population dynamics on 11A2 and 11B1 states of ozone molecule

In this work, we examine nonadiabatic population dynamics for 1¹B₁ and 1¹A₂ states of ozone molecule (O₃). In O₃, two lowest singlet excited states, ¹A₂ and ¹B₁, can be coupled. Thus, population transfer between them occurs through the seam involving these two states. At any point of the seam (conical intersection), the Born-Oppenheimer approximation breaks down, and it is necessary to investigate nonadiabatic dynamics. We consider a linear vibronic coupling Hamiltonian model and evaluate vibronic coupling constant, diabatic frequencies for three modes of O₃, bilinear and quadratic coupling constants for diabatic potentials, displacements, and Huang-Rhys coupling constants using ab initio calculations. The electronic structure calculations have been performed at the multireference configuration interaction and complete active space with second-order perturbation theory with a full-valence complete active space self-consistent field methods and augmented Dunning's standard correlation-consistent-polarized quadruple zeta basis set to determine ab initio potential energy surfaces for the ground state and first two excited states of O₃, respectively. We have chosen active space comprising 18 electrons distributed over 12 active orbitals. Our calculations predict the linear vibronic coupling constant 0.123 eV. We have obtained the population on the 1¹B₁ and 1¹A₂ excited electronic states for the first 500 fs after photoexcitation.     

High resolution vacuum ultraviolet emission spectrum of D(2) from 78 to 103 nm: The D (1)Pi(u)-->X (1)Sigma(g) (+) and D(') (1)Pi(u) (-)-->X (1)Sigma(g) (+) band systems

The emission spectrum of the D(2) molecule has been studied at high resolution in the vacuum ultraviolet region 78.5-102.7 nm. A detailed analysis of the two D (1)Pi(u)-->X (1)Sigma(g) (+) and D(') (1)Pi(u) (-)-->X (1)Sigma(g) (+) electronic band systems is reported. New and improved values of the level energies of the two upper states have been derived with the help of the program IDEN [V. I. Azarov, Phys. Scr. 44, 528 (1991); 48, 656 (1993)], originally developed for atomic spectral analysis. A detailed comparison is made between the observed energy levels and solutions of coupled equations using the newest ab initio potentials by Wolniewicz and co-workers [J. Chem. Phys. 103, 1792 (1995); 99, 1851 (1993); J. Mol. Spectros. 212, 208 (2002); 220, 45 (2003)] taking into account the nonadiabatic coupling terms for the D (1)Pi(u) state with the lowest electronic states B (1)Sigma(u) (+), C (1)Pi(u), and B(') (1)Sigma(u) (+). A satisfactory agreement has been found for most of the level energies belonging to the D and D(') states. The remaining differences between observation and theory are probably due to nonadiabatic couplings with other higher electronic states which were neglected in the calculations.     

Analysis of Diffusion in Hollow Geometries

The diffusion in hollow particles of solid adsorbent materials was analyzed based on analytical solutions to the basic diffusion equation. Three geometric shapes (plane sheet, cylinder, and sphere) of sorbent material were considered for two kinds of boundary conditions. The equations for determining the equivalent sizes compared to their corresponding solid particles were obtained directly from the theoretical expressions of sorption uptake curves. Among the three hollow particles of impermeable inner surface, the sphere gives the highest gain in effective diffusion rate compared to the corresponding solid particle. For permeable inner surface, at lower hollow volume fractions, the plane sheet shows the highest gain, while at higher hollow volume fractions, the sphere shows the highest gain in effective diffusion rate.     

Linearized non-equilibrium and non-isothermal two-dimensional model of liquid chromatography for studying thermal effects in cylindrical columns

A non-equilibrium and non-isothermal two-dimensional lumped kinetic model (2?D-LKM) is formulated and analytically solved to study the influence of temperature variations along the axial and radial coordinates of a liquid chromatographic column. The model includes convection-diffusion partial differential equations for mass and energy balances in the mobile phase coupled with differential equations for mass and energy in the stationary phase. The solutions are derived analytically through sequential implementation of finite Hankel and Laplace transformations using the Dirichlet inlet boundary conditions. The coupling between the thermal waves and concentration fronts is demonstrated through numerical simulations and important parameters are recognized that influence the column performance. For a more comprehensive study of the considered model, numerical temporal moments are obtained from the derived solutions. Several case studies are conducted and validity ranges of the derived analytical solutions are identified. The current analytical results will play a major role in the improvements of non-equilibrium and non-isothermal liquid chromatographic processes.     

An investigation of nonadiabatic interactions in Cl(2Pj) + D2 via crossed-molecular-beam scattering

We have determined limits on the cross section for both electronically nonadiabatic excitation and quenching in the Cl((2)P(j)) + D(2) system. Our experiment incorporates crossed-molecular-beam scattering with state-selective Cl((2)P(12,32)) detection and velocity-mapped ion imaging. By colliding atomic chlorine with D(2), we address the propensity for collisions that result in a change of the spin-orbit level of atomic chlorine either through electronically nonadiabatic spin-orbit excitation Cl((2)P(32)) + D(2)-->Cl(*)((2)P(12)) + D(2) or through electronically nonadiabatic spin-orbit quenching Cl(*)((2)P(12)) + D(2)-->Cl((2)P(32)) + D(2). In the first part of this report, we estimate an upper limit for the electronically nonadiabatic spin-orbit excitation cross section at a collision energy of 5.3 kcal/mol, which lies above the energy of the reaction barrier (4.9 kcal/mol). Our analysis and simulation of the experimental data determine an upper limit for the excitation cross section as sigma(NA)< or =0.012 A(2). In the second part of this paper we investigate the propensity for electronically nonadiabatic spin-orbit quenching of Cl(*) following a collision with D(2) or He. We perform these experiments at collision energies above and below the energy of the reaction barrier. By comparing the amount of scattered Cl(*) in our images to the amount of Cl(*) lost from the atomic beam we obtain the maximum cross section for electronically nonadiabatic quenching as sigma(NA)< or =15(-15) (+44) A(2) for a collision energy of 7.6 kcal/mol. Our experiments show the probability for electronically nonadiabatic quenching in Cl(*) + D(2) to be indistinguishable to that for the kinematically identical system of Cl(*) + He.     

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.     

On quantifying the vibronic interaction between electronic states

An accurate estimation of the interstate vibronic coupling strength is of particular relevance for the treatment of nonadiabatic dynamics. This is not a trivial task because direct interactions between electronic states have to be separated from intrinsic frequency shifts. Surprisingly, this issue has not been discussed in detail in the literature so far. An analysis of the error dependence is given for two formulas derived from linear vibronic coupling theory. The difficulty in estimating the interstate coupling parameters is shown to originate from the initially unknown contribution of the diagonal quadratic coupling coefficients to the total vibronic coupling. An interpretation of the error analysis including a numerical case study is followed by a more general discussion of the different mechanisms that can shape adiabatic electronic potential energy functions. Qualitative criteria are formulated for the differentiation between interstate and intrastate vibronic coupling effects based on electronic structure information. These ideas are then applied to investigate vibronic coupling problems in pyrazine as well as trans- and cis-hexatriene.     

Jahn-Teller effect in van der Waals complexes; Ar-C6H6 + and Ar-C6D6 +

The two asymptotically degenerate potential energy surfaces of argon interacting with the X (2)E(1g) ground state benzene(+) cation were calculated ab initio from the interaction energy of the neutral Ar-benzene complex given by Koch et al. [J. Chem. Phys. 111, 198 (1999)] and the difference of the geometry-dependent ionization energies of the complex and the benzene monomer computed by the outer valence Green's function method. Coinciding minima in the two potential surfaces of the ionic complex occur for Ar on the C(6v) symmetry axis of benzene(+) (the z axis) at z(e)=3.506 A. The binding energy D(e) of 520 cm(-1) is only 34% larger than the value for the neutral Ar-benzene complex. The higher one of the two surfaces is similar in shape to the neutral Ar-benzene potential, the lower potential is much flatter in the (x,y) bend direction. Nonadiabatic (Jahn-Teller) coupling was taken into account by transformation of the two adiabatic potentials to a two-by-two matrix of diabatic potentials. This transformation is based on the assumption that the adiabatic states of the Ar-benzene(+) complex geometrically follow the Ar atom. Ab initio calculations of the nonadiabatic coupling matrix element between the adiabatic states with the two-state-averaged CAS-SCF(5,6) method confirmed the validity of this assumption. The bound vibronic states of both Ar-C(6)H(6) (+) and Ar-C(6)D(6) (+) were computed with this two-state diabatic model in a basis of three-dimensional harmonic oscillator functions for the van der Waals modes. The binding energy D(0)=480 cm(-1) of the perdeuterated complex agrees well with the experimental upper bound of 485 cm(-1). The ground and excited vibronic levels and wave functions were used, with a simple model dipole function, to generate a theoretical far-infrared spectrum. Strong absorption lines were found at 10.1 cm(-1) (bend) and 47.9 cm(-1) (stretch) that agree well with measurements. The unusually low bend frequency is related to the flatness of the lower adiabatic potential in the (x,y) direction. The van der Waals bend mode of e(1) symmetry is quadratically Jahn-Teller active and shows a large splitting, with vibronic levels of A(1), E(2), and A(2) symmetry at 1.3, 10.1, and 50.2 cm(-1). The level at 1.3 cm(-1) leads to a strong absorption line as well, which could not be measured because it is too close to the monomer line. The level at 50.2 cm(-1) gives rise to weaker absorption. Several other weak lines in the frequency range of 10 to 60 cm(-1) were found.     

Jahn-Teller effect in tetraclusters of transition metals,taking into account vibronic coupling with E-vibrations

Conclusion From the results obtained it follows that the shape of the surface for the lower sheet of the adiabatic potential for the cluster in the E vibration space is similar to the surface for the E-e problem both in the linear and in the quadratic cases. Taking into account the extra modes of the E vibration which were discarded in this treatment does not qualitatively change the shape of the adiabatic potential surface (the number of minima and saddle points), which is a consequence of its symmetry [4]. The deformation pattern of the cluster corresponding to movement of the system along the bottom of the channel of one of the lower sheets is analogous to the case of the E-e problem for an octahedral and tetrahedral environment of a Jahn-Teller ion [10]. The average values of the pseudospin operators in the vibronic ground states are the same for all centers, which is evidence (taking into account the equality of the Jahn-Teller distortions at all centers) for ferrodistortion ordering of all the distortions, both in the linear and in the quadratic cases. If each vibronic center has spin s=1/2, all vibronic levels of the cluster have additional 16-fold spin degeneracy. Taking into account the effect of spin-orbit coupling on the shape of the adiabatic potential using perturbation theory shows that removal of spin degeneracy for the lower sheet does not occur, that the adiabatic potential does not change its shape but rather is only reduced in energy by , where is the spin-orbit coupling constant.The high orbital degeneracy of sheets 2, 3, 4 (Fig. 2a, b) is not removed by taking into account the quadratic vibronic coupling and is accidental. The indicated orbital and spin accidental degeneracy of the sheets is removed by taking into account the intercenter interaction.Translated from Teoreticheskaya i Éksperimental' naya Khimiya, Vol. 20, No. 1, pp. 1–9, January–February, 1984.     

Nonadiabatic transitions between lambda-doubling states in the capture of a diatomic molecule by an ion

The low-energy capture of a dipolar diatomic molecule in an adiabatically isolated electronic state with a good quantum number Omega (Hund's coupling case a) by an ion occurs adiabatically with respect to rotational transitions of the diatom. However, the capture dynamics may be nonadiabatic with respect to transitions between the pair of the Lambda-doubling states belonging to the same value of the intrinsic angular momentum j. In this work, nonadiabatic transition probabilities are calculated which define the Lambda-doubling j-specific capture rate coefficients. It is shown that the transition from linear to quadratic Stark effect in the ion-dipole interaction, which damps the T(-1/2) divergence of the capture rate coefficient calculated with vanishing Lambda-doubling splitting, occurs in the adiabatic regime with respect to transitions between Lambda-doubling adiabatic channel potentials. This allows one to suggest simple analytical expressions for the rate coefficients in the temperature range which covers the region between the sudden and the adiabatic limits with respect to the Lambda-doubling states.     

Substituent effects on the vibronic coupling for the phenoxyl/phenol self-exchange reaction

The impact of substituents on the vibronic coupling for the phenoxyl/phenol self-exchange reaction, which occurs by a proton-coupled electron transfer mechanism, is investigated. The vibronic couplings are calculated with a grid-based nonadiabatic method and a nuclear-electronic orbital nonorthogonal configuration interaction method. The quantitative agreement between these two methods for the unsubstituted phenoxyl/phenol system and the qualitative agreement in the predicted trends for the substituted phenoxyl/phenol systems provides a level of validation for both methods. Analysis of the results indicates that electron-donating groups enhance the vibronic coupling, while electron-withdrawing groups attenuate the vibronic coupling. Thus, if all other aspects of the reaction are the same, then electron-donating groups will increase the rate, while electron-withdrawing groups will decrease the rate. Correlations between the vibronic coupling and physical properties of the phenol are also analyzed. Negative Hammett constants correspond to higher vibronic couplings, while positive Hammett constants correspond to similar or slightly lower vibronic couplings relative to the unsubstituted phenoxyl/phenol system. In addition, lower bond dissociation enthalpies, ionization potentials, and redox potentials, as well as higher pKa values, tend to correspond to higher vibronic couplings relative to the unsubstituted phenoxyl/phenol system. The observed trends enable the prediction of the impact of general substituents on the vibronic coupling, and hence the rate, for the phenoxyl/phenol self-exchange reaction. The fundamental physical insights obtained from these studies are applicable to other proton-coupled electron transfer systems.     

Nuclear motion on the orbitally degenerate electronic ground state of fully deuterated triatomic hydrogen

Nuclear motion in the vicinity of conical intersections of the degenerate electronic ground state of fully deuterated triatomic hydrogen, D(3), is investigated with the aid of a time-dependent wavepacket approach in hyperspherical coordinates. Vibronic energy level spectra and the eigenfunctions are examined by including, for example, (1) geometric phase (GP) correction, (2) diagonal Born-Huang (BH) correction, and (3) both GP and BH corrections to the Born-Oppenheimer adiabatic Hamiltonian and finally by considering the nonadiabatic coupling between the two electronic surfaces explicitly. It emerges from this study that inclusion of both the GP and BH corrections is insufficient to explain the spectral features observed in the experiment. The latter are recovered by considering the complete two-states coupled Hamiltonian only. This study shows that both the GP and BH corrections constitute a minor part of the surface coupling effects, in particular, on the dynamics of the upper adiabatic sheet. Most importantly, we add that the experimental signature of the GP effect appears only in the observed shift of the eigenlevels of the electronic state when compared to those obtained from a completely Born-Oppenheimer Hamiltonian. The detail fine structure of the observed band of the electronic state is shaped by the off-diagonal derivative coupling elements of the nonadiabatic coupling operator.