Source author record

Angel Rubio

Angel Rubio appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.mtrl-sci cond-mat.mes-hall physics.chem-ph cond-mat.str-el quant-ph physics.optics cond-mat.other physics.comp-ph physics.atom-ph cond-mat.supr-con physics.atm-clus physics.plasm-ph cond-mat cond-mat.quant-gas cond-mat.soft cond-mat.stat-mech

Catalog footprint

What is connected

137works

16topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

Dynamical instability in a Floquet-Driven Dissipative System

We analyse the magnon spectrum and distribution function of the antiferromagnetic phase of the Floquet-driven Hubbard model. Above a critical drive strength, we find a dynamical instability, resulting from a change in sign of the magnon damping at a non-zero wavevector. The change in sign means that infinitesimal fluctuations grow with time, corresponding to an instability of the driven state. Implications for the nonequilibrium distribution function and the strong drive nonlinear dynamics are discussed.

preprint2026arXiv

Vacuum-dressed superconductivity in NbN observed in a high-$Q$ terahertz cavity

Emerging theoretical frameworks suggest that physical properties of matter can be altered within an optical cavity by harnessing quantum vacuum electromagnetic fluctuations, even in the total absence of external driving fields. Among the most intriguing predictions is the potential to noninvasively manipulate superconductivity. Here, we experimentally observe modified superconductivity in niobium nitride (NbN) thin films within high-quality-factor ($Q$) terahertz cavities. Using terahertz time-domain spectroscopy, we characterize the NbN response both in free space and within a high-$Q$ photonic-crystal cavity. Our analysis reveals significant cavity-induced modifications to the optical conductivity. A theoretical model indicates that these changes originate from a substantial ($\sim13\,\%$) reduction in the superfluid density and a minor ($\sim2\,\%$) reduction in the superconducting gap, driven by cavity vacuum fluctuations. These results demonstrate a platform for engineering ground states via vacuum--matter coupling, opening frontiers in cavity materials science.

preprint2024arXiv

Ab initio calculations of quantum light-matter interactions in general electromagnetic environments

The emerging field of strongly coupled light-matter systems has drawn significant attention in recent years due to the prospect of altering physical and chemical properties of molecules and materials. Because this emerging field draws on ideas from both condensed-matter physics and quantum optics, it has attracted attention from theoreticians from both fields. While the former employ accurate descriptions of the electronic structure of the matter the description of the electromagnetic environment is often oversimplified. Contrastingly, the latter often employs sophisticated descriptions of the electromagnetic environment, while using simple few-level approximations for the matter. Both approaches are problematic because the oversimplified descriptions of the electronic system are incapable of describing effects such as light-induced structural changes, while the oversimplified descriptions of the electromagnetic environments can lead to unphysical predictions because the light-matter interactions strengths are misrepresented. Here we overcome these shortcomings and present the first method which can quantitatively describe both the electronic system and general electromagnetic environments from first principles. We realize this by combining macroscopic QED (MQED) with Quantum Electrodynamical Density-functional Theory. To exemplify this approach, we consider an absorbing spherical cavity and study the impact of different parameters of both the environment and the electronic system on the transition from weak-to-strong coupling for different aromatic molecules. As part of this work, we also provide an easy-to-use tool to calculate the cavity coupling strengths for simple cavity setups. Our work is a step towards parameter-free ab initio calculations for strongly coupled quantum light-matter systems and will help bridge the gap between theoretical methods and experiments in the field.

preprint2022arXiv

A perspective on ab initio modeling of polaritonic chemistry: The role of non-equilibrium effects and quantum collectivity

This perspective provides a brief introduction into the theoretical complexity of polaritonic chemistry, which emerges from the hybrid nature of strongly coupled light-matter states. To tackle this complexity, the importance of ab initio methods is highlighted. Based on those, novel ideas and research avenues are developed with respect to quantum collectivity, as well as for resonance phenomena immanent in reaction rates under vibrational strong coupling. Indeed, fundamental theoretical questions arise about the mesoscopic scale of quantum-collectively coupled molecules, when considering the depolarization shift in the interpretation of experimental data. Furthermore, to rationalise recent findings based on quantum electrodynamical density-functional theory (QEDFT), a simple, but computationally efficient, Langevin framework is proposed, based on well-established methods from molecular dynamics. It suggests the emergence of cavity induced non-equilibrium nuclear dynamics, where thermal (stochastic) resonance phenomena could emerge in the absence of external periodic driving. Overall, we believe the latest ab initio results indeed suggest a paradigmatic shift for ground-state chemical reactions under vibrational strong coupling, from the collective quantum interpretation towards a more local, (semi)-classically and non-equilibrium dominated perspective. Finally, various extensions towards a refined description of cavity-modified chemistry are introduced in the context of QEDFT and future directions of the field are sketched.

preprint2022arXiv

Ab-initio cluster approach for high-harmonic generation in liquids

High harmonic generation (HHG) takes place in all phases of matter. In gaseous atomic and molecular media, it has been extensively studied and is very well-understood. In solids research is ongoing, but a consensus is forming for the dominant microscopic HHG mechanisms. In liquids on the other hand, no established theory yet exit and approaches developed for gases and solids are generally inapplicable, hindering our current understanding. We develop here a powerful and reliable ab-initio cluster-based approach for describing the nonlinear interactions between isotropic bulk liquids and intense laser pulses. The scheme is based on time-dependent density functional theory and utilizes several approximations that make it feasible yet accurate in realistic systems. We demonstrate our approach with HHG calculations in water, ammonia, and methane liquids, and compare the characteristic response of polar and non-polar liquids. We identify unique features in the HHG spectra of liquid methane that could be utilized for ultrafast spectroscopy of its chemical and physical properties: (i) a structural minima at 15-17eV, and (ii) a well-like shape in the perturbative region that is reminiscent of a shape resonance. Our results pave the way to accessible calculations of HHG in liquids and illustrate the unique nonlinear nature of liquid systems.

preprint2022arXiv

Common microscopic origin of the phase transitions in Ta$_2$NiS$_5$ and the excitonic insulator candidate Ta$_2$NiSe$_5$

The structural phase transition in Ta$_2$NiSe$_5$ has been envisioned as driven by the formation of an excitonic insulating phase. However, the role of structural and electronic instabilities on crystal symmetry breaking has yet to be disentangled. Meanwhile, the phase transition in its complementary material Ta$_2$NiS$_5$ does not show any experimental hints of an excitonic insulating phase. We present a microscopic investigation of the electronic and phononic effects involved in the structural phase transition in Ta$_2$NiSe$_5$ and Ta$_2$NiS$_5$ using extensive first-principles calculations. In both materials the crystal symmetries are broken by phonon instabilities, which in turn lead to changes in the electronic bandstructure also observed in experiment. A total energy landscape analysis shows no tendency towards a purely electronic instability and we find that a sizeable lattice distortion is needed to open a bandgap. We conclude that an excitonic instability is not needed to explain the phase transition in both Ta$_2$NiSe$_5$ and Ta$_2$NiS$_5$.

preprint2022arXiv

Exact Solution for A Real Polaritonic System Under Vibrational Strong Coupling in Thermodynamic Equilibrium: Absence of Zero Temperature and Loss of Light-Matter Entanglement

The first exact quantum simulation of a real molecular system (HD$^+$) under strong ro-vibrational coupling to a quantized optical cavity mode in thermal equilibrium is presented. Theoretical challenges in describing strongly coupled systems of mixed quantum statistics (Bosons and Fermions) are discussed and circumvented by the specific choice of our molecular system. Our exact simulations reveal the absence of a zero temperature for the strongly coupled matter and light subsystems, due to cavity induced non-equilibrium conditions. Furthermore, we explore the temperature dependency of light-matter quantum entanglement, which emerges for the groundstate, but is quickly lost already in the deep cryogenic regime, opposing predictions from phenomenological models (Jaynes-Cummings). Distillable molecular light-matter entanglement of ro-vibrational states may open interesting perspectives for quantum technological applications. Moreover, we find that the dynamics (fluctuations) of matter remains modified by the quantum nature of the thermal and vacuum field fluctuations for significant temperatures, e.g. at ambient conditions. These observations (loss of entanglement and coupling to quantum fluctuations) has far reaching consequences for the understanding and control of polaritonic chemistry and materials science, since a semi-classical theoretical description of light-matter interaction becomes feasible, but the typical canonical equilibrium assumption for the nuclear dynamics remains broken. This opens the door for quantum fluctuations induced stochastic resonance phenomena under vibrational strong coupling. A plausible theoretical mechanism to explain the experimentally observed resonance phenomena in absence of periodic driving, which have not yet been understood.

preprint2022arXiv

Floquet engineering the band structure of materials with optimal control theory

We demonstrate that the electronic structure of a material can be deformed into Floquet pseudo-bands with arbitrarily tailored shapes. We achieve this goal with a novel combination of quantum optimal control theory and Floquet engineering. The power and versatility of this framework is demonstrated here by utilizing the independent-electron tight-binding description of the $π$ electronic system of graphene. We show several prototype examples focusing on the region around the K (Dirac) point of the Brillouin zone: creation of a gap with opposing flat valence and conduction bands, creation of a gap with opposing concave symmetric valence and conduction bands -- which would correspond to a material with an effective negative electron-hole mass --, or closure of the gap when departing from a modified graphene model with a non-zero field-free gap. We employ time periodic drives with several frequency components and polarizations, in contrast to the usual monochromatic fields, and use control theory to find the amplitudes of each component that optimize the shape of the bands as desired. In addition, we use quantum control methods to find realistic switch-on pulses that bring the material into the predefined stationary Floquet band structure, i.e. into a state in which the desired Floquet modes of the target bands are fully occupied, so that they should remain stroboscopically stationary, with long lifetimes, when the weak periodic drives are started. Finally, we note that although we have focused on solid state materials, the technique that we propose could be equally used for the Floquet engineering of ultracold atoms in optical lattices, and to other non-equilibrium dynamical and correlated systems.

preprint2022arXiv

Frequency-dependent Sternheimer linear-response formalism for strongly coupled light-matter systems

The rapid progress in quantum-optical experiments especially in the field of cavity quantum electrodynamics and nanoplasmonics, allows to substantially modify and control chemical and physical properties of atoms, molecules and solids by strongly coupling to the quantized field. Alongside such experimental advances has been the recent development of ab-initio approaches such as quantum electrodynamical density-functional theory (QEDFT) that is capable of describing these strongly coupled systems from first-principles. To investigate response properties of relatively large systems coupled to a wide range of photon modes, ab-initio methods that scale well with system size become relevant. In light of this, we extend the linear-response Sternheimer approach within the framework of QEDFT to efficiently compute excited-state properties of strongly coupled light-matter systems. Using this method, we capture features of strong light-matter coupling both in the dispersion and absorption properties of a molecular system strongly coupled to the modes of a cavity. We exemplify the efficiency of the Sternheimer approach by coupling the matter system to the continuum of an electromagnetic field. We observe changes in the spectral features of the coupled system as Lorentzian line shapes turn into Fano resonances when the molecule interacts strongly with the continuum of modes. This work provides an alternative approach for computing efficiently excited-state properties of large molecular systems interacting with the quantized electromagnetic field.

