Next Article in Journal
Microstructure, Hardness and EIS Evaluation of Ti-15Zr-5Nb Dental Alloy
Previous Article in Journal
Thermal Expansion of Alkaline-Earth Borates
Previous Article in Special Issue
Amine-Terminated Modified Succinic Acid-Magnetite Nanoparticles for Effective Removal of Malachite Green Dye from Aqueous Environment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Temperature on Seebeck Coefficient of Nodal Line Semimetal in Molecular Conductor

by
Yoshikazu Suzumura
Department of Physics, Nagoya University, Nagoya 464-8602, Japan
Crystals 2024, 14(7), 601; https://doi.org/10.3390/cryst14070601
Submission received: 25 May 2024 / Revised: 22 June 2024 / Accepted: 26 June 2024 / Published: 28 June 2024
(This article belongs to the Special Issue Nanomaterials: Nanostructures and Magnetic Semiconductors)

Abstract

:
We examine the impact of temperature (T) on the Seebeck coefficient S, i.e., the T dependence of S for a single-component molecular conductor [Pd(dddt)2] (dddt = 5,6-dihydro-1,4-dithiin-2,3-dithiolate) with a half-filled band, where the coefficient is obtained from a ratio of the thermal conductivity to the electrical conductivity. The present paper demonstrates theoretically the novel result of large anisotropy in the Seebeck coefficient components of three-dimensional Dirac electrons in a molecular conductor. The conductor exhibits a nodal line with the energy variation around the chemical potential and provides the density of states (DOS) with a minimum. Using a threedimensional tight-binding (TB) model in the presence of both impurity and electron–phonon (e–p) scatterings, we study the Seebeck coefficient Sy for the molecular stacking and the most conducting direction. The impact of T on Sy exhibits a sign change, where Sy > 0 with a maximum at high temperatures and Sy < 0 with a minimum at low temperatures. The T dependence of Sy suggests that the contribution from the conduction (valence) band is dominant at low (high) temperatures. Further, it is shown that the the Seebeck coefficient components for perpendicular directions Sx and Sz are much smaller than Sy and present no sign change, in contrast to Sy. These results are analyzed in terms of the spectral conductivity as a function of the energy ϵ close to the chemical potential μ.

1. Introduction

