Exchange-only virial relation from the adiabatic connection.

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

Exchange energies with forces in density-functional theory.

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

High-harmonic spectroscopy of low-energy electron-scattering dynamics in liquids.

High-harmonic spectroscopy is an all-optical nonlinear technique with inherent attosecond temporal resolution. It has been applied to a variety of systems in the gas phase and solid state. Here we extend its use to liquid samples. By studying high-harmonic generation over a broad range of wavelengths and intensities, we show that the cut-off energy is independent of the wavelength beyond a threshold intensity and that it is a characteristic property of the studied liquid. We explain these observations with a semi-classical model based on electron trajectories that are limited by the electron scattering. This is further confirmed by measurements performed with elliptically polarized light and with ab-initio time-dependent density functional theory calculations. Our results propose high-harmonic spectroscopy as an all-optical approach for determining the effective mean free paths of slow electrons in liquids. This regime is extremely difficult to access with other methodologies, but is critical for understanding radiation damage to living tissues. Our work also indicates the possibility of resolving subfemtosecond electron dynamics in liquids offering an all-optical approach to attosecond spectroscopy of chemical processes in their native liquid environment.

Self-interaction correction schemes for non-collinear spin-density-functional theory.

We extend some of the well-established self-interaction correction (SIC) schemes of density-functional theory-the Perdew-Zunger SIC and the average-density SIC-to the case of systems with noncollinear magnetism. Our proposed SIC schemes are tested on a set of molecules and metallic clusters in combination with the widely used local spin-density approximation. As expected from the collinear SIC, we show that the averaged-density SIC works well for improving ionization energies but fails to improve more subtle quantities like the dipole moments of polar molecules. We investigate the exchange-correlation magnetic field produced by our extension of the Perdew-Zunger SIC, showing that it is not aligned with the local total magnetization, thus producing an exchange-correlation torque.

Phonon-enhanced nonlinearities in hexagonal boron nitride.

Polar crystals can be driven into collective oscillations by optical fields tuned to precise resonance frequencies. As the amplitude of the excited phonon modes increases, novel processes scaling non-linearly with the applied fields begin to contribute to the dynamics of the atomic system. Here we show two such optical nonlinearities that are induced and enhanced by the strong phonon resonance in the van der Waals crystal hexagonal boron nitride (hBN). We predict and observe large sub-picosecond duration signals due to four-wave mixing (FWM) during resonant excitation. The resulting FWM signal allows for time-resolved observation of the crystal motion. In addition, we observe enhancements of third-harmonic generation with resonant pumping at the hBN transverse optical phonon. Phonon-induced nonlinear enhancements are also predicted to yield large increases in high-harmonic efficiencies beyond the third.

High-harmonic generation in liquids with few-cycle pulses: effect of laser-pulse duration on the cut-off energy.

High-harmonic generation (HHG) in liquids is opening new opportunities for attosecond light sources and attosecond time-resolved studies of dynamics in the liquid phase. In gas-phase HHG, few-cycle pulses are routinely used to create isolated attosecond pulses and to extend the cut-off energy. Here, we study the properties of HHG in liquids, including heavy water, ethanol and isopropanol, by continuously tuning the pulse duration of a mid-infrared driver from the multi- to the two-cycle regime. Similar to the gas phase, we observe the transition from discrete odd-order harmonics to continuous extreme-ultraviolet emission. However, the cut-off energy is shown to be entirely independent of the pulse duration. These observations are confirmed by ab-initio simulations of HHG in large liquid clusters. Our results support the notion that the cut-off energy is a fundamental property of the liquid, independent of the driving-pulse properties. Our work implies that few-cycle mid-infrared laser pulses are suitable drivers for generating isolated attosecond pulses from liquids and confirm the capability of high-harmonic spectroscopy to determine the mean-free paths of slow electrons in liquids.

Ultrafast Spin Dynamics and Photoinduced Insulator-to-Metal Transition in α-RuCl3.

Laser-induced ultrafast demagnetization is a phenomenon of utmost interest and attracts significant attention because it enables potential applications in ultrafast optoelectronics and spintronics. As a spin-orbit coupling assisted magnetic insulator, α-RuCl3 provides an attractive platform to explore the physics of electronic correlations and unconventional magnetism. Using time-dependent density functional theory, we explore the ultrafast laser-induced dynamics of the electronic and magnetic structures in α-RuCl3. Our study unveils that laser pulses can introduce ultrafast demagnetizations, accompanied by an out-of-equilibrium insulator-to-metal transition in a few tens of femtoseconds. The spin response significantly depends on the laser wavelength and polarization on account of the electron correlations, band renormalizations, and charge redistributions. These findings provide physical insights into the coupling between the electronic and magnetic degrees of freedom in α-RuCl3 and shed light on suppressing the long-range magnetic orders and reaching a proximate spin liquid phase for two-dimensional magnets on an ultrafast time scale.

