lopscience

¡opscience.iop.org

Home Search Collections Journals About Contact us My IOPscience

Quantum beats in the polarization response of a dielectric to intense few-cycle laser pulses

This content has been downloaded from IOPscience. Please scroll down to see the full text. View the table of contents for this issue, or go to the journal homepage for more

Download details:

IP Address: 137.195.150.201

This content was downloaded on 28/10/2014 at 06:00

Please note that terms and conditions apply.

New Journal of Physics

The open access journal for physics

Quantum beats in the polarization response of a dielectric to intense few-cycle laser pulses

M Korbman1, S Yu Kruchinin1 and V S Yakovlev1,2,3

1 Max-Planck-Institut für Quantenoptik, Hans-Kopfermann-Straße 1, D-85748 Garching, Germany

2 Department für Physik der Ludwig-Maximilians-Universitat München, Am Coulombwall 1, D-85748 Garching, Germany

E-mail: vladislav.yakovlev@lmu.de

New Journal of Physics 15 (2013) 013006 (16pp)

Received 5 October 2012 Published 4 January 2013 Online at http://www.njp.org/

doi:10.1088/1367-2630/15/1/013006

Abstract. We have investigated the polarization response of a dielectric to intense few-cycle laser pulses with a focus on interband tunnelling. Once charge carriers are created in an initially empty conduction band, they make a significant contribution to the polarization response. In particular, the coherent superposition of conduction- and valence-band states results in quantum beats. This quantum-beat part of the polarization response is affected by the excitation dynamics and attosecond-scale motion of charge carriers in an intense laser field. Our analysis shows that, with the onset of Bloch oscillations or tunnelling, the nonlinear polarization response becomes sensitive to the carrier-envelope phase of the laser pulse.

Contents

1. Introduction 2

2. Theory 3

3. Results and discussion 8

4. Conclusions and outlook 14 Acknowledgments 15 References 15

3 Author to whom any correspondence should be addressed.

Content from this work may be used under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike 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.

New Journal of Physics 15 (2013) 013006 1367-2630/13/013006+16$33.00

© IOP Publishing Ltd and Deutsche Physikalische Gesellschaft

1. Introduction

Multiphoton strong-field excitation of valence-band electrons has been intensively investigated over the last several decades. Until recently, the main motivation for experimental research has been laser damage of optical materials [1-3]. Detailed time-resolved measurements of electron dynamics used to be out of reach, as the characteristic times of the most essential processes belong to the attosecond domain. The situation has recently changed due to the progress in the generation of few-cycle laser pulses and the development of attosecond science [4, 5]. In this context, several recent papers should be pointed out. Gertsvolf et al [6] observed that an elliptically polarized 40 fs pulse of near-infrared (NIR) radiation changes its ellipticity as it propagates through a few micrometres of fused silica. Based on their numerical calculations, these authors concluded that the change of ellipticity is due to subcycle, attosecond-scale excitation dynamics of valence-band electrons. Their conclusions were confirmed by Mitrofanov et al [7], who used a non-collinear pump-probe technique to observe an attosecond-scale modulation of the electron density created by a 5 fs NIR pulse in a SiO2 sample. A different optical strong-field phenomenon was observed by Ghimire et al [8] in ZnO exposed to high-power few-cycle mid-infrared laser pulses: they observed the generation of high harmonics with properties that cannot be explained by the conventional nonlinear optics. The origin of these harmonics was linked to a combination of anharmonic electron motion and multiple Bragg reflections at the boundaries of the first Brillouin zone, also known as 'Bloch oscillations' [9]. Another manifestation of Bloch oscillations observed by the same team was a redshift in the optical absorption in ZnO [10]. Most recently, it was demonstrated that an intense few-cycle laser pulse not only excites charge carriers in SiO2, but it also drives measurable currents that can be steered by controlling the carrier-envelope phase (CEP) [11] of the laser pulse [12]. The theoretical interpretation of that experiment was given in terms of the Wannier-Stark localization [13, 14] that can be viewed as a frequency-domain description of Bloch oscillations. These developments created a new perspective for studying electron dynamics in solids exposed to strong fields, where the emphasis is placed on attosecond-scale dynamics.

The foundations for the theoretical description of extremely nonlinear interaction with strong fields, where conventional perturbation theory breaks down, were laid by Keldysh in his seminal papers [15, 16] where multiphoton and tunnelling excitations [17], which are two extreme regimes of interband transitions, were distinguished by the Keldysh parameter

&>la /mE g

Y = \ g. (1)

