Publications

Condensates are a hallmark of emergence in quantum materials with superconductors and charge density wave as prominent examples. An excitonic insulator (EI) is an intriguing addition to this library, exhibiting spontaneous condensation of electron-hole pairs. However, condensate observables can be obscured through parasitic coupling to the lattice. Time-resolved terahertz (THz) spectroscopy can disentangle such obscurants through measurement of the quantum dynamics. We target Ta

The exchange-only virial relation due to Levy and Perdew is revisited. Invoking the adiabatic connection, we introduce the exchange energy in terms of the right-derivative of the universal density functional w.r.t. the coupling strength λ at λ=0. This agrees with the Levy-Perdew definition of the exchange energy as a high-density limit of the full exchange-correlation energy. By relying on v-representability for a fixed density at varying coupling strength, we prove an exchange-only virial relation without an explicit local-exchange potential. Instead, the relation is in terms of a limit (λ↘0) involving the exchange-correlation potential v

Density functional theory (DFT) offers a desirable balance between quantitative accuracy and computational efficiency in practical many-electron calculations. Its central component, the exchange-correlation energy functional, has been approximated with increasing levels of complexity ranging from strictly local approximations to nonlocal and orbital-dependent expressions with many tuned parameters. In this work, we formulate a general way of rewriting complex density functionals using deep neural networks in a way that allows for simplified computation of Kohn-Sham potentials as well as higher functional derivatives through automatic differentiation, enabling access to highly nonlinear response functions and forces. These goals are achieved by using a recently developed class of robust neural network models capable of modeling functionals, as opposed to functions, with explicitly enforced spatial symmetries. Functionals treated in this way are then called global density approximations and can be seamlessly integrated with existing DFT workflows. Tests are performed for a dataset featuring a large variety of molecular structures and popular meta-GGA density functionals, where we successfully eliminate orbital dependencies coming from the kinetic energy density, and discover a high degree of transferability to a variety of physical systems. The presented framework is general and could be extended to more complex orbital and energy dependent functionals as well as refined with specialized datasets.

Helical spin structures are expressions of magnetically induced chirality, entangling the dipolar and magnetic orders in materials1,2,3,4. The recent discovery of helical van der Waals multiferroics down to the ultrathin limit raises prospects of large chiral magnetoelectric correlations in two dimensions5,6. However, the exact nature and magnitude of these couplings have remained unknown so far. Here we perform a precision measurement of the dynamical magnetoelectric coupling for an enantiopure domain in an exfoliated van der Waals multiferroic. We evaluate this interaction in resonance with a collective electromagnon mode, capturing the impact of its oscillations on the dipolar and magnetic orders of the material with a suite of ultrafast optical probes. Our data show a giant natural optical activity at terahertz frequencies, characterized by quadrature modulations between the electric polarization and magnetization components. First-principles calculations further show that these chiral couplings originate from the synergy between the non-collinear spin texture and relativistic spin–orbit interactions, resulting in substantial enhancements over lattice-mediated effects. Our findings highlight the potential for intertwined orders to enable unique functionalities in the two-dimensional limit and pave the way for the development of van der Waals magnetoelectric devices operating at terahertz speeds.

Ultraclean graphene at charge neutrality hosts a quantum critical Dirac fluid of interacting electrons and holes. Interactions profoundly affect the charge dynamics of graphene, which is encoded in the properties of its electron-photon collective modes: surface plasmon polaritons (SPPs). Here, we show that polaritonic interference patterns are particularly well suited to unveil the interactions in Dirac fluids by tracking polaritonic interference in time at temporal scales commensurate with the electronic scattering. Spacetime SPP interference patterns recorded in terahertz (THz) frequency range provided unobstructed readouts of the group velocity and lifetime of polariton that can be directly mapped onto the electronic spectral weight and the relaxation rate. Our data uncovered prominent departures of the electron dynamics from the predictions of the conventional Fermi-liquid theory. The deviations are particularly strong when the densities of electrons and holes are approximately equal. The proposed spacetime imaging methodology can be broadly applied to probe the electrodynamics of quantum materials. The worldlines of plasmon polaritons in spacetime are visualized, and the electronic correlations in a Dirac fluid are unveiled.

We develop a comprehensive method to construct analytical continuum models for moiré systems directly from first-principle calculations without any parameter fitting. The core idea of this method is to interpret the terms in the continuum model as a basis, allowing us to determine model parameters as coefficients of this basis through Gram-Schmidt orthogonalization. We apply our method to twisted MoTe

