Correlated Band Structure of a Transition Metal Oxide ZnO Obtained from a Many-Body Wave Function Theory

Correlated Band Structure of a Transition Metal Oxide ZnO Obtained from a Many-Body Wave Function Theory

Masayuki Ochi    Ryotaro Arita    Shinji Tsuneyuki Department of Physics, Osaka University, Machikaneyama-cho, Toyonaka, Osaka 560-0043, Japan RIKEN Center for Emergent Matter Science (CEMS), Wako, Saitama 351-0198, Japan Department of Physics, The University of Tokyo, Hongo, Bunkyo-ku, Tokyo 113-0033, Japan Institute for Solid State Physics, The University of Tokyo, Kashiwa, Chiba 277-8581, Japan
July 14, 2019

Obtaining accurate band structures of correlated solids has been one of the most important and challenging problems in first-principles electronic structure calculation. There have been promising recent active developments of wave function theory for condensed matter, but its application to band-structure calculation remains computationally expensive. In this Letter, we report the first application of the bi-orthogonal transcorrelated (BiTC) method: self-consistent, free from adjustable parameters, and systematically improvable many-body wave function theory, to solid-state calculations with electrons: wurtzite ZnO. We find that the BiTC band structure better reproduces the experimental values of the gaps between the bands with different characters than several other conventional methods. This study paves the way for reliable first-principles calculations of the properties of strongly correlated materials.

71.10.-w, 71.15.-m, 71.20.-b

To reveal fertile and nontrivial physics in condensed matter, first-principles electronic-structure calculation has established itself as an indispensable tool in recent studies. For this purpose, density functional theory (DFT) HK (); KS () has played a leading role and has been applied to various materials; however, the limitations of this theoretical framework have come to light. One of the major problems is an inaccurate description of strong electron correlations, e.g., in transition metal oxides. The method GW1 (); GW2 (); GW3 () is a promising way to ameliorate the inaccuracy of the band structures and has been applied to several solids, including -electron systems. However, because the method is often applied without satisfying self-consistency, a nontrivial dependence on the initial DFT calculations is introduced. It has also been reported that the method sometimes exhibits severe difficulty in obtaining converged results GWconvergence (), owing to its perturbative nature. Another possible choice is to construct effective models from DFT that include correlation terms such as Hubbard and to solve the models using elaborated methodologies Models () such as dynamical mean-field theory DMFTreview (); DMFTreview2 (); DMFT2 (). However, the exact correspondence between the effective models and the first-principles Hamiltonian is a nontrivial problem.

Recently, wave function theory (WFT), which had been mainly applied to molecular systems and established itself as the gold standard in theoretical chemistry Szabo (), has become a promising alternative to DFT for accurate descriptions of electron correlation in solids FCIQMC (). Among the most powerful frameworks in WFT are first-principles quantum Monte Carlo (QMC) methods QMCreview (); QMCreview2 (), such as the variational Monte Carlo (VMC), diffusion Monte Carlo (DMC), auxiliary-field quantum Monte Carlo (AFQMC), and full-configuration-interaction (FCI) QMC methods. Other kinds of WFT, called post-Hartree-Fock (post-HF) methods, have also been applied to condensed matter in recent years MP2_Kresse_band (); MP2_Kresse_totE (); MP2_Shepherd (); CCD_Shepherd (); CCD_Shepherd2 (). However, their targets are in most cases limited to solids with small unit cells, owing to their expensive computational cost. In addition, the correlated band structure, which is quite useful in various kinds of theoretical analyses, is not easily obtained in many WFTs. For example, calculation of the band structure in the framework of VMC or DMC requires a large number of single-point calculations of the excited states, which is a clear difference from a mean-field-like approach such as DFT, whereby the whole band structure is obtained at once.

