---
abstract: 'Ultrafast electronic dynamics are typically studied using pulsed lasers. 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 twelve 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.'
author:
- Ruwan Senaratne
- 'Shankari V. Rajagopal'
- Toshihiko Shimasaki
- 'Peter E. Dotti'
- 'Kurt M. Fujiwara'
- Kevin Singh
- 'Zachary A. Geiger'
- 'David M. Weld'
title: Quantum Simulation of Ultrafast Dynamics Using Trapped Ultracold Atoms
---

[^1]

[^2]

The study of ultrafast-equivalent electronic and vibrational dynamics is a natural but largely unexplored application of cold-atom quantum simulation techniques [@lewenstein-wavepacketdynamics; @saenz-attosecondsimulatorPRA; @holthaus-strongfieldsim; @sengstock-moleculesimulator; @adppaper]. Quantum simulation experiments often rely on an analogy between trapped neutral atoms and electrons in matter [@lewenstein-review; @blochdalibardQSreview; @bloch2017review]. 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 [@greinerSFMI; @HuletAF], and dynamical phenomena from Bloch oscillations to many-body localization [@salomon-blochoscs; @bloch-mblscience]. 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 [@krausz-attosecond-review; @corkumkrauszreview; @clusters], testing approximate theories [@Keldysh1964; @Faisal1973; @Reiss1980; @corkum-model; @lewenstein-hhg], and pushing into experimentally unexplored regimes.

![**Quantum simulation of ultrafast dynamics.** **a:** Schematic bound states of Lennard-Jones, Gaussian, and $1/r$ potentials (offset for clarity). **b:** Diagram of optical trap (red), which is shaken in the $\mathbf{\hat{x}}$ direction to generate inertial forces on the condensate (blue). **c:** Measured trap position $\alpha(t)$ during a pulse. **d:** Response to a weak pulse. Colourmap shows density distribution after time-of-flight as a function of time. Pulse carrier frequency is 450 Hz, pulse envelope width is 3.76 ms, pulse amplitude is 0.6 $\upmu$m, and carrier-envelope phase is 0, as defined in Eq. \[eqn:pulse\]. **e:** Response to a stronger pulse. Unbinding occurs near 8 ms, after which the atoms propagate with constant velocity. Pulse amplitude is 2.4 $\upmu$m. All other pulse parameters are identical to those in **d**.[]{data-label="fig:summary"}](fig1.pdf){width="\columnwidth"}

