Gravitational caustics in an atom laser

Typically discussed in the context of optics, caustics are envelopes of classical trajectories (rays) where the density of states diverges, resulting in pronounced observable features such as bright points, curves, and extended networks of patterns. Here, we generate caustics in the matter waves of an atom laser, providing a striking experimental example of catastrophe theory applied to atom optics in an accelerated (gravitational) reference frame. We showcase caustics formed by individual attractive and repulsive potentials, and present an example of a network generated by multiple potentials. Exploiting internal atomic states, we demonstrate fluid-flow tracing as another tool of this flexible experimental platform. The effective gravity experienced by the atoms can be tuned with magnetic gradients, forming caustics analogous to those produced by gravitational lensing. From a more applied point of view, atom optics affords perspectives for metrology, atom interferometry, and nanofabrication. Caustics in this context may lead to quantum innovations as they are an inherently robust way of manipulating matter waves.

F rom light refracted by a sheet of glass to the light patterns seen on the ocean floor, rainbows, or the observation of gravitational lensing, caustics play a central role in the way optics presents itself in nature 1,2 . Unlike foci produced by optical instruments, caustics are generic in the sense that they do not need very specialized circumstances to exist, and are structurally stable 2 , leading to their widespread occurrence. They are formed, for example, when light is reflected (catacaustic) or refracted (diacaustic) from a curved surface [1][2][3][4] . While most visible in optics-prototypical caustics can readily be observed with polarized, coherent light 1,5-8 -the phenomenon of caustics and the underlying catastrophe theory 9,10 have found far-reaching interest. For example, caustics and catastrophe theory have been discussed in the context of generic two-mode quantum systems [11][12][13] , nuclear physics 14 , general relativity 15 , social sciences 16 , and robotics 17 . Caustics have also come into the focus of studies with electron microscopes where they may be exploited for advanced imaging techniques 18 .
Ultracold quantum gases provide a flexible platform for performing atom-optics experiments [19][20][21][22][23][24][25][26] where cold atoms, instead of photons, are used to generate atom-optical components with possible applications for fundamental science, atom interferometry, metrology, and new nano-fabrication approaches. Catastrophe atom optics-the formation of caustics by atom trajectories-has previously been discussed in specific settings, for example in the context of atoms being released from a magnetooptical trap 27,28 , atoms diffracting from a one-dimensional optical lattice 29,30 , or expanding Bose-Einstein condensates (BECs) with spatially varying initial phase 31 . To study caustics in matter waves, a particularly powerful tool and natural setting is an atom laser [32][33][34][35][36][37][38][39][40][41][42][43][44][45][46] , which is a coherent stream of collimated atoms outcoupled from a dilute-gas BEC. While caustics and catastrophe theory have been used to characterize the atom laser itself, in particular in terms of its transverse beam profile [47][48][49][50] , here we exploit the atom laser as a source of flow interacting with external potentials to generate a broad variety of caustic features. These features include individual fold and cusp caustics, and even complex caustic networks.
On large scales, the networks formed by caustics can be quite intricate. An example is shown in Fig. 1, where two repulsive optical Gaussian potentials are placed in the atom laser beam. Caustics arise from singularities in the continuous map (x i , t i ) ↦ (x, z) of atoms injected at time t = −t i and position (x = x i , z = 0) to the imaging plane (x, z) at the time of imaging t = 0. As shown in Fig. 1b, we can visualize this as a sheet embedded in three dimensions (x, z, t): the caustics occur where this sheet has vertical tangents, colored red in the figure. Despite the intricacy of the produced patterns, catastrophe theory reveals that all generic features can be categorized. In this case, we observe two stable types of singularities-folds and cusps.
To explore and classify the generated flow patterns, we apply a variety of flow-visualization techniques, including experimental fluid-flow tracing based on internal-state manipulation of the atoms, and simulated three-dimensional folded sheets. We achieve quantitative agreement between theory and experimental results. Atom optics differs from terrestrial light optics in a variety of ways, including the ability to easily introduce attractive as well as repulsive potentials, the power to perform internal state manipulation, e.g., for fluid flow tracing, and the ability to study non-negligible effects of gravity. The accelerated reference frame in our experiments, which is due to the presence of an effective gravity that can be adjusted using magnetic gradients, is natural for atom optics, and leads to features that would not exist otherwise, such as the transition from attached fan-like features to a detached crescent-shape caustic in the case of a repulsive potential. In contrast to previous work, the use of an atom laser combined with external potentials and an adjustable effective gravity provides a general platform for catastrophe atom optics investigations that can lead to rich observable dynamics with applications to cosmological effects. The direct imaging of sharply delineated features, a wealth of observable phenomena dependent on the shape, strength and sign of the potential(s), and the ability to trace atomic flow make this setting a very powerful platform.