Common wisdom dictates that physical systems become less ordered when heated to higher temperature. However, several systems display the opposite phenomenon and move to a more ordered state upon heating, e.g. at low temperature piezoelectric quartz is paraelectric and it only becomes piezoelectric when heated to sufficiently high temperature. The presence, or better, the re-entrance of unordered phases at low temperature is more prevalent than one might think. Although specific models have been developed to understand the phenomenon in specific systems, a universal explanation is lacking. Here we propose a universal simple microscopic theory which predicts the existence of two critical temperatures in inhomogeneous systems, where the lower one marks the re-entrance into the less ordered phase. We show that the re-entrant phase transition is caused by disorder-induced spatial localization of the zero-mode on a finite, i.e. sub-extensive, region of the system. Specifically, this trapping of the zero-mode disconnects the fluctuations of the order parameter in distant regions of the system, thus triggering the loss of long-range order and the re-entrance into the disordered phase. This makes the phenomenon quite universal and robust to the underlying details of the model, and explains its ubiquitous observation.

We adapt the topological quantum chemistry formalism to layer groups, and apply it to study the band topology of 8,872 entries from the computational two-dimensional (2D) materials databases C2DB and MC2D. In our analysis, we find 4,073 topologically non-trivial or obstructed atomic insulator entries, including 905 topological insulators, 602 even-electron number topological semimetals, and 1,003 obstructed atomic insulators. We thus largely expand the library of known topological or obstructed materials in two dimensions, beyond the few hundreds known to date. We additionally classify the materials into four categories: experimentally existing, stable, computationally exfoliated, and not stable. We present a detailed analysis of the edge states emerging in a number of selected new materials, and compile a Topological 2D Materials Database (2D-TQCDB) containing the band structures and detailed topological properties of all the materials studied in this work. The methodology here developed is implemented in new programs available to the public, designed to study the topology of any non-magnetic monolayer or multilayer 2D material.

We describe the evolution of low-temperature thermopower across Fermi-volume-changing quantum phase transitions in Kondo lattice models without translational symmetry breaking. This transition moves from a heavy Fermi liquid with a conventional Luttinger-volume large Fermi surface to a 'FL*' state, characterized by a small Fermi surface and a spin liquid with fractionalized excitations. The onset of the large Fermi surface phase is driven by the condensation of a Higgs field that carries a unit gauge charge under an emergent U(1) gauge field. We consider the case with spatially random Kondo exchange, as this leads to strange metal behavior in electrical transport. We find a large asymmetric thermopower in a 'skewed' marginal Fermi liquid, with similarities to the skewed non-Fermi liquid of Georges and Mravlje (arXiv:2102.13224). Our findings are consistent with recent observations in heavy fermion compounds (Z.-Y. Cao et al., arXiv:2408.13604), and describe an enhancement of thermopower on the large Fermi surface side as well as a non-monotonic behavior on the small Fermi surface side. Our results also apply to single-band Hubbard models and the pseudogap transition in the cuprates. In the ancilla framework, single-band models exhibit an inverted Kondo lattice transition: the small Fermi surface pseudogap state corresponds to the condensed Higgs state. This inversion results in an enhancement of thermopower on the pseudogap side in our theory, consistent with observations in the cuprates (C. Collignon et al., arXiv:2011.14927; A. Gourgout et al., arXiv:2106.05959). We argue that these observations support a non-symmetry-breaking Fermi-volume-changing transition as the underlying description of the onset of the pseudogap in the cuprates.

We demonstrate that collective vibrational strong coupling of molecules in thermal equilibrium can give rise to significant local electronic polarizations in the thermodynamic limit. We do so by first showing that the full non-relativistic Pauli-Fierz problem of an ensemble of strongly-coupled molecules in the dilute-gas limit reduces in the cavity Born-Oppenheimer approximation to a cavity-Hartree equation for the electronic structure. Consequently, each individual molecule experiences a self-consistent coupling to the dipoles of all other molecules, which amount to non-negligible values in the thermodynamic limit (large ensembles). Thus collective vibrational strong coupling can alter individual molecules strongly for localized "hotspots" within the ensemble. Moreover, the discovered cavity-induced polarization pattern possesses a zero net polarization, which resembles a continuous form of a spin glass (or better polarization glass). Our findings suggest that the thorough understanding of polaritonic chemistry, requires a self-consistent treatment of dressed electronic structure, which can give rise to numerous, so far overlooked, physical mechanisms.