From this viewpoint, the transcorrelated (TC) method BoysHandy (); Handy (); Ten-no (); Umezawa () is a fascinating WFT that can be applied to solids with reliable accuracy and moderate computational cost Sakuma (); TCaccel (); TCjfo (); TCCIS (); TCMP2 (); TCPW (). The TC method adopts the so-called Jastrow ansatz, which is based on a promising strategy often adopted in several WFTs such as QMC methods to describe strong electron correlations; i.e., the electron-electron distance is included into many-body wave functions. Explicitly correlated electronic-structure theory R12 () in quantum chemistry also adopts this strategy. In fact, the Gutzwiller- and Jastrow-correlation factors have often been used to describe strong electron correlation, including the Mott physics in systems such as the Hubbard model VMC (); MottVMC1 (); MottVMC2 (); MottVMC3 (). It is also important to note that, unlike several WFTs, the whole band structure is obtained at once by solving a one-body self-consistent-field (SCF) equation in the TC method, as described later in this Letter. Moreover, the TC method is deterministic, i.e., free from the statistical error, unlike the QMC methods. Accurate calculations for the Hubbard model TCHubbard (); LieAlgebra () and molecular systems LuoTC (); LuoVTC (); CanonicalTC () were also reported using the TC method or other theories that have a close relationship with the TC method. However, insofar as solid-state calculations are concerned, the TC method has so far been applied only to weakly correlated systems.

In this Letter, we present the first application of the TC method to the band-structure calculation of a -electron system: wurtzite ZnO. transition metal oxides have posed theoretical challenges for first-principles band-structure calculations, as it is well known that popular approximations such as the local-density approximation (LDA) fail to provide their band structures accurately, as we shall see later. We find that the TC method with the bi-orthogonal formulation (the BiTC method) BiTC (); TCMP2 () successfully reproduces the experimental band structure of ZnO. We also clarify how the Jastrow factor improves the first-principles description of the correlated electronic states through comparison of the band structures and electron densities among the BiTC and other methods.

The central concept of the TC method is to make use of the similarity transformation of the many-body Hamiltonian with the Jastrow factor ,


where the correlated wave function is represented as and is called the TC Hamiltonian. Here, denotes a pair of space and spin coordinates: . By adopting the so-called Slater-Jastrow ansatz, becomes a Slater determinant consisting of one-electron orbitals, : , and Eq. (1) yields an SCF equation for one-electron orbitals that experience the effective interaction described with the TC Hamiltonian:


where , , and are the external potential including the nucleus-electron interaction note_Vext () and the two- and three-body effective interactions in the TC Hamiltonian, defined as




respectively. As is evident, the HF method can be regarded as the TC method with . Owing to this effective one-body picture, it is possible to treat the many-body correlation with moderate computational cost. In addition, one can obtain the band structure of the quasiparticles by using the real part of the eigenvalues of the matrix on the right-hand side of Eq. (2). One of the authors proved in prior work Umezawa () that such a use of the matrix as quasiparticle energies is consistent with Koopmans’ theorem. We note that one can systematically improve the accuracy of the TC method by utilizing quantum chemical methodologies such as the coupled-cluster and configuration interaction methods to go beyond a single Slater determinant (e.g. Refs. Ten-no (); BiTC (); TCMP2 ()).

Here, the Jastrow function, , is set to the following simple form without adjustable parameters: Sakuma (); QMCreview (); Ceperley (); CeperleyAlder ()


where ( is the number of valence electrons in the simulation cell, is the volume of the simulation cell) and (spin parallel:), (spin antiparallel:). The long-range asymptotic form of this function describes the screening effect of the electron-electron Coulomb interaction BohmPines (). The short-range behavior of the exact Jastrow function should obey the cusp condition cusp (); cusp2 (); cusp_note (). The Jastrow ansatz adopted here works well for state-of-the-art QMC methods QMCreview (); QMCreview2 (). Although our choice of the Jastrow function is rather simple, we shall see that, nevertheless, it works well not only for weakly correlated systems Sakuma (); TCaccel () but also for the 3-electron system. Of course, it is possible to improve a quality of the many-body wave function by using a complicated Jastrow factor, but we adopted this simple trial wave function to realize moderate computational cost.

In this study calc_note (), we adopted the BiTC method BiTC (); TCMP2 (), in which the left one-electron orbitals, in the left Slater determinant X , replace the bra orbitals in the SCF equation (2), while the ket orbitals remain . Because the TC Hamiltonian is non-Hermitian, and become different. We do not show the band structure calculated using the TC method without the bi-orthogonal extension here, because of the large imaginary part of the eigenvalues Sppl ().

