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

Charge dynamics of correlated electrons formulated by variational description with inclusion of composite fermions

Kota Ido, Masatoshi Imada, Takahiro Misawa

Published Jul 9, 2019• 0 citations
OpenAlex

Abstract

We propose a method to calculate the charge dynamical structure factors for the ground states of correlated electron systems based on the variational Monte Carlo method. Our benchmarks for the one- and two-dimensional Hubbard models show that inclusion of composite-fermion excitations in the basis set greatly improves the accuracy, in reference to the exact charge dynamical structure factors for clusters. Together with examination for larger systems beyond tractable sizes by the exact diagonalization, our results indicate that the variational Monte Carlo method is a promising way for studies on the nature of charge dynamics in correlated materials such as the copper oxide superconductors, if the composite-fermion excitations are properly included in the restricted Hilbert space of intermediate states in the linear response theory. We also discuss the importance of incorporating nonlocal composite fermion for more accurate description.

Related papers

Charge dynamics of correlated electrons: Variational description with inclusion of composite fermions

Kota Ido, Masatoshi Imada, Takahiro Misawa

We propose a method to calculate the charge dynamical structure factors for the ground states of correlated electron systems based on the variational Monte Carlo method. Our benchmarks for the one- and two-dimensional Hubbard models show that inclusion of composite-fermion excitations in the basis set greatly improves the accuracy, in reference to the exact charge dynamical structure factors for clusters. Together with examination for larger systems beyond tractable sizes by the exact diagonalization, our results indicate that the variational Monte Carlo method is a promising way for studies on the nature of charge dynamics in correlated materials such as the copper oxide superconductors if the composite-fermion excitations are properly included in the restricted Hilbert space of intermediate states in the linear response theory. Our results are consistent with the particle-hole excitations inferred from the single-particle spectral function $A(\mathbit{k},\ensuremath{\omega})$ in the literature. We also discuss the importance of incorporating nonlocal composite fermion for a more accurate description. Future issues for further improvements are also discussed.

arXiv:1907.04123•Physical review. B./Physical review. B•Feb 18, 2020•14 citations

Variational Monte Carlo method for fermionic models combined with tensor networks and applications to the hole-doped two-dimensional Hubbard model

Hui‐Hai Zhao, Kota Ido, Satoshi Morita +1 more

The conventional tensor-network states employ real-space product states as reference wave functions. Here, we propose a many-variable variational Monte Carlo (mVMC) method combined with tensor networks by taking advantages of both to study fermionic models. The variational wave function is composed of a pair product wave function operated by real-space correlation factors and tensor networks. Moreover, we can apply quantum number projections, such as spin, momentum, and lattice symmetry projections, to recover the symmetry of the wave function to further improve the accuracy. We benchmark our method for one- and two-dimensional Hubbard models, which show significant improvement over the results obtained individually either by mVMC or by tensor network. We have applied the present method to a hole-doped Hubbard model on the square lattice, which indicates the stripe charge/spin order coexisting with a weak $d$-wave superconducting order in the ground state for the doping concentration of less than 0.3, where the stripe oscillation period gets longer with increasing hole concentration. The charge homogeneous and highly superconducting state also exists as a metastable excited state for the doping concentration less than 0.25.

arXiv:1703.03537•Physical review. B./Physical review. B•Aug 1, 2017•42 citations

Dynamical-quantum-cluster approach to two-particle correlation functions in the Hubbard model

S. Hochkeppel, Fakher F. Assaad, Werner R. Hanke

We investigate the charge- and spin-dynamical structure factors for the two-dimensional one-band Hubbard model in the strong coupling regime within the dynamical cluster approximation (DCA) that is explicitly extended to two-particle response functions. The full irreducible two-particle vertex with three momenta and frequencies is approximated by an effective vertex dependent on the momentum and frequency of the spin and/or charge excitation. In the spirit of the DCA, the effective vertex is calculated by using the quantum Monte Carlo methods on a finite cluster. On the basis of a comparison to high temperature auxiliary field quantum Monte Carlo data, we show that near and beyond optimal doping, our results provide a consistent overall picture of the interplay between charge, spin, and single-particle excitations.

arXiv:0801.3377•Physical Review B•May 13, 2008•11 citations

Time-dependent many-variable variational Monte Carlo method for nonequilibrium strongly correlated electron systems

Kota Ido, Takahiro Ohgoe, Masatoshi Imada

We develop a time-dependent variational Monte Carlo (t-VMC) method for quantum dynamics of strongly correlated electrons. The t-VMC method has been recently applied to bosonic systems and quantum spin systems. Here we propose a time-dependent trial wave function with many variational parameters, which is suitable for nonequilibrium strongly correlated electron systems. As the trial state, we adopt the generalized pair-product wave function with correlation factors and quantum-number projections. This trial wave function has been proven to accurately describe ground states of strongly correlated electron systems. To show the accuracy and efficiency of our trial wave function in nonequilibrium states as well, we present our benchmark results for relaxation dynamics during and after interaction quench protocols of fermionic Hubbard models. We find that our trial wave function well reproduces the exact results for the time evolution of physical quantities such as energy, momentum distribution, spin structure factor, and superconducting correlations. These results show that the t-VMC with our trial wave function offers an efficient and accurate way to study challenging problems of nonequilibrium dynamics in strongly correlated electron systems.

arXiv:1507.00274•Physical Review B•Dec 7, 2015•48 citations

Evaluation of the excitation spectra with diffusion Monte Carlo on an auxiliary bosonic ground state

Fernando Agustín Reboredo, Paul R. C. Kent, Jaron T. Krogel

We aim to improve upon the variational Monte Carlo (VMC) approach for excitations replacing the Jastrow factor by an auxiliary bosonic (AB) ground state and multiplying it by a fermionic component factor. The instantaneous change in imaginary time of an arbitrary excitation in the original interacting fermionic system is obtained by measuring observables via the ground-state distribution of walkers of an AB system that is subject to an auxiliary effective potential. The effective potential is used to (i) drive the AB system's ground-state configuration space toward the configuration space of the excitations of the original fermionic system and (ii) subtract from a diffusion Monte Carlo (DMC) calculation contributions that can be included in conventional approximations, such as mean-field and configuration interaction (CI) methods. In this novel approach, the AB ground state is treated statistically in DMC, whereas the fermionic component of the original system is expanded in a basis. The excitation energies of the fermionic eigenstates are obtained by sampling a fermion-boson coupling term on the AB ground state. We show that this approach can take advantage of and correct for approximate eigenstates obtained via mean-field calculations or truncated interactions. We demonstrate that the AB ground-state factor incorporates the correlations missed by standard Jastrow factors, further reducing basis truncation errors. Relevant parts of the theory have been tested in soluble model systems and exhibit excellent agreement with exact analytical data and CI and VMC approaches. In particular, for limited basis set expansions and sufficient statistics, AB approaches outperform CI and VMC in terms of basis size for the same systems. The implementation of this method in current codes, despite being demanding, will be facilitated by reusing procedures already developed for calculating ground-state properties with DMC and excitations with VMC.

The Journal of Chemical Physics•Sep 19, 2023•1 citation

Charge ordering in extended Hubbard models: Variational cluster approach

Markus Aichhorn, Hans Gerd Evertz, Wolfgang von der Linden +1 more

We present a generalization of the recently proposed variational cluster perturbation theory to extended Hubbard models at half-filling with repulsive nearest neighbor interaction. The method takes into account short-range correlations correctly by the exact diagonalization of clusters of finite size, whereas long-range order beyond the size of the clusters is treated on a mean-field level. For one dimension, we show that quantum Monte Carlo and density-matrix renormalization-group results can be reproduced with very good accuracy. Moreover we apply the method to the two-dimensional extended Hubbard model on a square lattice. In contrast to the one-dimensional case, a first order phase transition between spin density wave phase and charge density wave phase is found as function of the nearest-neighbor interaction at onsite interactions $U\ensuremath{\geqslant}3t$. The single-particle spectral function is calculated for both the one-dimensional and the two-dimensional system.

arXiv:cond-mat/0402580•Physical Review B•Dec 3, 2004•99 citations
We propose a method to calculate the charge dynamical structure factors for the ground states of correlated electron systems based on the variational Monte Carlo method. Our benchmarks for the one- and two-dimensional Hubbard models show that inclusion of composite-fermion excitations in the basis set greatly improves the accuracy, in reference to the exact charge dynamical structure factors for clusters. Together with examination for larger systems beyond tractable sizes by the exact diagonalization, our results indicate that the variational Monte Carlo method is a promising way for studies on the nature of charge dynamics in correlated materials such as the copper oxide superconductors, if the composite-fermion excitations are properly included in the restricted Hilbert space of intermediate states in the linear response theory. We also discuss the importance of incorporating nonlocal composite fermion for more accurate description.
We propose a method to calculate the charge dynamical structure factors for the ground states of correlated electron systems based on the variational Monte Carlo method. Our benchmarks for the one- and two-dimensional Hubbard models show that inclusion of composite-fermion excitations in the basis set greatly improves the accuracy, in reference to the exact charge dynamical structure factors for clusters. Together with examination for larger systems beyond tractable sizes by the exact diagonalization, our results indicate that the variational Monte Carlo method is a promising way for studies on the nature of charge dynamics in correlated materials such as the copper oxide superconductors if the composite-fermion excitations are properly included in the restricted Hilbert space of intermediate states in the linear response theory. Our results are consistent with the particle-hole excitations inferred from the single-particle spectral function $A(\mathbit{k},\ensuremath{\omega})$ in the literature. We also discuss the importance of incorporating nonlocal composite fermion for a more accurate description. Future issues for further improvements are also discussed.
The conventional tensor-network states employ real-space product states as reference wave functions. Here, we propose a many-variable variational Monte Carlo (mVMC) method combined with tensor networks by taking advantages of both to study fermionic models. The variational wave function is composed of a pair product wave function operated by real-space correlation factors and tensor networks. Moreover, we can apply quantum number projections, such as spin, momentum, and lattice symmetry projections, to recover the symmetry of the wave function to further improve the accuracy. We benchmark our method for one- and two-dimensional Hubbard models, which show significant improvement over the results obtained individually either by mVMC or by tensor network. We have applied the present method to a hole-doped Hubbard model on the square lattice, which indicates the stripe charge/spin order coexisting with a weak ddd-wave superconducting order in the ground state for the doping concentration of less than 0.3, where the stripe oscillation period gets longer with increasing hole concentration. The charge homogeneous and highly superconducting state also exists as a metastable excited state for the doping concentration less than 0.25.
We investigate the charge- and spin-dynamical structure factors for the two-dimensional one-band Hubbard model in the strong coupling regime within the dynamical cluster approximation (DCA) that is explicitly extended to two-particle response functions. The full irreducible two-particle vertex with three momenta and frequencies is approximated by an effective vertex dependent on the momentum and frequency of the spin and/or charge excitation. In the spirit of the DCA, the effective vertex is calculated by using the quantum Monte Carlo methods on a finite cluster. On the basis of a comparison to high temperature auxiliary field quantum Monte Carlo data, we show that near and beyond optimal doping, our results provide a consistent overall picture of the interplay between charge, spin, and single-particle excitations.
We develop a time-dependent variational Monte Carlo (t-VMC) method for quantum dynamics of strongly correlated electrons. The t-VMC method has been recently applied to bosonic systems and quantum spin systems. Here we propose a time-dependent trial wave function with many variational parameters, which is suitable for nonequilibrium strongly correlated electron systems. As the trial state, we adopt the generalized pair-product wave function with correlation factors and quantum-number projections. This trial wave function has been proven to accurately describe ground states of strongly correlated electron systems. To show the accuracy and efficiency of our trial wave function in nonequilibrium states as well, we present our benchmark results for relaxation dynamics during and after interaction quench protocols of fermionic Hubbard models. We find that our trial wave function well reproduces the exact results for the time evolution of physical quantities such as energy, momentum distribution, spin structure factor, and superconducting correlations. These results show that the t-VMC with our trial wave function offers an efficient and accurate way to study challenging problems of nonequilibrium dynamics in strongly correlated electron systems.
We aim to improve upon the variational Monte Carlo (VMC) approach for excitations replacing the Jastrow factor by an auxiliary bosonic (AB) ground state and multiplying it by a fermionic component factor. The instantaneous change in imaginary time of an arbitrary excitation in the original interacting fermionic system is obtained by measuring observables via the ground-state distribution of walkers of an AB system that is subject to an auxiliary effective potential. The effective potential is used to (i) drive the AB system's ground-state configuration space toward the configuration space of the excitations of the original fermionic system and (ii) subtract from a diffusion Monte Carlo (DMC) calculation contributions that can be included in conventional approximations, such as mean-field and configuration interaction (CI) methods. In this novel approach, the AB ground state is treated statistically in DMC, whereas the fermionic component of the original system is expanded in a basis. The excitation energies of the fermionic eigenstates are obtained by sampling a fermion-boson coupling term on the AB ground state. We show that this approach can take advantage of and correct for approximate eigenstates obtained via mean-field calculations or truncated interactions. We demonstrate that the AB ground-state factor incorporates the correlations missed by standard Jastrow factors, further reducing basis truncation errors. Relevant parts of the theory have been tested in soluble model systems and exhibit excellent agreement with exact analytical data and CI and VMC approaches. In particular, for limited basis set expansions and sufficient statistics, AB approaches outperform CI and VMC in terms of basis size for the same systems. The implementation of this method in current codes, despite being demanding, will be facilitated by reusing procedures already developed for calculating ground-state properties with DMC and excitations with VMC.
We present a generalization of the recently proposed variational cluster perturbation theory to extended Hubbard models at half-filling with repulsive nearest neighbor interaction. The method takes into account short-range correlations correctly by the exact diagonalization of clusters of finite size, whereas long-range order beyond the size of the clusters is treated on a mean-field level. For one dimension, we show that quantum Monte Carlo and density-matrix renormalization-group results can be reproduced with very good accuracy. Moreover we apply the method to the two-dimensional extended Hubbard model on a square lattice. In contrast to the one-dimensional case, a first order phase transition between spin density wave phase and charge density wave phase is found as function of the nearest-neighbor interaction at onsite interactions $U\ensuremath{\geqslant}3t$. The single-particle spectral function is calculated for both the one-dimensional and the two-dimensional system.