Scientific Research

An Academic Publisher

**Pseudo Jahn-Teller Effect in Puckering and Planarization of Heterocyclic Compounds** ()

_{12}H

_{8}S

_{2 }molecule up to the dianion (C

_{12}H

_{8}S

_{2})

^{2-}enhances the PJTE, followed by the S-S bond cleavage and significant structural distortions of the system. The change of the PJTE by chemical substitutions, accompanied either by puckering or by planarization of heterocyclic compounds, is discussed using as examples 1,4-ditinine and its S-oxygenated derivatives.

Keywords

Share and Cite:

*International Journal of Organic Chemistry*,

**8**, 142-159. doi: 10.4236/ijoc.2018.81010.

1. Introduction

From general space-symmetry considerations, one would expect that in an isotropic environment any molecular system should have the highest possible symmetry, but in fact we see that in many cases the symmetry of the system is much lower (see [1] ). In this respect, cyclic molecules (including the heterocyclic ones) should be at least planar, and any deviation from the planar structure (puckering) can be described through the instability of the reference planar high-symmetry nuclear configuration. A general approach to handle instabilities and structural changes in molecular systems in non-degenerate states is the pseudo Jahn-Teller effect (PJTE). It was proved that the PJTE is the only source of instability of high-symmetry configurations of molecules or solid in non-degenerate states (See, e.g., the monographs [2] [3] , and the reviews [1] [4] [5] and references therein). The phenomenon of puckering as due to the PJTE was first formulated in [3] , and then employed to explain the origin of the out-of-plane distortions of a variety of cyclic and heterocyclic systems [6] - [11] . Knowledge of the mechanism of puckering of cyclic molecules provides the key to understanding how one can influence the system to suppress this effect and to flatten the system. In a recent review of Bersuker [12] devoted to manipulation of structure of two-dimensional (2D) or quasi 2D systems, the methods of targeted external influence that suppresses the PJTE are demonstrated on a series of molecular and extended 2D systems including silicene, phosphorene, germanene, graphitic carbon nitride, etc.

In this paper, we first summarize the results of studies published over the last few years regarding the PJTE origin of the puckered structures of heterocyclic compounds. Then, the possibility of enhancement and suppression of the pseudo Jahn-Teller effect by oxidation, reduction and chemical substitutions and therefore the possibility of changing the structure of the systems is discussed in more detail using a series of heterocyclic compounds as examples.

Chart 1 lists the heterocyclic systems which were studied recently with regard to the PJTE in the origin of their puckering. The figures in the first column of Chart 1 show schematically the initial planar and equilibrium configurations of the compounds under study, with the indication of the symmetry of the configurations and the type of distortions. The second column indicates the PJTE problem, which was solved in the corresponding work.

The first group of compounds includes 1,4-dithiine molecules C_{4}S_{2}L_{4} [10] with a variety of the ligands (L = H, F, Cl, Br), S-oxygenated derivatives of 1,4-dithiine, two tricyclic systems (thianthrene [11] and C_{6}S_{8} with two S-atoms in the 1,4-position of the central six-membered ring [13] ), and 1,4-dihydropyrazine and its derivatives C_{4}N_{2}H_{4}L_{2} [15] in which two imide hydrogen atoms are replaced with halogens (L = F, Cl, Br). All these compounds are unstable in the high-symmetry planar nuclear configuration of D_{2h} symmetry, in the equilibrium configuration of C_{2v} symmetry they are bent at the S-S axis. The exception is the tricyclic C_{6}S_{8} molecule, which has C_{2h} symmetry in the planar configuration and C_{2} symmetry in the equilibrium.

The next group of molecules includes 1,2-dichalcogenins, 1,2-C_{4}X_{2}L_{4}, with X = O, S, Se, Te and L = H, F [10] , two tricyclic species C_{6}S_{8} and
${\text{C}}_{\text{6}}{\text{S}}_{8}^{2-}$ containing two S-atoms in the 1,2-positions of the central six-membered ring and one

Chart 1. List of heterocyclic systems studied with respect to the PJTE origin of their puckering.

thione (C=S) bond in the five-membered rings on its either side [11] , two tricyclic molecules C_{12}H_{8}S_{2} and (C_{12}H_{8}S_{2})^{2−} [this work], 3,6-pyridazinedione derivatives N_{2}C_{4}H_{2}Y_{2}Z_{2} (Y = O, S, Se; Z = H, F, Cl, Br) [16] , and 1,2-diazetes C_{2}N_{2}L_{4} with L = H, F, Cl, Br [17] . The planar configuration of these systems is unstable with respect to the twisting coordinate of a_{2} symmetry that leads to a puckered equilibrium structure of C_{2} symmetry.

Finally, one can distinguish the third group of five-membered heterocyclic molecules of the type C_{2}X_{3}Y_{2} (X = O, S, Se, Te; Y = H, F) [18] and C_{4}XY_{5} (X = N, P, As; Y = H, F, Cl) [19] , and the seven-membered rings C_{6}Y_{6}O and C_{6}Y_{7}N with Y = H, F, Cl, Br [20] . As in the previous cases, in these compounds the planar nuclear configuration of C_{2v} symmetry is distorted along the out-of-plane b_{1} coordinate of instability that leads to equilibrium puckered structure of Cs symmetry.

