Theoretical Model for a Novel Electronic State in a Dirac Electron System Close to Merging: An Imaginary Element between Sulphur and Selenium

Abstract

Topological materials with Dirac electron systems have been extensively studied. Organic crystalline materials form a unique group of such compounds with well-defined crystal structures. While most organic compounds require high pressures to exhibit Dirac-cone-type band structures, the title compound, α-STF₂I₃, has garnered increasing interest due to its Dirac-cone-type band structure under ambient pressure. Various experiments have been conducted under ambient pressure; their results can be compared with those of theoretical calculations to obtain insights into Dirac electron systems. However, structural disorder peculiar to the STF molecules in the solid-state has prevented any type of theoretical calculation of the states. In this study, we report a new method for calculating intermolecular interactions in disordered systems based on the extended Hückel approximation. This method enables band calculations, suggesting that this material is a rare example of a system close to merging. The obtained band structure indicates that the characteristic disorder in the STF solids distributed electrons equally on the sulphur and selenium atoms as if they belong to an imaginary element between sulphur and selenium and are arranged without disorder.

1. Introduction

Recently, Dirac electron systems (DESs), which are characterised by massless relativistic electrons with speeds of 1/100–1/1000th of the velocity of light, have been widely studied. The majority of the studies on DESs are theoretical because of the limited availability of DES materials [1,2,3,4,5,6,7,8,9]. Concerning the general molecular systems, the electronic structures including those under high pressure are calculated by DFT methods [10,11,12]. Although DESs were initially found in graphene [13,14] and some inorganic compounds [15], those found in organic compounds [16] possess clear advantages over their inorganic counterparts. For example, most organic DESs are found in bulk systems with well-defined crystal structures and chemical stoichiometries, whereas inorganic DESs often lack either of these two features. However, organic DESs have some disadvantages; for example, most become DES only under high pressure [17,18,19,20], unlike inorganic DESs. Performing detailed, accurate experiments and measurements (including crystal structure analyses) at high pressures is more complicated than those at ambient pressure. In this regard, the organic charge-transfer (CT) complex, α-STF₂I₃ (Figure 1; STF = bis(ethylenedithio)-diselenadithiafulvalene), is unique as it contains Dirac electrons at ambient pressure [21,22,23,24,25,26]; that is, it possesses the advantages of both organic and inorganic DESs.

All the electrical, magnetic, and optical properties of the DESs are governed by the intermolecular interactions between the STF molecules in α-STF₂I₃. In principle, intermolecular interactions, band structures, and various physical properties of crystalline materials can be calculated when all the atomic positions, namely atomic parameters, are known, assuming perfect three-dimensional periodic arrangements of the atoms in the framework of well-established band theories. However, it is challenging to calculate the intermolecular interactions, that is, the overlap and transfer integrals between the STF molecules in α-STF₂I₃, because all the STF molecules in this CT complex are randomly disordered between the two patterns shown in Figure 2. The calculations of the electronic structures and related quantities in the solid states require well-defined atomic parameters for all the atoms involved. Although DES materials with the disorder are seldom reported and have garnered considerable attention [27,28,29], the lack of the abovementioned calculation results has seriously hindered further studies on this CT complex and a clear understanding of the general DESs.

We have recently proposed a tight-binding band structure of α-STF₂I₃ based on the overlapping integrals S_STF, assuming statistically averaged structures between all the possible molecular arrangements at each site in the disordered crystal (Figure 3; see Appendix A for details) in addition to the assumption of t_STF/eV = −10S_STF, where t_STF represents the corresponding transfer integrals [21]. The band structure obtained in our previous study was characterised by less tilted and more isotropic Dirac cones compared to that in the present study.

The calculated band structure failed to reproduce the observed anisotropic temperature dependence of resistivity, despite being qualitatively consistent with the nearly temperature-independent electrical resistivity. In this paper, we propose a new method based on molecular orbitals (MOs) to estimate the overlap and transfer integrals in such disordered systems more clearly. More importantly, we found that solid-state STF molecules form unique electron configurations, which are qualitatively different from those in other disordered crystalline materials. The differences between the STF solids and other disordered systems include observed qualitative differences in electrical, magnetic, spectroscopic, and structural properties. The STF solids almost always behave as if there is no disorder in their crystal structures. In contrast, the slightest disorder frequently alters qualitative behaviours compared with those of related or corresponding materials without the disorder in common crystalline materials. To explain this fact, we propose a new idea of wavefunctions analogous to molecular orbitals and based on the original meanings of wavefunctions.

