arxiv: v1 [cond-mat.mtrl-sci] 13 Nov 2018

Transcription

1 Typeset with jpsj3.cls <ver..> Letter Analysis of Dirac Point in the Organic onductor α-(bedt-ttf) I 3 Yoshikazu Suzumura arxiv:8.564v [cond-mat.mtrl-sci] 3 Nov 8 Department of Physics, Nagoya University, hikusa-ku, Nagoya , Japan (Received March, 6) The Dirac electron in the organic conductor α-(bedt-ttf) I 3 under pressure is analyzed using a tight-binding model with nearest-neighbor transfer energies and four molecules per unit cell. By noting that the Dirac point between the first and second energy bands emerges or merges followed by the level crossing at a time-reversal invariant momentum (TRIM), an effective Hamiltonian is derived on the basis of these two wave functions at the TRIM, which have different parities associated with an inversion symmetry around the inversion center. We demonstrate that the Dirac point is determined by an intersection of two kinds of lines originating from the Hamiltonian described by symmetric and antisymmetric functions around the TRIM. The present method quantitatively gives a reasonable location of the Dirac point of α-(bedt-ttf) I 3 in a wide pressure range. It is well known that the organic conductor α-(bedt- TTF) I 3 exhibits the Dirac point on the Fermi energy under pressure. ) By utilizing the energy band with such a Dirac cone, several studies have been performed to explore the properties of the Dirac electron. ) However, in contrast to graphene, 3,4) it is not straightforward to obtain the Dirac point on the Brillouin zone owing to the accidental degeneracy at the Dirac point. 5) The existence of such a Dirac point is verified using a product of the parity at the time-reversal invariant momentum (TRIM), 6 9) and the merging or emergence of a pair ofdirac points occurs at the TRIM followed by level crossing. 8) The node of the component of the wave function ) connectsthe Diracpoint with the TRIM resulting in the Berry phase. ) Although the correspondence between the Dirac point and the TRIM has been clarified in these works, it is not yet successful to determine the location of the Dirac point except for the numerical diagonalization or semi-analytical calculation. ) Thus, further study of the role of the TRIM is needed to elucidate the mechanism of the formation of the Dirac point. In this present paper, a method of finding the Dirac point is demonstrated using an effective Hamiltonian based on the parity of the wave function at the TRIM. Figure shows the structure of α-(bedt-ttf) I 3 with four molecules A, A, B, and in the unit cell (the dot-dashed square), which gives a tight-binding model with seven kinds of nearest-neighbor transfer energies shown by a bond, a,,b 4. 3) There is an inversion symmetry around the sites of B,, and the middle of A and A. The tight-binding model of Fig. is written as Ĥ = l β > t α,β;l,l < lα, () l,l α,β where lα > denotes a state vector on the molecular site, and t α,β;l,l is the transfer energy between nearest-neighbor molecular sites. The quantities α and β denote the sites A, A, B, and, and l and l are the position vectors of the cell forming a square lattice with N sites. By using the Fourier transform y a a B b 3 b b b 4 Fig.. (olor online) Structure of α-(bedt-ttf) I 3 on the x-y plane with four molecules A, A, B, and per unit cell (dotdashed square). The bond (dotted line) shows nearest-neighbor transfer energies, a,,b 4. The cross denotes an inversion center located at the middle of A and A sites. 3) lα >= N / kexp[ikl] kα > with a wave vector k = (k x,k y ), Eq. () is rewritten as Ĥ = k Ψ(k) > H S < Ψ(k), where Ψ(k) >= ( / ( ka > + ka > ), / ( ka > ka >), kb >, k >) t with H S being a 4 x 4 matrix Hamiltonian. By using a unitary transformation, H S (k) is rewritten as a real Hamiltonian, ) A A a 3 a x a A a 3 A H(k) = P(k) / H S (k)p(k) / ε ε ε 3 ε 4 = ε ε ε 3 ε 4 ε 3 ε 3 ε 33 ε 34, () ε 4 ε 4 ε 43 ε 44 where ε = a 3 + a cosk y, ε = ε, ε = a sink y, ε 3 = (b + b 3 )cos(k x /), ε 4 = b cos( kx+ky ) + b 4 cos( kx ky ), ε 3 = (b 3 b )sin(k x /),ε 4 = b sin( kx+ky )+ b 4 sin( kx ky ), and ε 34 = a cos(k y /). ε 33 = h 44 = and ε ij = ε ji. b b 3 b 4 b B