In the next Section 2, we give the basic formulas of the PJTE theory, which will be needed in discussing the results. In Section 3, the possibility of flattening the 1,4-dithiine containing compounds by suppressing the PJTE via oxidation is discussed. The enhancement of the PJTE induced distortion by reduction is analyzed in Section 4 using the dibenzo[1,2]dithiine molecule as an example. In the last Section 5, the change of the PJTE by chemical substitution, accompanied by a change in the structure of heterocyclic compounds, is discussed using, as examples, 1,4-ditinine and its S-oxygenated derivatives. Geometry optimization and vibration frequency analysis for all the molecules considered in Sections 4 and 5 were performed at the B3LYP level of the DFT method [21] . The Pople’s 6-31+G(d,p) split basis sets [22] were utilized in all steps of the calculations. The potential energy profiles along the distortion coordinates were calculated with the CISD method. The applied active space for the C_{12}H_{8}S_{2} and (C_{12}H_{8}S_{2})^{2−} species was composed of eight electrons and nine active orbitals, while that for the 1,4-ditinine S-oxygenated derivatives included ten electrons and ten active states. All the calculations were carried out using the GAUSSIAN 09 program package [23] .The numerical values of the vibronic coupling constants were estimated by means of fitting the solutions of the secular equations to the ab initio calculated energy profiles.

2. Basic Formulas of the PJTE Theory

The theory of the PJTE is well developed (see, for example, [1] [2] [3] [4] ). In this approach the problem of the stability or instability of molecular nuclear configuration is reduced to estimation of the curvature K of the adiabatic potential energy surface (APES) in direction of any normal displacement Q, at the high-symmetry configuration Q_{0} for which the first derivatives are zero. At the minimum of the APES
$K={\left({\partial}^{2}E\left(Q\right)/\partial {Q}^{2}\right)}_{0}>0$ , while K < 0 if the high-symmetry configuration is unstable.

The exact expression for K of any molecular system in nondegenerate ground state Ψ_{0} can be obtained in the second order perturbation theory with respect to small nuclear displacements Q:

$K=\langle {\Psi}_{0}|{\left({\partial}^{2}H/\partial {Q}^{2}\right)}_{0}|{\Psi}_{0}\rangle -2\underset{i}{{\displaystyle \sum}}\frac{{\left|\langle {\Psi}_{0}|{\left(\partial H/\partial Q\right)}_{0}|{\Psi}_{i}\rangle \right|}^{2}}{\left({E}_{i}-{E}_{0}\right)}.$ (1)

The first term in Equation (1), K_{0}, is the so-called primary force constant. It determines the restoring force arising when the nuclei are displaced with respect to the “frozen” electron distribution. The second term, K_{v},

${K}_{v}=-2\underset{i}{{\displaystyle \sum}}\frac{{F}_{0i}^{2}}{{E}_{i}-{E}_{0}}$ (2)

is the vibronic contribution to the curvature coming from the vibronic interaction of the ground Ψ_{0} and excited Ψ_{i} states under the nuclear displacements, E_{0} and E_{i} are their energies. For the ground state this term is always negative. The matrix elements F_{0i} in Equations (1) and (2),

${F}_{0i}=\langle {\Psi}_{0}|{\left(\partial H/\partial Q\right)}_{0}|{\Psi}_{i}\rangle ,$ (3)

are the off-diagonal vibronic coupling constants. They are non-zero only if the product of irreducible representations of the wave functions of the ground and excited states ${\Gamma}_{{\Psi}_{0}}\times {\Gamma}_{{\Psi}_{i}}$ contains the irreducible representation ${\Gamma}_{Q}$ of the displacement Q.

It was proved analytically and confirmed by a series of numerical calculations (see in [3] [4] [5] ) that for any molecular system in the high-symmetry configuration K_{0} > 0. This means that the negative value of the curvature and, therefore, structural instabilities and distortions of high-symmetry nuclear configurations are only due to the vibronic contribution K_{v}. The instability takes place if the inequality

$\left|{K}_{v}\right|>{K}_{0}$ or $\Delta <2{F}^{2}/{K}_{0}$ (4)

holds.

In the simplest two-level problem, when only one nondegenerate excited state contributes significantly to the instability of the ground state, we get the following secular equation for the two states (the energy is read of the ground state level):

$\left|\begin{array}{cc}\left(1/2\right){K}_{0}{Q}^{2}-\epsilon & FQ\\ FQ& \left(1/2\right){K}_{1}{Q}^{2}-\epsilon +\Delta \end{array}\right|=0$ (5)

where K_{0} and K_{1} are the primary force constants in the ground and excited states, respectively.

The solution of this 2 × 2 equation is straightforward:

${\epsilon}_{1,2}=\frac{1}{4}\left({K}_{0}+{K}_{1}\right){Q}^{2}+\frac{\Delta}{2}\pm \frac{1}{2}\sqrt{{\left[\frac{1}{2}\left({K}_{0}-{K}_{1}\right){Q}^{2}-\Delta \right]}^{2}+4{F}^{2}{Q}^{2}}$ (6)

Direct calculation of the vibronic constants F and K_{i} involved in this PJT model is rather difficult mathematically. Usually, the values of the parameters are estimated by fitting the solutions of the secular equations to the ab initio calculated energy profiles along the coordinate of instability. In the more complex case of the three-level problem, when two excited states are active in the PJT mixing, the matrix equation is