preprint2022arXiv

Moiré Engineering of Nonsymmorphic Symmetries and Hourglass Superconductors

Moiré heterostructures hold the promise to provide platforms to tailor strongly correlated and topological states of matter. Here, we theoretically propose the emergence of an effective, rectangular moiré lattice in twisted bilayers of SnS with nonsymmorphic symmetry. Based on first-principles calculations, we demonstrate that strong intrinsic spin-orbit interactions render this tunable platform a moiré semimetal that hosts 2D hourglass fermions protected by time-reversal symmetry $\mathcal{T}$ and the nonsymmorphic screw rotation symmetry $\widetilde{\mathcal{C}}_{2y}$. We show that topological Fermi arcs connecting pairs of Weyl nodal points in the hourglass dispersion are preserved for weak electron-electron interactions, particularly in regions of superconducting order that emerge in the phase diagram of interaction strength and filling. Our work established moiré engineering as an inroad into the realm of correlated topological semimetals and may motivate further topology related researches in moiré heterostructures.

preprint2022arXiv

Moiré Engineering of Spin-Orbit Coupling in Twisted Platinum Diselenide

We study the electronic structure and correlated phases of twisted bilayers of platinum diselenide using large-scale ab initio simulations combined with the functional renormalization group. PtSe$_2$ is a group-X transition metal dichalcogenide, which hosts emergent flat bands at small twist angles in the twisted bilayer. Remarkably, we find that moiré engineering can be used to tune the strength of Rashba spin-orbit interactions, altering the electronic behavior in a novel manner. We reveal that an effective triangular lattice with a twist-controlled ratio between kinetic and spin-orbit coupling scales can be realized. Even dominant spin-orbit coupling can be accessed in this way and we discuss consequences for the interaction driven phase diagram, which features pronounced exotic superconducting and entangled spin-charge density waves.

preprint2022arXiv

New Class of Landau Levels and Hall Phases in a 2D Electron Gas Subject to an Inhomogeneous Magnetic Field: An Analytic Solution

An analytic closed form solution is derived for the bound states of electrons subject to a static, inhomogeneous ($1/r$-decaying) magnetic field, including the Zeeman interaction. The solution provides access to many-body properties of a two-dimensional, non-interacting, electron gas in the thermodynamic limit. Radially distorted Landau levels can be identified as well as magnetic field induced density and current oscillations close to the magnetic impurity. These radially localised oscillations depend strongly on the coupling of the spin to the magnetic field, which give rise to non-trivial spin currents. Moreover, the Zeeman interaction introduces a lowest flat band for $E_F=0^+$ assuming a spin $g_s$-factor of two. Surprisingly, in this case the charge and current densities can be computed analytically in the thermodynamic limit. Numerical calculations show that the total magnetic response of the electron gas remains diamagnetic (similar to Landau levels) independent of the Fermi energy. However, the contribution of certain, infinitely degenerate energy levels may become paramagnetic. Furthermore, numerical computations of the Hall conductivity reveal asymptotic properties of the electron gas, which are driven by the anisotropy of the vector potential instead of the magnetic field, i.e. become independent of spin. Eventually, the distorted Landau levels give rise to different Hall conductivity phases, which are characterized by sharp sign flips at specific Fermi energies. Overall, our work merges "impurity" with Landau-level physics, which provides novel physical insights, not only locally, but also in the asymptotic limit. This paves the way for a large number of future theoretical as well as experimental investigations.

preprint2022arXiv

Polaritonic Hofstadter Butterfly and Cavity-Control of the Quantized Hall Conductance

In a previous work [Phys. Rev. Lett. 123, 047202 (2019)] a translationally invariant framework called quantum-electrodynamical Bloch (QED-Bloch) theory was introduced for the description of periodic materials in homogeneous magnetic fields and strongly coupled to the quantized photon field in the optical limit. For such systems, we show that QED-Bloch theory predicts the existence of fractal polaritonic spectra as a function of the cavity coupling strength. In addition, for the energy spectrum as a function of the relative magnetic flux we find that a terahertz cavity can modify the standard Hofstadter butterfly. In the limit of no quantized photon field, QED-Bloch theory captures the well-known fractal spectrum of the Hofstadter butterfly and can be used for the description of 2D materials in strong magnetic fields, which are of great experimental interest. As a further application, we consider Landau levels under cavity confinement and show that the cavity alters the quantized Hall conductance and that the Hall plateaus are modified as $σ_{xy}=e^2ν/h(1+η^2)$ by the light-matter coupling $η$. Most of the aforementioned phenomena should be experimentally accessible and corresponding implications are discussed.

preprint2022arXiv

Probing phonon dynamics with multi-dimensional high harmonic carrier envelope phase spectroscopy

We explore pump-probe high harmonic generation (HHG) from monolayer hexagonal-Boron-Nitride, where a terahertz pump excites coherent optical phonons that are subsequently probed by an intense infrared pulse that drives HHG. We find, through state-of-the-art ab-initio calculations, that the structure of the emission spectrum is attenuated by the presence of coherent phonons, and is no longer comprised of discrete harmonic orders, but rather of a continuous emission in the plateau region. The HHG yield strongly oscillates as a function of the pump-probe delay, corresponding to ultrafast changes in the lattice such as bond compression or stretching. We further show that in the regime where the excited phonon period and the pulse duration are of the same order of magnitude, the HHG process becomes sensitive to the carrier-envelope-phase (CEP) of the driving field, even though the pulse duration is so long that no such sensitivity is observed in the absence of coherent phonons. The degree of CEP sensitivity vs. pump-probe delay is shown to be a highly selective measure for instantaneous structural changes in the lattice, providing a new approach for ultrafast multi-dimensional HHG-spectroscopy. Remarkably, the obtained temporal resolution for phonon dynamics is ~1 femtosecond, which is much shorter than the probe pulse duration because of the inherent sub-cycle contrast mechanism. Our work paves the way towards novel routes of probing phonons and ultrafast material structural changes and provides a mechanism for controlling the high harmonic response.

preprint2022arXiv

Quadrupolar excitons in a tunnel-coupled van der Waals heterotrilayer

Strongly bound excitons and many-body interactions between them determine light-matter interactions in van der Waals (vdW) heterostructures of 2D semiconductors. Unlike fundamental particles, quasiparticles in condensed matter, such as excitons, can be tailored to alter their interactions and realize emergent quantum phases. Here, using a WS$_2$/WSe$_2$/WS$_2$ heterotrilayer, we create a quantum superposition of oppositely oriented dipolar excitons - a quadrupolar exciton - wherein an electron is layer-hybridized in WS$_2$ layers while the hole localizes in WSe$_2$. In contrast to dipolar excitons, symmetric quadrupolar excitons only redshift in an out-of-plane electric field, consistent with ab initio calculations, regaining dipolar characteristics at higher fields. Electric field tunes the hybridization and allows for lifetime control through modification of the excitonic wavefunction. Lack of density-dependent blue shift of heterotrilayer excitons compared to dipolar excitons is consistent with quadrupolar interactions. Our results present vdW heterotrilayers as a field-tunable platform to engineer light-matter interactions and explore quantum phase transitions between spontaneously ordered many-exciton phases.

preprint2022arXiv

Soft-x-ray spectroscopy of monolayer transition metal dichalcogenides: valley selectivity?

Optical properties of semiconducting monolayer transition metal dichalcogenides have received a lot of attention in recent years, following the discovery of the valley selective optical population of either $K_{+}$ or $ K_{-}$ valleys at the direct band gap, depending on the polarisation of the incoming light. We use group theoretical selection rules, as well as ab initio DFT calculations, to investigate whether this valence/conduction valley selectivity effect is also present in x-ray optical transitions from a core level of the metal atom to the valence and conduction band $K$ valleys. Valley selectivity is predicted for $s$, $p_{1/2}$ and $p_{3/2}$ edges in transitions to and from the valence band edges with circularly polarised radiation. Possible novel applications to the diagnostics of valleytronic properties and intervalley dynamics are investigated and the feasibility of time-resolved pump-probe and high resolution experiments with suitable soft-x-ray free-electron laser sources is discussed.

preprint2022arXiv

Spin-orbit induced equilibrium spin currents in materials

The existence of spin-currents in absence of any driving external fields is commonly considered an exotic phenomenon appearing only in quantum materials, such as topological insulators. We demonstrate instead that equilibrium spin currents are a rather general property of materials with non negligible spin-orbit coupling (SOC). Equilibrium spin currents can be present at the surfaces of a slab. Yet, we also propose the existence of global equilibrium spin currents, which are net bulk spin-currents along specific crystallographic directions of materials. Equilibrium spin currents are allowed by symmetry in a very broad class of systems having gyrotropic point groups. The physics behind equilibrium spin currents is uncovered by making an analogy between electronic systems with SOC and non-Abelian gauge theories. The electron spin can be seen as the analogous of the color degree of freedom and equilibrium spin currents can then be identified with diamagnetic color currents appearing as the response to an effective non-Abelian magnetic field generated by SOC. Equilibrium spin currents are not associated with spin transport and accumulation, but they should nonetheless be carefully taken into account when computing transport spin currents. We provide quantitative estimates of equilibrium spin currents for several systems, specifically metallic surfaces presenting Rashba-like surface states, nitride semiconducting nanostructures and bulk materials, such as the prototypical gyrotropic medium tellurium. In doing so, we also point out the limitations of model approaches showing that first-principles calculations are needed to obtain reliable predictions. We therefore use Density Functional Theory computing the so-called bond currents, which represent a powerful tool to understand the relation between equilibrium currents, electronic structure and crystal point group.

preprint2022arXiv

Terahertz-induced high-order harmonic generation and nonlinear charge transport in graphene

We theoretically study the THz-induced high-order harmonic generation (HHG) and nonlinear electric transport in graphene based on the quantum master equation with the relaxation time approximation. To obtain microscopic insight into the phenomena, we compare the results of the fully dynamical calculations with those under a quasi-static approximation, where the electronic system is approximated as a nonequilibrium steady state. As a result, we find that the THz-induced electron dynamics in graphene can be accurately modeled with the nonequilibrium steady-state at each instance. The population distribution analysis further clarifies that the THz-induced HHG in graphene originates from the reduction of effective conductivity due to a large displacement of electrons in the Brillouin zone. By comparing the present nonequilibrium picture with a thermodynamic picture, we explore the role of the nonequilibrium nature of electron dynamics on the extremely nonlinear optical and transport phenomena in graphene.

preprint2022arXiv

Time- and angle-resolved photoelectron spectroscopy of strong-field light-dressed solids: prevalence of the adiabatic band picture

