## Abstract

Ultrafast electronic dynamics are typically studied using pulsed lasers. Here we demonstrate a complementary experimental approach: quantum simulation of ultrafast dynamics using trapped ultracold atoms. Counter-intuitively, this technique emulates some of the fastest processes in atomic physics with some of the slowest, leading to a temporal magnification factor of up to 12 orders of magnitude. In these experiments, time-varying forces on neutral atoms in the ground state of a tunable optical trap emulate the electric fields of a pulsed laser acting on bound charged particles. We demonstrate the correspondence with ultrafast science by a sequence of experiments: nonlinear spectroscopy of a many-body bound state, control of the excitation spectrum by potential shaping, observation of sub-cycle unbinding dynamics during strong few-cycle pulses, and direct measurement of carrier-envelope phase dependence of the response to an ultrafast-equivalent pulse. These results establish cold-atom quantum simulation as a complementary tool for studying ultrafast dynamics.

## Introduction

The study of ultrafast-equivalent electronic and vibrational dynamics is a natural but largely unexplored application of cold-atom quantum simulation techniques^{1,2,3,4,5}. Quantum simulation experiments often rely on an analogy between trapped neutral atoms and electrons in matter^{6,7,8}. Although these two systems have vastly different energy densities and constituents which differ in mass and charge, they can often be described by equivalent Hamiltonians, which give rise to equivalent physics. This analogy has been used to explore equilibrium solid-state phenomena from Mott insulators to antiferromagnets^{9, 10}, and dynamical phenomena from Bloch oscillations to many-body localization^{11, 12}. Here we extend this analogy to quantum simulation of ultrafast dynamics, with the aim of realizing an alternate experimental approach to open questions in a vibrant and expanding area of science^{13,14,15}, testing approximate theories^{16,17,18,19,20}, and pushing into experimentally unexplored regimes.

The quantum simulator we describe consists of an artificial atom or molecule made from a trapped quantum gas. The analog of the atomic or molecular binding potential is the tunable AC Stark potential of an optical trap (Fig. 1a), and the analog of the pulsed laser’s electric field is an inertial force arising from rapid trap translation (Fig. 1b). The time-dependent Gross–Pitaevskii equation describing the evolution of the condensate wavefunction *Ψ*(**r**, *t*) is^{1}