This parameter was introduced for a monochromatic laser field in the approximation of parabolic bands. In (1), &>L is the laser frequency, m = (1/me + 1/mh)-1 is the reduced mass of an electron and a hole, Eg = h&>g is the bandgap, e is the absolute value of the electron charge and F0 is the amplitude of the laser field in the medium. In the limit y ^ 1, where the field is weak and/or its frequency is large, an electron is excited from the valence band by absorbing an integer number of photons in a way that is well described by conventional perturbation theory. In the opposite extreme y ^ 1, where the field is strong and its frequency is sufficiently small, it is quantum tunnelling that promotes electrons to the conduction band. Obviously, this is the case if the external field is constant (&>L = 0), and interband tunnelling as a concept was first introduced by Zener [18] for time-independent fields. In the tunnelling regime, the excitation rate is often assumed to be determined by the instantaneous value of the electric field FL(t) [7].

Electron dynamics in the intermediate regime y ~ 1 are particularly complex—a laser field cannot be treated as a small perturbation in this regime and, at the same time, interband transitions cannot be described by pure tunnelling. A general expression for the cycle-averaged excitation rate was given by Keldysh in [16], but a completely satisfactory description of electron dynamics within an optical cycle of an ultrashort pulse is still missing. The dynamics can be modelled by solving the time-dependent Schrodinger equation (TDSE) in a suitable approximation, but there is no unambiguous way to distinguish between valence-and conduction-band electrons as long as the external field is present. This is analogous to the problem of determining the ionization rate of atoms and molecules in strong fields [19], for which many practical solutions exist, but there is no mathematically rigorous solution in the case of a few-cycle pulse. This motivated us to study the polarization response of a dielectric in the regime y ~ 1, as interband transitions directly affect polarization and, at the same time, polarization is a quantum-mechanical observable that determines all optical properties of a solid.

2. Theory

We are mainly interested in the strong-field regime where perturbation theory shows its shortcomings. In particular, the polarization response of a solid to a laser pulse is no longer described by a set of linear and nonlinear susceptibilities; therefore, quantum-mechanical simulations are necessary to model it. We solve the TDSE numerically in the velocity gauge, using the method described in [20]. For each value of the quasimomentum k, the electron wave function is represented in the basis of Bloch states |0k}, which are evaluated by solving the

single-electron stationary Schrodinger equation with an unperturbed Hamiltonian H0:

H o| 0k) = E'k\<Pl). (2)

Here, n is a band index. This equation is solved in the basis of Nmax plane waves:

0k (r) = cl j exp[i(k + K j)r], (3)

where K j denotes reciprocal lattice vectors, each of which can be expressed as a combination of the primitive vectors {b1, b2, b3}:

K = ^^ maba, ma e Z. (4)

Ansatz (3) ensures that all Bloch states satisfy 0k (r + R) = 0k (r) exp[ikR], where R is a vector in the direct Bravais lattice. This property immediately follows from exp[iKR] = 1.

The coefficients in expansion (3), as well as the energies Ek, are determined by diagonalizing the Hamiltonian H0 in the basis of plane waves exp[iKjr]. We assume that the interaction of a particular electron with lattice ions and other electrons is described by a pseudopotential U(r) and expand this pseudopotential in the plane-wave basis:

U (r) = Uj eiKjr, (5)

Uj = 1 f e-iKjr U (r ) dr, (6)

where Q is the volume of a unit cell. Writing the unperturbed Hamiltonian as

Ho = — t—V 2 + U (r ) (7)

and substituting ansatz (3) into (2), we obtain the following eigensystem for the expansion coefficients Ci i and energies E' :

( h2 2 \

^ {k + Kj) Sii + Uk'-K0 Clj = E"kCli■ (8)

Here, m0 is the electron rest mass, and the eigensystem can be solved independently for each quasimomentum k. The solutions of this eigensystem are normalized to satisfy

(0m 0) = Q j [€(r)]* 0 (r) d3r = Smn, (9)

mi _ I "I* ,43.

which implies that J^j ICn j I2 = 1. Here and in the following, a scalar product implies integration over one unit cell.

Having evaluated Bloch states |pf}, we use them as a basis to solve the TDSE and thus model the interaction of electrons with a laser pulse:

ih- №k(t)) = (H0 + Hint(t)Mk(t)), (10)

where Hn(t) is the interaction Hamiltonian. We perform our simulations in the velocity gauge, where an external electric field is described by the vector potential A(r, t). We also employ the dipole approximation: A(r, t) = A(t). In SI units, this corresponds to using the following form of the interaction Hamiltonian:

Hint (t) = — A(t) p. (11)

Note that in CGS units the right-hand side of (11) must be divided by the vacuum speed of light c. The substitution of the ansatz

|fk (t )) = £ ank (t )|(k} (12)

into (10) leads to the following system of coupled differential equations:

q /. \ Nmax

ih = El al (t) + - A(t) £ pf ak (t). (13)

d t m 0 ^

Nmax ^

pf = (Pl| P| Pl) = »£(l + Rj)(Clj) Clj (14)

are the matrix elements of the momentum operator p. An important advantage of the velocity gauge is that, as long as the dipole approximation holds and electron scattering is neglected, there is no coupling between different values of k. Mathematically, this is a consequence of

