Umberto De Giovannini

Real-space grids and the Octopus code as tools for the development of new simulation approaches for electronic systems

Real-space grids are a powerful alternative for the simulation of electronic systems. One of the main advantages of the approach is the flexibility and simplicity of working directly in real space where the different fields are discretized on a grid, combined with competitive numerical performance and great potential for parallelization. These properties constitute a great advantage at the time of implementing and testing new physical models. Based on our experience with the Octopus code, in this article we discuss how the real-space approach has allowed for the recent development of new ideas for the simulation of electronic systems. Among these applications are approaches to calculate response properties, modeling of photoemission, optimal control of quantum systems, simulation of plasmonic systems, and the exact solution of the Schrödinger equation for low-dimensionality systems.

Creating stable Floquet–Weyl semimetals by laser-driving of 3D Dirac materials

Tuning and stabilizing topological states, such as Weyl semimetals, Dirac semimetals or topological insulators, is emerging as one of the major topics in materials science. Periodic driving of many-body systems offers a platform to design Floquet states of matter with tunable electronic properties on ultrafast timescales. Here we show by first principles calculations how femtosecond laser pulses with circularly polarized light can be used to switch between Weyl semimetal, Dirac semimetal and topological insulator states in a prototypical three-dimensional (3D) Dirac material, Na3Bi. Our findings are general and apply to any 3D Dirac semimetal. We discuss the concept of time-dependent bands and steering of Floquet–Weyl points and demonstrate how light can enhance topological protection against lattice perturbations. This work has potential practical implications for the ultrafast switching of materials properties, such as optical band gaps or anomalous magnetoresistance.

Simulating Pump–Probe Photoelectron and Absorption Spectroscopy on the Attosecond Timescale with Time‐Dependent Density Functional Theory

Molecular absorption and photoelectron spectra can be efficiently predicted with real-time time-dependent density functional theory. We show herein how these techniques can be easily extended to study time-resolved pump-probe experiments, in which a system response (absorption or electron emission) to a probe pulse is measured in an excited state. This simulation tool helps with the interpretation of fast-evolving attosecond time-resolved spectroscopic experiments, in which electronic motion must be followed at its natural timescale. We show how the extra degrees of freedom (pump-pulse duration, intensity, frequency, and time delay), which are absent in a conventional steady-state experiment, provide additional information about electronic structure and dynamics that improve characterization of a system. As an extension of this approach, time-dependent 2D spectroscopy can also be simulated, in principle, for large-scale structures and extended systems.

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

We report the spin-selective optical excitation of carriers in inversion-symmetric bulk samples of the transition metal dichalcogenide (TMDC) WSe_{2}. Employing time- and angle-resolved photoelectron spectroscopy (trARPES) and complementary time-dependent density functional theory (TDDFT), we observe spin-, valley-, and layer-polarized excited state populations upon excitation with circularly polarized pump pulses, followed by ultrafast (<100  fs) scattering of carriers towards the global minimum of the conduction band. TDDFT reveals the character of the conduction band, into which electrons are initially excited, to be two-dimensional and localized within individual layers, whereas at the minimum of the conduction band, states have a three-dimensional character, facilitating interlayer charge transfer. These results establish the optical control of coupled spin-, valley-, and layer-polarized states in centrosymmetric materials with locally broken symmetries and suggest the suitability of TMDC multilayer and heterostructure materials for valleytronic and spintronic device concepts.

Modeling electron dynamics coupled to continuum states in finite volumes with absorbing boundaries

Absorbing boundaries are frequently employed in real-time propagation of the Schrödinger equation to remove spurious reflections and efficiently emulate outgoing boundary conditions. These conditions are a fundamental ingredient for the calculation of observables involving infinitely extended continuum states in finite volumes. In the literature, several boundary absorbers have been proposed. They mostly fall into three main families: mask function absorbers, complex absorbing potentials, and exterior complex-scaled potentials. To date none of the proposed absorbers is perfect, and all present a certain degree of reflections. Characterization of such reflections is thus a critical task with strong implications for time-dependent simulations of atoms and molecules. We introduce a method to evaluate the reflection properties of a given absorber and present a comparison of selected samples for each family of absorbers. Further, we discuss the connections between members of each family and show how the same reflection curves can be obtained with very different absorption schemes.

