Electronic Transitions of Molecules: Vibrating Lewis Structures
Since the conception of the electron pair bond, Lewis structures have been used to illustrate the electronic structure of a molecule in its ground state. But, for excited states, most descriptions rely on the concept of molecular orbitals. In this work we demonstrate a simple and intuitive description of electronic resonances in terms of localized electron vibrations. By partitioning the 3N -dimensional space of a manyelectron wavefunction into hyper-regions related by permutation symmetry, chemical structures naturally result which correspond closely to Lewis structures, with identifiable single and double bonds, and lone pairs. Here we demonstrate how this picture of electronic structure develops upon the admixture of electronic wavefunctions, in the spirit of coherent electronic transitions. We show that π − π * transitions correspond to double-bonding electrons oscillating along the bond axis, and n − π * transitions reveal lone-pairs vibrating out of plane. In butadiene and hexatriene, the double-bond oscillations combine with in-and out-of-phase combinations, revealing the correspondence 1 between electronic transitions, molecular normal mode vibrations, and molecular plasmonics. This analysis allows electronic excitations to be described by building upon ground state electronic structures, without the need for molecular orbitals.
Introduction<!>Methods Formalism<!>Wavefunctions<!>Results and Discussion<!>α β<!>Conclusions
<p>The chemist's picture of electronic structure is highly adaptable. For depicting ground state structures, it is usual to draw one or more Lewis structures, where valence electrons are assigned to bonds, or lone pairs. 1 The theoretical foundation of this picture is accounted for by valence bond theory, which is therefore more intuitive than the molecular orbital theory approach. 2 However, the understanding of excited states of molecules is almost always within the framework of excitations of electrons within molecular orbitals (MOs). The MO approach is particularly prevalent in the description of conjugated systems, where the delocalized electronic orbitals are separated by their symmetry in π and σ representations. 3 For accurate calculations of excited states, for reasons of computational simplicity one is left with no choice but to take an approach based on molecular orbitals. However, as is well known, the calculated wavefunctions are invariant to the choice of orbitals, providing they span the same Hilbert space. Also, interpretations of electronic structure based on orbitals alone ignore the effect of (anti)symmetry: In many cases the S 1 and T 1 states are largely described by the same orbital occupancy, but have entirely different wavefunctions. Furthermore, interpreting the excited-state wavefunction can be problematic when the transition is described by a multitude of single and double (or higher) excitations from the reference wavefunction. 4 It is desirable to arrive at a simpler view of excited states based on accurately calculated and correctly antisymmetrised wavefunctions, ideally providing insight in an intuition-free manner. Recently, we reported a partitioning of the 3N -dimensional ground-state wavefunction into "tiles". 5 The equivalent tiles, which do not overlap, are themselves 3N -dimensional objects and are related by permutation of like-spin electrons. Each tile contains all the information of the entire wavefunction, since other tiles are merely repeats related by the interchange of electrons of the same spin: Tiles are positive or negative depending on the sign of the permutation. We and others demonstrated that if the boundaries between likesigned tiles are defined by a Voronoi diagram, with sites given by the centre-of-mass of the tile, the tiles reproduce Lewis structures of molecules and the double-quartet structures of Linnett. [5][6][7][8][9][10] Further, by following the tiles as a function of a reaction coordinate, one obtains electron movements akin to the "curly arrows" beloved of organic chemists. 11 Our procedure to obtain these structures is designated dynamic Voronoi Metropolis sampling (DVMS). Figure 1 shows calculated DVMS structures for some simple molecules: Only one spin is shown, since they are both equivalent in these closed-shell species. The DVMS structures exhibit core electrons, single-bonding electrons, double "banana"-bonding electrons and non-bonding lone-pairs, corresponding closely to the Lewis structures. 9,10 The ground-state electronic wavefunction only contains nodes engendered by the antisymmetry required to satisfy Fermi-Dirac statistics. The tiles themselves are nodeless.</p><p>Excited-state wavefunctions also obey Fermi-Dirac statistics, and thus will also exhibit tiling.</p><p>However, they will exhibit more nodes, and/or different nodes to the ground state. 12,13 The excited-state wavefunction can thus also be explored by the DVMS procedure, but with nodes within the tiles.</p><p>In this work, we present an approach to describing molecular electronic transitions inspired by the time-dependence of coherent electronic transitions. If the time-dependent Schrödinger equation is solved for a molecule subject to an oscillating electric field, the electronic wavefunction, originally in the ground state, will evolve such that after the perturbation it can be described by a sum over electronic states. If the angular frequency, ω, of the oscillating electric field is well matched to the energy difference between the ground and excited state (resonance condition, ω ≈ ∆ε), then the resulting wavefunction will be an admixture of the ground and the excited state. This mixed state is time-dependent, and exhibits "quantum beats". Here we show that, within the ground state tile, the average electron positions can be plotted as a function of phase, ωt, for a wavefunction mixed between the ground and excited state of interest. The resulting electron motions reveal the connection between classical vibrational motion and electronic spectroscopy. This picture also clearly connects molecular electronic transitions with the plasmon resonances exhibited by metal nanoparticles, providing a clear framework for molecular plasmonics.</p><!><p>Energy eigenstates are "stationary" in that the observables are time-independent. However, a state that is written as an admixture of one or more energy eigenstates evolves in time, exhibiting "quantum beats". In order to maximize the magnitude of the quantum beats, we calculated the trajectory of the electron density for a maximally mixed wavefunction,</p><p>where Ψ i (x) and Ψ f (x) are respectively the initial and final electronic wavefunctions. The energy of the initial wavefunction is set arbitrarily to zero. The final-state energy is ε = ω. In this work we are interested in analysing wavefunction changes for excitations from the ground state, so Ψ i and Ψ f are the wavefunctions for the ground and excited states, respectively. We use Ψ i to define the tiling of Ψ(x, t) and for convenience we here use Ψ i and Ψ f that are normalised over a single tile.</p><p>The centroid of the tile (a 3N vector x) evolves as</p><p>where R is the region defining a single tile, xi and xf are the position expectation values of Ψ i and Ψ f over the tile and</p><p>The centroid thus oscillates about the average density centroids of the two states. The amplitude of this oscillation, x f i , may be partitioned in an analogous way to the DVMS centroid positions, as</p><p>were the r j are the N individual electron displacements. Their sum is the electronic transition dipole moment, µ f i , divided by the elemental charge, e.</p><p>As such, the oscillation direction and amplitude is directly related to the polarization and strength of the electronic transition.</p><!><p>All electronic wavefunctions were calculated using the FIREFLY program, 14 which is largely based on the GAMESS package. 15 Unless where stated, the 6-31G(d) basis was used. Ground and excited state wavefunctions were calculated at the CAS-CI level using ground state RHF orbitals, except in the case of ethylene.</p><!><p>The first excited state of ethylene, C 2 H 4 , has been the subject of controversy. The so-called While we cannot say how the position of a particular electron evolves (they are all equivalent), DVMS allows us to inspect how the average position of each electron evolves in time from the perspective of a single wavefunction tile. In DVMS, the σ − π description of the double bond is replaced by a "banana bond" picture, 17 with one electron of each spin lying above and below the bond axis (Figure 1).  3). The wavefunction is mixed between the ground 'N ' state and excited 'V ' state (Equation 1). Only one spin-set is shown, since both spins move identically. Right: The electron positions for the mid-point of the vibration (gold spheres), x = (x i + xf )/2, and the displacement vector x f i (blue arrows).</p><p>Since the motion of the centroid is described by an oscillation (Equation 3), it may be illustrated by vectors indicating the motions of particular positions ({r j }, Equation 4), as is routine for vibrational normal modes of molecules. This depiction is shown on the right of Figure 2, succinctly representing the correlated motion of the banana-bonding electrons in the C=C bond in the N − V transition. With this connection made, it is interesting to explore how other vibrational modes might correspond to different excited states.</p><p>In the MO picture, the singlet and triplet states derived from the same π − π * excitation are crudely described by the same determinants, but with a change of sign:</p><p>This is the case for many, if not most, chromophores. So, where the excited singlet (Figure 2) has the α and β electrons vibrating in-phase, the (m S = 0) triplet state has the α and β electrons vibrating out-of-phase, thereby exhibiting no electronic transition dipole moment.</p><p>The DVMS structure and displacement vector for the forbidden T 1 ← S 0 transition are shown in Figure 3.</p><!><p>Figure 3: The electron positions for the mid-point of the vibration, x = (x i + xf )/2, and the displacement vector x f i for the wavefunction mixed between the S 0 and T 1 states of ethylene. The α and β electrons oscillate out of phase about the same midpoint, so there is no transition dipole moment.</p><p>The present treatment may be extended to other chromophores. Formaldehyde is isoelectronic with ethylene, but exhibits a richer electronic spectrum on account of the involvement of the non-bonding, "lone pair" electrons. There are three low-lying excited states: π − π * , and two n − π * states. 18 From the standpoint of the MO approach, the latter two can be considered excitations from symmetry-adapted combinations of lone pairs giving rise to a 1 → b 1 and b 2 → b 1 transitions to 1 B 1 and 1 A 2 excited states respectively. The latter is of considerable photophysical interest. 19 The DVMS picture is shown in Figure 4.</p><p>Rotation about the z-axis of formaldehyde transforms as the a 2 representation in the C 2v point group. Since the ground state is 1 A 1 , the transition moment to the 1 1 A 2 state also transforms as a 2 and the transition is electric-dipole forbidden, but magnetic-dipole allowed.</p><p>The rotation of charge is evident in Figure 4: The lone-pair electrons rotate out of plane.</p><p>Surprisingly, at this level of theory ((6,4) CAS-CI) there is also considerable out-of-plane motion of the C-H bonding electrons.</p><p>The 2 1 A 1 ← 1 1 A 1 transition is analogous to the π − π * transition of ethylene, where the double-bonding electrons oscillate along the molecular symmetry axis. Finally, the 1 1 B 1 ← trans-butadiene is normally drawn with two conjugated double bonds, and its DVMS structure bears this out. Like for ethylene, accurate calculations for excited states of transbutadiene require attention to correlation and basis set. [20][21][22][23] The present wavefunctions were calculated at the (4,4) CAS-CI/6-31G(d) level, and we do not expect accurate energies or state-ordering. However, for ethylene, the qualitative nature of the excitation was not found to be sensitive to the level of theory.</p><p>The two lowest excited states are 1 1 B u and 2 1 A g (C 2h point group). The 2 1 A g state of trans-butadiene is dark, being described by two triplet configurations, one on each double bond. [24][25][26] We do not investigate this state here, since it is not amenable to a treatment rooted in time-dependent quantum mechanics. The lowest bright excited state is 1 1 B u .</p><p>The DVMS picture of the electron vibration arising from the 1 1 B u wavefunction mixed with the 1 1 A g ground state is shown in Figure 5. Here, eight electrons, four from each double bond, oscillate along the C=C bonds, all in-phase. Since the individual C=C bond transition dipoles add, the 1 1 B u ← 1 1 A g transition is seen to be strongly allowed.</p><p>The question naturally beckons as to whether there is an excitation described by the outof-phase excitation of the two double bonds. In trans-butadiene, such a transition would be forbidden, since the transition dipoles of the two double bonds cancel.</p><p>At the present level of theory this transition was found to be the 3</p><p>The corresponding electron vibration is illustrated in Figure 5. It stands to reason that this</p><p>The electron positions for the mid-point of the vibration, x = (x i + xf )/2, and the displacement vector x f i for the wavefunction mixed between the 1 1 A 1 (ground), and 1 1 A 2 , 2 1 A 1 and 1 1 B 1 excited states of formaldehyde. The vector length is doubled for clarity.</p><p>The electron positions for the mid-point of the vibration, x = (x i + xf )/2, and the displacement vectors x f i , for the 1 1 B u ← 1 1 A g and 3 1 A g ← 1 1 A g transitions of transbutadiene. The vector length is doubled for clarity.</p><p>state would lie higher in energy than the bright state, since its opposing dipoles lie end-to-end, analogous to a molecular exciton (J-aggregate). 27 This is not the 2 1 A g ← 1 1 A g transition, as one might expect, because the 2 1 A g state is the dark triplet-pair state mentioned above.</p><p>The exact nature of the 2 1 A g state is a subject of ongoing investigation. 22,28 Like in trans-butadiene, the lowest excited state of trans,trans-hexatriene is described as a dark, doubly excited triplet pair. 29 The next lowest-energy states calculated are 1 1 B u , 3 1 A g and 2 1 B u . The electron vibrations corresponding to excitations to these states from mixing CAS-CI wavefunctions are illustrated in Figure 6. The lowest-energy bright state is 1 1 B u , and the electron vibration is the in-phase combination of double-bond oscillations.</p><p>The length of the displacement vector x f i is larger for the central double bond, compared to the outer bonds.</p><p>The excitation to the 3 1 A g state is dipole forbidden, and the displacement vector shows the outer double bonds vibrating out-of-phase, such that the local dipoles cancel. The excitation to the 2 1 B u state is symmetry allowed, but the displacement vectors partially This picture of electronic transitions in conjugated systems can be extended to cyclic systems. Cyclic polyenes require special attention owing to the emergence of resonance structures. Indeed, in a forthcoming publication we will demonstrate that the most appropriate description of benzene is a superposition of Kekulé structures where α and β electrons are correlated: They occupy alternate Kekulé structures. 30 In this picture, the in-phase vibration of the double bonds generates a cyclic motion of charge and therefore no net dipole.</p><p>This corresponds to the transition to the 1 B 1u (S 2 ) state. Transitions to the bright degenerate E 1u states (S 3 ) are described by vectors which have magnitudes given by cos θ and sin θ, where θ is the angle describing the position of the double bond on the ring. The S 1 ← S 0 transition is forbidden: The S 1 state of benzene is related to S 0 by a change in relative phase of the comprising (staggered-spin) Kekulé stuctures, and as such the tile centroids are unmoved.</p><p>Plasmon resonances in metal nanoparticles are described in a classical way, as the concerted oscillation of conducting electrons: There is no succinct molecular orbital description. However, the term Molecular Plasmonics has gained traction in recent years. [31][32][33] In molecules, the electronic transitions analogous to plasmon resonances would be characterized by concerted, in-phase oscillations of the double-bonding electrons. As shown here, the bright 1 B u states of polyenes exhibit this property: The double-bonds oscillate in-phase, giving rise to a large transition dipole. In this picture, as polyenes increase in length, the dipole-dipole coupling of each double-bond oscillation will lower the frequency of the resonance. The naturally occuring β-carotene has 11 conjugated double bonds, and a "plasmon resonance" at 2.6 eV. For the polyenes, the 1 B u plasmon band converges to an energy near 2 eV for long chains (n 20). 34</p><!><p>Using the dynamic Voronoi Metropolis sampling procedure (DVMS), we defined wavefunction tiles for a range of chromophores and showed that the centroids corresponded closely with the canonical Lewis structures. We investigated the centroid of the electron density in the tile for time-dependent wavefunctions mixed between the ground and various excited states and found that the centroid oscillates with an amplitude proportional to the transition dipole moment. For polyenes, where the DVMS describes the structure as isolated double "banana bonds", it was found that transitions to various excited states may be described as the in-and out-of-phase vibrations of the individual double bonds. The most intense transition is characterized by an in-phase vibration of the double bonds, which provides a clear link to the phenomenon of plasmon resonances in metal nanoparticles. We anticipate that this method will allow a variety of electronic transitions to be interpreted in terms of electrodynamically-induced vibrations of the ground state electronic structure. Because this viewpoint does not require the invocation of molecular orbitals, it allows the essence of electronic spectroscopy to be described outside this paradigm.</p>