where *m* is the atomic mass, *g* = 4π*ħ*^{2}*a*_{s}/*m* parameterizes interactions among *N* atoms, *a*_{s} is the scattering length, and the optical potential *V*(**r**) is shaken in the \({\hat{\mathbf x}}\) direction with waveform *α*(*t*) (see Fig. 1c for a sample waveform). Crucially, the same equation also describes the evolution of the wavefunction of an atomic electron interacting with a linearly polarized laser field in the Kramers–Henneberger frame of reference^{1}, taking *m* to be the electron mass, *V* the nuclear potential, including screening effects, and *g* → 0. For dipolar excitations like those in the experiments we present, the impact of the atoms’ nonzero *g* is minimized due to Kohn’s theorem. Very similar dynamics have been theoretically predicted^{1} for the cold-atom and ultrafast realizations of Eq. (1). A closely related equivalence is described in ref. ^{2}. This equivalence between the evolution of condensate and electron wavefunctions motivates cold-atom quantum simulation of ultrafast dynamics, much as the Bose–Hubbard model motivated early quantum simulation of Mott insulators^{9}.

Though very little of the growing body of quantum simulation work has addressed ultrafast phenomena, a robust toolkit exists for controlling and measuring excitations in trapped gases. Collective excitations in Bose condensates were a major focus of early experimental and theoretical research^{21,22,23,24,25,26}, and the analogy between degenerate trapped gases and individual atoms was noted at that time^{1, 27, 28}. Ultrafast probes have recently been used to study many-body dynamics in Rydberg atoms^{29}, and recent theoretical proposals have suggested the use of cold atoms to simulate ultrafast dynamics in atoms^{2, 5}, molecules^{4}, and solids^{3}.

Cold gases offer unique capabilities for dynamical quantum simulation. Due to the extremely low energy scales, the dynamics are slowed, or magnified, with respect to atomic or molecular timescales by as much as 12 orders of magnitude, allowing the observation of ultrafast-equivalent processes in ultra-slow motion^{2}. This extreme temporal magnification—quantum gas chronoscopy—enables simple and complete control over all parameters of an applied force pulse, as well as straightforward measurement of the artificial atom’s or molecule’s response, with time resolution much faster than all relevant dynamics. The excitation spectrum itself can also be controlled by trap shaping. Using this toolkit of capabilities, we demonstrate experimentally that cold-atom quantum simulation can be used to probe complex phenomena of ultrafast science such as the effect of carrier-envelope phase and pulse intensity on unbinding dynamics. Fig. 1d, e shows examples of experimentally observed sub-cycle unbinding. These experiments demonstrate a new application for cold-atom quantum simulation and establish a potentially fruitful connection between ultrafast and ultracold atomic physics.

## Results

### Ultrafast-equivalent pulse synthesis

The experiments we describe use a Bose condensate of \(N \simeq 20\,000\) atoms of ^{84}Sr^{30}, with a scattering length \(a_{\mathrm{s}} \simeq 6.5{\kern 1pt} {\mathrm{nm}}\). Rapid trap translation gives rise to time-dependent inertial forces designed to have the same approximate functional form, and the same effect of driving dipole-allowed transitions, as the electric field of an ultrafast pulsed laser (see Fig. 1b for diagram of experimental setup). This is achieved by applying a trap, which depends on *x* and *t* as \(V(x,t) = - V_{{\mathrm{trap}}} \times {\mathrm{exp}}[ - 2(x - \alpha (t))^2{\mathrm{/}}w^2]\), where *w* is the 1/e^{2} trap waist and

Control over the pulse is effectively arbitrary; variable parameters include amplitude *A*, carrier frequency *f*, pulse full-width at half-maximum \(\tau = \left( {2{\kern 1pt} {\mathrm{ln}}\left( {2 + \sqrt 3 } \right)} \right){\mathrm{/}}\eta\), and carrier-envelope phase *ϕ* (CEP). The measured trap center translation as a function of time during a typical pulse is shown in Fig. 1c. All data reported here use pulse amplitudes well below the trap width. The effective Keldysh parameter in such an experiment is \(\gamma _{\mathrm{K}} = \sqrt {V_{{\mathrm{trap}}}{\mathrm{/}}2U_{\mathrm{p}}}\), where the optical trap depth *V*_{trap} corresponds to the ionization energy and the ponderomotive potential \(U_{\mathrm{p}} \simeq m\overline {\dot \alpha ^2} {\mathrm{/}}2\) is the time-averaged kinetic energy imparted to the atoms by the pulse. The use of inertial forces enables realization of Keldysh parameters of order unity and greater. Keldysh parameters much <1 could be straightforwardly attained by using a time-varying optical potential gradient rather than trap motion to apply the simulated electric field.

### Spectroscopy of tunable collective excitations

We performed initial spectroscopic characterization of our quantum simulator by applying pulses of constant length much greater than a drive period and variable carrier frequency *f*. After each pulse, the atoms that had not been unbound from the trap were counted with absorption imaging. The resulting plots of bound fraction versus pulse frequency characterize the collective excitation spectra of the trapped condensate. Nonlinear effects are straightforwardly probed by increasing the pulse intensity.

Excitation spectra for one particular trap are shown in Fig. 2a. The resonance at ~750 Hz corresponds to dipole oscillation in the trap and is at the same frequency as the resonance for a non-interacting gas. As the pulse amplitude is increased, higher modes are excited and power broadening is observed. Since our trap is deeply in the Thomas-Fermi regime, these modes are anharmonic and strongly collective. Fig. 2b shows a two-dimensional amplitude-frequency spectrogram.

The excitation spectrum can be tuned by adjusting the trap shape, enabling the study of ultrafast-equivalent dynamics in systems with specific spectral characteristics such as mode degeneracies. The results of such tuning of the excitation spectrum are presented in Fig. 3a–c. We observe good agreement with analytic predictions for dipole-allowed collective resonance positions in the broadened and unbroadened trap^{26}. Note that the frequencies of these complex anharmonic modes are not simply rescaled by broadening, but disperse at different rates; this enables tunable creation of mode degeneracies. This tunability of the collective excitation spectra is a key feature of cold-atom-based quantum simulation of ultrafast dynamics. Static adjustments like those demonstrated here enable the realization of desired spectral properties, and rapid tuning of mode degeneracies could enable the study of controllably diabatic or adiabatic dynamics. Future experiments could use this ability for quantum simulation of molecular energy relaxation mechanisms in the vicinity of tunable mode degeneracies similar to conical intersections^{31}.

### Momentum-resolved sub-cycle unbinding dynamics

Having demonstrated quasi-continuous-wave spectroscopy of bound states with tunable energy spectra, we turn to the use of this tool for quantum simulation of ultrafast dynamics during few-cycle pulses^{32}. In ultrafast streaking measurements, the electric field of a few-cycle femtosecond pulse deflects photoelectrons produced by an attosecond extreme ultraviolet pulse striking an atom, allowing characterization of both the pulses and the atom^{13, 33, 34}. In the quantum simulator, qualitatively similar techniques allow high-resolution measurement of sub-cycle quantum dynamics. Here, instead of using photoionization to terminate the dynamics, the experimenter can simply instantaneously turn off the trapping potential at any point before, during, or after the pulse. The atoms then propagate freely in space, and their instantaneous momenta at the time of trap removal are mapped onto their positions after some time of flight. Varying the time at which the trap is removed enables measurement of the time evolution of the bound quantum system with time resolution far below a drive period. This experimental technique, while commonplace in ultracold atomic physics, represents a powerful and general tool for the study of ultrafast-equivalent dynamics in our quantum simulator.

Fig. 4a–e presents the results of such measurements for both off-resonant and near-resonant pulses. The Bose-condensed atoms initially occupy mainly a single eigenstate of the transverse trapping potential. Quantum dynamics during and after the pulses shown in Fig. 4a, c can be tracked by direct momentum-space imaging of the atoms. Fig. 4b, d shows the density distribution after time of flight, integrated over the directions transverse to the excitation, as a function of time. For a pulse carrier frequency significantly below the dipole oscillation frequency *ν*_{
x
} in the dimension of driving (see Fig. 4a for pulse waveform), the momentum of the Bose–Einstein condensate (BEC) evolves coherently during and after the pulse, as shown in Fig. 4b. Incoherent heating due to the pulse is observed to be minimal on the few-cycle timescales we probe. The atoms respond to the pulse at *ν*_{
x
}—a higher frequency than the carrier—but remain bound. During a pulse with carrier frequency near *ν*_{
x
}, however, qualitatively different behavior is observed. Fig. 4d shows momentum evolution during a near-resonant pulse for an amplitude near the unbinding threshold (see Fig. 4c for pulse waveform). In this parameter regime, atoms do not leave the trap all at once, but do not incoherently heat either; instead, ejection starts at the time of the pulse peak, with additional bursts of atoms emitted during each subsequent half-cycle of the pulse. Fig. 4e shows one such burst. This unbinding process models ionization or molecular disintegration during an ultrafast laser pulse.

### Varying pulse amplitude and carrier-envelope phase

The ability to precisely measure the population and momenta of unbound states as a function of pulse parameters and time opens up the possibility of flexible quantum simulation of ultrafast unbinding dynamics. As an initial application of the quantum simulator presented herein we have measured the dependence of simulated ionization yield or photodissociation on both pulse amplitude and carrier-envelope phase. This represents a complementary method of testing the effects of two parameters central to numerous experimental and theoretical studies of ultrafast multiphoton ionization and bond-breaking processes^{16,17,18, 35,36,37}.

Both the precise unbinding time during an applied force pulse and the final unbound momentum depend sensitively and non-monotonically on pulse amplitude. In the quantum simulator, the amplitude of the pulse can be straightforwardly varied over a wide range, keeping the carrier-envelope phase, carrier frequency, and total pulse time constant. As the amplitude is increased from that used in Fig. 4d, the unbinding dynamics change drastically. Fig. 5a shows the momentum distribution of the atoms (measured by detecting the position distribution after 2 ms time of flight) after few-cycle pulses with amplitudes from 0 up to 6 μm. Figure 5b–e shows the full-time evolution of the momentum distribution during few-cycle pulses of selected amplitudes. Below a critical amplitude, no atoms are ejected from the trap. For some intermediate amplitudes, the behavior mirrors that shown in Fig. 4d, with bursts of atoms unbinding at different points during the pulse. Above that intermediate regime, all of the atoms unbind at one well-defined time and continue to move with constant momentum after unbinding. Strikingly, as the amplitude is increased further, the momentum of the unbound atoms reverses sign, as they unbind half a drive cycle earlier, in an oppositely directed simulated electric field.

Even for fixed pulse amplitude, the final state of the initially bound system after the force pulse depends sensitively and non-trivially on the CEP. In the pulsed-laser experimental context, the advent of few-cycle pulses with adjustable, stabilized CEP^{38} has enabled advances such as probes of the effects of CEP on ultrafast dynamics^{39, 40}, study of interference patterns in multiparticle ionization signals^{41}, and control of recollision processes in molecular ions^{42, 43}. The nearly arbitrary pulse shape control available in the cold-atom quantum simulator makes it a flexible tool for probing the dependence of ultrafast-equivalent dynamics on CEP. Figure 6a shows a measurement of post-pulse momentum distribution (again detected via time of flight) as a function of the carrier-envelope phase of a near-resonant applied pulse. Changing the CEP from 0 to π flips the sign of all forces during the pulse, and results in inversion of the momentum of the unbound atoms. As shown in Fig. 6b–d, the pulses at integer values of *ϕ*/π have sine-like character, possessing odd symmetry under reflection in time around the pulse center. Pulses with a CEP of 3π/2 have cosine-like character and give rise to very different unbinding dynamics at this pulse amplitude, populating more than one momentum class of unbound atoms. More complex dynamical phenomena are also visible in Fig. 6: the inward slope of the unbound momentum as CEP increases in the neighborhood of *ϕ* = π can be understood as the consequence of the force at the first unbinding peak sliding down the pulse envelope, and the observed asymmetry between *ϕ* = π/2 and *ϕ* = 3π/2 indicates a violation of inversion symmetry in the potential. The most likely cause of this symmetry-breaking is slight trap aberration; this points the way to future work elucidating the effects of potential shape on unbinding dynamics.