$\begin{array}{l}\left|\begin{array}{ccc}\frac{1}{2}{K}_{0}{Q}^{2}-\epsilon & {F}_{01}Q& {F}_{02}Q\\ {F}_{01}Q& \frac{1}{2}{K}_{1}{Q}^{2}+{\Delta}_{01}-\epsilon & 0\\ {F}_{02}Q& 0& \frac{1}{2}{K}_{2}{Q}^{2}+{\Delta}_{02}-\epsilon \end{array}\right|\\ =\left|\begin{array}{ccc}a-\epsilon & f& g\\ f& b-\epsilon & 0\\ g& 0& c-\epsilon \end{array}\right|=0\end{array}$ (7)

Here K_{0}, K_{1}, and K_{2} are the primary force constants for the ground and two excited states, respectively, Δ_{01}, and Δ_{02} are the energy gaps between the ground and the two excited states, and F_{01} and F_{02} are the corresponding vibronic coupling constants,

${F}_{01}=\langle 0|{\left(\partial H/\partial Q\right)}_{0}|1\rangle $ and ${F}_{02}=\langle {0}_{1}|{\left(\partial H/\partial Q\right)}_{0}|2\rangle $ . (8)

The constants a, b, c, f, and g are simplifying denotations. Solving for the 3×3 secular determinant, one can get the following equation for the energies of the three states (1 ground state and 2 excited states) along the distortion coordinate:

${\epsilon}^{3}-\left(a+b+c\right){\epsilon}^{2}+\left(ab+ac+bc-{f}^{2}-{g}^{2}\right)\epsilon -abc+c{f}^{2}+b{g}^{2}=0$ (9)

A similar equation can also be obtained for the four-level PJT problem. The procedure for estimating the vibronic constants in these cases can be found in the works [10] [13] .

Note that Equations (5)-(7) are correct for small Q values only; they characterize the changes in the K values induced by the term interactions under the Q displacements.

For a qualitative interpretation of phenomena related to the PJTE, it is worthwhile to use a more visual orbital scheme that provides information on the role of individual orbitals in the origin of instability. First of all we note that in Equation (1) ∂H/∂Q = ∂V/∂Q where V is the part of the Hamiltonian that depends on the nuclear coordinates

$V=-\underset{\alpha =1}{\overset{N}{{\displaystyle \sum}}}\underset{i=1}{\overset{n}{{\displaystyle \sum}}}\frac{{Z}_{\alpha}}{\left|{r}_{i}-{R}_{\alpha}\right|}+\frac{1}{2}\underset{\alpha >\beta}{\overset{N}{{\displaystyle \sum}}}\frac{{Z}_{\alpha}{Z}_{\beta}}{\left|{R}_{\alpha}-{R}_{\beta}\right|}$ (10)

Here N is the number of atoms, Z_{α} are the nuclear charges, and R_{α} and r_{i} are the nuclear and electron coordinates respectively. Further, assume that the wavefunction of the excited state Ψ_{i} is the Slater determinant which differs from that of the ground state just by one-electron excitation α → m from the double occupied MO |φ_{α}ñ to the unoccupied MO |φ_{m}ñ. Then, the vibronic coupling constant F_{0i} in Equation (3) can be simply expressed by means of the off-diagonal orbital vibronic coupling constant (OVCC)
${f}_{Q}^{\alpha m}$ :

${F}_{0i}=\langle {\Psi}_{0}|{\left(\partial H/\partial Q\right)}_{0}|{\Psi}_{i}\rangle =\sqrt{2}{f}^{\alpha m}$ , ${f}^{\alpha m}=\langle {\phi}_{\alpha}|{\left(\partial V/\partial Q\right)}_{0}|{\phi}_{m}\rangle $ , (11)

and
${f}_{Q}^{\alpha m}$ are nonzero only if the product of symmetries of the two MOs |φ_{α}ñ and |φ_{m}ñ contains the symmetry of distortion Q. If |φ_{α}ñ is the single occupied MO then

${F}_{0i}=\langle {\Psi}_{0}|{\left(\partial H/\partial Q\right)}_{0}|{\Psi}_{i}\rangle ={f}^{\alpha m}$ (12)

Thus, a complicated PJT problem involving many excited states can be reduced to the study of the orbital pairs which are mixed by distortion thereby destabilizing the ground state.

3. Suppression of the PJTE by Oxidation. Possibility of Planarization of 1,4-Dithiine Containing Compounds

The reference nuclear configuration of all the considered systems from this group of compounds (see Introduction) is a planar structure of D_{2h} symmetry. In this configuration they have one imaginary frequency corresponding to the out-of-plane distortion of b_{1u} symmetry that leads to the C_{2v} butterfly-like equilibrium geometry [10] [11] [15] . This means that the planar structure is unstable with respect to the symmetrized b_{1u} displacements (note that this distortion is of b_{1u}-type if the symmetry axis of the second order is perpendicular to the plane of the molecule [11] , but if it lies in the molecular plane, then the same distortion is qualified as of b_{3u}-type [10] [15] ). The ground electronic state of all the molecules is the ^{1}A_{g} one. According to the PJTE (Equation (3)) the excited states that cause the instability of the ground state should have the ^{1}B_{1u} symmetry (A_{g} × B_{1u} = b_{1u}). In all the systems there is only one low-lying ^{1}B_{1u} excited state. Therefore, in the studied compounds we have a two-level (A_{g} + B_{1u}) Ä b_{1u} PJTE problem. An exception is the C_{6}S_{8} molecule in which three excited states of A_{u} symmetry contribute to the instability of the ground state trough the four-level (A_{g} + 1A_{u} + 2A_{u} + 3A_{u}) Ä a_{u} PJTE problem [13] .