In recent years, strong-field physics in condensed-matter was pioneered as a novel approach for controlling material properties through laser-dressing, as well as for ultrafast spectroscopy via nonlinear light-matter interactions (e.g. harmonic generation). A potential controversy arising from these advancements is that it is sometimes vague which band-picture should be used to interpret strong-field experiments: the field-free bands, the adiabatic (instantaneous) field-dressed bands, Floquet bands, or some other intermediate picture. We here try to resolve this issue by performing 'theoretical experiments' of time- and angle-resolved photoelectron spectroscopy (Tr-ARPES) for a strong-field laser-pumped solid, which should give access to the actual observable bands of the irradiated material. To our surprise, we find that the adiabatic band-picture survives quite well, up to high field intensities (~10^12 W/cm^2), and in a wide frequency range (driving wavelengths of 4000 to 800nm, with Keldysh parameters ranging up to ~7). We conclude that to first order, the adiabatic instantaneous bands should be the standard blueprint for interpreting ultrafast electron dynamics in solids when the field is highly off-resonant with characteristic energy scales of the material. We then discuss weaker effects of modifications of the bands beyond this picture that are non-adiabatic, showing that by using bi-chromatic fields the deviations from the standard picture can be probed with enhanced sensitivity. Our work outlines a clear band picture for the physics of strong-field interactions in solids, which should be useful for designing and analyzing strong-field experimental observables and also to formulate simpler semi-empirical models.

preprint2022arXiv

Unconventional Excitonic States with Phonon Sidebands in Layered Silicon Diphosphide

Many-body interactions between quasiparticles (electrons, excitons, and phonons) have led to the emergence of new complex correlated states and are at the core of condensed matter physics and material science. In low-dimensional materials, unique electronic properties for these correlated states could significantly affect their optical properties. Herein, combining photoluminescence, optical reflection measurements and theoretical calculations, we demonstrate an unconventional excitonic state and its bound phonon sideband in layered silicon diphosphide (SiP$_2$), in which the bound electron-hole pair is composed of electrons confined within one-dimensional phosphorus$-$phosphorus chains and holes extended in two-dimensional SiP$_2$ layers. The excitonic state and the emergent phonon sideband show linear dichroism and large energy redshifts with increasing temperature. Within the $GW$ plus Bethe$-$Salpeter equation calculations and solving the generalized Holstein model non-perturbatively, we confirm that the observed sideband feature results from the correlated interaction between excitons and optical phonons. Such a layered material provides a new platform to study excitonic physics and many-particle effects.

preprint2022arXiv

Wavefunction Embedding for Molecular Polaritons

Polaritonic chemistry relies on the strong light-matter interaction phenomena for altering the chemical reaction rates inside optical cavities. To explain and to understand these processes, the development of reliable theoretical models is essential. While computationally efficient quantum electrodynamics self-consistent field (QED-SCF) methods, such as quantum electrodynamics density functional theory (QEDFT) needs accurate functionals, quantum electrodynamics coupled cluster (QED-CC) methods provide a systematic increase in accuracy but at much greater cost. To overcome this computational bottleneck, herein we introduce and develop the QED-CC-in-QED-SCF projection-based embedding method that inherits all the favorable properties from the two worlds, computational efficiency and accuracy. The performance of the embedding method is assessed by studying some prototypical but relevant reactions, such as methyl transfer reaction, proton transfer reaction, as well as protonation reaction in a complex environment. The results obtained with the new embedding method are in excellent agreement with more expensive QED-CC results. The analysis performed on these reactions indicate that the strong light-matter interaction is very local in nature and that only a small region should be treated at the QED-CC level for capturing important effects due to cavity. This work sets the stage for future developments of polaritonic quantum chemistry methods and it will serve as a guideline for development of other polaritonic embedding models.

preprint2021arXiv

Approximations based on density-matrix embedding theory for density-functional theories

Recently a novel approach to find approximate exchange-correlation functionals in density-functional theory (DFT) was presented (U. Mordovina et. al., JCTC 15, 5209 (2019)), which relies on approximations to the interacting wave function using density-matrix embedding theory (DMET). This approximate interacting wave function is constructed by using a projection determined by an iterative procedure that makes parts of the reduced density matrix of an auxiliary system the same as the approximate interacting density matrix. If only the diagonal of both systems are connected this leads to an approximation of the interacting-to-non-interacting mapping of the Kohn-Sham approach to DFT. Yet other choices are possible and allow to connect DMET with other DFTs such as kinetic-energy DFT or reduced density-matrix functional theory. In this work we give a detailed review of the basics of the DMET procedure from a DFT perspective and show how both approaches can be used to supplement each other. We do so explicitly for the case of a one-dimensional lattice system, as this is the simplest setting where we can apply DMET and the one that was originally presented. Among others we highlight how the mappings of DFTs can be used to identify uniquely defined auxiliary systems and auxiliary projections in DMET and how to construct approximations for different DFTs using DMET inspired projections. Such alternative approximation strategies become especially important for DFTs that are based on non-linearly coupled observables such as kinetic-energy DFT, where the Kohn-Sham fields are no longer simply obtainable by functional differentiation of an energy expression, or for reduced density-matrix functional theories, where a straightforward Kohn-Sham construction is not feasible.

preprint2021arXiv

Correlation-driven sub-3 fs charge migration in ionised adenine

Sudden ionisation of a relatively large molecule can initiate a correlation-driven process dubbed charge migration, where the electron density distribution is expected to rapidly change. Capturing this few-femtosecond/attosecond charge redistribution represents the real-time observation of the electron correlation in the molecule. So far, there has been no experimental evidence of this process. Here we report on a time-resolved study of the correlation-driven charge migration process occurring in the bio-relevant molecule adenine after ionisation by a 15-35 eV attosecond pulse . We find that, the production of intact doubly charged adenine - via a shortly-delayed laser-induced second ionisation event - represents the signature of a charge inflation mechanism resulting from the many-body excitation. This conclusion is supported by first-principles time-dependent simulations. Our findings opens new important perspectives for the control of the molecular reactivity at the electronic timescale.

preprint2021arXiv

Detecting multiple chirality centers in chiral molecules with high harmonic generation

Characterizing chirality is highly important for applications in the pharmaceutical industry, as well as in the study of dynamical chemical and biological systems. However, this task has remained challenging, especially due to the ongoing increasing complexity and size of the molecular structure of drugs and active compounds. In particular, large molecules with many active chirality centers are today ubiquitous, but remain difficult to structurally analyze due to their high number of stereoisomers. Here we theoretically explore the sensitivity of high harmonic generation (HHG) to the chirality of molecules with a varying number of active chiral centers. We find that HHG driven by bi-chromatic non-collinear lasers is a sensitive probe for the stereo-configuration of a chiral molecule. We first show through calculations (from benchmark chiral molecules with up to three chirality centers) that the HHG spectrum is imprinted with information about the handedness of each chirality center in the driven molecule. Next, we show that using both classical- and deep-learning-based reconstruction algorithms, the composition of an unknown mixture of stereoisomers can be reconstructed with high fidelity by a single-shot HHG measurement. Our work illustrates how the combination of non-linear optics and machine learning might open routes for ultra-sensitive sensing in chiral systems.

preprint2021arXiv

Moiré heterostructures as a condensed matter quantum simulator

Twisted van der Waals heterostructures have latterly received prominent attention for their many remarkable experimental properties, and the promise that they hold for realising elusive states of matter in the laboratory. We propose that these systems can, in fact, be used as a robust quantum simulation platform that enables the study of strongly correlated physics and topology in quantum materials. Among the features that make these materials a versatile toolbox are the tunability of their properties through readily accessible external parameters such as gating, straining, packing and twist angle; the feasibility to realize and control a large number of fundamental many-body quantum models relevant in the field of condensed-matter physics; and finally, the availability of experimental readout protocols that directly map their rich phase diagrams in and out of equilibrium. This general framework makes it possible to robustly realize and functionalize new phases of matter in a modular fashion, thus broadening the landscape of accessible physics and holding promise for future technological applications.

preprint2021arXiv

Nanometer-scale lateral p-n junctions in graphene/$α$-RuCl$_3$ heterostructures

The ability to create high-quality lateral p-n junctions at nanometer length scales is essential for the next generation of two-dimensional (2D) electronic and plasmonic devices. Using a charge-transfer heterostructure consisting of graphene on $α$-RuCl$_3$, we conduct a proof-of-concept study demonstrating the existence of intrinsic nanoscale lateral p-n junctions in the vicinity of graphene nanobubbles. Our multi-pronged experimental approach incorporates scanning tunneling microscopy (STM) and spectroscopy (STS) and scattering-type scanning near-field optical microscopy ($\textit{s}$-SNOM) in order to simultaneously probe both the electronic and optical responses of nanobubble p-n junctions. Our STM and STS results reveal that p-n junctions with a band offset of more than 0.6 eV can be achieved over lateral length scale of less than 3 nm, giving rise to a staggering effective in-plane field in excess of 10$^8$ V/m. Concurrent $\textit{s}$-SNOM measurements confirm the utility of these nano-junctions in plasmonically-active media, and validate the use of a point-scatterer formalism for modeling surface plasmon polaritons (SPPs). Model $\textit{ab initio}$ density functional theory (DFT) calculations corroborate our experimental data and reveal a combination of sub-angstrom and few-angstrom decay processes dictating the dependence of charge transfer on layer separation. Our study provides experimental and conceptual foundations for the use of charge-transfer interfaces such as graphene/$α$-RuCl$_3$ to generate p-n nano-junctions.

preprint2021arXiv

Photoionization and transient Wannier-Stark ladder in silicon: First principle simulations versus Keldysh theory

Nonlinear photoionization of dielectrics and semiconductors is widely treated in the frames of the Keldysh theory whose validity is limited to small photon energies compared to the band gap and relatively low laser intensities. The time-dependent density functional theory (TDDFT) simulations, which are free of these limitations, enable to gain insight into non-equilibrium dynamics of the electronic structure. Here we apply the TDDFT to investigate photoionization of silicon crystal by ultrashort laser pulses in a wide range of laser wavelengths and intensities and compare the results with predictions of the Keldysh theory. Photoionization rates derived from the simulations considerably exceed the data obtained with the Keldysh theory within the validity range of the latter. Possible reasons of the discrepancy are discussed and we provide fundamental data on the photoionization rates beyond the limits of the Keldysh theory. By investigating the features of the Stark shift as a function of photon energy and laser field strength, a manifestation of the transient Wannier-Stark ladder states have been revealed which become blurred with increasing laser field strength. Finally, it is shown that the TDDFT simulations can potentially provide reliable data on the electron damping time that is of high importance for large-scale modeling.

preprint2021arXiv

The Free Electron Gas in Cavity Quantum Electrodynamics