## Discussion

The results presented here open the door to a broad class of quantum simulation experiments investigating ultrafast nonequilibrium phenomena, with numerous possible scientific targets. Emulation of pump-probe experiments, multichromatic light fields, and non-physical (for example, half-cycle) pulse shapes impossible to create with lasers would require no techniques beyond those demonstrated here apart from changing the form of *α*(*t*).

Additional possibilities require only modest extensions of the experimental approach reported here. The simplest such extension would be to replace the inertial forces used to emulate electric fields with time-varying Zeeman or Stark potential gradients. This would allow ultrafast quantum simulation in the regime of Keldysh parameter *γ*_{K} <1, and enable the direct experimental investigation of additional open questions of current interest in ultrafast science. Potential scientific targets include the creation of photoelectron vortices with circularly polarized pulses^{44}, the demonstration and study of strong-field stabilization, wherein the ionization probability becomes a decreasing function of pulse amplitude^{45}, and detailed quantitative measurement of sub-cycle tunnel ionization timing effects^{46,47,48,49}.

Other extensions to the basic technique are also possible. The use of traps with multiple minima could enable modeling of more complex molecular configurations^{4}. Ultrafast quantum simulation could also be pursued with small numbers of trapped fermions^{50}, making a more direct analog of atomic electrons. However, the use of Bose condensates and the analogy of Eq. (1) does greatly magnify the signal, making experiments feasible with bosons that would be very challenging with fermions. Finally, an expansion of the analogy underlying these quantum simulation experiments beyond atoms and molecules could enable the study of ultrafast-equivalent dynamical phenomena relevant to nuclear excitations^{51} and strong-field dynamics in solids^{3}.