Monitoring Electron-Photon Dressing in WSe2

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 pumpprobe 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 trARPES 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 spinresolved PES can be used to predict asymmetric dichroic response linked to the valley selective optical excitations in monolayer transition metal dichalcogenides (TMDs). 1 ar X iv :1 60 9. 03 21 8v 2 [ co nd -m at .m tr lsc i] 1 3 Se p 20 16

A First-Principles Time-Dependent Density Functional Theory Framework for Spin and Time-Resolved Angular-Resolved Photoelectron Spectroscopy 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 analyzing surface. The resulting approach, that we named t-SURFFP, allows us 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, monolayer, and bilayer WSe2, and hexagonal BN (hBN) under different laser configurations.Angel Rubio‡ University of the Basque Country UPV/EHU, Nano-Bio Spectroscopy Group, Avenida de Tolosa 72, 20018 San Sebastian, Spain and Max Planck Institute for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany. Center for Free-Electron Laser Science and Department of Physics, University of Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany (Dated: September 13, 2016)

Stark Ionization of Atoms and Molecules within Density Functional Resonance Theory

We show that the energetics and lifetimes of resonances of finite systems under an external electric field can be captured by Kohn–Sham density functional theory (DFT) within the formalism of uniform complex scaling. Properties of resonances are calculated self-consistently in terms of complex densities, potentials, and wave functions using adapted versions of the known algorithms from DFT. We illustrate this new formalism by calculating ionization rates using the complex-scaled local density approximation and exact exchange. We consider a variety of atoms (H, He, Li, and Be) as well as the H2 molecule. Extensions are briefly discussed.

Efficient and accurate modeling of electron photoemission in nanostructures with TDDFT

Abstract We derive and extend the time-dependent surface-flux method introduced in [L. Tao, A. Scrinzi, New J. Phys. 14, 013021 (2012)] within a time-dependent density-functional theory (TDDFT) formalism and use it to calculate photoelectron spectra and angular distributions of atoms and molecules when excited by laser pulses. We present other, existing computational TDDFT methods that are suitable for the calculation of electron emission in compact spatial regions, and compare their results. We illustrate the performance of the new method by simulating strong-field ionization of C60 fullerene and discuss final state effects in the orbital reconstruction of planar organic molecules.

Phonon-driven spin-Floquet magneto-valleytronics in MoS2

Two-dimensional materials equipped with strong spin–orbit coupling can display novel electronic, spintronic, and topological properties originating from the breaking of time or inversion symmetry. A lot of interest has focused on the valley degrees of freedom that can be used to encode binary information. By performing ab initio time-dependent density functional simulation on MoS2, here we show that the spin is not only locked to the valley momenta but strongly coupled to the optical E″ phonon that lifts the lattice mirror symmetry. Once the phonon is pumped so as to break time-reversal symmetry, the resulting Floquet spectra of the phonon-dressed spins carry a net out-of-plane magnetization (≈0.024μB for single-phonon quantum) even though the original system is non-magnetic. This dichroic magnetic response of the valley states is general for all 2H semiconducting transition-metal dichalcogenides and can be probed and controlled by infrared coherent laser excitation.In 2H semiconducting transition-metal dichalcogenides the valley-selective excitation has been achieved with circularly polarized photons. Here, the authors show that circularly polarized phonons produce a valley-dependent dynamic spin state as a result of strong spin-phonon coupling.The valley degree of freedom and the possibility of spin-valley coupling of solid materials have attracted growing interest, and the relaxation dynamics of spin- and valley-polarized states has become an important focus of recent studies. In spin-orbit-coupled inversion-asymmetric two-dimensional materials, such as MoS_{2} it has been found that the spin randomization is characteristically faster than the time scales for inter- and intra-valley scatterings. In this study, we examined the ultrafast non-collinear spin dynamics of an electron valley in monolayer MoS_{2} by using real-time propagation time-dependent density functional theory. We found that the spin precession of an electron in the valley is sharply coupled with the lowest-lying optical phonon that release the in-plane mirror symmetry. This indicates that the spin randomization of MoS_{2} is mainly caused by spin-phonon interaction. We further suggest that flipping of spins in a spin-orbit-coupled system can be achieved by the control over phonons.