Cavity modification of material properties and phenomena is a novel research field largely motivated by the advances in strong light-matter interactions. Despite this progress, exact solutions for extended systems strongly coupled to the photon field are not available, and both theory and experiments rely mainly on finite-system models. Therefore a paradigmatic example of an exactly solvable extended system in a cavity becomes highly desireable. To fill this gap we revisit Sommerfeld's theory of the free electron gas in cavity quantum electrodynamics (QED). We solve this system analytically in the long-wavelength limit for an arbitrary number of non-interacting electrons, and we demonstrate that the electron-photon ground state is a Fermi liquid which contains virtual photons. In contrast to models of finite systems, no ground state exists if the diamagentic $\textbf{A}^2$ term is omitted. Further, by performing linear response we show that the cavity field induces plasmon-polariton excitations and modifies the optical and the DC conductivity of the electron gas. Our exact solution allows us to consider the thermodynamic limit for both electrons and photons by constructing an effective quantum field theory. The continuum of modes leads to a many-body renormalization of the electron mass, which modifies the fermionic quasiparticle excitations of the Fermi liquid and the Wigner-Seitz radius of the interacting electron gas. Lastly, we show how the matter-modified photon field leads to a repulsive Casimir force and how the continuum of modes introduces dissipation into the light-matter system. Several of the presented findings should be experimentally accessible.

preprint2021arXiv

TMDs as a platform for spin liquid physics: A strong coupling study of twisted bilayer WSe$_2$

The advent of twisted moiré heterostructures as a playground for strongly correlated electron physics has led to a plethora of experimental and theoretical efforts seeking to unravel the nature of the emergent superconducting and insulating states. Amongst these layered compositions of two dimensional materials, transition metal dichalcogenides (TMDs) are by now appreciated as highly-tunable platforms to simulate reinforced electronic interactions in the presence of low-energy bands with almost negligible bandwidth. Here, we focus on the twisted homobilayer WSe$_2$ and the insulating phase at half-filling of the flat bands reported therein. More specifically, we explore the possibility of realizing quantum spin liquid (QSL) physics on the basis of a strong coupling description, including up to second nearest neighbor Heisenberg couplings $J_1$ and $J_2$, as well as Dzyaloshinskii-Moriya (DM) interactions. Mapping out the global phase diagram as a function of an out-of-plane displacement field, we indeed find evidence for putative QSL states, albeit only close to SU$(2)$ symmetric points. In the presence of finite DM couplings and XXZ anisotropy, long-range order is predominantly present, with a mix of both commensurate and incommensurate magnetic phases.

preprint2021arXiv

Ultra-Strong Spin-Orbit Coupling and Topological Moiré Engineering in Twisted ZrS2 Bilayers

We predict that twisted bilayers of 1T-ZrS$_2$ realize a novel and tunable platform to engineer two-dimensional topological quantum phases dominated by strong spin-orbit interactions. At small twist angles, ZrS$_2$ heterostructures give rise to an emergent and twist-controlled moiré Kagomé lattice, combining geometric frustration and strong spin-orbit coupling to give rise to a moiré quantum spin Hall insulator with highly controllable and nearly-dispersionless bands. We devise a generic pseudo-spin theory for group-IV transition metal dichalcogenides that relies on the two-component character of the valence band maximum of the 1T structure at $Γ$, and study the emergence of a robust quantum anomalous Hall phase as well as possible fractional Chern insulating states from strong Coulomb repulsion at fractional fillings of the topological moiré Kagomé bands. Our results establish group-IV transition metal dichalcogenide bilayers as a novel moiré platform to realize strongly-correlated topological phases in a twist-tunable setting.

preprint2021arXiv

Ultrafast Dynamical Lifshitz Transition

Fermi surface is at the heart of our understanding of metals and strongly correlated many-body systems. An abrupt change in the Fermi surface topology, also called Lifshitz transition, can lead to the emergence of fascinating phenomena like colossal magnetoresistance and superconductivity. While Lifshitz transitions have been demonstrated for a broad range of materials by equilibrium tuning of macroscopic parameters such as strain, doping, pressure and temperature, a non-equilibrium dynamical route toward ultrafast modification of the Fermi surface topology has not been experimentally demonstrated. Combining time-resolved multidimensional photoemission spectroscopy with state-of-the-art TDDFT+$U$ simulations, we introduce a novel scheme for driving an ultrafast Lifshitz transition in the correlated type-II Weyl semimetal T$\mathrm{_{d}}$-MoTe$_{2}$. We demonstrate that this non-equilibrium topological electronic transition finds its microscopic origin in the dynamical modification of the effective electronic correlations. These results shed light on a novel ultrafast scheme for controlling the Fermi surface topology in correlated quantum materials.

preprint2020arXiv

Chemistry in Quantum Cavities: Exact Results, the Impact of Thermal Velocities and Modified Dissociation

In recent years tremendous progress in the field of light-matter interactions has unveiled that strong coupling to the modes of an optical cavity can alter chemistry even at room temperature. Despite these impressive advances, many fundamental questions of chemistry in cavities remain unanswered. This is also due to a lack of exact results that can be used to validate and benchmark approximate approaches. In this work we provide such reference calculations from exact diagonalisation of the Pauli-Fierz Hamiltonian in the long-wavelength limit with an effective cavity mode. This allows us to investigate the reliability of the ubiquitous Jaynes-Cummings model not only for electronic but also for the case of ro-vibrational transitions. We demonstrate how the commonly ignored thermal velocity of charged molecular systems can influence chemical properties, while leaving the spectra invariant. Furthermore, we show the emergence of new bound polaritonic states beyond the dissociation energy limit.

preprint2020arXiv

Effect of Many Modes on Self-Polarization and Photochemical Suppression in Cavities

The standard description of cavity-modified molecular reactions typically involves a single (resonant) mode, while in reality the quantum cavity supports a range of photon modes. Here we demonstrate that as more photon modes are accounted for, physico-chemical phenomena can dramatically change, as illustrated by the cavity-induced suppression of the important and ubiquitous process of proton-coupled electron-transfer. Using a multi-trajectory Ehrenfest treatment for the photon-modes, we find that self-polarization effects become essential, and we introduce the concept of self-polarization-modified Born-Oppenheimer surfaces as a new construct to analyze dynamics. As the number of cavity photon modes increases, the increasing deviation of these surfaces from the cavity-free Born-Oppenheimer surfaces, together with the interplay between photon emission and absorption inside the widening bands of these surfaces, leads to enhanced suppression. The present findings are general and will have implications for the description and control of cavity-driven physical processes of molecules, nanostructures and solids embedded in cavities.

preprint2020arXiv

Electron-Phonon-Driven Three-Dimensional Metallicity in an Insulating Cuprate

The role of the crystal lattice for the electronic properties of cuprates and other high-temperature superconductors remains controversial despite decades of theoretical and experimental efforts. While the paradigm of strong electronic correlations suggests a purely electronic mechanism behind the insulator-to-metal transition, recently the mutual enhancement of the electron-electron and the electron-phonon interaction and its relevance to the formation of the ordered phases have also been emphasized. Here, we combine polarization-resolved ultrafast optical spectroscopy and state-of-the-art dynamical mean-field theory to show the importance of the crystal lattice in the breakdown of the correlated insulating state in an archetypal undoped cuprate. We identify signatures of electron-phonon coupling to specific fully-symmetric optical modes during the build-up of a three-dimensional metallic state that follows charge photodoping. Calculations for coherently displaced crystal structures along the relevant phonon coordinates indicate that the insulating state is remarkably unstable toward metallization despite the seemingly large charge-transfer energy scale. This hitherto unobserved insulator-to-metal transition mediated by fully-symmetric lattice modes can find extensive application in a plethora of correlated solids.

preprint2020arXiv

Giant exciton Mott density in anatase TiO2

Elucidating the carrier density at which strongly bound excitons dissociate into a plasma of uncorrelated electron-hole pairs is a central topic in the many-body physics of semiconductors. However, there is a lack of information on the high-density response of excitons absorbing in the near-to-mid ultraviolet, due to the absence of suitable experimental probes in this elusive spectral range. Here, we present a unique combination of many-body perturbation theory and state-of-the-art ultrafast broadband ultraviolet spectroscopy to unveil the interplay between the ultraviolet-absorbing two-dimensional excitons of anatase TiO$_2$ and a sea of electron-hole pairs. We discover that the critical density for the exciton Mott transition in this material is the highest ever reported in semiconductors. These results deepen our knowledge of the exciton Mott transition and pave the route toward the investigation of the exciton phase diagram in a variety of wide-gap insulators.

preprint2020arXiv

High-Harmonic generation from spin-polarised defects in solids

Generation of high-order harmonics in gases enabled to probe the attosecond electron dynamics in atoms and molecules with unprecedented resolution. Extending the techniques developed originally for atomic and molecular gases to solid state materials requires a fundamental understanding of the physics that has been partially addressed theoretically. Here we employ time-dependent density-functional theory to investigate how the electron dynamics resulting in high-harmonic emission in monolayer hexagonal boron nitride is affected by the presence of vacancies. We show how these realistic spin-polarised defects modify the harmonic emission, and demonstrate that important differences exist between harmonics from a pristine solid and a defected-solid. In particular, we found that the different spin channels are affected differently because of the presence of the spin-polarized point defect, and that localisation of the wavefunction, the geometry of the defect and the electron-electron interaction are all important ingredients to describe high-harmonic generation in defected-solids. We show that different vacancies lead to qualitatively different effects, thus opening the door to the high-harmonic imaging of spin-polarised defects in solids.

preprint2020arXiv

High-harmonics and isolated attosecond pulses from MgO

On the basis of real-time ab initio calculations, we study the non-perturbative interaction of two-color laser pulses with MgO crystal in the strong field regime to generate isolated attosecond pulse from high-harmonic emissions from MgO crystal. In this regard, we examine the impact of incident pulse characteristics such as its shape, intensity, and ellipticity as well as the consequence of the crystal anisotropy on the emitted harmonics and their corresponding isolated attosecond pulses. Our calculations predict the creation of isolated attosecond pulses with a duration of ~ 300 attoseconds; in addition, using elliptical driving pulses, the generation of elliptical isolated attosecond pulses is shown. Our work prepares the path for all solid-state compact optical devices offering perspectives beyond traditional isolated attosecond pulse emitted from atoms.

preprint2020arXiv

High-order harmonic generation in graphene: nonlinear coupling of intra and interband transitions

We investigate high-order harmonic generation (HHG) in graphene with a quantum master equation approach. The simulations reproduce the observed enhancement in HHG in graphene under elliptically polarized light [N. Yoshikawa et al, Science 356, 736 (2017)]. On the basis of a microscopic decomposition of the emitted high-order harmonics, we find that the enhancement in HHG originates from an intricate nonlinear coupling between the intraband and interband transitions that are respectively induced by perpendicular electric field components of the elliptically polarized light. Furthermore, we reveal that contributions from different excitation channels destructively interfere with each other. This finding suggests a path to potentially enhance the HHG by blocking a part of the channels and canceling the destructive interference through band-gap or chemical potential manipulation.

preprint2020arXiv

Light-matter hybrid-orbital-based first-principles methods: the influence of the polariton statistics