Then the energy profiles (cross-sections of the APES) along the coordinate of instability b_{1u} were calculated for the ground and the B_{1u} excited states. The vibronic coupling constants were estimated by means of fitting the solutions of the 2 × 2 secular equations (Equation (6)) to the corresponding energy profiles. Thus obtained numerical values of the PJTE parameters allow one to verify the condition of PJTE instability (4). Indeed, for small values of Q the expression for the energy can be approximately written as

${\epsilon}_{gr}=\frac{1}{2}K{Q}^{2}=\frac{1}{2}\left({K}_{0}-\frac{2{F}^{2}}{\Delta}\right){Q}^{2}$ (13)

Substituting in this formula the numerical data from the works [10] [11] [15] , one can see that in all the cases the resulting values of the curvature K is negative. Thus, we can conclude that the instability of the high-symmetry D_{2h} nuclear configuration of the 1,4-dithiins with respect to the b_{1u} puckering modes and their distortions are due to the PJT mixing of the ground ^{1}A_{g} and the first excited ^{1}B_{1u} states.

In all the considered 1,4-dithiins these ^{1}B_{1u} excited states are mainly formed by one-electron excitations from the HOMO b_{1u} to the unoccupied a_{g} MO. Then the removal of an electron from this orbital reduces by half the PJTE negative contribution K_{v} to the curvature of the APES compared with a neutral molecule. Indeed, as follows from Equations (11) (12), the vibronic contribution to the curvature for the neutral molecules is approximately equal to
${K}_{v}\approx -\text{4}{f}^{\text{2}}/\Delta $ , whereas for the radical cations
${K}_{v}\approx -2{f}^{\text{2}}/\Delta $ where
$f=\langle {b}_{1u}|{\left(\partial H/\partial {Q}_{b1u}\right)}_{0}|{a}_{g}\rangle $ .

The decrease in the vibronic contribution leads to a positive value for the curvature K of the APES, indicating that the cations of 1,4-dithiine containing compounds should have a planar nuclear configurations. To verify this statement, we have performed calculations of the oxidized molecules, which have shown the absence of any imaginary frequency in the ground state of their planar configuration, that is, the planar geometry of the cations is an equilibrium one.

Thus, one can conclude that bending of all the considered 1,4-dithiine containing molecules is due to the PJT coupling between the ground A_{g} and the excited B_{1u} states and restoration of their planar nuclear configurations upon oxidation is directly related to the decrease in the orbital vibronic coupling between these states, and hence, suppression of the PJTE.

4. Enhancement of the PJTE Distortions by Reduction. S-S Bond Cleavage in 1,2-Dithiine Containing Molecules

In [24] 3,8-diiodo-dibenzo[1,2]dithiine was investigated both theoretically and experimentally as a potential molecular system which can act as a memory storage bank. It was shown that the neutral molecule is slightly twisted around the central axis whereas two-electron reduction of this compound causes S-S bond cleavage and significant structural rearrangement [24] . Earlier we have studied two tricyclic species, the neutral carbon sulfide C_{6}S_{8} and its dianion C_{6}S_{8}^{2−} containing two S-atoms in the 1,2-positions of the central six-membered ring and one thione (C=S) bond in the five-membered rings on its either side. We have shown [11] that in both C_{6}S_{8} and
${\text{C}}_{\text{6}}{\text{S}}_{8}^{2-}$ systems the out-of-plane distortions are due to the PJT coupling between the ground ^{1}A_{1} and two excited electronic states of ^{1}A_{2} symmetry and that reduction of C_{6}S_{8} up to
${\text{C}}_{\text{6}}{\text{S}}_{8}^{2-}$ leads to enhancement of the PJTE and as a result, to a much stronger distortion of the dianion compared with the neutral molecule. In the present paper we present the results of the electronic structure and PJTE investigation of dibenzo[1,2]dithiine molecule, C_{12}H_{8}S_{2}.

For this purpose the electronic structure calculations and vibrational frequency analysis for C_{12}H_{8}S_{2} and (C_{12}H_{8}S_{2})^{2−} molecules were carried out in their reference planar nuclear configuration of C_{2v} symmetry. In both cases, the calculations have shown the presence of one imaginary frequency of a_{2} symmetry with the values of 146.3 cm^{−1} and 236.7 cm^{−1} for C_{12}H_{8}S_{2} and (C_{12}H_{8}S_{2})^{2−}, respectively. This means that the planar configuration is unstable with respect to the symmetrized out-of-plane a_{2} displacement of the atoms. The calculations give also the directions and relative values of displacements (Figure 1 ).