The quantum simulator we describe consists of an artificial atom or molecule made from a trapped quantum gas. The analogue of the atomic or molecular binding potential is the tunable AC Stark potential of an optical trap, and the analogue of the pulsed laser’s electric field is an inertial force arising from rapid trap translation. The time-dependent Gross-Pitaevskii equation describing the evolution of the condensate wavefunction $\mathit{\Psi}(\boldsymbol{\mathrm{r}},t)$ is [@lewenstein-wavepacketdynamics] $$\left[-\mathrm{i}\hbar\partial_t -\frac{\hbar^2\nabla^2}{2m}+V(\boldsymbol{\mathrm{r}}+\alpha(t)\boldsymbol{\mathrm{\hat{x}}})+gN|\mathit{\Psi|}^2\right]\mathit{\Psi}=0,
\label{GPE}$$ where $m$ is the atomic mass, $g=4\mathrm{\pi}\hbar^2a_{\mathrm{s}}/m$ parameterizes interactions among $N$ atoms, $a_{\mathrm{s}}$ is the scattering length, and the optical potential $V(\boldsymbol{\mathrm{r}})$ is shaken in the $\boldsymbol{\mathrm{\hat{x}}}$ direction with waveform $\alpha(t)$. 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 [@lewenstein-wavepacketdynamics], taking $m$ to be the electron mass, $V$ the nuclear potential including screening effects, and $g\rightarrow 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 [@lewenstein-wavepacketdynamics] for the cold-atom and ultrafast realizations of Eq. \[GPE\]. A closely related equivalence is described in [@saenz-attosecondsimulatorPRA]. 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 [@greinerSFMI].

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 [@jin-collectiveexc; @ketterle-collectiveexcitations; @inguscioexcitations; @stringari-orig; @shlyapnikov-excitations; @graham-anisotropictrapexc], and the analogy between degenerate trapped gases and individual atoms was noted at that time [@esry97; @walsworth97; @lewenstein-wavepacketdynamics]. Ultrafast probes have recently been used to study many-body dynamics in Rydberg atoms [@ultrafastrydberg], and recent theoretical proposals have suggested the use of cold atoms to simulate ultrafast dynamics in atoms [@saenz-attosecondsimulatorPRA; @adppaper], molecules [@sengstock-moleculesimulator], and solids [@holthaus-strongfieldsim].

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 twelve orders of magnitude, allowing the observation of ultrafast-equivalent processes in ultra-slow-motion [@saenz-attosecondsimulatorPRA]. 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. 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 [@strontium84BEC], with a scattering length $a_{\mathrm{s}}\!\simeq\!6.5~\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. This is achieved by applying a trap which depends on $x$ and $t$ as , where $w$ is the 1/$\mathrm{e}^2$ trap waist and $$\alpha(t) = A\  {\rm sech} \left[ \eta \left( t - t_0 \right) \right]  \sin \left[ 2 \pi f \left( t - t_0 \right) + \phi  + \mathrm{\pi} \right].
\label{eqn:pulse}$$ Control over the pulse is effectively arbitrary; variable parameters include amplitude $A$, carrier frequency $f$, pulse full-width at half-maximum $\tau=(2 \ln{(2+\sqrt{3})})/\eta$, and carrier-envelope phase $\phi$. The measured trap centre translation as a function of time during a typical pulse is shown in Fig. \[fig:summary\]c. 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_\text{trap}/2U_{\mathrm{p}}}$, where the optical trap depth $V_\text{trap}$ corresponds to the ionization energy and the ponderomotive potential $U_{\mathrm{p}} \simeq  m \overline{\dot{\alpha}^2}/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 less than 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 the quantum simulator.** **a:** Remaining bound fraction as a function of carrier frequency for pulses with $\tau =250$ ms, $\phi=\mathrm{\pi}$ and amplitudes of 0.9 $\upmu$m (dotted) and 2.4 $\upmu$m (solid). Note the emergence of higher-order peaks and power broadening at larger amplitudes. **b:** Unbound fraction as a function of applied pulse frequency and amplitude, for a 250 ms pulse. Lines indicate cuts plotted in panel **a**. []{data-label="fig:spec"}](fig2.pdf){width="0.95\linewidth"}

**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 the top panel of Fig. \[fig:spec\]. The resonance at $\sim$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. The bottom panel of Fig. \[fig:spec\] shows a 2D amplitude-frequency spectrogram.

![**Tunable excitation spectra via trap shaping.** **a:** Bound fraction after a 1-second pulse as a function of excitation frequency for an unmodified trap. Bold dotted lines are theoretically predicted frequencies of collective resonances expected to couple to our drive. Labels on theory lines indicate the quantum numbers $(k,\beta,\gamma)$, using the notation of [@graham-anisotropictrapexc]. $\beta$ and $\gamma$ are parity quantum numbers, and $k$ indicates the form of the nodal surface for the excitation. The quantum number $m$ is 1 for all resonances plotted. The only inputs to this theory are the three trap frequencies. The resonance at half the fundamental frequency is believed to be due to parametric excitation of a dipole oscillation in the direction of gravity. Pulse amplitudes were increased from 0.6 $\upmu$m at low frequency to 3 $\upmu$m at the highest frequency to maximize peak visibility. **b:** Evolution of predicted resonances under continuously increasing trap broadening. Thinner dotted lines represent resonances which are not dipole-allowed for this drive polarization. **c:** Bound fraction after a 1-second pulse as a function of excitation frequency for a trap broadened in one direction as described in the methods section. []{data-label="fig:broadening"}](fig3.pdf){width="0.95\linewidth"}

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. \[fig:broadening\]. We observe good agreement with analytic predictions for dipole-allowed collective resonance positions in the broadened and unbroadened trap [@graham-anisotropictrapexc]. 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 [@CI-2001].

