Renner-Teller nonadiabatic coupling terms: An ab-initio study of the HNH molecule

In this Communication we present the first theoretical/numerical treatment of nonadiabatic coupling terms (NACT) that originate from the Renner-Teller (RT) model, namely, those that follow from the splitting of an electronic level of a linear molecule when it becomes bent. These two newly formed states are characterized by different symmetries and are designated as A and B. Our main findings: (1) The RT NACTs are quantized as long as they are calculated close enough to collinear configuration of the molecule (in this case HNH). Their value is tau = 1 (the Jahn-Teller values in similar situations, are tau = (1/2)). (2) Calculation of RT NACTs at bent configurations (i.e., at a distance from the linear axis) yield decreased values, sometimes by more than 50%. This last finding implies that in strongly bent configurations the two-state Hilbert subspace (formed by the above mentioned A and B states) is affected by upper states, most likely via Jahn-Teller conical intersections. (3) This study has also important practical implications. The fact that the RT NACTs decrease in (strongly) bent situations implies that analyzing spectroscopic data employing only the two Pi-states may not be sufficient in order to achieve the required accuracy.     

Curl equations as substratum for the derivation of the electronic nonadiabatic coupling terms

The idea to derive the nonadiabatic coupling terms by solving the Curl equations (Avery, J.; Baer, M.; Billing, G. D. Mol Phys 2002, 100, 1011) is extended to a three‐state system where the first and second states form one conical intersection, i.e., τ₁₂ and the second and the third states form another conical intersection, i.e., τ₂₃. Whereas the two‐state Curl equations form a set of linear differential equations, the extension to a three‐state system not only increases the number of equations but also leads to nonlinear terms. In the present study, we developed a perturbative scheme, which guarantees convergence if the overlap between the two interacting conical intersections is not too strong. Among other things, we also revealed that the nonadiabatic coupling term between the first and third states, i.e., τ₁₃ (such interactions do not originate from conical intersection) is formed due to the interaction between τ₁₂ and τ₂₃. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

Assigning signs to the electronic nonadiabatic coupling terms: the H2,O system as a case study

This paper is devoted to a specific difficulty related to the electronic nonadiabatic coupling terms (NACT), namely, how to determine correctly their signs. It is well known that correct NACTs, including their signs, are crucial for any numerical treatment of the nuclear Schrodinger equation [see, i.e., A. Kuppermaan and R. Abrol, Adv. Chem. Phys. 124, 283 (2003)]. In most cases the derivation of the correct sign of the nonadiabatic coupling matrix (NACM) is done employing various continuity procedures. However, there are cases where these procedures do not suffice and for these cases we suggest to apply an additional procedure based on a mathematical lemma which asserts that the exponentiated line integral which yields the D matrix is invariant with respect to the initial point of the integration [M. Baer, J. Phys. Chem. A 104, 3181 (2000)]. In the numerical study we apply this lemma to determine the signs of the 3x3 NACM elements for the three excited states of the {H(2),O} system (some of these NACTs are presented here for the first time). It turns out that the ab initio treatment yields results from which one can form eight different 3x3 NACMs. However the application of this lemma (which does not require any significant additional numerical effort) reduces this number to two. The final selection is done by an enhanced numerical study which requires more accurate calculations.     

Organic functional group analysis at the micromole level. I

Microchimica Acta - The N-oxides of pyridine and of its homologs and cyano, nitro and halogenated derivatives have been studied by differential pulse polarography. The sensitivity is approximately...     

Interactions between frenkel excitons. I. Formalism and preliminary results

For the tightest-bound Frenkel excitons a many-particle hamiltonian is adopted which includes terms representing interactions between the excitons. The parameters of the model are the excitation energy of an isolated localized excitation, the excitation-transfer matrix elements, and the interaction energies. Several special cases, involving particular relations between the latter two sets of parameters, are treated qualitatively. Equations are translated into the language of spin-$\text{1}\text{2}$ lattices, so that use may be made of the results of that theory. Favorable conditions for observation of polyexcitons and of phase transition to a liquid of excitons are discussed. A standard formalism for the determination of absorption and emission spectra and of their moments is adapted to the present problem. Possible generalizations of the model are briefly discussed.     

Ab initio electron propagator theory of molecular wires. I. Formalism

Ab initio electron propagator methodology may be applied to the calculation of electrical current through a molecular wire. A new theoretical approach is developed for the calculation of the retarded and advanced Green functions in terms of the electron propagator matrix for the bridge molecule. The calculation of the current requires integration in a complex half plane for a trace that involves terminal and Green's-function matrices. Because the Green's-function matrices have complex poles represented by matrices, a special scheme is developed to express these "matrix poles" in terms of ordinary poles. An expression for the current is derived for a terminal matrix of arbitrary rank. For a single terminal orbital, the analytical expression for the current is given in terms of pole strengths, poles, and terminal matrix elements of the electron propagator. It is shown that Dyson orbitals with high pole strengths and overlaps with terminal orbitals are most responsible for the conduction of electrical current.     

Reclassifying exciton-phonon coupling in molecular aggregates: evidence of strong nonadiabatic coupling in oligothiophene crystals