Both systems have a nondegenerate ground electronic state ^{1}A_{1}. Therefore, only the excited states of ^{1}A_{2} symmetry can produce the instability of the ground state (Equation (3)). In Figure 2 the calculated energy profiles (cross sections of APES) for the ground ^{1}A_{1} and the low-lying excited ^{1}A_{2} states of C_{12}H_{8}S_{2} and (C_{12}H_{8}S_{2})^{2−} molecules along the instability coordinate a_{2} are presented. In the neutral C_{12}H_{8}S_{2} compound the first 1^{1}A_{2} state (∆ = 6.36 eV) is formed by one-electron excitation
${a}_{\text{2}}\left(\text{HOMO}\right)\to {a}_{\text{1}}\left(\text{LUMO}+\text{3}\right)$ , while the 2^{1}A_{2} state corresponds to the excitation
${b}_{\text{1}}\left(\text{HOMO}-\text{2}\right)\to {b}_{\text{2}}\left(\text{LUMO}+\text{2}\right)$ . This 2^{1}A_{2} state gives only a slight negative contribution to the curvature of the APES firstly due to the large energy gap (8.05 eV) and, secondly, because of small value of the vibronic constant (F_{02} = 0.51 eV/Å). In the dianion (C_{12}H_{8}S_{2})^{2−} the first excited 1^{1}A_{2} state
${b}_{\text{1}}\left(\text{HOMO}\right)\to {b}_{\text{2}}\left(\text{LUMO}\right)$ does not contribute to the instability of the ground state. Although the energy gap is only 2.5 eV, the vibronic constant is very small (F_{01} = 0.23 eV/Å) due to the fact that the mixed orbitals are localized in different parts of the molecule (Figure 2). Therefore, in both cases we can assume that the two-level PJTE problem is adequate for describing the instability of the ground state.

By fitting the solutions of the secular Equations (5) to the ab initio calculated energy profiles along the coordinate of instability we can estimate the values of the PJTE parameters. As a result we obtain: for the neutral molecule K_{0} = 0.69 eV/Å^{2}, F_{01} = 1.56 eV/Å, ∆ = 6.36 eV, and according to Equation (13) the resulting value of the curvature K = −0.075 eV/Å^{2}; for the dianion K_{0} = 0.63 eV/Å^{2},

Figure 1. Out-of-plane a_{2} normal mode in planar C_{2v} configurations for (a) the neutral C_{12}H_{8}S_{2} and (b) the dianion form (C_{12}H_{8}S_{2})^{2−}; (c) and (d) are the optimized twisted C_{2}^{ }structures for C_{12}H_{8}S_{2} and (C_{12}H_{8}S_{2})^{2−}, respectively. The arrow lengths show expected relative (not absolute) displacements.

(a) (b)

Figure 2. Ab initio calculated energy profiles (cross sections of the APES) along the coordinate of instability a_{2} for the ground A_{1} and the excited A_{2} electronic terms of (a) the neutral C_{12}H_{8}S_{2} molecule and (b) its reduced form (C_{12}H_{8}S_{2})^{2−}.

F_{02} = 1.51 eV/Å, ∆ = 5.52 eV, and K = −0.196 eV/Å^{2}. We see that the absolute value of the curvature for the reduced form estimated via the PJTE is almost two and a half times larger than this value for the neutral molecule. Thus, the reduction of the C_{12}H_{8}S_{2} molecule enhances significantly the PJTE and, as a consequence, enhances the distortions induced by it.

It should be noted that reduction does not always lead to enhancement of the PJTE and to an increase of the PJTE-induced distortions. In those cases when the excited state which is active in the PJTE is formed by excitation to the lowest unoccupied molecular orbital, the population of the latter in the process of reduction leads to suppression of the PJT effect and to the planarization of the distorted molecule. One such example is a cyclic molecule P_{6} which is distorted in a free state but it becomes planar being coordinated in the triple-decker sandwich complexes CpMoP_{6}MoCp due to the transfer of the electron density from the metal to the lowest unoccupied MO of P_{6} molecule [7] .

5. Change in the PJTE by Chemical Substitutions. Puckering and Planarization of 1,4-Dithiine and Its S-Oxygenated Derivatives

In [11] , it was shown that chemical substitution of any atoms in heterocyclic compounds can lead both to suppression of the PJTE and planarization of heterocycles and to its enhancement resulting in the PJTE induced symmetry breaking. So, for example, the anthracene molecule (An) in its ground electronic state has a planar nuclear configuration of D_{2h} symmetry. In this configuration it does not have any low-lying excited states of B_{1u} symmetry which produce the PJTE instability with respect to the folding of the molecule. In comparison, the substituted 9,10-dihydroanthracene (9, 10-H_{2}An) has a folded geometry of C_{2v} symmetry. The addition of two hydrogen atoms in this system leads to the appearance of a new excited B_{1u} state, the vibronic mixing of which with the ground state results in the enhancement of the PJTE and the distortion [11] . On the other hand, substitution of carbon atoms at the 1, 4, 6, and 9 positions in thianthrene molecule by more electronegative nitrogen atoms suppresses the PJTE due to increasing the energy gap between the ground and the excited electronic state responsible for the distortion that leads to planarization of the C_{12}H_{8}S_{2}N_{4} molecule [11] .

A similar suppression of the PJTE was shown to take place in the S-oxygenated derivative of 1,4-dithiine, C_{4}H_{4}(SO_{2})_{2}, which are planar [13] . The absence of the butterfly-type puckered structure in this molecule is explained by a significant increase of the energy gap between the ground and the PJT active excited states and by very small value of the vibronic coupling constant. The condition of instability is not satisfied and the system remains planar [13] . At the same time, DFT calculations have shown that the S-oxygenated derivatives of 1,4-dithiine with one and two SO groups are puckered and have a shadow boat conformation [25] . Below we demonstrate that the puckering of these molecules and planarization of the systems with the SO_{2} groups are due to enhancement and suppression of the PJTE.

