Researcher profile

Richard Scalettar

Richard Scalettar contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

5 published item(s)

preprint2022arXiv

Dynamical tuning of the chemical potential to achieve a target particle number in grand canonical Monte Carlo simulations

We present a method to facilitate Monte Carlo simulations in the grand canonical ensemble given a target mean particle number. The method imposes a fictitious dynamics on the chemical potential, to be run concurrently with the Monte Carlo sampling of the physical system. Corrections to the chemical potential are made according to time-averaged estimates of the mean and variance of the particle number, with the latter being proportional to thermodynamic compressibility. We perform a variety of tests, and in all cases find rapid convergence of the chemical potential -- inexactness of the tuning algorithm contributes only a minor part of the total measurement error for realistic simulations.

preprint2022arXiv

Fast and scalable quantum Monte Carlo simulations of electron-phonon models

We introduce methodologies for highly scalable quantum Monte Carlo simulations of electron-phonon models, and report benchmark results for the Holstein model on the square lattice. The determinant quantum Monte Carlo (DQMC) method is a widely used tool for simulating simple electron-phonon models at finite temperatures, but incurs a computational cost that scales cubically with system size. Alternatively, near-linear scaling with system size can be achieved with the hybrid Monte Carlo (HMC) method and an integral representation of the Fermion determinant. Here, we introduce a collection of methodologies that make such simulations even faster. To combat "stiffness" arising from the bosonic action, we review how Fourier acceleration can be combined with time-step splitting. To overcome phonon sampling barriers associated with strongly-bound bipolaron formation, we design global Monte Carlo updates that approximately respect particle-hole symmetry. To accelerate the iterative linear solver, we introduce a preconditioner that becomes exact in the adiabatic limit of infinite atomic mass. Finally, we demonstrate how stochastic measurements can be accelerated using fast Fourier transforms. These methods are all complementary and, combined, may produce multiple orders of magnitude speedup, depending on model details.

preprint2022arXiv

Phase Diagram of the Su-Schrieffer-Heeger-Hubbard model on a square lattice

The Hubbard and Su-Schrieffer-Heeger Hamiltonians (SSH) are iconic models for understanding the qualitative effects of electron-electron and electron-phonon interactions respectively. In the two-dimensional square lattice Hubbard model at half filling, the on-site Coulomb repulsion, $U$, between up and down electrons induces antiferromagnetic (AF) order and a Mott insulating phase. On the other hand, for the SSH model, there is an AF phase when the electron-phonon coupling $λ$ is less than a critical value $λ_c$ and a bond order wave when $λ> λ_c$. In this work, we perform numerical studies on the square lattice optical Su-Schrieffer-Heeger-Hubbard Hamiltonian (SSHH), which combines both interactions. We use the determinant quantum Monte Carlo (DQMC) method which does not suffer from the fermionic sign problem at half filling. We map out the phase diagram and find that it exhibits a direct first-order transition between an antiferromagnetic phase and a bond-ordered wave as $λ$ increases. The AF phase is characterized by two different regions. At smaller $λ$ the behavior is similar to that of the pure Hubbard model; the other region, while maintaining long range AF order, exhibits larger kinetic energies and double occupancy, i.e. larger quantum fluctuations, similar to the AF phase found in the pure SSH model.

preprint2022arXiv

Quantum Membrane Phases in Synthetic Lattices of Cold Molecules or Rydberg Atoms

We calculate properties of dipolar interacting ultracold molecules or Rydberg atoms in a semi-synthetic three-dimensional configuration -- one synthetic dimension plus a two-dimensional real space optical lattice or periodic microtrap array -- using the stochastic Green function Quantum Monte Carlo method. Through a calculation of thermodynamic quantities and appropriate correlation functions, along with their finite size scalings, we show that there is a second order transition to a low temperature phase in which two-dimensional `sheets' form in the synthetic dimension of internal rotational or electronic states of the molecules or Rydberg atoms, respectively. Simulations for different values of the interaction $V$, which acts between atoms or molecules that are adjacent both in real and synthetic space, allow us to compute a phase diagram. We find a finite-temperature transition at sufficiently large $V$, as well as a quantum phase transition -- a critical value $V_c$ below which the transition temperature vanishes.

preprint2021arXiv

$π$-Phase shift across stripes in a charge density wave system

Many strongly correlated materials are characterized by deeply intertwined charge and spin order. Besides their high superconducting transition temperatures, one of the central features of these complex patterns in cuprates is a phase shift which occurs across lines of decreased hole density. That is, when doped away from their AF phase, the additional charge is not distributed uniformly, but rather in `stripes'. The sublattices preferentially occupied by up and down spin are reversed across these stripes, a phenomonenon referred to as a `$π$-phase shift'. Many of the spin-charge patterns, including the $π$-phase shift, are reproduced by Density Matrix Renormalization Group and Quantum Monte Carlo calculations of simplified tight binding (repulsive Hubbard) models. In this paper we demonstrate that this sublattice reversal is generic by considering the corresponding phenomenon in the attractive Hubbard Hamiltonian, where a charge density wave phase forms at half-filling. We introduce charge stripes via an appropriate local chemical potential; measurements of charge correlation across the resulting lines of lowered density reveal a clear $π$ phase.