Source author record

Bryan Ostdiek

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

hep-ph hep-ex astro-ph.CO astro-ph.GA physics.data-an Machine Learning astro-ph.IM

Catalog footprint

Research overview

18works

7topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

preprint2022arXiv

Chasing Accreted Structures within Gaia DR2 using Deep Learning

In previous work, we developed a deep neural network classifier that only relies on phase-space information to obtain a catalog of accreted stars based on the second data release of Gaia (DR2). In this paper, we apply two clustering algorithms to identify velocity substructure within this catalog. We focus on the subset of stars with line-of-sight velocity measurements that fall in the range of Galactocentric radii $r \in [6.5, 9.5]$ kpc and vertical distances $|z| < 3$ kpc. Known structures such as Gaia Enceladus and the Helmi stream are identified. The largest previously-unknown structure, Nyx, is a vast stream consisting of at least 90 stars in the region of interest. This study displays the power of the machine learning approach by not only successfully identifying known features, but also discovering new kinematic structures that may shed light on the merger history of the Milky Way.

preprint2022arXiv

Creating Simple, Interpretable Anomaly Detectors for New Physics in Jet Substructure

Anomaly detection with convolutional autoencoders is a popular method to search for new physics in a model-agnostic manner. These techniques are powerful, but they are still a "black box," since we do not know what high-level physical observables determine how anomalous an event is. To address this, we adapt a recently proposed technique by Faucett et al., which maps out the physical observables learned by a neural network classifier, to the case of anomaly detection. We propose two different strategies that use a small number of high-level observables to mimic the decisions made by the autoencoder on background events, one designed to directly learn the output of the autoencoder, and the other designed to learn the difference between the autoencoder's outputs on a pair of events. Despite the underlying differences in their approach, we find that both strategies have similar ordering performance as the autoencoder and independently use the same six high-level observables. From there, we compare the performance of these networks as anomaly detectors. We find that both strategies perform similarly to the autoencoder across a variety of signals, giving a nontrivial demonstration that learning to order background events transfers to ordering a variety of signal events.

preprint2022arXiv

Evidence for a Vast Prograde Stellar Stream in the Solar Vicinity

Massive dwarf galaxies that merge with the Milky Way on prograde orbits can be dragged into the disk plane before being completely disrupted. Such mergers can contribute to an accreted stellar disk and a dark matter disk. We present evidence for Nyx, a vast new stellar stream in the vicinity of the Sun, that may provide the first indication that such an event occurred in the Milky Way. We identify about 500 stars that have coherent radial and prograde motion in this stream using a catalog of accreted stars built by applying deep learning methods to the second Gaia data release. Nyx is concentrated within $\pm 2$ kpc of the Galactic midplane and spans the full radial range studied (6.5-9.5 kpc). The kinematics of Nyx stars are distinct from those of both the thin and thick disk. In particular, its rotational speed lags the disk by $\sim 80$ km/s and its stars follow more eccentric orbits. A small number of Nyx stars have chemical abundances or inferred ages; from these, we deduce that Nyx stars have a peak metallicity of [Fe/H] $\sim -0.5$ and ages $\sim$ 10-13 Gyr. Taken together with the kinematic observations, these results strongly favor the interpretation that Nyx is the remnant of a disrupted dwarf galaxy. To further justify this interpretation, we explicitly demonstrate that metal-rich, prograde streams like Nyx can be found in the disk plane of Milky Way-like galaxies using the FIRE hydrodynamic simulations. Future spectroscopic studies will be able to validate whether Nyx stars originate from a single progenitor.

preprint2022arXiv

Extracting the Subhalo Mass Function from Strong Lens Images with Image Segmentation

