Publications

Time-dependent Wannier functions were initially proposed as a means for calculating the polarization current in crystals driven by external fields. In this work, we present a simple gauge where Wannier states are defined based on the maximally localized functions at the initial time, and are propagated using the time-dependent Bloch states obtained from established first-principles calculations, avoiding the costly Wannierization at ech time step. We show that this basis efficiently describes the time-dependent polarization of the laser driven system through the analysis of the motion of Wannier centers. We use this technique to analyze highly nonlinear and non-perturbative responses such as high harmonic generation in solids, using the hexagonal boron nitride as an illustrative example, and we show how it provides an intuitive picture for the physical mechanisms.

We study the electronic structure of bulk 1T-TaSe

We present the theoretical derivation and numerical implementation of the linear response equations for relativistic quantum electrodynamical density functional theory (QEDFT). In contrast to previous works based on the Pauli-Fierz Hamiltonian, our approach describes electrons interacting with photonic cavity modes at the four-component Dirac-Kohn-Sham level, derived from fully relativistic QED through a series of established approximations. Moreover, we show that a new type of spin-orbit-like (SO) cavity-mediated interaction appears under the relativistic description of the coupling of matter with quantized cavity modes. Benchmark calculations performed for atoms of group 12 elements (Zn, Cd, Hg) demonstrate how a relativistic treatment enables the description of exciton polaritons which arise from the hybridization of formally forbidden singlet-triplet transitions with cavity modes. For atoms in cavities tuned on resonance with a singlet-triplet transition we discover a significant interplay between SO effects and coupling to an off-resonant intense singlet-singlet transition. This dynamic relationship highlights the crucial role of ab initio approaches in understanding cavity quantum electrodynamics. Finally, using the mercury porphyrin complex as an example, we show that relativistic linear response QEDFT provides computationally feasible first-principles calculations of polaritonic states in large heavy element-containing molecules of chemical interest.

Many-body systems which saturate the quantum bound on chaos are attracting interest across a wide range of fields. Notable examples include the Sachdev-Ye-Kitaev model and its variations, all characterised by some form or randomness and all to all couplings. Here we study many-body quantum chaos in a quantum impurity model showing Non-Fermi-Liquid physics, the overscreened multichannel SU(N) Kondo model. We compute exactly the low-temperature behavior of the out-of time order correlator in the limit of large N and large number of channels K, at fixed ratio γ=K/N. Due to strong correlations at the impurity site the spin fractionalizes in auxiliary fermions and bosons. We show that all the degrees of freedom of our theory acquire a Lyapunov exponent which is linear in temperature as T→0, with a prefactor that depends on γ. Remarkably, for N=K the impurity spin displays maximal chaos, while bosons and fermions only get up to half of the maximal Lyapunov exponent. Our results highlights two new features: a non-disordered model which is maximally chaotic due to strong correlations at its boundary and a fractionalization of quantum chaos.

Quantum-enhanced sensors, which surpass the standard quantum limit (SQL) and approach the fundamental precision limits dictated by quantum mechanics, are finding applications across a wide range of scientific fields. This quantum advantage becomes particularly significant when a large number of particles are included in the sensing circuit. Achieving such enhancement requires introducing and preserving entanglement among many particles, posing significant experimental challenges. In this work, we integrate concepts from Floquet theory and quantum information to design an entangler capable of generating the desired entanglement between two paths of a quantum interferometer. We demonstrate that our path-entangled states enable sensing beyond the SQL, reaching the fundamental Heisenberg limit (HL) of quantum mechanics. Moreover, we show that a decoding parity measurement maintains the HL when specific conditions from Floquet theory are satisfied–particularly those related to the periodic driving parameters that preserve entanglement during evolution. We address the effects of a priori phase uncertainty and imperfect transmission, showing that our method remains robust under realistic conditions. Finally, we propose a superconducting-circuit implementation of our sensor in the microwave regime, highlighting its potential for practical applications in high-precision measurements.

