Time-resolved photoemission of infinitely periodic atomic arrangements: correlation-dressed excited states of solids

A theory of the time-resolved photoemission spectroscopy (TRPES) is developed, which enables to explore the real-time electron dynamics of infinitely periodic crystalline solids. In the strongly correlated electron systems NiO and CuO, the early-stage dynamics of the valence band edge are found to be sharply contrasted between those in the spectra of TRPES. This provides a new dynamical insight to the Zaanen–Sawatzky–Allen (ZSA) classification scheme of correlated insulators and makes us assert that NiO dynamically behaves like the Mott–Hubbard insulator (MHI) and CuO like the charge transfer insulator (CTI). In the two-dimensional carbon layer graphene, the real-time electron dynamics of quantum-phase-dressed excited states, i.e., due to the Berry phase and the pseudospin correlation, are investigated in an unprecedented way through the time-resolved angle-resolved photoemission spectroscopy (TR-ARPES). In particular, the dephasing dynamics of optically doped electrons and holes in the massless Dirac band, accompanying a field-induced gliding of the Dirac cone, are discovered.


INTRODUCTION
Ultrafast phenomena of the optically driven transient changes in charge, spin, lattice, and orbital degrees of freedom and their crosslinks provide benchmarks for the photoinduced dynamics of a given material. In particular, the dynamical instabilities induced by the cooperative interplays among those degrees of freedom would be classified as the photoinduced phase transitions (PIPTs) 1,2 , which are in fact observed in a wide range of materials. An optical doping through the pulse laser pumping triggers the photoinduced deployments of simultaneously introduced electrons and holes differently from the chemical doping, which carries all the essence of the photoinduced dynamics and related PIPTs. Hence, instead of conventional approaches based on the time-resolved macroscopic electro-optic measurements [3][4][5] , one may attempt to explore the photoinduced ultrafast dynamics through an explicit tracking of the creation, development, and destruction of optically doped electrons and holes. For the purpose, the most promising framework would be the time-resolved photoemission spectroscopy (TRPES). TRPES, a kind of the pump-probe spectroscopy, takes a snap shot of the photoemission using the delayed probe pulse under the dynamics initiated by the pump pulse and monitors the time evolution in the electronic structures of a material out of equilibrium in a manifold of energy and time 6-16 , i.e., conjugate variables in the Heisenberg uncertainty principle.
Insulator-metal transition (IMT) is one of the most exciting phenomena in condensed matter physics 17 . One of the possible routes to drive IMT is to simultaneously introduce electrons and holes by irradiating an insulator, i.e., the photoinduced IMT (PIIMT) 1,2 . It is TRPES that is a powerful tool to probe PIIMT in a real-time domain through a direct observation of temporal changes of the electronic structure caused by optically doped electrons and holes [6][7][8][9][10][11]13,16,18,19 . In the experimental aspect, TRPES requires high standards of quantum optics and electronics so that it has only recently been implemented actively. Meanwhile, in the theoretical side, TRPES has been less explored due to inherent difficulties in describing the real-time dynamics. Although the nonequilibrium Green's function method was proposed, an actual treatment of the method has been still in a sort of model calculation [20][21][22] . Instead, for a theoretical description of TRPES for real solids (i.e., infinitely periodic solids beyond the model consideration), we attempt to reformulate the theory of TRPES based on the many-body time-dependent Schrödinger equation 23 such that the first-principles calculation can be applied. In fact, no established theoretical formulation is available yet especially for infinitely periodic crystalline solids in the firstprinciples framework 24,25 . For the numerical encoding of TRPES, a complicated engagement of many electrons system in the photoemission process should be managed at an appropriate level of approximation.
We develop an original theory of TRPES within the three-step model of photoemission and address optically induced transient changes in the electronic structure of transition metal oxides NiO and CuO in the angle-integrated mode and two-dimensional gapless semiconductor graphene in the angle-resolved mode. First, we choose NiO and CuO due to the following reasons. (i) Although NiO and CuO are strongly correlated electron systems, the modified Becke-Johnson local density approximation (mBJ-LDA) 26 potential is found to describe their ground states quite successfully. This gives a possibility to manage the electron dynamics of the two materials in a new formalism of TRPES with time-dependent density functional theory (TDDFT) based on mBJ-LDA. (2) NiO and CuO have been studied for a long time as typical transition metal oxides but rarely in respects of the real-time electron dynamics. Hence, critical scientific insights may be expected in their real-time electron dynamics. Next, we choose graphene because it may be a very good test bed for a new formulation of TRPES based on the first-principles calculation and also a very interesting system especially in the dynamic respects. Now we examine the time evolution of the valence band edges of NiO and CuO and provide a new dynamical vision to the Zaanen-Sawatzky-Allen (ZSA) scheme 27 . It is found that a strong optical pump results in a reconstruction of valence band edge to create incoherent in-gap states in NiO, whereas a simple weight reduction in the valence band edge of CuO. With this finding, we conclude that the dynamical response of NiO is like the Mott-Hubbard insulator (MHI) and that of CuO like the charge transfer insulator (CTI). For a two-dimensional hexagonal lattice of carbon atoms graphene, we investigate the dynamics of kdependent excited states, which are dressed with the dynamical quantum phase determined by the Berry phase and the pseudospin correlation depending on the polarizations of pump and probe pulses. Those excited states in the massless Dirac band, incorporating a gliding of the Dirac cone according to an acceleration of electrons along the easy direction, are further explored by tracing optically doped electrons and holes in the real-time domain.