Detecting substructure within strongly lensed images is a promising route to shed light on the nature of dark matter. However, it is a challenging task, which traditionally requires detailed lens modeling and source reconstruction, taking weeks to analyze each system. We use machine-learning to circumvent the need for lens and source modeling and develop a neural network to both locate subhalos in an image as well as determine their mass using the technique of image segmentation. The network is trained on images with a single subhalo located near the Einstein ring across a wide range of apparent source magnitudes. The network is then able to resolve subhalos with masses $m\gtrsim 10^{8.5} M_{\odot}$ . Training in this way allows the network to learn the gravitational lensing of light, and remarkably, it is then able to detect entire populations of substructure, even for locations further away from the Einstein ring than those used in training. Over a wide range of the apparent source magnitude, the false-positive rate is around three false subhalos per 100 images, coming mostly from the lightest detectable subhalo for that signal-to-noise ratio. With good accuracy and a low false-positive rate, counting the number of pixels assigned to each subhalo class over multiple images allows for a measurement of the subhalo mass function (SMF). When measured over three mass bins from $10^9M_{\odot}$ -- $10^{10} M_{\odot}$ the SMF slope is recovered with an error of 36% for 50 images, and this improves to 10% for 1000 images with Hubble Space Telescope-like noise.

preprint2022arXiv

Image segmentation for analyzing galaxy-galaxy strong lensing systems

The goal of this paper is to develop a machine learning model to analyze the main gravitational lens and detect dark substructure (subhalos) within simulated images of strongly lensed galaxies. Using the technique of image segmentation, we turn the task of identifying subhalos into a classification problem, where we label each pixel in an image as coming from the main lens, a subhalo within a binned mass range, or neither. Our network is only trained on images with a single smooth lens and either zero or one subhalo near the Einstein ring. On an independent test set with lenses with large ellipticities, quadrupole and octopole moments, and for source apparent magnitudes between 17-25, the area of the main lens is recovered accurately. On average, only 1.3% of the true area is missed and 1.2% of the true area is added to another part of the lens. In addition, subhalos as light as $10^{8.5}M_{\odot}$ can be detected if they lie in bright pixels along the Einstein ring. Furthermore, the model is able to generalize to new contexts it has not been trained on, such as locating multiple subhalos with varying masses or more than one large smooth lens.

preprint2022arXiv

Substructure Detection Reanalyzed: Dark Perturber shown to be a Line-of-Sight Halo

Observations of structure at sub-galactic scales are crucial for probing the properties of dark matter, which is the dominant source of gravity in the universe. It will become increasingly important for future surveys to distinguish between line-of-sight halos and subhalos to avoid wrong inferences on the nature of dark matter. We reanalyze a sub-galactic structure (in lens JVAS B1938+666) that has been previously found using the gravitational imaging technique in galaxy-galaxy lensing systems. This structure has been assumed to be a satellite in the halo of the main lens galaxy. We fit the redshift of the perturber of the system as a free parameter, using the multi-plane thin-lens approximation, and find that the redshift of the perturber is $z_\mathrm{int} = 1.42\substack{+0.10 \\ -0.15}$ (with a main lens redshift of $z=0.881$ ). Our analysis indicates that this structure is more massive than the previous result by an order of magnitude. This constitutes the first dark perturber shown to be a line-of-sight halo with a gravitational lensing method.

preprint2021arXiv

Challenges for Unsupervised Anomaly Detection in Particle Physics

Anomaly detection relies on designing a score to determine whether a particular event is uncharacteristic of a given background distribution. One way to define a score is to use autoencoders, which rely on the ability to reconstruct certain types of data (background) but not others (signals). In this paper, we study some challenges associated with variational autoencoders, such as the dependence on hyperparameters and the metric used, in the context of anomalous signal (top and $W$ ) jets in a QCD background. We find that the hyperparameter choices strongly affect the network performance and that the optimal parameters for one signal are non-optimal for another. In exploring the networks, we uncover a connection between the latent space of a variational autoencoder trained using mean-squared-error and the optimal transport distances within the dataset. We then show that optimal transport distances to representative events in the background dataset can be used directly for anomaly detection, with performance comparable to the autoencoders. Whether using autoencoders or optimal transport distances for anomaly detection, we find that the choices that best represent the background are not necessarily best for signal identification. These challenges with unsupervised anomaly detection bolster the case for additional exploration of semi-supervised or alternative approaches.

preprint2021arXiv

Deep Set Auto Encoders for Anomaly Detection in Particle Physics