2 J. Phys. Soc. Jpn. Letter Author Name P(k) denotes a matrix for the transformation of the base by a π-rotation around the inversion center, which is taken at the middle between A and A. The elements of P(k) are given by [P(k)] γ,γ = p γ (k)δ γ,γ with p γ (k) =,, e ikx, and e i(kx+ky) for γ =,, 3, and 4, respectively. The equation for P(k) is written as P(k) u γ >= p γ (k) u γ >, (3) which gives p γ (k) =,, e ikx, e i(kx+ky), and u γ >= (,,,) t, (,,,) t, (,,,) t, (,,,) t, respectively. The TRIM is given by k = G/ with G being the reciprocal lattice vector, where G/ is written as (k x /π,k y /π) = (,) (Γ point), (,) (X point), (,) (Y point), and (,) (M point). We mention two factors for deriving the effective Hamiltonian. One is the parity of the wave function at the TRIM, which is obtained by P(G/). 8) Noting p γ (G/) = ± in Eq. (3) and a relation one obtains [P(G/),H(G/)] =, (4) P(G/) j(g/) > = p Ej (G/) j(g/) >, (5) H(G/) j(g/) > = E j (G/) j(g/) >, (6) where p Ej (G/) = ±, and E > E > E 3 > E 4. Since the +( ) sign gives the even (odd) parity for both u γ > and j(g/) >, the wave function j(g/) >= γ d jγ(g/) u γ >, with the even (odd) parity, is described by a linear combination of u γ > with even ( odd) parity. The other is the relation between the parity and the matrix element of Eq. (). When the parity of i(g/) > is different from that of j(g/) >, leads to < i(g/) H(G/) j(g/) > = γ,γ d iγ (G/)ε γ,γ d jγ (G/) =, (7) ε γ,γ (G/) =, (8) for p γ (G/)p γ (G/) =. In general, ε γ,γ (G/) for p γ (G/)p γ (G/) =. Thus, for p γ (G/)p γ (G/) = (+), the matrix element ε γ,γ (k) with k = G/+q is antisymmetric(symmetric) as a function of q owing to E j (G/+q) = E j (G/ q), i.e., the time-reversal symmetry. In the case of the Γ point (G/=), the antisymmetric ( symmetric) function with respect to k is given by ε (k), ε 3 (k), and ε 4 (k) (the others) in Eq. (). When i(g/) > and j(g/) > have different parities, the quantity < i(g/) H(G/ + q) j(g/) > becomes antisymmetric as a function of q, and then < i(g/) H(G/+q) j(g/) >=, (9) gives a line passing through G/ owing to ε γ,γ being real. Now, we examine the Dirac point between E (k) and E (k), which emerges at k = G/ followed by a level crossing of E (G/) and E (G/) with different parities. Taking these > and > as the unperturbed states where E j [ E j (G/)] and j > [ j(g/) >], we calculate E (k) and E (k) at an arbitrary k by the perturbation of V k defined as V k = H k H(G/), () and H k = H(k). By extending our previous works, 4,5) in which the base is taken at the Dirac point or the M point in the charge-ordered state, the effective Hamiltonian up to the second order is written as ( ) < Hk > < Ṽk > H eff (k) = < Ṽk > < H, () k > where H k = H(G/)+Ṽk and Ṽk is given by < i Ṽk j >=< i V k j > + < i V k n > n=3,4 ( ) + < n V k j >. E i E n E j E n () From Eq. (), the energies Ẽ(k) and Ẽ(k) and the gap function k are estimated as Ẽ j (k) = < H k > + < H k > ( ) j k, (3) k = Ẽ(k) Ẽ(k) = f(k) +g(k), (4) where f(k) = < H k >, (5) g(k) = < H k > < H k >. (6) The coupling between two bands of < H k > and < H k >vanishesonalinegivenbyf(k) = [Eq.(9)], while < H k > becomes equal to < H k > on a line of g(k) =, as will be shown later. The Dirac point ±k is obtained from k =, (7) which is equivalent to f(k ) = and g(k ) =. The effective Hamiltonian [Eq. ()] consisting of the antisymmetric f(k) and symmetric g(k) has a common feature with that of a 3 PbO 6) in which the Dirac point is obtained using the a-d x y and Pb-p orbitals with different symmetries. We numerically examine the case of α-(bedt- TTF) I 3. The transfer energy t = a,,b 4 (ev) at the uniaxial pressure P (kbar) is estimated as,7) t(p) = t()(+k t P), (8) where a (),,b 4 () =.8,.48,.,.3,.4,.6,.5, and the corresponding K t is given by.89,.67,.5,,., and.3, respectively. The calculation of the Dirac point in terms of the states of the TRIMs (G/ = Γ, X, Y, and M) is performed as follows. (i) At a given pressure, verify the existence of the Dirac point using the condition p E (Γ)p E (X)p E (Y)p E (M) =. 8) (ii) Vary (increase

3 J. Phys. Soc. Jpn. Letter Author Name 3 (a) (b) E i (Γ).4.4 E even odd E d(γ) d γ d γ E i (Y).5.5 E E even odd d(x) d(y) d γ d γ γ P P γ (d). Fig.. (olor online) Pressure (P) dependence of E j (G/) (j= and ) with the even or odd parity where G/ corresponds to the Γ point (a) and Y point (b). The inset shows components of d jγ (G/) for E (open square) and E (closed square) at P=38 kbar in (a), and at P =4 kbar in (b) where the open and closed circles in (b) denote d jγ (X) (j= and, respectively) at 4 kbar. and decrease) pressure to obtain the TRIM at which the level crossing of E (G/) and E (G/) with different parities occurs. In this case, the Dirac point moves to the respective TRIM owing to the merging of a pair of Dirac points. (iii) alculate the Dirac point by substituting the state of the TRIM into Eq. (7). One of these TRIMs is chosen depending on the location of the Dirac point. The parity of the wave function in the present case is found as[(p E (G/),p E (G/)]=(+, ),(, ),(+, ), and (+, ) at the Γ, X, Y, and M points respectively. Then the present method of Eq. () may be applied to the Γ, Y, and M points owing to p E (G/)p E (G/) =. In Figs. (a) and (b), the pressure dependence of E j (G/) (j= and ) is shown for the Γ and Y points wherethe levelcrossingofe (k) ande (k) occursat the Γ point for P = P (= 39. kbar), and at the Y point for P = -.8 kbar(from the extrapolation), respectively. The Dirac point for P P is examined using the wave function of the Γ point, while the wave function of the Y point is applied to the Dirac point being far away from the Γ point. The inset denotes the corresponding d jγ (Γ) at P = 38 kbar (a) and d jγ (Y) at P = 4 kbar (b) [and d jγ (X) for comparison]. Each component has either the even parity or the odd parity. For the Γ point, the wave functions of E, E 3, and E 4 (E ) have the even (odd) parity. For the Y point, the wave functions of E and E 3, (E and E 4 ) have the even (odd) parity. First,weshowthemergingbehaviorbytakingEq.() up to the first order for simplicity. Figure 3(a) shows f(k) at P=38 kbar where f(k) = (the bright line) gives the merging direction, k y.7k x. Figure 3(b) shows the corresponding g(k) where the bright line denotes g(k) =, and g() > at the Γ point. For P > P, the bright line moves from the horizontal axis to the verticalaxisowingtog() <.Notethatg(G/) < forthe Xpoint andg(g/) > for the YandMpoints. The gap function k [Eq. (4)] is obtained in Fig. 3(c) representing the contour just before the merging, which is slightly larger than that of the numerical diagonalization. The contour of k close to the Dirac point is elliptic and the k y /π..... Fig. 3. (olor online) ontours of f(k) (a) and g(k) (b) at P = 38 kbar where the center k = denotes the Γ point. The bright line shows f(k) = and g(k) =. With respect to the Γ point, f(k) is antisymmetric and g(k) is symmetric. The corresponding gap function of k (c) shows a pair of Dirac points (in the bright region) just before the merging. In (d), the solid(dot-dashed) line denotes g(k) = at P= 38 kar (4 kbar) and the dashed line denotes f(k) = for P P. The dotted line denotes g(k) = at P = P (= 39. kbar). principal axis is almost parallel to the merging direction. The merging at G/ is determined by g(g/) =, while f(g/) = is the identity. Figure 3(d) shows the lines of g(k) = for P= 38 kar (solid line) and 4 kbar (dotdashed line), and f(k) = for P P (dashed line). The dotted line denotes g(k) = at P = P, which is given by k y (.±.45)k x. The Dirac point exists only for P < P owing to the presence of an intersection between f(k) = and g(k) =, which is shown by the closed circle with (k x /π,k y /π) ±(.,.8). The absence of the intersection for P < P leads to the formation of the gap in k. Figure 3(d) is also understood by the following analysis with G/ = (the Γ point). For the pressure close to P = P, Eqs. (5) and (6) are expanded with respect to q(= k G/) as g(k) = D(P) x q x + xy q x q y + y q y,(9) and f(k) x q x + y q y, where D(P) D (P P). Figure 3(d) is well described by the parameters estimated as D.39, x.386, xy.4, and y.95, x.4, and y.56. Next, weexamine the Diracpoint by applyingeq. () up to the second-order perturbation. With decreasing pressure, the difference in k between the first- and second-order perturbations increases, and the choice of the Γ point becomes invalid at low pressure, e.g., the obtained Dirac point at P= kbar exhibits a different behavior. Thus, the Dirac point for lower pressures is calculated by taking the Y point where the matrix element with the antisymmetric function of q[= k (,π)] is given by ε, ε 4, ε 3, and ε 34. Figure 4(a) shows f(q) for P= 4 kbar corresponding to the pressure of the Dirac P=4 P=P k x /π P=38

4 4 J. Phys. Soc. Jpn. Letter Author Name k y /π. P= Γ.8.8. k x /π Fig. 4. (olor online) q[= k (,π)] dependence of the contour for f(q) (a), g(q) (b), and q (c) at P=4 kbar, which are obtained from Eqs. (5), (6), and (4) with G/ = (,π) (Y point). The center (q = ) denotes the Y point where f(q) [g(q)] is antisymmetric (symmetric) with respect to q. These (a) and (b) give q in (c) and are compared with q in (d), which is obtained from the numerical diagonalization of H(k). The Dirac point in the bright region is given by (q x/π,q y/π) = ±(.58,.67) for (c) and ±(.6,.68) for (d). electron found in α-(bedt-ttf) I 3. ) The bright line obtained from f(q) = exists also for the second-order perturbation since the intermediate state n > with n = 3 and 4 in Eq. () is orthogonal to either i > or j > with a different parity. Note that such a line exists in the presence and absence of the Dirac point, suggesting an intrinsic property that originates from the states i > and j > with a different parity. Figure 4(b) shows the corresponding g(q) where g(q) = (the bright line) intersects with f(q) = in Fig. 4(a). The existence of the line of g(q) = is understood as follows. At the Y point(center),g(q) > fromthedefinition, andg(q) < at the X point since < H k > becomes smaller than < H k > at the X point, as seen from the comparison of each component of the Y and X points in the inset of Fig. (b). Figure 4(c) shows q [Eq. (4)], which is calculated from Figs. 4(a) and 4(b). Figure 4(d) shows E (q) E (q), which is numerically calculated by the diagonalization of Eq. (). The behavior around the Dirac point shows a good correspondence between Figs. 4(c) and 4(d), suggesting the validity of Eqs. () and () for finding the Dirac point. Figure 5 shows the Dirac point on the plane of k x and k y with some choices of P (kbar). The closed circle is calculated by the numerical diagonalization of Eq. (), where the Dirac point moves from the Γ point (P = P ) to the Y point (P= -.8 kbar) with decreasing pressure. The open square is obtained from Eq. () with the unperturbed state at the Γ point. lose to the Γ point, the difference between the closed circle and the open square is small, but increases with decreasing P. The open circle is obtained from Eq. () with the unperturbed state at the Y point. The difference between the closed circle Fig. 5. Dirac point on the plane of k x and k y with some choices of P (kbar). The closed circle is the exact one and the open square (open circle) is obtained using the unperturbed state of the Γ point (Y point) in Eq. (). and the open circle is small but increases with increasing P from kbar. As shown in Figs. 4(c) and 4(d), the Dirac point and the gap function q at P = 4 kbar are well described by Eq. (). Thus, the Dirac point is attributable to two ingredients given by Figs. 4(a) and 4(b), which exhibit the antisymmetric and symmetric functions with respect to the TRIM (of the Y point) and the resultant lines of f(q) = and g(q) =. Finally, we note a range for using Eq. (3) where the region with > Ẽ(k) E (k) >. covers a wide range including the Y and Γ points, and the Dirac point. However, although the location of the Dirac point is successfully obtained, the present method gives the overtilted Dirac cone for Fig. 4(c) owing to the first term of Eq. (3), and remains to be explored as the next step. In summary, using an effective Hamiltonian based on two states with even and odd parities at the TRIM, we obtained the Dirac point as the intersection of two lines. One comes from the diagonal element with the symmetric function relevant to the crossing of two energy bands, and the other originates from the off-diagonal element with the antisymmetric function leading to the vanishing of the gap. Acknowledgements The author thanks R. Kato for useful discussions. This work was supported by a Grant-in-Aid for Scientific Research (A) (No. 5H8) and () (No ) from the Ministry of Education, ulture, Sports, Science, and Technology, Japan, ) S. Katayama, A. Kobayashi, and Y. Suzumura, J. Phys. Soc. Jpn. 75, 5475 (6). ) K. Kajita, Y. Nishio, N. Tajima, Y. Suzumura, and A. Kobayashi, J. Phys. Soc. Jpn. 83, 7 (4). 3) K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and A. A. Firsov, Nature 438, 97 (5). 4) For example, see the review by T. Ando, J. Phys. Soc. Jpn 74, 777 (5).