For this purpose we studied a series of six molecules: 1,4-dithiine and five its oxygenated derivatives. Starting with the planar configuration, ab initio calculations of the electronic structure and vibrational frequencies of the systems in their high-symmetry nuclear configurations were carried out. One can see from Figure 3 that molecules with one and two SO_{2} groups have a planar conformation (no imaginary frequency), while the other four compounds have an imaginary frequency of b_{1u}/b_{1} symmetry. This means that the planar nuclear configurations of these molecules are unstable and undergo symmetry breaking along the normal coordinates Q(b_{1u}/b_{1}) transforming the planar D_{2h}/C_{2v} structures to the bending C_{2v}/C_{2} ones. The calculations give also the directions and relative values of these displacements (Figure 3).

All the considered molecules have a nondegenerate ground electronic states, ^{1}A_{g} in systems with D_{2h} reference configuration and ^{1}A_{1} for molecules with C_{2v} planar nuclear configuration (Figure 4). According to Equation (3), the excited

Figure 3. Puckering b_{1u}/b_{1} normal modes in the series of 1,4-dithiine and S-oxygenated derivatives; the arrow lengths show expected relative (not absolute) displacements. Optimized structures (C_{2v}/C_{s} configurations) are shown in two projections. In C_{4}H_{4}SS and C_{4}H_{4}SO_{2}SO_{2} with the reference configuration of D_{2h} symmetry, the z-axis is perpendicular to the plane of the molecules. In other systems with planar reference configurations C_{2v}, the z axis passes through the sulfur atoms, and the molecules lie in the yz plane. The designation of molecules is the same as in [25] .

Figure 4. Relative energies of the ground and excited states (in eV) of 1,4-dithiine (C_{4}H_{4}SS) and its S-oxygenated derivatives obtained by ab initio calculations at CAS (10, 10)/6-31+G(d, p) level in the optimized D_{2h} (C_{2v}) geometry.

states which produce the instability of the ground state should have the B_{1u}/B_{1} symmetry. It can be seen from Figure 4 that the addition of oxygen atoms to the systems significantly changes the number and energies of these excited states.

In 1,4-dithiine molecule, the main vibronic contribution to the instability comes from only one low-lying excited B_{1u} state which is formed by one-electron excitation from the HOMO b_{1u} to the (LUMO + 4) a_{g} (Figure 4 and Figure 5). The values of the PJT parameters we obtained are: K_{0} = 0.37 eV/Å^{2}, ∆ = 6.77 eV, and F = 1.32 eV/Å. They agree rather well with those calculated in [10] . At small distortions, the resulting value of the curvature of the AP along the b_{1u} puckering coordinate is negative, K = −0.14 eV/Å^{2}. Hence the condition of puckering instability (Equation (4)) is fully satisfied.

Passing to the C_{4}H_{4}(SO_{2})S and C_{4}H_{4}(SO_{2})_{2} molecules, one can see that the energy gaps between the ground and the excited B_{1}/B_{1u} states significantly increase (Figure 4). The second reason for suppression of the PJTE is a significant decrease in the vibronic coupling constants (F = 0.76 eV/Å and F = 0.04 eV/Å for the first and the second molecules, respectively). From Equation (11) it follows that
$F=\sqrt{2}{f}^{{b}_{1}{a}_{1}}$ for the C_{4}H_{4}(SO_{2})S molecule and
$F=\sqrt{2}{f}^{{b}_{1u}{a}_{g}}$ for the C_{4}H_{4}(SO_{2})_{2} one. Since the highest occupied b_{1} MO in the first molecule and b_{1u} MO in the second one have a zero value at one and, respectively, at both sulfur atoms (Figure 5 and Figure 6), this leads to a decrease in the values of the

Figure 5. MO energy levels schemes for 1,4-dithiine and its S-oxygenated derivatives in the ground A_{g} (A_{1}) states in optimized D_{2h} (C_{2v}) geometry with indication of the one-electron excitations to the B_{1u} (B_{1}) terms.

Figure 6. Combinations of the occupied and unoccupied molecular orbitals for all the compounds, which give the excited B_{1u} (B_{1}) states. H, H − 1, H − 2, H − 3, H − 4, L, L + 1, L + 2, L + 3 and L + 4 represent the respective HOMO, HOMO − 1, HOMO − 2, HOMO − 3, HOMO − 4, LUMO, LUMO + 1, LUMO + 2, and LUMO + 3.

orbital vibronic constants ${f}^{{b}_{1}{a}_{1}}/{f}^{{b}_{1u}{a}_{1g}}$ . Both these factors (an increase of the energy gaps and a decrease of the vibronic coupling constants) lead to suppression of the PJTE, the condition of Equation (4) is not satisfied, so these molecules have planar conformations.