There is an increased interest in model agnostic search strategies for physics beyond the standard model at the Large Hadron Collider. We introduce a Deep Set Variational Autoencoder and present results on the Dark Machines Anomaly Score Challenge. We find that the method attains the best anomaly detection ability when there is no decoding step for the network, and the anomaly score is based solely on the representation within the encoded latent space. This method was one of the top-performing models in the Dark Machines Challenge, both for the open data sets as well as the blinded data sets.

preprint2020arXiv

Cataloging Accreted Stars within Gaia DR2 using Deep Learning

The goal of this study is to present the development of a machine learning based approach that utilizes phase space alone to separate the Gaia DR2 stars into two categories: those accreted onto the Milky Way from those that are in situ. Traditional selection methods that have been used to identify accreted stars typically rely on full 3D velocity, metallicity information, or both, which significantly reduces the number of classifiable stars. The approach advocated here is applicable to a much larger portion of Gaia DR2. A method known as "transfer learning" is shown to be effective through extensive testing on a set of mock Gaia catalogs that are based on the FIRE cosmological zoom-in hydrodynamic simulations of Milky Way-mass galaxies. The machine is first trained on simulated data using only 5D kinematics as inputs and is then further trained on a cross-matched Gaia/RAVE data set, which improves sensitivity to properties of the real Milky Way. The result is a catalog that identifies around 767,000 accreted stars within Gaia DR2. This catalog can yield empirical insights into the merger history of the Milky Way and could be used to infer properties of the dark matter distribution.

preprint2020arXiv

On the ATLAS Top Mass Measurements and the Potential for Stealth Stop Contamination

The discovery of the stop - the Supersymmetric partner of the top quark - is a key goal of the physics program enabled by the Large Hadron Collider. Although much of the accessible parameter space has already been probed, all current searches assume the top mass is known. This is relevant for the "stealth stop" regime, which is characterized by decay kinematics that force the final state top quark off its mass shell; such decays would contaminate the top mass measurements. We investigate the resulting bias imparted to the template method based ATLAS approach. A careful recasting of these results shows that effect can be as large as 2.0 GeV, comparable to the current quoted uncertainty on the top mass. Thus, a robust exploration of the stealth stop splinter requires the simultaneous consideration of the impact on the top mass. Additionally, we explore the robustness of the template technique, and point out a simple strategy for improving the methodology implemented for the semi-leptonic channel.

preprint2019arXiv

Mass Agnostic Jet Taggers

Searching for new physics in large data sets needs a balance between two competing effects---signal identification vs background distortion. In this work, we perform a systematic study of both single variable and multivariate jet tagging methods that aim for this balance. The methods preserve the shape of the background distribution by either augmenting the training procedure or the data itself. Multiple quantitative metrics to compare the methods are considered, for tagging 2-, 3-, or 4-prong jets from the QCD background. This is the first study to show that the data augmentation techniques of Planing and PCA based scaling deliver similar performance as the augmented training techniques of Adversarial NN and uBoost, but are both easier to implement and computationally cheaper.

preprint2016arXiv

750 GeV Diphotons from Supersymmetry with Dirac Gauginos

Motivated by the recent excess in the diphoton invariant mass near 750 GeV, we explore a supersymmetric extension of the Standard Model that includes the minimal set of superpartners as well as additional Dirac partner chiral superfields in the adjoint representation for each gauge group. The bino partner pseudoscalar is identified as the 750 GeV resonance, while superpotential interactions between it and the gluino (wino) partners yield production via gluon fusion (decay to photon pairs) at one-loop. The gauginos and these additional adjoint superpartners are married by a Dirac mass and must also have Majorana masses. While a large wino partner Majorana mass is necessary to explain the excess, the gluino can be approximately Dirac-like, providing benefits consistent with being both "supersoft" (loop corrections to the scalar masses from Dirac gauginos are free of logarithmic enhancements) and "supersafe" (the experimental limits on the squark/gluino masses can be relaxed due to the reduced production rate). Consistency with the measured Standard Model-like Higgs boson mass is imposed, and a numerical exploration of the parameter space is provided. Models that can account for the diphoton excess are additionally characterized by having couplings that can remain perturbative up to very high scales, while remaining consistent with experimental constraints, the Higgs boson mass, and an electroweak scale which is not excessively fine tuned.