f (4(r))* p 4>lk(r) d3r = (4 Ip|4) £ ei(k-k'>Rj a <W

J R3 ,

where the summation is performed over all the vectors of the Bravais lattice R j, and we assume that both k and kf belong to the first Brillouin zone. Due to this property, equations (13) can be solved independently for each quasimomentum k. According to our experience, the downside of the velocity gauge is a larger number of bands required for convergence, as compared to length-gauge simulations, where different values of k are coupled to each other [21, 22], but even the two-band approximation may be sufficient [23].

At the beginning of a simulation, all electrons are supposed to reside in the valence bands (VB), while all the conduction bands are empty. We solve equations (13) for each of the electrons independently and then add single-electron contributions together to evaluate physical observables such as current density or polarization. Obviously, our approach is only applicable if (i) lattice vibrations have no significant effect on the investigated dynamics, (ii) effects related to electron correlation and electron-hole interaction are not important and (iii) distortions of the pseudopotential U (r) upon electronic excitations are negligible. While these approximations obviously have their limitations, single-electron models proved to be useful in studying basic physics related to phenomena occurring on the attosecond and few-femtosecond time scales [20, 24, 25].

The initial condition for solving equations (13) for a particular k is

ak (tmin) = V - (15)

Here, n is a band where the electron was before the interaction with a laser pulse.

In this paper, we are most interested in the macroscopic electric current density j (t) induced by the laser pulse

j (t) = £/ jkn (t) d3k, (16)

n eVB ^

where the integral is taken over the first Brillouin zone and jk n (t) represents single-electron contributions to the current density averaged over a unit cell:

jkn(t) = - —1 f d3r {Re [rkn(r, t)pfk,n(r, t)] + eA(t)\fk,n(r, t)|2}. (17)

mo ^ Jq

Here and in the following, we add index n to quantities related to single-electron wave functions in order to indicate the initial band. Thus, jk n (t) represents the contribution from all the bands to j (t) in a calculation where the electron initially occupied band n. The term containing the vector potential is required in the velocity gauge.

With the aid of (12), the single-electron current density can be expressed via probability amplitudes a\ (t):

jk.n (t) =--I Re

EKn (t)) *«in (t) pk

+ eA(t) I . (18)

Three different physical effects determine the induced currents: polarization due to the displacement of valence-band electrons, the light-driven motion of conduction-band electrons and quantum beats due to coherent superpositions of valence- and conduction-band states. In the presence of a strong external electric field, these contributions cannot be fully separated from each other, but the overall induced current density j (t) is defined unambiguously. It is this induced current that determines the change of the refractive index and the generation of new frequency components. However, it is conventional to describe optical response in terms of polarization, which we define as