Polarization spectroscopy of high-order harmonic generation in gallium arsenide.

An interesting property of high harmonic generation in solids is its laser polarization dependent nature which in turn provides information about the crystal and band structure of the generation medium. Here we report on the linear polarization dependence of high-order harmonic generation from a gallium arsenide crystal. Interestingly, we observe a significant evolution of the anisotropic response of above bandgap harmonics as a function of the laser intensity. We attribute this change to fundamental microscopic effects of the emission process comprising a competition between intraband and interband dynamics. This intensity dependence of the anisotropic nature of the generation process offers the possibility to drive and control the electron current along preferred directions of the crystal, and could serve as a switching technique in an integrated all-solid-state petahertz optoelectronic device.

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 exists, 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 nonpolar 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, including a structural minimum at 15-17 eV that is associated solely with the liquid phase. Our results pave the way to accessible calculations of HHG in liquids and illustrate the unique nonlinear nature of liquid systems.

Erratum: Identification of the Mott Insulating Charge Density Wave State in 1T-TaS_{2} [Phys. Rev. Lett. 126, 196406 (2021)].

This corrects the article DOI: 10.1103/PhysRevLett.126.196406.

Light-Driven Extremely Nonlinear Bulk Photogalvanic Currents.

We predict the generation of bulk photocurrents in materials driven by bichromatic fields that are circularly polarized and corotating. The nonlinear photocurrents have a fully controllable directionality and amplitude without requiring carrier-envelope-phase stabilization or few-cycle pulses, and can be generated with photon energies much smaller than the band gap (reducing heating in the photoconversion process). We demonstrate with ab initio calculations that the photocurrent generation mechanism is universal and arises in gaped materials (Si, diamond, MgO, hBN), in semimetals (graphene), and in two- and three-dimensional systems. Photocurrents are shown to rely on sub-laser-cycle asymmetries in the nonlinear response that build-up coherently from cycle to cycle as the conduction band is populated. Importantly, the photocurrents are always transverse to the major axis of the co-circular lasers regardless of the material's structure and orientation (analogously to a Hall current), which we find originates from a generalized time-reversal symmetry in the driven system. At high laser powers (∼10^{13} W/cm^{2}) this symmetry can be spontaneously broken by vast electronic excitations, which is accompanied by an onset of carrier-envelope-phase sensitivity and ultrafast many-body effects. Our results are directly applicable for efficient light-driven control of electronics, and for enhancing sub-band-gap bulk photogalvanic effects.

Identification of the Mott Insulating Charge Density Wave State in 1T-TaS_{2}.

We investigate the low-temperature charge density wave (CDW) state of bulk TaS_{2} with a fully self-consistent density-functional theory with the Hubbard U potential, over which the controversy has remained unresolved regarding the out-of-plane metallic band. By examining the innate structure of the Hubbard U potential, we reveal that the conventional use of atomic-orbital basis could seriously misevaluate the electron correlation in the CDW state. By adopting a generalized basis, covering the whole David star, we successfully reproduce the Mott insulating nature with the layer-by-layer antiferromagnetic order. Similar consideration should be applied for description of the electron correlation in molecular solid.

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 nonequilibrium 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 scheme for driving an ultrafast Lifshitz transition in the correlated type-II Weyl semimetal T d-MoTe2 We demonstrate that this nonequilibrium topological electronic transition finds its microscopic origin in the dynamical modification of the effective electronic correlations. These results shed light on a previously unexplored ultrafast scheme for controlling the Fermi surface topology in correlated quantum materials.

Enhanced tunable second harmonic generation from twistable interfaces and vertical superlattices in boron nitride homostructures.

Broken symmetries induce strong even-order nonlinear optical responses in materials and at interfaces. Unlike conventional covalently bonded nonlinear crystals, van der Waals (vdW) heterostructures feature layers that can be stacked at arbitrary angles, giving complete control over the presence or lack of inversion symmetry at a crystal interface. Here, we report highly tunable second harmonic generation (SHG) from nanomechanically rotatable stacks of bulk hexagonal boron nitride (BN) crystals and introduce the term twistoptics to describe studies of optical properties in twistable vdW systems. By suppressing residual bulk effects, we observe SHG intensity modulated by a factor of more than 50, and polarization patterns determined by moiré interface symmetry. Last, we demonstrate greatly enhanced conversion efficiency in vdW vertical superlattice structures with multiple symmetry-broken interfaces. Our study paves the way for compact twistoptics architectures aimed at efficient tunable frequency conversion and demonstrates SHG as a robust probe of buried vdW interfaces.

Time-Resolved Exciton Wave Functions from Time-Dependent Density-Functional Theory.