preprint2015arXiv

Constraining the minimal dark matter fiveplet with LHC searches

Adding a fermion to the standard model particle content which is a fiveplet under $SU(2)_L$ gives a dark matter candidate. The lightest state in the multiplet is neutral and automatically stable. The charged component of the multiplet is only slightly heavier and can travel a finite distance in the LHC detectors before decaying, leaving a charged track which disappears before the edge of the detector. We use the recent ATLAS and CMS searches to exclude a Majorana fiveplet with a mass up to 267 GeV. We estimate that with 3 ab $^{-1}$ of $\sqrt{s}=14$ TeV data this could be pushed to a mass 520 GeV. These exclusions are $\sim 10\%$ greater than what is achieved for wino-like dark matter. We also discuss how the doubly charged states could be used to distinguish a disappearing track signal from that given by a triplet such as the pure wino.

preprint2015arXiv

Dark Matter from the Supersymmetric Custodial Triplet Model

The Supersymmetric Custodial Triplet Model (SCTM) adds to the particle content of the MSSM three $SU(2)_L$ triplet chiral superfields with hypercharge $Y=(0,\pm1)$ . At the superpotential level the model respects a global $SU(2)_L \otimes SU(2)_R$ symmetry only broken by the Yukawa interactions. The pattern of vacuum expectation values of the neutral doublet and triplet scalar fields depends on the symmetry pattern of the Higgs soft breaking masses. We study the cases where this symmetry is maintained in the Higgs sector, and when it is broken only by the two doublets attaining different vacuum expectation values. In the former case, the symmetry is spontaneously broken down to the vectorial subgroup $SU(2)_V$ and the $ρ$ parameter is protected by the custodial symmetry. However in both situations the $ρ$ parameter is protected at tree level, allowing for light triplet scalars with large vacuum expectation values. We find that over a large range of parameter space, a light neutralino can supply the correct relic abundance of dark matter either through resonant s-channel triplet scalar funnels or well tempering of the Bino with the triplet fermions. Direct detection experiments have trouble probing these model points because the custodial symmetry suppresses the coupling of the neutralino and the $Z$ and a small Higgsino component of the neutralino suppresses the coupling with the Higgs. Likewise the annihilation cross sections for indirect detection lie below the Fermi-LAT upper bounds for the different channels.

preprint2015arXiv

The Relic Neutralino Surface at a 100 TeV collider

We map the parameter space for MSSM neutralino dark matter which freezes out to the observed relic abundance, in the limit that all superpartners except the neutralinos and charginos are decoupled. In this space of relic neutralinos, we show the dominant dark matter annihilation modes, the mass splittings among the electroweakinos, direct detection rates, and collider cross-sections. The mass difference between the dark matter and the next-to-lightest neutral and charged states is typically much less than electroweak gauge boson masses. With these small mass differences, the relic neutralino surface is accessible to a future 100 TeV hadron collider, which can discover inter-neutralino mass splittings down to 1 GeV and thermal relic dark matter neutralino masses up to 1.5 TeV with a few inverse attobarns of luminosity. This coverage is a direct consequence of the increased collider energy: the Standard Model events with missing transverse momentum in the TeV range have mostly hard electroweak radiation, distinct from the soft radiation shed in compressed electroweakino decays. We exploit this kinematic feature in final states including photons and leptons, tailored to the 100 TeV collider environment.

preprint2015arXiv

Towards the Final Word on Neutralino Dark Matter

We present a complete phenomenological prospectus for thermal relic neutralinos. Including Sommerfeld enhancements to relic abundance and halo annihilation calculations, we obtain direct, indirect, and collider discovery prospects for all neutralinos with mass parameters $M_1,M_2,|μ| < 4$ TeV, that freeze out to the observed dark matter abundance, with scalar superpartners decoupled. Much of the relic neutralino sector will be uncovered by the direct detection experiments Xenon1T and LZ, as well as indirect detection with CTA. We emphasize that thermal relic higgsinos will be found by next-generation direct detection experiments, so long as $M_{1,2} < 4$ TeV. Charged tracks at a 100 TeV hadron collider complement indirect searches for relic winos. Thermal relic bino-winos still evade all planned experiments, including disappearing charged-track searches. However, they can be discovered by compressed electroweakino searches at a 100 TeV collider, completing the full coverage of the relic neutralino surface.

