Quasi-one-dimensional (Q1D) organic crystals of tetrathiotetracene-iodide, TTT2I3, were synthesized independently and nearly simultaneously     with the aim to find superconductivity in a low dimensional conductor. However, these crystals with rather high electrical conductivity near room temperature, at low temperature showed a transition into a dielectric state. Such transition has firstly observed in the Q1D charge transfer compound TTF-TCNQ (tetrathiafulvalene-tetracyanoquinodimethane)   and was the first experimental confirmation of the structural transition, predicted earlier by Peierls  in 1D conductors.
TTT2I3 is also a charge transfer compound. The orthorhombic crystal structure consists of segregated chains or stacks of plane TTT molecules and of iodine chains. The lattice constants are a = 18.40 Å, b = 4.96 Å and c =18.32 Å, which demonstrates a very pronounced crystal quasi-one-dimensionality. The highly conducting direction is along b. The compound is of mixed valence. Two molecules of TTT give one electron to iodine chain formed of ions that play the role of acceptors. Only TTT chains are conductive and the carriers are holes. The electrons on iodine ions are in a rather localized states and do not participate in the transport. In the crystals grown by sublimation of TTT and iodine in an inert gas flow  the room electrical conductivity σ along b direction achieves (103 - 104) Ω−1∙cm−1, but in those grown from solution   σ ~ (800 - 103) Ω−1∙cm−1. Such variation in σ of crystals, grown in different laboratories, shows that the conductivity properties of TTT stacks are highly sensitive to defects and impurities. It is caused by the purity of initial materials and the conditions of crystal growth. In all crystals, with the lowering of temperature the conductivity firstly grows, reaches a maximum after that falls. The temperature of the maximum, Tmax, and the value of the ratio σmax/σ300 depends on the iodine content. Crystals with a surplus of iodine, TTT2I3.1, have Tmax ~ (34 - 35) K and very sharp fall of σ(T) after the maximum.
The aim of present paper is to demonstrate that this sharp decrease of σ(T) is determined by the Peierls structural transition in the TTT chains. At our knowledge, the Peierls transition in TTT2I3 was not studied neither theoretically, nor experimentally. It is known that the Peierls structural transition is connected by the competition of two processes that take place, when the temperature is decreased. From one side, it is favorable that the lattice distorts because this diminishes the electronic energy of the crystal, lowering the Fermi energy. However, this distortion increases the lattice elastic energy. At some temperature, named the Peierls critical temperature Tp, when the first process prevails over the second one, a Peierls structural transition takes place.
The Peierls structural transition has been studied in many Q1D crystals  -  . In TTF-TCNQ crystals, the transition takes place at 54 K into TCNQ stacks and at 38 K into TTF stacks with the opening of respective band gaps in the electronic spectrum above the Fermi energy and a strong diminution of electrical conductivity. We have also studied the transition at 54 K in a more complete physical model   . With the lowering temperature, some modifications in the phonon spectrum take place   , and at some temperature, the renormalized phonon frequency becomes equal to zero for a given value of the phonon wave vector. At this temperature the Peierls transition occurs.
Note that earlier the crystals of TTT2I3 have been investigated as good candidates for thermoelectric applications   . It was predicted   that after optimization of the carrier concentration in such crystals values of dimensionless thermoelectric figure of merit~4 could be realized. However, not all parameters of these crystals are determined experimentally. Other aim of present paper is to use the investigation of the Peierls transition in order to determine more precisely some parameters of TTT2I3.
We will apply a more complete crystal model   , which takes into account two the most important hole-phonon interactions. The first interaction is similar to that of deformation potential and is determined by the variation of the transfer energy of a carrier from one molecule to the nearest one, caused by acoustic lattice vibrations. Other interaction is of polaron type and is determined by the variation of the polarization energy of molecules surrounding the conduction electron caused by the same acoustic vibrations. The scattering on structural defects is also taken into account. We show that the Peierls structural transition explains the sharp decrease of electrical conductivity in TTT2I3 at low temperature. The dispersion of renormalized phonons and the Peierls critical temperature are determined. For the simplicity, we consider the 2D physical model.
2. The Physical Model
We apply the two-dimensional physical model described in  for TTF-TCNQ crystals, considering the interchain interaction in the plane of TTT stacks small. The Hamiltonian of the crystal in the tight binding and nearest neighbor approximations has the form:
where the first term is the energy operator of free holes in the periodic field of the lattice, are the creation and annihilation operators of such hole with a two-dimensional wave vector k and projections (kx, ky). The energy of the hole, measured from the top of conduction band, has the form:
Here and are the transfer energies of a hole from one molecule to another along the chain (x direction) and in perpendicular direction (y direction). In Equation (1) are creation and annihilation operators of an acoustic phonon with two-dimensional wave vector q and frequency. The second term in the Equation (1) is the energy operator of longitudinal acoustic phonons,
where and are limit frequencies for oscillations in x and y directions. The third term in Equation (1) represents the hole-phonon interactions. As it was mentioned above, two interaction mechanisms are considered: the first is similar to that of deformation potential and the second is of polaron type. The coupling constants of the first interaction are proportional to the derivatives and of and with respect to the intermolecular distances. The coupling constant of second interaction is proportional to the average polarizability of the molecule. This interaction is important for crystals composed of large molecules as TTT, so as is roughly proportional to the volume of molecule.
The square module of matrix element from Equation (1) can be written in the form:
where M is the mass of the TTT molecule, N is the number of molecules in the basic region of the crystal, , the parameters γ1 and γ2 have the sense of the amplitudes ratio of the second hole-phonon interaction to the first one along chains and in the transversal direction:
The analysis shows that the Hamiltonian from the Equation (1) can not explain the sharp decrease of electrical conductivity for temperatures lower than Tmax = 35 K, even, when we vary the crystal parameters and consider only the first interaction mechanism. It is necessary to take into account also the dynamical interaction of carriers with the defects. The static interaction will give contribution to the renormalization of hole spectrum. The defects in TTT2I3 crystals are created due to different coefficients of dilatation of TTT and iodine chains. The Hamiltonian of this interaction is presented in the form:
Here is the matrix element of the hole interaction with a defect, where is the Fourier transformation of
the derivative with respect to intermolecular distance from the energy of interaction of a carrier with a defect, numbers the defects, which are considered linear along x-direction of TTT chains and distributed randomly.
where the constant D = 1.05 and determines the intensity of hole interaction with a defect.
In order to investigate the Peierls transition, the method of temperature dependent retarded Green functions is applied  . The retarded Green functions
for lattice displacements are determined as:
where is abbreviated notation of retarded Green function, indicates an average over a grand canonical ensemble from Equation (1), -are operators in Heisenberg representation, when and when.
The equation of motion for the operator uq is deduced as follows:
On the base of Equation (8), one can obtain the first equation for the Green function:
Further, one can obtain the equation of motion for the new Green function, which will contain new Green functions of higher order of the type etc. We obtain an infinite chain of equa-
tions. In order to cut up the chain, let’s consider that the hole-phonon interaction is weak and express the three-particle Green function through the one-par- ticle Green function as follows
Thus, a closed equation for the function is obtained that corresponds to the random phase approximation.
Now it is conveniently to pass to Fourier transformation of the function after:
As a result, it follows:
From Equations (12) and (13) it results the expression for the Fourier transformation of the lattice displacement Green function
In order to distinguish the retarded Green function it is needed to put, , then the pole of determines the real and im-
aginary part of the renormalized lattice frequency.
The real part of Equation (15) will determine the renormalized lattice frequency Ω(q), as the solution of the transcendent equation
where the principal value of the dimensionless polarization operator takes the form:
Here, and are, respectively, the square module of ma-
trix elements of the hole-phonon interaction from Equation (4), and of hole interaction with defects from Equation (6), the is the Fermi distribution function. The Equation (16) can be solved only numerically.
3. Results and Discussions
Computer simulations are performed for the following parameters  : M = 6.5 ´ 105 me (me is the mass of the free electron), w1 = 0.16 eV, = 0.26 eV∙Å−1, a = 18.35 Å, b = 4.96 Å, c = 18.46 Å. The sound velocity at low temperatures is vs1 = 1.5・105 cm/s along chains (in b direction), d = 0.015, γ1 = 1.7, and γ2 is determined from the relations: γ2 = γ1b5/(a5d). For vs2 in a transversal (in a direction) we have taken 1.35 × 105 cm/s.
Figures 1-4 present the dependences of renormalized phonon frequencies Ω(qx) as functions of qx for different temperatures and different values of qy . In the same graphs, the dependences for initial phonon frequency ω(qx) are presented too. It is seen that the values of Ω(qx) are diminished in comparison with those of frequency ω(qx) in the absence of hole-phonon interaction. This means that the hole-phonon interaction and structural defects diminish the values of lattice elastic constants. Also, one can observe that with a decrease of temperature T the curves change their form, and in dependencies Ω(qx) a minimum appears. This minimum becomes more pronounced at lower temperatures.
Figure 1 shows the case, when qy = 0. In this case the interaction between TTT chains is neglected. The Peierls transition begins at T = 35 K. At this temperature, the electrical conductivity is strongly diminished, so as a gap in the carrier spectrum is fully opened just above the Fermi energy. In addition, it is
Figure 1. Renormalized phonon spectrum Ω(qx) for γ1 = 1.7 and different temperatures. The dashed line is for the spectrum of free phonons. In this case qy = 0.
Figure 2. Renormalized phonon spectrum Ω(qx) for γ1 = 1.7 and different temperatures. The dashed line is for the spectrum of free phonons. The case of qy = π/4.
Figure 3. Renormalized phonon spectrum Ω(qx) for γ1 = 1.7 and different temperatures. The dashed line is for the spectrum of free phonons. The case of qy = 2kF.
Figure 4. Renormalized phonon spectrum Ω(qx) for γ1 = 1.7 and different temperatures. The dashed line is for the spectrum of free phonons. The case of qy = π.
seen that the slope of Ω(qx) at small qx is diminished in comparison with that of ω(qx). This means that the hole-phonon interaction and structural defects have reduced also the sound velocity in a large temperature interval. When the interaction between transversal chains is taken into account (qy ≠ 0), the temperature when Ω(qx) = 0 is diminished.
Figure 2 shows Ω(qx) for qy = π/4 and different temperatures. One can see that Ω(qx) attains zero at T ~ 30 K. Figures 2-4 correspond to the 2D physical model, qy ≠ 0.
When qy = 2kF (Figure 3), the temperature, when Ω(qx) = 0, decreases additionally and has the value of T = 21 K.
Figure 4 shows the dependences of Ω(qx) on qx for qy = π and different temperatures. It is observed that the temperature, when Ω(qx) = 0, decreases still more and equals T = 19 K, thus, our calculations show that at this temperature the Peierls transition is finished. A new superstructure must appear. Unfortunately, at our knowledge, such experiments were not realized. It would be interesting to verify experimentally our conclusions. As it is seen from  , at T = 19 K, the electrical conductivity is strongly reduced, but achieves zero at T = 10 K. This can be explained by further increase of the gap above the Fermi level, when the temperature decreases from 19 K to 10 K. The existence of an energy gap above the Fermi energy at temperatures higher than that of phase transition have been observed in crystals of tetramethyl tetrathiafulvalene  , which have similar to TTT2I3 temperature dependence of electrical conductivity.
Figure 5 and Figure 6 show the dependencies of the real part of dimensionless polarization operator as function of qx for different values of qy and different temperatures at Ω = 0. In Figure 5, it is presented the case, when qy = 0, and the interaction between TTT chains is not taken into account. It is observed a peak near the value of unity. This means that the Peierls transition begins at T = 35 K.
In Figure 6 it is presented the same dependence of polarization operator as function of qx, but for qy = π. From this graph it is observed that, when the interaction between TTT chains is taken into account, the Peierls critical temperature decreases and transition is finished at T = 19 K.
We have investigated the behavior of phonons near Peierls structural transition in quasi-one-dimensional organic crystals of TTT2I3 (tetrathiotetracene iodide)
Figure 5. Polarization operator as function of qx for qy = 0 and T = 35 K.
Figure 6. Polarization operator as function of qx for qy = π and T = 19 K.
type in 2D approximation. A more complete crystal model is applied which takes into account two the most important hole-phonon interactions. One interaction is of deformation potential type and the other is similar to that of polaron. The ratios of amplitudes of second hole-phonon interaction to the first one along chains and in transversal direction are noted by γ1 and γ2, respectively. The interaction of holes with the structural defects in direction of TTT chains is taken into account too. Analytical expression for the polarization operator was obtained in random phase approximation. The method of retarded temperature dependent Green function is applied. The numerical calculations for renormalized phonon spectrum, Ω(qx), for different temperatures are presented in two cases: 1) when qy = 0 and the interaction between transversal chains is neglected and 2) when qy ≠ 0 and interactions between the adjacent chains are considered. It has been established that Peierls transition begins at T ~ 35 K in TTT chains and reduces considerably the electrical conductivity. Due to interchain interaction the transition is finished at T ~ 19 K. It is demonstrated that the hole-phonon interaction and the interactions with the structural defects diminish Ω(qx) and reduce the sound velocity in a large temperature interval.