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₄S₂L₄ [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₆S₈ with two S-atoms in the 1,4-position of the central six-membered ring [13] ), and 1,4-dihydropyrazine and its derivatives C₄N₂H₄L₂ [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₆S₈ molecule, which has C_2h symmetry in the planar configuration and C₂ symmetry in the equilibrium.

The next group of molecules includes 1,2-dichalcogenins, 1,2-C₄X₂L₄, with X = O, S, Se, Te and L = H, F [10] , two tricyclic species C₆S₈ and C 6 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₁₂H₈S₂ and (C₁₂H₈S₂)²⁻ [this work], 3,6-pyridazinedione derivatives N₂C₄H₂Y₂Z₂ (Y = O, S, Se; Z = H, F, Cl, Br) [16] , and 1,2-diazetes C₂N₂L₄ with L = H, F, Cl, Br [17] . The planar configuration of these systems is unstable with respect to the twisting coordinate of a₂ symmetry that leads to a puckered equilibrium structure of C₂ symmetry.

Finally, one can distinguish the third group of five-membered heterocyclic molecules of the type C₂X₃Y₂ (X = O, S, Se, Te; Y = H, F) [18] and C₄XY₅ (X = N, P, As; Y = H, F, Cl) [19] , and the seven-membered rings C₆Y₆O and C₆Y₇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₁ 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₁₂H₈S₂ and (C₁₂H₈S₂)²⁻ 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.

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₀ for which the first derivatives are zero. At the minimum of the APES K = ( ∂ 2 E ( Q ) / ∂ Q 2 ) 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 Ψ₀ can be obtained in the second order perturbation theory with respect to small nuclear displacements Q:

K = 〈 Ψ 0 | ( ∂ 2 H / ∂ Q 2 ) 0 | Ψ 0 〉 − 2 ∑ i | 〈 Ψ 0 | ( ∂ H / ∂ Q ) 0 | Ψ i 〉 | 2 ( E i − E 0 ) . (1)

The first term in Equation (1), K₀, 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 ∑ i F 0 i 2 E i − E 0 (2)

is the vibronic contribution to the curvature coming from the vibronic interaction of the ground Ψ₀ and excited Ψ_i states under the nuclear displacements, E₀ 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 0 i = 〈 Ψ 0 | ( ∂ H / ∂ Q ) 0 | Ψ i 〉 , (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 Γ Ψ 0 × Γ Ψ i contains the irreducible representation Γ 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. 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

| K v | > K 0 or Δ < 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):

| ( 1 / 2 ) K 0 Q 2 − ε F Q F Q ( 1 / 2 ) K 1 Q 2 − ε + Δ | = 0 (5)

where K₀ and K₁ are the primary force constants in the ground and excited states, respectively.

The solution of this 2 × 2 equation is straightforward:

ε 1 , 2 = 1 4 ( K 0 + K 1 ) Q 2 + Δ 2 ± 1 2 [ 1 2 ( K 0 − K 1 ) Q 2 − Δ ] 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

| 1 2 K 0 Q 2 − ε F 01 Q F 02 Q F 01 Q 1 2 K 1 Q 2 + Δ 01 − ε 0 F 02 Q 0 1 2 K 2 Q 2 + Δ 02 − ε | = | a − ε f g f b − ε 0 g 0 c − ε | = 0 (7)

Here K₀, K₁, and K₂ are the primary force constants for the ground and two excited states, respectively, Δ₀₁, and Δ₀₂ are the energy gaps between the ground and the two excited states, and F₀₁ and F₀₂ are the corresponding vibronic coupling constants,

F 01 = 〈 0 | ( ∂ H / ∂ Q ) 0 | 1 〉 and F 02 = 〈 0 1 | ( ∂ H / ∂ Q ) 0 | 2 〉 . (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:

ε 3 − ( a + b + c ) ε 2 + ( a b + a c + b c − f 2 − g 2 ) ε − a b c + 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 = − ∑ α = 1 N ∑ i = 1 n Z α | r i − R α | + 1 2 ∑ α > β N Z α Z β | R α − R β | (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 α m :

F 0 i = 〈 Ψ 0 | ( ∂ H / ∂ Q ) 0 | Ψ i 〉 = 2 f α m , f α m = 〈 φ α | ( ∂ V / ∂ Q ) 0 | φ m 〉 , (11)

and f Q α 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 0 i = 〈 Ψ 0 | ( ∂ H / ∂ Q ) 0 | Ψ i 〉 = f α 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.

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 ¹A_g one. According to the PJTE (Equation (3)) the excited states that cause the instability of the ground state should have the ¹B_1u symmetry (A_g × B_1u = b_1u). In all the systems there is only one low-lying ¹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₆S₈ 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

ε g r = 1 2 K Q 2 = 1 2 ( K 0 − 2 F 2 Δ ) 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 ¹A_g and the first excited ¹B_1u states.

In all the considered 1,4-dithiins these ¹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 ≈ − 4 f 2 / Δ , whereas for the radical cations K v ≈ − 2 f 2 / Δ where f = 〈 b 1 u | ( ∂ H / ∂ Q b 1 u ) 0 | a g 〉 .

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.

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₆S₈ and its dianion C₆S₈²⁻ 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₆S₈ and C 6 S 8 2 − systems the out-of-plane distortions are due to the PJT coupling between the ground ¹A₁ and two excited electronic states of ¹A₂ symmetry and that reduction of C₆S₈ up to C 6 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₁₂H₈S₂.

For this purpose the electronic structure calculations and vibrational frequency analysis for C₁₂H₈S₂ and (C₁₂H₈S₂)²⁻ 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₂ symmetry with the values of 146.3 cm⁻¹ and 236.7 cm⁻¹ for C₁₂H₈S₂ and (C₁₂H₈S₂)²⁻, respectively. This means that the planar configuration is unstable with respect to the symmetrized out-of-plane a₂ 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 ¹A₁. Therefore, only the excited states of ¹A₂ 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 ¹A₁ and the low-lying excited ¹A₂ states of C₁₂H₈S₂ and (C₁₂H₈S₂)²⁻ molecules along the instability coordinate a₂ are presented. In the neutral C₁₂H₈S₂ compound the first 1¹A₂ state (∆ = 6.36 eV) is formed by one-electron excitation a 2 ( HOMO ) → a 1 ( LUMO + 3 ) , while the 2¹A₂ state corresponds to the excitation b 1 ( HOMO − 2 ) → b 2 ( LUMO + 2 ) . This 2¹A₂ 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₀₂ = 0.51 eV/Å). In the dianion (C₁₂H₈S₂)²⁻ the first excited 1¹A₂ state b 1 ( HOMO ) → b 2 ( LUMO ) 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₀₁ = 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.69 eV/Ų, F₀₁ = 1.56 eV/Å, ∆ = 6.36 eV, and according to Equation (13) the resulting value of the curvature K = −0.075 eV/Ų; for the dianion K₀ = 0.63 eV/Ų,

F₀₂ = 1.51 eV/Å, ∆ = 5.52 eV, and K = −0.196 eV/Ų. 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₁₂H₈S₂ 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₆ which is distorted in a free state but it becomes planar being coordinated in the triple-decker sandwich complexes CpMoP₆MoCp due to the transfer of the electron density from the metal to the lowest unoccupied MO of P₆ molecule [7] .

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₂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₁₂H₈S₂N₄ molecule [11] .

A similar suppression of the PJTE was shown to take place in the S-oxygenated derivative of 1,4-dithiine, C₄H₄(SO₂)₂, 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₂ 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₂ groups have a planar conformation (no imaginary frequency), while the other four compounds have an imaginary frequency of b_1u/b₁ 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₁) transforming the planar D_2h/C_2v structures to the bending C_2v/C₂ 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, ¹A_g in systems with D_2h reference configuration and ¹A₁ for molecules with C_2v planar nuclear configuration (Figure 4). According to Equation (3), the excited

states which produce the instability of the ground state should have the B_1u/B₁ 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.37 eV/Ų, ∆ = 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/Ų. Hence the condition of puckering instability (Equation (4)) is fully satisfied.

Passing to the C₄H₄(SO₂)S and C₄H₄(SO₂)₂ molecules, one can see that the energy gaps between the ground and the excited B₁/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 = 2 f b 1 a 1 for the C₄H₄(SO₂)S molecule and F = 2 f b 1 u a g for the C₄H₄(SO₂)₂ one. Since the highest occupied b₁ 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

orbital vibronic constants f b 1 a 1 / f b 1 u a 1 g . 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₁ MOs (HOMO-1 in Figure 5 and Figure 6). This leads to the appearance of additional excited B₁ 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₁ + 1B₁ + 2B₁) Ä b₁. 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₄H₄SOS molecule we obtain K₀ = 0.44 eV/Ų, F₀₁ = 1.26 eV/Å, Δ₀₁ = 5.14 eV, F₀₂ = 1.39 eV/Å, Δ₀₂ = 8.57 eV. At small values of Q_b1 the resulting value of the curvature of the APES is equal to K ≈ K 0 − 2 F 01 2 / Δ 01 − 2 F 02 2 / Δ 02 = − 0.63 eV/Ų. For the C₄H₄SOSO molecule these values are: K₀ = 0.56 eV/Ų, F₀₁ = 1.20 eV/Å, Δ₀₁ = 4.16 eV, F₀₂ = 1.48 eV/Å, Δ₀₂ = 6.42 eV, and K ≈ K 0 − 2 F 01 2 / Δ 01 − 2 F 02 2 / Δ 02 = − 0.81 eV/Ų. 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₄H₄SO₂SO molecule, in comparison with the C₄H₄(SO₂)S one, the energy gap between the ground and the excited B₁ states decreases, while the value of the vibronic coupling constant F = 〈 A 1 | ( ∂ H / ∂ Q ) 0 | 1 B 1 〉 becomes quite large (F = 1.38 eV/Å) due to the fact that the mixed orbitals (HOMO b₁ and (LUMO + 1) a₁) are localized on the OSC₂ 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₂ groups.

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₂ groups.

Gorinchoy, N. (2018) Pseudo Jahn-Teller Effect in Puckering and Planarization of Heterocyclic Compounds. International Journal of Organic Chemistry, 8, 142-159. https://doi.org/10.4236/ijoc.2018.81010