![**Sub-cycle dynamics during off-resonant and near-resonant pulses.** **a:** Trap minimum position as a function of time for a pulse with $\tau = 3.76$ ms, $f = 200$ Hz, $A = 3 \upmu$m and $\phi = \mathrm{\pi}$. **b:** Post-time-of-flight integrated spatial density distribution versus trap turn-off time during the off-resonant pulse depicted in panel **a**. Here $\nu_{x}=450$ Hz. **c:** Trap minimum position as a function of time for a pulse with $\tau = 3.76$ ms, $f = 550$ Hz, $A = 1.5 \upmu$m and $\phi = \frac{3 \mathrm{\pi}}{2}$. **d:** Post-time-of-flight integrated spatial density distribution versus trap turn-off time during the near-resonant pulse depicted in panel **c**. Here $\nu_{x}=600$ Hz. **e:** Density distribution at the time indicated by the dashed line in panel **d**. Peaks from bound and ejected atoms are visible.[]{data-label="fig:streak"}](fig4.pdf){width="0.95\linewidth"}

![image](fig5.pdf){width="0.75\linewidth"}

**Momentum-resolved sub-cycle unbinding dynamics.** Having demonstrated quasi-CW 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 [@fewcycleionization-zerorange]. 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 [@Hentschel2001; @Itatani2002; @Krausz2009]. 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. \[fig:streak\] 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 pulse can be tracked by direct momentum-space imaging of the atoms. Panels **b** and **d** of Fig. \[fig:streak\] show 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 $\nu_{x}$ in the dimension of driving, the momentum of the BEC evolves coherently during and after the pulse, as shown in Fig. \[fig:streak\]b. Incoherent heating due to the pulse is observed to be minimal on the few-cycle time scales we probe. The atoms respond to the pulse at $\nu_{x}$ — a higher frequency than the carrier — but remain bound. During a pulse with carrier frequency near $\nu_{x}$, however, qualitatively different is observed. Fig. \[fig:streak\]d shows momentum evolution during a near-resonant pulse for an amplitude near the unbinding threshold. 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 \[fig:streak\]e 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 [@Keldysh1964; @Faisal1973; @Reiss1980; @Chin1970; @Mainfray1984; @LHuillier1983].

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. \[fig:streak\]d, the unbinding dynamics change drastically. Fig. \[fig:amp\]a 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 $\upmu$m. The bottom panels show 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 behaviour mirrors that shown in Fig. \[fig:streak\]d, 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 carrier-envelope phase $\phi$ (CEP). In the pulsed-laser experimental context, the advent of few-cycle pulses with adjustable, stabilized CEP [@Baltuska2003] has enabled advances such as probes of the effects of CEP on ultrafast dynamics [@Peng2007; @Kling2008], study of interference patterns in multiparticle ionization signals [@Kruger2011], and control of recollision processes in molecular ions [@Kling2006; @Rathje2013]. 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. Fig. \[fig:CEP1\] 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 $\mathrm{\pi}$ flips the sign of all forces during the pulse, and results in inversion of the momentum of the unbound atoms. As shown in the bottom panel of Fig. \[fig:CEP1\], the pulses at integer values of $\phi/\mathrm{\pi}$ have sine-like character, possessing odd symmetry under reflection in time around the pulse centre. Pulses with a CEP of $3\mathrm{\pi}/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. \[fig:CEP1\]: the inward slope of the unbound momentum as CEP increases in the neighborhood of $\phi=\mathrm{\pi}$ can be understood as the consequence of the force at the first unbinding peak sliding down the pulse envelope, and the observed asymmetry between $\phi=\mathrm{\pi}/2$ and $\phi=3\mathrm{\pi}/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.