Figure 1: Calculated band structures using the LDA (solid lines), all-electron  GWZnO () (green circles), BiTC (solid line), DMC ZnO_DMC () (orange circles), and HF (solid line) methods. Experimental data taken from Ref. ZnOexpt_dots () are shown with black dots in the middle figure. Blue broken lines show the other experimental data for the positions of the conduction-band minimum, the O- bottom ZnOexpt_old (), and the Zn- peak position ZnO3d (), which might correspond to the averaged position of the Zn- bands. The valence-band maximum energy is set to zero.
Band gap [Error] O- bottom [Error] Zn- averaged Zn- bottom
DFT LDA 0.7 [2.7] 5.8
HSE03 2.1 [1.3] 4.9 [0.4] 6.5
(LDA) 2.4 [1.0] 5.2 [0.1] 6.5
(HSE03) 3.2, 3.46 [0.2, 0.06] 6.21 7.2
(GGA+) 2.94 [0.46] 5.6 [0.3] 6.33 7.1
( 6 eV)
+ (GGA+) 3.30 [0.1] 5.5 [0.2] 7.45 8.0
( 6 eV, 1.5 eV)
QS 3.87 [0.47] 5.3 [ 0] 7.2
sc (RPA) 3.8 [0.4] 6.4
sc (-) 3.2 [0.2] 6.7
WFT AFQMC 3.26(16) [0.14]
DMC 3.8(2) [0.4]
BiTC 3.1 [0.3] 5.1 [0.2] 9.1 9.7
HF 11.4 [8.0] 5.7 [0.4] 9.3 9.9
Expt. 3.4 5.3, 5.2(3) 7.5, 7.5(2), 8.5(4),
8.6(2), 8.81(15)
Table 1: Some characteristic values in the band structures of ZnO, as calculated by several methods Tab1_note (). The bottoms of the Zn- and O- bands are evaluated at the point where the valence-band top is set to zero for each method. The errors between the calculated and experimental values shown in Ref. ZnOexpt_old () are presented in parentheses. For LDA, the O- bottom and Zn--averaged levels are not presented here because of the overlap among the O- and Zn- bands. For (HSE03), the Zn- bottom position was read from the density of states presented in Ref. ZnOHSEGW (). All values are in eV. Ref. ZnOHSEGW (). Ref. GWZnO (). Ref. ZnOARPES (). Ref. ZnOQSGW (). Ref. ZnOscGW (). Ref. ZnOAFQMC (). Ref. ZnO_DMC (). Ref. ZnOexpt_old (). Ref. ZnO3d (). Ref. ZnO3d_3 (). Ref. ZnO3d_2 ().

Figure 1 presents the band structures of ZnO calculated with the LDA, all-electron starting from LDA using the LAPW LAPW () method, BiTC, and HF methods. The characteristic energy values in these band structures, as evaluated at the point, are listed in Table 1. Table 1 also lists calculated values with various other methods. By focusing on the gaps between the bands with different characters as listed in Table 1, the BiTC band structure exhibits better accuracy than many other conventional methods, including the method starting from LDA, and accuracy comparable to that for the most accurate varieties of the scheme, i.e., the method starting from HSE03 and the self-consistent methods such as the QS method QSGW (). It is noteworthy that both the and BiTC methods yield such successful results despite being based on conceptually different formulations. It is important that the BiTC method is based on the self-consistent formulation and thus is independent of DFT calculations, whereas the method strongly depends upon the unperturbed DFT calculations, as seen in Table 1. We should also note that, while some methods employ parameters and , which are difficult to determine in the ab initio way, the BiTC method does not use such parameters.

Consistency of the calculated band gaps among WFTs that use similar trial wave functions (AFQMC, DMC, and BiTC) is also remarkable. We again stress that the whole band structure is not easily obtained in QMC simulations and requires many single-point excited-state calculations. We can see that the band structures calculated with the BiTC and HF methods are very similar, except for the band gap. Therefore, the main role of the Jastrow factor used in this study on the band structure seems to be improvement of the size of the band gap through the screening effect of the electron–electron interaction, which is described with the long-ranged asymptotic behavior of the Jastrow factor. To obtain a more accurate band structure, e.g., with respect to the depth of the Zn-3 bands, more elaborated Jastrow factors, as used in QMC studies QMCreview (); QMCreview2 (), will be necessary. Because the Zn- bands are almost flat, a key point might be accurate description of the atomic states, which is an important issue for future investigation.

For the method, it was pointed out that the shallow Zn- bands can strengthen the - hybridization, and thus can result in underestimation of the band gap ZnOARPES (); ZnOscGW (). A similar situation might also be realized in the BiTC band structure, whereas the Zn- bands are rather deep note_GW () and so the band gap can be overestimated. However, the BiTC band gap is also affected by the parameter in the Jastrow factor, as mentioned in the above comparison between the BiTC and HF band structures. Because the parameter used in this study was determined by RPA analysis of the uniform electron gas, it can cause overscreening in the insulator TCjfo (), thereby decreasing the band gap. One possibility is that these two factors are canceled here, but more detailed investigation on other materials is also an important future issue note_GW2 ().

