Electron wavefunction in strong and ultra-strong laser field

1. Electron wavefunction in strong and ultra-strong laser field One- and two-dimensional ab initio simulations and models Jacek Matulewski Division of Collision Theory and Nonlinear Systems Faculty of Physics, Astronomy and Informatics Nicolaus Copernicus University in Toruń, Poland Marburg, 20 I 2004

2. Outline 1. Stabilisation and related phenomena in 1D strong field ionisation 2. Beyond the dipole approximation - magnetic drift in 2D ultra-strong field 3. Brief report on other projects (control of wavefunction by pulses) Authors: Andrzej Raczyński, Jarosław Zaremba and Jacek Matulewski

3. Outline 1. Stabilisation and related phenomena in 1D strong field ionisation 2. Beyond the dipole approximation - magnetic drift in 2D ultra-strong field 3. Brief report on other projects (control of wavefunction by pulses)

4. Language Space built with : Initial state = ground bound state of quantum system Time depended quantities we’re interested in

5. Simulation in 1D - typical result for strong laser field Evolution of wavefunction (fast oscillations with frequency of laser field was removed) Potential well: 2a = 0.244 a.u. = 1.3·10-11 m V0 = 2.049 a.u. = 55.7 eV Electric field: w = 1 a.u. = 6.6·1015 s-1 e = 1 a.u. = 5.1·1011 V/m (I = e 2 = 1016 W/cm2) Phase and pulse shape must be chosen carefully because of possibility of electron escape due to fast drift(Newton’s equations)

6. Simulation in 1D - typical result for strong laser field Population of initial state: fast oscillations Photoelectron spectra: ATI phenomena

7. Simulation in 1D - slow drift and stabilisation Evolution of wavefunction in e = 5 a.u. (precisely: of its part located near potential well) • Wavefunction properties • slowly changes position • remains its shape

8. Simulation in 1D - slow drift and stabilisation Some photos of wavefunction in e = 5 a.u. (precisely: of its part located near potential well) • Wavefunction properties • slowly changes position • remains its shape The idea is to find the potential which keeps shape of the wavefunction unchanged and model its movement due to perturbation of free oscillating electron by well potential

9. Transformation to Kramers frame , where and is putted together fast oscillations and slow drift Spreading of oscillating potential around (x0 is average value of b0(t)) Simulation in 1D - Generalised KH model Phys Rev A 61, 043402 (2000) One dimensional Schrödinger equation for one electron in laser field (dipole approximation):

10. Replacing time depended potential with HK well which is zeroth term of its Fourier expansion: In our parameters HK well has only one bound state Approximation: Function describing slow drift is well fitted by solution of Ehrenfest equation The reason of slow drift is interaction of wavefunction with the well edges Simulation in 1D - Generalised KH modelPhys Rev A 61, 043402 (2000)

11. Simulation in 1D - Generalised KH model (results) Phys Rev A 61, 043402 (2000) One period of oscillations - marks of the interaction with the potential well

12. Simulation in 1D - Generalised KH model (results) Phys Rev A 61, 043402 (2000) Position of the electron

13. Simulation in 1D - Generalised KH model (results) Phys Rev A 61, 043402 (2000) Occupation of the potential well ground state (initial state)

14. Transformed back to laboratory frame wavefunction of HK eigenstate The way of switching on the pulse Slow drift and fast oscillations Time of live in HK-well(Volkova et all, Zh. Eksp. Teor. Fiz. 106, 1360 (1994)) Stabilisation is not permanent because of HK state is ionised (higher terms of Fourier expansion) Simulation in 1D - Generalised KH model (results) Phys Rev A 61, 043402 (2000) What else influence the final population of the initial state?

15. Simulation in 1D - Summary Phys Rev A, 61, 043402 (2000) 1. Stabilisation in Kramers- Henneberger well 2. Slow drift of almost free electrons (interaction with well)

16. Outline 1. Stabilisation and related phenomena in 1D strong field ionisation 2. Beyond the dipole approximation - magnetic drift in 2D ultra-strong field 3. Brief report on other projects (control of wavefunction by pulses)

17. Typical simulation set in 3D: wave propagation direction: y electric field polarisation direction: x magnetic field polarisation direction: z

18. dipole approximation Omitting free evolution in z direction Simplifying the y dependence of hamiltonian Atom + pole lasera w 1D: stabilisation - HK well slow drift Atom + laser field in 2D: stabilisation - HK well in 2D magnetic drift Possible description approaches Atom + laser field in 3D without dipole approximation: classical description of free electron (Lorentz force) quantum description of atom in classical external field

19. Electromagnetic field: Lorentz force acting on free electron: Equations of motion: and trivial equation for z(t) Classical description in 3D (de facto 2D)

20. Equations of motion: Solution: x(t) and y(t) Classical description in 3D (de facto 2D)

21. Dipole approximation + unitary transformation Two-dimensional calculation pA (no dipole approximation!!!) One-dimensional calculations dE (oscillations, stabilisation, slow drift) Quantum description in 3D in Coulomb gauge Schrödinger equation (omitting evolution in z direction):