In summary, we present experimental results from a cold-atom quantum simulator of ultrafast phenomena, including nonlinear spectroscopy of collective excitations, control of the energy spectra of the bound states of the simulator, imaging of sub-cycle dynamics during an unbinding process similar to ultrafast ionization, and measurement of the effects on unbinding dynamics of pulse amplitude and carrier-envelope phase. Such cold-atom quantum simulation of ultrafast dynamical phenomena has the potential to enable benchmarking of relevant theories and explorations of experimentally challenging regimes, in an approach complementary to both ultrafast theory and pulsed-laser experiments.

## Methods

### Preparation of a BEC in an optical trap

During the computer-controlled experimental sequence^{52}, atoms from an effusive source are collimated by a nozzle^{53}, Zeeman-slowed, trapped, and pre-cooled by sequential magneto-optical traps using the 461 and 689 nm ground-state transitions, and evaporatively cooled to degeneracy in a crossed-beam 1064 nm optical dipole trap (ODT)^{30}. In the emulation stage of the experiment, the resulting condensate of 2 × 10^{4 84}Sr atoms is adiabatically loaded into a single-beam ODT with a waist of 15 μm. The trap depth can be varied across a wide range; a typical value used in the work presented here is 10 μK. The s-wave scattering length of ^{84}Sr is 6.5 nm.