A detailed understanding of strong matter-photon interactions requires first-principle methods that can solve the fundamental Pauli-Fierz Hamiltonian of non-relativistic quantum electrodynamics efficiently. A possible way to extend well-established electronic-structure methods to this situation is to embed the Pauli-Fierz Hamiltonian in a higher-dimensional light-matter hybrid auxiliary configuration space. In this work we show the importance of the resulting hybrid Fermi-Bose statistics of the polaritons, which are the new fundamental particles of the ``photon-dressed'' Pauli-Fierz Hamiltonian for systems in cavities. We show that violations of these statistics can lead to unphysical results. We present an efficient way to ensure the proper symmetry of the underlying wave functions by enforcing representability conditions on the dressed one-body reduced density matrix. We further present a general prescription how to extend a given first-principles approach to polaritons and as an example introduce polaritonic Hartree-Fock theory. While being a single-reference method in polariton space, polaritonic Hartree-Fock is a multi-reference method in the electronic space, i.e. it describes electronic correlations. We also discuss possible applications to polaritonic QEDFT. We apply this theory to a lattice model and find that the more delocalized the bound-state wave function of the particles is, the stronger it reacts to photons. The main reason is that within a small energy range many states with different electronic configurations are available as opposed to a strongly bound (and hence energetically separated) ground-state wave function. This indicates that under certain conditions coupling to the quantum vacuum of a cavity can indeed modify ground state properties.

preprint2020arXiv

Moiré metrology of energy landscapes in van der Waals heterostructures

The emerging field of twistronics, which harnesses the twist angle between two-dimensional materials, represents a promising route for the design of quantum materials, as the twist-angle-induced superlattices offer means to control topology and strong correlations. At the small twist limit, and particularly under strain, as atomic relaxation prevails, the emergent moiré superlattice encodes elusive insights into the local interlayer interaction. Here we introduce moiré metrology as a combined experiment-theory framework to probe the stacking energy landscape of bilayer structures at the 0.1 meV/atom scale, outperforming the gold-standard of quantum chemistry. Through studying the shapes of moiré domains with numerous nano-imaging techniques, and correlating with multi-scale modelling, we assess and refine first-principle models for the interlayer interaction. We document the prowess of moiré metrology for three representative twisted systems: bilayer graphene, double bilayer graphene and H-stacked $MoSe_2/WSe_2$. Moiré metrology establishes sought after experimental benchmarks for interlayer interaction, thus enabling accurate modelling of twisted multilayers.

preprint2020arXiv

Nature of symmetry breaking at the excitonic insulator transition: Ta$_2$NiSe$_5$

Ta$_2$NiSe$_5$ is one of the most promising materials for hosting an excitonic insulator ground state. While a number of experimental observations have been interpreted in this way, the precise nature of the symmetry breaking occurring in Ta$_2$NiSe$_5$, the electronic order parameter, and a realistic microscopic description of the transition mechanism are, however, missing. By a symmetry analysis based on first-principles calculations, we uncover the \emph{discrete} lattice symmetries which are broken at the transition. We identify a purely electronic order parameter of excitonic nature that breaks these discrete crystal symmetries and contributes to the experimentally observed lattice distortion from an orthorombic to a monoclinic phase. Our results provide a theoretical framework to understand and analyze the excitonic transition in Ta$_2$NiSe$_5$ and settle the fundamental questions about symmetry breaking governing the spontaneous formation of excitonic insulating phases in solid-state materials.

preprint2020arXiv

Relevance of the quadratic diamagnetic and self-polarization terms in cavity quantum electrodynamics

Experiments at the interface of quantum-optics and chemistry have revealed that strong coupling between light and matter can substantially modify chemical and physical properties of molecules and solids. While the theoretical description of such situations is usually based on non-relativistic quantum electrodynamics, which contains quadratic light-matter coupling terms, it is commonplace to disregard these terms and restrict to purely bilinear couplings. In this work we clarify the physical origin and the substantial impact of the most common quadratic terms, the diamagnetic and self-polarization terms, and highlight why neglecting them can lead to rather unphysical results. Specifically we demonstrate its relevance by showing that neglecting it leads to the loss of gauge invariance, basis-set dependence, disintegration (loss of bound states) of any system in the basis set-limit, unphysical radiation of the ground state and an artificial dependence on the static dipole. Besides providing important guidance for modeling strongly coupled light-matter systems, the presented results do also indicate under which conditions those effects might become accessible.

preprint2020arXiv

Role of electron scattering on the high-order harmonic generation from solids

We extend the semi-classical trajectory description for the high-order harmonic generation~(HHG) from solids by integrating the effect of electron-scattering. Comparing the extended semi-classical trajectory model with a one-dimensional quantum mechanical simulation, we find that the multi-plateau feature of the HHG spectrum is formed by Umklapp scattering under the electron-hole acceleration dynamics by laser fields. Furthermore, by tracing the scattered trajectories in real-space, the model fairly describes the emitted photon energy and the emission timing of the HHG even in the higher plateau regions. We further consider the loss of trajectories by scattering processes with a mean-free-path approximation and evaluate the HHG cutoff energy as a function of laser wavelength. As a result, we find that the trajectory loss by scattering causes the wavelength independence of the HHG from solids.

preprint2020arXiv

Room Temperature Terahertz Electroabsorption Modulation by Excitons in Monolayer Transition Metal Dichalcogenides

The interaction between off-resonant laser pulses and excitons in monolayer transition metal dichalcogenides is attracting increasing interest as a route for the valley-selective coherent control of the exciton properties. Here, we extend the classification of the known off-resonant phenomena by unveiling the impact of a strong THz field on the excitonic resonances of monolayer MoS$_2$. We observe that the THz pump pulse causes a selective modification of the coherence lifetime of the excitons, while keeping their oscillator strength and peak energy unchanged. We rationalize these results theoretically by invoking a hitherto unobserved manifestation of the Franz-Keldysh effect on an exciton resonance. As the modulation depth of the optical absorption reaches values as large as 0.05 dB/nm at room temperature, our findings open the way to the use of semiconducting transition metal dichalcogenides as compact and efficient platforms for high-speed electroabsorption devices.

preprint2020arXiv

Self-consistent potential correction for charged periodic systems

Supercell models are often used to calculate the electronic structure of local perturbations from the ideal periodicity in the bulk or on the surface of a crystal or in wires. When the defect or adsorbent is charged, a jellium counter charge is applied to maintain overall neutrality, but the interaction of the artificially repeated charges has to be corrected, both in the total energy and in the one-electron eigenvalues and eigenstates. This becomes paramount in slab or wire calculations, where the jellium counter charge may induce spurious states in the vacuum. We present here a self-consistent potential correction scheme and provide successful tests of it for bulk and slab calculations.

preprint2020arXiv

Setting the photoelectron clock through molecular alignment

The interaction of strong laser fields with matter intrinsically provides powerful tools to image transient dynamics with an extremely high spatiotemporal resolution. Here, we study strong-field ionisation of laser-aligned molecules and show a full real-time picture of the photoelectron dynamics in the combined action of the laser field and the molecular interaction. We demonstrate that the molecule has a dramatic impact on the overall strong-field dynamics: it sets the clock for the emission of electrons with a given rescattering kinetic energy. This result represents a benchmark for the seminal statements of molecular-frame strong-field physics and has strong impact on the interpretation of self-diffraction experiments. Furthermore, the resulting encoding of the time-energy relation in molecular-frame photoelectron momentum distributions shows the way of probing the molecular potential in real-time and accessing a deeper understanding of electron transport during strong-field interactions.

preprint2020arXiv

The spontaneous symmetry breaking in Ta$_2$NiSe$_5$ is structural in nature

The excitonic insulator is an electronically-driven phase of matter that emerges upon the spontaneous formation and Bose condensation of excitons. Detecting this exotic order in candidate materials is a subject of paramount importance, as the size of the excitonic gap in the band structure establishes the potential of this collective state for superfluid energy transport. However, the identification of this phase in real solids is hindered by the coexistence of a structural order parameter with the same symmetry as the excitonic order. Only a few materials are currently believed to host a dominant excitonic phase, Ta$_2$NiSe$_5$ being the most promising. Here, we test this scenario by using an ultrashort laser pulse to quench the broken-symmetry phase of this transition metal chalcogenide. Tracking the dynamics of the material's electronic and crystal structure after light excitation reveals surprising spectroscopic fingerprints that are only compatible with a primary order parameter of phononic nature. We rationalize our findings through state-of-the-art calculations, confirming that the structural order accounts for most of the electronic gap opening. Not only do our results uncover the long-sought mechanism driving the phase transition of Ta$_2$NiSe$_5$, but they also conclusively rule out any substantial excitonic character in this instability.

preprint2020arXiv

Unravelling the intertwined atomic and bulk nature of localised excitons by attosecond spectroscopy

The electro-optical properties of most semiconductors and insulators of technological interest are dominated by the presence of electron-hole quasiparticles called excitons. The manipulation of these hydrogen-like quasi-particles in dielectrics, has received great interest under the name excitonics that is expected to be of great potential for a variety of applications, including optoelectronics and photonics. A crucial step for such exploitation of excitons in advanced technological applications is a detailed understanding of their dynamical nature. However, the ultrafast processes unfolding on few-femtosecond and attosecond time scales, of primary relevance in view of the desired extension of electronic devices towards the petahertz regime, remain largely unexplored. Here we apply attosecond transient reflection spectroscopy in a sequential two-foci geometry and observe sub-femtosecond dynamics of a core-level exciton in bulk MgF$_2$ single crystals. With our unique setup, we can access absolute phase delays which allow for an unambiguous comparison with theoretical calculations based on the Wannier-Mott model. Our results show that excitons surprisingly exhibit a dual atomic- and solid-like character which manifests itself on different time scales. While the former is responsible for a femtosecond optical Stark effect, the latter dominates the attosecond excitonic response and originates by the interaction with the crystal. Further investigation of the role of exciton localization proves that the bulk character persists also for strongly localised quasi-particles and allows us to envision a new route to control exciton dynamics in the close-to-petahertz regime.

preprint2020arXiv

Vibrational coherent control of localized d-d electronic excitation

Addressing the role of quantum coherence in the interplay between the different matter constituents (electrons, phonons and spin) is a critical step towards understanding transition metal oxides and design complex materials with new functionalities. Here we use coherent vibrational control of onsite d-d electronic transitions in a model edge-sharing insulating transition metal oxide (CuGeO3) to single-out the effects of vibrational coherence in electron-phonon coupling. By comparing time domain experiments based on high and low frequency ultrashort pumps with a fully quantum description of phonon assisted absorption, we could distinguish the processes associated to incoherent thermal lattice fluctuations from those driven by the coherent motion of the atoms. In particular, while thermal fluctuation of the phonon bath uniformly increases the electronic absorption, the resonant excitation of phonon modes results also in light-induced transparency which is coherently controlled by the vibrational motion.

preprint2019arXiv

Benchmarking Semiclassical and Perturbative Methods for Real-time Simulations of Cavity-Bound Emission and Interference

We benchmark a selection of semiclassical and perturbative dynamics techniques by investigating the correlated evolution of a cavity-bound atomic system to assess their applicability to study problems involving strong light-matter interactions in quantum cavities. The model system of interest features spontaneous emission, interference, and strong coupling behaviour, and necessitates the consideration of vacuum fluctuations and correlated light-matter dynamics. We compare a selection of approximate dynamics approaches including fewest switches surface hopping, multi-trajectory Ehrenfest dynamics, linearized semiclasical dynamics, and partially linearized semiclassical dynamics. Furthermore, investigating self-consistent perturbative methods, we apply the Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy in the second Born approximation. With the exception of fewest switches surface hopping, all methods provide a reasonable level of accuracy for the correlated light-matter dynamics, with most methods lacking the capacity to fully capture interference effects.

