Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:cond-mat/9209026

Extracting Infinite System Properties from Finite Size Clusters: `Phase Randomization/Boundary Condition Averaging'

J. Tinka Gammel, D. K. Campbell, Jr Loh

Published Sep 19, 1992• 0 citations
arXivHTMLPDF

Abstract

When electron-electron correlations are important, it is often necessary to use exact numerical methods, such as Lanczos diagonalization, to study the full many-body Hamiltonian. Unfortunately, such exact diagonalization methods are restricted to small system sizes. We show that if the Hubbard $U$ term is replaced by a ``periodic Hubbard" term, the full many body Hamiltonian may be numerically exactly solved, even for very large systems (even $>$100 sites), though only for low fillings. However, for half-filled systems and large $U$ this approach is not only no longer exact, it no longer improves extrapolation to larger systems. We discuss how generalized ``randomized variable averaging'' (RVA) or ``phase randomization'' schemes can be reliably employed to improve extrapolation to large system sizes in this regime. This general approach can be combined with any many-body method and is thus of broad interest and applicability.

DOI
10.48550/arxiv.cond-mat/9209026

Related papers

Ground-state energy of the Hubbard model: Cluster-perturbative results

Maciej M. Maśka

When electron correlations are important it is often necessary to use numerical methods to solve the Hamiltonian for a finite system (cluster) ``exactly.'' Unfortunately, such methods are restricted to small systems. We propose to combine the ``exact'' numerical diagonalization for small clusters with the perturbative calculations to take into account the intracluster as well as intercluster interactions.

arXiv:cond-mat/9802171•Physical review. B, Condensed matter•Apr 15, 1998•7 citations

Massively parallel exact diagonalization of strongly correlated systems

Andreas Dolfen

The physics of strongly correlated materials poses one of the most challenging problems in condensed-matter sciences. Standard approximations applicable to wide classes of materials such as the local density approximation fail, due to the importance of the Coulomb repulsion between localized electrons. Instead, we resort to non-perturbative many-body methods. The calculations are, however, only feasible for rather small model systems. The full Hamiltonian of a real material is approximated by a model Hamiltonian comprising only the most important electronic degrees of freedom, while the effect of all other electrons is included in an average way by renormalizing the parameters. Realistic calculations of strongly correlated materials need to include sufficiently many of these electronic degrees of freedom. The new generation of massively parallel supercomputers allows for these realistic calculations. However, exploiting their computational power requires newly devised algorithms. As a solver we use the Lanczos method, which needs the full many-body state of the correlated system. It is thus limited by the available main memory. The foremost problem for a distributed-memory implementation is that the multiplication of the Hamiltonian to the many-body state leads to highly non-local memory access patterns. A simple yet important observation leads to an efficient solution: in the kinetic term of the Hamiltonian the electron-spin is conserved. Thus, writing the many-body vector as a matrix v(iup,idn), where the indices label spin-configurations, we find that the hopping term only connects vector elements that differ in one index. Hence, storing entire slices v(iup,:) on one node, the kinetic term for the spin-down electrons is local to that process. After transposing v, the same is true for the hopping of the spin-up electrons. With an efficient matrix transpose, implemented via MPI collectives, we thus obtain a version of the Lanczos method that efficiently well on JUBL, Jülich's Blue Gene/L. Using the code we treated two-orbital Hubbard Models which show interesting Mott-band insulator transitions. Kohn's criterium which determines the phase of a system by ground-state properties alone, namely the response of the ground-state energies to a vector potential, was employed. This vector potential can be equivalently regarded as a change in boundary conditions, requiring complex wave vectors. We found that half-filled systems become metallic in the region between a band and Mott insulating phase. To extrapolate to infinite system sizes and therefore getting access to an arbitrarily high resolution in reciprocal space, we use cluster perturbation theory (CPT). It works by diagonalizing a finite cluster exactly and then treating the hopping between identical clusters perturbatively in strong coupling perturbation theory. This leads to an effectively infinite chain. Due to our efficient Lanczos implementation it is possible to calculate relatively large clusters, minimizing finite-size effects. The combination of our Lanczos code and CPT allowed for a treatment of the the organic metal TTF-TCNQ. We were able to solve a long-standing problem in the interpretation of its experimental spectra. Being an interesting one-dimensional metal TTF-TCNQ has been studied thoroughly for more than 30 years. Theoretical estimates as well as experimental measurements suggest values for the parameters of the Hubbard model t and U. With these values Hubbard model calculations show qualitatively the same features as APRES experiments, namely the signatures of spin-charge separation at the Gamma-point. But in order to fit the experimental data to the numerical calculations t has to be doubled, which translates into a broadening of the spectra. This ad-hoc adjustment is not only unsatisfactory but also the temperature dependence of the spectral function is incorrectly reproduced. We resolved this problem by properly including nearest-neighbor Coulomb interactions. Recent DFT calculations show that indeed the commonly used values of t and U are correct. But Coulomb interactions between electrons on neighboring sites are also significant and must not be neglected. Thus, we studied the effect of the next-neighbor interaction parameter V being at about U/2. Including V in our Lanczos calculations for obtaining the angular-integrated spectral function shows that its effect is to broaden the spectra, comparable to increasing the value of t. Moreover, we found that doubling t indeed mimics the effect of Vapprox U/2. Surprisingly, the broadening of the spectrum as a result of V can be understood in first-order Rayleigh-Schrödinger perturbation theory.

•Jan 1, 2011•3 citations

Exact Diagonalization Approach for the infinite D Hubbard Model

Michel Caffarel, Werner Krauth

We present a powerful method for calculating the thermodynamic properties of the Hubbard model in infinite dimensions, using an exact diagonalization of an Anderson model with a finite number of sites. At finite temperatures, the explicit diagonalization of the Anderson Hamiltonian allows the calculation of Green's functions with a resolution far superior to that of Quantum Monte Carlo calculations. At zero temperature, the Lanczòs method is used and yields the essentially exact zero-temperature solution of the model, except in a region of very small frequencies. Numerical results for the half-filled case in the paramagnetic phase (quasi-particle weight, self-energy, and also real-frequency spectral densities) are presented.

arXiv:cond-mat/9306057•Jun 28, 1993•0 citations

Magnetic correlations in the Anderson lattice: An exact-diagonalization study

Paolo Maria Santini, Lucio Claudio Andreani, Hans Beck

The ground-state properties of the periodic Anderson model on one-dimensional chains from three to six sites are studied by exact diagonalization. An appropriate average over boundary conditions is found to be necessary in order to eliminate finite-size oscillations. In the weak-coupling limit, the results for the short-range correlations suggest a smooth transition from a ferromagnetic state at low electron concentration to an antiferromagnetic state at half-filling, rather than a sudden jump at a critical concentration.

Physical review. B, Condensed matter•Jan 1, 1993•18 citations

Density-matrix functional theory of the Hubbard model: An exact numerical study

Román López-Sandoval, G. M. Pastor

A density-functional theory for many-body lattice models is considered in which the single-particle density matrix ${\ensuremath{\gamma}}_{\mathrm{ij}}$ is the basic variable. Eigenvalue equations are derived for solving Levy's constrained search of the interaction energy functional $W[{\ensuremath{\gamma}}_{\mathrm{ij}}].$ $W[{\ensuremath{\gamma}}_{\mathrm{ij}}]$ is expressed as the sum of Hartree-Fock energy ${E}_{\mathrm{HF}}[{\ensuremath{\gamma}}_{\mathrm{ij}}]$ and the correlation energy ${E}_{C}[{\ensuremath{\gamma}}_{\mathrm{ij}}].$ Exact results are obtained for ${E}_{C}({\ensuremath{\gamma}}_{12})$ of the Hubbard model on various periodic lattices, where ${\ensuremath{\gamma}}_{\mathrm{ij}}={\ensuremath{\gamma}}_{12}$ for all nearest neighbors i and j. The functional dependence of ${E}_{C}({\ensuremath{\gamma}}_{12})$ is analyzed by varying the number of sites ${N}_{a},$ band filling ${N}_{e},$ and lattice structure. The infinite one-dimensional chain and one-, two-, or three-dimensional finite clusters with periodic boundary conditions are considered. The properties of ${E}_{C}({\ensuremath{\gamma}}_{12})$ are discussed in the limits of weak $({\ensuremath{\gamma}}_{12}\ensuremath{\simeq}{\ensuremath{\gamma}}_{12}^{0})$ and strong $({\ensuremath{\gamma}}_{12}\ensuremath{\simeq}{\ensuremath{\gamma}}_{12}^{\ensuremath{\infty}})$ electronic correlations, and in the crossover region $({\ensuremath{\gamma}}_{12}^{\ensuremath{\infty}}<~{\ensuremath{\gamma}}_{12}<~{\ensuremath{\gamma}}_{12}^{0}).$ Using an appropriate scaling we observe that ${\ensuremath{\varepsilon}}_{C}{(g}_{12}{)=E}_{C}{/E}_{\mathrm{HF}}$ has a pseudo-universal behavior as a function of ${g}_{12}=({\ensuremath{\gamma}}_{12}\ensuremath{-}{\ensuremath{\gamma}}_{12}^{\ensuremath{\infty}})/({\ensuremath{\gamma}}_{12}^{0}\ensuremath{-}{\ensuremath{\gamma}}_{12}^{\ensuremath{\infty}}).$ The fact that ${\ensuremath{\varepsilon}}_{C}{(g}_{12})$ depends weakly on ${N}_{a},$ ${N}_{e},$ and lattice structure suggests that the correlation energy of extended systems could be obtained quite accurately from finite-cluster calculations. Finally, the behaviors of ${E}_{C}({\ensuremath{\gamma}}_{12})$ for repulsive $(U>0)$ and attractive $(U<0)$ interactions are contrasted.

arXiv:cond-mat/9910264•Physical review. B, Condensed matter•Jan 15, 2000•45 citations

Theoretical study of interacting low dimensional electron systems

Jianzhong Zheng

This thesis addresses the question of a strongly interacting many-body problem: the Hubbard model on one- and two-dimensional lattices. The technique employed is the small-cluster method wherein the full many-body Hamiltonian for finite clusters with periodic boundary conditions (PBC's) is solved exactly. We first apply this technique to an extended one-dimensional chain and a two-dimensional square-lattice Hubbard model. We then construct a two-dimensional overlayer/substrate Hubbard model with first- and second-nearest-neighbor hopping, on-site and nearest-neighbor interactions, and overlayer-site orbital energy. The ground-state properties are investigated and discussed in a many-body picture.

•May 29, 2019•0 citations
When electron-electron correlations are important, it is often necessary to use exact numerical methods, such as Lanczos diagonalization, to study the full many-body Hamiltonian. Unfortunately, such exact diagonalization methods are restricted to small system sizes. We show that if the Hubbard UUU term is replaced by a periodic Hubbard" term, the full many body Hamiltonian may be numerically exactly solved, even for very large systems (even $>$100 sites), though only for low fillings. However, for half-filled systems and large $U$ this approach is not only no longer exact, it no longer improves extrapolation to larger systems. We discuss how generalized randomized variable averaging'' (RVA) or ``phase randomization'' schemes can be reliably employed to improve extrapolation to large system sizes in this regime. This general approach can be combined with any many-body method and is thus of broad interest and applicability.
The physics of strongly correlated materials poses one of the most challenging problems in condensed-matter sciences. Standard approximations applicable to wide classes of materials such as the local density approximation fail, due to the importance of the Coulomb repulsion between localized electrons. Instead, we resort to non-perturbative many-body methods. The calculations are, however, only feasible for rather small model systems. The full Hamiltonian of a real material is approximated by a model Hamiltonian comprising only the most important electronic degrees of freedom, while the effect of all other electrons is included in an average way by renormalizing the parameters. Realistic calculations of strongly correlated materials need to include sufficiently many of these electronic degrees of freedom. The new generation of massively parallel supercomputers allows for these realistic calculations. However, exploiting their computational power requires newly devised algorithms. As a solver we use the Lanczos method, which needs the full many-body state of the correlated system. It is thus limited by the available main memory. The foremost problem for a distributed-memory implementation is that the multiplication of the Hamiltonian to the many-body state leads to highly non-local memory access patterns. A simple yet important observation leads to an efficient solution: in the kinetic term of the Hamiltonian the electron-spin is conserved. Thus, writing the many-body vector as a matrix v(iup,idn), where the indices label spin-configurations, we find that the hopping term only connects vector elements that differ in one index. Hence, storing entire slices v(iup,:) on one node, the kinetic term for the spin-down electrons is local to that process. After transposing v, the same is true for the hopping of the spin-up electrons. With an efficient matrix transpose, implemented via MPI collectives, we thus obtain a version of the Lanczos method that efficiently well on JUBL, Jülich's Blue Gene/L. Using the code we treated two-orbital Hubbard Models which show interesting Mott-band insulator transitions. Kohn's criterium which determines the phase of a system by ground-state properties alone, namely the response of the ground-state energies to a vector potential, was employed. This vector potential can be equivalently regarded as a change in boundary conditions, requiring complex wave vectors. We found that half-filled systems become metallic in the region between a band and Mott insulating phase. To extrapolate to infinite system sizes and therefore getting access to an arbitrarily high resolution in reciprocal space, we use cluster perturbation theory (CPT). It works by diagonalizing a finite cluster exactly and then treating the hopping between identical clusters perturbatively in strong coupling perturbation theory. This leads to an effectively infinite chain. Due to our efficient Lanczos implementation it is possible to calculate relatively large clusters, minimizing finite-size effects. The combination of our Lanczos code and CPT allowed for a treatment of the the organic metal TTF-TCNQ. We were able to solve a long-standing problem in the interpretation of its experimental spectra. Being an interesting one-dimensional metal TTF-TCNQ has been studied thoroughly for more than 30 years. Theoretical estimates as well as experimental measurements suggest values for the parameters of the Hubbard model t and U. With these values Hubbard model calculations show qualitatively the same features as APRES experiments, namely the signatures of spin-charge separation at the Gamma-point. But in order to fit the experimental data to the numerical calculations t has to be doubled, which translates into a broadening of the spectra. This ad-hoc adjustment is not only unsatisfactory but also the temperature dependence of the spectral function is incorrectly reproduced. We resolved this problem by properly including nearest-neighbor Coulomb interactions. Recent DFT calculations show that indeed the commonly used values of t and U are correct. But Coulomb interactions between electrons on neighboring sites are also significant and must not be neglected. Thus, we studied the effect of the next-neighbor interaction parameter V being at about U/2. Including V in our Lanczos calculations for obtaining the angular-integrated spectral function shows that its effect is to broaden the spectra, comparable to increasing the value of t. Moreover, we found that doubling t indeed mimics the effect of Vapprox U/2. Surprisingly, the broadening of the spectrum as a result of V can be understood in first-order Rayleigh-Schrödinger perturbation theory.
We present a powerful method for calculating the thermodynamic properties of the Hubbard model in infinite dimensions, using an exact diagonalization of an Anderson model with a finite number of sites. At finite temperatures, the explicit diagonalization of the Anderson Hamiltonian allows the calculation of Green's functions with a resolution far superior to that of Quantum Monte Carlo calculations. At zero temperature, the Lanczòs method is used and yields the essentially exact zero-temperature solution of the model, except in a region of very small frequencies. Numerical results for the half-filled case in the paramagnetic phase (quasi-particle weight, self-energy, and also real-frequency spectral densities) are presented.
The ground-state properties of the periodic Anderson model on one-dimensional chains from three to six sites are studied by exact diagonalization. An appropriate average over boundary conditions is found to be necessary in order to eliminate finite-size oscillations. In the weak-coupling limit, the results for the short-range correlations suggest a smooth transition from a ferromagnetic state at low electron concentration to an antiferromagnetic state at half-filling, rather than a sudden jump at a critical concentration.
A density-functional theory for many-body lattice models is considered in which the single-particle density matrix ${\ensuremath{\gamma}}_{\mathrm{ij}}$ is the basic variable. Eigenvalue equations are derived for solving Levy's constrained search of the interaction energy functional $W[{\ensuremath{\gamma}}_{\mathrm{ij}}].$ $W[{\ensuremath{\gamma}}_{\mathrm{ij}}]$ is expressed as the sum of Hartree-Fock energy ${E}_{\mathrm{HF}}[{\ensuremath{\gamma}}_{\mathrm{ij}}]$ and the correlation energy ${E}_{C}[{\ensuremath{\gamma}}_{\mathrm{ij}}].$ Exact results are obtained for ${E}_{C}({\ensuremath{\gamma}}_{12})$ of the Hubbard model on various periodic lattices, where ${\ensuremath{\gamma}}_{\mathrm{ij}}={\ensuremath{\gamma}}_{12}$ for all nearest neighbors i and j. The functional dependence of ${E}_{C}({\ensuremath{\gamma}}_{12})$ is analyzed by varying the number of sites Na,{N}_{a},Na​, band filling Ne,{N}_{e},Ne​, and lattice structure. The infinite one-dimensional chain and one-, two-, or three-dimensional finite clusters with periodic boundary conditions are considered. The properties of ${E}_{C}({\ensuremath{\gamma}}_{12})$ are discussed in the limits of weak $({\ensuremath{\gamma}}_{12}\ensuremath{\simeq}{\ensuremath{\gamma}}_{12}^{0})$ and strong $({\ensuremath{\gamma}}_{12}\ensuremath{\simeq}{\ensuremath{\gamma}}_{12}^{\ensuremath{\infty}})$ electronic correlations, and in the crossover region $({\ensuremath{\gamma}}_{12}^{\ensuremath{\infty}}<~{\ensuremath{\gamma}}_{12}<~{\ensuremath{\gamma}}_{12}^{0}).$ Using an appropriate scaling we observe that ${\ensuremath{\varepsilon}}_{C}{(g}_{12}{)=E}_{C}{/E}_{\mathrm{HF}}$ has a pseudo-universal behavior as a function of ${g}_{12}=({\ensuremath{\gamma}}_{12}\ensuremath{-}{\ensuremath{\gamma}}_{12}^{\ensuremath{\infty}})/({\ensuremath{\gamma}}_{12}^{0}\ensuremath{-}{\ensuremath{\gamma}}_{12}^{\ensuremath{\infty}}).$ The fact that ${\ensuremath{\varepsilon}}_{C}{(g}_{12})$ depends weakly on Na,{N}_{a},Na​, Ne,{N}_{e},Ne​, and lattice structure suggests that the correlation energy of extended systems could be obtained quite accurately from finite-cluster calculations. Finally, the behaviors of ${E}_{C}({\ensuremath{\gamma}}_{12})$ for repulsive (U>0)(U>0)(U>0) and attractive (U<0)(U<0)(U<0) interactions are contrasted.
This thesis addresses the question of a strongly interacting many-body problem: the Hubbard model on one- and two-dimensional lattices. The technique employed is the small-cluster method wherein the full many-body Hamiltonian for finite clusters with periodic boundary conditions (PBC's) is solved exactly. We first apply this technique to an extended one-dimensional chain and a two-dimensional square-lattice Hubbard model. We then construct a two-dimensional overlayer/substrate Hubbard model with first- and second-nearest-neighbor hopping, on-site and nearest-neighbor interactions, and overlayer-site orbital energy. The ground-state properties are investigated and discussed in a many-body picture.