Search

Instabilities, where small fluctuations seed the formation of large-scale structures, govern dynamics in a variety of fluid systems. The Rayleigh-Taylor instability (RTI), present from tabletop to astronomical scales, is an iconic example characterized by mushroom-shaped incursions appearing when immiscible fluids are forced together. Despite its ubiquity, RTI experiments are challenging; here, we report the observation of the RTI in an immiscible binary superfluid consisting of a two-component Bose-Einstein condensate. We force these components together to initiate the instability, and observe the growth of mushroom-like structures. The interface can also be stabilized, allowing us to spectroscopically measure the "ripplon" interface modes. Lastly, we use matter-wave interferometry to transform the superfluid velocity field at the interface into a vortex chain. These results-in agreement with our theory-demonstrate the close connection between the RTI in classical and quantum fluids., Comment: 10 pages (including SM), 4 figures + 3 figures in SM

We study the problem of detecting or recovering a planted ranked subgraph from a directed graph, an analog for directed graphs of the well-studied planted dense subgraph model. We suppose that, among a set of $n$ items, there is a subset $S$ of $k$ items having a latent ranking in the form of a permutation $\pi$ of $S$, and that we observe a fraction $p$ of pairwise orderings between elements of $\{1, \dots, n\}$ which agree with $\pi$ with probability $\frac{1}{2} + q$ between elements of $S$ and otherwise are uniformly random. Unlike in the planted dense subgraph and planted clique problems where the community $S$ is distinguished by its unusual density of edges, here the community is only distinguished by the unusual consistency of its pairwise orderings. We establish computational and statistical thresholds for both detecting and recovering such a ranked community. In the log-density setting where $k$, $p$, and $q$ all scale as powers of $n$, we establish the exact thresholds in the associated exponents at which detection and recovery become statistically and computationally feasible. These regimes include a rich variety of behaviors, exhibiting both statistical-computational and detection-recovery gaps. We also give finer-grained results for two extreme cases: (1) $p = 1$, $k = n$, and $q$ small, where a full tournament is observed that is weakly correlated with a global ranking, and (2) $p = 1$, $q = \frac{1}{2}$, and $k$ small, where a small "ordered clique" (totally ordered directed subgraph) is planted in a random tournament., Comment: 79 pages, 1 figure. This paper subsumes most of the results of the earlier arXiv:2407.16597 by a subset of the authors

Weak measurement enables the extraction of targeted information from a quantum system while minimizing decoherence due to measurement backaction. However, in many-body quantum systems backaction can have unexpected effects on wavefunction collapse. We theoretically study a minimal many-particle model consisting of weakly measured non-interacting fermions in a one dimensional lattice. Repeated measurement of on-site occupation number with single-site resolution stochastically drives the system toward a Fock state, regardless of the initial state. This need not be the case for measurements that do not, even in principle, have single-site spatial resolution. We numerically show for systems with up to 16 sites that decreasing the spatial resolution strongly affects both the rate of stochastic evolution for each quantum trajectory and the allowed final states. The full Hilbert space can be partitioned into backaction-free subspaces (BFSs) the elements of which are indistinguishable to these measurements. Repeated measurements will drive any initial state into a single BFS, leading to a steady state that is a fixed point of the measurement process. We exactly calculate the properties of these BFSs for systems up to 32 sites and find that even for moderate reductions in measurement resolution they yield non-trivial steady state entanglement and coherence.

The Pound-Drever-Hall (PDH) technique is routinely used to stabilize the frequency of a laser to a reference cavity. The electronic sideband (ESB) locking scheme, a PDH variant, helps bridge the frequency difference between the quantized frequencies enforced by the cavity and the laser frequency of interest. Here we use quadrature amplitude modulation (QAM), a technique used in digital signal communication, to engineer the high-quality phase-modulated radio-frequency (rf) signal required for ESB locking scheme. We develop a theoretical framework to analyze the effects of in-phase/quadrature-phase (I/Q) impairments on the ESB error signal for ultra-narrow linewidth lasers. We design and implement two baseband-sampling software-defined radio variants for implementing QAM that compensate for these I/Q impairments. Using these variants, we engineer high-quality phase-modulated radio-frequency (rf) signals with a large phase modulation index of 1.01 radians, a maximum modulation frequency of 3 MHz, a tunable carrier wave frequency range of 450 MHz to 4 GHz, and I/Q errors of less than 2.25 % over the entire carrier wave frequency range., Comment: 11 pages, 7 figures