Theoretical formulation of TRPES
In the present formulation, TRPES is formally based on the threestep model of the photoemission: the photoelectron (i) creation, (ii) propagation and penetration through the surface, and (iii) escape to the vacuum, and additionally the free electron approximation for the final state 28 . Inspired from the first-order time-dependent perturbation theory, we propose the photoemission amplitude To build up a theoretical formulation of TRPES from Eq. (1), we employ the first-principles TDDFT calculation implemented in ELK (see http://elk.sourceforge.net). We obtain the time-evolving Kohn-Sham (KS) wave function |n, q;τ〉 from the time-dependent KS equation for the occupied electron with the momentum q and the band index n, whose time evolution is governed entirely by an optical pump A pump (τ) through TDDFT. |k(τ)〉 is a time-evolved final state of a plane wave to be overlapped with |n, q;τ〉 through p, where p is the linear momentum operator. The equation provides a powerful chance for directly calculating TRPES within the sudden approximation even though the KS wave function is in an infinitely periodic atomic system. Now c n;q k ðτ d Þ stands for a probability amplitude for finding the photoelectron from |n, q;τ〉, which can then be computed through the first-order timedependent perturbation theory due to the probe pulse A probe (τ − τ d ), i.e., extreme ultraviolet pulse, with the time delay τ d with respect to the pump pulse A pump (τ), In Eq. (2), we simply write |k(τ)〉 as k j ie Ài R dτ k 2 2 þEF À Á and insert a completeness identity for the ground KS states. This increases an efficiency for calculating spectra of TRPES since only the timedependent momentum matrices 〈n′, q|p|n, q;τ′〉 are required for the calculation process, that is, |n, q;τ′〉 does not have to be expanded into |k(τ) basis for every time step. Vector potentials of pump and probe pulses are given by Here Θ(x) is the Heaviside step function and τ is adopted as 4.9 fs for the pump pulse and 2.4 fs for the probe pulse. Each pulse has then the duration of 2τ. Furthermore, A 0 probe A 0 pump À1 ( 1 is assumed. Hence, the spectrum of TRPES gets to be J k τ d ð Þ ¼ P n;q c n;q k τ d ð Þ 2 . By using Eq. (2), we demonstrate photon-dressed Floquet states of Ar atomic orbitals in Supplementary Fig. 1, which guarantees that Eq. (2) can provide proper spectra of nonequilibrium states. Based on the method, in the present study, we mainly concentrate on the angle-integrated spectra of TRPES for the strongly correlated systems NiO and CuO, whereas on the angle-resolved spectra of TRPES by keeping the wave vectors k for the two-dimensional semiconductor graphene. In our approach of the three-step model, step (ii) is actually ignored by neglecting the surface scattering, which confirms that the TDDFT calculation, i.e., the time evolution of |n, q;τ〉, could be carried out with the full translational symmetry of the material. It is noted that the formulation of Eqs. (1) and (2) does not include the higher-order electron emission process, such as the Auger process.
Correlation induced band mixing in the excited states of NiO and CuO Over the past few decades, one of the most intensively discussed issues at the heart of condensed matter physics is the strongly correlated electron system. Historically, the fundamental test bed of correlated solid electrons has been the 3d transition metal oxide, which includes NiO that has then been long considered to be MHI. In MHI, a band gap opened between the upper Hubbard band (UHB) and the lower Hubbard band (LHB) is scaled by the d-d Coulomb correlation U 29 , which is described by the Hubbard model 30 . Later, however, the X-ray photoemission spectroscopy 31,32 and the cluster calculation 33 of NiO have suggested that, differently from the previous understanding, ligand p orbitals from oxygen atoms are strongly associated with the valence edge and the band gap is determined between the p-like ligand band and the d-like UHB. This indicates another kind of insulator, CTI, according to the ZSA scheme 27 , which classifies transition metal oxides into MHI and CTI according to a characteristic of the valence band edge, i.e., whether the valence band edge is made of d or p orbital. Later, there have followed more investigations that NiO should be possibly classified as an intermediate one rather than either MHI or CTI because the p-d hybridization nontrivially affects the electronic structure [34][35][36] . One may expect better understanding of the insulating nature of NiO and solve the aforementioned long-time controversy through a different angle of view, for instance, the real-time response of the electronic structure. Unlike NiO, CuO has been regarded as CTI due to a relatively simple electronic structure with the Coulomb correlation U larger than the charge transfer strength 37 . Here we demonstrate photoinduced electron dynamics of NiO and CuO through an angle integrated mode of TRPES, i.e., Moreover, for the calculation, we adopt the Tran-Blaha modified Becke-Johnson potential and the Perdew-Wang LDA (TB-mBJLDA) as the exchange and the correlation potential 26 . TB-mBJLDA is shown to predict correctly band gaps and magnetic moments for NiO and CuO, whose reliability would be of the same order as the hybrid functionals or the GW method (Supplementary Table 1).
Density of states (DOS) of NiO and CuO are displayed in Fig. 1a, b. Comparing between those, p-and d-band characters are mixed at the valence band edges of both materials, but the p character is found relatively stronger in CuO than in NiO.
to a band gap of the material. ω probe is taken as 25 eV. In Fig. 1c, d, we provide the spectra of TRPES at time delays τ d = −9.8 and 13.4 fs for NiO and CuO, respectively. τ d is the delay time between the two pulses. Let us be aware that spectra at τ d = −9.8 fs (i.e., before the pump pulse arrives) corresponds to the conventional PES, which directly reflect the occupied DOS in the ground state. However, spectra at τ d = 13.4 fs are indeed interesting in that snapshots of dynamical responses of the electronic structure under the optical pumping are taken. It is of great importance to note that both valence and conduction band edges are reconstructed to shift inwardly so that in-gap states are created in NiO of Fig. 1c. The resultant gap is found much smaller than the original gap E 0 g . In CuO of Fig. 1d, however, such a shift of the valence band edge is not shown. Instead, just a slight decrease in the weight near the edge is observed. Besides, the reconstructed edge of the conduction band looks likely spreading over the gap. It is remarkable to find that, although p and d characters at the valence band edges are similarly mixed in both materials in the ground state, the edge behaviors are obviously different from each other in TRPES. Schematics of band reconstructions of NiO and CuO are sketched in Fig. 1e, f.
Real-time dynamics are demonstrated through the spectra of TRPES at τ d = −4.9, −2.4, 0, 2.4, 4.9 fs shown (Fig. 2a, b). We keep track of a reconstruction of the valence and conduction band edges for NiO and CuO, which are marked by V and C, respectively. According to Fig. 2c, as time goes, both edges V and C of NiO are clearly shown to proceed to the inside of the gap and make in-gap states, as previously discussed in Fig. 1. Experiments of 1T-TaS 2 6 , TbTe 3 7 , 1T-TiSe2 8 , VO 2 11 , and cuprates 16 have actually reported inward shifts of LHBs, which result in a formation of forbidden mid-gap states. In the sense, a reduction of the gap could be regarded as intrinsic in correlated d bands under the optical pumping. In Fig. 2d, on the other hand, V of CuO in fact shows no noticeable shift even though C moves similarly to NiO. Moreover, the two-dimensional plots of spectra of TRPES clearly identify the retardation time τ ret to reach a reconstruction of the electronic structure (Fig. 2e, f). Hence, those dynamics of the band edges make us stride toward that NiO is dynamically MHI-like and CuO is dynamically CTI-like.
Different behaviors in the valence band edge between NiO and CuO are then further analyzed through the orbital-resolved spectra of TRPES ( Supplementary Fig. 2). Interestingly, the orbital-resolved spectra tell us that, due to the probable p-d hybridization, the oxygen p band of NiO is d-like and the Cu d band of CuO is p-like. In addition, exact diagonalization studies of the many-body model Hamiltonian evidently support our conclusion. We now propose a model Hamiltonian H comprising H 0 and H 1 , where H 0 describes the one-dimensional crystalline solid and H 1 the dynamical process of TRPES Here c þ iσ (c iσ ) is an electron creation (annihilation) operator at the ith site with its spin σ and similarly c þ kσ or c kσ a photoelectron operator with its energy ε k = k 2 /2 and spin σ. n iσ (=c þ iσ c iσ ) is the electron number operator. t implies the nearest neighbor hopping and U the on-site Coulomb repulsion. μ is the chemical potential  Table S1), respectively. e, f Schematics of the band reconstruction induced by the optical pump in NiO and CuO. In the calculation of the angle-integrated spectra, a slow variation of the matrix element to remain after completing an angle integration is neglected. and Δ the ionic potential. Pump and probe terms of Here Θ(x) is the Heaviside step function and τ is adopted as 4.9 and 2.4 fs for the pump and probe pulses, respectively. Further, V 0 probe =V 0 pump ( 1 and V 0 pump ¼ 0:5 is assumed. In Fig. 3a-d, Δ = 0 and U/t = 5 is taken and H 0 then falls to the N-site (N = 6) Hubbard model. Figure 3a-c provide the spectral functions of N (i.e., half-filling), N + 1, and N − 1 electrons, respectively. At half-filling, the spectral function shows a Hubbard gap E 0 g (=3.6 eV) opened between the UHB and the LHB (Fig. 3a). This defines a correlated insulator corresponding to the prototype of the MHI. The (N + 1)-electron spectral function shows that, when adding an electron to the N-electron system, the electron does not jump simply into the N-electron conduction band but into a reconstructed state (marked by "e") just below the conduction band (Fig. 3b). Similarly, when adding a hole to the N-electron state, the hole does not stay in the N-electron valence band but instead prefers a reconstructed state (marked by "h") just above the valence band (Fig. 3c). TRPES incorporates an optical doping (i.e., a simultaneous creation of electron and hole due to the optical pumping) and naturally involves both features of (N + 1)-and (N − 1)-electron spectral functions, i.e., both of "e" and "h", which clearly explains inward shifts of the conduction and valence edges into the gap as schematically drawn in Fig. 2e of the main text. Figure 3d gives the spectra of TRPES and PES, whose calculation can be done by extending the many-body Hilbert space to include the photoelectron under the total Hamiltonian H (=H 0 + H 1 ) 23 . Spectra of TRPES show well the conjugated feature of (N + 1)-and (N − 1)-electron spectra, which confirms the discussion above. In a sense, on the other hand, since nonexcited electron states (i.e., N-electron states) still remain, TRPES may be approximately viewed as an addition of (N + 1)-and (N − 1)-electron spectral functions on an existence of N-electron function. Besides, it is seen that the configuration interaction between (N + 1)-and (N − 1)-electron states leading to a formation of exciton would be more or less negligible in a small size Hubbard model. Figure 3e-h treat H 0 with U = 0 and 2Δ/t = 5, which describes an ordinary semiconducting material. As illustrated in Fig. 3e-g, it is seen that, when adding either an electron or a hole to the N-electron system, the band edges behave rigidly Quantum-phase-dressed excited states of graphene During the past decades, the angle-resolved photoemission spectroscopy (ARPES) has been the most direct and powerful method for an investigation of two-dimensional solids in their reciprocal lattices. An extension of ARPES along the time dimension is the time-resolved ARPES (TR-ARPES). Because TR-ARPES can track the q-dependent electronic structure in the manifold of energy and time and provide a wealth of real-time quantum information on the electronic dynamics, it should be a powerful tool especially for the two-dimensional semiconducting systems. Here we demonstrate the photoinduced real-time electron dynamics of the Dirac cone of graphene through an angle-resolved mode of TRPES, i.e., J k τ d ð Þ ¼ P n;q c n;q k τ d ð Þ 2 . In the calculation, the LDA is adopted for the exchange and correlation potential.
ARPES has a merit in observing the quantum phase absorbed in the matrix element, for instance, the Berry phase 38 and nonequilibrium carrier dynamics 39 for graphene. Taking ω probe = 35 eV, we provide the momentum-energy distribution curves (MEDCs) of ARPES, i.e., the band mapping, of graphene (Fig. 4a, b) and electron-doped (n-type) graphene (Fig. 4c, d) along the K′ − Γ − K direction and theq y -direction through the K point. In MEDCs, not only the massless Dirac bands cutting the Fermi level but also low-lying 2s-like bands are shown (Fig. 4a-d). In Fig. 4e, f, schematics of two Dirac cones centered at K′ and K points with Matrix elements for ARPES, which contain the information on the Berry phase 38 , make a characteristic selection for the Dirac cones such that only half (blue-colored conical section) of the lower or the upper Dirac cone along theq x -direction through the K′ or K point survives for thex-polarized probe pulse (Fig. 4a, c, e, f). Different phases of p z orbitals from two carbon atoms characteristically cancel the matrix element for the photoelectron state |k〉 40,41 depending on whether it is emitted from the lower or the upper Dirac cone. Along theq y -direction through the K point, however, the full Dirac cones are signaled in Fig. 4b, d. Moreover, taking ω pump = 3 eV and E 0 pump ¼ ω pump A 0 pump ¼ 0:1 VÅ À1 , we provide MEDCs of TR-ARPES of graphene at a fixed delay time τ d = 7 fs along theq y -direction through the K point for thex-polarized pump (Fig. 4g) and the K′ − Γ − K direction for the b y-polarized pump (Fig. 4i) and the constant energy mappings (at E B = 1.5 eV (>0), i.e., in the upper Dirac cone) in the crystal momentum space for thex-polarized (Fig. 4h) and the b y-polarized pump (Fig. 4j). TR-ARPES captures the physics of dynamical quantum phases of the excited electron or hole states through the pseudospin correlation. Pseudospins of the Dirac cones are definitely coupled with the optical transition (Fig. 4k, l) 42 .
Obviously, the pseudospin correlation should be considered together with the selection due to the matrix elements in determining the spectra of TR-ARPES. Consequently, vanishing nodal lines result along q y = 0 and q x = K around the K point for thex-polarized and b y-polarized pumps, respectively, in the constant energy mappings of TR-ARPES (Fig. 4h, j). This can be directly understood through a compact form of the model Hamiltonian for the Dirac cone, H 0 = v F q ⋅ σ, where q is the electron wave vector, σ the Pauli matrix representing the pseudospin, and v F the Fermi velocity, is cast to  . e, f Half of the Dirac cone for graphene (e) and halves of the Dirac cones for the n-type graphene (f) with nonzero matrix elements (blue-colored conical surfaces), which contain the information on the Berry phase, are captured by ARPES. Red arrows denote the pseudospin configuration. g-j Spectra of TR-ARPES of graphene at τ d = 7 fs for thex-polarized (g, h) and the b y-polarized pump (i, j), i.e., MEDCs (g, i) and the constant energy mapping at E B = 1.5 eV (h, j). Arrows in g, i indicate the weight transfers to the conduction bands due to the optical pumping. k, l Pesudospin selection under the interband transition due to thex-polarized (k) and the b y-polarized pump (l).
at q y = 0 on the Dirac cone because σ x should be reserved in the optical transition (Fig. 4k). A case for the b y-polarized pump could be understood in the same fashion (Fig. 4l). In Fig. 5, the time-dependent constant energy mappings of TR-ARPES of graphene focusing on the K-centered Dirac cone are displayed at τ d = −7, −2.4, 0, 2.4, 7 fs. First, it is interesting to compare the spectra of ARPES of the n-type graphene with the Dirac point near −1.5 eV (Fig. 5a) with TR-ARPES of graphene from thex-polarized (Fig. 5b) and the b y-polarized pump (Fig. 5c). Spectra of TR-ARPES at E B < 0 are from the lower Dirac cone with optically doped holes at ω pump = 3 eV are found to be almost same as ARPES of the n-type graphene at E B < −1.5 eV, which indicates that the valence band is in fact rigid with respect to the optical doping but a tiny reduction in the spectral weight. As previously stated, the anisotropy selecting a half of the Dirac cone with q x < K (q x > K) in the spectra of TR-ARPES at E B < 0 (E B > 0) under thex-polarized probe pulse comes from the Berry phase involved in the matrix element. Meanwhile, the spectra at E B > 0 are from the upper Dirac cone with optically doped electrons, i.e., electrons excited to the conduction band by the pump pulse, which are inevitably carried with the pseudospin correlation. As the time delay increases, nodal lines become clear along q y = 0 and q x = K for thex-polarized (Fig. 5b) and the b y-polarized pump (Fig. 5c), respectively, which results from the dynamical quantum phase due to the pseudospin correlation. Spectra of ARPES of the n-type graphene would be more or less reproduced by making a summation of both of those of TR-ARPES from thex-polarized and b y-polarized pumps.
Needless to say, the dynamical quantum phase should be also incorporated in a part of the optically doped holes. We attain the difference spectra as J k (τ 1 ) − J k (τ 2 ) between the two time delays τ 1 and τ 2 (τ 1 > τ 2 ) in Fig. 6. For a case of thex-polarized pump, the spectral differences of doped electrons (red) and holes (blue) in the constant energy mapping at E B = 1.5 eV and E B = −1.5 eV, respectively, are reasonably dispersed according to the spectra of TR-ARPES. Notably, those of doped holes have a nodal line along q y = 0 just like electrons 43 . Both of doped electrons and holes undergo dephasing monotonically as time goes and get to reside on well-defined Dirac cones through the energy conservation (Fig.  6b). Thus a silhouette of the full Dirac cone with the nodal line is acquired by adding the spectral differences at E B = 1.5 eV (electrons) and E B = −1.5 eV (holes) in J k (7 fs) − J k (2.4 fs) (Fig. 6c). For the b y-polarized pump, doped electrons and hole-like negative-E B signals (i.e., instead of "doped holes," due to a reason to be cleared soon) are definitely unbalanced in their spectral densities ( Fig. 6d) even if, for thex-polarized pump, the charge neutrality of the doped carriers is reasonably maintained. Furthermore, the hole-like negative-E B signals are found not to have a nodal line along q x = K, either. These indicate that such signals at the b y-polarized pump should not be due solely to the optical doping. Let us note that, in our geometry, the b y-polarized pump is along the armchair direction of graphene so that it could efficiently accelerate p z -orbitals belonging to a pair of neighboring carbon atoms along the b y-direction. In this consideration, the b y-polarized pump breaks the C 3 rotational symmetry of the system and gives rise to a real-time moving of the Dirac cone 44 . In fact, a gliding shift of the Dirac cone outward the Brillouin zone along thê x-direction is found for the static field along the b y-direction ( Fig.  6e and Supplementary Fig. 8). Therefore, broad hole-like negative-E B signals around the edge of the hexagon (Fig. 6d) indicate a decrease (blue) and an increase (red) in the electronic states in the momentum space at E B = −1.5 eV according to a continuous fieldinduced moving of the Dirac cone beyond the K-point along thê x-direction (Fig. 6f), which would be clear in the difference spectra of MEDCs in Fig. 7a-c (also see Supplementary Fig. 5).
Difference spectra of MEDCs from TR-ARPES (Fig. 7a-c) show the overall dynamics of Dirac cones in a clearer fashion. At the first stage under the b y-polarized pump pulse (Fig. 7a), together with the optical doping between conduction and valence bands, Dirac cones are found to glide toward the outward direction of K′ and K. In the intermediate time delays (Fig. 7b), Dirac cones seem to remain temporarily at the displaced positions. At the final stage (Fig. 7c), Dirac cones return to their original positions at K′ and K. A gliding shift of the Dirac cone is discussed previously in relation to a breaking of the C 3 rotational symmetry. This can be reproduced by a simple tight-binding Hamiltonian for graphene given by where Δ k ¼ Àt 0 ðe ikÁa1 þ e ikÁa2 þ e ikÁa3 Þ and a n is the unit vector to the nearest neighbor C atom and t 0 (=3 eV) is the hopping (or hybridization) parameter. Here, assuming a 3 ∥ b y, a strong pump field along the b y-direction would immediately render a hybridization strength between C atoms along the a 3 -direction to increase. That is, in the new Hamiltonian, Δ k may be replaced by Δ k ¼ Àt 0 e ikÁa1 þ e ikÁa2 À Á À γt 0 e ikÁa3 , where γ will be regarded as a symmetry breaking parameter in the hybridization between the nearest neighbor C atoms along b y-direction, that is, at γ = 1, graphene recovers the C 3 rotational symmetry, i.e., an equilibrium state, whereas at γ > 1, the symmetry is broken under the b y-polarized pump pulse. In Fig. 7d, e, γ is found to linearly result in outward glidings of Dirac cones of K′ and K, respectively, and γ~1 .02 to induce δκ consistent with TR-ARPES of Figs 6e, f and 7a-c.
The present calculation has been performed based on the single-body equation with the Born-Oppenheimer approximation, Fig. 7 Difference spectra of MEDCs of TR-ARPES and gliding shifts of Dirac cones. a-c Different spectra of MEDCs, i.e., J k (τ 1 ) − J k (τ 2 ) (τ 1 > τ 2 ) for the b y-polarized pump. Schematics on right panels represent time-dependent movings of Dirac cones at K′ and K. Black and red arrows indicate the field-induced moving of Dirac cones and the interband optical transitions, respectively. d, e Tight-binding model calculation of graphene. d A model calculation for the outward gliding shift of the Dirac cone. Red and Blue band structures represent symmetry conserved (i.e., γ = 1; equilibrium) and symmetry broken (i.e., γ = 1.1; nonequilibrium) states, respectively. Here, γ is a symmetry breaking parameter. The inset displays schematics of the gliding shifts of Dirac cones. e Estimation of δκ (Fig. 6f) with respect to values of γ. b Difference spectra of constant energy mappings of TR-ARPES as J k (τ 1 ) − J k (τ 2 ) between the two time delays τ 1 and τ 2 (τ 1 > τ 2 ) for thex-polarized pump. c Silhouette of the full Dirac cone with the nodal line of k y = 0 from an addition of difference spectra at E B = 1.5 eV and E B = −1.5 eV for thex-polarized pump. d Difference spectra for the b y-polarized pump. e Difference spectra between with and without the static field E field b y at E B = 0 eV ( Supplementary Fig. 8). Dirac cone of graphene undergoes an outward shift along thex-direction as well as a gliding shift along the static field (∝b y). For the b y-polarized optical pump, only the outward shift along thex-direction would remain in an approximate sense. f Field-induced moving (δκ) of the Dirac cone with respect to the time delay at the b y-polarized pump.
but the adiabatic Born-Oppenheimer approximation is recognized not to be valid in graphene 45 . Here it is stressed that a gliding shift of the Dirac cone is originated by a breaking of the sixfold rotation symmetry of the sp 3 hybridization due to the optical pumping, not an instantaneous change of atomic positions. Therefore, we still expect that the gliding would be observed at least during the optical pumping, i.e., the external optical perturbation may overwhelm the phonon-related underlying process. What would happen after the optical pump field has passed will be an interesting question due to the nonadiabatic electron-phonon coupling.

DISCUSSION
In the first-principles framework, we have explored an original theory of TRPES or TR-ARPES for infinitely periodic solids and examined the real-time electron spectra focusing on the excited valence band edges of strongly correlated electron systems and the quantum-phase-dressed excited states of a two-dimensional semiconductor graphene. In strongly correlated systems, dynamics of valence band edges under a strong optical pump have been found to be very different between NiO and CuO. This could provide a critical insight to the dynamical aspect of the ZSA scheme, from which we claim that NiO dynamically behaves like the Mott-Hubbard system and CuO like the charge transfer system. In graphene, a single layer of carbon atoms, excited states dressed with the Berry curvature, and the pseudospin correlation have been thoroughly tracked in the angle-resolved mode. Particularly, the dephasing dynamics of optically doped electrons and holes on the Dirac cone, more intricately on a field-induced moving of the Dirac cone, have been discovered in the real-time domain.

First-principles calculations
Calculations for initial ground states were computed by using the allelectron linearized augmented planewave (APW) method implemented in ELK (see http://elk.sourceforge.net) code. First order of extra radial function is added to the local orbitals and APWs. The maximum length for the momentums multiplying averaged muffin-tin radius is fixed to be 8.0. The reciprocal lattice vector cutoff for the electron density and potential is 12 atomic unit. For the calculation of CuO and NiO, we performed the spinpolarized calculation and TB-mBJLDA DFT with sampling of the first Brillouin zone of 4 × 4 × 4. We confirmed the convergence of the time evolution until 8 × 8 × 8. For graphene, the LDA is adopted together with the first Brillouin zone sampling 27 × 27 × 1. The time-step interval for the time-evolution calculation is fixed to 0.1 atomic unit.