### Temporal magnification of the quantum simulator

Varying the power and shape of the trap beam yields transverse trap frequencies *ν*_{
x
} and *ν*_{
z
} between 300 and 1000 Hz, and an axial trap frequency *ν*_{
y
} between 5 and 15 Hz. Drawing the analogy between the optical trap and a single hydrogen atom, the energy difference *hν*_{
x
} between the ground and first relevant excited states (typically \(\simeq\) 500 Hz) is analogous to the 2.47 PHz Ly-*α* line. Although the detailed energy spectrum of the ODT differs from that of hydrogen due to the different potential shapes, comparison of these two frequency scales indicates an approximate temporal magnification factor of 5 × 10^{12}. Emulation of molecular vibrational excitations or cluster dynamics leads to a ratio of up to 10^{11} between characteristic timescales of the emulator and emulated system.

### Application of time-varying inertial forces

An acousto-optic modulator (AOM) can translate the trap center in the *x*-direction (see Fig. 1b for axis definitions) at amplitudes up to 6 μm and frequencies from DC up to hundreds of kHz. Because of the direction of translation and the wide separation between transverse and axial frequency scales, the axial degree of freedom is irrelevant to the results we present. Under the assumption that the atoms remain near the center of the trap during the pulse, which we observe to be true until unbinding, the effective applied force is \(m\ddot x = F(t) = \left. { - {\mathrm{d}}V(x,t){\mathrm{/}}{\mathrm{d}}x} \right|_{x = 0}\). \(\ddot x\) can be specified as desired; for the particular functional form chosen, \(\ddot x(t)\) has a similar enveloped-pulse shape to *x*(*t*).

### Tuning the trap geometry

Shaping of the trap is achieved most simply by periodic translation on much faster timescales than the dynamics of the BEC, so that the atoms experience a time-averaged potential. For this purpose, the trap AOM’s radio-frequency drive frequency was sinusoidally modulated at 500 kHz, giving rise to a maximum trap translation amplitude of 4.2 ± 0.3 μm.

### Data availability

The data that support the findings in this study are available from the corresponding author upon reasonable request.

## References

- 1.
Dum, R. et al. Wave packet dynamics with Bose-Einstein condensates.

*Phys. Rev. Lett.***80**, 3899–3902 (1998). - 2.
Sala, S., Förster, J. & Saenz, A. Ultracold-atom quantum simulator for attosecond science.

*Phys. Rev. A***95**, 011403 (2017). - 3.
Arlinghaus, S. & Holthaus, M. Driven optical lattices as strong-field simulators.

*Phys. Rev. A***81**, 063612 (2010). - 4.
Lühmann, D.-S., Weitenberg, C. & Sengstock, K. Emulating molecular orbitals and electronic dynamics with ultracold atoms.

*Phys. Rev. X***5**, 031016 (2015). - 5.
Rajagopal, S. V. et al. Quantum emulation of extreme non-equilibrium phenomena with trapped atoms.

*Ann. Phys.***529**, 1700008 (2017). - 6.
Lewenstein, M. et al. Ultracold atomic gases in optical lattices: mimicking condensed matter physics and beyond.

*Adv. Phys.***56**, 243–379 (2007). - 7.
Bloch, I., Dalibard, J. & Nascimbene, S. Quantum simulations with ultracold quantum gases.

*Nat. Phys.***8**, 267–276 (2012). - 8.
Gross, C. & Bloch, I. Quantum simulations with ultracold atoms in optical lattices.

*Science***357**, 995–1001 (2017). - 9.
Greiner, M., Mandel, O., Esslinger, T., Hänsch, T. W. & Bloch, I. Quantum phase transition from a superfluid to a mott insulator in a gas of ultracold atoms.

*Nature***415**, 39–44 (2002). - 10.
Hart, R. A. et al. Observation of antiferromagnetic correlations in the hubbard model with ultracold atoms.

*Nature***519**, 211–214 (2015). - 11.
Ben Dahan, M., Peik, E., Reichel, J., Castin, Y. & Salomon, C. Bloch oscillations of atoms in an optical potential.

*Phys. Rev. Lett.***76**, 4508–4511 (1996). - 12.
Schreiber, M. et al. Observation of many-body localization of interacting fermions in a quasirandom optical lattice.