2. Calculation Methods

Single-crystal X-ray structural analyses demonstrated that all four chalcogen atoms (X: Figure 2) in the solid-state STF molecules possessed equal electron densities [21,22,23,24,25,26]. This trend is not only the case with conducting STF CT complexes, but also the case with insulating STF CT complexes and even with the neutral STF molecular crystals, irrespective of the crystal structures and the molecular symmetries required for the positions in the unit cells [26]. This trend is peculiar to the STF molecule in the solid-state, which cannot be explained by the known quantum stabilisation effects, such as delocalisation or intermolecular interactions. Such electron configurations require that all the atomic positions and the electron densities on each X atom are fixed in each molecule, regardless of the packing pattern. Thus, only the degrees of freedom remain in the assignment of S and Se atoms at fixed X positions with equal electron densities. In other words, the electron densities were averaged for all four X atoms as if the X atoms belonged to imaginary element between sulphur and selenium. This is probably because the differences in the atomic energy levels are sufficiently small between S and Se atoms, favouring such averaged electron configuration. Ultimately, this configuration would lead to the reduction of Coulombic repulsion between electrons, in addition to stabilisation due to the increased kinetic energy of electrons. As wavefunctions for electrons are associated with electron densities, we assumed that the wavefunction of an STF molecule (Equation (1)) in the solid state (ϕ_STF (MO of STF)) should be described using the two wavefunctions corresponding to the average of the two patterns in Figure 2.

$${\varphi }_{\mathrm{STF}}=\frac{1}{\sqrt{2}}{\varphi }_1+\frac{1}{\sqrt{2}}{\varphi }_2$$

(1)

where suffixes 1 and 2 denote the different patterns shown in Figure 2. To equally distribute the electrons on the X atoms in the STF molecules, both coefficients must be equal, that is, $1/\sqrt{2}$ if Equation (1) is normalised and if the cross-terms are ignored. The cross-terms between $\varphi$₁ and $\varphi$₂ are ignored because the two orientations ($\varphi$₁ and $\varphi$₂) never coincide at a given site for a single STF molecule. In other words, such an averaged electron configuration does not involve the exchange or oscillating movement between the S and Se nuclei in the STF molecules. Equation (1) describes the electronic configuration of the X atoms in Figure 2, corresponding to the imaginary element between S and Se. Notably, it appears that Equation (1) describes the interaction between the two states or the fluctuation originating from the quantum interference between the two states. However, this is not the case. As we will see by substituting atomic orbitals in Equation (1), a “conjugate” equation, that is, a state corresponding to Equation (2) does not describe the electronic configuration of the solid STF molecules (Equation (6)).

$${{\varphi }_{\mathrm{STF}}}^{\prime }=\frac{1}{\sqrt{2}}{\varphi }_1-\frac{1}{\sqrt{2}}{\varphi }_2$$

(2)

All the wavefunctions below are assumed to be normalised. We use the same coefficients even after altering the expression, except when it is misleading.

The molecular orbital ${\varphi }_{\mathrm{i}}$ (i = 1, 2) was substituted into Equation (1) through a linear combination of atomic orbitals, ${\phi }_{\mathrm{j}}$,

$${\varphi }_{\mathrm{i}}={\displaystyle \sum }_{\mathrm{j}}{c}_{\mathrm{j}}{\phi }_{\mathrm{j}}$$

(3)

where c_j and ${\phi }_{\mathrm{j}}$ denote the coefficient and atomic orbital of atom j, respectively. For example, ${\phi }_{\mathrm{H}}$ designates the 1s orbital of a hydrogen atom. For j = C (carbon), ${\phi }_{\mathrm{j}}$ should describe one of the C 1s, C 2s, C 2p_x, C 2p_y, and C 2p_z orbitals. Note that all the atomic positions, except for the X atoms, $\overrightarrow{r_{\mathrm{j}}}$, and electron densities on all the atoms in an STF molecule are common between ${\varphi }_1$ and ${\varphi }_2$.