Figure 2: Calculated electron densities on the blue line shown in the crystal structure are presented for the LDA, BiTC, and HF methods with solid black, broken blue, and broken red lines, respectively. The crystal structure was depicted using the VESTA software VESTA ().

Figure 2 presents the electron densities calculated with the LDA, HF, and BiTC methods on the line shown in the crystal structure. The electron density obtained by the BiTC method is defined as Re, where the condition is satisfied due to the bi-orthonormalization condition . We can see that the electron densities of these methods are almost the same. However, a slight increase of the electron density at the atomic sites is observed for the BiTC method compared with the others. Such a tendency is consistent with the fact that the strong divergence of the electron-electron Coulomb repulsion is alleviated for the effective two-body interaction in the similarity-transformed TC Hamiltonian, because the Jastrow factor satisfies the cusp condition. More concretely, the and terms in Eq. (3) yield divergence with a different sign than the electron-electron Coulomb repulsion. As can be seen in the proof of the cusp condition cusp (); cusp2 (), the true many-body wave function should exhibit deformation described with the two-body degrees of freedom near the electron-electron coalescence point, which cannot be represented solely with one-body degrees of freedom. This is a characteristic advantage of the Slater-Jastrow-type wave function for the description of localized electronic states, such as in strongly correlated systems. It is noteworthy that the atomic calculations of the TC method also exhibit a similar tendency for localization Umezawa ().

Finally, we mention the computational effort required for the BiTC calculation. Computation takes place on time scales given by , where , , and are the numbers of -points, occupied bands, and plane waves, respectively note_computation (). This is the same order as that for the HF or hybrid DFT calculations with a prefactor about 20 to 40 note_cmptime (). The BiTC calculation involves neither the frequency index nor the convergence with respect to the number of conduction bands, unlike some perturbative methodologies such as the method. As can be seen from the fact that the hybrid DFT calculations have now been applied to various periodic systems, the computational cost of the BiTC method is reasonable for solid-state calculations. One remaining obstacle for wide application of the BiTC method is that we use the norm-conserving pseudopotential with a very high cutoff energy to handle the semicore states at present. However, this problem is not inherent to the BiTC method and can be overcome in principle by the development of a pseudopotential formalism such as the PAW method PAW1 (); PAW2 () adapted to the TC method, which is an important future issue note_PAW (). We also note that the deformation of the electron density near atomic sites shown in Fig. 2 also implies the importance of careful treatment for the core states. This can also be a common problem for QMC calculations using the Jastrow correlation factor.

To conclude, we apply the bi-orthogonal version of the TC method to wurtzite ZnO and find that it well reproduces the experimental band structure. Our study encourages further investigation of other strongly correlated materials using the BiTC method.