P (t) = i t j (t') dt'. (19)

Note that the polarization defined this way is due to both bound and free electrons. Ordinarily, P (t) would be assigned to bound electrons, while j (t) would describe the motion of free electrons. We make no attempt to make this distinction because, to the best of our knowledge, there is no rigorous way to distinguish between bound and free electrons as long as the external field is present [19]. It is only after the laser pulse that projecting a wave function onto conduction-band Bloch states yields physically relevant excitation probabilities.

For simplicity, we perform our numerical simulations in one spatial dimension x, which we assume to be the polarization direction of the laser field. Most of the equations in this section can be transformed to their one-dimensional (1D) forms by replacing all vector quantities and 3D integrals with scalar quantities and 1D integrals: r ^ x, l ^ k. The only exception is (16), where this procedure would result in j (t) being measured in wrong units. Reducing the dimensionality of our original problem, we essentially assume that physical observables do not depend on the coordinates y and z, so that the integral over the first Brillouin zone in (16) reduces to

j (t) = n £ / jk , n (t) dk , (20)

neVB ^

where the factor n, measured in m-2, accounts for the integration over ky and kz. In practice, the value of n should be chosen to approximate some known optical properties of the solid, such as its refractive index or absorption.

From this point on, we will write our equations in atomic units (au), unless stated otherwise. In atomic units, h = e = m0 = 1, a unit of energy is equal to 27.21 eV, and a unit of electric field is equal to 5.142 x 1011 Vm-1. Because h = 1, we will interchangeably use the words 'energy' and 'frequency'.

We adjusted our model potential U(x) to obtain a bandgap close to that of SiO2: 8.95 eV [26]. To this end, we set the lattice constant to a0 = 0.5 nm = 9.45 au, which is one of the lattice constants in a-quartz, and use the following expression for the central unit cell:

U(x) = -0.7 (1 +tanh[x + 0.8]) (1 +tanh[-x + 0.8]), |x| ^ a0/2. (21)

Outside the central unit cell, the potential is continued periodically: U (x + a0) = U (x). The band structure that corresponds to this potential is shown in figure 1(a). At k = 0 the energies E% of the first four bands are -35.48, -8.33, 0.62 and 13.13 eV. We regard the first two bands with negative values of E% as valence bands (the lowest valence band is not shown in figure 1(a)), while all the bands with E% > 0 are regarded as conduction bands.

Figure 1. (a) The uppermost valence band and the two lowest conduction bands of the model potential (21). (b) A comparison of the refractive index evaluated according to (24) with the refractive index of fused silica [27].

We used the following expression for the vector potential of the (linearly polarized) laser field acting on electrons in the solid:

A(t) = -0(tl -It |) — cos4 ( JM sin(^Lt + ^cep). (22)

wl \2tlj

Here, F0 is the amplitude of the electric field, wL is the central frequency of the laser pulse, <pCEP is its CEP, 0(x) is the Heaviside step function and tl is related to the full-width at halfmaximum (FWHM) of |A(t)|2 by FWHM = 4arccos(2-1/8)tl/n ^ 0.5224tl. The FWHM of | A(t) |2 is very close, although not precisely equal, to the FWHM of | F(t) |2. The external electric field acting on electrons is, by definition,

F (t) = - A (t). (23)

Note that, in this work, we do not make any attempt to evaluate the screening field created by the displacement of electrons—this field is assumed to be a part of A(t) and F(t).

For a weak pulse, the polarization response is linear, so that, in the frequency domain, P (w) = x(w) F (w), and the refractive index can be evaluated from the linear susceptibility

n (w) = 71 + x(w). (24)

In figure 1(b), we compare the refractive index of fused silica with n(w) evaluated using our model. The refractive index is plotted against the laser wavelength X = 2nc/w. For these simulations, we used several pulses with F0 = 10-5 au = 5 x 106 V m-1, FWHM = 4n/wL, and values of wL that allowed us to cover the spectral range presented in figure 1(b). Normalizing the polarization response, we set n = 0.111 au in (20) and use this value henceforth. Given the simplicity of our model, we find that the agreement with the measured refractive index is very satisfactory.

The next section reports on simulations with much more intense fields, where multiphoton excitations populate conduction bands. For those simulations, we had to use 15 bands in order to achieve numerical convergence in solving the TDSE. A smaller number of bands results in discrepancies that are most visible in the spectral range occupied by low-order harmonics, while

? 2 "c

1 0 ? -1

3 -2 -3

12 11.5 11 10.5 10 9.5 9 8.5

-2 0 2 4 time (optical cycles)

0 (atom -1 ?

-2 210 -3

1 2 3 4 5 6 7 time (optical cycles)

-2 0 2 4 time (optical cycles)

Figure 2. (a) The polarization P(t) induced by a laser pulse with aL = ag/4 (XL = 557 nm), a peak intensity of 3 x 1013 Wcm-2, and ^CEP = 0. The electric field F(t) is shown by the dotted blue line; the FWHM of the pulse corresponds to the time range from -1.5 to 1.5 optical cycles. (b) S(t, a), which is the result of applying the wavelet transform to P (t). The dashed white line marks the energy Eg = hag. (c) The solid red line is the sum of the bandgap energy Eg and the ponderomotive energy Up(t). For each time t, the dashed blue line depicts the above-bandgap frequency where S (t , a) has a maximum. (d) The quantum-beat part of the polarization response F-1[wW(a) P (a)] (grey) and its envelope (red). The spectral window used for the analysis is given by (28).

increasing the number of bands to 20 has a negligible effect on the polarization response even for the most intense pulses that we used in our modelling.

3. Results and discussion

For all the simulations in this section, we used pulses with the FWHM equal to three optical oscillations: FWHM = 3 x 2n/aL. The electric field F(t) of such a pulse is shown in figure 2(a) by the dotted blue line plotted against the number of optical cycles aLt/(2n). The solid red line in this figure shows the polarization response P (t) evaluated according to (19) for the case where the central frequency is equal to a quarter of the bandgap (aL = ag/4, XL = 2nc/aL = 557 nm), and the peak intensity is equal to IL = 3 x 1013 Wcm-2, which corresponds to an amplitude of the electric field of F0 = 0.029 au = 1.5 V A-1. The scales chosen for the figure emphasize that, at the leading edge of the laser pulse, P (t) is proportional to F (t). At later times, the polarization response becomes increasingly nonlinear and high-frequency oscillations appear, which persist even after the laser pulse ends. The frequency of these oscillations is close to the bandgap

frequency wg = Eg/h. Thus, these are quantum beats appearing due to the presence of coherent superpositions of valence- and conduction-band states. This part of the polarization response will be the main topic of our further discussion. In order to distinguish it from conventional harmonics due to x(3), X(5), X(7\ ... nonlinearities, we chose the central laser frequency to be an even divisor of the bandgap, taking advantage of the fact that harmonics of even orders are absent if the potential U(x) possesses the inversion symmetry U(-x) = U(x), like our model potential (21). Note that, due to the same symmetry, the four-photon transition from the uppermost valence band to the lowest conduction band at k = 0 is forbidden, while interband tunnelling can efficiently populate conduction-band states in the middle of the Brillouin zone. Indirectly, this is confirmed by the fact that the fast oscillations in figure 2(a) only appear at very

high intensities. At a peak intensity of IL = 10 W cm-2, which is just three times smaller than the one used for figure 2, the amplitude of quantum-beat oscillations is smaller by approximately a factor of 400, which corresponds to a drop of their intensity by five orders of magnitude.

We investigate the quantum-beat part of the polarization response by means of the wavelet analysis using the Morlet wavelet [28]:

5 (t, œ)

= œ |£

P(tf)W(œ(tf- t)) dt7

W (x ) =

exp ix

The parameter tw determines the size of the temporal window, and it must be chosen to provide an optimal compromise between the temporal and spectral resolution. Our choice was tw = 3wg/Wl.

The outcome of this analysis, S(t, w), is shown as a false-colour diagram in figure 2(b), the bandgap energy being indicated by a dashed white line. For these parameters, the quantum-beat signal completely dominates the third-order harmonic, which is barely visible in the plot. In this time-frequency analysis, the third-harmonic signal is centred at 6.7 eV, and it is temporally confined to the FWHM of the light pulse: wL|t I/(2n) < 1.5.

One of the most striking features observed in figure 2(b) is the fact that the quantumbeat signal initially appears at frequencies exceeding wg by roughly 2eV/h. The instantaneous frequency of quantum beats then gradually decreases, approaching wg at the tailing edge of the laser pulse. In other words, the quantum-beat signal is negatively chirped—its lower-frequency components are delayed with respect to the higher-frequency ones. This cannot be explained by the presence of the fifth-order harmonic—even though the fifth harmonic occupies the relevant spectral range (5hwL = 11.2 eV), it must be confined to an even smaller temporal range than the third harmonic, so that S(t, w) can be considered free from its contribution for wLt/(2n) > 1, where most of the quantum-beat signal is observed.

We seek an explanation for the chirp of the quantum-beat signal in the laser-driven motion of charge carriers. The average kinetic energy of an electron-hole pair is the ponderomotive energy. For a monochromatic field with an amplitude F0, the ponderomotive energy is given, in the approximation of parabolic bands, by

F 2 F0

4 m œ>L

u 10 8

i \ №<0 -

i i \ \

-2 0 time (optical

cycles)

Eg+Up(t) -

time (optical cycles)

Figure 3. The wavelet analysis of the polarization induced by three-cycle pulses with a peak intensity of 3 x 1013 Wcm-2, ^CEP = 0, and the following central frequencies: (a) aL = ag/5 (XL = 697 nm); (b) aL = ag/6 (XL = 836 nm); (c) aL = ag/8 (Al = 1108 nm); (d) aL = ag/10 (kL = 1385 nm). The units in colour coding are arbitrary, but they are the same as those used for figure 2(b).

Here, m is the reduced mass of an electron and a hole. We obtain the effective masses of electrons and holes by fitting parabolas to the lowest conduction and uppermost valence bands in the region |k| < n/(2a0), which yields me = 0.34 au and mh = 3.61 au. These values are in good agreement with those that can be found in the literature [29, 30]. Thus, the reduced effective mass m in (27) is equal to m = (1/me + 1/mh)-1 = 0.35 au. The cycle-averaged total energy of an electron-hole pair is equal to Eg + Up (t), where Up (t) is to be evaluated by replacing F0 in (27) with the pulse envelope. This quantity is plotted as a red line in figure 2(c). It is fairly close to the energies where S(t, a) has local maxima with respect to a. We observed this kind of negative chirp in numerous simulations with different parameters of the laser pulse, some of which are shown in figure 3. Therefore, we conclude that the frequency of quantum beats exceeds ag by approximately the ponderomotive energy.

Figure 2(d) presents a different kind of analysis that highlights the quantum-beat part of the polarization response. In order to suppress the third- and fifth-order harmonics, we multiplied the Fourier transform of P (t) with a soft spectral window that cuts all frequency components that are further than aL/2 from the bandgap frequency:

(œ) = e(

^ cos2 (n

/o>l )

The grey line in figure 2(d) shows the inverse Fourier transform of w (a) P (a). In other words, it represents the part of the polarization response that oscillates at frequencies close to the bandgap frequency. The envelope of these fast oscillations is shown as a solid red line, and the dotted blue line depicts the square of the laser field. From this analysis, it is obvious that the quantum-beat part of the polarization response experiences its most rapid increase at a time that lies about two optical cycles after the peak of the laser pulse. This is counter-intuitive, as the quantum-beat signal is due to interband excitations, and the rate of the interband excitations is expected to reach its peak at the peak of the laser pulse. Another prominent feature apparent from figure 2(d) is the decay of quantum beats that begins at the trailing edge of the laser pulse. As our model does not account for any relaxation phenomena (such as electron-phonon interaction or spontaneous emission of radiation), the origin of this decay lies in the dephasing of dipole oscillators associated with each pair of coherently populated valence- and conduction-band states. Each such coherent superposition creates a current jk(t) described by (18). If the

q (t) = aqk (ti) exp[-i(t - ti)El

laser pulse ends at a time t1, then al(t) = al(t1) exp[-i(t —11)El], so that

jk (t) = —Re

\\* i(t—ti)(El — E[)

J>l (ti))* ak(ti)plle- —k ^

For brevity, we have omitted the index specifying the initial band of the electron. For a given pair of bands l = l, (29) describes a current oscillating with a k-dependent frequency El — Elk. Because of this dependence, the net current density (20) attenuates as most currents jk (t) get out of phase with each other.

Figure 3 presents the wavelet analysis of polarization response functions for several other values of the central laser frequency aL. The dashed white lines depict Eg + Up(t). In figures 3(a) and (b), these lines closely follow local maxima of S(t, a). For smaller laser frequencies, the ponderomotive potential Up significantly exceeds the width of the lowest conduction band (see figure i(a)). This is the regime where Bloch oscillations [3i] must play an important role [8].

In figure 3(a), where valence-band electrons can be excited by five-photon absorption, we made an exception from our choice to use those values of aL that are even divisors of ag. Still, the time-frequency analysis looks very similar to that for aL = ag/4, shown in figure 2(b), and the contribution from the fifth harmonic generation is so weak that it is practically invisible in our colour scheme. For all the cases shown in figure 3, we observe the same qualitative features as those in figure 2(b): the frequency components of the polarization response that lie above the bandgap appear delayed with respect to the peak of the laser pulse, and the instantaneous frequency of these oscillations decreases with time, eventually approaching ag. We also observe a rapid decrease in the magnitude of S(t, a).

To further understand the origin of these trends, we plot the sum of all conduction-band populations at the end of our simulations in figure 4. These are the probabilities of finding an electron with a given quasimomentum k in one of the conduction bands at a time t > tl where the laser field is absent (since we do not account for relaxation phenomena, populations do not change in the absence of external fields). The distribution of excited electrons dramatically changes as aL decreases from ag/4 to ag/i0: the pronounced peak at k = 0 disappears, and many irregular peaks appear across the Brillouin zone. At the same time, the distributions become increasingly sensitive to the CEP of the laser pulse. We consider two physical phenomena that may be responsible for this trend: interband tunnelling and Bloch oscillations.

a0 k / n a0 k / n

Figure 4. The probability for an electron that initially occupies a state in the uppermost valence band with a certain quasimomentum k to be excited to one of the conduction bands by the end of a three-cycle laser pulse with a peak intensity of 3 x 1013 W cm-2. Results obtained for cosine pulses (^CEP = 0, solid red lines) are compared with those for sine pulses (^CEP = n/2, dotted blue lines). The central laser frequency aL, the Keldysh parameter y and the ratio of the Bloch frequency aB to the laser frequency are specified in the plots.

The first phenomenon, interband tunnelling, represents an extreme regime of interband excitations where the Keldysh parameter y is much smaller than 1. In this regime, interband transitions should be viewed not as a result of absorbing a certain number of photons, but as a result of quantum tunnelling [18]. Also, in the tunnelling regime, the CEP is an important parameter because it is the field of the laser pulse, rather than its envelope, that controls the tunnelling rate. According to figure 4 and all the other similar simulations we have performed, the onset of the CEP dependence coincides with y decreasing below a value approximately equal to 0.7. Furthermore, the tunnelling picture provides an intuitive explanation for the multitude of peaks in figures 4(c) and (d): they may be attributed to the interference among electron wave packets launched to conduction bands within different half-cycles of the laser pulse, which is analogous to the above-threshold ionization spectra of atoms exposed to intense laser pulses [32].

What the tunnelling picture alone cannot explain is the disappearance of the peak at k = 0. Indeed, all analytical models for the tunnelling rate predict that it should rapidly decrease as the bandgap increases [16, 18, 33], and our model potential has the smallest bandgap at k = 0.

Therefore, the motion of excited electrons in conduction bands must play a significant role [23], Bloch oscillations being the most prominent manifestation of this motion. An electron excited to the bottom of a band at a zero-crossing of the vector potential (A(t0) = 0) may reach the boundary of the first Brillouin zone if, at some later moment ti, the condition | A(ti)| = n/a0 is satisfied. Thus, Bloch oscillations should be considered important if the amplitude of the vector potential F0/(L is comparable to n/a0. Introducing the Bloch frequency

ea0 F0

(b = —— (SI units), (30)

this condition can be expressed as (B/(L = n. The labels in figures 4(a)-(d) provide values of (B/(L in units of n for all the four simulations. The appearance of the CEP dependence and the complex structures in the distribution of conduction-band electrons coincide not only with Y becoming smaller than 1, but also with (B /(L becoming larger than n. In fact, this is not a coincidence, as

(B a0x ImE g

— = -g (SI units), (31)

(l hi Y

as follows from equations (1) and (30). For our parameters, (B/(L ^ 1.02n/Y, so that the onset of tunnelling coincides with the onset of Bloch oscillations, irrespective of the laser pulse parameters.

Without interband excitations, Bloch oscillations obviously cannot occur, but is it the nature of interband transitions or the Bloch oscillations themselves that are mainly responsible for the observed CEP dependences? While we cannot answer this question with certainty, we can argue that interband tunnelling alone is not sufficient for a polarization response to be CEP dependent. Indeed, tunnelling must be the dominant excitation mechanism in the parameter regime of figure 4(a) because, in spite of the forbidden four-photon transitions at k = 0, the bottom of the conduction band is strongly populated after the laser pulse. However, the polarization response is almost CEP independent. Therefore, our simulations with (L = (g/4 provide an example where the dominance of interband tunnelling over (perturbative) multiphoton excitations does not lead to a significant CEP dependence.

The quantum-mechanical observables that we have discussed so far cannot be directly measured in a realistic experiment. The most accessible quantity for experiments is the spectrum of emitted radiation. In order to evaluate such spectra, one needs to account for absorption and phase-matching [9], which is beyond the scope of this work. Nevertheless, important conclusions about these spectra can be made by analysing the Fourier transform of the polarization response P (t) that we evaluate with our model. In figure 5, we plot polarization spectra for the same simulations as those presented in figure 4. We see that the spectra are also sensitive to the CEP of the laser pulse. Similarly to excitation probabilities, the CEP dependence is very weak for (L = (g/4, but it becomes very pronounced already at (L = (g/6, and it affects not only the quantum-beat part of the polarization response, but also the conventional harmonics of third, fifth, etc orders.

Consistent with the results shown in figure 3, the quantum-beat signal rapidly decreases with decreasing (L. At the same time, according to figure 4, the excitation probabilities are comparable in all the four cases. Consequently, the decrease in the polarization response is due to dephasing. Most importantly, we used longer pulses for longer wavelengths in order to have the same number of optical cycles in all the pulses. As a result, dephasing played a more

14 (a)

100 10"1 10"2 10"3 10-4 10"5 10"6

1 1 1 PCEP " 0 A 1 1 fflL = œg/4 = 2.2 eV (D)

" PCEP " n/2 Y = 0.93 -

V. fflg/fflL = 1.1 n 'cn

\ _ 'E

/ N. .Q

_ y-S \ ' \ ÄL

Jr \ '/' ty

if \ /1 \ 1 W

" :y \/Oil \ l ' ten

!:i i i , .\ir c

6 8 10 photon energy (eV)

10"° -

10"° -

6 8 10 photon energy (eV)

10° 10" 10" 10" 10" 10" 10"

9CEP = 0 9cep = n/2

l-;-1-1-

ffl|_ = œg/6 = 1.5 eV

Y = 0.62 fflß/ffl|_ = 1.6 n

6 8 10 photon energy (eV)

(a10 ty

ten 10

6 8 10 photon energy (eV)

Figure 5. The power spectra of the polarization response for ^CEP = 0 (solid red lines) and <pCEP = n/2 (dotted blue lines). The parameters of the laser pulse are the same as those used for figures 4(a)-(d); the values of the laser frequency &>L are specified in plot labels. The double-dashed black vertical lines mark the bandgap energy. Note that a part of the quantum-beat spectrum appears below the bandgap.

significant role for the 14 fs pulse with = &>g/10 than for the 5.5 fs pulse with &>L = /4. In addition to that, the dephasing time is inversely proportional to the spectral width of the excitation, and excited electrons occupy a larger part of the Brillouin zone in the tunnelling regime (see figure 4).

4. Conclusions and outlook

We have analysed the polarization response of a model dielectric resembling SiO2 to few-cycle laser pulses that are strong enough to efficiently excite electrons from valence to conduction bands. These interband transitions create coherent superpositions of states that manifest themselves as quantum beats in the polarization. Obviously, the dynamics of interband excitations and dephasing determine the properties of the quantum-beat part of the polarization response. It is less obvious that the instantaneous frequency of quantum-beat oscillations changes with time, approaching the bandgap frequency as the laser field attenuates. We find that, during the laser pulse, the quantum-beat frequency (averaged over a laser cycle) is

larger than the bandgap frequency by Up/h, unless the ponderomotive energy Up exceeds the width of the lowest conduction band. Thus, the emitted burst of high-frequency radiation is negatively chirped. We have also observed that quantum beats appear in the polarization response significantly delayed with respect to the peak of a laser pulse, even though most electron-hole pairs must appear when the laser field is strongest. Apparently, the effect of a strong field on charge carriers suppresses single-photon transitions from conduction- to valence-band states. A satisfactory explanation of this phenomenon probably requires the development of an analytical theory, which may be triggered by this work.

For a dielectric such as SiO2, the onset of interband tunnelling (y ~ i) coincides with the onset of Bloch oscillations (aB ~ naL). In the tunnelling regime, we have observed that the nonlinear polarization and the final distribution of charge carriers are sensitive to the CEP of a laser pulse. In particular, the spectrum of the emitted high-frequency radiation depends on the CEP, which should be possible to observe in experiments. Such a measurement would have to overcome obstacles related to the strong absorption above the bandgap and the phase mismatch below the bandgap [9]. Also, it must be mentioned that the coherent few-femtosecond pulse radiated by quantum beats spectrally overlaps with the incoherent radiation due to fluorescence, which our model does not account for. Some ultrashort temporal gating may therefore be necessary to isolate the radiation due to quantum beats. Still, in spite of all these obstacles, such measurements should be feasible at the current stage of technology.

Acknowledgments

The authors gratefully acknowledge fruitful discussions with Professor F Krausz and Professor K Yabana. This work was supported by the DFG Cluster of Excellence: Munich-Centre for Advanced Photonics.

References

[1] Stuart B C, Feit M D, Herman S, Rubenchik A M, Shore B W and Perry M D i996 Phys. Rev. B 53 i749-6i

[2] Lenzner M, Kriiger J, Sartania S, Cheng Z, Spielmann C, Mourou G, Kautek W and Krausz F i998 Phys.

Rev. Lett. 80 4076-9

[3] Mao S S, Quere F, Guizard S, Mao X, Russo R E, Petite G and Martin P 2004 Appl. Phys. A 79 i695-709

[4] Hentschel M, Kienberger R, Spielmann C, Reider G and Milosevic N 200i Nature 414 509-i3

[5] Krausz F and Ivanov M 2009 Rev. Mod. Phys. 81 i63-234

[6] Gertsvolf M, Spanner M, Rayner D M and Corkum P B 20i0 J. Phys. B: At. Mol. Opt. Phys. 43 i3i002

[7] Mitrofanov A V, Verhoef A J, Serebryannikov E E, Lumeau J, Glebov L, Zheltikov A M and Baltuska A 20ii

Phys. Rev. Lett. 106 i4740i

[8] Ghimire S, DiChiara A D, Sistrunk E, Agostini P, DiMauro L F and Reis D A 20ii Nature Phys. 7 i38-4i

[9] Ghimire S, DiChiara A D, Sistrunk E, Ndabashimiye G, Szafruga U B, Mohammad A, Agostini P,

DiMauro L F and Reis D A 20i2 Phys. Rev. A 85 043836

[10] Ghimire S, DiChiara A D, Sistrunk E, Szafruga U B, Agostini P, DiMauro L F and Reis D A 20ii Phys. Rev.

Lett. 107 i67407

[11] Baltuska A, Udem T, Uiberacker M, Hentschel M and Goulielmakis E 2003 Nature 421 6ii-5

[12] Schiffrin A et al Optical-field-induced current in dielectrics Nature at press (doi:i0.i038/natureii567)

[13] Wannier G H i960 Phys. Rev. 117 432-9

[14] Durach M, Rusina A, Kling M F and Stockman M I 20ii Phys. Rev. Lett. 107 086602

[15] Keldysh L V i958 Sov. Phys.—JETP 7 788-90

[16] Keldysh L V i965 Sov. Phys.—JETP 20 i307-i4

[17] Wegener M 2004 Extreme Nonlinear Optics: An Introduction (Advanced Texts in Physics) (Berlin: Springer)

[18] Zener C i934 Proc. R. Soc. Lond. A 145 523-9

[19] Yakovlev V S, Korbman M and Scrinzi A 20i2 Dressed bound states for attosecond dynamics in strong laser

fields Chem. Phys. at press (doi:i0.i0i6/j.chemphys.20i2.0i.005)

[20] Bachau H, Belsky A N, Martin P, Vasil'ev A N and Yatsenko B N 2006 Phys. Rev. B 74 2352i5

[21] Blount E11962 Formalisms of Band Theory (Solid State Physics vol i3) ed F Seitz and D Turnbull (New York:

Academic) pp 305-73

[22] Glutsch S 2004 Phys. Rev. B 69 2353i7

[23] Golde D, Meier T and Koch S W 2008 Phys. Rev. B 77 075330

[24] Kazansky A K and Echenique P M 2009 Phys. Rev. Lett. 102 i7740i

[25] Brabec T 2008 Strong Field Laser Physics (Springer Series in Optical Sciences) (Berlin: Springer)

[26] Laughlin R B i980 Phys. Rev. B 22 302i-9

[27] Malitson I H i965 J. Opt. Soc. Am. 55 i205-8

[28] Chui C K i992 An Introduction to Wavelets (San Diego, CA: Academic)

[29] Schneider P M and Fowler W B i976 Phys. Rev. Lett. 36 425-8

[30] Chelikowsky J R and Schluter M i977 Phys. Rev. B 15 4020-9

[31] Bloch F i929 Z. Phys. A 52 555-600

[32] Delone N B and Krainov V P 2000 Multiphoton Processes in Atoms (Springer Series on Atoms + Plasmas)

(Berlin: Springer)

[33] Kane E O i960 J. Phys. Chem. Solids 12 i8i-8