In the derivatives with one and two SO groups, oxygen atoms bring their 2p_{π} electrons to the π-electron system of the heterocycle, thereby providing the appearance of additional occupied b_{1} MOs (HOMO-1 in Figure 5 and Figure 6). This leads to the appearance of additional excited B_{1} states (Figure 4) which also contribute to the instability of the ground state of planar nuclear configuration. The PJTE problem in these cases becomes the three-level one (A_{1} + 1B_{1} + 2B_{1}) Ä b_{1}. Following the procedure outlined in Section 1, and using calculated APES cross-sections for these compounds, we can estimate the values of the PJT parameters and the resulting values of the curvature of the APES in these cases. Such, for the C_{4}H_{4}SOS molecule we obtain K_{0} = 0.44 eV/Å^{2}, F_{01} = 1.26 eV/Å, Δ_{01} = 5.14 eV, F_{02} = 1.39 eV/Å, Δ_{02} = 8.57 eV. At small values of Q_{b}_{1} the resulting value of the curvature of the APES is equal to
$K\approx {K}_{0}-2{F}_{01}^{2}/{\Delta}_{01}-2{F}_{02}^{2}/{\Delta}_{02}=-0.63$ eV/Å^{2}. For the C_{4}H_{4}SOSO molecule these values are: K_{0} = 0.56 eV/Å^{2}, F_{01} = 1.20 eV/Å, Δ_{01} = 4.16 eV, F_{02} = 1.48 eV/Å, Δ_{02} = 6.42 eV, and
$K\approx {K}_{0}-2{F}_{01}^{2}/{\Delta}_{01}-2{F}_{02}^{2}/{\Delta}_{02}=-0.81$ eV/Å^{2}. We see that the absolute values of the curvature of the APES for these molecules estimated via the PJTE are much larger than that for the 1,4-dithiine molecule. At last, in the C_{4}H_{4}SO_{2}SO molecule, in comparison with the C_{4}H_{4}(SO_{2})S one, the energy gap between the ground and the excited B_{1} states decreases, while the value of the vibronic coupling constant
$F=\langle {A}_{1}|{\left(\partial H/\partial Q\right)}_{0}|1{B}_{1}\rangle $ becomes quite large (F = 1.38 eV/Å) due to the fact that the mixed orbitals (HOMO b_{1} and (LUMO + 1) a_{1}) are localized on the OSC_{2} fragment of the molecule which undergoes the most significant distortion (see Figure 5 and Figure 6). This leads to the large negative contribution K_{v} and to negative value of the curvature of the APES, so the heterocycle is distorted.

Thus, the oxygenation of the 1,4-dithiine molecule can lead to both the enhancement of the PJTE accompanied by puckering of the heterocycles (if oxygen atoms donate their 2p_{π} electrons to the π-system of the heterocycle) and to the PJTE suppression and subsequent flattening of derivatives with one and two SO_{2} groups.

6. Conclusions

Analysis of the results of published works and the results obtained in this paper allows us to draw the following conclusions:

1) The instability of planar nuclear configurations of all the considered heterocyclic molecules and their out-of-plane distortions are due to the pseudo Jahn-Teller effect.

2) If the instability of the ground state is provided mainly by the excited state formed by one-electron excitation from the HOMO to the appropriate by symmetry unoccupied MO, the oxidation of the systems by removing electrons from this MO leads to the suppression of the PJTE. In the case of 1,4-dithiine containing molecules, this results in the restoration of planar configuration.

3) In 1,2-dithiin containing molecules, reduction leads to enhancement of the PJTE followed by S-S bond cleavage and significant structural rearrangements of the systems.

4) Changes of the PJTE in the series of 1,4-ditinin and its S-oxygenated derivatives are accompanied either by out-of-plane distortions of their heterocycles (if oxygen atoms donate their 2p_{π} electrons to the π-system of a heterocycle) or by their flattening in derivatives with one and two SO_{2} groups.

Conflicts of Interest

The authors declare no conflicts of interest.

[1] | Bersuker, I.B. (2016) Spontaneous Symmetry Breaking in Matter Induced by Degeneracy and Pseudodegeneracy. In: Rice, S.A. and Dinner, A.R., Eds., Advances in Chemical Physics, Vol. 160, Wiley, Hoboken, NJ, 159-208. |

[2] |
Bersuker, I.B. (2006) The Jahn-Teller Effect. Cambridge University Press, Cambridge, UK. https://doi.org/10.1017/CBO9780511524769 |

[3] |
Bersuker, I.B. and Polinger, V.Z. (1989) Vibronic Interactions in Molecules and Crystals. Springer-Verlag, Berlin-Heidelberg-New York. https://doi.org/10.1007/978-3-642-83479-0 |

[4] |
Bersuker, I.B. (2001) Modern Aspects of the Jahn-Teller Effect Theory and Applications to Molecular Problems. Chemical Reviews, 101, 1067-1114. https://doi.org/10.1021/cr0004411 |

[5] |
Bersuker, I.B. (2013) Pseudo Jahn-Teller Effect—A Two-State Paradigm in Formation, Deformation, and Transformation of Molecular Systems and Solids. Chemical Reviews, 113, 1351-1390. https://doi.org/10.1021/cr300279n |

[6] |
Liu., Y., Bersuker, I.B. and Boggs, J.E. (2013) Pseudo Jahn-Teller Origin of Puckering in C4H42 + , Si4H42 + , and C4F42 + Dications. Chemical Physics, 417, 26-29. https://doi.org/10.1016/j.chemphys.2013.02.033 |

[7] |
Ivanov, A.S., Bozhenko, K.V. and Boldyrev, A.I. (2012) On the Suppression Mechanism of the Pseudo-Jahn-Teller Effect in Middle E6 (E = P, As, Sb) Rings of Triple-Decker Sandwich Complexes. Inorganic Chemistry, 51, 8868-8872. https://doi.org/10.1021/ic300786w |

[8] |
Sergeeva, A.P. and Boldyrev, A.I. (2010) Flattening a Puckered Pentasilacyclopentadienide Ring by Suppression of the Pseudo Jahn-Teller Effect. Organometallics, 29, 3951-3954. https://doi.org/10.1021/om1006038 |

