The ferroelectric photo ground state of SrTiO3: Cavity materials engineering.

Optical cavities confine light on a small region in space, which can result in a strong coupling of light with materials inside the cavity. This gives rise to new states where quantum fluctuations of light and matter can alter the properties of the material altogether. Here we demonstrate, based on first-principles calculations, that such light-matter coupling induces a change of the collective phase from quantum paraelectric to ferroelectric in the [Formula: see text] ground state, which has thus far only been achieved in out-of-equilibrium strongly excited conditions [X. Li et al., Science 364, 1079-1082 (2019) and T. F. Nova, A. S. Disa, M. Fechner, A. Cavalleri, Science 364, 1075-1079 (2019)]. This is a light-matter hybrid ground state which can only exist because of the coupling to the vacuum fluctuations of light, a photo ground state The phase transition is accompanied by changes in the crystal structure, showing that fundamental ground state properties of materials can be controlled via strong light-matter coupling. Such a control of quantum states enables the tailoring of materials properties or even the design of novel materials purely by exposing them to confined light.

Reconstruction of ultrafast exciton dynamics with a phase-retrieval algorithm.

The first step to gain optical control over the ultrafast processes initiated by light in solids is a correct identification of the physical mechanisms at play. Among them, exciton formation has been identified as a crucial phenomenon which deeply affects the electro-optical properties of most semiconductors and insulators of technological interest. While recent experiments based on attosecond spectroscopy techniques have demonstrated the possibility to observe the early-stage exciton dynamics, the description of the underlying exciton properties remains non-trivial. In this work we propose a new method called extended Ptychographic Iterative engine for eXcitons (ePIX), capable of reconstructing the main physical properties which determine the evolution of the quasi-particle with no prior knowledge of the exact relaxation dynamics or the pump temporal characteristics. By demonstrating its accuracy even when the exciton dynamics is comparable to the pump pulse duration, ePIX is established as a powerful approach to widen our knowledge of solid-state physics.

Simulating Terahertz Field-Induced Ferroelectricity in Quantum Paraelectric SrTiO_{3}.

Recent experiments have demonstrated that light can induce a transition from the quantum paraelectric to the ferroelectric phase of SrTiO_{3}. Here, we investigate this terahertz field-induced ferroelectric phase transition by solving the time-dependent lattice Schrödinger equation based on first-principles calculations. We find that ferroelectricity originates from a light-induced mixing between ground and first excited lattice states in the quantum paraelectric phase. In agreement with the experimental findings, our study shows that the nonoscillatory second harmonic generation signal can be evidence of ferroelectricity in SrTiO_{3}. We reveal the microscopic details of this exotic phase transition and highlight that this phenomenon is a unique behavior of the quantum paraelectric phase.

Unraveling materials Berry curvature and Chern numbers from real-time evolution of Bloch states.

Materials can be classified by the topological character of their electronic structure and, in this perspective, global attributes immune to local deformations have been discussed in terms of Berry curvature and Chern numbers. Except for instructional simple models, linear response theories have been ubiquitously used in calculations of topological properties of real materials. Here we propose a completely different and versatile approach to obtain the topological characteristics of materials by calculating physical observables from the real-time evolving Bloch states: The cell-averaged current density reveals the anomalous velocities that lead to the conductivity quantum. Results for prototypical cases are shown, including a spin-frozen valley Hall and a quantum anomalous Hall insulator. The advantage of this method is best illustrated by the example of a quantum spin Hall insulator: The quantized spin Hall conductivity is straightforwardly obtained irrespective of the non-Abelian nature in its Berry curvature. Moreover, the method can be extended to the description of real observables in nonequilibrium states of topological materials.

Survival of Floquet-Bloch States in the Presence of Scattering.

Floquet theory has spawned many exciting possibilities for electronic structure control with light, with enormous potential for future applications. The experimental demonstration in solids, however, remains largely unrealized. In particular, the influence of scattering on the formation of Floquet-Bloch states remains poorly understood. Here we combine time- and angle-resolved photoemission spectroscopy with time-dependent density functional theory and a two-level model with relaxation to investigate the survival of Floquet-Bloch states in the presence of scattering. We find that Floquet-Bloch states will be destroyed if scattering-activated by electronic excitations-prevents the Bloch electrons from following the driving field coherently. The two-level model also shows that Floquet-Bloch states reappear at high field intensities where energy exchange with the driving field dominates over energy dissipation to the bath. Our results clearly indicate the importance of long scattering times combined with strong driving fields for the successful realization of various Floquet phenomena.

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 MoS2. 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.