Massless Dirac fermions exhibit characteristic properties of electrons, which originate from the energy band with a linear dispersion [1]. The Dirac electron in a bulk system has been found in the following two kinds of molecular conductors [2,3].
One is the organic conductor α -(BEDT-TTF)2I3 under pressures [(BEDT-TTF=bis (ethylenedithio)tetrathiafulvalene], in which the Dirac cone is discovered [4,5] using the tight-binding (TB) model with the transfer energy estimated by the extended Hückel method [6,7,8]. The conductor exhibits a zero-gap state (ZGS) with the density of states (DOS) vanishing linearly at the Fermi energy. The Dirac cone was confirmed by the first-principles DFT calculation [9]. Further, a two-band model [10,11] was proposed to examine the Dirac electron in an organic conductor. Several properties of the Dirac cone appear in the temperature (T) dependence of physical quantities. The reversal of the sign of the Hall coefficient, which occurs for the chemical potential being equal to the energy of the Dirac point, was proposed theoretically [12] and was also observed experimentally in the Hall conductivity [13]. The conductivity of Dirac electrons has been examined using a two-band model with the conduction and valence bands. The static conductivity at absolute zero temperature remains finite with a universal value, i.e., independent of the magnitude of the impurity scattering due to a quantum effect [14]. At finite temperatures, the conductivity depends on the magnitude of the impurity scattering, Γ 0 , which is proportional to the inverse of the life-time by the disorder. With increasing temperature, the conductivity increases for Γ 0 T [15]. Although a monotonic increase in the conductivity is expected, the measured conductivity (or resistivity) on the above organic conductor shows an almost constant behavior at high temperatures [16,17,18,19,20]. Noting the electron–phonon (e–p) interaction enhances the resistivity of the conventional metal at high temperatures; a nearly constant conductivity at high temperatures is explained by adding an acoustic phonon scattering to a two-band model with the Dirac cone [21].
Another Dirac electron system was found in a single-component molecular conductor, [Pd(dddt)2], (dddt = 5,6-dihydro-1,4-dithiin-2,3-dithiolate) [3], which exhibits a three-dimensional Dirac electron system under a high pressure with almost temperature-independent resistivity [22]. First-principles calculations [23] show a three-dimensional Dirac electron system, consisting of HOMO (Highest Occupied Molecular Orbital) and LUMO (Lowest Unoccupied Molecular Orbital) bands, and a TB model exhibits a loop of Dirac points [24]. This system is one of a type of system called a nodal line semimetal, i.e., a loop of Dirac points [25,26,27]. There are several studies on the effective Hamiltonian, where a general two-band model is introduced [28] and the explicit calculation is performed for the nodal line semimetal [29]. The conductivity was also examined [30] to comprehend the almost temperature-independent conductivity. Recently, a TB model was improved using the crystal structure, which was obtained under high pressure [31]. This band calculation shows the DOS, which depends linearly on a wide region of the energy being compatible with the temperature corresponding to almost constant conductivity. Thus, we recalculated the almost constant resistivity in [Pd(dddt)2] using the newly found TB model and by taking account of both impurity and acoustic phonon scatterings [32]. It was shown that an interplay of two kinds of scattering explains the T dependence of the resistivity obtained by the experiment. In addition to two-dimensional systems with Dirac points, three-dimensional systems with nodal line semimetal are of interest due to anisotropic conductivity and a common feature of the almost constant conductivity.
The Seebeck coefficient, which is proportional to a ratio of the thermal conductivity to the electrical conductivity, is a quantity displaying a competition between the conduction and valence bands in Dirac electrons. It does not depend on details of the e–p interaction and impurity scattering. The general formula has been established in terms of linear response theory [33,34,35]. For the Seebeck coefficient, S ν ( ν = x and y) of the organic conductor α -(BEDT-TTF)2I3, there are experiments on the T dependence under hydrostatic pressures [36,37], where the sign of S x (perpendicular to the molecular stacking axis) is positive except for low temperatures depending on samples. There are theoretical studies for S y at ambient pressure with a correlation [38] and for S x under a uniaxial pressure [39]. Using a TB model [40], which was derived from the DFT calculation with the experimental lattice structure, it has been shown that the Seebeck coefficient under hydrostatic pressures is positive, i.e., S x , S y > 0 , at finite temperatures. However, the Seebeck coefficient of the nodal line semimetal has not yet been clarified.
In the present paper, we examine the Seebeck effect of [Pd(dddt)2] in addition to the conductivity [32], where the T dependence of μ takes a crucial role. In Section 2, the model and formulation for the Seebeck coefficient of [Pd(dddt)2] with 1/2-filled band are given. In Section 3, the electronic states are examined by calculating the energy band, the nodal line of Dirac points, DOS, and the T dependence of the chemical potential. In Section 4, we present the Seebeck coefficient and the electric conductivity for both directions being parallel and perpendicular to the molecular stacking (y) axis. They are analyzed in terms of the spectral conductivity [39]. Section 5 is devoted to summary and discussion.

2. Formulation

2.1. Model

We consider a two-dimensional Dirac electron system given by
H = H TB + H p + H e p + H imp .
Here, H TB denotes a TB model of the single-component molecular conductor consisting of four molecules per unit cell, which is shown below. The second term denotes the harmonic phonon given by H p = q ω q b q b q with ω q = v s | q | and = 1. The third term is the electron–phonon (e–p) interaction with a coupling constant g q [41],
H e p = k , γ q g q c γ ( k + q ) c γ ( k ) ( b q + b q ) ,
with c γ ( k ) = α d α γ a α ( k ) , which is obtained by the diagonalization of H TB , as shown later. The e–p scattering is considered within the same band (i.e., intraband) due to the energy conservation with v v s , where v 0.05 [10] denotes the averaged velocity of the Dirac cone. The last term of Equation (1), H imp , denotes a normal impurity scattering.
Figure 1 displays the crystal structure of [Pd(dddt)2] [3,24], which consists of four molecules (1, 2, 3, and 4) with HOMO and LUMO per unit cell providing eight molecular orbitals. These molecules are located on two kinds of layers with the xy plane, where layer 1 consists of molecules 1 and 3, and layer 2 consists of molecules 2 and 4. The z axis is perpendicular to the xy plane of layers 1 and 2, forming a three-dimensional system.
A revised TB model corresponding to Figure 1 has been recently obtained using the crystal structure observed under pressure [31]. There are several kinds of transfer energies between two molecular orbitals, which are listed in Table 1. The interlayer energies in the z direction are given by a (1 and 2 molecules, and 3 and 4 molecules), and c (1 and 4 molecules, and 2 and 3 molecules). The intralayer energies in the xy plane are given by p (1 and 3 molecules), q (2 and 4 molecules), and b (along the molecular stacking y axis). Further, these energies are classified by three kinds of transfer energies given by HOMO-HOMO (H), LUMO-LUMO (L), and HOMO-LUMO (HL).
The TB model Hamiltonian in Equation (1) is expressed as
H TB = i , j = 1 N α , β t i , j ; α , β | i , α j , β | ,
where t i , j ; α , β are transfer energies between nearest-neighbor sites and < i , α | is a state vector. The spin degree of freedom is discarded. In addition, α , β = H1, H2, ⋯, L3, and L4.
Using a Fourier transform | α ( k )   = N 1 / 2 j exp [ i k r j ] | j , α with a wave vector k = ( k x , k y , k z ) , Equation (3) is rewritten as
H TB = k | Φ ( k ) H ^ ( k ) Φ ( k ) | ,
where Φ ( k ) | = ( H 1 | , H 2 | , H 3 | , H 4 | , L 1 | , L 2 | , L 3 | , L 4 | ) . We take the lattice constant as unity and then 0 < | k x | , | k y | , | k z | < π in the first Brillouin zone. The expression H ^ ( k ) is an 8 × 8 matrix Hamiltonian, where h α , β = H ^ ( k ) α , β ; h α , β ( k ) denotes a Fourier transform of t i , j ; α , β with a complex conjugate relation h α , β ( k ) = h β , α ( k ) ¯ ; k = k x a * + k y b * + k z c * ( k x , k y , k z ) , where the k y corresponds to the molecular staking axis, and the lattice constant is taken as unity. Matrix elements h α , β ( k ) are given in the previous work [32]. These energies in the unit of eV are listed in Table 1, where the gap between the energy of HOMO and that of LUMO is taken as Δ E = 0.696 eV to reproduce the energy band of the first principle calculation [24].
The energy band E j ( k ) and the wave function | Ψ j ( k ) , ( j = 1 , 2 , , 8 ) are calculated from
H ^ ( k ) | Ψ j ( k ) = E j ( k ) | Ψ j ( k ) ,
where E 1 > E 2 > > E 8 and
| Ψ j ( k ) = α d α j ( k ) | α ,
with α = H1, H2, H3, H4, L1, L2, L3, and L4.

2.2. Dirac Points and DOS

Since the electron close to the chemical potential is relevant for the electron-hole excitation, we consider only E 4 ( k ) and E 5 ( k ) , i.e., the valence and conduction bands for the present calculation. Thus, E 4 ( k ) and E 5 ( k ) are replaced by E + ( k ) and E ( k ) for the calculation of the transport, while E ± ( k ) represents not only the Dirac cone but also full dispersion of E 4 ( k ) and E 5 ( k ) in the first Brillouin zone. The present energy bands E ± ( k ) provide a nodal line, i.e., a loop of the Dirac point k D , which is obtained from
E + ( k D ) = E ( k D ) .
The chemical potential μ = μ ( T ) is determined self-consistently from
1 N k j = 1 8 f ( E j ( k ) μ ( T ) ) = d ω D ( ω ) f ( ω μ ) = 4 ,
where f ( ω ) = 1 / ( exp [ ω / T ] + 1 ) , with T being temperature in units of eV, and k B = 1 . Equation (8) is the condition of the half-filled band due to the HOMO and LUMO bands. The expression D ( ω ) denotes a density of states (DOS) per spin and per unit cell, which is given by
D ( ω ) = 1 N k γ = ± δ ( ω E γ ( k ) ) ,
where d ω D ( ω ) = 8 .

2.3. Electric Transport

The conductivity is given by [42]
σ ν ( T ) = d ϵ ( f ( ϵ μ ) ϵ ) × σ ν ( ϵ , T ) .
The quantity σ ν ( ϵ , T ) denotes the spectral conductivity [39] for ν = x , y , and z, which is calculated as ( γ , γ = 4 and 5)
σ ν ( ϵ , T ) = γ , γ σ ν γ γ ( ϵ , T ) ,
σ ν γ γ ( ϵ , T ) = e 2 π N k v γ γ ν ( k ) * v γ γ ν ( k ) × Γ γ ( ϵ E γ ( k ) ) 2 + Γ γ 2 × Γ γ ( ϵ E γ ( k ) ) 2 + Γ γ 2 ,
v γ γ ν ( k ) = α β d α γ ( k ) * h α β k ν d β γ ( k ) .
Here, h = 2 π and e denote Planck’s constant and electric charge, respectively. The spectral conductivity depends on T due to the e–p interaction; Γ γ denotes the damping of the electron of the γ band given by
Γ γ = Γ + Γ ph γ ,
where the first term comes from the impurity scattering and the second term corresponding to the phonon scattering is given by [21,43]
Γ ph γ = C 0 R × T | ξ γ , k | ,
R = λ λ 0 ,
λ = | g q | 2 / ω q , ξ γ , k = E γ ( k ) μ , C 0 = 6.25 λ 0 / ( 2 π v 2 ) , v = 0.05 and λ 0 / 2 π v = 0.1 ; λ 0 corresponds to λ for an organic conductor [44,45], and λ becomes independent of | q | for small | q | . For parameters Γ and R, we take mainly Γ = 0.0003 and R = 2 in the present numerical calculation.
In linear response theory [33,34,35], the electric current density j = ( j x , j y , j z ) is obtained by the electric field E = ( E x , E y , E z ) , and the temperature gradient T , i.e., the ν (= x , y , and z) component of the current density is expressed as
j ν = L 11 ν E ν L 12 ν ν T / T ,
where L 11 ν is the electrical conductivity σ ν [42] as shown in Equation (10), and L 12 ν is the thermal conductivity. From (17), the Seebeck coefficient S ν ( T ) is obtained by
S ν ( T ) = L 12 ν T L 11 ν = γ , γ S ν γ γ ( T ) ,
where L 11 ν and L 12 ν are calculated as
L 11 ν = σ ν ( T ) ,
L 12 ν = 1 e d ϵ ( f ( ϵ μ ) ϵ ) ( ϵ μ ) σ ν ( ϵ , T ) ,
and S ν γ γ ( T ) is given by
S ν γ γ = 1 e d ϵ ( f ( ϵ μ ) ϵ ) × ( ϵ μ ) σ ν γ γ ( ϵ , T ) × 1 T L 11 ν .
Noting that f ( ϵ μ ) / ϵ is the even function of ϵ μ , Equations (19) and (21) are calculated as
σ ν ( ϵ , T ) = σ ν ( μ , T ) + σ ν ( μ , T ) ( ϵ μ ) + 1 2 σ ν ( μ , T ) ( ϵ μ ) 2 + ,
e L 12 ν ( T ) = π 2 3 σ ν ( μ , T ) T 2 7 π 4 90 σ ν ( μ , T ) T 4 + .
At low temperatures, the sign change of S ν ( T ) with decreasing T comes from the first term of Equation (23).

3. Electronic States

3.1. Energy Band

Figure 2a shows energy bands E j (j = 4 and 5) for the fixed k z = 0 , where E 4 (upper band) and E 5 (lower band) correspond to the conduction and valence bands, respectively. It is found that E 4 ( k D ) = E 5 ( k D ) = ϵ D and E 4 ( k ) > ϵ D > E 5 ( k ) . A one-dimensional band is seen along the k y direction. Two bands correspond to LUMO and HOMO, which are convex downward and upward, respectively. When HOMO–LUMO coupling is absent, there is an intersection due to overlapping. When HOMO–LUMO coupling is present, the intersection disappears due to a gap except for a Dirac point k D . LUMO (HOMO) corresponds to E 4 ( E 5 ), when LUMO band is larger than HOMO band. The relation is reversed when the LUMO band is smaller than the HOMO band. Figure 2b shows energy bands E j (j = 4 and 5) on the k z k x plane for the fixed k y = 0 . The Dirac point exists between E 4 and E 5 , which correspond to HOMO or LUMO in similar way to those of Figure 2a. Note that ϵ D of Figure 2a is smaller than that of Figure 2b.

3.2. Nodal Line and DOS

Figure 3a shows Dirac points in 3D momentum space forming a nodal line, where A (B) on the line corresponds to a minimum (maximum) of the energy and exists on a plane of k z = 0 ( k y = 0 ). The width of the energy variation along the nodal line is given by 0.002 . The chemical potential is located between A and B. The energy is symmetric with respect to the Γ point, ( k x , k y , k z ) = ( 0 , 0 , 0 ) and k x = 0 . Figure 3b shows DOS as function of ω μ 0 , where μ 0 denotes the chemical potential at T = 0. It is found that D ( μ 0 ) 0 due to the nodal line, where the energy varies around μ 0 . There is an asymmetry of D ( ω ) with respect to ω = μ 0 , which shows D ( μ 0 ω ˜ ) < D ( μ 0 + ω ˜ ) for ω ˜ = ω μ 0 > 0.002 , and D ( μ 0 ω ˜ ) > D ( μ 0 + ω ˜ ) for ω ˜ < 0.002 [32]. The dashed line denotes D ˜ ( μ 0 + ω ˜ ) , which is defined by D ( μ 0 ω ˜ ) . Noting that D ( ω ) | ω μ 0 | without tilting of the Dirac cone, which is proportional to the inverse of the averaged velocity of the Dirac cone [42], the velocity of the valence band ( ω ˜ < 0 ) is larger than that of the conduction band ( ω ˜ > 0 ) except for the momentum around the the nodal line. The inset denotes the T dependence of μ μ 0 , which increases slightly for μ close to the nodal line and decreases with increasing T for μ being away from the nodal line.

4. Seebeck Coefficients

Since the present salt of [Pd(dddt)2] shows the largest conductivity along the y-axis, we first examine the Seebeck coefficient S y . Next, we examine S ν ( ν = x and z) with the direction perpendicular to the y-direction. Further, the mutual relation between the conductivity and the Seebeck coefficient is clarified, where both quantities are determined by the spectral conductivity.

4.1. Coefficient for the y-Axis Direction

Figure 4 shows the conductivity σ y and the band components σ y γ γ as a function of T for some choices of Γ 0 = 0.0003 and 0.0005 with R = 2 ; E 5 ( E 4 ) corresponds to the valence band (conduction band), which is given by ϵ < μ 0 ( ϵ > μ 0 ). However, such a correspondence is invalid in a small region of | ϵ μ 0 |   < 0.0002 due to the variation of the energy on nodal line. It is found that σ y 44 > σ y 55 at low temperatures ( T < 0.002 ), whereas σ y 44 < σ y 55 at high temperatures ( T > 0.002 ). Thus, the contribution from the valence band is larger (smaller) than that from the conduction band at high (low) temperatures. Such a relation can be understood from DOS in Figure 3b, where the average velocity of the Dirac cone of the conduction band is larger (smaller) than that of the valence band for the electron with 0 < ϵ μ 0 < 0.002 ( ϵ μ 0 > 0.002 ) . The band component normalized by the total σ y shows a small dependence on Γ 0 (impurity scattering), while the total one decreases clearly by the increase in Γ 0 .
In Figure 5, the Seebeck coefficient S y ( Γ 0 = 0.0003 and 0.0005) is examined with some choices of R = 1, 2, and 4. It turns out that there is a sign change of S y around T 0.006 , where S y ( > 0 ) at high temperatures takes a maximum and S y ( < 0 ) at low temperatures takes a minimum. At high temperatures, S y decreases with the increase in R due to the effect of the e–p interaction, which is enhanced by the increase in T, as seen from Equation (15). Note that S y for T 0.008 is almost independent of R. At low temperatures, S y reduces to zero in the limit of T 0 , as seen from Equations (22) and (23) [39]. The Seebeck coefficient S y with Γ 0 = 0.0003 is compared with Γ 0 = 0.0005 . The increase in Γ 0 reduces S y ( > 0 ) , implying that the reduction of L 12 y is larger than that of L 11 y . The temperature corresponding to the sign change ( S y = 0 ) decreases for increasing R, but it remains almost the same for the increase in Γ 0 .
Figure 6 shows that the band components of the Seebeck coefficients S y γ γ and S y for Γ 0 = 0.0003 with R = 2 . S y 44 ( < 0 ) comes from the conduction band with ϵ > μ 0 , which is obtained from the LUMO band outside the nodal line and from the HOMO band inside the nodal line; S y 55 ( > 0 ) comes from the valence band with ϵ < μ 0 , which is given by the HOMO band outside the nodal line and by the LUMO band inside the nodal line. The off-diagonal component S y 45 ( = S y 54 ) is negligibly small compared with the diagonal components S y 44 and S y 55 . The total contribution S y ( = S y 44 + S y 55 ) is shown by the right hand axis, where S y < 0 at low temperatures, S y > 0 at high temperatures, and S y = 0 , i.e., the sign change occurs at T 0.006 . The dotted line denotes S y 44 and is compared with S y 55 , where their intersection gives S y = 0 .

4.2. Coefficients for the x and z-Axes Directions

Figure 7 shows S ν ( T ) and σ ν ( T ) for ν = x and z. The conductivity σ ν with ν = x and z is much smaller than σ y due to the large anisotropy of the velocity v x , v z < < v y , since the transfer energy for the y axis is largest due to the molecular stacking direction. Both S x ( > 0 ) and S z ( > 0 ) show no sign change. The fact that σ z ( T ) > σ x ( T ) and S x > S z at high temperatures suggests that the effect of the conductivity is larger than that of the thermal conductivity. At high temperatures, S x ( T ) is much larger than S z , and S z decreases rapidly. Note that S x takes a minimum at low T and S x 0 at T 0.002 . Such a minimum suggests that S x is close to the sign change, as discussed later.
Figure 8 shows band components of the normalized conductivity σ ν γ γ / σ ν with ν = x (solid line) and z (dashed line) for Γ 0 = 0.0003 and R = 2 . The total conductivity is given by σ ν (= σ ν 44 + σ ν 55 + σ ν 54 + σ ν 45 ). Quantities σ ν 44 / σ ν and σ ν 55 / σ ν ( σ ν 45 / σ ν and σ ν 54 / σ ν ) correspond to intraband (interband) contribution, where the interband contribution is negligibly small. The difference between σ ν 55 (the valence band) and σ ν 44 (conduction band) increases by the increase in T. The difference between σ x γ γ / σ x and σ z γ γ / σ z is small compared with that between σ x and σ z . Thus, the main contribution for the conductivity is given by the valence band at high temperatures.
Figure 9 shows band components of the Seebeck coefficient S ν γ γ ( T ) for ν = x and z, which are mainly determined by the intraband contributions, i.e., S ν 55 and S ν 44 . The contribution from the valence band S ν 55 has a positive sign and that of the conduction band S ν 44 has a negative sign. The total Seebeck coefficient has a positive sign due to S ν 55 > | S ν 44 | . The difference between S x 44 and S z 44 is negligibly small, while S x 55 is larger than S z 55 except for low temperatures. The minimum of S x at T 0.002 in Figure 7 is compatible with the fact that S x 55 S x 44 at low temperatures in Figure 9 (the solid line is compared with the dotted line).

4.3. Spectral Conductivity

The Seebeck coefficient is obtained from the spectral conductivity, which is expanded in terms of ϵ μ . From Equations (22) and (23), the Seebeck coefficient is written as
S ν ( T ) = π 2 T 3 e × σ ν ( μ , T ) + ( 7 π 2 / 30 ) σ ν ( μ , T ) T 2 + σ ν ( μ , T ) + ( π 2 / 6 ) σ ν ( μ , T ) T 2 + .
In the limit of low T, the sign of S y ( T ) is determined by that of σ y ( μ , T ) due to S y ( T ) ( π 2 / 3 e ) T σ y ( μ , T ) / σ y ( μ , T ) . Since the sign change of S ν is determined by zero of the numerator, the temperature for the sign change deviates slightly from that of σ ( μ , T ) = 0 . In Figure 10a,b, spectral conductivities σ ν ( ϵ , T ) are shown as a function of ϵ μ 0 .
Figure 10a presents σ y ( ϵ , T ) , which takes a minimum at ϵ = ϵ min . The variation of σ y ( ϵ , T ) at low temperatures is small for | ϵ μ 0 | < 0.0002 , which corresponds to the energy region of the nodal line. A linear increase in σ y ( ϵ , T ) for | ϵ μ 0 | > 0.0002 can be understood from the increase in DOS (Figure 3b). When T increases, ϵ min ( < 0 ) increases to zero, and σ y ( ϵ , T ) at ϵ = ϵ min increases slowly. However, for ϵ being away from min ϵ min , the opposite behavior of the T dependence of σ y ( ϵ , T ) is found, i.e., σ y ( ϵ , T ) decreases with increasing T. These two behaviors can be understood from Equation (12). The T dependence of σ y ( ϵ , T ) is determined by that of Γ γ , which increases by T as shown in Equations (14) and (15). From Equation (12), the numerator gives the increase in σ y ( ϵ , T ) at ϵ = ϵ min , whereas the denominator gives the decrease in σ y ( ϵ , T ) at ϵ being far away from ϵ min . Such a mechanism also explains the T dependence of the conductivity σ y in Figure 4, which increases at low T but decreases at high T. The vertical lines denote the corresponding μ ( T ) , where μ ( 0 ) is shown by the solid line; μ ( T ) increases for 0 < T < 0.001 and decreases for 0.001 < T , as shown in the inset of Figure 3b. Since σ y ( μ , T ) > 0 for 0 < T < 0.002 and σ y ( μ , T ) < 0 at T = 0.004, T corresponding to σ y ( μ , T ) = 0 is lower than that of the sign change of S y = 0 (Figure 5), i.e., T 0.006 . Such a discrepancy suggests that the second term of Equation (23) becomes relevant with increasing T. Figure 10b presents the ϵ μ 0 dependence of σ ν ( ϵ , T ) for ν = x and z, which takes a minimum at ϵ = ϵ min . Since ϵ min > μ ( T ) is found at any T, σ ν ( μ , T ) < 0 at finite temperatures suggests S ν ( T ) > 0 , i.e., no sign change of the Seebeck coefficient. At T 0.002 , the unequality of | σ x ( μ , T ) | < | σ z ( μ , T ) | suggests S x ( T ) < S z ( T ) , while the effect of the higher order of Equation (23) may give rise to a minimum of S x ( T ) in Figure 7.

5. Summary and Discussion

In summary, we calculated the T dependence on the Seebeck coefficient S ν ( T ) ( ν = y , x , and z) of the molecular conductor [Pd(dddt)2] under a high pressure and examined in terms of the spectral conductivity σ ν ( μ , T ) . The conductor exhibits the largest conductivity along the y direction corresponding to the molecular stacking, where the z axis is perpendicular to a two-dimensional yx plane.
Noticeable behavior is found in the T dependence of S y ( T ) . With decreasing T, S y ( T ) changes the sign from positive to negative, as seen in Figure 5. This implies the crossover of the dominant contribution from the hole of the valence band ( E 5 ) to the electron of the conduction band ( E 4 ) (Figure 6). The Seebeck coefficient S y in Figure 5 is determined by the k dependence of the velocity v γ γ y . The sign change is also understood from the energy dependence of DOS being proportional to the inverse of the averaged velocity (Figure 3b), where the average velocity of the conduction (valence) band is larger than that of the valence (conduction) band close to (away from) ω = μ 0 . This is quantitatively understood from the crossover of the band components of the conductivity (Figure 4), where σ y 44 > σ y 55 at lower temperatures and σ y 55 > σ y 44 at higher temperatures. The result of S y ( T ) < 0 at low temperature is consistent with the spectral conductivity at low temperature (Figure 10a), i.e., σ y ( μ , T ) > 0 at T = 0.001. Note that such a sign change has been also found for α -(BEDT-TTF)2I3 under a uniaxial pressure [39], which is understood from the band components of the conductivity and the spectral conductivity. However, the relevance of the Seebeck coefficient to the asymmetry of the Dirac cone is unclear in the case of the uniaxial pressure.
Finally, we note the T dependence of S ν with ν = x and z. There is a large anisotropy of the conductivity, where σ x and σ z are much smaller than σ y , suggesting that the velocity for x and z are much smaller than that for y. In this case, it is complicated to discuss the Seebeck coefficient in terms of the ω dependence of DOS. In fact, there is no change of the sign at low temperatures, i.e., S x > 0 and S z > 0 (Figure 7), implying that the contribution from the valence band is always larger than that from the conduction band (Figure 8 and Figure 9). No sign change at low temperature is also understood from the spectral function, where σ y ( μ , T ) < 0 at low temperatures (Figure 10b). This result shares a common feature with the Seebeck coefficient of the α -(BEDT-TTF)2I3 under a hydrostatic pressure, where the hole-like behavior at finite temperatures with the zero doping is obtained [40].
Here, we note a possible experiment on the Seebeck coefficent. The nodal line semimetal in the single component molecular conductor [Pd(dddt)2] was discovered under pressures of P 12 GPa [31], which are much higher than P 2 GPa corresponding to the experiment on the Seebeck coefficient in organic conductors [36,37]. Since the ambient or low pressures are needed for the measurement of the thermal conductivity, it is a future problem to find such a new molecular conductor that displays the nodal line semimetal with the chemical potentials located between the HOMO and LUMO bands in addition to the small T dependence of resistivity [3].

Funding

This research received no external funding.

Data Availability Statement

Data supporting the findings of this study are available within the paper or from the author.

Acknowledgments

The author thanks M. Ogata for useful discussions on the present work.

Conflicts of Interest

The author declares no conflicts of interest.

References

  1. Novoselov, K.S.; Geim, A.K.; Morozov, S.V.; Jiang, D.; Katsnelson, M.I.; Grigorieva, I.V.; Dubonos, S.V.; Firsov, A.A. Two-dimensional gas of massless Dirac fermions in graphene. Nature 2005, 438, 197–200. [Google Scholar] [CrossRef]
  2. Kajita, K.; Nishio, Y.; Tajima, N.; Suzumura, Y.; Kobayashi, A. Molecular Dirac Fermion Systems—Theoretical and Experimental Approaches. J. Phys. Soc. Jpn. 2014, 83, 072002. [Google Scholar] [CrossRef]
  3. Kato, R.; Cui, H.B.; Tsumuraya, T.; Miyazaki, T.; Suzumura, Y. Emergence of the Dirac Electron System in a Single-Component Molecular Conductor under High Pressure. J. Am. Chem. Soc. 2017, 139, 1770–1773. [Google Scholar] [CrossRef]
  4. Katayama, S.; Kobayashi, A.; Suzumura, Y. Pressure-Induced Zero-Gap Semiconducting State in Organic Conductor α-(BEDT-TTF)2I3. J. Phys. Soc. Jpn. 2006, 75, 054705. [Google Scholar]
  5. Kobayashi, A.; Katayama, S.; Noguchi, K.; Suzumura, Y. Superconductivity in Charge Ordered Organic Conductor –α-(ET)2I3 Salt–. J. Phys. Soc. Jpn. 2004, 73, 3135. [Google Scholar] [CrossRef]
  6. Mori, T.; Kobayashi, A.; Sasaki, Y.; Kobayashi, H.; Saito, G.; Inokuchi, H. Band Structure of two types of (BEDT-TTF)2I3. Chem. Lett. 1984, 13, 957–960. [Google Scholar] [CrossRef]
  7. Kondo, R.; Kagoshima, S.; Harada, J. Crystal structure analysis under uniaxial strain at low temperature using a unique design of four-axis X-ray diffractometer with a fixed sample. Rev. Sci. Instrum. 2005, 76, 093902. [Google Scholar] [CrossRef]
  8. Kondo, R.; Kagoshima, S.; Tajima, N.; Kato, R. Crystal and Electronic Structures of the Quasi-Two-Dimensional Organic Conductor α-(BEDT-TTF)2I3 and Its Selenium Analogue α-(BEDT-TSeF)2I3 under Hydrostatic Pressure at Room Temperature. J. Phys. Soc. Jpn. 2009, 78, 114714. [Google Scholar] [CrossRef]
  9. Kino, H.; Miyazaki, T. First-Principles Study of Electronic Structure in α-(BEDT-TTF)2I3 at Ambient Pressure and with Uniaxial Strain. J. Phys. Soc. Jpn 2006, 75, 034704. [Google Scholar] [CrossRef]
  10. Kobayashi, A.; Katayama, S.; Suzumura, Y.; Fukuyama, H. Massless Fermions in Organic Conductor. J. Phys. Soc. Jpn. 2007, 76, 034711. [Google Scholar] [CrossRef]
  11. Goerbig, M.O.; Fuchs, J.-N.; Montambaux, G.; Piéchon, F. Tilted anisotropic Dirac cones in quinoid-type graphene and α-(BEDT-TTF)2I3. Phys. Rev. B 2008, 78, 045415. [Google Scholar] [CrossRef]
  12. Kobayashi, A.; Suzumura, Y.; Fukuyama, H. Hall Effect and Orbital Diamagnetism in Zerogap State of Molecular Conductor α-(BEDT-TTF)2I3. J. Phys. Soc. Jpn. 2008, 77, 064718. [Google Scholar]
  13. Tajima, N.; Kato, R.; Sugawara, S.; Nishio, Y.; Kajita, K. Interband effects of magnetic field on Hall conductivity in the multilayered massless Dirac fermion system α-(BEDT-TTF)2I3. Phys. Rev. B 2012, 85, 033401. [Google Scholar]
  14. Shon, N.H.; Ando, T. Quantum Transport in Two-Dimensional Graphite System. J. Phys. Soc. Jpn. 1998, 67, 2421. [Google Scholar] [CrossRef]
  15. Peres, N.M.R.; Guinea, F.; Castro Neto, A.H. Electronic properties of disordered two-dimensional carbon. Phys. Rev. B 2006, 83, 125411. [Google Scholar] [CrossRef]
  16. Kajita, K.; Ojiro, T.; Fujii, H.; Nishio, Y.; Kobayashi, H.; Kobayashi, A.; Kato, R. Magnetotransport Phenomena of α-Type (BEDT-TTF)2I3 under High Pressures. J. Phys. Soc. Jpn. 1992, 61, 23. [Google Scholar] [CrossRef]
  17. Tajima, N.; Tamura, M.; Nishio, Y.; Kajita, K.; Iye, Y. Transport Property of an Organic Conductor α-(BEDT-TTF)2I3 under High Pressure –Discovery of a Novel Type of Conductor–. J. Phys. Soc. Jpn. 2000, 69, 543–551. [Google Scholar] [CrossRef]
  18. Tajima, N.; Ebina-Tajima, A.; Tamura, M.; Nishio, Y.; Kajita, K. Effects of Uniaxial Strain on Transport Properties of Organic Conductor α-(BEDT-TTF)2I3 and Discovery of Superconductivity. J. Phys. Soc. Jpn. 2002, 71, 1832. [Google Scholar] [CrossRef]
  19. Tajima, N.; Sugawara, S.; Tamura, M.; Kato, R.; Nishio, Y.; Kajita, K. Transport properties of massless Dirac fermions in an organic conductor α-(BEDT-TTF)2I3 under pressure. EPL 2007, 80, 47002. [Google Scholar] [CrossRef]
  20. Liu, D.; Ishikawa, K.; Takehara, R.; Miyagawa, K.; Tanuma, M.; Kanoda, K. Insulating Nature of Strongly Correlated Massless Dirac Fermions in an Organic Crystal. Phys. Rev. Lett. 2016, 116, 226401. [Google Scholar] [CrossRef]
  21. Suzumura, Y.; Ogata, M. Role of acoustic phonons in exotic conductivity of two-dimensional Dirac electrons. Phys. Rev. B 2018, 98, 161205. [Google Scholar] [CrossRef]
  22. Cui, H.B.; Tsumuraya, T.; Kawasugi, Y.; Kato, R. Pressure Induced Superconductivity and Dirac Cone Formation in Single-Component Molecular Conductors. In Proceedings of the 17th International Conference on High Pressure in Semiconductor Physics (HPSP-17), Tokyo, Japan, 7–11 August 2016. [Google Scholar]
  23. Tsumuraya, T.; Cui, H.B.; Miyazaki, T.; Kato, R. First-Principles Study of Single Component Molecular Crystals under Pressure II. Read at the Meeting of the Physical Society Japan at Hiroshima on March 28, 2013 (28pXN-2). Available online: https://doi.org/10.11316/jpsgaiyo.68.1.4.0_950_2 (accessed on 24 May 2024).
  24. Kato, R.; Suzumura, Y. Novel Dirac Electron in Single-Component Molecular Conductor [Pd(dddt)2] (dddt =5,6-dihydro-1,4-dithiin-2,3-dithiolate). J. Phys. Soc. Jpn. 2017, 86, 064705. [Google Scholar] [CrossRef]
  25. Murakami, S. Phase transition between the quantum spin Hall and insulator phases in 3D: Emergence of a topological gapless phase. New J. Phys. 2007, 9, 356. [Google Scholar] [CrossRef]
  26. Hirayama, M.; Okugawa, R.; Murakami, S. Topological Semimetals Studied by Ab Initio Calculations. J. Phys. Soc. Jpn. 2018, 87, 041002. [Google Scholar] [CrossRef]
  27. Bernevig, A.; Weng, H.; Fang, Z.; Dai, X. Recent Progress in the Study of Topological Semimetals. J. Phys. Soc. Jpn. 2018, 87, 041001. [Google Scholar] [CrossRef]
  28. Liu, Z.; Wang, H.; Wang, Z.F.; Yang, J.; Liu, F. Pressure-induced organic topological nodal-line semimetal in the three-dimensional molecular crystal Pd(dddt)2. Phys. Rev. B 2018, 97, 155138. [Google Scholar]
  29. Tsumuraya, T.; Kato, R.; Suzumura, Y. Effective Hamiltonian of Topological Nodal Line Semimetal in Single-Component Molecular Conductor [Pd(dddt)2] from First-Principles. J. Phys. Soc. Jpn. 2018, 87, 113701. [Google Scholar] [CrossRef]
  30. Suzumura, Y.; Cui, H.B.; Kato, R. Conductivity and Resistivity of Dirac Electrons in Single-Component Molecular Conductor [Pd(dddt)2]. J. Phys. Soc. Jpn. 2018, 87, 084702. [Google Scholar] [CrossRef]
  31. Kato, R.; Cui, H.; Minamidate, T.; Yeung, H.H.-M.; Suzumura, Y. Electronic Structure of a Single-Component Molecular Conductor [Pd(dddt)2] (dddt =5,6-dihydro-1,4-dithiin-2,3-dithiolate) under High Pressure. J. Phys. Soc. Jpn. 2020, 89, 124706. [Google Scholar] [CrossRef]
  32. Suzumura, Y.; Kato, R.; Ogata, M. Electric Transport of Nodal Line Semimetal in Single-Component Molecular Conductors. Crystals 2020, 10, 862. [Google Scholar] [CrossRef]
  33. Kubo, R. Statistical-Mechanical Theory of Irreversible Processes: I. General Theory and Simple Applications to Magnetic and Conduction Problems. J. Phys. Soc. Jpn. 1957, 12, 570. [Google Scholar] [CrossRef]
  34. Luttinger, J.M. Theory of Thermal Transport Coefficients. Phys. Rev. A 1964, 135, 1505. [Google Scholar] [CrossRef]
  35. Ogata, M.; Fukuyama, H. Range of Validity of Sommerfeld–Bethe Relation Associated with Seebeck Coefficient and Phonon Drag Contribution. J. Phys. Soc. Jpn. 2019, 88, 074703. [Google Scholar] [CrossRef]
  36. Kitamura, R.; Tajima, N.; Kajita, K.; Reizo, R.; Tamura, M.; Naito, T.; Nishio, Y. Thermoelectric Power of Multilayered Massless Dirac Fermion System α-(BEDT-TTF)2I3—Charge Ordering and Zero-Gap States–. JPS Conf. Proc. 2014, 1, 012097. [Google Scholar]
  37. Konoike, T.; Sato, M.; Uchida, K.; Osada, T. Anomalous Thermoelectric Transport and Giant Nernst Effect in Multilayered Massless Dirac Fermion System. J. Phys. Soc. Jpn. 2013, 82, 073601. [Google Scholar] [CrossRef]
  38. Ohki, D.; Omori, Y.; Kobayashi, A. Effect of Coulomb interactions on the Seebeck coefficient of the organic Dirac electron system α-(BEDT-TTF)2I3. Phys. Rev. B 2020, 101, 245201. [Google Scholar] [CrossRef]
  39. Suzumura, Y.; Ogata, M. Seebeck Coefficient of Two-dimensional Electrons in Organic Conductor under Pressure. Phys. Rev. 2023, B107, 195416. [Google Scholar] [CrossRef]
  40. Suzumura, Y.; Tsumuraya, T.; Ogata, M. Seebeck Effect of Dirac Electrons in Organic Conductors under Hydrostatic Pressure using a Tight-Binding Model Derived from First Principles. J. Phys. Soc. Jpn. 2024, 93, 054704. [Google Scholar] [CrossRef]
  41. Fröhlich, H. On the theory of superconductivity: The one-dimensional case. Proc. Phys. Soc. A 1954, 223, 296. [Google Scholar]
  42. Katayama, S.; Kobayashi, A.; Suzumura, Y. Electric Conductivity of the Zero-Gap Semiconducting State in α-(BEDT-TTF)2I3 Salt. J. Phys. Soc. Jpn. 2006, 75, 023708. [Google Scholar] [CrossRef]
  43. Abrikosov, A.A.; Gorkov, L.P.; Dzyaloshinskii, I.E. Methods of Quantum Field Theory in Statistical Physics; Prentice-Hall: Englewood Cliffs, NJ, USA, 1963. [Google Scholar]
  44. Rice, M.J.; Pietronero, L.; Brüesh, P. Phase phonons and intramolecular electron-phonon coupling in the organic linear chain semiconductor TEA(TCNQ)2. Solid State Commun. 1977, 21, 757. [Google Scholar] [CrossRef]
  45. Gutfreund, H.; Hartzstein, C.; Weger, M. The electron-phonon coupling strength in TTF-TCNQ. Solid State Commun. 1980, 36, 647. [Google Scholar] [CrossRef]
Figure 1. Crystal structure of [Pd(dddt)2] shown in the xz plane [24], where the molecule is stacked along the y direction perpendicular to the plane. Layer 1 (molecules 1 and 3) and layer 2 (molecules 2 and 4) are parallel to the xy plane and alternated along the z direction. The notations x ,   y , and z correspond to a ,   b , and c of the conventional crystallography.
Figure 1. Crystal structure of [Pd(dddt)2] shown in the xz plane [24], where the molecule is stacked along the y direction perpendicular to the plane. Layer 1 (molecules 1 and 3) and layer 2 (molecules 2 and 4) are parallel to the xy plane and alternated along the z direction. The notations x ,   y , and z correspond to a ,   b , and c of the conventional crystallography.
Crystals 14 00601 g001
Figure 2. (a) Energy band at k z = 0 , where the Dirac point is given by k D = ± ( 0 , 0.075 , 0 ) π ; (b) Energy band at k y = 0 , where the Dirac point is given by k D = ± ( 0.344 , 0 , 0.275 ) π .
Figure 2. (a) Energy band at k z = 0 , where the Dirac point is given by k D = ± ( 0 , 0.075 , 0 ) π ; (b) Energy band at k y = 0 , where the Dirac point is given by k D = ± ( 0.344 , 0 , 0.275 ) π .
Crystals 14 00601 g002
Figure 3. Nodal line (a) and DOS (b) [32]. (a) Closed circle denotes Dirac point in 3D momentum space ( k x , k y , k z ), which gives a nodal line. The points A and B on the line correspond to a Dirac point in Figure 2a and Figure 2b, respectively. They provide a minimum and a maximum of the energy on the nodal line. The chemical potential exists on the line between A and B. (b) DOS [ D ( ω )] is shown as a function of ω μ 0 , where μ 0 denotes the chemical potential at T = 0. The dashed line is drawn to compare with the blue line, where the dashed line and the red line are symmetric around ω = μ 0 . Inset denotes the T dependence of μ μ 0 with μ 0 = 0.5053 .
Figure 3. Nodal line (a) and DOS (b) [32]. (a) Closed circle denotes Dirac point in 3D momentum space ( k x , k y , k z ), which gives a nodal line. The points A and B on the line correspond to a Dirac point in Figure 2a and Figure 2b, respectively. They provide a minimum and a maximum of the energy on the nodal line. The chemical potential exists on the line between A and B. (b) DOS [ D ( ω )] is shown as a function of ω μ 0 , where μ 0 denotes the chemical potential at T = 0. The dashed line is drawn to compare with the blue line, where the dashed line and the red line are symmetric around ω = μ 0 . Inset denotes the T dependence of μ μ 0 with μ 0 = 0.5053 .
Crystals 14 00601 g003
Figure 4. T dependence of the conductivity σ y (right axis) and the band components σ y γ γ (left axis) for Γ 0 = 0.0003 (solid line) and 0.0005 (dotted line) with R = 2 , where σ y = σ y 44 + σ y 55 + σ y 54 + σ y 45 . Quantities σ y 44 / σ y and σ y 55 / σ y ( σ y 45 / σ y and σ y 54 / σ y ) correspond to intraband (interband) contribution.
Figure 4. T dependence of the conductivity σ y (right axis) and the band components σ y γ γ (left axis) for Γ 0 = 0.0003 (solid line) and 0.0005 (dotted line) with R = 2 , where σ y = σ y 44 + σ y 55 + σ y 54 + σ y 45 . Quantities σ y 44 / σ y and σ y 55 / σ y ( σ y 45 / σ y and σ y 54 / σ y ) correspond to intraband (interband) contribution.
Crystals 14 00601 g004
Figure 5. T dependence of the Seebeck coefficient S y with R = 1 (dot-dashed line), 2 (solid line), and 4 (dashed line) for Γ 0 = 0.0003 and 0.0005 (dot).
Figure 5. T dependence of the Seebeck coefficient S y with R = 1 (dot-dashed line), 2 (solid line), and 4 (dashed line) for Γ 0 = 0.0003 and 0.0005 (dot).
Crystals 14 00601 g005
Figure 6. T dependence of the band components of the Seebeck coefficients S y γ γ and S y for Γ 0 = 0.0003 with R = 2 ; S y 44 is negative due to the conduction band with ϵ > μ 0 , whereas S y 55 is positive due to the valence band with ϵ < μ 0 . The interband component S y 45 is negligibly small compared with the intraband components S y 44 and S y 55 . The total contribution S y ( = S y 44 + S y 55 ) is shown by the right hand axis, where S y < 0 at low temperatures and S y > 0 at high temperatures leading to S y = 0 at T 0.006 . The S y 44 (dotted line) is shown to compare with S y 55 .
Figure 6. T dependence of the band components of the Seebeck coefficients S y γ γ and S y for Γ 0 = 0.0003 with R = 2 ; S y 44 is negative due to the conduction band with ϵ > μ 0 , whereas S y 55 is positive due to the valence band with ϵ < μ 0 . The interband component S y 45 is negligibly small compared with the intraband components S y 44 and S y 55 . The total contribution S y ( = S y 44 + S y 55 ) is shown by the right hand axis, where S y < 0 at low temperatures and S y > 0 at high temperatures leading to S y = 0 at T 0.006 . The S y 44 (dotted line) is shown to compare with S y 55 .
Crystals 14 00601 g006
Figure 7. Seebeck coefficient S ν ( T ) and conductivity σ ν ( T ) for ν = x and z with fixed Γ 0 = 0.0003 and R = 2 .
Figure 7. Seebeck coefficient S ν ( T ) and conductivity σ ν ( T ) for ν = x and z with fixed Γ 0 = 0.0003 and R = 2 .
Crystals 14 00601 g007
Figure 8. Components of the normalized conductivity σ ν γ γ / σ ν with ν = x (solid line) and z (dashed line) for Γ 0 = 0.0003 and R = 2 , where σ ν = σ ν 44 + σ ν 55 + σ ν 54 + σ ν 45 . Quantities σ ν 44 / σ ν and σ ν 55 / σ ν ( σ ν 45 / σ ν and σ ν 54 / σ ν ) correspond to intraband (interband) contribution.
Figure 8. Components of the normalized conductivity σ ν γ γ / σ ν with ν = x (solid line) and z (dashed line) for Γ 0 = 0.0003 and R = 2 , where σ ν = σ ν 44 + σ ν 55 + σ ν 54 + σ ν 45 . Quantities σ ν 44 / σ ν and σ ν 55 / σ ν ( σ ν 45 / σ ν and σ ν 54 / σ ν ) correspond to intraband (interband) contribution.
Crystals 14 00601 g008
Figure 9. Components of Seebeck coefficient S ν γ γ ( T ) for ν = x and z as a function of T.
Figure 9. Components of Seebeck coefficient S ν γ γ ( T ) for ν = x and z as a function of T.
Crystals 14 00601 g009
Figure 10. Spectral conductivities (a) σ y ( ϵ , T ) and (b) σ ν ( ϵ , T ) as a function of ϵ μ 0 with Γ 0 = 0.0003 and R = 2 for T = 0.001(1), 0.002 (2), 0.004, and 0.008 (4). The vertical lines denote μ μ 0 for the corresponding T; T = 0 and 0.002.
Figure 10. Spectral conductivities (a) σ y ( ϵ , T ) and (b) σ ν ( ϵ , T ) as a function of ϵ μ 0 with Γ 0 = 0.0003 and R = 2 for T = 0.001(1), 0.002 (2), 0.004, and 0.008 (4). The vertical lines denote μ μ 0 for the corresponding T; T = 0 and 0.002.
Crystals 14 00601 g010
Table 1. Transfer energies for P = 5.9 GPa [31], which are multiplied by 10 3 eV. The energy difference between the HOMO and LUMO is taken as Δ E = 0.696 eV.
Table 1. Transfer energies for P = 5.9 GPa [31], which are multiplied by 10 3 eV. The energy difference between the HOMO and LUMO is taken as Δ E = 0.696 eV.
H H L L H L
b 1 209.3 1.9 51.2 (stacking)
p 1 ( p ) 28.1 12.4 19.9 Layer 1
p 2 17.1
b 2 49.9 80.4 67.2 (stacking)
q 1 ( q ) 10.8 8.1 9.3 Layer 2
q 2 9.2
a 1 28.2 14.6 20.1
a 2 2.2 1.3 1.7 Interlayer
c 1 15.4 12.7 14.1
c 2 3.9 15.8 11.8
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Suzumura, Y. Impact of Temperature on Seebeck Coefficient of Nodal Line Semimetal in Molecular Conductor. Crystals 2024, 14, 601. https://doi.org/10.3390/cryst14070601

AMA Style

Suzumura Y. Impact of Temperature on Seebeck Coefficient of Nodal Line Semimetal in Molecular Conductor. Crystals. 2024; 14(7):601. https://doi.org/10.3390/cryst14070601

Chicago/Turabian Style

Suzumura, Yoshikazu. 2024. "Impact of Temperature on Seebeck Coefficient of Nodal Line Semimetal in Molecular Conductor" Crystals 14, no. 7: 601. https://doi.org/10.3390/cryst14070601

APA Style

Suzumura, Y. (2024). Impact of Temperature on Seebeck Coefficient of Nodal Line Semimetal in Molecular Conductor. Crystals, 14(7), 601. https://doi.org/10.3390/cryst14070601

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop