Search for: All records

We introduce nested gausslet bases, an improvement on previous gausslet bases that can treat systems containing atoms with much larger atomic numbers. We also introduce pure Gaussian distorted gausslet bases, which allow the Hamiltonian integrals to be performed analytically, as well as hybrid bases in which the gausslets are combined with standard Gaussian-type bases. All these bases feature the diagonal approximation for the electron–electron interactions so that the Hamiltonian is completely defined by two Nb × Nb matrices, where Nb ≈ 104 is small enough to permit fast calculations at the Hartree–Fock level. In constructing these bases, we have gained new mathematical insight into the construction of one-dimensional diagonal bases. In particular, we have proved an important theorem relating four key basis set properties: completeness, orthogonality, zero-moment conditions, and diagonalization of the coordinate operator matrix. We test our basis sets on small systems with a focus on high accuracy, obtaining, for example, an accuracy of 2 × 10−5 Ha for the total Hartree–Fock energy of the neon atom in the complete basis set limit. 

Abstract A key component of the phase diagram of many iron-based superconductors and electron-doped cuprates is believed to be a quantum critical point (QCP), delineating the onset of antiferromagnetic spin-density wave order in a quasi-two-dimensional metal. The universality class of this QCP is believed to play a fundamental role in the description of the proximate non-Fermi liquid behavior and superconducting phase. A minimal model for this transition is the O(3) spin-fermion model. Despite many efforts, a definitive characterization of its universal properties is still lacking. Here, we numerically study the O(3) spin-fermion model and extract the scaling exponents and functional form of the static and zero-momentum dynamical spin susceptibility. We do this using a Hybrid Monte Carlo (HMC) algorithm with a novel auto-tuning procedure, which allows us to study unprecedentedly large systems of 80 × 80 sites. We find a strong violation of the Hertz-Millis form, contrary to all previous numerical results. Furthermore, the form that we do observe provides good evidence that the universal scaling is actually governed by the analytically tractable fixed point discovered near perfect “hot-spot’ nesting, even for a larger nesting window. Our predictions can be directly tested with neutron scattering. Additionally, the HMC method we introduce is generic and can be used to study other fermionic models of quantum criticality, where there is a strong need to simulate large systems.