The molecular orbitals ${\varphi }_{\mathrm{i}}$ (i = 1, 2) obtained in this way are as follows:

$${\varphi }_1={\displaystyle \sum }_{\mathrm{j}}{a}_{\mathrm{Y}}{\phi }_{\mathrm{Y}}\left(\overrightarrow{r_{\mathrm{J}}}\right)+{\displaystyle \sum }_{\mathrm{k}}{b}_{\mathrm{S}}{\phi }_{\mathrm{S}}\left(\overrightarrow{r_{\mathrm{k}}}\right)+{\displaystyle \sum }_{\mathrm{l}}{c}_{\mathrm{Se}}{\phi }_{\mathrm{Se}}\left(\overrightarrow{r_{\mathrm{l}}}\right)$$

(4)

$${\varphi }_2={\displaystyle \sum }_{\mathrm{j}}{a}_{\mathrm{Y}}{\phi }_{\mathrm{Y}}\left(\overrightarrow{r_{\mathrm{J}}}\right)+{\displaystyle \sum }_{\mathrm{k}}{b}_{\mathrm{Se}}{\phi }_{\mathrm{Se}}\left(\overrightarrow{r_{\mathrm{k}}}\right)+{\displaystyle \sum }_{\mathrm{l}}{c}_{\mathrm{S}}{\phi }_{\mathrm{S}}\left(\overrightarrow{r_{\mathrm{l}}}\right)$$

(5)

where suffixes j, k, and l indicate the atomic positions in the STF molecule. Suffix Y indicates atoms other than X atoms in the STF molecule. Notably, the first terms in Equations (4) and (5) are identical. Meanwhile, the coordinates $\overrightarrow{r_{\mathrm{k}}}$ and $\overrightarrow{r_{\mathrm{l}}}$ for S and Se atoms indicate the X atoms on the left and right sides of each STF molecule, respectively. Accordingly, the two equations of ${\varphi }_1$ and ${\varphi }_2$ are identical, except for the coordinates $\overrightarrow{r_{\mathrm{k}}}$ and $\overrightarrow{r_{\mathrm{l}}}$ for S and Se atoms (the position vectors in real space for X atoms), respectively. By substituting Equations (4) and (5) in Equation (2), we obtain

$$\begin{array}{ll}{{\varphi }_{\mathrm{STF}}}^{\prime }& =\frac{1}{\sqrt{2}}{\varphi }_1-\frac{1}{\sqrt{2}}{\varphi }_2\\ & =\frac{1}{\sqrt{2}}\left\{{\displaystyle {\displaystyle \sum }_{\mathrm{k}}}{b}_{\mathrm{S}}{\phi }_{\mathrm{S}}\left(\overrightarrow{r_{\mathrm{k}}}\right)+{\displaystyle {\displaystyle \sum }_{\mathrm{l}}}{c}_{\mathrm{Se}}{\phi }_{\mathrm{Se}}\left(\overrightarrow{r_{\mathrm{l}}}\right)\right\}\\ & -\frac{1}{\sqrt{2}}\left\{{\displaystyle {\displaystyle \sum }_{\mathrm{k}}}{b}_{\mathrm{Se}}{\phi }_{\mathrm{Se}}\left(\overrightarrow{r_{\mathrm{k}}}\right)+{\displaystyle {\displaystyle \sum }_{\mathrm{l}}}{c}_{\mathrm{S}}{\phi }_{\mathrm{S}}\left(\overrightarrow{r_{\mathrm{l}}}\right)\right\}\end{array}$$

(6)

which evidently does not describe the electronic configuration required for the solid STF molecules.

By substituting Equations (4) and (5) in Equation (1), we obtain