Much of our knowledge of quantum systems is encapsulated in the expectation value of Hermitian operators, experimentally obtained by averaging projective measurements. However, dynamical properties are often described by products of operators evaluated at different times; such observables cannot be measured by individual projective measurements, which occur at a single time. For example, the dynamical structure factor describes the propagation of density excitations, such as phonons, and is derived from the spatial density operator evaluated at different times. Conventionally, this is measured by first exciting the system at a specific wavevector and frequency, then measuring the response. Here, we describe an alternative approach using a pair of time-separated weak measurements, and analytically show that their cross-correlation function directly recovers the dynamical structure factor. We provide numerical confirmation of this technique with a matrix product states simulation of the one-dimensional Bose-Hubbard model, weakly measured by phase contrast imaging. We explore the limits of the method and demonstrate its applicability to real experiments with limited imaging resolution., Comment: 6 pages, 4 figures

Geometrical frustration and long-range couplings are key contributors to create quantum phases with different properties throughout physics. We propose a scheme where both ingredients naturally emerge in a Raman induced subwavelength lattice. We first demonstrate that Raman-coupled multicomponent quantum gases can realize a highly versatile frustrated Hubbard Hamiltonian with long-range interactions. The deeply subwavelength lattice period leads to strong long-range interparticle repulsion with tunable range and decay. We numerically demonstrate that the combination of frustration and long-range couplings generates many-body phases of bosons, including a range of density-wave and superfluid phases with broken translational and time reversal symmetries, respectively. Our results thus represent a powerful approach for efficiently combining long-range interactions and frustration in quantum simulations., Comment: 17 pages, 10 figures

We investigated turbulence in atomic Bose-Einstein condensates (BECs) using a minimally destructive, impurity injection technique analogous to particle image velocimetry in conventional fluids. Our approach transfers small regions of the BEC into a different hyperfine state, and tracks their displacement ultimately yielding the velocity field. This allows us to quantify turbulence in the same way as conventional in fluid dynamics in terms of velocity-velocity correlation functions called velocity structure functions that obey a Kolmogorov scaling law. Furthermore the velocity increments show a clear fat-tail non-Gaussian distribution that results from intermittency corrections to the initial "K41" Kolmogorov theory. Our observations are fully consistent with the later "KO62" description. These results are validated by a 2D dissipative Gross-Pitaevskii simulation., Comment: 6 pages, 4 figures

A static electric field of a few V/cm shifts the energy levels of ultracold Rydberg atoms in a magneto-optical trap. For a given principle quantum number, most of the energy levels are nearly degenerate at zero field and fan out with increasing field to form a manifold. We excite Rydberg atoms to energy levels near the center of the manifold, where the spacing is nearly harmonic, and allow them to exchange energy via resonant dipole-dipole interactions. We measure the time evolution as energy spreads away from the center of the manifold, which reveals that the system fails to thermalize for long interaction times. A computational model that includes only a few essential features of the system qualitatively agrees with this result., Comment: 6 pages, 4 figures

We consider the problem of inferring an unknown ranking of $n$ items from a random tournament on $n$ vertices whose edge directions are correlated with the ranking. We establish, in terms of the strength of these correlations, the computational and statistical thresholds for detection (deciding whether an observed tournament is purely random or drawn correlated with a hidden ranking) and recovery (estimating the hidden ranking with small error in Spearman's footrule or Kendall's tau metric on permutations). Notably, we find that this problem provides a new instance of a detection-recovery gap: solving the detection problem requires much weaker correlations than solving the recovery problem. In establishing these thresholds, we also identify simple algorithms for detection (thresholding a degree 2 polynomial) and recovery (outputting a ranking by the number of "wins" of a tournament vertex, i.e., the out-degree) that achieve optimal performance up to constants in the correlation strength. For detection, we find that the above low-degree polynomial algorithm is superior to a natural spectral algorithm. We also find that, whenever it is possible to achieve strong recovery (i.e., to estimate with vanishing error in the above metrics) of the hidden ranking, then the above "Ranking By Wins" algorithm not only does so, but also outputs a close approximation of the maximum likelihood estimator, a task that is NP-hard in the worst case., Comment: 34 pages

