Abstract
We investigate the effect of broken inversion symmetry on the generation and control of ultrafast currents in a transparent dielectric (SiO2) by strong femtosecond optical laser pulses. Ab initio simulations based on time-dependent density functional theory predict ultrafast direct currents that can be viewed as a nonlinear photogalvanic effect. Most surprisingly, the direction of the current undergoes a sudden reversal above a critical threshold value of laser intensity of about W cm−2. We trace this switching to the transition from nonlinear polarisation currents to the tunnelling excitation regime. The latter is found to be sensitive to the relative orientation between laser polarisation and chemical bonds. We demonstrate control of the ultrafast currents by the time delay between two laser pulses. While two temporally separated laser pulses lead to currents along one direction their temporal overlap can reverse the current. We find the ultrafast current control by the nonlinear photogalvanic effect to be remarkably robust and insensitive to the laser-pulse shape and the carrier-envelope phase.

Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
1. Introduction
In the last decade, ultrafast few-cycle laser pulses with well-defined carrier-envelope phase (CEP) have become available providing novel opportunities to explore the ultrafast and nonlinear response of matter to strong optical fields. The study of the induced electronic motion and of the highly nonlinear optical response have focussed on rare gas atoms [1], molecules [2], and, more recently, on nanostructures, surfaces, and bulk matter [3–5]. The driven electron dynamics can be monitored through optical signals [6–10] and through emitted electrons [11–16]. Very recently, Schiffrin et al [17] have demonstrated directed electron currents generated inside transparent dielectrics by carefully tailored laser pulses. In turn, the ultrafast response can characterise the impinging laser field [18]. Currently, avenues are explored to exploit such ultrafast modulation of electric currents for petahertz-scale signal processing [19] enabled by the short intrinsic time scale of the electron motion (∼1 fs), orders of magnitude faster than semiconductor electronics.
In this work, we explore a novel channel for the ultrafast electronic response that is unique to dielectrics with a non-centrosymmetric crystallographic structure: the generation of direct currents (dc) induced by strong optical laser pulses. Fully three-dimensional ab-initio simulations based on time-dependent density functional theory (TD-DFT) predict the generation of strongly nonlinear currents in α-quartz that are, in contrast to previously observed currents [17, 18], independent of the details of the laser pulse shape. The direction of the currents is found to be sensitive to the instantaneous laser intensity. Analysis of the spatio-temporal charge dynamics on the atomic length and time scales allows us to link this to the transition from nonlinear polarisation currents to directional tunnelling excitation, the latter being highly sensitive to the alignment between the laser polarisation and the chemical bonds in the crystal. We demonstrate that this transition may be investigated in a pump-probe setup leaving its marks as a change of the direction of the current as a function of the pump-probe delay.
2. Method
Theoretical exploration of ultrafast processes in solids faces the challenge to tackle the time-dependent many-body problem. TD-DFT has emerged as a versatile tool allowing for an ab-initio description of a variety of nonlinear and strong field processes [20–24] including in the solid state [9, 14, 16, 25–27]. Here, we employ a real-space, real-time formulation of TD-DFT [28–33] for the electronic dynamics induced by strong few-cycle laser pulses in α-SiO2 (α-quartz). Briefly, we solve the time-dependent Kohn–Sham equations (atomic units are used unless stated otherwise)
where i runs over the occupied Kohn–Sham orbitals The Hamiltonian
describes the system under the influence of a homogeneous time-dependent electric field of amplitude F0 along with vector potential in the velocity gauge. The periodic lattice potential is given by norm-conserving pseudopotentials of the Troullier–Martins form [34] representing the ionic cores (O(1s2) and Si(1s22s22p6)). The valence electron density is For the exchange and correlation potential we employ the adiabatic Tran–Blaha modified Becke–Johnson (TB-mBJ) meta-GGA functional [35–37]. The validity of the adiabatic approximation is still an open question [38]. For example, representing driven few-state dynamics, in particular Rabi oscillations, by adiabatic functionals has been found to be difficult [39, 40]. However, implementation of functionals that include effects non-local in time [41] into realistic three-dimensional simulations have remained a challenge. Currently, the applicability of adiabatic functionals to ultrafast phenomena can only be assessed by comparison with experiment. Good agreement for a wide array of observables including high-harmonic generation [20], electronic dynamics in metal clusters [22] and semiconductors [9], as well as reflectance and ablation [27] point to its usefulness and approximate validity. One likely reason for its success is that for these cases of many-particle systems the coherent dynamics involves averaging over a large number of states. The TB-mBJ functional accurately reproduces the band gap for SiO2 and yields good agreement with the experimental dielectric function over the spectral range of interest including at optical frequencies [42].
The time-dependent Kohn–Sham equations (1) are solved in the transverse geometry [43] to treat the bulk polarisation response of the infinitely extended system along the polarisation direction. We use a Cartesian grid with discretisation ∼0.25 a.u. in laser polarisation direction and ∼0.45 a.u. perpendicular to the polarisation direction in a cuboid cell of dimensions 9.28 × 16.05 × 10.21 a.u.3 employing a nine-point stencil for the kinetic energy operator and a Bloch-momentum grid of 43-points. The time evolution is performed with a 4th-order Taylor approximation to the Hamiltonian with a time step of 0.02 a.u. including a predictor-corrector step. The solution of equations (1) and (2) allows to analyze the time and space dependent microscopic vectorial current density
as well as the mean current density J(t) along the laser polarisation direction averaged over the unit cell of volume Ω,
The polarisation density [44] gives the charge density D(t) transferred by the pulse. The experimentally observed total charge Q will depend also on the details of the geometry of the laser focus and of the collection volume not explicitly treated in the following.
While dephasing due to elastic scattering is self-consistently included in TD-DFT, inelastic electron–phonon and electron–electron scattering is not. For low excitation densities and large band gap insulators electron-electron scattering can be neglected [45]. Electron–phonon scattering can be included on a phenomenological (i.e., non-self consistent) level via damping or relaxation rates. The possible influence of relaxation on the anisotropy of the induced nonlinear current are estimated by applying the quantum friction approach of Neuhauser and Lopata [46] and a purely phenomenological approach employing damping constants.
3. Nonlinear direct currents
First studies of the short-pulse induced current and charge transfer in polycrystalline SiO2 [17, 18, 33] found a sinusoidal dependence on the CEP, of the few-cycle electromagnetic field with -envelope and the carrier frequency of the IR laser. Sub-cycle control and steering of electrons required exquisite control over the instantaneous electric field The direction of the induced current was found to be determined by the CE phase and, thus, by the anisotropy of the few-cycle laser pulse. Here, we explore an alternate route towards steering, controlling and switching ultrafast currents that does not rely on control of the instantaneous field but on the instantaneous intensity dependence of a (dc) that is applicable for longer laser pulses where CEP effects become negligible. We propose to exploit the non-inversion symmetric lattice structure of the target (in this case polarization parallel to the direction of SiO2) to cause an intensity dependence of the direction of the charge transfer which can be tested in a simple pump-probe setting described in section 4.
Starting point is the observation that the total charge density transferred at the conclusion of the pulse can be split into a CEP dependent part with amplitude which vanishes upon averaging over and a residual or direct part tends to decrease with increasing pulse length while the magnitude of increases with the pulse length (figure 1). For pulse lengths exceeding a few optical cycles ( fs, nm), is approximately proportional to the pulse duration and dominates the signal exceeding DCEP by about one order of magnitude. The directionality of D0 is the focus of the present investigation.
Figure 1. Pulse length dependence of transferred charge at intensities (a) 1013 W cm−2, (b) W cm−2 and (c) 1014 W cm−2, each with pulse shape. Carrier-envelope phase dependent part DCEP (blue open squares); CEP independent part D0 (green full squares); dashed lines: linear slope going through the origin.
Download figure:
Standard image High-resolution imageThe charge transferred by the induced dc current, features a strongly nonlinear scaling with intensity or, equivalently, field strength (figure 2). The origin of this highly nonlinear response lies in the broken centrosymmetry of the SiO2 crystal along the direction. In general, generation of a directed flow of charge by a laser field requires a broken inversion symmetry. For few-cycle laser pulses with well-defined CEP, inversion symmetry is violated by a suitable choice of In the present case, it is not the temporal shape of the laser electric field but the electronic and crystallographic structure of the target the laser interacts with that causes ultrafast currents. This novel mechanism does not rely on delicate CEP control yet offers sub-cycle response and switching. The effect is also robust against damping and relaxation included into the simulations by adding a Drude-friction term
to the vector potential in equation (2) with damping constant α.
Figure 2. Carrier-envelope phase independent transferred charge density D0 (green full squares) as function of laser intensity ( pulse with full duration fs, eV) and absolute value (open squares). Red circles show the values of D0 for selected laser intensities including strong damping ( corresponding to a damping time of ∼1 fs). Power law (dashed line). Electrons move along for W cm−2 (upper panel) and along for .
Download figure:
Standard image High-resolution imageThe appearance of a direct current in a homogeneous medium under illumination, independent of the CEP, and linearly increasing with pulse duration, can be viewed as a nonlinear analogue to the well-known photogalvanic (PG) effect [47–51] as first qualitatively discussed by Alon [52]. Conventionally, the lowest order photogalvanic effect is described by
JPG is quadratic in the electric field components and linear in the time-averaged laser intensity For linearly polarized light, the photogalvanic tensor associated with the two-wave mixing in the second-order susceptibility is non-zero only in non-centrosymmetric crystals [53]. Microscopically, a variety of mechanisms may contribute to the PG effect such as asymmetric excitation, scattering, or recombination of electrons and electronic defects [48]. One important realisation is the so-called 'shift current' [53–55] due to the shift between the centre of charge of the valence electrons and the excited electrons in the conduction band. This shift current has been predicted to be important in several semi-conductors [53, 56, 57] and has been first experimentally verified for ferroelectrics [55].
The present nonlinear generalisation of the photogalvanic effects is obviously a signature of strong-field interaction with matter. This is underscored by the surprising observation of current reversal as a function of laser intensity (figure 2). We find a critical value of current reversal at W cm−2. Electrons move preferentially along the direction for lower intensities while they propagate along direction for higher intensities We have checked that this effect does not sensitively depend on the XC functional employed and is also obtained with the local density approximation [58]. We have also checked that the sign reversal persists when relaxation is included on a phenomenological level. We find the magnitude of the charge transfer to be somewhat reduced, its intensity dependence and, most notably, the sign inversion above a critical intensity Ic (figure 2) remain unchanged. We expect the sign reversal of the transferred charge to be experimentally observable.
We elucidate the microscopic mechanism for this reversal by analysis of the spatio-temporal charge dynamics. At lower intensities the multi-photon driven nonlinear polarisation current leads to a localized accumulation of charge in between the Si-O bonds as displayed in the time-averaged density fluctuations at the conclusion of the laser pulse (inset figure 3(e)). This implies the formation of an induced atomic-scale dipole around the oxygen atoms, i.e. a displacement of the centre of charge by vertical excitation, resembling the shift current mechanism of the standard photogalvanic effect but generalized to higher order reflected in the nonlinear intensity scaling of We have verified that for the direction of charge transfer agrees with that of the linear photogalvanic effect observed at lower intensities but higher photon energies. The nonlinear charge transfer along the direction can thus be viewed as the strong-field (or multi-photon) realisation of the shift current. For higher laser intensities the dominant charge transfer mechanism is excitation of the tilted conduction band by tunnelling. Tunnelling significantly depends on the local potential landscape in the tunnelling direction. We find tunnelling is enhanced when the bond direction is aligned with the laser field as illustrated by a strongly asymmetric current density at times near the maxima of the electric field (figures 3(c) and (d)). Excitations along the Si–O–Si bond chain give rise to a current after the conclusion of the laser pulse (figure 3(f)). Tunnelling excitation is more efficient along where the O–Si bond is more closely aligned with the laser field (bond-alignment angle ) while in direction tunnelling is suppressed because of the larger angle () between the bond axis and the laser polarisation (along ). Following tunnelling excitation to the conduction band, the current is mainly driven along the helical channel formed by the O–Si–O–Si–O chain along the direction (see labelled atoms in figure 3(c) and light contours in (d) and (e)). The reversal of the charge transfer and current direction is therefore most likely associated with the increased weight of tunnelling excitation, consistent with the onset of a sub-cycle time structure of charge transfer for The transition to tunnelling excitation is therefore accompanied by a reversal of the charge transfer and current direction. As tunnelling rates scale exponentially with the peak intensity the transition is quite abrupt suggesting its potential for femtosecond current switching. It is worth noting that despite the sudden switch in direction the overall intensity dependence of the transferred charge density is comparatively smooth. This observation is in line with other strong-field phenomena such as ionization that also display a smooth dependence on peak-field intensity of ultrashort pulses when crossing from the multi-photon to the tunnelling regime (e.g., [59]).
Figure 3. (a) Time-dependent polarisation density P(t) (solid green), laser field F(t) (black dashed, scaled) for fs, intensity W cm−2. (b) Time-dependent nonlinear current (equation (7), red solid). (c), (d) Snapshot of the current density in an - cut plane going through the central oxygen at a time (c), (d). Angles between the electric field and the O-Si bond are , . (e), (f) Current density after conclusion of the laser pulse for W cm−2 (e) and W cm−2 (f). While a current flows along the Si–O–Si bonds for W cm−2, we find no current for W cm−2 on the same colour scale. (inset panel e) Time-averaged density modulation after the laser pulse in plane cut through the central oxygen atom on a magnified colour scale.
Download figure:
Standard image High-resolution image4. Pump-probe protocol for dc currents
For all laser intensities, the dominant part of the CEP-independent dc charge transfer happens during the laser pulse (figure 3(a)), in contrast to the CEP-controlled ac charge transfer [33]. For high intensities the charge transfer shows sub-cycle time structure. The time-dependence of the tunnelling current can be conveniently analysed by the nonlinear response contribution after subtracting the linear-response current scaled to the instantaneous field,
with the conductivity determined for low intensity [43]. During the rise time of the pulse (figure 3(b)) is still as linear response prevails. However, once a field strength sufficient for tunnelling between neighbouring atoms is reached around fs, the nonlinear current shows strong spikes. While the linear response current is, to a good approximation, 90° out of phase with the electric field and the polarisation is in phase with the current spikes are in phase with F(t) as expected for tunnelling excitation. At later times (from fs on), remains in phase with, and becomes proportional to the laser field, indicative of a conductor-like linear response with a Drude (free carrier-like) conductivity for the tunnelling-induced electron population in the conduction band.
The present analysis of the nonlinear photogalvanic dc current suggests that the key control parameter is the instantaneous intensity I(t) rather than the cycle averaged intensity in the conventional photogalvanic effect or the instantaneous value of the field F(t) in the CEP controlled ac current. This sensitivity to I(t) can be explored in a pump–probe setting, in which the instantaneous intensity can by manipulated by the delay between pump and probe pulses. The pump–probe delay may therefore serve as knob for fast charge transfer by the nonlinear photogalvanic effect. To demonstrate this control we choose the intensity of both pump and probe pulses to be separately subcritical () with pump intensity W cm−2 and probe intensity W cm−2. However, the superimposed fields give rise to a maximum intensity of W cm−2 above Ic. The sign and amplitude of the induced current is controlled by the time delay between the laser pulses (figure 4).
Figure 4. (a) Laser-induced transferred charge density as a function of the peak-peak delay ( W cm−2, W cm−2, fs). Red squares, dashed line: TD-DFT simulations; solid line: model equation (8). Inset: magnified data around to 2 fs (dashed–dotted arrow). (b) Maximum instantaneous laser intensity for a given pump-probe delay. Insets: pulse shape after superposition of pump and probe pulse F(t) at fs (weak pulse before strong pulse) and at +10 fs.
Download figure:
Standard image High-resolution imageFor large positive and negative delays, the transferred charge saturates at the same positive value. In contrast, for near-zero delay where the maximum intensity is attained, the dc current switches direction and the transferred charge becomes negative. Remarkably, during the period of strong overlap the modulation of the dc current occurs on the sub-fs time scale resulting from the strongly varying maximum instantaneous laser intensity as a function of pump–probe delay (figure 4(b)).
Assuming that the charge transfer is governed by the central peak of the combined laser pulse, a simple estimate in analogy to equation (6) predicts
where denotes the sign function and is the maximum instantaneous laser intensity for pump–probe delay (figure 4(b)). In equation (8), we denote the nonlinear generalisation of the photogalvanic tensor by This simple model reproduces the temporal variation of in the full TD-DFT calculations remarkably well, underlining that the maximum instantaneous laser field drives the nonlinear photogalvanic effect through tunnelling near the field maximum.
5. Conclusions
We predict a nonlinear extension of the photogalvanic effect into the strong-field regime giving rise to ultrafast dc currents in insulators illuminated by multi-femtosecond laser pulses. We observe a strongly nonlinear intensity dependence and even a reversal of the induced currents above a critical intensity Ic associated with the transition from nonlinear polarisation currents to tunnelling excitation. The charge transfer is rather insensitive to details of the laser pulse shape and CEP but strongly dependent on the maximum instantaneous field strength. The latter may be controlled by the pump–probe delay in a two-pulse setup giving rise to a distinct sign change in the transferred charge as function of the pump–probe delay. The nonlinear photogalvanic effect opens up opportunities for light-field controlled femtosecond charge separation with relatively modest requirements on the driving laser. Even many-cycle pulses without CEP stabilisation can be used as the lattice structure instead of the CEP is employed to break the inversion symmetry along the laser polarisation axis. The nonlinear photogalvanic effect is conceptually simpler than the CEP dependent charge transfer since no elaborate steering of the conduction band electrons is necessary. Therefore, the effect is robust against changes in the laser pulse parameters. We envision the nonlinear photogalvanic effect may be useful for ultrafast signal processing as the sign of the current may be controlled by the time delay between two laser pulses only, and does not rely on stable and custom-tailored wave shapes nor short pulses with only a few cycles. This may be advantageous in particular for optical interconnects based on surface plasmon propagation [19] where the pulse shape and duration of a surface plasmon wave packet are difficult to control. The relatively sharp threshold intensity Ic for reversal of the current may provide a simple route towards femtosecond current switching and, moreover, a sensitive intensity calibration for laser pulses that directly measures the maximum electric field strength in the material. Finally, the photogalvanic effect may also be investigated by associated terahertz emission [60–62].
Acknowledgments
This work was supported by the Austrian Science Fund (FWF): P21141-N16, special research programmes SFB-041 ViCoM, SFB-049 Next Lite and doctoral college W1243. GW thanks the IMPRS-APS for financial support. X-MT was supported by a Grant-in-Aid for Scientific Research (C24540421) from the JSPS. KY acknowledges support by the Grants-in-Aid for Scientific Research Nos. 23340113 and 25104702. Calculations were performed using the Vienna Scientific Cluster (VSC) and the supercomputer at the Institute of Solid State Physics, University of Tokyo.