*Science***349**, 842–845 (2015). - 13.
Krausz, F. & Ivanov, M. Attosecond physics.

*Rev. Mod. Phys.***81**, 163–234 (2009). - 14.
Corkum, P. B. & Krausz, F. Attosecond science.

*Nat. Phys.***3**, 381–387 (2007). - 15.
Fennel, T. et al. Laser-driven nonlinear cluster dynamics.

*Rev. Mod. Phys.***82**, 1793–1842 (2010). - 16.
Keldysh, L. V. Ionization in the field of a strong electromagnetic wave.

*J. Exp. Theor. Phys.***47**, 1945–1957 (1964). - 17.
Faisal, F. H. M. Multiphoton transitions. IV. Bound-free transition integrals in compact forms.

*J. Phys. B: At. Mol. Phys.***6**, 553 (1973). - 18.
Reiss, H. R. Gauges for intense-field electrodynamics.

*Phys. Rev. A***22**, 770–772 (1980). - 19.
Corkum, P. B. Plasma perspective on strong field multiphoton ionization.

*Phys. Rev. Lett.***71**, 1994–1997 (1993). - 20.
Lewenstein, M., Balcou, P., Ivanov, M. Y., L’Huillier, A. & Corkum, P. B. Theory of high-harmonic generation by low-frequency laser fields.

*Phys. Rev. A***49**, 2117–2132 (1994). - 21.
Jin, D. S., Ensher, J. R., Matthews, M. R., Wieman, C. E. & Cornell, E. A. Collective excitations of a Bose-Einstein condensate in a dilute gas.

*Phys. Rev. Lett.***77**, 420–423 (1996). - 22.
Mewes, M.-O. et al. Collective excitations of a Bose-Einstein condensate in a magnetic trap.

*Phys. Rev. Lett.***77**, 988–991 (1996). - 23.
Fort, C. et al. Collective excitations of a

^{87}Rb Bose condensate in the Thomas-Fermi regime.*EPL (Europhys. Lett.)***49**, 8 (2000). - 24.
Stringari, S. Collective excitations of a trapped Bose-condensed gas.

*Phys. Rev. Lett.***77**, 2360–2363 (1996). - 25.
Öhberg, P. et al. Low-energy elementary excitations of a trapped Bose-condensed gas.

*Phys. Rev. A***56**, R3346–R3349 (1997). - 26.
Csordás, A. & Graham, R. Collective excitations in Bose-Einstein condensates in triaxially anisotropic parabolic traps.

*Phys. Rev. A***59**, 1477–1487 (1999). - 27.
Esry, B. D. Hartree-Fock theory for Bose-Einstein condensates and the inclusion of correlation effects.

*Phys. Rev. A***55**, 1147–1159 (1997). - 28.
Walsworth, R. & You, L. Selective creation of quasiparticles in trapped Bose condensates.

*Phys. Rev. A***56**, 555–559 (1997). - 29.
Takei, N. et al. Direct observation of ultrafast many-body electron dynamics in an ultracold Rydberg gas.

*Nat. Commun.***7**, 13449 (2016). - 30.
Stellmer, S., Tey, M. K., Huang, B., Grimm, R. & Schreck, F. Bose-Einstein condensation of strontium.

*Phys. Rev. Lett.***103**, 200401 (2009). - 31.
Yarkony, D. R. Conical intersections: the new conventional wisdom.

*J. Phys. Chem. A***105**, 6277–6293 (2001). - 32.
Milošević, D. B., Paulus, G. G., Bauer, D. & Becker, W. Above-threshold ionization by few-cycle pulses.

*J. Phys. B: At. Mol. Opt. Phys.***39**, R203 (2006). - 33.
Hentschel, M. et al. Attosecond metrology.

*Nature***414**, 509–513 (2001). - 34.
Itatani, J. et al. Attosecond streak camera.

*Phys. Rev. Lett.***88**, 173903 (2002). - 35.
Chin, S. L. & Isenor, N. R. Multiphoton ionization in atomic gases with depletion of neutral atoms.

*Can. J. Phys.***48**, 1445–1447 (1970). - 36.
Mainfray, G. & Manus, G. Multiphoton ionization of atoms.

*Rep. Progress Phys.***54**, 1333 (1991). - 37.
l’Huillier, A., Lompre, L. A., Mainfray, G. & Manus, C. Multiply charged ions induced by multiphoton absorption in rare gases at 0.53