A part of the calculations was performed using supercomputers at the Supercomputer Center, Institute of Solid State Physics, The University of Tokyo. This study was supported by Grant-in-Aid for young scientists (B) (No. 15K17724) from the Japan Society for the Promotion of Science, MEXT Element Strategy Initiative to Form Core Research Center, CREST, JST, and Computational Materials Science Initiative, Japan.


  • (1) P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964).
  • (2) W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965).
  • (3) L. Hedin, Phys. Rev. 139, A796 (1965).
  • (4) M. S. Hybertsen and S. G. Louie, Phys. Rev. Lett. 55, 1418 (1985).
  • (5) M. S. Hybertsen and S. G. Louie, Phys. Rev. B 34, 5390 (1986).
  • (6) B.-C. Shih, Y. Xue, P. Zhang, M. L. Cohen, and S. G. Louie, Phys. Rev. Lett. 105, 146401 (2010).
  • (7) M. Imada and T. Miyake, J. Phys. Soc. Jpn. 79, 112001 (2010).
  • (8) A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, Rev. Mod. Phys. 68, 13 (1996).
  • (9) G. Kotliar, S. Y. Savrasov, K. Haule, V. S. Oudovenko, O. Parcollet, and C. A. Marianetti, Rev. Mod. Phys. 78, 865 (2006).
  • (10) K. Held, Adv. Phys. 56, 829 (2007).
  • (11) A. Szabo and N. S. Ostlund, Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory, (Macmillan, New York, 1982).
  • (12) G. H. Booth, A. Grüneis, G. Kresse, and A. Alavi, Nature (London) 493, 365 (2013).
  • (13) W. M. C. Foulkes, L. Mitas, R. J. Needs, and G. Rajagopal, Rev. Mod. Phys. 73, 33 (2001).
  • (14) J. Kolorenč and L. Mitas, Rep. Prog. Phys. 74, 026502 (2011).
  • (15) M. Marsman, A. Grüneis, J. Paier, G. Kresse, J. Chem. Phys. 130, 184103 (2009).
  • (16) A. Grüneis, M. Marsman, G. Kresse, J. Chem. Phys. 133, 074107 (2010).
  • (17) J. J. Shepherd, A. Grüneis, Phys. Rev. Lett. 110, 226401 (2013).
  • (18) J. J. Shepherd, T. M. Henderson, G. E. Scuseria, Phys. Rev. Lett. 112, 133002 (2014).
  • (19) J. J. Shepherd, T. M. Henderson, G. E. Scuseria, J. Chem. Phys. 140, 124102 (2014).
  • (20) S. F. Boys and N. C. Handy, Proc. R. Soc. London Ser. A 309, 209 (1969); ibid. 310, 43 (1969); ibid. 310, 63 (1969); ibid. 311, 309 (1969).
  • (21) N. C. Handy, Mol. Phys. 21, 817 (1971).
  • (22) S. Ten-no, Chem. Phys. Lett. 330, 169 (2000); ibid. 175 (2000).
  • (23) N. Umezawa and S. Tsuneyuki, J. Chem. Phys. 119, 10015 (2003).
  • (24) R. Sakuma and S. Tsuneyuki, J. Phys. Soc. Jpn. 75, 103705 (2006).
  • (25) M. Ochi, K. Sodeyama, R. Sakuma, and S. Tsuneyuki, J. Chem. Phys. 136, 094108 (2012).
  • (26) M. Ochi, K. Sodeyama, and S. Tsuneyuki, J. Chem. Phys. 140, 074112 (2014).
  • (27) M. Ochi, Y. Yamamoto, R. Arita, and S. Tsuneyuki, J. Chem. Phys. 144, 104109 (2016).
  • (28) M. Ochi and S. Tsuneyuki, J. Chem. Theory Comput. 10, 4098 (2014).
  • (29) M. Ochi and S. Tsuneyuki, Chem. Phys. Lett. 621, 177 (2015).
  • (30) W. Kutzelnigg, Theor. Chim. Acta 68, 445 (1985).
  • (31) H. Yokoyama and H. Shiba, J. Phys. Soc. Jpn. 56, 1490 (1987).
  • (32) M. Capello, F. Becca, M. Fabrizio, S. Sorella, and E. Tosatti, Phys. Rev. Lett. 94, 026406 (2005).
  • (33) H. Yokoyama, T. Miyagawa, and M. Ogata, J. Phys. Soc. Jpn. 80, 084607 (2011).
  • (34) T. Miyagawa and H. Yokoyama, J. Phys. Soc. Jpn. 80, 084705 (2011).
  • (35) S. Tsuneyuki, Prog. Theor. Phys. Suppl. 176, 134 (2008).
  • (36) J. M. Wahlen-Strothman, C. A. Jiménez-Hoyos, T. M. Henderson, and G. E. Scuseria, Phys. Rev. B 91, 041114(R) (2015).
  • (37) H. Luo, W. Hackbusch, and H.-J. Flad, Mol. Phys. 108, 425 (2010).
  • (38) H. Luo, J. Chem. Phys. 133, 154109 (2010).
  • (39) T. Yanai and T. Shiozaki, J. Chem. Phys. 136, 084107 (2012).
  • (40) O. Hino, Y. Tanimura, and S. Ten-no, J. Chem. Phys. 115, 7865 (2001).
  • (41) Note that the external potential does not include the DFT potential unlike some other post-Kohn-Sham treatments.
  • (42) D. M. Ceperley, Phys. Rev. B 18, 3126 (1978).
  • (43) D. M. Ceperley and B. J. Alder, Phys. Rev. Lett. 45, 566 (1980).
  • (44) D. Bohm and D. Pines, Phys. Rev. 92, 609 (1953).
  • (45) T. Kato, Commun. Pure Appl. Math. 10, 151 (1957).
  • (46) R. T. Pack and W. B. Brown, J. Chem. Phys. 45, 556 (1966).
  • (47) To fully satisfy the cusp condition for the singlet and triplet spin pairs, one should adopt the operator representation for the Jastrow function Ten-nocusp (), which introduces the troublesome non-terminating series of the effective interaction in the TC Hamiltonian. Thus, we used this approximated cusp condition in the same manner as in many QMC studies QMCreview (); QMCreview2 ().
  • (48) S. Ten-no, J. Chem. Phys. 121, 117 (2004).
  • (49) The following HF and (Bi)TC calculations were performed using the tc++ code Sakuma (); TCaccel (), which makes use of the LDA PZ81 () orbitals obtained using the tapp code tapp1 (); tapp2 () as an initial guess of the HF and (Bi)TC orbitals. Singularities of the electron-electron Coulomb repulsion term and the Jastrow function in the -space were handled using a method proposed by Gygi and Baldereschi GygiBaldereschi () using an auxiliary function of the same form as Ref. auxfunc (). For calculations of wurtzite ZnO, we used a -mesh, a cutoff energy of 361 Ry, and 150 bands to expand the subspace for block-Davidson diagonalization in the HF and (Bi)TC methods ( with the notation used in Ref. [TCPW, ]). Lattice parameters were extracted from the experiment: Å, Å, and  ZnOstructure (). The same nonlocal norm-conserving pseudopotentials TM (), where the semicore and orbitals of zinc atoms were treated as valence electrons, were used for our LDA, HF, and (Bi)TC calculations. Sppl ()
  • (50) J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981).
  • (51) J. Yamauchi, M. Tsukada, S. Watanabe, and O. Sugino, Phys. Rev. B 54, 5586 (1996).
  • (52) O. Sugino and A. Oshiyama, Phys. Rev. Lett. 68, 1858 (1992).
  • (53) F. Gygi and A. Baldereschi, Phys. Rev. B 34, 4405 (1986).
  • (54) S. Massidda, M. Posternak, and A. Baldereschi, Phys. Rev. B 48, 5058 (1993).
  • (55) H. Karzel, W. Potzel, M. Köfferlein, W. Schiessl, M. Steiner, U. Hiller, G. M. Kalvius, D. W. Mitchell, T. P. Das, P. Blaha, K. Schwarz, and M. P. Pasternak, Phys. Rev. B 53, 11425 (1996).
  • (56) N. Troullier and J. L. Martins, Phys. Rev. B 43, 1993 (1991).
  • (57) See Supplemental Material for more information.
  • (58) O. K. Andersen, Phys. Rev. B 12, 3060 (1975).
  • (59) S. V. Faleev, M. van Schilfgaarde, and T. Kotani, Phys. Rev. Lett. 93, 126406 (2004).
  • (60) M. Usuda, N. Hamada, T. Kotani, and M. van Schilfgaarde, Phys. Rev. B 66, 125101 (2002).
  • (61) G. Zwicker and K. Jacobi, Solid State Commun. 54, 701 (1985).
  • (62) Numerical Data and Functional Relationships in Science and Technology, edited by K. H. Hellwege and O. Madelung, Landolt-Börnstein, New Series, Group III (Springer, Berlin, 1982), Vols. 17a and 22a.
  • (63) L. Y. Lim, S. Lany, Y. J. Chang, E. Rotenberg, A. Zunger, and M. F. Toney, Phys. Rev. B 86, 235113 (2012).
  • (64) A. R. H. Preston, B. J. Ruck, L. F. J. Piper, A. DeMasi, K. E. Smith, A. Schleife, F. Fuchs, F. Bechstedt, J. Chai, and S. M. Durbin, Phys. Rev. B 78, 155114 (2008).
  • (65) T. Kotani, M. van Schilfgaarde, and S. V. Faleev, Phys. Rev. B 76, 165106 (2007).
  • (66) M. Shishkin, M. Marsman, and G. Kresse, Phys. Rev. Lett. 99, 246403 (2007).
  • (67) F. Ma, S. Zhang, and H. Krakauer, New J. Phys. 15, 093017 (2013).
  • (68) J. Yu, L. K. Wagner, and E. Ertekin, J. Chem. Phys. 143, 224707 (2015).
  • (69) L. Ley, R. A. Pollak, F. R. McFeely, S. P. Kowalczyk, and D. A. Shirley, Phys. Rev. B 9, 600 (1974).
  • (70) R. A. Powell, W. E. Spicer, and J. C. McMenamin, Phys. Rev. Lett. 27, 97 (1971).
  • (71) C. J. Vesely, R. L. Hengehold, and D. W. Langer, Phys. Rev. B 5, 2296 (1972).
  • (72) Also for weakly correlated systems, we have observed that localized states (e.g., the 3 states of bulk silicon) tend to become deeper in their band structures TCaccel (); TCjfo (); TCPW ().
  • (73) In our observation, the similar trend is also realized for other semiconductors and insulators with electrons, which include both the ionic and covalent crystals.
  • (74) Refs. GWZnO () and ZnOARPES () used the lattice parameters that are a bit different from ours. However, such difference yields a difference in the band gap and the -band position for less than 0.1 eV in LDA calculations. Therefore, we can conclude that the difference of the lattice parameters among several references has a minor effect on the comparison presented in Table 1.
  • (75) K. Momma and F. Izumi, J. Appl. Crystallogr. 44, 1272 (2011).
  • (76) This scaling is valid only when a single Slater determinant is used (i.e. multi-determinant methodologies are not used) and the Jastrow function is represented as (i.e. not a general two-body function).
  • (77) The ratio of computational time required for the TC and HF calculations is about 15 to 20 TCaccel (). As for the BiTC method, this ratio becomes about 20 to 40 owing to the nonequivalency of the bra and ket orbitals.
  • (78) P. E. Blöchl, Phys. Rev. B 50, 17953 (1994).
  • (79) G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999).
  • (80) Our present implementation requires about 20,000 core hours for solving the BiTC-SCF equation, which is determined by a very large cutoff energy of plane waves. If we can make the cutoff energy an order of magnitude lower, which is a usual demand for other methods, by development of efficient core treatment, this cost will be reduced also by an order of magnitude.