5 J. Phys. Soc. Jpn. Letter Author Name 5 5). Herring, Phys. Rev. 5, 365 (937). 6) L. Fu and. L. Kane, Phys. Rev. B 76, 453 (7). 7) T. Mori, J. Phys. Soc. Jpn. 8, 347 (3). 8) F. Piéchon and Y. Suzumura, J. Phys. Soc. Jpn. 8, 3373 (3). 9) T. Kariyado and Y. Hatsugai, Phys. Rev. B 88, 456 (3). ) A. Kobayashi and Y. Suzumura, J. Phys. Soc. Jpn. 8, 5475 (3). ) F. Piéchon and Y. Suzumura, J. Phys. Soc. Jpn. 8, 373 (3). ) Y. Suzumura, T. Morinari, and F. Piéchon, J. Phys. Soc. Jpn. 8, 378 (3). 3) T.Mori, A.Kobayashi, T. Sasaki, H.Kobayashi, G. Saito, and H. Inokuchi, hem. Lett. 3, 957 (984). 4) S. Katayama, A. Kobayashi, and Y. Suzumura, Eur. Phys. J. B. 67, 39 (9). 5) A. Kobayashi, Y. Suzumura, F. Piéchon, and G. Montambaux, Phys. Rev. B 84, 7545 (). 6) T. Kariyado and M. Ogata, J. Phys. Soc. Jpn. 8, 8374 (); T.Kariyado and M.Ogata, J.Phys.Soc.Jpn.8, 647 (). 7) R. Kondo, S. Kagoshima, and J. Harada, Rev. Sci. Instrum. 76, 939 (5).