Direct Measurement of Electron-Phonon Coupling with Time-Resolved ARPES.

Time- and angular- resolved photoelectron spectroscopy is a powerful technique to measure electron dynamics in solids. Recent advances in this technique have facilitated band and energy resolved observations of the effect that excited phonons, have on the electronic structure. Here, we show with the help of ab initio simulations that the Fourier analysis of the time-resolved measurements of solids with excited phonon modes enables the determination of the band- and mode-resolved electron-phonon coupling directly from the experimental data without any additional input from theory. Such an observation is not restricted to regions of strong electron-phonon coupling and does not require strongly excited or hot phonons, but can be employed to monitor the dynamical renormalization of phonons in driven phases of matter.

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).

Nonlinear polarization evolution using time-dependent density functional theory.

We propose a theoretical and computational approach to investigate temporal behavior of a nonlinear polarization in a perturbative regime induced by an intense and ultrashort pulsed electric field. First-principles time-dependent density functional theory is employed to describe the electron dynamics. Temporal evolution of third-order nonlinear polarization is extracted from a few calculations of electron dynamics induced by pulsed electric fields with the same time profile but different amplitudes. We discuss characteristic features of the nonlinear polarization evolution as well as an extraction of nonlinear susceptibilities and time delays by fitting the polarization. We also carry out a decomposition of temporal and spatial changes in the electron density in power series with respect to the field amplitude. It helps to get insight into the origin of the nonlinear polarization in atomic scale.

Nonlinear electronic excitations in crystalline solids using meta-generalized gradient approximation and hybrid functional in time-dependent density functional theory.

We develop methods to calculate electron dynamics in crystalline solids in real-time time-dependent density functional theory employing exchange-correlation potentials which reproduce band gap energies of dielectrics; a meta-generalized gradient approximation was proposed by Tran and Blaha [Phys. Rev. Lett. 102, 226401 (2009)] (TBm-BJ) and a hybrid functional was proposed by Heyd, Scuseria, and Ernzerhof [J. Chem. Phys. 118, 8207 (2003)] (HSE). In time evolution calculations employing the TB-mBJ potential, we have found it necessary to adopt the predictor-corrector step for a stable time evolution. We have developed a method to evaluate electronic excitation energy without referring to the energy functional which is unknown for the TB-mBJ potential. For the HSE functional, we have developed a method for the operation of the Fock-like term in Fourier space to facilitate efficient use of massive parallel computers equipped with graphic processing units. We compare electronic excitations in silicon and germanium induced by femtosecond laser pulses using the TB-mBJ, HSE, and a simple local density approximation (LDA). At low laser intensities, electronic excitations are found to be sensitive to the band gap energy: they are close to each other using TB-mBJ and HSE and are much smaller in LDA. At high laser intensities close to the damage threshold, electronic excitation energies do not differ much among the three cases.

Ab initio simulation of electrical currents induced by ultrafast laser excitation of dielectric materials.