preprint2019arXiv

Berry Curvature Engineering by Gating Two-Dimensional Antiferromagnets

Recent advances in tuning electronic, magnetic, and topological properties of two-dimensional (2D) magnets have opened a new frontier in the study of quantum physics and promised exciting possibilities for future quantum technologies. In this study, we find that the dual-gate technology can well tune the electronic and topological properties of antiferromagnetic (AFM) even septuple-layer (SL) MnBi$_2$Te$_4$ thin films. Under an out-of-plane electric field that breaks $\mathcal{PT}$ symmetry, the Berry curvature of the thin film could be engineered efficiently, resulting in a huge change of anomalous Hall (AH) signal. Beyond the critical electric field, the double-SL MnBi$_2$Te$_4$ thin film becomes a Chern insulator with a high Chern number of 3. We further demonstrate that such 2D material can be used as an AFM switch via electric-field control of the AH signal. These discoveries inspire the design of low-power memory prototype for future AFM spintronic applications.

preprint2019arXiv

Exact exchange-correlation potential of effectively interacting Kohn-Sham systems

Aiming to combine density functional theory (DFT) and wavefunction theory, we study a mapping from the many-body interacting system to an effectively-interacting Kohn-Sham system instead of a non-interacting Kohn-Sham system. Because a ground state of effectively-interacting systems requires having a solution for the correlated many-body wavefunctions, this provides a natural framework to many-body wavefunction theories such as the configuration interaction and the coupled cluster method in the formal theoretical framework of DFT. Employing simple one-dimensional two-electron systems -- namely, the one-dimensional helium atom, hydrogen molecule and heteronuclear diatomic molecule -- we investigate properties of many-body wavefunctions and exact exchange-correlation potentials of effectively-interacting Kohn-Sham systems. As a result, we find that the asymptotic behavior of the exact exchange-correlation potential can be controlled by optimizing that of the effective interaction. Furthermore, the typical features of the exact non-interacting Kohn-Sham system, namely a spiky feature and a step feature in the exchange-correlation potential for the molecular dissociation limit can be suppressed by a proper choice of the effective interaction. These findings open a possibility to construct numerically robust and efficient exchange-correlation potentials and functionals based on the effectively-interacting Kohn-Sham scheme.

preprint2019arXiv

Local Berry curvature signatures in dichroic angle-resolved photoelectron spectroscopy

Topologically nontrivial two-dimensional materials hold great promise for next-generation optoelectronic applications. However, measuring the Hall or spin-Hall response is often a challenge and practically limited to the ground state. An experimental technique for tracing the topological character in a differential fashion would provide useful insights. In this work, we show that circular dichroism angle-resolved photoelectron spectroscopy (ARPES) provides a powerful tool which can resolve the topological and quantum-geometrical character in momentum space. In particular, we investigate how to map out the signatures of the local Berry curvature by exploiting its intimate connection to the orbital angular momentum. A spin-resolved detection of the photoelectrons allows to extend the approach to spin-Chern insulators. Our predictions are corroborated by state-of-the art \emph{ab initio} simulations employing time-dependent density functional theory, complemented with model calculations. The present proposal can be extended to address topological properties in materials out of equilibrium in a time-resolved fashion.

preprint2019arXiv

Magic continuum in twisted bilayer WSe2

Emergent quantum phases driven by electronic interactions can manifest in materials with narrowly dispersing, i.e. "flat", energy bands. Recently, flat bands have been realized in a variety of graphene-based heterostructures using the tuning parameters of twist angle, layer stacking and pressure, and resulting in correlated insulator and superconducting states. Here we report the experimental observation of similar correlated phenomena in twisted bilayer tungsten diselenide (tWSe2), a semiconducting transition metal dichalcogenide (TMD). Unlike twisted bilayer graphene where the flat band appears only within a narrow range around a "magic angle", we observe correlated states over a continuum of angles, spanning 4 degree to 5.1 degree. A Mott-like insulator appears at half band filling that can be sensitively tuned with displacement field. Hall measurements supported by ab initio calculations suggest that the strength of the insulator is driven by the density of states at half filling, consistent with a 2D Hubbard model in a regime of moderate interactions. At 5.1 degree twist, we observe evidence of superconductivity upon doping away from half filling, reaching zero resistivity around 3 K. Our results establish twisted bilayer TMDs as a model system to study interaction-driven phenomena in flat bands with dynamically tunable interactions.

preprint2019arXiv

Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems

Over the last years extraordinary advances in experimental and theoretical tools have allowed us to monitor and control matter at short time and atomic scales with a high-degree of precision. An appealing and challenging route towards engineering materials with tailored properties is to find ways to design or selectively manipulate materials, especially at the quantum level. To this end, having a state-of-the-art ab initio computer simulation tool that enables a reliable and accurate simulation of light-induced changes in the physical and chemical properties of complex systems is of utmost importance. The first principles real-space-based Octopus project was born with that idea in mind, providing an unique framework allowing to describe non-equilibrium phenomena in molecular complexes, low dimensional materials, and extended systems by accounting for electronic, ionic, and photon quantum mechanical effects within a generalized time-dependent density functional theory framework. The present article aims to present the new features that have been implemented over the last few years, including technical developments related to performance and massive parallelism. We also describe the major theoretical developments to address ultrafast light-driven processes, like the new theoretical framework of quantum electrodynamics density-functional formalism (QEDFT) for the description of novel light-matter hybrid states. Those advances, and other being released soon as part of the Octopus package, will enable the scientific community to simulate and characterize spatial and time-resolved spectroscopies, ultrafast phenomena in molecules and materials, and new emergent states of matter (QED-materials).

preprint2019arXiv

Polaritonic Coupled-Cluster Theory

We develop coupled-cluster theory for systems of electrons strongly coupled to photons, providing a promising theoretical tool in polaritonic chemistry with a perspective of application to all types of fermion-boson coupled systems. We show benchmark results for model molecular Hamiltonians coupled to cavity photons. By comparing to full configuration interaction results for various ground-state properties and optical spectra, we demonstrate that our method captures all key features present in the exact reference, including Rabi splittings and multi-photon processes. Further, a path on how to incorporate our bosonic extension of coupled-cluster theory into existing quantum chemistry programs is given.

preprint2019arXiv

Ultrafast Transient Absorption Spectroscopy of the Charge-Transfer Insulator NiO: Beyond the Dynamical Franz-Keldysh Effect

We demonstrate that a dynamical modification of the Hubbard U in the model charge-transfer insulator NiO can be observed with state-of-the-art time-resolved absorption spectroscopy. Using a self-consistent time-dependent density functional theory plus U computational framework, we show that the dynamical modulation of screening and Hubbard U significantly changes the transient optical spectroscopy. Whereas we find the well-known dynamical Franz-Keldysh effect when the U is frozen, we observe a dynamical band-gap renormalization for dynamical U. The renormalization of the optical gap is found to be smaller than the renormalization of U. This work opens up the possibility of driving a light-induced transition from a charge-transfer into a Mott insulator phase.

preprint2018arXiv

Real-time solutions of coupled Ehrenfest-Maxwell-Pauli-Kohn-Sham equations: fundamentals, implementation, and nano-optical applications

We present the theoretical foundations and the implementation details of a density-functional approach for coupled photons, electrons, and effective nuclei in non-relativistic quantum electrodynamics. Starting point of the formalism is a generalization of the Pauli-Fierz field theory for which we establish a one-to-one correspondence between external fields and internal variables. Based on this correspondence, we introduce a Kohn-Sham construction which provides a computationally feasible approach for ab-initio light-matter interactions. In the mean-field limit for the effective nuclei the formalism reduces to coupled Ehrenfest-Maxwell-Pauli-Kohn-Sham equations. We present an implementation of the approach in the real-space real-time code Octopus. For the implementation we use the Riemann-Silberstein formulation of classical electrodynamics and rewrite Maxwell's equations in Schrödinger form. This allows us to use existing time-evolution algorithms developed for quantum-mechanical systems also for Maxwell's equations. Overall, our approach is ideally suited for applications in nano-optics, nano-plasmonics, (photo) electrocatalysis, light-matter coupling in 2D materials, cases where laser pulses carry orbital angular momentum, or light-tailored chemical reactions in optical cavities to name but a few.

preprint2017arXiv

Generation and evolution of spin-, valley- and layer-polarized excited carriers in inversion-symmetric WSe2

Manipulation of spin and valley degrees of freedom is a key step towards realizing novel quantum technologies, for which atomically thin transition metal dichalcogenides (TMDCs) have been established as promising candidates. In monolayer TMDCs, the lack of inversion symmetry gives rise to a spin-valley correlation of the band structure allowing for valley-selective electronic excitation with circularly polarized light. Here we show that, even in centrosymmetric samples of 2H-WSe2, circularly polarized light can generate spin-, valley- and layer-polarized excited states in the conduction band. Employing time- and angle-resolved photoemission spectroscopy (trARPES) with spin-selective excitation, the dynamics of valley and layer pseudospins of the excited carriers are investigated. Complementary time-dependent density functional theory (TDDFT) calculations of the excited state populations reveal a strong circular dichroism of the spin-, valley- and layer-polarizations and a pronounced 2D character of the excited states in the K valleys. We observe scattering of carriers towards the global minimum of the conduction band on a sub-100 femtosecond timescale to states with three-dimensional character facilitating inter-layer charge transfer. Our results establish the optical control of coupled spin-, valley- and layer-polarized states in centrosymmetric materials and suggest the suitability of TMDC multilayer materials for valleytronic and spintronic device concepts.

preprint2016arXiv

A first principles TDDFT framework for spin and time-resolved ARPES in periodic systems

We present a novel theoretical approach to simulate spin, time and angular-resolved photoelectron spectroscopy (ARPES) from first principles that is applicable to surfaces, thin films, few layer systems, and low-dimensional nanostructures. The method is based on a general formulation in the framework of time-dependent density functional theory (TDDFT) to describe the real time-evolution of electrons escaping from a surface under the effect of any external (arbitrary) laser field. By extending the so called t-SURFF method to periodic systems one can calculate the final photoelectron spectrum by collecting the flux of the ionization current trough an analysing surface. The resulting approach, that we named t-SURFFP, allows to describe a wide range of irradiation conditions without any assumption on the dynamics of the ionization process allowing for pump-probe simulations on an equal footing. To illustrate the wide scope of applicability of the method we present applications to graphene, mono- and bi-layer WSe$_2$, and hexagonal BN under different laser configurations.

preprint2016arXiv

Correlated Electron-Nuclear Dynamics with Conditional Wave Functions