A beyond electric-dipole light-matter theory is needed to describe emerging X-ray and THz applications for characterization and control of quantum materials but inaccessible as nondipole lattice-aperiodic terms impede on the use of Bloch's theorem. To circumvent this, we derive a formalism that captures dominant nondipole effects in intense electromagnetic fields while conserving lattice translational symmetry. Our approach enables the first accurate nondipole first-principles microscopic simulation of nonperturbative harmonic generation in Si. We reveal nondipole-induced transverse currents generating perturbative even-ordered harmonics and display the onset of nondipole high harmonic generation near the laser damage threshold.

Recently, neural networks (NNs) have become a powerful tool for detecting quantum phases of matter. Unfortunately, NNs are black boxes and only identify phases without elucidating their properties. Novel physics benefits most from insights about phases, traditionally extracted in spin systems using spin correlators. Here, we combine two approaches and design TetrisCNN, a convolutional NN with parallel branches using different kernels that detects the phases of spin systems and expresses their essential descriptors, called order parameters, in a symbolic form based on spin correlators. We demonstrate this on the example of snapshots of the one-dimensional transverse-field Ising model taken in various bases. We show also that TetrisCNN can detect more complex order parameters using the example of two-dimensional Ising gauge theory. This work can lead to the integration of NNs with quantum simulators to study new exotic phases of matter.

The tensor cross interpolation (TCI) algorithm is a rank-revealing algorithm for decomposing low-rank, high-dimensional tensors into tensor trains/matrix product states (MPS). TCI learns a compact MPS representation of the entire object from a tiny training data set. Once obtained, the large existing MPS toolbox provides exponentially fast algorithms for performing a large set of operations. We discuss several improvements and variants of TCI. In particular, we show that replacing the cross interpolation by the partially rank-revealing LU decomposition yields a more stable and more flexible algorithm than the original algorithm. We also present two open source libraries, xfac in Python/C++ and this http URL in Julia, that implement these improved algorithms, and illustrate them on several applications. These include sign-problem-free integration in large dimension, the superhigh-resolution quantics representation of functions, the solution of partial differential equations, the superfast Fourier transform, the computation of partition functions, and the construction of matrix product operators.

We propose exchanging the energy functionals in ground-state DFT with physically equivalent exact force expressions as a new promising route towards approximations to the exchange-correlation potential and energy. In analogy to the usual energy-based procedure, we split the force difference between the interacting and auxiliary Kohn-Sham system into a Hartree, an exchange, and a correlation force. The corresponding scalar potential is obtained by solving a Poisson equation, while an additional transverse part of the force yields a vector potential. These vector potentials obey an exact constraint between the exchange and correlation contribution and can further be related to the atomic-shell structure. Numerically, the force-based local-exchange potential and the corresponding exchange energy compare well with the numerically more involved optimized effective-potential method. Overall, the force-based method has several benefits when compared to the usual energy-based approach and opens a route towards numerically inexpensive non-local and (in the time-dependent case) non-adiabatic approximations.

We report the first ab initio, non-relativistic QED method that couples light and matter self-consistently beyond the electric dipole approximation and without multipolar truncations. This method is based on an extension of the Maxwell-Pauli-Kohn-Sham approach to a full minimal coupling Hamiltonian, where the space- and time-dependent vector potential is coupled to the matter system, and its back-reaction to the radiated fields is generated by the full current density. The implementation in the open-source Octopus code is designed for massively-parallel multiscale simulations considering different grid spacings for the Maxwell and matter subsystems. Here, we show the first applications of this framework to simulate renormalized Cherenkov radiation of an electronic wavepacket, magnetooptical effects with non-chiral light in non-chiral molecular systems, and renormalized plasmonic modes in a nanoplasmonic dimer. We show that in some cases the beyond-dipole effects can not be captured by a multipolar expansion Hamiltonian in the length gauge. Finally, we discuss further opportunities enabled by the framework in the field of twisted light and orbital angular momentum, inelastic light scattering and strong field physics.