Publications

The intersection of electronic topology and strong correlations offers a rich platform to discover exotic quantum phases of matter and unusual materials. An overarching challenge that impedes the discovery is how to diagnose topology in strongly correlated settings, as exemplified by Mott insulators. Here, we develop a general framework to address this outstanding question and illustrate its power in the case of Mott insulators. The concept of Green's function Berry curvature -- which is frequency dependent -- is introduced. We apply this notion in a system that contains symmetry-protected nodes in its noninteracting bandstructure; strong correlations drive the system into a Mott insulating state, creating contours in frequency-momentum space where the Green's function vanishes. The Green's function Berry flux of such zeros is found to be quantized, and is as such direct probe of the system's topology. Our framework allows for a comprehensive search of strongly correlated topological materials with Green's function topology.

The Matsubara Green's function formalism stands as a powerful technique for computing the thermodynamic characteristics of interacting quantum many-particle systems at finite temperatures. In this manuscript, our focus centers on introducing MatsubaraFunctions.jl, a Julia library that implements data structures for generalized n-point Green's functions on Matsubara frequency grids. The package's architecture prioritizes user-friendliness without compromising the development of efficient solvers for quantum field theories in equilibrium. Following a comprehensive introduction of the fundamental types, we delve into a thorough examination of key facets of the interface. This encompasses avenues for accessing Green's functions, techniques for extrapolation and interpolation, as well as the incorporation of symmetries and a variety of parallelization strategies. Examples of increasing complexity serve to demonstrate the practical utility of the library, supplemented by discussions on strategies for sidestepping impediments to optimal performance.

A major challenge in the field of correlated electrons is the computation of dynamical correlation functions. For comparisons with experiment, one is interested in their real-frequency dependence. This is difficult to compute, as imaginary-frequency data from the Matsubara formalism require analytic continuation, a numerically ill-posed problem. Here, we apply quantum field theory to the single-impurity Anderson model (AM), using the Keldysh instead of the Matsubara formalism with direct access to the self-energy and dynamical susceptibilities on the real-frequency axis. We present results from the functional renormalization group (fRG) at one-loop level and from solving the self-consistent parquet equations in the parquet approximation. In contrast to previous Keldysh fRG works, we employ a parametrization of the four-point vertex which captures its full dependence on three real-frequency arguments. We compare our results to benchmark data obtained with the numerical renormalization group and to second-order perturbation theory. We find that capturing the full frequency dependence of the four-point vertex significantly improves the fRG results compared to previous implementations, and that solving the parquet equations yields the best agreement with the NRG benchmark data, but is only feasible up to moderate interaction strengths. Our methodical advances pave the way for treating more complicated models in the future.

Kohn's theorem places strong constraints on the cyclotron response of Fermi liquids. Recent observations of a doping dependence in the cyclotron mass of La

We present a method of calculating the energy gap of a charge-neutral excitation using only ground-state calculations. We report Quantum Monte Carlo calculations of Γ→Γ and Γ→X particle-hole excitation energies in diamond carbon. We analyze the finite-size effect and find the same 1/L decay rate as that in a charged excitation, where L is the linear extension of the supercell. This slow decay is attributed to the delocalized nature of the excitation in supercells too small to accommodate excitonic binding effects. At larger system sizes, the apparent 1/L decay crosses over to a 1/L

Variational wave function ansatze are an invaluable tool to study the properties of strongly correlated systems. We propose such a wave function, based on the theory of auxiliary fields and combining aspects of auxiliary-field quantum Monte Carlo and modern variational optimization techniques including automatic differentiation. The resulting ansatz, consisting of several slices of optimized projectors, is highly expressive and systematically improvable. We benchmark this form on the two-dimensional Hubbard model, using both cylindrical and large, fully periodic supercells. The computed ground-state energies are competitive with the best variational results. Moreover, the optimized wave functions predict the correct ground-state order with near full symmetry restoration (i.e. translation invariance) despite initial states with incorrect orders. The ansatz can become a tool for local order prediction, leading to a new paradigm for variational studies of bulk systems. It can also be viewed as an approach to produce accurate and systematically improvable wave functions in a convenient form of non-orthogonal Slater determinants (e.g., for quantum chemistry) at polynomial computational cost.

Semiconducting transition-metal dichalcogenides (TMDs) exhibit high mobility, strong spin-orbit coupling, and large effective masses, which simultaneously leads to a rich wealth of Landau quantizations and inherently strong electronic interactions. However, in spite of their extensively explored Landau levels (LL) structure, probing electron correlations in the fractionally filled LL regime has not been possible due to the difficulty of reaching the quantum limit. Here, we report evidence for fractional quantum Hall (FQH) states at filling fractions 4/5 and 2/5 in the lowest LL of bilayer MoS

Synthetic dimension platforms offer unique pathways for engineering quantum matter. We compute the phase diagram of a many-body system of ultracold atoms (or polar molecules) with a set of Rydberg states (or rotational states) as a synthetic dimension, where the particles are arranged in real space in optical microtrap arrays and interact via dipole-dipole exchange interaction. Using mean-field theory, we find three ordered phases - two are localized in the synthetic dimension, predicted as zero-temperature ground states in Refs. [Sci. Rep., 8, 1 (2018) and Phys. Rev. A 99, 013624 (2019)], and a delocalized phase. We characterize them by identifying the spontaneously broken discrete symmetries of the Hamiltonian. We also compute the phase diagram as a function of temperature and interaction strength, for both signs of the interaction. For system sizes with more than six synthetic sites and attractive interactions, we find that the thermal phase transitions can be first or second order, which leads to a tri-critical point on the phase boundary. By examining the dependence of the tri-critical point and other special points of the phase boundary on the synthetic dimension size, we shed light on the physics for thermodynamically large synthetic dimension.

Strongly interacting electronic systems possess rich phase diagrams resulting from the competition between different quantum ground states. A general mechanism that relieves this frustration is the emergence of microemulsion phases, where regions of different phase self-organize across multiple length scales. The experimental characterization of these phases often poses significant challenges, as the long-range Coulomb interaction microscopically mingles the competing states. Here, we use cryogenic reflectance and magneto-optical spectroscopy to observe the signatures of the mixed state between an electronic Wigner crystal and an electron liquid in a MoSe2 monolayer. We find that the transit into this 'microemulsion' state is marked by anomalies in exciton reflectance, spin susceptibility, and Umklapp scattering, establishing it as a distinct phase of electronic matter. Our study of the two-dimensional electronic microemulsion phase elucidates the physics of novel correlated electron states with strong Coulomb interactions.

Over the last two decades, following the early developments on maximally-localized Wannier functions, an ecosystem of electronic-structure simulation techniques and software leveraging the Wannier representation has flourished. This environment includes codes to obtain Wannier functions and interfaces with first-principles simulation software, as well as an increasing number of related post-processing packages. Wannier functions can be obtained for isolated or extended systems (both crystalline and disordered), and can be used to understand chemical bonding, to characterize polarization, magnetization, and topology, or as an optimal basis set, providing very accurate interpolations in reciprocal space or large-scale Hamiltonians in real space. In this review, we summarize the current landscape of techniques, materials properties and simulation codes based on Wannier functions that have been made accessible to the research community, and that are now well integrated into what we term a