The radiative decay of excited states of the negative ion of rhodium, Rh$^-$, has been investigated experimentally and theoretically. The experiments were conducted at the Double ElectroStatic Ion Ring Experiment (DESIREE) facility at Stockholm University using selective photodetachment from a stored ion beam to monitor the time evolution of the excited state populations. The lifetimes of the Rh$^-$ $^3F_{3}$ and $^3F_{2}$ fine structure levels were measured to be 3.2(6)~s and 21(4)~s, respectively. An additional, previously unreported, higher-lying bound state of mixed $^1D_2+^3P_2+(4d^95s)^1D_2+^3F_2$ composition was observed and found to have a lifetime of 10.9(8)s. The binding energy of this state was determined to be in the interval $0.1584(2) $ eV $ < E_b < 0.2669(2)$ eV, using laser photodetachment threshold (LPT) spectroscopy. An autodetaching state with a lifetime of 480(10) microseconds was also observed. Theoretical calculations of the excited-state compositions, energies, and magnetic-dipole transition lifetimes were performed using the multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction methods. The calculated lifetimes of the $^3F_{3}$ and $^3F_{2}$ fine structure levels are in excellent agreement with the measured values. The present study should provide valuable insights into electron correlation effects in negative ions and forbidden radiative transitions.

The purpose of this style guide is to assist developers in designing effective and consistent-looking user interfaces for accelerator control rooms. A similar purpose is to help developers avoid the creation of user interfaces that needlessly stray from the accepted standard set forth in this document. This way, all interfaces combined will look congruous. This is especially beneficial for development that spans multiple years by many different contributors. This document is intended as a ready reference source for all user interface design for the Fermilab accelerator complex.

Relevance labels, which indicate whether a search result is valuable to a searcher, are key to evaluating and optimising search systems. The best way to capture the true preferences of users is to ask them for their careful feedback on which results would be useful, but this approach does not scale to produce a large number of labels. Getting relevance labels at scale is usually done with third-party labellers, who judge on behalf of the user, but there is a risk of low-quality data if the labeller doesn't understand user needs. To improve quality, one standard approach is to study real users through interviews, user studies and direct feedback, find areas where labels are systematically disagreeing with users, then educate labellers about user needs through judging guidelines, training and monitoring. This paper introduces an alternate approach for improving label quality. It takes careful feedback from real users, which by definition is the highest-quality first-party gold data that can be derived, and develops an large language model prompt that agrees with that data. We present ideas and observations from deploying language models for large-scale relevance labelling at Bing, and illustrate with data from TREC. We have found large language models can be effective, with accuracy as good as human labellers and similar capability to pick the hardest queries, best runs, and best groups. Systematic changes to the prompts make a difference in accuracy, but so too do simple paraphrases. To measure agreement with real searchers needs high-quality "gold" labels, but with these we find that models produce better labels than third-party workers, for a fraction of the cost, and these labels let us train notably better rankers.

A majority of ultracold atom experiments utilize resonant absorption imaging techniques to obtain the atomic density. To make well-controlled quantitative measurements, the optical intensity of the probe beam must be precisely calibrated in units of the atomic saturation intensity ${I_{\rm sat}}$. In quantum gas experiments, the atomic sample is enclosed in an ultra-high vacuum system that introduces loss and limits optical access; this precludes a direct determination of the intensity. Here, we use quantum coherence to create a robust technique for measuring the probe beam intensity in units of ${I_{\rm sat}}$ via Ramsey interferometry. Our technique characterizes the ac Stark shift of the atomic levels due to an off-resonant probe beam. Furthermore, this technique gives access to the spatial variation of the probe intensity at the location of the atomic cloud. By directly measuring the probe intensity just before the imaging sensor our method in addition yields a direct calibration of imaging system losses as well as the quantum efficiency of the sensor., Comment: 15 pages, 7 figures