We theoretically investigate the generation of ultrafast currents in insulators induced by strong few-cycle laser pulses. Ab initio simulations based on time-dependent density functional theory give insight into the atomic-scale properties of the induced current signifying a femtosecond-scale insulator-metal transition. We observe the transition from nonlinear polarization currents during the laser pulse at low intensities to tunnelinglike excitation into the conduction band at higher laser intensities. At high intensities, the current persists after the conclusion of the laser pulse considered to be the precursor of the dielectric breakdown on the femtosecond scale. We show that the transferred charge sensitively depends on the orientation of the polarization axis relative to the crystal axis, suggesting that the induced charge separation reflects the anisotropic electronic structure. We find good agreement with very recent experimental data on the intensity and carrier-envelope phase dependence [A. Schiffrin et al., Nature (London) 493, 70 (2013).

High-Order Harmonic Generation in Solids: The Role of Intraband Transitions in Extreme Nonlinear Optics.

High-order harmonic generation (HHG) in gases is frequently used nowadays to produce attosecond pulses and coherent radiation in the visible-to-soft X-ray spectral range. HHG in solids is a natural extension of the idea of HHG in gases, and its first observation about ten years ago opened the door to investigations on attosecond electron dynamics in solids and the development of solid-state attosecond light sources. The common process in both types of HHG is nonlinear photocarrier generation, and thus, transitions between different bands (interband transitions) are always important for HHG. As well, in the case of solids, the transitions within a band (intraband transitions) also need to be considered, because efficient carrier acceleration is possible due to them. This Perspective focuses on experimental findings that show how intraband transitions can be controlled because such an understanding will be essential in the development of unique optoelectronics that can operate at petahertz frequencies.

Controlling Floquet states on ultrashort time scales.

The advent of ultrafast laser science offers the unique opportunity to combine Floquet engineering with extreme time resolution, further pushing the optical control of matter into the petahertz domain. However, what is the shortest driving pulse for which Floquet states can be realised remains an unsolved matter, thus limiting the application of Floquet theory to pulses composed by many optical cycles. Here we ionized Ne atoms with few-femtosecond pulses of selected time duration and show that a Floquet state can be observed already with a driving field that lasts for only 10 cycles. For shorter pulses, down to 2 cycles, the finite lifetime of the driven state can still be explained using an analytical model based on Floquet theory. By demonstrating that the amplitude and number of Floquet-like sidebands in the photoelectron spectrum can be controlled not only with the driving laser pulse intensity and frequency, but also by its duration, our results add a new lever to the toolbox of Floquet engineering.

Two-step Brillouin zone sampling for efficient computation of electron dynamics in solids.

We develop a numerical Brillouin-zone integration scheme for real-time propagation of electronic systems with time-dependent density functional theory. This scheme is based on the decomposition of a large simulation into a set of small independent simulations. The performance of the decomposition scheme is examined in both linear and nonlinear regimes by computing the linear optical properties of bulk silicon and high-order harmonic generation. The decomposition of a large simulation into a set of independent simulations can improve the efficiency of parallel computation by reducing communication and synchronization overhead and enhancing the portability of simulations across a relatively small cluster machine.

Conditional Wave Function Theory: A Unified Treatment of Molecular Structure and Nonadiabatic Dynamics.

We demonstrate that a conditional wave function theory enables a unified and efficient treatment of the equilibrium structure and nonadiabatic dynamics of correlated electron-ion systems. The conditional decomposition of the many-body wave function formally recasts the full interacting wave function of a closed system as a set of lower-dimensional (conditional) coupled "slices". We formulate a variational wave function ansatz based on a set of conditional wave function slices and demonstrate its accuracy by determining the structural and time-dependent response properties of the hydrogen molecule. We then extend this approach to include time-dependent conditional wave functions and address paradigmatic nonequilibrium processes including strong-field molecular ionization, laser-driven proton transfer, and nuclear quantum effects induced by a conical intersection. This work paves the road for the application of conditional wave function theory in equilibrium and out-of-equilibrium ab initio molecular simulations of finite and extended systems.

Simulating Vibronic Spectra without Born-Oppenheimer Surfaces.

We show how linear vibronic spectra in molecular systems can be simulated efficiently using first-principles approaches without relying on the explicit use of multiple Born-Oppenheimer potential energy surfaces. We demonstrate and analyze the performance of mean-field and beyond-mean-field dynamics techniques for the H2 molecule in one dimension, in the later case capturing the vibronic structure quite accurately, including quantum Franck-Condon effects. In a practical application of this methodology we simulate the absorption spectrum of benzene in full dimensionality using time-dependent density functional theory at the multitrajectory Ehrenfest level, finding good qualitative agreement with experiment and significant spectral reweighting compared to commonly used single-trajectory Ehrenfest dynamics. These results form the foundation for nonlinear spectral calculations and show promise for future application in capturing phenomena associated with vibronic coupling in more complex molecular and potentially condensed phase systems.

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 quasi-particles, called excitons. The manipulation of excitons in dielectrics has recently received great attention, with possible applications in different fields including optoelectronics and photonics. 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 MgF2 single crystals. Furthermore, we access absolute phase delays, which allow for an unambiguous comparison with theoretical calculations. 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. Further theoretical investigation reveals a link with the exciton sub-femtosecond nanometric motion and allows us to envision a new route to control exciton dynamics in the close-to-petahertz regime.