![**Carrier-envelope phase dependence of final momentum.** **a:** Integrated spatial density distribution after application of a near-resonant 450 Hz pulse with $\tau=3.76$ ms and $A=2.4$ $\upmu$m followed by 3 ms of time-of-flight, versus CEP. **b-d:** Pulse waveforms (force versus time) for CEP values indicated in the inset.[]{data-label="fig:CEP1"}](fig6.pdf){width="0.99\columnwidth"}

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 $\alpha(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 $\gamma_\mathrm{K}$ less than one, and enable the direct experimental investigation of open questions of current interest in ultrafast science. Potential scientific targets include the creation of photoelectron vortices with circularly-polarized pulses [@starace-polarization], the demonstration and study of strong-field stabilization, wherein the ionization probability becomes a decreasing function of pulse amplitude [@stronfieldstabilization], and detailed quantitative measurement of sub-cycle tunnel ionization timing effects [@dudovich-tunnelingtime; @attoclock-tunnelingtiming; @krausz-tunnelingtime; @lindrothhuillierphotoionization].

Other extensions to the basic technique are also possible. The use of traps with multiple minima could enable modelling of more complex molecular configurations [@sengstock-moleculesimulator]. Ultrafast quantum simulation could also be pursued with small numbers of trapped fermions [@jochim-fewfermion], making a more direct analogue of atomic electrons. However, the use of Bose condensates and the analogy of Eq. \[GPE\] 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 [@zinner-nuclearQS] and strong-field dynamics in solids [@holthaus-strongfieldsim].

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 Bose-Einstein condensate in an optical trap.** During the computer-controlled experimental sequence [@Cicero], atoms from an effusive source are collimated by a nozzle [@nozzleRSI], Zeeman-slowed, trapped and pre-cooled by sequential magneto-optical traps using the 461 nm and 689 nm ground-state transitions, and evaporatively cooled to degeneracy in a crossed-beam 1064 nm optical dipole trap (ODT) [@strontium84BEC]. In the emulation stage of the experiment, the resulting condensate of $2\times10^4$ $^{84}$Sr atoms is adiabatically loaded into a single-beam ODT with a waist of 15 $\upmu$m. The trap depth can be varied across a wide range; a typical value used in the work presented here is 10 $\upmu$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 $\nu_x$ and $\nu_z$ between $300$ and $1000$ Hz, and an axial trap frequency $\nu_y$ between 5 and 15 Hz. Drawing the analogy between the optical trap and a single hydrogen atom, the energy difference $h\nu_x$ between the ground and first relevant excited states (typically $\simeq$500 Hz) is analogous to the 2.47 PHz Ly-$\alpha$ 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$\times 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 centre in the $x$-direction (see Fig. \[fig:summary\]c for axis definitions) at amplitudes up to 6 $\upmu$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 centre of the trap during the pulse, which we observe to be true until unbinding, the effective applied force is . $\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 time scales than the dynamics of the BEC, so that the atoms experience a time-averaged potential. For this purpose, the trap AOM’s RF drive frequency was sinusoidally modulated at 500 kHz, giving rise to a maximum trap translation amplitude of 4.2$\pm$0.3 $\upmu$m.

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

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 contributions {#author-contributions .unnumbered}
====================

All authors contributed extensively to the work presented in this paper. RS, SVR, TS, and PED performed the measurements, with experimental assistance from KMF, KS, and ZAG. RS, SVR, TS, PED and DMW analysed the data and wrote the manuscript. All authors contributed insights and discussions on the experiment. DMW supervised the project.

Additional Information {#additional-information .unnumbered}
======================

**Competing financial interests:** The authors declare no competing financial interests.

[54]{}ifxundefined \[1\][ ifx[\#1]{} ]{}ifnum \[1\][ \#1firstoftwo secondoftwo ]{}ifx \[1\][ \#1firstoftwo secondoftwo ]{}““\#1””@noop \[0\][secondoftwo]{}sanitize@url \[0\][‘\
12‘\$12 ‘&12‘\#12‘12‘\_12‘%12]{}@startlink\[1\]@endlink\[0\]@bib@innerbibempty **, [****, ()](\doibase 10.1103/PhysRevLett.80.3899) **, [****, ()](\doibase 10.1103/PhysRevA.95.011403) **, [****, ()](\doibase 10.1103/PhysRevA.81.063612) **, [****, ()](\doibase
  10.1103/PhysRevX.5.031016) **, [****, ()](\doibase 10.1002/andp.201700008) **, [****, ()](\doibase 10.1080/00018730701223200) **, [****, ()](http://dx.doi.org/10.1038/nphys2259) **, [****, ()](\doibase 10.1126/science.aal3837) **, [****, ()](\doibase 10.1038/415039a) **, [****, ()](http://dx.doi.org/10.1038/nature14223) **, [****, ()](\doibase
  10.1103/PhysRevLett.76.4508) **, [****, ()](\doibase 10.1126/science.aaa7432) **, [****, ()](\doibase
  10.1103/RevModPhys.81.163) **, [****, ()](http://dx.doi.org/10.1038/nphys620) **, [****, ()](\doibase 10.1103/RevModPhys.82.1793) **, [****, ()](http://www.jetp.ac.ru/cgi-bin/dn/e{_}020{_}05{_}1307.pdf) **, [****, ()](\doibase
  10.1088/0022-3700/6/3/023) **, [****, ()](\doibase 10.1103/PhysRevA.22.770) **, [****, ()](\doibase 10.1103/PhysRevLett.71.1994) **, [****, ()](\doibase 10.1103/PhysRevA.49.2117) **, [****, ()](\doibase 10.1103/PhysRevLett.77.420) **, [****, ()](\doibase 10.1103/PhysRevLett.77.988) **, [****, ()](\doibase 10.1209/epl/i2000-00112-5) **, [****, ()](\doibase 10.1103/PhysRevLett.77.2360) **, [****, ()](\doibase 10.1103/PhysRevA.56.R3346) **, [****, ()](\doibase
  10.1103/PhysRevA.59.1477) **, [****, ()](\doibase
  10.1103/PhysRevA.55.1147) **, [****, ()](\doibase 10.1103/PhysRevA.56.555) **, [****, ()](http://dx.doi.org/10.1038/ncomms13449) **, [****, ()](\doibase
  10.1103/PhysRevLett.103.200401) **, [****, ()](http://dx.doi.org/10.1021/jp003731u) **, [****, ()](\doibase 10.1088/0953-4075/39/14/R01) **, [****, ()](\doibase 10.1038/35107000) **, [****, ()](\doibase 10.1103/PhysRevLett.88.173903) **, [****, ()](\doibase
  10.1103/RevModPhys.81.163) **, [****, ()](\doibase 10.1139/p70-183) in @noop [**]{},  (, ) pp.  **, [****, ()](\doibase
  10.1103/PhysRevA.27.2503) **, [****, ()](\doibase 10.1038/nature01414) **, [****, ()](\doibase 10.1103/PhysRevA.76.043401) **, [****, ()](\doibase 10.1088/1367-2630/10/2/025024) **, [****, ()](\doibase 10.1038/nature10196) **, [****, ()](\doibase 10.1126/science.1126259) **, [****, ()](\doibase
  10.1103/PhysRevLett.111.093002) **, [****, ()](\doibase 10.1103/PhysRevLett.115.113004) **, [****, ()](\doibase
  10.1103/PhysRevLett.65.2362) **, [****, ()](http://dx.doi.org/10.1038/nature11025) **, [****, ()](http://dx.doi.org/10.1038/nphys2125) **, [****, ()](http://dx.doi.org/10.1038/nature05648) **, [****, ()](\doibase 10.1126/science.aao7043) **, [****, ()](http://dx.doi.org/10.1126/science.1201351) **, [****, ()](\doibase 10.1088/0954-3899/40/5/053101) **, [****, ()](\doibase 10.1063/1.4773536) **, [****,  ()](\doibase http://dx.doi.org/10.1063/1.4907401)

[^1]: These authors contributed equally.

[^2]: These authors contributed equally.