*μ*m.*Phys. Rev. A***27**, 2503–2512 (1983). - 38.
Baltuška, A. et al. Attosecond control of electronic processes by intense light fields.

*Nature***421**, 611–615 (2003). - 39.
Peng, L.-Y. & Starace, A. F. Attosecond pulse carrier-envelope phase effects on ionized electron momentum and energy distributions.

*Phys. Rev. A***76**, 043401 (2007). - 40.
Kling, M. F. et al. Imaging of carrier-envelope phase effects in above-threshold ionization with intense few-cycle laser fields.

*New J. Phys.***10**, 25024–25024 (2008). - 41.
Krüger, M., Schenk, M. & Hommelhoff, P. Attosecond control of electrons emitted from a nanoscale metal tip.

*Nature***475**, 78–81 (2011). - 42.
Kling, M. F. et al. Control of electron localization in molecular dissociation.

*Science***312**, 246–8 (2006). - 43.
Rathje, T. et al. Coherent control at its most fundamental: carrier-envelope-phase-dependent electron localization in photodissociation of a H molecular ion beam target.

*Phys. Rev. Lett.***111**, 093002 (2013). - 44.
Ngoko Djiokap, J. M. et al. Electron vortices in photoionization by circularly polarized attosecond pulses.

*Phys. Rev. Lett.***115**, 113004 (2015). - 45.
Pont, M. & Gavrila, M. Stabilization of atomic hydrogen in superintense, high-frequency laser fields of circular polarization.

*Phys. Rev. Lett.***65**, 2362–2365 (1990). - 46.
Shafir, D. et al. Resolving the time when an electron exits a tunnelling barrier.

*Nature***485**, 343–346 (2012). - 47.
Pfeiffer, A. N. et al. Attoclock reveals natural coordinates of the laser-induced tunnelling current flow in atoms.

*Nat. Phys.***8**, 76–80 (2012). - 48.
Uiberacker, M. et al. Attosecond real-time observation of electron tunnelling in atoms.

*Nature***446**, 627–632 (2007). - 49.
Isinger, M. et al. Photoionization in the time and frequency domain.

*Science***358**, 893–896 (2017). - 50.
Serwane, F. et al. Deterministic preparation of a tunable few-fermion system.

*Science***332**, 336–338 (2011). - 51.
Zinner, N. T. & Jensen, A. S. Comparing and contrasting nuclei and cold atomic gases.

*J. Phys. G: Nucl. Part. Phys.***40**, 053101 (2013). - 52.
Keshet, A. & Ketterle, W. A distributed, graphical user interface based, computer control system for atomic physics experiments.

*Rev. Sci. Instrum.***84**, 015105 (2013). - 53.
Senaratne, R. et al. Effusive atomic oven nozzle design using an aligned microcapillary array.

*Rev. Sci. Instrum.***86**, 023105 (2015).

## Acknowledgements

The authors thank Vyacheslav Lebedev, Mikhail Lipatov, Ethan Q. Simmons, Yi Zeng, Jacob Hines, Ian Harley-Trochimczyk, James Chow, and Bryance Oyang for experimental assistance; thank Andrew Jayich and Alejandro Saenz for critical readings of the manuscript; and acknowledge support from the National Science Foundation (CAREER 1555313), Air Force Office of Scientific Research (YIP FA9550-12-1-0305), Army Research Office (PECASE W911NF1410154 and DURIP W911NF1510436), and a President’s Research Catalyst Award (CA-15-327861) from the UC Office of the President.

## Author information

### Affiliations

### Contributions

All authors contributed extensively to the work presented in this paper. R.S., S.V.R., T.S., and P.E.D. performed the measurements, with experimental assistance from K.M.F., K.S., and Z.A.G. R.S., S.V.R., T.S., P.E.D., and D.M.W. analyzed the data and wrote the manuscript. All authors contributed insights and discussions on the experiment. D.M.W. supervised the project.

### Corresponding author

## Ethics declarations

### Competing interests

The authors declare no competing interests.

## Additional information

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

## About this article

### Cite this article

Senaratne, R., Rajagopal, S.V., Shimasaki, T. *et al.* Quantum simulation of ultrafast dynamics using trapped ultracold atoms.
*Nat Commun* **9, **2065 (2018). https://doi.org/10.1038/s41467-018-04556-3

Received:

Accepted:

Published:

## Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.