$$\begin{array}{ll}{\varphi }_{\mathrm{STF}}& =\frac{1}{\sqrt{2}}\left\{{\displaystyle {\displaystyle \sum }_{\mathrm{j}}}{a}_{\mathrm{Y}}{\phi }_{\mathrm{Y}}\left(\overrightarrow{r_{\mathrm{J}}}\right)+{\displaystyle {\displaystyle \sum }_{\mathrm{k}}}{b}_{\mathrm{S}}{\phi }_{\mathrm{S}}\left(\overrightarrow{r_{\mathrm{k}}}\right)+{\displaystyle {\displaystyle \sum }_{\mathrm{l}}}{c}_{\mathrm{Se}}{\phi }_{\mathrm{Se}}\left(\overrightarrow{r_{\mathrm{l}}}\right)\right\}\\ & +\frac{1}{\sqrt{2}}\left\{{\displaystyle {\displaystyle \sum }_{\mathrm{j}}}{a}_{\mathrm{Y}}{\phi }_{\mathrm{Y}}\left(\overrightarrow{r_{\mathrm{J}}}\right)+{\displaystyle {\displaystyle \sum }_{\mathrm{k}}}{b}_{\mathrm{Se}}{\phi }_{\mathrm{Se}}\left(\overrightarrow{r_{\mathrm{k}}}\right)+{\displaystyle {\displaystyle \sum }_{\mathrm{l}}}{c}_{\mathrm{S}}{\phi }_{\mathrm{S}}\left(\overrightarrow{r_{\mathrm{l}}}\right)\right\}\\ & =\frac{1}{\sqrt{2}}\left\{{\displaystyle {\displaystyle \sum }_{\mathrm{j}}}{a}_{\mathrm{Y}}{\phi }_{\mathrm{Y}}\left(\overrightarrow{r_{\mathrm{J}}}\right)+{\displaystyle {\displaystyle \sum }_{\mathrm{k}}}{b}_{\mathrm{S}}{\phi }_{\mathrm{S}}\left(\overrightarrow{r_{\mathrm{k}}}\right)+{\displaystyle {\displaystyle \sum }_{\mathrm{l}}}{c}_{\mathrm{S}}{\phi }_{\mathrm{S}}\left(\overrightarrow{r_{\mathrm{l}}}\right)\right\}\\ & +\frac{1}{\sqrt{2}}\left\{{\displaystyle {\displaystyle \sum }_{\mathrm{j}}}{a}_{\mathrm{Y}}{\phi }_{\mathrm{Y}}\left(\overrightarrow{r_{\mathrm{J}}}\right)+{\displaystyle {\displaystyle \sum }_{\mathrm{k}}}{b}_{\mathrm{Se}}{\phi }_{\mathrm{Se}}\left(\overrightarrow{r_{\mathrm{k}}}\right)+{\displaystyle {\displaystyle \sum }_{\mathrm{l}}}{c}_{\mathrm{Se}}{\phi }_{\mathrm{Se}}\left(\overrightarrow{r_{\mathrm{l}}}\right)\right\}\end{array}$$

(7)

Therefore,

$${\varphi }_{\mathrm{STF}}=\frac{1}{\sqrt{2}}{\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}+\frac{1}{\sqrt{2}}{\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}$$

(8)

where ${\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}$ and ${\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}$ represent the following MOs, respectively.

$${\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}={\displaystyle \sum }_{\mathrm{j}}{a}_{\mathrm{Y}}{\phi }_{\mathrm{Y}}\left(\overrightarrow{r_{\mathrm{J}}}\right)+{\displaystyle \sum }_{\mathrm{k}}{b}_{\mathrm{S}}{\phi }_{\mathrm{S}}\left(\overrightarrow{r_{\mathrm{k}}}\right)+{\displaystyle \sum }_{\mathrm{l}}{c}_{\mathrm{S}}{\phi }_{\mathrm{S}}\left(\overrightarrow{r_{\mathrm{l}}}\right)$$

(9)

$${\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}={\displaystyle \sum }_{\mathrm{j}}{a}_{\mathrm{Y}}{\phi }_{\mathrm{Y}}\left(\overrightarrow{r_{\mathrm{J}}}\right)+{\displaystyle \sum }_{\mathrm{k}}{b}_{\mathrm{Se}}{\phi }_{\mathrm{Se}}\left(\overrightarrow{r_{\mathrm{k}}}\right)+{\displaystyle \sum }_{\mathrm{l}}{c}_{\mathrm{Se}}{\phi }_{\mathrm{Se}}\left(\overrightarrow{r_{\mathrm{l}}}\right)$$

