A practical algorithm for many-electron systems based on the path-integral renormalization group (PIRG) method is proposed in the real-space finite-difference (RSFD) approach. The PIRG method, developed for investigating strongly correlated electron systems, has been successfully applied to some models such as Hubbard models. However, to apply this method to more realistic systems of electrons with long-range Coulomb interactions within the RSFD formalism, the one-body Green's function, which requires large computational resources, is to be replaced with an alternative. For the same reason, an efficient algorithm for computing the Fock matrix is needed. The newly proposed algorithm is free of the one-body Green's function and enables us to compute the Fock matrix efficiently. Our result shows a significant reduction in CPU time and the possibility of using the present algorithm as a practical numerical tool.
KEYWORDS: ab initio calculation; strongly correlated electrons; long-range Coulomb interactions
I TRODUCTIOUp to now, many numerical algorithms for strongly correlated electron systems have been proposed and applied to various systems, e.g., quantum Monte Carlo method [1] , density matrix renormalization group (DMRG) method [2] , configuration interaction method [3] , and coupled cluster method [4] . However, even now, the nature of the ground state still remains a challenge because of the immaturity of numerical tools.The path-integral renormalization group (PIRG) method [5,6] was proposed to obtain the many-body ground state. Unlike the quantum Monte Carlo method, the PIRG method does not suffer from the sign problem [7] . Furthermore, unlike the DMRG method, the PIRG method does not limit the dimensionality of systems because numerical renormalization is carried out in an imaginary time space. In the PIRG method, the ground-state wave function is expressed by a linear combination of basis states, e.g., Slater determinants, in a truncated Hilbert space. While retaining the size of the truncated Hilbert space, the optimized basis states and the ground state are projected out numerically.To make the PIRG method applicable to more realistic systems, we extend the PIRG method with the real-space finite-difference (RSFD) approach in which every physical quantity is defined only on grid points in the discretized space [8][9][10][11] . In this endeavor, The process of "choosing more preferable basis states" becomes the main drawback with respect to the computational cost. In particular, because more basis states tend to be required in the RSFD scheme, one-body Green's functions and the Fock matrix dominate the computational cost and prevent us from applying this method to realistic systems.The aim of this work is to introduce a new algorithm within the framework of the RSFD approach to overcome the above-mentioned problems and to show its applicability to a realistic quantum system.
METHODOLOGY A D UMERICAL APPLICATIO Wave Function RepresentationLet us expand the many-body wave function