Results
Experimental setup. The experimental setup for these investigations is schematically depicted in Fig. 2. A 87 Rb BEC is initially confined in an elongated harmonic trap. The trap is formed by a combination of an attractive focused dipole laser and an additional magnetic field gradient (see "Methods" section for details). In this setup, the x axis is oriented along the weakly confined direction of the trap, y along the imaging axis, and z vertically. The BEC is prepared in the F; m F ¼ 1; À1 j i spin state, which is supported against gravity by the trap. A 7 ms long microwave (mw) tone is used to coherently transfer atoms to the F; m F ¼ 2; À2 j i spin state, which is magnetically expelled out of the trap in the  3)). The right barrier is offset from this by (Δx, Δz) = (39, −42) μm (ε = 0.64). b A numerical rendering of the classical trajectories as a sheet ðx i ; t i Þ7 !ðx; z; t i À ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi À2z=a z p Þ, the singularities of which appear as caustics when projected down into the (x, z) imaging plane where the experimental data is shown again. Here, t is the vertical direction, while x increases to the left, and z decreases into the image. Red shading denotes regions where det J À1 is large [Eq. (1)], corresponding to the caustics when projected onto the imaging plane. c A numerical rendering of the classical trajectories as a sheet (x i , t i ) ↦ (x, z, t i ), which now includes the free-fall background.
downward direction with an acceleration of a z = 26.9(3) μm/ ms 2 ≈ 2.75g. This accelerated, collimated stream of atoms forms the atom laser. A laser propagating in the positive y direction generates an attractive or repulsive potential. Imaging is performed along the negative y direction, and a dichroic mirror is used to overlay the potential beam path with the imaging laser.
We note that our observation geometry differs from that typically used in optics where a wavefront passes through a refracting surface and then propagates to a screen on which it is observed. In our case, the images include the propagation direction as one axis of the two-dimensional observation plane.
Folds and cusps. To prepare for the discussion of our experimental observations, we begin with a brief theoretical review of caustics.
As we shall show below, the prominent features in our experiment can be well described by analyzing the classical dynamics of the atoms. In classical mechanics (optics), caustics occur on the envelope of classical trajectories q(t) where the trajectories (rays) and their tangents converge. These classical trajectories are stationary points of the action S 0 ½q ¼ 0, and the caustics arise when the hessian S″(q) is degenerate in some directions. These overlapping trajectories thus lead to an infinite density of states-classical divergences apparent to anyone who has accidentally started a fire from a curved mirror-that are softened by quantum mechanics where they are characterized by a breakdown in the Wentzel-Kramers-Brillouin (WKB) approximation 13,51,52 .
The atom laser continuously injects a thin line of atoms at height z i = 0, which then fall under a constant acceleration a z and are scattered by the potential. The injection region has limited spatial extent (~1 μm in the y and z directions), and we assume z i = y i ≈ 0 in our analysis.
Geometrically, this input can be described by a uniform sheet in the two-dimensional state space (x i , t i ) spanned by the initial injection sites of the atom laser x = x i , z = 0, at time t = − t i . The atoms then follow classical trajectories, falling under the acceleration a z , scattering from the optical potential(s), and ending at a final location (x, z) at the time of imaging, which we take as t = 0 so that t i has the interpretation of the time over which the atoms have fallen. Assuming a thin initial stream of atoms, the images represent a continuous mapping of ( Þ which we show in Fig. 1b. This representation removes the effect of the background acceleration a z from the visualization and accentuates the observed features.
Caustics correspond to singularities in this mapping. Assuming a constant injection rate-a good approximation for our experiments-the observed local density is inversely proportional to the determinant of the Jacobian of the mapping which becomes zero at the points corresponding to the caustics. These classical divergences are softened by quantum mechanics 51,52 , resulting in an Airy-function interference pattern.
The size of these features ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 4_ 2 =m 2 a z 3 q $ 0:5 μm is too small to observe in the current experiment, but ripe for future study.
Whitney 53 proved that continuous mappings from a plane into a plane can result in only two stable types of singularity-folds and cusps, both of which are observed here. Folds appear as curves where the surface in our three-dimensional embedding (x, z, t) has vertical tangents along t, and cusps appear as singular points where these folds meet. These are the only singularities that are stable, in the sense that they will persist, for example, even if the imaging angle is changed slightly. Additional nongeneric singularities can in principle be seen with this geometry, but these must be artificially tuned, for example by intentionally aligning cusp caustics. Arnold 10 provides a complete classification of these singularities.
We demonstrate here both fold and cusp caustics by inserting a Gaussian optical potential at position z = −h: where σ is the Gaussian waist and U 0 is the central strength of the potential. After scattering through this potential, classical particles will eventually move on parabolic trajectories. Qualitatively, the nature of the scattering and the associated caustics will be largely governed by the ratio The caustic structure changes dramatically at ε ≈ 1. In the experiments, the sign and magnitude of ε can be varied over a wide range, for example by adjusting the wavelength and intensity of the laser generating the optical potential.
Typical results are shown in Fig. 3. The top row ( Fig. 3a-d) was generated by inserting a repulsive potential (ε > 0) in the atom laser, while the bottom row ( Fig. 3e-h) was generated by inserting an attractive potential (ε < 0).
Repulsive potential. The results presented in Fig. 3a-d have been obtained using a repulsive Gaussian potential. The potentials were generated by a laser with a wavelength of 660 nm and a beam waist of σ ≈ 11.3(5) μm located h ≈ 78(3) μm below the trapped BEC. The ratio ε was varied by changing the laser intensity. For low values of ε, pronounced fan-like features are seen to emanate from both sides of the repulsive potential (Fig. 3a). The observed edge steepness of these features appears to be limited only by our imaging resolution of~3 μm. Our analysis presented below in the context of Fig. 4 identifies these edges as fold caustics. As the strength of the repulsive potential is increased, the attached fanlike features increase in width (Fig. 3b), ultimately detaching from the potential at ε ≈ 1 (Fig. 3c) and forming a half-ring-shaped detached caustic as ε increases above unity (Fig. 3d). This crossover occurs at a point where, at zero impact parameter, a falling classical atom would slow to rest at the top of the potential, signifying the reflection/transmission threshold.
The detached caustic here is a unique feature of atom optics in a sloped potential (such as the one generated by a constant downward acceleration) and would not exist in the absence of such a slope. Unlike the fan-like feature, the shape of this detached feature does not change as the potential height of the barrier is further increased. This is consistent with the analysis presented in the "Methods" section, which shows that any dependence should appear only weakly through effects related to the finite size of the potential. The observed features presented here are independent of the atomic density; the density of the atom laser is sufficiently low that mean-field effects can be neglected.
These features are reminiscent of the shockwaves created by a supersonic object moving through a fluid (see refs. 54,55 for classic examples). In particular, for ε < 1, the caustics look like an attached oblique shock (Fig. 3a, b), while for ε > 1, the shape of the pronounced caustic appearing above the potential (Fig. 3d) resembles that of a detached bow shock. The transition between these occurs for ε ≈ 1, and is shown in Fig. 3c, which shows faint signatures of both types of caustic. While the features in our experiment can be well-described by classical free-particle dynamics and, unlike shocks, do not involve non-linear selfsteepening, this analogy is intriguing. In our experiments, the atoms scatter from the potential with an impact velocity of about 6.5 cm/s. For comparison, even in the dense region of the trapped e-h An attractive potential is centered h = 46(1) μm below the trapped BEC. Individual panels are labeled by the energy ratio ϵ as defined in the main text. All images have been averaged over 6 independent experimental images with the same parameters. A faint grid is overlaid at 100 μm increments. The scale for square root of the density has been set so that pure white corresponds to the 99.95th percentile over all images to emphasize the caustic structure. The slight fringing seen in the images is due to optical effects in the imaging system, not matter-wave interference. to emphasize the scattering, overlaid on the experimental data. For clarity, a reversed color scheme has been used to plot the experimental images. e-h Corresponding map ðx i ; t i Þ 7 !ðx; z; t i À ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi À2z=a z p Þ with the vertical time axis scaled by a factor of 5 for frame (g). These numerics are generated for a single repulsive (a, b, e, f) or attractive (c, d, g, h) potential with the energy ratio noted at the top of each panel.
BEC, the bulk speed of sound c s ≈ 3 mm/s is more than an order of magnitude smaller. Thus, while our experiments operate in a regime where mean-field effects play no role, one can envision designing a similar procedure to explore supersonic, or even hypersonic, shockwaves 56 . It should be noted that our observations of caustics are distinct from the observation of Bogoliubov-Cherenkov radiation 57 .
To provide a clearer view on the physics behind the observed features, Fig. 4a, b show a comparison to numerical simulation for two repulsive potential strengths, along with corresponding visualizations in the form of folded sheets in Fig. 4e, f, respectively. As in Fig. 1b, the projection of the sheets onto the imaging plane reveal the caustics: Caustics occur along singularities of the map [Eq. (1)], which correspond to portions of this surface with vertical slopes. The repulsive potential causes the sheet to fold back on itself without intersections, forming multiples of four-fold caustics as one descends (Fig. 4e, f). Cusp caustics occur where the number of fold caustics changes.
For 0 < ε < 1, the sheet overlaps at most three times between the fan-like caustics (Fig. 4a-e). In the absence of an external acceleration, the maximally scattered trajectory always scatters less than 90°for Gaussian potentials. This implies that the outermost fan caustic is not an envelope, but is instead a single parabolic trajectory (see "Methods" section). In the transition region ε ≈ 1, the maximal scattering angle rapidly increases from 90°to 180°. For Gaussian potentials, this region exists only because of the finite acceleration a z ≠ 0, and for our parameters is limited to 0.997 ≲ ε < 1 (see "Methods" section).
As ε approaches unity, the abrupt change of the dynamics observed in the experiments corresponds to a drastic change in the sheet structure (Fig. 4f). In terms of the classical trajectories, for ε ≥ 1, the central particle with zero impact parameter x i = 0 bounces infinitely many times. This leads to an infinite series of overlapping sheets as shown in Fig. 4f. In the experimental image Fig. 3d, these collapse to a single feature that looks like a detached bow shock. Quantum mechanics softens this structure since particles can tunnel through the barrier, but on a scale that we cannot resolve in this experiment.
To provide a quantitative analysis, we measure the scattering angle of the caustic θ c (ε) which is dependent on the strength of the potential ε. As shown in the "Methods" section, the outer caustics of these fans correspond to a particular parabolic trajectory: where the maximum of the parabola (x Ã , z Ã ) depends on details of the scattering. In Fig. 5, we compare θ c (ε) extracted by fitting Eq. (4) to the outer caustics of the experimental images, with the values computed from the classical scattering problem. As discussed in the "Methods" section, the fact that |θ c | < 90°in this range is a peculiar feature of Gaussian potentials. The quantitative agreement seen between the experiment and classical scattering allows one to use caustics in the spirit of inverse classical scattering theory. Solving the inverse scattering problem provides sensitive input to calibrate properties of the potentials used in an experiment. For example, slight deviations from the theory, such as an asymmetry between left and right caustics, likely indicate slight asymmetries in the optical potentials. In this way, catastrophe atom optics can be used as a tool to precisely measure spatial properties of experimental potentials when designing atomtronic devices. For more information about the analysis of Fig. 5, see the "Methods" section.
Attractive potential. The great flexibility afforded by atom optical techniques allows one to not only change the strength of the potential but also its sign. In our experiment, we introduce an attractive potential by using a laser with a wavelength of 850 nm and a Gaussian beam waist of σ ≈ 27 μm located h ≈ 46 μm below the BEC. This attractive potential is an analog of a lens that focuses the atom laser. For very low potential depths (Fig. 3e), the focusing is weak and the potential appears to produce a nearly collimated, narrow stream of increased density below the potential. As the potential depth is increased, the parabolic trajectories from either side of the potential are seen to cross (Fig. 3f-h), and the structure is that of two-fold caustics emerging from a cusp (Fig. 4g, h). As with the repulsive potential, for maximal scattering less than 90°-which is the case for −2.59 ⪅ ε < 1 (see "Methods" section)-these caustics are also parabolic trajectories [Eq. (4)]. Figure 4c, d shows a comparison with classical trajectory calculations with corresponding sheet visualizations in Fig. 4g, h. Here, the attractive potential draws the sheet downward, causing it to eventually self-intersect. This allows for the formation of zero or two-fold caustics as one descends. A central cusp caustic occurs where the number of fold caustics changes from zero to two. As for the repulsive potentials, excellent quantitative agreement is found between numerics and the experimental results.
Fluid flow tracing. Going beyond the imaging of the overall fluid flow pattern, internal state manipulation of the atoms affords further powerful ways to visualize and analyze the flow: flow tracers can be created that indicate the evolution of a wavefront as it propagates along the atom laser stream. This technique is demonstrated in Fig. 6 where a horizontal line of atoms, located in the region between the trapped BEC and the potential, is transferred into a state that appears dark in the absorption images. The transfer is effected by a 100 μs mw pulse on the 2; À2 j i to 1; À1 j i transition. Spatial selectivity is possible due to the magnetic gradient in which the experiments are performed. The dark line flows with the atom laser, tracing slices of specific evolution times. In Fig. 6, a series of such lines has been injected into the stream at a fixed position between the trapped BEC and the potential in 0.5 ms time intervals. (4) to the outer left (orange squares) and right (green circles) caustics in the experimental images. The standard deviation at each ϵ provides an estimate of the uncertainty based on three data points. The solid line is the classical result θ c (ε) for a particle with kinetic energy K = mv 2 /2 = εU 0 scattering off of a Gaussian potential. The small differences between the left and right θ c above ε ≈ 0.5 is likely due to slight asymmetries in optical potential. Upon scattering from a strong repulsive potential (Fig. 6a), full dark rings are observed to propagate away from the potential and connect to the bow-shaped caustic, indicating the wavefront after the scattering event. Once the particles have fallen far enough, they are no longer influenced by the potential, and these bands are approximately circular arcs of radius t i ffiffiffiffiffiffiffiffiffi 2a z h p centered at height Àh þ ða z À gÞt 2 I =2 À a z t 2 i =2, extending over θ ∈ [−θ c , θ c ] (see "Methods" section). This provides a very visual explanation for the emergence of the caustic as an envelope of rays.
Fluid flow tracing in the case of an attractive potential is shown in Fig. 6b. The image clearly reveals how the initially horizontal fluid tracer lines are drawn into the region of the potential, from which they emerge as loop structures. The left and right apexes of these loops reveal the position of the fold caustics, providing an independent and direct experimental visualization of the formation of the caustics. Red dotted lines in Fig. 6 show the results of classical trajectories, which are in agreement with the experiment.

Discussion
In contrast to catastrophe optics of light, catastrophe atom optics has unusual features and provides unique opportunities. For example, the pronounced half-ring-shaped caustic appearing for sufficiently strong repulsive potentials would not exist without the influence of gravity, and a similar feature would be difficult to see with light outside of cosmological contexts. The experiments presented here demonstrate the ability to generate and directly image complex caustics in an accelerated reference frame. These experiments build upon and contribute to the field of catastrophe atom optics by introducing a powerful experimental setting for investigations.
This system suggests an intuitive approach for visualizing the origin of the generic cusp and fold caustics proved by Whitney 53 . Real-time dynamics provide a natural embedding of the mapping whose singularities define the caustics, and our technique of fluidflow tracing generates direct experimental images of slices through this embedding, aiding the interpretation of such atom optics experiments.
The complexity of the observed phenomenology suggests many interesting future extensions of this work, including the construction of more complex caustic networks and the study of departures from classical catastrophe theory as inter-atomic interactions and quantum interference effects start to appear. Furthermore, our demonstration of fluid-flow tracing might find interesting applications in the study of quantum turbulence or quantum shocks, where following fluid tracers can help to reveal underlying flow patterns.
Our theoretical investigations complement these experiments by using a range of techniques. The results presented here are based on classical trajectories, providing a concrete demonstration of the principles behind classical catastrophe theory 9,10 in a context beyond the usual setting in optics. Extending these to include quantum effects can proceed through a semi-classical expansion, with leading order corrections given by the WKB approximation 52,58-60 , and full quantum effects included numerically (see, e.g., refs. 13,61-63 ). As a future direction, the excellent agreement with the experiment allows us to validate these different approaches, and thus to design future experiments that will probe higher-order corrections. This work also lays the foundation to further explore the quantum ramifications of caustics: Appearing as singularities in the classical action 13,51,52 , caustics play a similar role to quantum scars 64 , giving rise to robust features in the presence of classical chaos. Atom lasers have a distinct advantage over optics here in that non-linear interactions can also be manipulated.
From a point of view of applications, it is worth noting that caustics also occur in electron microscopes 18 . One can imagine an analogous setup for an atom laser, where a specimen under study is placed into the atom laser stream by an optical tweezer operated at the tune-out wavelength 65 for Rb so that the trap itself does not perturb the atom laser. In this context, the experimental platform presented here provides a flexible model system for studying the effects of advanced beam shaping techniques.

Methods
Experimental setup. To investigate the dynamics of an atom laser scattered by a Gaussian potential, a 87 Rb Bose-Einstein condensate composed of 5.5 × 10 5 atoms is initially prepared in the F; m F ¼ 1; À1 j ispin state. The condensate is confined in a hybrid trap consisting of a single dipole trap with a waist of 20 μm and a magnetic quadrupole field which has been vertically shifted above the dipole trap. This leads to a hybrid trap with trap frequencies of {ω x , ω y , ω z } = 2π × {7.1, 167, 180} Hz. The initial condensate in the F; m F ¼ 1; À1 j ispin state has a Thomas-Fermi (TF) radius in the x direction of 93 μm.
Atoms are then ejected from the trap by resonantly exciting atoms from the 1; À1 j ispin state to the 2; À2 j istate using a 7 ms mw pulse. Atoms in the 2; À2 j i state accelerate downwards at a z = 26.9(3) μm/ms 2 in the negative z direction.
The repulsive (attractive) Gaussian potential is created using a 660 nm (850 nm) laser focused to a σ = 11.3(5) μm (27(1) μm) Gaussian waist located h = 78(3) μm (46(3) μm) below the trapped BEC. In the case of Fig. 1a, where two repulsive potentials were used, orthogonal polarizations of light were used to prevent interference effects between the beams. Absorption imaging is performed along the −y direction, using the F ¼ 2 ! 3 0 cycling transition after a brief 0.5 ms timeof-flight expansion.
Classical trajectories. Here we present the analysis of the classical trajectories scattered by the Gaussian potential Eq. (2). Sufficiently far from the potential, U(x, z) ≈ 0 and the trajectories will be parabolic due to the constant downward acceleration a z : qðtÞ ¼ xðtÞ Here is the speed of particles falling from the injection site without the potential. In the limit of a zero-range potential, σ → 0, one has x 0 (x i ) → 0, z 0 (x i ) → 0. In this limit, θ(x i ) is the scattering angle at time t = 0 when the particle hits the potential, and is the only parameter that depends on the impact parameter Fig. 6 Fluid-flow tracing in the presence of a repulsive or attractive potential. a A repulsive potential with waist σ ≈ 12 μm is located h = 68.0(5) μm below the trapped BEC. A small horizontal stripe of atoms (located in the region between the trapped BEC and the potential) is converted back to the 1; À1 j istate by a series of brief 100 μs mw pulses in 0.5 ms intervals, leaving a void of atoms in the imaged 2; À2 j istate. Overlaying the experimental image are the corresponding classical trajectories at these fixed t i , shown as dotted (red) curves for positive initial impact parameters 0 < x 0 < 100 μm. These correspond to slices through the embedded sheet (x i , t i ) ↦ (x, z, t i ) (see for example Fig. 1c), demonstrating that these fluid-flow tracers provide a direct way to visualize the folding of this sheet which gives rise to the caustics. b Same as a, but for an attractive potential with waist σ ≈ 25 μm located h = 50(5) μm below the trapped BEC. Both of the experimental images here have been averaged over 101 runs.
x i . Deviations from these values characterize the finite-size effects of the potential, and must be calculated numerically.
The effects of imaging after a time-of-flight expansion of t I = 0.5 ms can be included by extending the trajectories from time t to time t f = t + t I under the reduced acceleration g: The effect of imaging is, thus, to shift the trajectories after scattering vertically. The classical trajectories can be found by eliminating t f , and are inverted parabola with maxima at (x Ã , z Ã ): Note that if the scattering is such that the maximum scattering angle is less than 90°(|θ c | < π/2), then the trajectory where θ = θ c will have the widest parabola. This trajectory will then eventually overtake all other trajectories and corresponds to the limit of trajectories as θ → θ c , and will lie along a caustic, corresponding to a singularity in det J ¼ 0 [Eq. (1)] with vanishing partials ∂/∂x i . To perform the analysis in Fig. 5 of the case where 0 < ϵ < 1, we locate the outer caustic in each image and then fit the parabolic scattering trajectory to these, extracting the angle θ c (ε). This analysis was repeated for three experimental runs of the same parameters, where the resulting average and estimated error are reported in Fig. 5.
All data has been corrected for a small camera tilt of θ cam = 0.81°, which was measured using a falling BEC from a tightly confined trap.
Neglecting the effects of acceleration during the scattering, θ c (ε) depends only on the dimensionless ratio ε. In particular, θ c is independent of the waist σ of the potential, which only changes which impact parameter x i (ε, σ) scatters maximally. Interestingly, as shown in Fig. 5, scattering from a Gaussian potential is somewhat peculiar in that the maximum scattering angle for −2.59 ⪅ ε ⪅ 1 is less than 90°. This is the case for all examples considered here except for Fig. 3d which has ϵ > 1 and hence θ c = 0 for the bouncing trajectories. In this limit, the transition from ε < 1 to ε > 1 is discontinuous with the sudden disappearance of the attached oblique-shock-like caustic and the appearance of the detached bow-shock-like caustic.
This discontinuity is slightly softened by the fact that, during the scattering, the acceleration a z is still present. In this case, θ c (ε, ς) depends on the two dimensionless ratios The effective potential for the central trajectory x i = 0 is thus: VðzÞ mah ¼ṼðzÞ ¼z þ εe À2z 2 =ς 2 ;z ¼ z þ h h : This particle will bounce ifṼ max ðzÞ ≥Ṽð1Þ % 1, where the latter approximation for all values of ς we consider for the repulsive potentials. The explicit solution to this optimization problem is: With our parameters, ς = 0.14, corresponding to a transition region of 0.997 ≤ε < 1. Frame Fig. 3c, for example, shows faint signatures of both features, and thus sits right at this transition.
We can also consider the slices of constant t f corresponding to the dark bands in Fig. 6: In the zero-range limit x 0 , z 0 → 0 no longer depend on x i so we can solve for θ(x i ) to obtain: These are circular arcs of radius v 0 t f centered at height Àh þ ða z À gÞt 2 I =2À a z t 2 f =2, extending through θ ∈ [−θ c , θ c ]. Our numerical results do not make these approximations. The classical trajectories are found by integrating the classical equations of motion with the physical potential, including the expansion time. The agreement between these, however, allows us to assert that the quantitative corrections from the acceleration during scattering are small (sub percent).

Data availability
The data used to generate the figures in the manuscript and accompanying code are available under osf.io/kdm9s/ 66 . Additional experimental and numerical data sets in this work will be made available from the corresponding authors upon reasonable request.