(10)

In short, ${\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}$ and ${\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}$ are molecular orbitals, assuming that all the X atoms in Figure 2 should be either S or Se atoms with identical atomic positions, respectively. As such, the wavefunctions (MOs) of disordered solid-state STF molecules can be described by the equal contributions of the MOs of ${\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}$ and ${\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}$, as shown in Equation (8).

Next, we discuss the electron densities of the STF sites in the solid state using Equation (8). The contribution of MO, ${\varphi }_{\mathrm{STF}}$, to the electron densities of a given STF site in a solid state is described as follows:

$$\begin{array}{ll}2{\displaystyle \int }{\left|{\varphi }_{\mathrm{STF}}\right|}^2\mathrm{d}V& =2{\displaystyle \int }{\left(\frac{1}{\sqrt{2}}{\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}+\frac{1}{\sqrt{2}}{\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}\right)}^{*}\left(\frac{1}{\sqrt{2}}{\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}+\frac{1}{\sqrt{2}}{\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}\right)\mathrm{d}V\\ & =\left({\displaystyle \int }{\left|{\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}\right|}^2\mathrm{d}V+{\displaystyle \int }{\left|{\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}\right|}^2\mathrm{d}V\right)\\ & +\left({\displaystyle \int }{\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}^{*}{\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}^{}\mathrm{d}V+{\displaystyle \int }{\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}^{*}{\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}^{}\mathrm{d}V\right)\end{array}$$

(11)

where the integral is calculated over the entire space. The coefficient “2” is required by spin multiplicity. Thus, by noting that ${\varphi }_{\mathrm{STF}}$, ${\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}$, and ${\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}$ are normalised,

$$2=2+\left({\displaystyle \int }{\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}^{*}{\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}\mathrm{d}V+{\displaystyle \int }{\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}^{*}{\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}\mathrm{d}V\right)$$

Then,

$${\displaystyle \int }{\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}^{*}{\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}\mathrm{d}V+{\displaystyle \int }{\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}^{*}{\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}\mathrm{d}V=0$$

(12)

In other words,

$${\displaystyle \int }{\left|{\varphi }_{\mathrm{STF}}\right|}^2\mathrm{d}V=\frac{1}{2}\left({\displaystyle \int }{\left|{\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}\right|}^2\mathrm{d}V+{\displaystyle \int }{\left|{\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}\right|}^2\mathrm{d}V\right)$$

(13)

Equation (13) indicates that the electron densities of disordered STF molecules can be described by the average of the electron densities of the ${\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}$ and ${\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}$ eigenstates, consistent with the observations.

The discussion thus far indicates that intermolecular interactions between disordered STF molecules in the solid state can be described by replacing the wavefunction of STF, ${\varphi }_{\mathrm{STF}}$, with an average of the wavefunctions of ${\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}$ and ${\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}$ (Equation (8)). In disordered STF CT complexes, one can similarly calculate the overlap (S) and transfer (t) integrals between STF molecules based on Equation (8). The overlap ($S_{\mathrm{STF}}$) and transfer ($t_{\mathrm{STF}}$) integrals between two STF molecules located at $\overrightarrow{R_1}$ and $\overrightarrow{R_2}$, respectively, are:

$$\begin{array}{ll}{S}_{\mathrm{STF}}& ={\displaystyle \int }{\varphi }_{\mathrm{STF}}\ast \left(\overrightarrow{R_1}\right){\varphi }_{\mathrm{STF}}\left(\overrightarrow{R_2}\right)\mathrm{d}V\\ & ={\displaystyle \int }{\left(\frac{1}{\sqrt{2}}{\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}\left(\overrightarrow{R_1}\right)+\frac{1}{\sqrt{2}}{\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}\left(\overrightarrow{R_1}\right)\right)}^{*}(\frac{1}{\sqrt{2}}{\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}\left(\overrightarrow{R_2}\right)\\ & +\frac{1}{\sqrt{2}}{\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}\left(\overrightarrow{R_2}\right))\mathrm{d}V \left(\because \mathrm{Equation}\left(8\right)\right)\\ & \approx \frac{1}{2}\{{\displaystyle \int }{\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}\ast \left(\overrightarrow{R_1}\right){\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}\left(\overrightarrow{R_2}\right)\mathrm{d}V+{\displaystyle \int }{\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}\ast \left(\overrightarrow{R_1}\right){\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}\left(\overrightarrow{R_2}\right)\mathrm{d}V\}\\ & =\frac{1}{2}\left\{S_{\left(\mathrm{X}=\mathrm{S}\right)}+S_{\left(\mathrm{X}=\mathrm{Se}\right)}\right\},\end{array}$$