Exciton-phonon (EP) coupling in molecular aggregates is reexamined in cases where extended intermolecular interactions result in low-energy excitons with high effective masses. The analysis is based on a single intramolecular vibrational mode with frequency omega0 and Huang-Rhys factor lambda2. When the curvature Jc at the exciton band bottom is much smaller than the free-exciton Davydov splitting W, the strength of the EP coupling is determined by comparing the nuclear relaxation energy lambda2omega0 with the curvature. In this way, weak (lambda2omega0<4piJc), intermediate I (lambda2omega0 approximately 4piJc), and strong I (lambda2omega0>4piJc) coupling regimes are introduced. The conventional intermediate (lambda2omega0 approximately W) and strong (lambda2omega0>W) EP coupling regimes originally defined by Simpson and Peterson [J. Chem. Phys. 26, 588 (1957)] are based solely on the Davydov splitting and are referred to here as intermediate II and strong II regimes, respectively. Within the intermediate I and strong I regimes the near degeneracy of the low-energy excitons allows efficient nonadiabatic coupling, resulting in a spectral splitting between the b- and ac-polarized first replicas in the vibronic progression characterizing optical absorption. Such spectral signatures are clearly observed in OT4 thin films and crystals, where splittings for the lowest energy mode with omega0=161 cm(-1) are as large as 30 cm(-1) with a small variation due to sample disorder. Numerical calculations using a multiphonon BO basis set and a Hamiltonian including linear EP coupling yield excellent agreement with experiment.     

Spheroidal multipolar expansions for intermolecular energy in crystal theory. I. Formalism

In view of further applications to molecular crystal theory, spheroidal expansion of intermolecular potential involving. anisotropic molecules is investigated, including multipole moments up to octupole. The spheroidal formalism is worked out for the electrostatic potential outside a non-spherica charge distribution as well as for interactions between two anisotropic molecules in the non-overlap assumption. The case where either molecule is axia symmetric is presented as a reduced form of the general formulas displayed.     

Computation of expectation values from vibrational coupled-cluster at the two-mode coupling level

In this work we show how the vibrational coupled-cluster method at the two-mode coupling level can be used to calculate zero-point vibrational averages of properties. A technique is presented, where any expectation value can be calculated using a single set of Lagrangian multipliers computed solving iteratively a single linear set of equations. Sample calculations are presented which show that the resulting algorithm scales only with the third power of the number of modes, therefore making large systems accessible. Moreover, we present applications to water, pyrrole, and para-nitroaniline.     

Derivation of the electronic nonadiabatic coupling field in molecular systems: an algebraic-vectorial approach

In this Communication it is suggested that various elements of the nonadiabatic coupling matrix, tau(jk)(s) are created by the singular nonadiabatic coupling terms of the system. Moreover, given the spatial distribution of these coupling terms in the close vicinity of their singularity points yields, according to this approach, the integrated intensity of the field at every point in the region of interest. To support these statements we consider the conical intersections of the three lower states of the H+H(2) system: From an ab initio treatment we obtain the nonadiabatic coupling terms around each conical intersection separately (at its close vicinity) and having those, create the field at every desired point employing vector-algebra. This approach is also used to calculate the intensity of the Curl of those matrix elements that lack their own sources [tau(13)(s) in the present case]. The final results are compared with relevant ab initio calculations.     

Study of nonadiabatic effects in the Li-Li 2 + system. Location of nonadiabatic regions

Nonadiabatic coupling between the lowest two singlet potential energy surfaces of the Li-Li ₂ ⁺ system is calculated using the diatomics-in-molecules method. Location of nonadiabatic regions in the configuration space of Li-Li ₂ ⁺ and their analysis is used to estimate those inner and translational states of the reactants which can lead to nonadiabatic behavior.     

Density functional perturbational orbital theory of spin polarization in electronic systems. I. Formalism

A perturbational approach is presented for the general analysis of spin-polarization effect on electronic structures and energies within spin-density functional formalism. Explicit expressions for the changes in Kohn-Sham [Phys. Rev. 140, 1133 (1965)] orbital energies and coefficients as well as for the change in total electronic energy are derived upon using the local spin density and self-interaction-corrected exchange-correlation functionals. The application of the method for atoms provides analytical expressions for the exchange splitting energy and spin-polarization energy. The atomic exchange parameters are obtained from the expressions for the elements with Z=1-92 and they match well with Stoner exchange parameters for 3d metal elements.     

Direct perturbation theory in terms of energy derivatives: fourth-order relativistic corrections at the Hartree-Fock level

In this work, the quantum-chemical treatment of relativistic effects by means of direct perturbation theory is extended from its lowest order, DPT2, to the next higher order, DPT4. The required theory is given in terms of energy derivatives with the DPT4 energy correction defined as the corresponding second derivative with respect to the relativistic perturbation parameter λ(rel) = c(2) and c as the speed of light. To facilitate the implementation in standard quantum-chemical program packages, a general formulation of DPT starting from a nonrelativistic Lagrangian is developed, thereby expanding both wave function and operators in terms of λ(rel). The corresponding expressions, which incorporate in an additive manner scalar-relativistic and spin-orbit contributions, are given at the Hartree-Fock level and have been implemented in the CFOUR program package using the available analytic second-derivative techniques. The accuracy of the DPT4 corrections at the HF level is investigated by comparison with rigorous four-component calculations. Scalar-relativistic and spin-orbit contributions are analyzed individually and the importance of the various terms to those corrections is discussed. Furthermore, the basis-set dependence of the computed DPT4 corrections is investigated.     

