Publications

We present an efficient method for propagating the time-dependent Kohn–Sham equations in free space, based on the recently introduced Fourier contour deformation (FCD) approach. For potentials which are constant outside a bounded domain, FCD yields a high-order accurate numerical solution of the time-dependent Schrödinger equation directly in free space, without the need for artificial boundary conditions. Of the many existing artificial boundary condition schemes, FCD is most similar to an exact nonlocal transparent boundary condition, but it works directly on Cartesian grids in any dimension, and runs on top of the fast Fourier transform rather than fast algorithms for the application of nonlocal history integral operators. We adapt FCD to time-dependent density functional theory (TDDFT), and describe a simple algorithm to smoothly and automatically truncate long-range Coulomb-like potentials to a time-dependent constant outside of a bounded domain of interest, so that FCD can be used. This approach eliminates errors originating from the use of artificial boundary conditions, leaving only the error of the potential truncation, which is controlled and can be systematically reduced. The method enables accurate simulations of ultrastrong nonlinear electronic processes in molecular complexes in which the interference between bound and continuum states is of paramount importance. We demonstrate results for many-electron TDDFT calculations of absorption and strong field photoelectron spectra for one and two-dimensional models, and observe a significant reduction in the size of the computational domain required to achieve high quality results, as compared with the popular method of complex absorbing potentials.

The quest for primordial B-modes in the cosmic microwave background has emphasized the need for refined models of the Galactic dust foreground. Here we aim at building a realistic statistical model of the multifrequency dust emission from a single example. We introduce a generic methodology relying on microcanonical gradient descent models conditioned by an extended family of wavelet phase harmonic (WPH) statistics. To tackle the multichannel aspect of the data, we define cross-WPH statistics, quantifying non-Gaussian correlations between maps. Our data-driven methodology could apply to various contexts, and we have updated the software PyWPH, on which this work relies, accordingly. Applying this to dust emission maps built from a magnetohydrodynamics simulation, we construct and assess two generative models: (1) a (I, E, B) multi-observable input, and (2) a {I

We investigate the real-space profile of effective Coulomb interactions in correlated kagome materials. By particularizing to KV

The Sliced-Wasserstein distance (SW) is a computationally efficient and theoretically grounded alternative to the Wasserstein distance. Yet, the literature on its statistical properties -- or, more accurately, its generalization properties -- with respect to the distribution of slices, beyond the uniform measure, is scarce. To bring new contributions to this line of research, we leverage the PAC-Bayesian theory and a central observation that SW may be interpreted as an average risk, the quantity PAC-Bayesian bounds have been designed to characterize. We provide three types of results: i) PAC-Bayesian generalization bounds that hold on what we refer as adaptive Sliced-Wasserstein distances, i.e. SW defined with respect to arbitrary distributions of slices (among which data-dependent distributions), ii) a principled procedure to learn the distribution of slices that yields maximally discriminative SW, by optimizing our theoretical bounds, and iii) empirical illustrations of our theoretical findings.

We introduce what we believe to be a novel method to perform linear optical random projections without the need for holography. Our method consists of a computationally trivial combination of multiple intensity measurements to mitigate the information loss usually associated with the absolute-square non-linearity imposed by optical intensity measurements. Both experimental and numerical findings demonstrate that the resulting matrix consists of real-valued, independent, and identically distributed (i.i.d.) Gaussian random entries. Our optical setup is simple and robust, as it does not require interference between two beams. We demonstrate the practical applicability of our method by performing dimensionality reduction on high-dimensional data, a common task in randomized numerical linear algebra with relevant applications in machine learning.

The spectral and transport properties of strongly correlated metals, such as SrVO

The spectral and transport properties of strongly correlated metals, such as SrVO

Combining the complementary capabilities of two of the most powerful modern computational methods, we find superconductivity in both the electron- and hole-doped regimes of the two-dimensional Hubbard model (with next nearest neighbor hopping). In the electron-doped regime, superconductivity is weaker and is accompanied by antiferromagnetic Néel correlations at low doping. The strong superconductivity on the hole-doped side coexists with stripe order, which persists into the overdoped region with weaker hole density modulation. These stripe orders, neither filled as in the pure Hubbard model (no next nearest neighbor hopping) nor half-filled as seen in previous state-of-the-art calculations, vary in fillings between 0.6 and 0.8. The resolution of the tiny energy scales separating competing orders requires exceedingly high accuracy combined with averaging and extrapolating with a wide range of system sizes and boundary conditions. These results validate the applicability of this iconic model for describing cuprate high-T

Two-dimensional moiré materials have emerged as the most versatile platforms for realizing quantum phases of electrons. Here, we explore the stability origins of correlated states in WSe2/WS2 moiré superlattices. We find that ultrafast electronic excitation leads to melting of the Mott states on time scales five times longer than predictions from the charge hopping integrals and the melting rates are thermally activated, with activation energies of 18 and 13 meV for the one- and two-hole Mott states, respectively, suggesting significant electron-phonon coupling. DFT calculation of the one-hole Mott state confirms polaron formation and yields a hole-polaron binding energy of 16 meV. These findings reveal a close interplay of electron-electron and electron-phonon interactions in stabilizing the polaronic Mott insulators at transition metal dichalcogenide moiré interfaces.

Controlling edge states of topological magnon insulators is a promising route to stable spintronics devices. However, to experimentally ascertain the topology of magnon bands is a challenging task. Here we derive a fundamental relation between the light-matter coupling and the quantum geometry of magnon states. This allows to establish the two-magnon Raman circular dichroism as an optical probe of magnon topology in honeycomb magnets, in particular of the Chern number and the topological gap. Our results pave the way for interfacing light and topological magnons in functional quantum devices.