22. Reiss, Phys. Rev 63 013409 (2000) Magnetic drift in propagation direction Stabilisation (HK well) Solutions of two-dimensional Schrödinger equations in pA gauge, no dipole approximation, ultra-strong laser field atom: radial well a = 1 a.u., V0 = 2 a.u. field: w = 1 a.u., E0 = 15 a.u.

23. Wavefunction - phenomenology of stabilisation Phys Rev A 68, 013408 (2003) atom: radial well a = 1 a.u., V0 = 2 a.u. field: w = 1 a.u. 4T, 7T, 10T, E0 = 20 a.u. E0 = 15 a.u. two scattering centres in x = 0 and x = 30 a.u.

24. Wavefunction - phenomenology of stabilisation Phys Rev A 68, 013408 (2003) Characteristics of wavefunction 1. Regardless the phase of laser field magnetic drift is always directed to +y 2. Drift is linear in time + 2w oscillations 3. Constant drift velocity depends on E02 4. Wave of electron finding probability has double frequency in y direction E0 = 20 a.u.

25. spreading of vector potential in series Patel et al., Phys.Rev.A 64 013411 (2) Lowest order approximation of nondipol Schrödinger equation (additional term proportional to y in hamiltonian describing magnetic field) transformation to frame of oscillating electron transformation removing A2 term (3) Equation in Kramers frame with additional magnetic field averaging of oscillating potential (4) Electron motion equation in HK well and coupling with continuum by magnetic field Quantum model of stabilisation and magnetic drift in 2D Simplifying the hamiltonian (1) Schrödinger equation with no dipole approximation (pA gauge)

26. In high frequency regime the distribution of wavepacket in momentum space is concentrated around the zero Term describing magnetic field (2) Quantum model of stabilisation and magnetic drift in 2D Looking for hamiltonian with simplier dependence on y (1)

27. Transformation to frame moving with an electron Transformation removing the A2 term (3) In Kramers frame potential oscillates Quantum model of stabilisation and magnetic drift in 2D Transformation to the Kramers frame (2)

28. where and (4) „typical” time independed part of hamiltonian with bounding potential Term similar to dE (but here all in pA) Quantum model of stabilisation and magnetic drift in 2D Fourier expansion of the time depended potential of the oscillating well In high frequency regime one can replace fast oscillating potential with time averaged potential of Kramers-Henneberger well:

29. Quantum model of stabilisation and magnetic drift in 2D Bounded states of Kramers-Henneberger well VHK(x, y) in 2D Eigenvalue problem: E0 = 15 a.u. Potential and ground state (E1 = – 0.0192 a.u.). Only excited state (E2 = – 0.0158 a.u.)

30. Classical model in 2D/3D Dipole approxination (regarding y dependence in Ex(y,t) and Bz(y,t)) Characteristics of wavefunction 1. Regardless the phase of laser field magnetic drift is always directed to +y 2. Drift is linear in time + 2w oscillations 3. Constant drift velocity depends on E02 4. Wave of electron finding probability has double frequency in y direction Quantum model of stabilisation and magnetic drift in 2D HK well bound states coupling with continuum Approximated hamiltonian for linear polarised laser field propagating in y direction: (4)

31. Stabilisation and magnetic drift in ultra-strong laser fieldSummary of 2D simulations Probability of finding electron inside the o box with size 100 a.u. around the well with dipole app. without dipole app. Field intensity Ryabikin, Sergeev, Optics Express 417, 7 12 (2000)

32. Saving the stabilisation using constant magnetic field for quantum model and details see Phys Rev A 68, 045401 (2003) Wavefunction without and with constant magnetic field

33. Outline 1. Stabilisation and related phenomena in 1D strong field ionisation 2. Beyond the dipole approximation - magnetic drift in 2D ultra-strong field 3. Brief report on other projects (control of wavefunction by pulses)

34. Outline 1. Stabilisation and related phenomena in 1D strong field ionisation 2. Beyond the dipole approximation - magnetic drift in 2D ultra-strong field 3. Brief report on other projects (control of wavefunction by pulses)

35. Other projects (quantum engineering)

36. Control of ionisation of hydrogen atom by two-colour pulsePhysics Letters A 25 205-211 (1999) detuning

37. Control of Rydberg atom by chirped pulsePhys Rev A 57, 4561 (1998) Rochester atom: Chirped pulse: w0 = w13 = 0.518 a.u. w1 = w34 = 0.0587 a.u. e = 0.05 a.u. By changing the depth of modulation one can control the coupling amplitudes depending on

38. (no modulation) Control of Rydberg atom by chirped pulsePhys Rev A 57, 4561 (1998) Chirped pulse: Photoelectron spectra

39. Control of Rydberg atom by chirped pulsePhys Rev A 57, 4561 (1998) Chirped pulse: Bound states population

40. Outline 1. Stabilisation and related phenomena in 1D strong field ionisation 2. Beyond the dipole approximation - magnetic drift in 2D ultra-strong field 3. Brief report on other projects (control of wavefunction by pulses)