Appendix A Supplemental Material

Appendix B Convergence with respect to the number of -points

We have performed band-structure calculations using a -mesh (32 atoms in the unit cell) for the HF and BiTC methods to confirm the convergence with respect to the number of -points. For the BiTC calculations, the band gap and the relative energy of the Zn--band bottom to the valence-band top are 2.92 and eV, respectively, for a -mesh, whereas they are 3.09 and eV, respectively, for a -mesh (108 atoms in the unit cell). For the HF calculations, the corresponding values are 11.53 and 10.07 eV, respectively, for a -mesh, whereas they are 11.42 and eV, respectively, for a -mesh. The difference between the values calculated using the two -meshes is less than 0.2 eV for all of these values.

In our previous study on the TC method TCaccel (), we saw that the calculated band gap with -points obeys an approximate relation:


where is a constant. Using this approximate relation for the band energies, the residual finite-size error in our calculations using a -mesh is estimated to be less than 0.1 eV. Therefore, we can say that the convergence with respect to the number of -points is sufficiently achieved. Although it is possible that the approximate relation that is used above does not hold for this case ZnO_DMC (), we can expect that the convergence will be sufficient for our discussion because the finite-size error is expected to be significantly reduced in a -mesh from its value for a one.

Appendix C Imaginary part of the TC eigenvalues without the bi-orthogonal formulation

We do not show the band structure calculated using the TC method without the bi-orthogonal extension in the main text. This is because we found that the imaginary part of the eigenvalues of the Zn- bands obtained by solving the one-body SCF equation became as large as 0.1 eV in the TC method, which is more than 1,000 times larger than that in the BiTC method. Because the non-Hermiticity in the (Bi)TC formulation originates from the similarity-transformation, the presumed equivalency between the TC solution and the true many-body eigenstate is suggested to break when the imaginary parts of the eigenvalues become too large. The difference between the band structures obtained with the TC and BiTC methods was small for weakly correlated systems except LiF that is rather localized-electron system TCMP2 ().

Appendix D Electron-density differences among the LDA, HF, and BiTC methods

The differences in the electron densities calculated with the LDA, HF, and BiTC methods are shown in Fig. S1. In Fig. S1, we can verify the tendencies explained in the main text, i.e., that the BiTC electron density tends to be localized near the atomic sites. In addition, we can see an increase of the BiTC electron density along the Zn-O bonds near the oxygen atoms.

Figure 1: The differences between the HF and LDA electron densities and between the BiTC and LDA electron densities are presented in panels (b) and (c), respectively. The electron densities on the plane shown in panel (a) are depicted here.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description