The molecular Schrödinger equation is rewritten in terms of non-unitary equations of motion for the nuclei (or electrons) that depend parametrically on the configuration of an ensemble of generally defined electronic (or nuclear) trajectories. This scheme is exact and does not rely on the tracing-out of degrees of freedom. Hence, the use of trajectory-based statistical techniques can be exploited to circumvent the calculation of the computationally demanding Born-Oppenheimer potential-energy surfaces and non-adiabatic coupling elements. The concept of potential-energy surface is restored by establishing a formal connection with the exact factorization of the full wave function. This connection is used to gain insight from a simplified form of the exact propagation scheme.

preprint2016arXiv

Disentangling Vacancy Oxidation on Metallicity-Sorted Carbon Nanotubes

Pristine single-walled carbon nanotubes (SWCNTs) are rather inert to O$_2$ and N$_2$, which for low doses chemisorb only on defect sites or vacancies of the SWCNTs at the ppm level. However, very low doping has a major effect on the electronic properties and conductivity of the SWCNTs. Already at low O$_2$ doses (80 L), the X-ray photoelectron spectroscopy (XPS) O 1s signal becomes saturated, indicating nearly all the SWCNT's vacancies have been oxidized. As a result, probing vacancy oxidation on SWCNTs via XPS yields spectra with rather low signal-to-noise ratios, even for metallicity-sorted SWCNTs. We show that, even under these conditions, the first principles density functional theory calculated Kohn-Sham O 1s binding energies may be used to assign the XPS O 1s spectra for oxidized vacancies on SWCNTs into its individual components. This allows one to determine the specific functional groups or bonding environments measured. We find the XPS O 1s signal is mostly due to three O-containing functional groups on SWCNT vacancies: epoxy (C$_2$$>$O), carbonyl (C$_2$$>$C$=$O), and ketene (C$=$C$=$O), as ordered by abundance. Upon oxidation of nearly all the SWCNT's vacancies, the central peak's intensity for the metallic SWCNT sample is 60\% greater than for the semiconducting SWCNT sample. This suggests a greater abundance of O-containing defect structures on the metallic SWCNT sample. For both metallic and semiconducting SWCNTs, we find O$_2$ does not contribute to the measured XPS O~1s spectra.

preprint2016arXiv

Monitoring in real time the photon-dressing and undressing of quasiparticles from first principles time-resolved photoelectron spectroscopy

Optical pumping of solids creates a non-equilibrium electronic structure where electrons and photons combine to form quasiparticles of dressed electronic states. The resulting shift of electronic levels is known as the optical Stark effect, visible as a red shift in the optical spectrum. Here we show that in a pump-probe setup we can uniquely define a non-equilibrium quasiparticle bandstructure that can be directly measurable with photo-electron spectroscopy. The dynamical photon-dressing (and undressing) of the many-body electronic states can be monitored by pump-probe time and angular resolved photoelectron spectroscopy (tr-ARPES) as the photon-dressed bandstructure evolves in time depending on the pump-probe pulse overlap. The computed tr-ARPES spectrum agrees perfectly with the quasi-energy spectrum of Floquet theory at maximum overlap and goes to the the equilibrium bandstructure as the pump-probe overlap goes to zero. Additionally, we show how this time-dependent non-equilibrium quasiparticle structure can be understood to be the bandstructure underlying the optical Stark effect. The extension to spin-resolved PES can be used to predict asymmetric dichroic response linked to the valley selective optical excitations in monolayer transition metal dichalcogenides (TMDs).

preprint2016arXiv

Photoinduced Absorption within Single-Walled Carbon Nanotube Systems

We study the photoabsorption properties of photoactive bulk polymer/ fullerene/nanotube heterojunctions in the near-infrared region. By combining pump-probe spectroscopy and linear response time-dependent density functional theory within the random phase approximation (TDDFT-RPA) we elucidate the excited state dynamics of the $E_{11}$ transition within (6,5) and (7,5) single-walled carbon nanotubes (SWNTs) and combined with poly(3-hexylthiophene-2,5-diyl) (P3HT) and [6,6]-phenyl-C$_{61}$-butyric acid methyl ester (PCBM) in P3HT/PCBM/SWNT blended samples. We find the presence of a photoinduced absorption (PA) peak is related mainly to the width of the photobleach (PB) peak and the charge carrier density of the SWNT system. For mixed SWNT samples, the PB peak is too broad to observe the PA peak, whereas within P3HT/PCBM/SWNT blended samples P3HT acts as a hole acceptor, narrowing the PB peak by exciton delocalization, which reveals a PA peak. Our results suggest that the PA peak originates from a widening of the band gap in the presence of excited electrons and holes. These results have important implications for the development of new organic photovoltaic heterojunctions including SWNTs.

preprint2016arXiv

Square selenene and tellurene: novel group VI elemental 2D semi-Dirac materials and topological insulators

With first principles calculations, we predict a novel stable 2D layered structure for group VI elements Se and Te that we call square selenene and square tellurene, respectively. They have chair-like buckled structures similar to other layered materials such as silicene and germanene but with a square unit cell rather than hexagonal. This special structure gives rise to anisotropic band dispersions near the Fermi level that can be described by a generalized semi-Dirac Hamiltonian. We show that the considerably large band gap ($\sim$0.1 eV) opened by spin-orbit coupling makes square selenene and tellurene topological insulators, hosting non-trivial edge states. Therefore, square selenene and tellurene are promising materials for novel electronic and spintronic applications. Finally, we show that this new type of 2D elemental material can potentially be grown on proper substrates, such as a Au(100) surface.

preprint2016arXiv

Stable Dirac semi-metal in the allotrope of IV elements

Three dimensional topological Dirac semi-metals represent a novel state of quantum matter with exotic electronic properties, in which a pair of Dirac points with the linear dispersion along all momentum directions exist in the bulk. Herein, by using the first principles calculations, we discover a new metastable allotrope of Ge and Sn in the staggered layered dumbbell structure, named as germancite and stancite, to be Dirac semi-metals with a pair of Dirac points on its rotation axis. On the surface parallel to the rotation axis, a pair of topologically non-trivial Fermi arcs are observed and a Lifshitz transition is found by tuning the Fermi level. Furthermore, the quantum thin film of germancite is found to be an intrinsic quantum spin Hall insulator. These discoveries suggest novel physical properties and future applications of the new metastable allotrope of Ge and Sn.

preprint2016arXiv

Systematic construction of density functionals based on matrix product state computations

We propose a systematic procedure for the approximation of density functionals in density functional theory that consists of two parts. First, for the efficient approximation of a general density functional, we introduce an efficient ansatz whose non-locality can be increased systematically. Second, we present a fitting strategy that is based on systematically increasing a reasonably chosen set of training densities. We investigate our procedure in the context of strongly correlated fermions on a one-dimensional lattice in which we compute accurate training densities with the help of matrix product states. Focusing on the exchange-correlation energy, we demonstrate how an efficient approximation can be found that includes and systematically improves beyond the local density approximation. Importantly, this systematic improvement is shown for target densities that are quite different from the training densities.

preprint2016arXiv

Tailored pump-probe transient spectroscopy with time-dependent density-functional theory: controlling absorption spectra

Recent advances in laser technology allow us to follow electronic motion at its natural time-scale with ultra-fast time resolution, leading the way towards attosecond physics experiments of extreme precision. In this work, we assess the use of tailored pumps in order to enhance (or reduce) some given features of the probe absorption (for example, absorption in the visible range of otherwise transparent samples). This type of manipulation of the system response could be helpful for its full characterization, since it would allow us to visualize transitions that are dark when using unshaped pulses. In order to investigate these possibilities, we perform first a theoretical analysis of the non-equilibrium response function in this context, aided by one simple numerical model of the Hydrogen atom. Then, we proceed to investigate the feasibility of using time-dependent density-functional theory as a means to implement, theoretically, this absorption-optimization idea, for more complex atoms or molecules. We conclude that the proposed idea could in principle be brought to the laboratory: tailored pump pulses can excite systems into light-absorbing states. However, we also highlight the severe numerical and theoretical difficulties posed by the problem: large-scale non-equilibrium quantum dynamics are cumbersome, even with TDDFT, and the shortcomings of state-of-the-art TDDFT functionals may still be serious for these out-of-equilibrium situations.

preprint2016arXiv

Taming singularities of the diagrammatic many-body perturbation theory

In a typical scenario the diagrammatic many-body perturbation theory generates asymptotic series. Despite non-convergence, the asymptotic expansions are useful when truncated to a finite number of terms. This is the reason for popularity of leading-order methods such as $GW$ approximation in condensed matter, molecular and atomic physics. Emerging higher-order implementations suffer from the appearance of nonsimple poles in the frequency-dependent Green's functions and negative spectral densities making self-consistent determination of the electronic structure impossible. Here a method based on the Padé approximation for overcomming these difficulties is proposed and applied to the Hamiltonian describing a core electron coupled to a single plasmonic excitation. By solving the model purely diagrammatically, expressing the self-energy in terms of combinatorics of chord diagrams, and regularizing the diverging perturbative expansions using the Padé approximation the spectral function is determined self-consistently using 3111 diagrams up to the sixth order.

preprint2015arXiv

Comparing quasiparticle H$_2$O level alignment on anatase and rutile TiO$_2$

Knowledge of the molecular frontier levels' alignment in the ground state can be used to predict the photocatalytic activity of an interface. The position of the adsorbate's highest occupied molecular orbital (HOMO) levels relative to the substrate's valence band maximum (VBM) in the interface describes the favorability of photogenerated hole transfer from the VBM to the adsorbed molecule. This is a key quantity for assessing and comparing H$_2$O photooxidation activities on two prototypical photocatalytic TiO$_2$ surfaces: anatase (A)-TiO$_2$(101) and rutile (R)-TiO$_2$(110). Using the projected density of states (DOS) from state-of-the-art quasiparticle (QP) $G_0W_0$ calculations, we assess the relative photocatalytic activity of intact and dissociated H$_2$O on coordinately unsaturated (Ti$_{\textit{cus}}$) sites of idealized stoichiometric A-TiO$_2$(101)/R-TiO$_2$(110) and bridging O vacancies (O$_{\textit{br}}^{\textit{vac}}$) of defective A-TiO$_{2-x}$(101)/R-TiO$_{2-x}$(110) surfaces ($x=\frac{1}{4},\frac{1}{8}$) for various coverages. Such a many-body treatment is necessary to correctly describe the anisotropic screening of electron-electron interactions at a photocatalytic interface, and hence obtain accurate interfacial level alignments. The more favorable ground state HOMO level alignment for A-TiO$_2$(101) may explain why the anatase polymorph shows higher photocatalytic activities than the rutile polymorph. Our results indicate that (1) hole trapping is more favored on A-TiO$_2$(101) than R-TiO$_2$(110) and (2) HO@Ti$_{\textit{cus}}$ is more photocatalytically active than intact H$_2$O@Ti$_{\textit{cus}}$.

preprint2015arXiv

Exact Maps in Density Functional Theory for Lattice Models