preprint2014arXiv

Catching sparks from well-forged neutralinos

In this paper we present a new search technique for electroweakinos, the superpartners of electroweak gauge and Higgs bosons, based on final states with missing transverse energy, a photon, and a dilepton pair, $\ell^+\,\ell^- + γ+ \displaystyle{\not} E_T$ . Unlike traditional electroweakino searches, which perform best when $m_{\widetildeχ^0_{2,3}} - m_{\widetildeχ^0_1}, m_{\widetildeχ^{\pm}} - m_{\widetildeχ^0_1} > m_Z$ , our search favors nearly degenerate spectra; degenerate electroweakinos typically have a larger branching ratio to photons, and the cut $m_{\ell\ell} \ll m_Z$ effectively removes on-shell Z boson backgrounds while retaining the signal. This feature makes our technique optimal for `well-tempered' scenarios, where the dark matter relic abundance is achieved with inter-electroweakino splittings of $\sim 20 - 70\,\text{GeV}$ . Additionally, our strategy applies to a wider range of scenarios where the lightest neutralinos are almost degenerate, but only make up a subdominant component of the dark matter -- a spectrum we dub `well-forged'. Focusing on bino-Higgsino admixtures, we present optimal cuts and expected efficiencies for several benchmark scenarios. We find bino-Higgsino mixtures with $m_{\widetildeχ^0_{2,3}} \lesssim 190\,\text{GeV}$ and $m_{\widetildeχ^0_{2,3}} - m_{\widetildeχ^0_1} \cong 30\,\text{GeV}$ can be uncovered after roughly $600\,\text{fb}^{-1}$ of luminosity at the 14 TeV LHC. Scenarios with lighter states require less data for discovery, while scenarios with heavier states or larger mass splittings are harder to discriminate from the background and require more data. Unlike many searches for supersymmetry, electroweakino searches are one area where the high luminosity of the next LHC run, rather than the increased energy, is crucial for discovery.

preprint2012arXiv

LHC Signals of Non-Custodial Warped 5D Models

We study the implications at the LHC for a recent class of non-custodial warped extra-dimensional models where the AdS_5 metric is modified near the infrared brane. Such models allow for TeV Kaluza-Klein excitations without conflict with electroweak precision tests. We discuss both the production of electroweak and strong Kaluza-Klein gauge bosons. As we will show, only signals involving the third generation of quarks seem to be feasible in order to probe this scenario.

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Bryan Ostdiek

Research overview

Is this your work?

18 publications

Chasing Accreted Structures within Gaia DR2 using Deep Learning

Creating Simple, Interpretable Anomaly Detectors for New Physics in Jet Substructure

Evidence for a Vast Prograde Stellar Stream in the Solar Vicinity

Extracting the Subhalo Mass Function from Strong Lens Images with Image Segmentation

Image segmentation for analyzing galaxy-galaxy strong lensing systems

Substructure Detection Reanalyzed: Dark Perturber shown to be a Line-of-Sight Halo

Challenges for Unsupervised Anomaly Detection in Particle Physics

Deep Set Auto Encoders for Anomaly Detection in Particle Physics

Cataloging Accreted Stars within Gaia DR2 using Deep Learning

On the ATLAS Top Mass Measurements and the Potential for Stealth Stop Contamination

Mass Agnostic Jet Taggers

750 GeV Diphotons from Supersymmetry with Dirac Gauginos

Constraining the minimal dark matter fiveplet with LHC searches

Dark Matter from the Supersymmetric Custodial Triplet Model

The Relic Neutralino Surface at a 100 TeV collider

Towards the Final Word on Neutralino Dark Matter

Catching sparks from well-forged neutralinos

LHC Signals of Non-Custodial Warped 5D Models

See the researcher in context

Building this map preview