(14)

where

$$S_{\left(\mathrm{X}=\mathrm{S}\right)}={\displaystyle \int }{\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}\ast \left(\overrightarrow{R_1}\right){\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}\left(\overrightarrow{R_2}\right)\mathrm{d}V$$

(15)

$$S_{\left(\mathrm{X}=\mathrm{Se}\right)}={\displaystyle \int }{\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}\ast \left(\overrightarrow{R_1}\right){\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}\left(\overrightarrow{R_2}\right)\mathrm{d}V$$

(16)

If the cross terms are ignored (Equation (17)):

$${\displaystyle \int }{\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}^{*}\left(\overrightarrow{R_1}\right){\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}^{}\left(\overrightarrow{R_2}\right)\mathrm{d}V+{\displaystyle \int }{\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}^{*}\left(\overrightarrow{R_1}\right){\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}^{}\left(\overrightarrow{R_2}\right)\mathrm{d}V=0$$

(17)

The validity for Equation (17) is discussed after deriving the transfer integrals. Similarly,

$$\begin{array}{ll}{t}_{\mathrm{STF}}& ={\displaystyle \int }{\varphi }_{\mathrm{STF}}\ast \left(\overrightarrow{R_1}\right)\widehat{H}{\varphi }_{\mathrm{STF}}\left(\overrightarrow{R_2}\right)\mathrm{d}V\\ & =\frac{1}{2}\left\{E_{\left(\mathrm{X}=\mathrm{S}\right)}{S}_{\left(\mathrm{X}=\mathrm{S}\right)}+E_{\left(\mathrm{X}=\mathrm{Se}\right)}{S}_{\left(\mathrm{X}=\mathrm{Se}\right)}\right\}\\ & +\frac{1}{2}\{E_{\left(\mathrm{X}=\mathrm{Se}\right)}{\displaystyle \int }{\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}\ast \left(\overrightarrow{R_1}\right){\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}\left(\overrightarrow{R_2}\right)\mathrm{d}V\\ & +E_{\left(\mathrm{X}=\mathrm{S}\right)}{\displaystyle \int }{\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}\ast \left(\overrightarrow{R_1}\right){\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}\left(\overrightarrow{R_2}\right)\mathrm{d}V\} \left(\because \mathrm{Equation}\left(14\right)\right)\end{array}$$

(18)

where E_{(X = S)} and E_{(X = Se)} are the energies of the orbitals of interest, for example, the HOMO. Below, we limit our discussion to the HOMO of an STF molecule, which dominates the physical properties of α-STF₂I₃.

$$\widehat{H}{\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}=E_{\left(\mathrm{X}=\mathrm{S}\right)}{\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}, \widehat{H}{\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}=E_{\left(\mathrm{X}=\mathrm{Se}\right)}{\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)},$$

(19)

Therefore, Equation (18) is further altered to

$$t_{\mathrm{STF}}=\frac{1}{2}\left\{t_{\left(\mathrm{X}=\mathrm{S}\right)}+t_{\left(\mathrm{X}=\mathrm{Se}\right)}\right\},$$

(20)

where

$$t_{\left(\mathrm{X}=\mathrm{S}\right)}=E_{\left(\mathrm{X}=\mathrm{S}\right)}{S}_{\left(\mathrm{X}=\mathrm{S}\right)},t_{\left(\mathrm{X}=\mathrm{Se}\right)}=E_{\left(\mathrm{X}=\mathrm{Se}\right)}{S}_{\left(\mathrm{X}=\mathrm{Se}\right)}$$

(21)

and we approximated