[9] |
Jose, D. and Datta, A. (2012) Understanding of the Buckling Distortions in Silicene. The Journal of Physical Chemistry C, 116, 24639-24648. https://doi.org/10.1021/jp3084716 |

[10] | Ilkhani, A.R., Hermoso, W. and Bersuker, I.B. (2015) Pseudo Jahn-Teller Origin of Instability of Planar Configurations of Hexa-Heterocycles. Application to Compounds with C4X2 Skeletons (X = O, S, Se, Te). Chemical Physics, 460, 75-82. |

[11] |
Gorinchoy, N.N. and Bersuker, I.B. (2017) Pseudo Jahn-Teller Effect in Control and Rationalization of Chemical Transformations in Two-Dimensional Compounds. Journal of Physics: Conference Series, 833, 012010. https://doi.org/10.1088/1742-6596/833/1/012010 |

[12] |
Bersuker, I.B. (2017) Manipulation of Structure and Properties of Two-Dimensional Systems Employing the Pseudo Jahn-Teller Effect. FlatChem, 6, 11-27. https://doi.org/10.1016/j.flatc.2017.10.001 |

[13] |
Pratik, S.M., Chowdhury, C., Bhattacharjee, R., Jahiruddin, Sk. and Datta, A. (2015) Pseudo Jahn-Teller Distortion for a Tricyclic Carbon Sulfide (C6S8) and Its Suppression in S-Oxygenated Dithiine (C4H4(SO2)2). Chemical Physics, 460, 101-105. https://doi.org/10.1016/j.chemphys.2015.04.009 |

[14] |
Pratik, S.M. and Datta, A. (2015) 1,4-Dithiine-Puckered in the Gas Phase But Planar in Crystals: Role of Cooperativity. The Journal of Physical Chemistry C, 119, 15770-15776. https://doi.org/10.1021/acs.jpcc.5b04908 |

[15] |
Hermoso, W., Ilkhani, A.R. and Bersuker, I.B. (2014) Pseudo Jahn-Teller Origin of Instability of Planar Configurations of Hexa-Heterocycles C4N2H4X2 (X = H, F, Cl, Br). Computational and Theoretical Chemistry, 1049, 109-114. https://doi.org/10.1016/j.comptc.2014.10.007 |

[16] |
Ilkhani, A.R. (2015) Pseudo Jahn-Teller Origin of Twisting in 3,6-Pyridazinedione Derivatives; N2C4H2Y2Z2 (Y = O, S, Se, Z = H, F, Cl, Br) Compounds. Journal of Theoretical and Computational Chemistry, 14, Article ID: 1550045. https://doi.org/10.1142/S0219633615500455 |

[17] |
Ilkhani, A.R., Gorinchoy, N.N. and Bersuker, I.B. (2015) Pseudo Jahn-Teller Effect in Distortion and Restoration of Planar Configurations of Tetra-Heterocyclic 1,2-diazetes C2N2E4, E = H, F, Cl, Br. Chemical Physics, 460, 106-110. https://doi.org/10.1016/j.chemphys.2015.07.015 |

[18] |
Ilkhani, A.R. (2017) The Symmetry Breaking Phenomenon in 1,2,3-Trioxolene and C2Y3Z2 (Z = O, S, Se, Te, Z = H, F) Compounds: A Pseudo Jahn-Teller Origin Study. Química Nova, 40, 491-495. https://doi.org/10.21577/0100-4042.20170029 |

[19] |
Ilkhani, A.R. and Monajjemi, M. (2015) The Pseudo Jahn-Teller Effect of Puckering in Pentatomic Unsaturated Rings C4AE5, A = N, P, As, E = H, F, Cl. Computational and Theoretical Chemistry, 1074, 19-25. https://doi.org/10.1016/j.comptc.2015.10.006 |

[20] |
Ilkhani, A.R. (2017) Pseudo Jahn-Teller Effect in Oxepin, Azepin, and Their Halogen Substituted Derivatives. Russian Journal of Physical Chemistry A, 9, 1743-1751. https://doi.org/10.1134/S0036024417090035 |

[21] |
Polly, R., Werner, H.J., Manby, F.R. and Knowles, P.J. (2004) Fast Hartree-Fock Theory using Local Density Fitting Approximations. Molecular Physics, 102, 2311-2321. https://doi.org/10.1080/0026897042000274801 |

[22] |
Hehre, W.J., Ditchfield, R. and Pople, J.A. (1972) Self-Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian-Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules. The Journal of Chemical Physics, 56, 2257-2261. https://doi.org/10.1063/1.1677527 |

[23] | Frisch, M.J., Trucks, G.W., Schlege, H.B., et al. (2009) Gaussian 09. Revision B.01, Gaussian, Inc., Wallingford. |

[24] |
Benniston, A.C., Allen, B.D., Harriman, A., Llarena, I., Rostron, J.P. and Stewart, B. (2009) Accessing Molecular Memory via a Disulfide Switch. New Journal of Chemistry, 33, 417-427. https://doi.org/10.1039/B814676D |

[25] |
Vessally, E. (2008) DFT Calculations on 1,4-dithiine and S-oxygenated Derivatives. Bulletin of the Chemical Society of Ethiopia, 22, 465-468. https://doi.org/10.4314/bcse.v22i3.61259 |

Copyright © 2020 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.