Time-dependent density-functional theory (TDDFT) is a computationally efficient first-principles approach for calculating optical spectra in insulators and semiconductors including excitonic effects. We show how exciton wave functions can be obtained from TDDFT via the Kohn-Sham transition density matrix, both in the frequency-dependent linear-response regime and in real-time propagation. The method is illustrated using one-dimensional model solids. In particular, we show that our approach provides insight into the formation and dissociation of excitons in real time. This opens the door to time-resolved studies of exciton dynamics in materials by means of real-time TDDFT.

The CECAM electronic structure library and the modular software development paradigm.

First-principles electronic structure calculations are now accessible to a very large community of users across many disciplines, thanks to many successful software packages, some of which are described in this special issue. The traditional coding paradigm for such packages is monolithic, i.e., regardless of how modular its internal structure may be, the code is built independently from others, essentially from the compiler up, possibly with the exception of linear-algebra and message-passing libraries. This model has endured and been quite successful for decades. The successful evolution of the electronic structure methodology itself, however, has resulted in an increasing complexity and an ever longer list of features expected within all software packages, which implies a growing amount of replication between different packages, not only in the initial coding but, more importantly, every time a code needs to be re-engineered to adapt to the evolution of computer hardware architecture. The Electronic Structure Library (ESL) was initiated by CECAM (the European Centre for Atomic and Molecular Calculations) to catalyze a paradigm shift away from the monolithic model and promote modularization, with the ambition to extract common tasks from electronic structure codes and redesign them as open-source libraries available to everybody. Such libraries include "heavy-duty" ones that have the potential for a high degree of parallelization and adaptation to novel hardware within them, thereby separating the sophisticated computer science aspects of performance optimization and re-engineering from the computational science done by, e.g., physicists and chemists when implementing new ideas. We envisage that this modular paradigm will improve overall coding efficiency and enable specialists (whether they be computer scientists or computational scientists) to use their skills more effectively and will lead to a more dynamic evolution of software in the community as well as lower barriers to entry for new developers. The model comes with new challenges, though. The building and compilation of a code based on many interdependent libraries (and their versions) is a much more complex task than that of a code delivered in a single self-contained package. Here, we describe the state of the ESL, the different libraries it now contains, the short- and mid-term plans for further libraries, and the way the new challenges are faced. The ESL is a community initiative into which several pre-existing codes and their developers have contributed with their software and efforts, from which several codes are already benefiting, and which remains open to the community.

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

Over the last few 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 toward 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, i.e., to provide a unique framework that allows us 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. This 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, such as the new theoretical framework of quantum electrodynamics density-functional formalism for the description of novel light-matter hybrid states. Those advances, and others being released soon as part of the Octopus package, will allow the scientific community to simulate and characterize spatial and time-resolved spectroscopies, ultrafast phenomena in molecules and materials, and new emergent states of matter (quantum electrodynamical-materials).

Multiflat Bands and Strong Correlations in Twisted Bilayer Boron Nitride: Doping-Induced Correlated Insulator and Superconductor.

Two-dimensional materials, obtained by van der Waals stacking of layers, are fascinating objects of contemporary condensed matter research, exhibiting a variety of new physics. Inspired by the breakthroughs of twisted bilayer graphene (TBG), we demonstrate that twisted bilayer boron nitride (TBBN) is an even more exciting novel system that turns out to be an excellent platform to realize new correlated phases and phenomena; exploration of its electronic properties shows that in contrast to TBG in TBBN multiple families of 2,4, and 6-fold degenerate flat bands emerge without the need to fine tune close to a "magic angle", resulting in dramatic and tunable changes in optical properties and exciton physics, and providing an additional platform to study strong correlations. Upon doping, unforeseen new correlated phases of matter (insulating and superconducting) emerge. TBBN could thus provide a promising experimental platform, insensitive to small deviations in the twist angle, to study novel exciton condensate and spatial confinement physics, and correlations in two dimensions.

All-optical nonequilibrium pathway to stabilising magnetic Weyl semimetals in pyrochlore iridates.

Nonequilibrium many-body dynamics is becoming a central topic in condensed matter physics. Floquet topological states were suggested to emerge in photodressed bands under periodic laser driving. Here we propose a viable nonequilibrium route without requiring coherent Floquet states to reach the elusive magnetic Weyl semimetallic phase in pyrochlore iridates by ultrafast modification of the effective electron-electron interaction with short laser pulses. Combining ab initio calculations for a time-dependent self-consistent light-reduced Hubbard U and nonequilibrium magnetism simulations for quantum quenches, we find dynamically modified magnetic order giving rise to transiently emerging Weyl cones that can be probed by time- and angle-resolved photoemission spectroscopy. Our work offers a unique and realistic pathway for nonequilibrium materials engineering beyond Floquet physics to create and sustain Weyl semimetals. This may lead to ultrafast, tens-of-femtoseconds switching protocols for light-engineered Berry curvature in combination with ultrafast magnetism.