Publications

A saddle point plus fluctuation analysis of the periodically driven half-filled two-dimensional Hubbard model is performed. For drive frequencies below the equilibrium gap, we find discontinuous transitions to time-dependent solutions. A highly excited, generically nonthermal distribution of magnons occurs even for drive frequencies far above the gap. Above a critical drive amplitude, the low-energy magnon distribution diverges as the frequency tends to zero and antiferromagnetism is destroyed, revealing the generic importance of collective mode excitations arising from a nonequilibrium drive.

The nonuniform fast Fourier transform (NUFFT) generalizes the FFT to off-grid data. Its many applications include image reconstruction, data analysis, and the numerical solution of differential equations. We present FINUFFT, an efficient parallel library for type 1 (nonuiform to uniform), type 2 (uniform to nonuniform), or type 3 (nonuniform to nonuniform) transforms, in dimensions 1, 2, or 3. It uses minimal RAM, requires no precomputation or plan steps, and has a simple interface to several languages. We perform the expensive spreading/interpolation between nonuniform points and the fine grid via a simple new kernel---the `exponential of semicircle' $e^{\beta \sqrt{1-x^2}}$ in $x\in[-1,1]$---in a cache-aware load-balanced multithreaded implementation. The deconvolution step requires the Fourier transform of the kernel, for which we propose efficient numerical quadrature. For types 1 and 2, rigorous error bounds asymptotic in the kernel width approach the fastest known exponential rate, namely that of the Kaiser--Bessel kernel. We benchmark against several popular CPU-based libraries, showing favorable speed and memory footprint, especially in three dimensions when high accuracy and/or clustered point distributions are desired.

Resistance in superconductors arises from the motion of vortices driven by flowing supercurrents or external electromagnetic fields and may be strongly affected by thermal or quantum fluctuations. The common expectation borne out in previous experiments is that as the temperature is lowered, vortex motion is suppressed, leading to a decreased resistance. A new generation of materials provides access to the previously inaccessible regime of clean-limit superconductivity in atomically thin superconducting layers. We show experimentally that for few-layer 2H-NbSe2 the resistance below the superconducting transition temperature may be non-monotonic, passing through a minimum and then increasing again as temperature is decreased further. The effects exists over a wide range of current and magnetic fields, but is most pronounced in monolayer devices at intermediate currents. Analytical and numerical calculations confirm that the findings can be understood in a two-fluid vortex model, in which a fraction of vortices flow in channels while the rest are pinned but thermally fluctuating in position. We show theoretically that the pinned, fluctuating vortices effectively control the mobility of the free vortices. The findings provide a new perspective on fundamental questions of vortex mobility and dissipation in superconductors.

Lecture notes of the Autumn School on Correlated Electrons 2019
1 Xavier Blase: Introduction to Density Functional Theory
2 Xinguo Ren: The Random Phase Approximation and its Application to Real Materials
3 Cyrus Umrigar: Introduction to Variational and Projector Monte Carlo
4 Arne Lüchow: Optimized Quantum Monte Carlo Wave Functions
5 Federico Becca: Variational Wave Functions for Strongly Correlated Fermionic Systems
6 Shiwei Zhang: Auxiliary-Field Quantum Monte Carlo at Zero- and Finite-Temperature
7 Erik Koch: Exact Diagonalization and Lanczos Method
8 Miles Stoudenmire: Quantum Chemistry DMRG in a Local Basis
9 Karen Hallberg: Density Matrix Renormalization
10 Marcelo Rozenberg: Dynamical Mean-Field Theory and Mott Transition
11 Eva Pavarini: Dynamical Mean-Field Theory for Materials
12 Robert Eder: Analytic Properties of Self-Energy and Luttinger-Ward Functional
13 James Freericks: Introduction to Many-Body Green Functions In and Out Of Equilibrium
14 Andrea Donarini: Electronic Transport in Correlated Single Molecule Junctions
15 Nikolay Prokof'ev: Diagrammatic Monte Carlo
16 Anders Sandvik: Stochastic Series Expansion Methods
17 Gerardo Ortiz: Algebraic Methods in Many-Body Physics

We propose a hierarchical approach to testing general relativity with multiple gravitational wave detections. Unlike existing strategies, our method does not assume that parameters quantifying deviations from general relativity are either common or completely unrelated across all sources. We instead assume that these parameters follow some underlying distribution, which we parametrize and constrain. This can be then compared to the distribution expected from general relativity, i.e. no deviation in any of the events. We demonstrate that our method is robust to measurement uncertainties and can be applied to theories of gravity where the parameters beyond general relativity are related to each other, as generally expected. Our method contains the two extremes of common and unrelated parameters as limiting cases. We apply the hierarchical model to the population of 10 binary black hole systems so far detected by LIGO and Virgo. We do this for a parametrized test of gravitational wave generation, by modeling the population distribution of beyond-general-relativity parameters with a Gaussian distribution. We compute the mean and the variance of the population and show that both are consistent with general relativity for all parameters we consider. In the best case, we find that the population properties of the existing binary signals are consistent with general relativity at the ~1% level. This hierarchical approach subsumes and extends existing methodologies, and is more robust at revealing potential subtle deviations from general relativity with increasing number of detections.

We study the magnetic susceptibility in the normal state of Sr2RuO4 using dynamical mean-field theory including dynamical vertex corrections. Besides the well known incommensurate response, our calculations yield quasi-local spin fluctuations which are broad in momentum and centered around the Γ point, in agreement with recent inelastic neutron scattering experiments [P. Steffens, et al., Phys. Rev. Lett. 122, 047004 (2019)]. We show that these quasi-local fluctuations are controlled by the Hund's coupling and account for the dominant contribution to the momentum-integrated response. While all orbitals contribute equally to the incommensurate response, the enhanced Γ point response originates from the planar xy orbital.