Polybenzopinacols. I. Photolytic coupling of aromatic diketones

Low molecular weight polybenzopinacols were obtained by the photolytic coupling of m- and p-dibenzoylbenzene and 4,4′-dibenzoyldiphenyl ether in isopropanol, tetrahydrofuran–isopropanol, benzene–isopropanol, and benzene–ethanol solutions. The polypinacols were soluble in common organic solvents such as tetrahydrofuran, ether, and benzene. The inherent viscosities ranged from 0.06 to 0.14. Average molecular weight (M?_n) data indicated that the polymers were mostly dimers and trimers.     

Rovibrational spectra of ammonia. I. Unprecedented accuracy of a potential energy surface used with nonadiabatic corrections

In this work, we build upon our previous work on the theoretical spectroscopy of ammonia, NH(3). Compared to our 2008 study, we include more physics in our rovibrational calculations and more experimental data in the refinement procedure, and these enable us to produce a potential energy surface (PES) of unprecedented accuracy. We call this the HSL-2 PES. The additional physics we include is a second-order correction for the breakdown of the Born-Oppenheimer approximation, and we find it to be critical for improved results. By including experimental data for higher rotational levels in the refinement procedure, we were able to greatly reduce our systematic errors for the rotational dependence of our predictions. These additions together lead to a significantly improved total angular momentum (J) dependence in our computed rovibrational energies. The root-mean-square error between our predictions using the HSL-2 PES and the reliable energy levels from the HITRAN database for J = 0-6 and J = 7∕8 for (14)NH(3) is only 0.015 cm(-1) and 0.020∕0.023 cm(-1), respectively. The root-mean-square errors for the characteristic inversion splittings are approximately 1∕3 smaller than those for energy levels. The root-mean-square error for the 6002 J = 0-8 transition energies is 0.020 cm(-1). Overall, for J = 0-8, the spectroscopic data computed with HSL-2 is roughly an order of magnitude more accurate relative to our previous best ammonia PES (denoted HSL-1). These impressive numbers are eclipsed only by the root-mean-square error between our predictions for purely rotational transition energies of (15)NH(3) and the highly accurate Cologne database (CDMS): 0.00034 cm(-1) (10 MHz), in other words, 2 orders of magnitude smaller. In addition, we identify a deficiency in the (15)NH(3) energy levels determined from a model of the experimental data.     

Chemometric evaluation of surface water quality at regional level

Summary To reduce the volume of information needed for the evaluation of the quality of some surface waters designated for purification, data gathered in a previous study on three rivers (Simeto, Alcantara and Oreto) in Sicily have been subjected to statistical multivariate analysis. The data base comprised twentyfive variables from eight sampling points at monthly intervals for one year. Such a complex matrix cannot easily be interpreted and a statistical package have been used to put the data into a more useful form. The three adopted methods (KNN, LDA and SIMCA) gave partially different results as regarding confidence limits and basically convergent results as regarding the geographic classification of the surface waters examined. This analysis also shows that the measured variables cluster into different groups, according to their information content, and provide a reliable criterion for the choice of optimal parameters.     

A computationally efficient exact pseudopotential method. I. Analytic reformulation of the Phillips-Kleinman theory

Even with modern computers, it is still not possible to solve the Schrodinger equation exactly for systems with more than a handful of electrons. For many systems, the deeply bound core electrons serve merely as placeholders and only a few valence electrons participate in the chemical process of interest. Pseudopotential theory takes advantage of this fact to reduce the dimensionality of a multielectron chemical problem: the Schrodinger equation is solved only for the valence electrons, and the effects of the core electrons are included implicitly via an extra term in the Hamiltonian known as the pseudopotential. Phillips and Kleinman (PK) [Phys. Rev. 116, 287 (1959)]. demonstrated that it is possible to derive a pseudopotential that guarantees that the valence electron wave function is orthogonal to the (implicitly included) core electron wave functions. The PK theory, however, is expensive to implement since the pseudopotential is nonlocal and its computation involves iterative evaluation of the full Hamiltonian. In this paper, we present an analytically exact reformulation of the PK pseudopotential theory. Our reformulation has the advantage that it greatly simplifies the expressions that need to be evaluated during the iterative determination of the pseudopotential, greatly increasing the computational efficiency. We demonstrate our new formalism by calculating the pseudopotential for the 3s valence electron of the Na atom, and in the subsequent paper, we show that pseudopotentials for molecules as complex as tetrahydrofuran can be calculated with our formalism in only a few seconds. Our reformulation also provides a clear geometric interpretation of how the constraint equations in the PK theory, which are required to obtain a unique solution, are themselves sufficient to calculate the pseudopotential.