In the present work, we employ exact diagonalization for model systems on a real-space lattice to explicitly construct the exact density-to-potential and for the first time the exact density-to-wavefunction map that underly the Hohenberg-Kohn theorem in density functional theory. Having the explicit wavefunction-to- density map at hand, we are able to construct arbitrary observables as functionals of the ground-state density. We analyze the density-to-potential map as the distance between the fragments of a system increases and the correlation in the system grows. We observe a feature that gradually develops in the density-to-potential map as well as in the density-to-wavefunction map. This feature is inherited by arbitrary expectation values as functional of the ground-state density. We explicitly show the excited-state energies, the excited-state densities, and the correlation entropy as functionals of the ground-state density. All of them show this exact feature that sharpens as the coupling of the fragments decreases and the correlation grows. We denominate this feature as intra-system steepening. We show that for fully decoupled subsystems the intra-system steepening transforms into the well-known inter-system derivative discontinuity. An important conclusion is that for e.g. charge transfer processes between localized fragments within the same system it is not the usual inter-system derivative discontinuity that is missing in common ground-state functionals, but rather the differentiable intra-system steepening that we illustrate in the present work.

preprint2015arXiv

Kohn-Sham Approach to Quantum Electrodynamical Density Functional Theory: Exact Time-Dependent Effective Potentials in Real Space

The density-functional approach to quantum electrodynamics is extending traditional density-functional theory and opens the possibility to describe electron-photon interactions in terms of effective Kohn-Sham potentials. In this work, we numerically construct the exact electron-photon Kohn-Sham potentials for a prototype system which consists of a trapped electron coupled to a quantized electromagnetic mode in an optical high-Q cavity. While the effective current that acts on the photons is known explicitly, the exact effective potential that describes the forces exerted by the photons on the electrons is obtained from a fixed-point inversion scheme. This procedure allows us to uncover important beyond-mean-field features of the effective potential which mark the breakdown of classical light-matter interactions. We observe peak and step structures in the effective potentials, which can be attributed solely to the quantum nature of light, i.e., they are real-space signatures of the photons. Our findings show how the ubiquitous dipole interaction with a classical electromagnetic field has to be modified in real-space in order to take the quantum nature of the electromagnetic field fully into account.

preprint2015arXiv

Orbitals from local RDMFT: Are they Kohn-Sham or Natural Orbitals?

Recently, an approximate theoretical framework was introduced, called local reduced density matrix functional theory (local-RDMFT), where functionals of the one-body reduced density matrix (1-RDM) are minimized under the additional condition that the optimal orbitals satisfy a single electron Schrödinger equation with a local potential. In the present work, we focus on the character of these optimal orbitals. In particular, we compare orbitals obtained by local-RDMFT with those obtained with the full minimization (without the extra condition) by contrasting them against the exact NOs and orbitals from a density functional calculation using the local density approximation (LDA). We find that the orbitals from local-RMDFT are very close to LDA orbitals, contrary to those of the full minimization that resemble the exact NOs. Since local RDMFT preserves the good quality of the description of strong static correlation, this finding opens the way to a mixed density/density matrix scheme, where Kohn-Sham orbitals obtain fractional occupations from a minimization of the occupation numbers using 1-RDM functionals. This will allow for a description of strong correlation at a cost only minimally higher than a density functional calculation.

preprint2015arXiv

Static correlation and electron localization in molecular dimers from the self-consistent RPA and GW approximation

We investigate static correlation and delocalization errors in the self-consistent GW and random-phase approximation (RPA) by studying molecular dissociation of the H_2 and LiH molecules. Although both approximations contain topologically identical diagrams, the non-locality and frequency dependence of the GW self-energy crucially influence the different energy contributions to the total energy as compared to the use of a static local potential in the RPA. The latter leads to significantly larger correlation energies which allow for a better description of static correlation at intermediate bond distances. The substantial error found in GW is further analyzed by comparing spin-restricted and spin-unrestricted calculations. At large but finite nuclear separation their difference gives an estimate of the so-called fractional spin error normally determined only in the dissociation limit. Furthermore, a calculation of the dipole moment of the LiH molecule at dissociation reveals a large delocalization error in GW making the fractional charge error comparable to the RPA. The analyses are supplemented by explicit formulae for the GW Green's function and total energy of a simplified two-level model providing additional insights into the dissociation limit.

preprint2015arXiv

Time-dependent thermal transport theory

Understanding thermal transport in nanoscale systems presents important challenges to both theory and experiment. In particular, the concept of local temperature at the nanoscale appears difficult to justify. Here, we propose a novel theoretical approach where we replace the temperature gradient with controllable external black-body radiations. The theory recovers known physical results, for example the linear relation between the thermal current and the temperatures difference of two black-bodies. Furthermore, our theory is not limited to the linear regime and goes beyond accounting for non linear effects and transient phenomena. In the strong coupling and large temperature gradients, we show that there is a maximum energy current that the system can sustain and we recover the Kramers' turnover. Since the present theory is general and can be adapted to describe both electron and phonon dynamics, it provides a first step towards a unified formalism for investigating thermal and electronic transport.

Angel Rubio

What is connected

Connect this record

See the researcher in context

Building this map preview

137 published item(s)

Dynamical instability in a Floquet-Driven Dissipative System

Vacuum-dressed superconductivity in NbN observed in a high-$Q$ terahertz cavity

Ab initio calculations of quantum light-matter interactions in general electromagnetic environments

A perspective on ab initio modeling of polaritonic chemistry: The role of non-equilibrium effects and quantum collectivity

Ab-initio cluster approach for high-harmonic generation in liquids

Common microscopic origin of the phase transitions in Ta$_2$NiS$_5$ and the excitonic insulator candidate Ta$_2$NiSe$_5$

Exact Solution for A Real Polaritonic System Under Vibrational Strong Coupling in Thermodynamic Equilibrium: Absence of Zero Temperature and Loss of Light-Matter Entanglement

Floquet engineering the band structure of materials with optimal control theory

Frequency-dependent Sternheimer linear-response formalism for strongly coupled light-matter systems

Moiré Engineering of Nonsymmorphic Symmetries and Hourglass Superconductors

Moiré Engineering of Spin-Orbit Coupling in Twisted Platinum Diselenide

New Class of Landau Levels and Hall Phases in a 2D Electron Gas Subject to an Inhomogeneous Magnetic Field: An Analytic Solution

Polaritonic Hofstadter Butterfly and Cavity-Control of the Quantized Hall Conductance

Probing phonon dynamics with multi-dimensional high harmonic carrier envelope phase spectroscopy

Quadrupolar excitons in a tunnel-coupled van der Waals heterotrilayer

Soft-x-ray spectroscopy of monolayer transition metal dichalcogenides: valley selectivity?

Spin-orbit induced equilibrium spin currents in materials

Terahertz-induced high-order harmonic generation and nonlinear charge transport in graphene

Time- and angle-resolved photoelectron spectroscopy of strong-field light-dressed solids: prevalence of the adiabatic band picture

Unconventional Excitonic States with Phonon Sidebands in Layered Silicon Diphosphide

Wavefunction Embedding for Molecular Polaritons

Approximations based on density-matrix embedding theory for density-functional theories

Correlation-driven sub-3 fs charge migration in ionised adenine

Detecting multiple chirality centers in chiral molecules with high harmonic generation

Moiré heterostructures as a condensed matter quantum simulator

Nanometer-scale lateral p-n junctions in graphene/$α$-RuCl$_3$ heterostructures

Photoionization and transient Wannier-Stark ladder in silicon: First principle simulations versus Keldysh theory

The Free Electron Gas in Cavity Quantum Electrodynamics

TMDs as a platform for spin liquid physics: A strong coupling study of twisted bilayer WSe$_2$

Ultra-Strong Spin-Orbit Coupling and Topological Moiré Engineering in Twisted ZrS2 Bilayers

Ultrafast Dynamical Lifshitz Transition

Chemistry in Quantum Cavities: Exact Results, the Impact of Thermal Velocities and Modified Dissociation

Effect of Many Modes on Self-Polarization and Photochemical Suppression in Cavities

Electron-Phonon-Driven Three-Dimensional Metallicity in an Insulating Cuprate

Giant exciton Mott density in anatase TiO2

High-Harmonic generation from spin-polarised defects in solids

High-harmonics and isolated attosecond pulses from MgO

High-order harmonic generation in graphene: nonlinear coupling of intra and interband transitions

Light-matter hybrid-orbital-based first-principles methods: the influence of the polariton statistics

Moiré metrology of energy landscapes in van der Waals heterostructures

Nature of symmetry breaking at the excitonic insulator transition: Ta$_2$NiSe$_5$

Relevance of the quadratic diamagnetic and self-polarization terms in cavity quantum electrodynamics

Role of electron scattering on the high-order harmonic generation from solids

Room Temperature Terahertz Electroabsorption Modulation by Excitons in Monolayer Transition Metal Dichalcogenides

Self-consistent potential correction for charged periodic systems

Setting the photoelectron clock through molecular alignment

The spontaneous symmetry breaking in Ta$_2$NiSe$_5$ is structural in nature

Unravelling the intertwined atomic and bulk nature of localised excitons by attosecond spectroscopy

Vibrational coherent control of localized d-d electronic excitation

Benchmarking Semiclassical and Perturbative Methods for Real-time Simulations of Cavity-Bound Emission and Interference

Berry Curvature Engineering by Gating Two-Dimensional Antiferromagnets

Exact exchange-correlation potential of effectively interacting Kohn-Sham systems

Local Berry curvature signatures in dichroic angle-resolved photoelectron spectroscopy

Magic continuum in twisted bilayer WSe2

Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems

Polaritonic Coupled-Cluster Theory

Ultrafast Transient Absorption Spectroscopy of the Charge-Transfer Insulator NiO: Beyond the Dynamical Franz-Keldysh Effect

Real-time solutions of coupled Ehrenfest-Maxwell-Pauli-Kohn-Sham equations: fundamentals, implementation, and nano-optical applications

Generation and evolution of spin-, valley- and layer-polarized excited carriers in inversion-symmetric WSe2

A first principles TDDFT framework for spin and time-resolved ARPES in periodic systems

Correlated Electron-Nuclear Dynamics with Conditional Wave Functions

Disentangling Vacancy Oxidation on Metallicity-Sorted Carbon Nanotubes

Monitoring in real time the photon-dressing and undressing of quasiparticles from first principles time-resolved photoelectron spectroscopy

Photoinduced Absorption within Single-Walled Carbon Nanotube Systems

Square selenene and tellurene: novel group VI elemental 2D semi-Dirac materials and topological insulators

Stable Dirac semi-metal in the allotrope of IV elements

Systematic construction of density functionals based on matrix product state computations

Tailored pump-probe transient spectroscopy with time-dependent density-functional theory: controlling absorption spectra

Taming singularities of the diagrammatic many-body perturbation theory

Comparing quasiparticle H$_2$O level alignment on anatase and rutile TiO$_2$

Exact Maps in Density Functional Theory for Lattice Models

Kohn-Sham Approach to Quantum Electrodynamical Density Functional Theory: Exact Time-Dependent Effective Potentials in Real Space

Orbitals from local RDMFT: Are they Kohn-Sham or Natural Orbitals?

Static correlation and electron localization in molecular dimers from the self-consistent RPA and GW approximation