$$E_{\left(\mathrm{X}=\mathrm{Se}\right)}{\displaystyle \int }{\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}\ast \left(\overrightarrow{R_1}\right){\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}\left(\overrightarrow{R_2}\right)\mathrm{d}V+E_{\left(\mathrm{X}=\mathrm{S}\right)}{\displaystyle \int }{\varphi }_{\left(\mathrm{X}=\mathrm{Se}\right)}\ast \left(\overrightarrow{R_1}\right){\varphi }_{\left(\mathrm{X}=\mathrm{S}\right)}\left(\overrightarrow{R_2}\right)\mathrm{d}V\approx 0$$

(22)

$$E_{\left(\mathrm{X}=\mathrm{Se}\right)}\approx E_{\left(\mathrm{X}=\mathrm{S}\right)}=-10 \left(\mathrm{eV}\right)$$

(23)

3. Results and Discussion

The obtained S_STF overlap integrals are listed in

Table 1. Estimated overlap integrals (S_STF × 10³/dimensionless) in α-STF₂I₃. HOMO-HOMO overlaps between adjacent STF radical cations.

Pairs (a)	X = Se	X = S	X = Se or S (50%)
g1	+9.0	−19.7	−5.35
g2	−8.50	−17.9	−13.2
g3	+2.2	−11.7	−4.75
h1	+31.8	−25.9	+2.95
h2	−33.7	−25.3	−29.5
h3	−15.0	−13.3	−14.15
h4	+6.8	−8.6	−0.9

^(a) For the interacting STF pairs in the unit cell, see Figure 4. The extended Hückel parameters utilised in the calculations are reported in [21]. For reference, the values of S_{(X = Se)} and S_{(X = S)} are listed in Table 1.

. Figure 4 shows the interacting pairs of the STF molecules corresponding to the S_STF values in Table 1.

The validities of Equations (17) and (22) were not proven in this general discussion. However, the derived parameters, S_STF (Table 1) and t_STF, semi-quantitatively reproduced the observed electrical properties. In contrast, those based on our former estimation [21] failed to reproduce them even qualitatively. Namely, the anisotropy in conductivity in the ab-plane was observed as σ_b > σ_a [22], which probably agrees with the contour plot of E₁ − E₂ in this study (Figure 5b), exhibiting an oval elongated approximately in the k_a-direction. Note that the ratio between the anisotropic conductivities in the ab-plane is proportional to the ratio between the square of the velocities in the corresponding directions [30]. In our previous method, we did not consider the details of the molecular orbitals shown in Equations (4)–(7). We simply averaged all the terms appearing in Equations (14) and (18) between the four interacting patterns (Figure 3). Consequently, the obtained overlap S_STF and transfer t_STF integrals in our previous study gave the contour plot of E₁ − E₂, exhibiting an oval elongated approximately in the k_b-direction [21], which qualitatively contradict the observed anisotropy in conductivity in the ab-plane. Once the overlap S_STF and transfer t_STF integrals are obtained in this study, we can calculate a tight-binding band structure using the model in our previous study [21]. The band structure obtained (Figure 5) is unique and essential because it indicates that α-STF₂I₃ has a pair of Dirac cones close to merging. The term “merging” implies that the Dirac electrons are close to transforming into normal electrons. Such cases have been extensively investigated theoretically [7,31,32,33,34]. However, to the best of our knowledge, actual examples have never been reported. The Dirac cones could not exhibit such a large anisotropy in resistivity, as observed, except when they are close to merging. Accordingly, the anisotropy observed in the resistivity of α-STF₂I₃ was exclusively explained by this band structure.

4. Conclusions

We proposed a calculation method for intermolecular interactions in disordered systems based on the extended Hückel approximation. Despite the generally negative impression of the effects of disorder on cooperative properties, the physical properties and band structures of α-STF₂I₃ indicated that the fine tuning of intermolecular interactions using disorder in organic CT complexes is an effective new method. This method enables electron densities between S and Se atoms as if there was no disorder in the solid state, which is otherwise impossible. The reconciliation of the established band theories with the disordered systems has remained uninvestigated for a long time. The calculation method proposed in this study was used to investigate the band structures of such systems in detail. It has not only revealed that α-STF₂I₃ is a rare example of a system close to merging, but also that there is a novel electronic state requiring a wavefunction of the type of Equation (1) and the relevant assumptions.