We describe a two-dimensional optical lattice for ultracold atoms with spatial features below the diffraction limit created by a bichromatic optical standing wave. At every point in space these fields couple the internal atomic states in a three-level Lambda coupling configuration. Adiabatically following the local wavefunction of the resulting dark state yields a spatially uniform Born-Oppenheimer potential augmented by geometric scalar and vector potentials appearing due to spatially rapid changes of the wavefunction. Depending on system parameters, we find that the geometric scalar potential can interpolate from a 2D analogue of the Kronig-Penney lattice, to an array of tubes with a zig-zag shaped barrier. The geometric vector potential induces a spatially periodic effective magnetic field (the Berry's curvature) that can be tuned to cause destructive interference between neighboring tubes, thereby decoupling them at a critical point in parameter space. We numerically investigate the energy spectrum including decay from the excited state, and find that the adiabatic approximation is sound for strong coupling strengths, leading to negligible loss in the dark state manifold. Furthermore, the spectrum is well-described by a non-Hermitian tight binding model with on-site losses, and hopping characterized by both loss and, surprisingly, gain., Comment: 11 pages, 7 figures. For supplementary numerical codes, see http://web.vu.lt/ff/g.juzeliunas/2D_Lambda_Supplementary.zip

We introduce a new algorithm and software for solving linear equations in symmetric diagonally dominant matrices with non-positive off-diagonal entries (SDDM matrices), including Laplacian matrices. We use pre-conditioned conjugate gradient (PCG) to solve the system of linear equations. Our preconditioner is a variant of the Approximate Cholesky factorization of Kyng and Sachdeva (FOCS 2016). Our factorization approach is simple: we eliminate matrix rows/columns one at a time and update the remaining matrix using sampling to approximate the outcome of complete Cholesky factorization. Unlike earlier approaches, our sampling always maintains a connectivity in the remaining non-zero structure. Our algorithm comes with a tuning parameter that upper bounds the number of samples made per original entry. We implement our algorithm in Julia, providing two versions, AC and AC2, that respectively use 1 and 2 samples per original entry. We compare their single-threaded performance to that of current state-of-the-art solvers Combinatorial Multigrid (CMG), BoomerAMG-preconditioned Krylov solvers from HyPre and PETSc, Lean Algebraic Multigrid (LAMG), and MATLAB's with Incomplete Cholesky Factorization (ICC). Our evaluation uses a broad class of problems, including all large SDDM matrices from the SuiteSparse collection and diverse programmatically generated instances. Our experiments suggest that our algorithm attains a level of robustness and reliability not seen before in SDDM solvers, while retaining good performance across all instances. Our code and data are public, and we provide a tutorial on how to replicate our tests. We hope that others will adopt this suite of tests as a benchmark, which we refer to as SDDM2023. Our solver code is available at: https://github.com/danspielman/Laplacians.jl/ Our benchmarking data and tutorial are available at: https://rjkyng.github.io/SDDM2023/

We experimentally and theoretically investigate the anisotropic speed of sound of an atomic superfluid (SF) Bose-Einstein condensate in a 1D optical lattice. Because the speed of sound derives from the SF density, this implies that the SF density is itself anisotropic. We find that the speed of sound is decreased by the optical lattice, and the SF density is concomitantly reduced. This reduction is accompanied by the appearance of a normal fluid in the purely Bose condensed phase. The reduction in SF density -- first predicted [A. J. Leggett, Phys. Rev. Lett. 1543--1546 (1970)] in the context of supersolidity -- results from the coexistence of superfluidity and density modulations, but is agnostic about the origin of the modulations. We additionally measure the moment of inertia of the system in a scissors mode experiment, demonstrating the existence of rotational flow. As such we shed light on some supersolid properties using imposed, rather than spontaneously formed, density-order.

We consider solitary wave excitations above the ground state of $F=1$ spin-orbit coupled Bose-Einstein condensates (SOBECs). The low energy properties of SOBECs in any of the three branches of the single particle dispersion relation can be described by suitable scalar nonlinear Schr\"odinger (NLS) equations which we obtain using multiple-scale expansions. This enables us to examine a variety of different configurations, such as dark solitary waves associated with higher energy branches, as well as dark and bright structures in the lowest branch. The lowest branch can also exhibit a ``superstripe'' phase that supports solitary waves. In all cases, we provide explicit expressions for the NLS coefficients, and confirm their validity with full numerical simulations of the SOBEC system including a harmonic confining potential., Comment: 16 figures

Ultracold atoms are an ideal platform for understanding system-reservoir dynamics of many-body systems. Here, we study quantum back-action in atomic Bose-Einstein condensates, weakly interacting with a far-from resonant, i.e., dispersively interacting, probe laser beam. The light scattered by the atoms can be considered as a part of quantum measurement process whereby the change in the system state derives from measurement back-action. We experimentally quantify the resulting back-action in terms of the deposited energy. We model the interaction of the system and environment with a generalized measurement process, leading to a Markovian reservoir. Further, we identify two systematic sources of heating and loss: a stray optical lattice and probe-induced light assisted collisions (an intrinsic atomic process). The observed heating and loss rates are larger for blue detuning than for red detuning, where they are oscillatory functions of detuning with increased loss at molecular resonances and reduced loss between molecular resonances., Comment: 13 pages, 8 figures. arXiv admin note: text overlap with arXiv:2209.04400

Recent experiments demonstrated deeply subwavelength lattices using atoms with $N$ internal states Raman-coupled with lasers of wavelength $\lambda$. The resulting unit cell was $\lambda/2N$ in extent, an $N$-fold reduction compared to the usual $\lambda/2$ periodicity of an optical lattice. For resonant Raman coupling, this lattice consists of $N$ independent sinusoidal potentials (with period $\lambda/2$) displaced by $\lambda/2N$ from each other. We show that detuning from Raman resonance induces tunneling between these potentials. Periodically modulating the detuning couples the $s$- and $p$-bands of the potentials, creating a pair of coupled subwavelength Rice--Mele chains. This operates as a novel topological charge pump that counter-intuitively can give half the displacement per pump cycle of each individual Rice--Mele chain separately. We analytically describe this behavior in terms of infinite-system Chern numbers, and numerically identify the associated finite-system edge states., Comment: 12 pages, 8 figures

A fundamental tenet of quantum mechanics is that measurements change a system's wavefunction to that most consistent with the measurement outcome, even if no observer is present. Weak measurements produce only limited information about the system, and as a result only minimally change the system's state. Here, we theoretically and experimentally characterize quantum back-action in atomic Bose-Einstein condensates interacting with a far-from resonant laser beam. We theoretically describe this process using a quantum trajectories approach where the environment measures the scattered light and present a measurement model based on an ideal photodetection mechanism. We experimentally quantify the resulting wavefunction change in terms of the contrast of a Ramsey interferometer and control parasitic effects associated with the measurement process. The observed back-action is in good agreement with our measurement model; this result is a necessary precursor for achieving true quantum back-action limited measurements of quantum gases., Comment: 15 pages including the SI; 5 figures plus 2 figures in the SI

Here we revisit the topic of stationary and propagating solitonic excitations in self-repulsive three-dimensional Bose-Einstein condensates by quantitatively comparing theoretical analysis and associated numerical computations with our experimental results. Using fully 3d numerical simulations, we explore the existence, stability, and evolution dynamics of planar dark solitons, as well as their instability-induced decay products including solitonic vortices and vortex rings. In the trapped case and with no adjustable parameters, our numerical findings are in correspondence with experimentally observed coherent structures. Without a longitudinal trap, we identify numerically exact traveling solutions and quantify how their transverse destabilization threshold changes as a function of the solitary wave speed.

We simulate the dynamics of Rydberg atoms resonantly exchanging energy via two-, three-, and four-body dipole-dipole interactions in a one-dimensional array. Using simplified models of a realistic experimental system, we study the initial state survival probability, mean level spacing, spread of entanglement, and properties of the energy eigenstates. By exploring a range of disorders and interaction strengths, we find regions in parameter space where the three- and four-body dynamics either fail to thermalize or do so slowly. The interplay between the stronger hopping and weaker field-tuned interactions gives rise to quantum many-body scar states, which play a critical role in slowing the dynamics of the three- and four-body interactions., Comment: 21 pages, 19 figures

Continuously measured interacting quantum systems almost invariably heat, causing loss of quantum coherence. Here, we study Bose-Einstein condensates (BECs) subject to repeated weak measurement of the atomic density and describe several protocols for generating a feedback signal designed to remove excitations created by measurement backaction. We use a stochastic Gross-Pitaevskii equation to model the system dynamics and find that a feedback protocol utilizing momentum dependant gain and filtering can effectively cool both 1D and 2D systems. The performance of these protocols is quantified in terms of the steady state energy, entropy, and condensed fraction. These are the first feedback cooling protocols demonstrated in 2D, and in 1D our optimal protocol reduces the equilibrium energy by more than a factor of 100 as compared with a previous cooling protocol developed using the same methodology. We also use this protocol to quench-cool 1D BECs from non-condensed highly excited states and find that they rapidly condense into a far from equilibrium state with energy orders of magnitude higher than the equilibrium ground state energy for that condensate fraction. We explain this in terms of the near-integrability of our 1D system, whereby efficiently cooled low momentum modes are effectively decoupled from the energetic `reservoir' of the higher momentum modes. We observe that the quench-cooled condensed states can have non-zero integer winding numbers described by quantized supercurrents., Comment: 14 pages, 7 figures

We establish a dataset of over $1.6\times10^4$ experimental images of Bose--Einstein condensates containing solitonic excitations to enable machine learning (ML) for many-body physics research. About $33~\%$ of this dataset has manually assigned and carefully curated labels. The remainder is automatically labeled using SolDet -- an implementation of a physics-informed ML data analysis framework -- consisting of a convolutional-neural-network-based classifier and OD as well as a statistically motivated physics-informed classifier and a quality metric. This technical note constitutes the definitive reference of the dataset, providing an opportunity for the data science community to develop more sophisticated analysis tools, to further understand nonlinear many-body physics, and even advance cold atom experiments., Comment: 16 pages, 4 figures

Marcus, Spielman and Srivastava (Annals of Mathematics 2014) solved the Kadison--Singer Problem by proving a strong form of Weaver's conjecture: they showed that for all $\alpha > 0$ and all lists of vectors of norm at most $\sqrt{\alpha}$ whose outer products sum to the identity, there exists a signed sum of those outer products with operator norm at most $\sqrt{8 \alpha} + 2 \alpha.$ We prove that it is NP-hard to distinguish such a list of vectors for which there is a signed sum that equals the zero matrix from those in which every signed sum has operator norm at least $\kappa \sqrt{\alpha}$, for some absolute constant $\kappa > 0.$ Thus, it is NP-hard to construct a signing that is a constant factor better than that guaranteed to exist. For $\alpha = 1/4$, we prove that it is NP-hard to distinguish whether there is a signed sum that equals the zero matrix from the case in which every signed sum has operator norm at least $1/4$.

Nontrivial topology in lattices is characterized by invariants--such as the Zak phase for one dimensional (1D) lattices--derived from wave functions covering the Brillouin zone. We realized the 1D bipartite Rice-Mele (RM) lattice using ultracold $^{87}$Rb and focus on lattice configurations possessing various combinations of chiral, time-reversal and particle-hole symmetries. We quenched between configurations and used a form of quantum state tomography, enabled by diabatically tuning lattice parameters, to directly follow the time evolution of the Zak phase as well as a chiral winding number. The Zak phase evolves continuously; however, when chiral symmetry transiently appears in the out-of-equilibrium system, the chiral winding number is well defined and can take on different integer values. When quenching between two configurations obeying all three symmetries the Zak phase is time independent; we confirm the contrasting prediction of [M. McGinley and N. R.Cooper, PRL 121 090401 (2018)] that chiral symmetry is periodically restored, at which times the winding number changes by $\pm 2$, yielding values that are not present in the native RM Hamiltonian.

We experimentally realized a time-periodically modulated 1D lattice for ultracold atoms featuring a pair of linear bands, each associated with a Floquet winding number: a topological invariant. These bands are spin-momentum locked and almost perfectly linear everywhere in the Brillouin zone (BZ), making this system a near-ideal realization of the 1D Dirac Hamiltonian. We characterized the Floquet winding number using a form of quantum state tomography, covering the BZ and following the micromotion through one Floquet period. Lastly, we altered the modulation timing to lift the topological protection, opening a gap at the Dirac point that grew in proportion to the deviation from the topological configuration.

In ultracold-atom experiments, data often comes in the form of images which suffer information loss inherent in the techniques used to prepare and measure the system. This is particularly problematic when the processes of interest are complicated, such as interactions among excitations in Bose-Einstein condensates (BECs). In this paper, we describe a framework combining machine learning (ML) models with physics-based traditional analyses to identify and track multiple solitonic excitations in images of BECs. We use an ML-based object detector to locate the solitonic excitations and develop a physics-informed classifier to sort solitonic excitations into physically motivated subcategories. Lastly, we introduce a quality metric quantifying the likelihood that a specific feature is a longitudinal soliton. Our trained implementation of this framework, SolDet, is publicly available as an open-source python package. SolDet is broadly applicable to feature identification in cold-atom images when trained on a suitable user-provided dataset., Comment: 13 pages, 9 figures