Journal of Physics: Condensed Matter
TOPICAL REVIEW Related content
- Foundations of the AdS5 times S5 Continuous vibronic symmetries in Jahn–Teller superstring: I Gleb Arutyunov and Sergey Frolov models - Cooperative Jahn–Teller phase transition of icosahedral molecular units Seyed H Nasrollahi and Dimitri D To cite this article: Raphael F Ribeiro and Joel Yuen-Zhou 2018 J. Phys.: Condens. Matter 30 333001 Vvedensky
- The many-nucleon theory of nuclear collective structure and its macroscopic limits: an algebraic perspective D J Rowe, A E McCoy and M A Caprio View the article online for updates and enhancements.
This content was downloaded from IP address 169.228.107.189 on 30/01/2019 at 01:23 Journal of Physics: Condensed Matter
J. Phys.: Condens. Matter 30 (2018) 333001 (20pp) https://doi.org/10.1088/1361-648X/aac89e
Topical Review Continuous vibronic symmetries in Jahn–Teller models
Raphael F Ribeiro and Joel Yuen-Zhou
Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA 92093, United States of America
E-mail: [email protected]
Received 11 October 2017, revised 26 April 2018 Accepted for publication 30 May 2018 Published 24 July 2018
Abstract Explorations of the consequences of the Jahn–Teller (JT) effect remain active in solid- state and chemical physics. In this topical review we revisit the class of JT models which exhibit continuous vibronic symmetries. A treatment of these systems is given in terms of their algebraic properties. In particular, the compact symmetric spaces corresponding to JT models carrying a vibronic Lie group action are identified, and their invariants used to reduce their adiabatic potential energy surfaces into orbit spaces of the corresponding Lie groups. Additionally, a general decomposition of the molecular motion into pseudorotational and radial components is given based on the behavior of the electronic adiabatic states under the corresponding motions. We also provide a simple proof that the electronic spectrum for the space of JT minimum-energy structures (trough) displays a universality predicted by the epikernel principle. This result is in turn used to prove the topological equivalence between bosonic (fermionic) JT troughs and real (quaternionic) projective spaces. The relevance of the class of systems studied here for the more common case of JT systems with only discrete point group symmetry, and for generic asymmetric molecular systems with conical intersections involving more than two states is likewise explored. Finally, we show that JT models with continuous symmetries present the simplest models of conical intersections among an arbitrary number of electronic state crossings, and outline how this information may be utilized to obtain additional insight into generic dynamics near conical intersections. Keywords: Jahn–Teller, vibronic, symmetry (Some figures may appear in colour only in the online journal)
1. Introduction vibrational mode [10], the so-called E ⊗ e system [11] (we employ the standard convention that the irreducible repre- Jahn–Teller (JT) models [1, 2] explain a rich variety of phe- sentation (irrep) corresponding to the electronic multiplet is nomena in condensed matter and chemical physics [3–6]. labeled by an upper case letter, while that of the vibrations is Modern studies have explored the role of JT distortions in given by a lower case). This model has been employed, for e.g. possible mechanisms for unconventional superconductiv- example, to describe the distortion of CuII and MnIII in an ity [6], colossal magnetoresistance [7], multiferroics [8], and octahedral environment [12] (figure1 ), and properties of tria- single-molecule transport [9]. The unifying feature of these tomic systems near the equilateral configuration [13]. A well- phenomena is that they involve significant coupling of orbital known characteristic of this system is that when only linear and vibrational degrees of freedom. vibronic couplings are included, it displays circular symmetry One of the simplest JT models is that consisting of a degen- [11] (figure 2). The reason is the linear E ⊗ e Hamiltonian is erate pair of electronic states coupled to a similarly degenerate invariant under simultaneous rotations of the electronic states
1361-648X/18/333001+20$33.00 1 © 2018 IOP Publishing Ltd Printed in the UK J. Phys.: Condens. Matter () 333001 Topical Review
Figure 1. Eg vibrational modes and electronic orbitals for a system with Oh symmetry. and vibrational coordinates (see equation (7)). The space of rotations of the plane [SO(2)] is topologically equivalent to the one-dimensional circle S1. Therefore, the symmetry group of linear E ⊗ e is continuous (as opposed to the discrete point groups). This has extreme implications, e.g. there exists a con- Figure 2. Ground and excited-state branches of the E ⊗ e APES. tinuous family of minima (trough) in the ground-state adiaba- tic potential energy surface (APES) (figure 2), the vibronic of a Berry phase [19–21]. On the other hand, if the warp- ground-state is doubly degenerate, and the vibronic (pseudo) ing generates new electronic-state intersections (esis, which angular momentum is quantized in odd half-integral units, include the ubiquitous conical intersections (cis) [22–24]) on thus indicating a vibronic motion with spinorial character the electronic ground-state APES in regions accessible at low- [3, 11]. All of these are surprising, as none are generic proper- energies, then the warped JT system becomes topologically ties of finite molecular systems irrespective of the existence of different [19, 25]. In other words, the symmetry and degen- point group symmetry. They are also intrinsically quant um- eracy of the vibronic ground-state changes. While it is well- mechanical (even though a semiclassical treatment leads to an known that warping necessarily generates additional esis on effective classical Hamiltonian which describes the properties the JT APES [19, 26, 27], these are brought from the vibra- of the system when → 0 [14]). tional configuration space infinity (where they coalesce when Continuous symmetries in JT systems have sometimes nonlinear vibronic couplings vanish) (see figure 3). Thus, if been described as accidental or emergent, since the molecular the warping is weak enough, and the new esis are as a result Hamiltonian is only constrained to be a molecular point (dou- far away, then at low-energies the molecular system will be ble) group scalar. Deeper mathematical analysis has revealed confined to a region that excludes esis other than that which the continuous invariance properties of JT models follows defines the JT model; in this case, homotopic invariants of the from the interplay between the representations of molecular continuously-symmetric electronic ground-state APES will be point (double) groups adopted by the electronic and nuclear preserved [19, 20]. degrees of freedom [15–17]. In particular, Pooler laid out the A less known example of the robustness of the properties of JT representation theory underlying linear JT problems with con- models with continuous symmetries was given by Markiewicz tinuous symmetries [15, 16]. Roughly speaking, for the lat- [28]. To understand, recall that continuous-symmetry breaking ter to occur, the JT distortions need to be isotropic, and the in the dynamical JT problem without quadratic or higher-order space of electronic Hamiltonians must be equivalent to the vibronic couplings will happen under one of the following vibrational configuration space as irreducible representations three conditions: (a) different vibrational frequencies for the (irreps) of a Lie group G. In this way, the molecular Hilbert JT active modes irreps, but equal JT stabilization energies, space H≡Hvib ⊗Hel carries a G-action. (b) equal vibrational frequencies, but different JT stabiliza- Due to the constraints on the fundamental parameters of tion energies, and (c) different vibrational frequencies and JT continuously-symmetric JT models, these have been only stabilization energies. In [28], the linear E ⊗ (b1 ⊕ b2) model rarely employed to extract quantitative information about (which can be understood as arising from symmetry break- physical systems (though some examples are given below). ing of the vibrational modes of E ⊗ e) was employed to study Quadratic or non-degenerate linear vibronic couplings are case (a). This has a continuous set of ground-state minima, known to break continuous symmetries [3, 18]. For instance, but pseudoangular momentum is not conserved (the b1 and in the presence of quadratic couplings, the continuous set of b2 vibrations have different frequencies). However, the aver- E ⊗ e minima imposed by SO(2) symmetry becomes a discrete age wavepacket pseudoangular momentum was numerically set separated by saddle points [3, 18, 19]. The APES associ- verified to be quantized for different values of vibrational fre- ated to this model is said to be warped (figure3 ). However, if quency anisotropy. While the dynamical continuous symmetry the distorted JT surfaces can be obtained by a gap-preserving was broken, it left clear signatures. continuous deformation (homotopy) of the electronic ground- There have been also some experimental studies of state APES, then basic features of the continuously-symmetric JT centers which benefited from an analysis based on an JT system will remain relevant (figure 3) [19]. These include ideal model admitting continuous symmetries. For exam- the vibronic ground-state degeneracy and symmetry, as well ple, O’Brien [29] employed the SO(3)-invariant version of + as any other non-trivial effects originating from the existence T ⊗ (e ⊕ t2) to investigate the spectra of F centers (crystal
2 J. Phys.: Condens. Matter () 333001 Topical Review
Figure 3. Left: weakly warped APES for E ⊗ e; right: strongly-warped APES for the same model. vacancies occupied by holes) in CaO, and obtained quantita- 2 tively accurate absorption bandshapes. The same model also H(Q, P)=Hvib(Q, P)+HJT(Q), Q ∈ R , (1) provided accurate fits to the absorption spectrum of F-centers where Q =(Q1, Q2), and P =(P1, P2) are the displacements (electron-occupied vacancies) in CsF. Another example is from the JT center (the molecular shape with Q = 0, hosting ’ ’ ’ P2 → P2 Pooler s and O Brien s [30] study of the 1/2 3/2 line- the electronically degenerate multiplet) and their canonically shapes of Tl atoms in halides. Tl is a heavy-atom and its conjugate momenta, respectively (figure 1). The vibrational energy levels are strongly influenced by spin–orbit coupling. contribution is given by the 2D isotropic harmonic oscillator The JT problem for this model is Γ8 ⊗ (e ⊕ t2), where Γ8 (also Hamiltonian, denoted by G3/2) refers to the quartet double group representa- 2 2 P1 P2 1 2 2 2 vib(Q, P)= + + ω ( + ), (2) tion of cubic systems [31]. With the assumption of degener- H 2 2 2 Q1 Q2 ate vibronic couplings and harmonic frequencies, qualitative agreement with exper imental data was also obtained. where ω is the vibrational frequency for the e modes. The JT In summary, while strong constraints need to be satisfied Hamiltonian HJT(Q) contains the interaction between nuclear ⊗ for a JT system to feature continuous symmetries, i.e. to carry and electronic degrees of freedom. Given that E e is a TRI a Lie group action, their study adds valuable insight to the spinless model, the time-reversal symmetry operator T satis- 2 generic problem of cis, as they are minimal models featuring fies T = 1, and we can take the electronic Hilbert space Hel universal properties of systems that contains them. Moreover, to be a real vector space with basis functions invariant under ⊗ the high symmetry of such systems imply that they form T. Thus, the vibronic E e Hamiltonian can be written as: convenient starting points for the investigation of the com- −i σ2 φ i σ2 φ HJT(Q)=FQe 2 σ3e 2 , (3) plex interplay between orbital, structural, and charge order in extended systems with JT active sites. where F ∈ R −{0} is the reduced vibronic coupling constant The connection between JT models and Lie groups has 2 2 + (from the Wigner–Eckart theorem), Q = Q1 + Q2 ∈ R , been explored before [15, 16, 18, 32–38]. In this work we tanφ = / , φ ∈ 1 σ , = 1, 2, 3, review and present some novel perspectives on this relation- Q2 Q1 S , and j j are the Pauli matri- H ship based on the theory of symmetric spaces [39–41]. The ces acting on el. A change of electronic frame (basis) that Q general relevance of symmetric spaces to quantum mechanics preserves the reality of the electronic eigenvectors for all can be illustrated with the Threefold way classification of non- can be parametrized by −i σ2 θ 1 relativistic Hamiltonian ensembles [42]. Dyson proved that, M(θ)=e 2 ∈ SO(2), θ ∈ S . (4) in the absence of anticommuting symmetries, any quantum- mechanical system belongs to either one of three symmetry Let R(θ) be the rotation of vibrational configuration space classes: real orthogonal, real symplectic and complex unitary. defined by The real orthogonal class contains all time-reversal invariant R(θ) : (Q, φ) → (Q, φ + θ). (5) (TRI) bosonic systems, while the TRI fermionic Hamiltonians φ → φ + θ correspond to the real symplectic, and the complex unitary It is equivalent to mapping in equation (3). Hence, ensemble contains models with broken time-reversal symme- a change of electronic frame is equivalent to a rotation of the try. Dyson’s classification (and its generalization 43[ , 44]) is vibrational configuration space, the basis for the application of random matrix models in con- −1 −1 M(θ)HJT (Q) M (θ)=HJT R (θ)Q . (6) densed matter, nuclear [45], and chemical physics [46, 47]. In fact, as we will see later, there exist several connections As we now see, the JT Hamiltonian (and therefore, also the between topological phenomena arising in JT models and molecular) is invariant under the simultaneous action of SO(2) condensed matter physics (section 3.3). on the electronic and vibrational Hilbert spaces defined by Let us illustrate some of the previous remarks with the −1 M(θ)HJT [R(θ)Q] M (θ)=HJT(Q). (7) simplest JT model with continuous symmetry, the linear E ⊗ e system. The molecular Hamiltonian can be written as The equation above explains the circular symmetry of a sum of a purely vibrational part and a vibronic component, figure 2: the SO(2) action on the space of molecular (nuclear)
3 J. Phys.: Condens. Matter () 333001 Topical Review
Figure 4. Lhs: E ⊗ e contours with equal adiabatic electronic spectrum; Rhs: reduced orbit space.
Figure 5. Pseudorotational motion in vibrational and electronic configuration space (the orbital at the center is the adiabatic electronic ground-state). geometries maps a nuclear configuration into another with the We also find that every molecular geometry of E ⊗ e lives same Q, but that is equivalent to a change of (electronic) basis in an SO(2) orbit labeled by the radial coordinate (figure 4). of Hel. Therefore, all invariant properties of the electronic Formally, to each Q = 0 we assign the space of molecular Hamiltonian are preserved (e.g. its eigenvalues) under a vibra- structures tional space rotation (pseudorotation). This explains the choice T T O(Q)={φ ∈ [0, 2π)|R(φ)(Q,0) =(Q1, Q2) }, of displacement coordinates for JT distortions (Q = 0): φ is (9) changed under SO(2) action on the vibrational configuration containing all Q = Q(cos φ, sin φ) with equal electronic space (pseudorotation), while variation of Q corresponds to energy spectrum. Therefore, the APES can be decomposed (radial) motion that modifies the adiabatic electronic spectrum. into a one-dimensional space of SO(2) orbits, i.e. Equation (7) demonstrates the invariance of the linear E ⊗ e Q ∈ O(Q ). Hamiltonian with respect to the continuous SO(2) action on (10) 0 electronic and nuclear degrees of freedom. It also implies a Q conserved quantity by Noether’s theorem [48], the vibronic Each point O(Q) of the orbit space, with the exception of pseudoangular momentum J (setting = 1 hereafter), O(0), has an internal space corresponding to a circle S1 of molecular geometries which can be mapped into each other ∂ σ3 = −i + , by free pseudorotational motion (figure5 ). J ∂φ 2 (8) These simple considerations on the SO(2) invariance of where the first term is a linear operator on the vibrational E ⊗ e have allowed us to quickly understand a variety of its Hilbert space Hvib, while the second acts on the electronic Hel non-trivial features. In particular, (i) the APES was simpli- [11, 12] (see equation (3)). Equation (8) is a simple, but exotic fied by decomposing it into a space of orbits of the 2D proper result: it implies fractionalization of the vibronic pseudoangu- rotation group SO(2), (ii) the existence and meaning of radial lar momentum. and pseudoangular coordinates were explained, and (iii) the
4 J. Phys.: Condens. Matter () 333001 Topical Review
anomalous vibronic angular momentum conservation law of QΛ =( Λ λ , Λ λ ..., Λ λ ) where i Q i i1 Q i i2 Q i i|Λ | (with similar nota- E ⊗ e was quickly obtained. In more complex JT systems P ω i tion for Λi), mΛi and Λi are the mass and harmonic frequency carrying a Lie group action, the above points are going to be Λ of the nuclear displacements i, FΛi is the reduced vibronic generalized. For instance, there will be more than one angu- coupling constant associated with the same modes (from the lar and one radial coordinate, so that orbits will require more – V – Wigner Eckart theorem), and ΛiΓ is a vector of Clebsch than a single number to be specified uniquely. There will also VΛ Γ =( Λ λ Γ, Λ λ Γ..., Λ λ Γ)T Gordan matrices, i.e. i V i i1 V i i2 V i i|Λ | , be orbits of different types with qualitatively distinct internal i with e.g. (pseudorotational) motion. In the remainder of this topical review we will employ the VΛ λ Γ = Γγ |Λ λ Γγ |Γγ Γγ |. i i1 k i i1 l k l (13) invariance properties of symmetric spaces to investigate sev- γk∈Γ γl∈Γ eral aspects of the higher-dimensional generalizations of the Pooler has given the formal theory underlying the SO(2)-invariant E ⊗ e model discussed above. In the process construction of molecular JT Hamiltonians with continuous we will present simple proofs of established facts about JT symmetries, i.e. invariant under the action of a Lie group systems, as well as obtain some new fundamental results on [15, 16]. First, it was shown that the space of electronic tensor the topological properties of JT APESs. The content is organ- operators living in Γ ⊗ Γ can be decomposed into even and ized as follows: in section 2 we discuss the general relation- odd subspaces, which we will denote by M and P, respec- ship between linear JT problems with continuous symmetries tively. The operators in M belong to the symmetric part of and symmetric spaces. Section 3 applies the framework pre- Γ ⊗ Γ − A1 (where A is the totally symmetric irrep) if the sented in the previous to characterize local and global proper- 1 system is spinless (we discuss the fermionic case later). The ties of the JT APES. In section 4 we discuss the connection odd live in the antisymmetric. Therefore, P is a Lie subal- between the investigated models and molecular systems with gebra of Γ ⊗ Γ − A1 (since the commutator of antisymmetric accidental cis involving more than two APESs. A summary of matrices is antisymmetric). Conversely, the commutator of the presented material and a discussion of future directions is even operators is an odd element, so the latter do not span a given in section 5. Lie algebra. For example, in E ⊗ eM is spanned by the Pauli matrices σ1 and σ3, while σ2 is the generator of P. 2. Symmetric spaces underlying JT models with P generates the |Γ|-dimensional special orthogonal group, continuous symmetries SO(|Γ|). If M is isomorphic to an irrep Λ of SO(|Γ|), then H(Q, P) is invariant under the action of the same group on We start this section by revisiting the conditions satisfied by both electronic and nuclear degrees of freedom, if and only if JT systems admitting continuous symmetries. These are used the vibrational configuration space can also be embedded in to determine the symmetric spaces corresponding to each ele- (i.e. transforms like) the irrep Λ of SO(|Γ|) [16]. In E ⊗ e, σ2 ment of the set of JT models studied in this article. We also generates SO(2), as the elements of the latter can be written as −iφσ /2 discuss some general properties of symmetric spaces which U(φ)=e 2 . Conjugation of σ1 and σ3 with U(φ) shows will be useful for establishing several universal properties of that M transforms in the vector irrep of SO(2), and so do JT systems carrying a Lie group action. the JT coordinates Q1 and Q2 (see equations (2)–(6)). Hence, Consider a JT problem defined by the coupling between linear E ⊗ e shows vibronic SO(2) invariance. Γ electronic states belonging to the irrep of a point group S In other words, sufficient conditions for invariance under (S is the symmetry group of the nuclear geometry at the JT continuous transformations of electronic and vibrational Q = 0 {γ } center ) spanned by the basis vectors i i=1,2,...,|Γ|, degrees of freedom in JT models are that: (a) only linear and JT active vibrations transforming like the vectors of the vibronic couplings are non-vanishing, and (b) the space of Λ=⊕ Λ Λ (generically) reducible representation i i, where i is traceless electronic Hamiltonians and the vibrational configu- {λ , λ ...λ } a non-totally symmetric irrep of S with basis 1 2 |Λ | . i i i i ration space are isomorphic as irreps of a Lie group G. The electronic multiplet is only defined to belong to a single The above conditions require that the electronic multiplet irrep, for none of the studied models contain reducible elec- Γ = is equally coupled to every JT active mode, while mΛi m tronic representations (Wiseman showed how to generalize ω = ω Λ ⊂ Λ and Λi Λ for all i . Physically, they imply equal the theory of [16] to the case of direct sum electronic Hilbert JT stabilization energies for molecular distortions along each – spaces [37]). Generalizing equations (1) (3), the molecular Λi, and (pseudo)rotational symmetry in the purely vibrational Hamiltonian of the system including only linear vibronic cou- part of H(Q, P), respectively. plings is given by [3] The embeddings of the vibrational and electronic point Λ=⊕λ Γ P2 1 group irreps i and into Lie group irreps will not be Λi 2 2 H(Q, P)= + mΛ ωΛ QΛ + HJT(Q), subject of future discussion. The interested reader may consult 2 2 i i i (11) mΛi Λi Λi [15, 16, 49] for details. Henceforth, we take the conditions for continuous invariance of JT Hamiltonians to be fulfilled and HJT(Q)= FΛ ΓQΛ · VΛ Γ, study its consequences. For the sake of simplicity, we will not i i i (12) Λi make a distinction between the continuous and point group
5 J. Phys.: Condens. Matter () 333001 Topical Review irreps anymore, so that from now on Λ and Γ correspond to Table 1. JT models with continuous symmetries investigated in the appropriate continuous group irreps of SO(N) or USp(2N) this paper. The last column gives examples of discrete symmetry groups from which each model may be obtained by imposing the (in the fermionic case discussed below) depending on the total appropriate constraints [55]. S∗ refers to the double group obtained spin of the considered electronic JT multiplet. from a point group S. Hel P M The relations satisfied by the operators in and 2 can be summarized by JT model T Symmetric Space Examples E ⊗ e [11] 1 SU(2)/SO(2) C3 , O [M, M]=P, [P, M]=M, [P, P]=P. (14) v h T ⊗ (e ⊕ t2) [50] 1 SU(3)/SO(3) Td, O , I, I L = M⊕P h h They imply that is a Lie algebra [41]. In other G ⊗ (g ⊕ h) [35] 1 SU(4)/SO(4) I,Ih L words, corresponds to the tangent space of a Lie group L at H ⊗ (g ⊕ 2h) [36] 1 SU(5)/SO(5) I,Ih the (group) identity. The odd operators define a subalgebra P. Γ ( ) ⊗ ( ⊕ ) − ∗, ∗, ∗ 8 G3/2 e t2 [30] 1 SU(4)/USp(4) Td Oh Ih They generate the Lie subgroup P ⊂ L. Importantly, the elec- Γ ( ) ⊗ ( ⊕ 2 ) − ∗, ∗ 9 I5/2 g h [38] 1 SU(6)/USp(6) I Ih tronic matrices in M define tangent vectors for the symmetric space L/P [39–41]. In fact, the algebraic definition of symmetric 2 where kΛ = mω . Note the matrices Vα ⊂M and the spaces is encapsulated by the relations shown in equation (14). Λ so(|Γ|) generators are both real and traceless, so their direct Some of their relevant properties will be discussed below. sum is the algebra sl(|Γ|, R), i.e. L = sl(|Γ|, R). As a con- Time-reversal symmetry implies the existence of a basis sequence, the underlying symmetric spaces for the spinless where the electronic Hamiltonian matrix elements are real. class of JT Hamiltonians with continuous symmetries is Thus, only the electronic tensor operators of L⊂sl(|Γ|, R) SL(|Γ|, R)/SO(|Γ|) (or SU(|Γ|)/SO(|Γ|) if we take the skew- (the space of traceless real matrices) which are in the symmet- hermitian Hel time-evolution generators iVα as the generators ric part of Γ ⊗ Γ − A1 lead to allowed quantum Hamiltonians. of the symmetric space). Hence, HJT(Q) belongs to M = L/P. Moreover, the A fundamental property of symmetric spaces is their rank antisymmetric electronic tensor operators generate the special [39, 41]: the number of anisotropic spatial directions on these orthogonal group SO(|Γ|). Its generators define the Lie alge- manifolds. Physically, the rank of a given JT model is equal bra so(|Γ|): the space of infinitesimal rotations of Hel. Thus, to the number of linearly independent JT active nuclear dist- so(|Γ|) acts on the electronic Hilbert space operator HJT(Q) ortions unrelated by SO(|Γ|) (or USp(|Γ|) in the fermionic by infinitesimal rotations. As we show below, while the case) transformations on electronic or vibrational degrees defined so(N) action occurs on Hel, it is equivalent to infini- of freedom. In other words, motion along anisotropic direc- tesimal rotation of the internal displacements Q → Q + δQ in tions induces changes in the eigenvalues of the electronic JT models with continuous symmetries. Hamiltonian. For instance, in E ⊗ e only geometries with the In the spinless case, M⊕P spans either (a) sl(|Γ|, R) same radial coordinate Q are related by nuclear pseudorotation (see below for the T2 = −1 case), or (b) one of its proper (figure 4). Therefore, in the space of time-evolution operators subsets [15]. Examples of class (a) are the Lie group invari- of E ⊗ e, there exists a single direction which is anisotropic. ant for mulations of E ⊗ e [11], T ⊗ (e ⊕ t2) [50], T ⊗ h [51], This conforms with the fact that the rank of SU(2)/SO(2) is G ⊗ (g ⊕ h) [35], H ⊗ (g ⊕ 2h) [36, 49, 52]. These include JT one [39, 41]. active dist ortions representing all independent ways to split the Another view on the rank of a symmetric space is that it electronic degeneracy at Q = 0. In these models, the continu- specifies the dimensionality of its maximal abelian (commu- ous symmetry is maximal, since the symmetry of the molecular tative) subalgebra (also commonly denoted as Cartan subal- Hamiltonian cannot be increased without changing the number gebra) C [39]. By a theorem of symmetric spaces theory, any of dimensions of Hel. Conversely, class (b) contains H ⊗ h m ∈L/P is related to an element m of the Cartan subalgebra [36, 51], H ⊗ (g ⊕ h) [36], etc [15]. Their vibronic coupling C C by conjugation with an element of the Lie group P (with matrices only span a proper subset of sl(N, R)/so(N, R), and tangent space P at the identity of P) [39], the molecular Hamiltonian symmetry can be increased by −1 including the remaining symmetry-allowed independent JT m = gmCg , g ∈ P, mC ∈C. (16) active couplings. Thus, the latter models do not include one The rank of the symmetric space SL(|Γ|, R)/SO(|Γ|) (or or more of the possible JT active distortions. In this paper we SU(|Γ|)/SO(|Γ|)), is equal to |Γ|−1 [41]. focus on the JT models with maximal continuous symmetries We can always take the Cartan subalgebra of (see table 1), for they share a variety of deeply related proper- sl(|Γ|, R|)/so(|Γ|) to be given by |Γ|−1 linearly independ- ties. Thus, from now on, anytime we mention JT models car- ent traceless diagonal matrices of sl(|Γ|, R) [41]. In this case, rying a Lie group action, it should be understood that we are equation (16) is the simple statement that any real traceless referring to models in class (a), unless otherwise noted. symmetric matrix m can be diagonalized by an orthogonal The electronic Hamiltonian can only be constructed with transformation. We note that conjugation preserves com- the even tensor operators Vα ∈M. Therefore, the constraints mutation relations, so the choice of Cartan subalgebra is not given above imply the most general JT Hamiltonian for a unique. Thus, a set of diagonal matrices may be the simplest, Γ ⊗ Λ system with a continuous symmetry can be written as: but it is not the only. P2 1 (Q, P)= Λ + Q2 + , It follows that the electronic part of the molecular H 2 2kΛ Λ FΛ QαVα (15) Hamiltonian (equation (15)) can be rewritten as m α∈M
6 J. Phys.: Condens. Matter () 333001 Topical Review HJT(Q)=FΛ QαVα 3. Local and global structure of JT orbits α∈M We now take advantage of the framework introduced in the = (Q) −1(Q), previous section to quickly obtain insight into the APESs of U FΛ QαC VαC U JT models carrying a Lie group action. As we show in sec- αC∈MC tion 4, the results which we present here also have some U(Q) ∈ SO(|Γ|), (17) significance for the treatment of dynamics in the neighbor- where MC ⊂M is a choice of Cartan subalgebra for hood of electronic degeneracies of generic molecular systems M≡sl(|Γ|, R)/so(|Γ|). By the invariance of H(Q, P) under beyond JT models. Therefore, they display a degree of uni- simultaneous SO(|Γ|) action on electronic and vibrational versality which is perhaps unanticipated in view of the non- degrees of freedom, there exists an embedding of SO(|Γ|) generic constraints satisfied by the continuous symmetries of in SO(|Λ|) [since proper rotations (orthogonal transforma- the investigated models. tions with determinant equal to one) of the vibrational con- figuration space are described by the latter group]. Hence, 3.1. Symmetry-adapted coordinates and orbit spaces there is a continuous injective map, R :SO(|Γ|) → SO(|Λ|) satisfying Consider first the problem of choosing coordinates for vibra- −1 −1 tional motion, which are adapted to the invariance proper- U(Q)Qα Vα U (Q)= [R (β)Q ]αVα, C C C (18) ties of the underlying symmetric space. These have various α ∈M α C C advantages compared to cartesian coordinate systems, since where β denotes an SO(|Λ|) point, R(β) specifies the they reduce the complexity of molecular potential energy sur- −1 SO(|Λ|) rotation satisfying R (β)QC = Q, and QC is the faces (e.g. figure 4). For instance, in E ⊗ e [10, 11] and T ⊗ h nuclear displacement with nonvanishing components only [56, 57], the identification of radial and angular coordinates along the directions which are dual to the Cartan subalgebra introduces significant simplification to treatments of the
matrices VαC. Equation (18) is the mathematical represen- dynamical JT problem. They are also useful in describing the tation of the statement that in JT models with continuous warped ground-state APES [18]. symmetries, a rotation of the electronic frame is equivalent As we showed in section 1, coordinates adapted to the to a rotation of the space of JT distorted structures. It is the circular symmetry of the SO(2) E ⊗ e model can be trivially generalization of equation (6) for N-dimensional electronic obtained. However, this is not the case for higher-dimensional Hilbert spaces. models with SO(N), N > 2 or USp(2N) symmetries. The In cases where spin–orbit coupling is strong and time-rever- situation is ameliorated if only the space of minima of the JT sal symmetry is implemented by T satisfying T2 = −1, the JT APES is of interest, for the ground-state APES is homemor- active modes live in the antisymmetric part of the Γ ⊗ Γ − A1 phic with the real (if T2 = 1 [18]) or quaternionic (when space, where Γ is now a spinorial double group irrep [2, 53]. T2 = −1) projective space (see sections 3.2.2 and 3.2.3). The The electronic Hilbert space for a fixed geometry can thus be real projective space (RPN ) can be obtained from the sphere given the structure of a quaternionic vector space [54]. The SN by identifying antipodal points of the latter [58]. Therefore, unitary transformations which preserve this structure form hyperspherical coordinates are a natural choice for the study ∼ the unitary symplectic group U(|Γ|/2, H) = USp(|Γ|) [41]. of vibrational motion on the spinless ground-state trough In this case, the molecular Hamiltonian is invariant under the [35, 36] of the continuously-symmetric JT models. While simultaneous action of the symplectic group on the electronic profitable in discerning properties of the space of minima and nuclear degrees of freedom whenever the symmetric part of the electronic ground-state APES, this approach gives no of Γ ⊗ Γ generates the usp(|Γ|) algebra, the traceless antisym- insight into the motion which is normal to the extremal sub- ∗ ∗ ∼ metric part lives in su (|Γ|)/usp(|Γ|) (su (|Γ|) = sl(|Γ|/2, H) space of the APES, nor does it explain the spectral flow of the is the space of traceless |Γ|×|Γ| complex matrices invariant Born–Oppenheimer JT Hamiltonian in non-stationary regions under conjugation by an antiunitary operator T satisfying of the molecular vibrational configuration space. Thus, hyper- T2 = −1, or simply the space of traceless |Γ|/2 ×|Γ|/2 matri- spherical coordinates do not take full advantage of the invari- ces with quaternionic entries [41]), and the vibrational degrees ance properties of the systems studied here to maximally of freedom transform in the same irrep of USp(|Γ|) as the simplify the description of the molecular APES. traceless antisymmetric electronic tensor operators [15, 16]. We will employ the isomorphism between traceless elec- Therefore, in cases with strong spin–orbit coupling, the under- tronic Hamiltonian operators and vibrational displacements lying symmetric spaces for JT models carrying a Lie group as irreps of a Lie group G, to choose vibrational coordinates action consist of SU∗(|Γ|)/USp(|Γ|) (or SU(|Γ|)/USp(|Γ|), if adapted to the corresponding SO(N) (or USp(2N)) action. instead of taking the Hamiltonians as generators, we consider Consider first the specific case of the cubic JT problem the corresponding time-evolution generators iH). A list of all T1 ⊗ (e ⊕ t2) [50], or equivalently the icosahedral T ⊗ h (the JT models carrying a Lie group action encompassed by this generalization to more complex models will be made later) study is given in table 1, along with their corresponding sym- [49]. The continuous symmetry of this system in the presence metric spaces, and examples of point groups giving rise to of degenerate couplings was originally investigated by O’Brien these models. [50, 59]. The JT active displacements are characterized by
7 J. Phys.: Condens. Matter () 333001 Topical Review
' G &L ' K
Figure 6. D3d, D4h epikernel and Ci kernel structures for the T1 ⊗ (eg ⊕ t2g) JT model. The first two have the same spectrum (in the presence of SO(3) symmetry) including a non-degenerate electronic ground-state and a doubly-degenerate excited state, while the latter is a lower symmetry structure with the same r1(Q) as the first two, but different r2(Q), reflecting the presence of three non-degenerate electronic states in its spectrum.
the vector Q =(Qθ, Q, Qxy, Qzx, Qyz) [3]. A global dia- and r2 are SO(3) invariants (no electronic SO(3) rotation batic basis for Hel is defined by Hel = span{|x , |y , |z } will correspond to nuclear motion normal to O(r1, r2)). It Q ∈ ( , ) (which may be thought of as px, py and pz atomic orbitals). follows that given a representative O r1 r2 , any ele- The traceless Hermitian electronic tensor operators are in the ment Q ∈ O(r1, r2) may be obtained by an SO(3) action symmetric part of T1 ⊗ T1 − A1 and will be labeled by Vα, on Q, this rotation in turn corresponding to the electronic α ∈{θ, , xy, zx, yz}, Hilbert space SO(3) transformation U(Q , Q) satisfying ⎛ ⎞ ⎛ √ ⎞ (Q, Q ) (Q) −1(Q, Q )= (Q ) 1 00 − 3 00 U HJT U HJT . Thus, the 5-dimen- 2 2 √ ⎝ 1 ⎠ ⎜ ⎟ sional APES of T ⊗ (e ⊕ t2) may be reduced into a 2-dimen- Vθ = 0 0 , V = ⎝ 0 3 0⎠ , 2 2 sional orbit space. 00−1 000 √ The remaining three vibrational degrees of free- 3 dom are expected to parametrize the internal space (V ) = − (δ δ + δ δ ) , ij ab 2 ia jb ja ib (19) of the orbits O(r1, r2). Thus, SO(3)-adapted coor- dinates for the molecular displacements are given where i, j ∈{x, y, z}, i = j, δab = 1 if a = b and 0 otherwise. Q = Q(r1, r2, β, α, φ), ∀ r1 ∈ R −{0}, r2 ∈ [−1, 1] The JT active distortions and the space of electronic by (see (β, α, φ) Hamiltonians are isomorphic to the space of quadrupoles [54] below), where parametrize SO(3) rotations. They can (the 2 irrep of SO(3)), which is also isomorphic to the tangent be understood as Euler angles in the space of JT active dis- space of SL(3, R)/SO(3) at the identity, i.e. sl(3, R)/so(3). torted structures [33, 60]. We employ the zyz convention for the parametrization of SO(3) elements, Its Cartan subalgebra is two-dimensional. It may be taken to −iφ −iα −iβ (β, α, φ)=e Jz e Jy e Jz , be for example, span{V, Vθ} or alternatively span{Vθ, Vxy}. U (20) Importantly, this implies there exists two SO(3) invariants in where U(β, α, φ) acts on the 1 (vector) irrep of SO(3), the Ji the universal enveloping algebra of sl(3, R)/so(3, R) [39]. are its Hermitian generators, (Ji)jk = i ijk, φ, β ∈ [0, 2π), and From the SO(3)-equivalence between symmetric electronic α ∈ [0, π]. tensor operators (minus the identity) and vibrational displace- In order to understand the internal space of each orbit, (Q) ments, it follows that there exists also two functions r1 note that the JT distortion-induced electronic multiplet split- (Q) and r2 of the nuclear degrees of freedom which are invari- ting can only happen in two different ways: either HJT(Q) ant under the action of SO(3). They will be denoted radial (with Q = 0) has three distinct eigenvalues or two degener- coordinates. Let r1 denote the quadratic SO(3) invariant of 2 ate and a single non-degenerate. The matrices representing the sl(3, R)/so(3), r1(Q)=|Q| = Q, and r2(Q) be the cubic. vibronic coupling due to t2 or e displacements demonstrate The latter may be obtained from the secular determinant of this point (equation (19)). HJT(Q) is related by an orthogo- (Q) 0 HJT [39], or as the totally symmetric SO(3) irrep in the nal transformation to an element of the Cartan subalgebra 2 ⊗ 2 ⊗ 2 decomposition of the tensor product , since the cubic of sl(3, R)/so(3) (section 2). In particular, if we choose Vθ invariant is a homogeneous polynomial of third degree in the and V to be a basis for the Cartan subalgebra, we find from 2 vibrational coordinates (which transform like the irrep of equation (17), SO(3)) [33]. We provide a simple derivation of r2(Q) below. (Q)= (β, α, φ) + θ −1(β, α, φ), This simple analysis indicates that each molecular geom- HJT FvU QcV Qc Vθ U (21) etry of T ⊗ (e ⊕ t2) belongs to a single SO(3) orbit charac- HJT(Q)=−Fv|Q|U(β, α, φ) terized by two radial coordinates, r1 and r2 (compare to the ⎛ ⎞ cos[γ(Q ) − 2π ] 00 E ⊗ e case which has orbits described by a single real num- c 3 × ⎝ 0 cos[γ(Q )+ 2π ] 0 ⎠ −1(β, α, φ), ber (figure4 , equations (9) and (10))). We denote a specific c 3 U 00cos[γ(Q )] orbit of this model by O(r1, r2). It only includes molecu- c (22) lar geometries with the same electronic spectrum, since r1
8 J. Phys.: Condens. Matter () 333001 Topical Review
where Fv denotes the vibronic coupling constant, electronically degenerate doublet iff γ = 0 or γ = π/3, i.e. Q =( θ, ,0,0,0) c Qc Qc defines an element of the Cartan sub- for r2 = ±1. For any other values of γ ∈ [0, π/3], the corre- algebra of sl(3, R)/so(3) specifying the eigenvalues of the JT sponding electronic JT Hamiltonian belongs to a 3D orbit of −1 Hamiltonian for a given geometry Q = R Qc, R ∈ SO(3), SO(3). γ(Q)=tan−1( / θ) and Qc Qc . These considerations drasti- Our analysis also reveals the topology of each type cally simplify the task of finding the coordinate mapping of orbit. The 2D orbits are parametrized by a 3D rota- 2 Q → r2(Q) (the cubic invariant of sl(3, R)/so(3)). It is pro- tion axis. Thus, each 2D orbit is isomorphic to RP . The 3 ∼ 3 portional to the O(λ ) term of det[HJT(Q) − Iλ] [39], which 3D orbits are copies of SO(3)/V = RP /V [58], where may now be easily calculated from equation (22), V = {U(0, 0, 0), U(0, 0, π), U(0, π,0) U(π/2, π, π/2)}. The latter follows from the fact that all elements of V are diagonal r2(Q) = cos[3γ(Qc)]. (23) matrices, and both Vε and Vθ are left invariant under conjuga- For any Q, HJT(Q) is related by an orthogonal transformation tion by an element of the subgroup V. Therefore, electronic (Q ) γ ∈ [0, π/3] ∈ [−1, 1] to a unique HJT c with [54], so r2 . Hilbert space transformations by elements of V correspond to According to the above, motion along the Qθ and Q direc- trivial action on the vibrational configuration space. tions provides all distinct possibilities for the splitting of the The qualitative analysis carried to this point indicates the = 0 electronic multiplet when Q . This may also be seen by APES of the SO(3)-invariant T ⊗ (e ⊕ t2) can be decomposed checking that the eigenvalues of Vij are equal to those of V, into a 2D orbit space. Only two topologically non-equivalent which in turn are different from those of Vθ. The orbit obtained types of orbits exist in agreement with the two linearly inde- by SO(3) action on Vθ defines a 2D subspace of the vibrational pendent ways of lifting the degeneracy of triply-degenerate [ , ]=0 configuration space, since Jz Vθ . In more detail, suppose electronic states. Q = e =( ,0,0,0,0) Q = Q Q θ Q . Then, c , and By applying equations (19)–(21) to (22), we obtain −1 the SO(3)-adapted parametrization of JT distortions [50], U(0, 0, φ)QVθU (0, 0, φ)=QVθ. (24) Q = Q(Q, γ, β, α, φ),