Showing 1–50 of 160 results for author: Green, S

Volumetric Primitives for Modeling and Rendering Scattering and Emissive Media

Authors: Jorge Condor, Sebastien Speierer, Lukas Bode, Aljaz Bozic, Simon Green, Piotr Didyk, Adrian Jarabo

Abstract: We propose a volumetric representation based on primitives to model scattering and emissive media. Accurate scene representations enabling efficient rendering are essential for many computer graphics applications. General and unified representations that can handle surface and volume-based representations simultaneously, allowing for physically accurate modeling, remain a research challenge. Inspi… ▽ More We propose a volumetric representation based on primitives to model scattering and emissive media. Accurate scene representations enabling efficient rendering are essential for many computer graphics applications. General and unified representations that can handle surface and volume-based representations simultaneously, allowing for physically accurate modeling, remain a research challenge. Inspired by recent methods for scene reconstruction that leverage mixtures of 3D Gaussians to model radiance fields, we formalize and generalize the modeling of scattering and emissive media using mixtures of simple kernel-based volumetric primitives. We introduce closed-form solutions for transmittance and free-flight distance sampling for 3D Gaussian kernels, and propose several optimizations to use our method efficiently within any off-the-shelf volumetric path tracer by leveraging ray tracing for efficiently querying the medium. We demonstrate our method as an alternative to other forms of volume modeling (e.g. voxel grid-based representations) for forward and inverse rendering of scattering media. Furthermore, we adapt our method to the problem of radiance field optimization and rendering, and demonstrate comparable performance to the state of the art, while providing additional flexibility in terms of performance and usability. △ Less

Submitted 24 May, 2024; originally announced May 2024.

Comments: 17 pages, 10 figures

ACM Class: I.3.2; I.3.3; I.3.6; I.3.5; I.3.7

Investigating the {Origin} of the Absorption-Line Variability in Narrow-Line Seyfert 1 Galaxy WPVS 007

Authors: Kaylie S. Green, Sarah C. Gallagher, Karen M. Leighly, Hyunseop Choi, Dirk Grupe, Donald M. Terndrup, Gordon T. Richards, S. Komossa

Abstract: Broad Absorption Line Quasars (BALQs) are actively accreting supermassive black holes that have strong outflows characterized by broad absorption lines in their rest-UV spectra. Variability in these absorption lines occurs over months to years depending on the source. WPVS 007, a low-redshift, low-luminosity Narrow-line Seyfert 1 (NLS1) shows strong variability over shorter timescales, providing a… ▽ More Broad Absorption Line Quasars (BALQs) are actively accreting supermassive black holes that have strong outflows characterized by broad absorption lines in their rest-UV spectra. Variability in these absorption lines occurs over months to years depending on the source. WPVS 007, a low-redshift, low-luminosity Narrow-line Seyfert 1 (NLS1) shows strong variability over shorter timescales, providing a unique opportunity to study the driving mechanism behind this variability that may mimic longer scale variability in much more massive quasars. We present the first variability study using {the} spectral synthesis code SimBAL, which provides velocity-resolved changes in physical conditions of the gas using constraints from multiple absorption lines. Overall, we find WPVS 007 to have a highly ionized outflow with a large mass-loss rate and kinetic luminosity. We determine the primary cause of the absorption-line variability in WPVS 007 to be a change in covering fraction of the continuum by the outflow. This study is the first SimBAL analysis where multiple epochs of observation were fit simultaneously, demonstrating the ability of SimBAL to use the time-domain as an additional constraint in spectral models. △ Less

Submitted 9 May, 2024; originally announced May 2024.

Comments: 26 pages, 3 tables, 11 figures

Journal ref: 2023, ApJ, 953, 186

Evidence for eccentricity in the population of binary black holes observed by LIGO-Virgo-KAGRA

Authors: Nihar Gupte, Antoni Ramos-Buades, Alessandra Buonanno, Jonathan Gair, M. Coleman Miller, Maximilian Dax, Stephen R. Green, Michael Pürrer, Jonas Wildberger, Jakob Macke, Bernhard Schölkopf

Abstract: Binary black holes (BBHs) in eccentric orbits produce distinct modulations the emitted gravitational waves (GWs). The measurement of orbital eccentricity can provide robust evidence for dynamical binary formation channels. We analyze 57 GW events from the first, second and third observing runs of the LIGO-Virgo-KAGRA (LVK) Collaboration using a multipolar aligned-spin inspiral-merger-ringdown wave… ▽ More Binary black holes (BBHs) in eccentric orbits produce distinct modulations the emitted gravitational waves (GWs). The measurement of orbital eccentricity can provide robust evidence for dynamical binary formation channels. We analyze 57 GW events from the first, second and third observing runs of the LIGO-Virgo-KAGRA (LVK) Collaboration using a multipolar aligned-spin inspiral-merger-ringdown waveform model with two eccentric parameters: eccentricity and relativistic anomaly. This is made computationally feasible with the machine-learning code DINGO which accelerates inference by 2-3 orders of magnitude compared to traditional inference. First, we find eccentric aligned-spin versus quasi-circular aligned-spin $\log_{10}$ Bayes factors of 1.84 to 4.75 (depending on the glitch mitigation) for GW200129, 3.0 for GW190701 and 1.77 for GW200208_22. We measure $e_{\text{gw}, 10Hzへるつ}$ to be $0.27_{-0.12}^{+0.10}$ to $0.17_{-0.13}^{+0.14}$ for GW200129, $0.35_{-0.11}^{+0.32}$ for GW190701 and $0.35_{-0.21}^{+0.18}$ for GW200208_22. Second, we find $\log_{10}$ Bayes factors between the eccentric aligned-spin versus quasi-circular precessing-spin hypothesis between 1.43 and 4.92 for GW200129, 2.61 for GW190701 and 1.23 for GW200208_22. Third, our analysis does not show evidence for eccentricity in GW190521, which has an eccentric aligned-spin against quasi-circular aligned-spin $\log_{10}$ Bayes factor of 0.04. Fourth, we estimate that if we neglect the spin-precession and use an astrophysical prior, the probability of one out of the 57 events being eccentric is greater than 99.5% or $(100 - 8.4 \times 10^{-4})$% (depending on the glitch mitigation). Fifth, we study the impact on parameter estimation when neglecting either eccentricity or higher modes in eccentric models. These results underscore the importance of including eccentric parameters in the characterization of BBHs for GW detectors. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Comments: 33 pages, 13 figures

EgoLifter: Open-world 3D Segmentation for Egocentric Perception

Authors: Qiao Gu, Zhaoyang Lv, Duncan Frost, Simon Green, Julian Straub, Chris Sweeney

Abstract: In this paper we present EgoLifter, a novel system that can automatically segment scenes captured from egocentric sensors into a complete decomposition of individual 3D objects. The system is specifically designed for egocentric data where scenes contain hundreds of objects captured from natural (non-scanning) motion. EgoLifter adopts 3D Gaussians as the underlying representation of 3D scenes and… ▽ More In this paper we present EgoLifter, a novel system that can automatically segment scenes captured from egocentric sensors into a complete decomposition of individual 3D objects. The system is specifically designed for egocentric data where scenes contain hundreds of objects captured from natural (non-scanning) motion. EgoLifter adopts 3D Gaussians as the underlying representation of 3D scenes and objects and uses segmentation masks from the Segment Anything Model (SAM) as weak supervision to learn flexible and promptable definitions of object instances free of any specific object taxonomy. To handle the challenge of dynamic objects in ego-centric videos, we design a transient prediction module that learns to filter out dynamic objects in the 3D reconstruction. The result is a fully automatic pipeline that is able to reconstruct 3D object instances as collections of 3D Gaussians that collectively compose the entire scene. We created a new benchmark on the Aria Digital Twin dataset that quantitatively demonstrates its state-of-the-art performance in open-world 3D segmentation from natural egocentric input. We run EgoLifter on various egocentric activity datasets which shows the promise of the method for 3D egocentric perception at scale. △ Less

Submitted 26 March, 2024; originally announced March 2024.

Comments: Preprint. Project page: https://egolifter.github.io/

Ultralight vector dark matter search using data from the KAGRA O3GK run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi , et al. (1778 additional authors not shown)

Abstract: Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we prese… ▽ More Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we present the result of a search for $U(1)_{B-L}$ gauge boson DM using the KAGRA data from auxiliary length channels during the first joint observation run together with GEO600. By applying our search pipeline, which takes into account the stochastic nature of ultralight DM, upper bounds on the coupling strength between the $U(1)_{B-L}$ gauge boson and ordinary matter are obtained for a range of DM masses. While our constraints are less stringent than those derived from previous experiments, this study demonstrates the applicability of our method to the lower-mass vector DM search, which is made difficult in this measurement by the short observation time compared to the auto-correlation time scale of DM. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 20 pages, 5 figures

Report number: LIGO-P2300250

Gravitational wave populations and cosmology with neural posterior estimation

Authors: Konstantin Leyde, Stephen R. Green, Alexandre Toubiana, Jonathan Gair

Abstract: We apply neural posterior estimation for fast-and-accurate hierarchical Bayesian inference of gravitational wave populations. We use a normalizing flow to estimate directly the population hyper-parameters from a collection of individual source observations. This approach provides complete freedom in event representation, automatic inclusion of selection effects, and (in contrast to likelihood esti… ▽ More We apply neural posterior estimation for fast-and-accurate hierarchical Bayesian inference of gravitational wave populations. We use a normalizing flow to estimate directly the population hyper-parameters from a collection of individual source observations. This approach provides complete freedom in event representation, automatic inclusion of selection effects, and (in contrast to likelihood estimation) without the need for stochastic samplers to obtain posterior samples. Since the number of events may be unknown when the network is trained, we split into sub-population analyses that we later recombine; this allows for fast sequential analyses as additional events are observed. We demonstrate our method on a toy problem of dark siren cosmology, and show that inference takes just a few minutes and scales to $\sim 600$ events before performance degrades. We argue that neural posterior estimation therefore represents a promising avenue for population inference with large numbers of events. △ Less

Submitted 6 April, 2024; v1 submitted 20 November, 2023; originally announced November 2023.

Comments: 16 + 5 pages, 7 + 1 figures. Small changes according to the published version

Journal ref: Phys. Rev. D 109, 064056 (2024)

VLT/MUSE Characterisation of Dimorphos Ejecta from the DART Impact

Authors: Brian P. Murphy, Cyrielle Opitom, Colin Snodgrass, Matthew M. Knight, Jian-Yang Li, Nancy L. Chabot, Andrew S. Rivkin, Simon F. Green, Paloma Guetzoyan, Daniel Gardener, Julia de León

Abstract: We have observed the Didymos-Dimorphos binary system with the MUSE integral field unit spectrograph mounted at the Very Large Telescope (VLT) pre and post-DART impact, and captured the ensuing ejecta cone, debris cloud, and tails at sub-arcsecond resolutions. We targeted the Didymos system over 11 nights from 26 September to 25 October 2022, and utilized both narrow and wide-field observations wit… ▽ More We have observed the Didymos-Dimorphos binary system with the MUSE integral field unit spectrograph mounted at the Very Large Telescope (VLT) pre and post-DART impact, and captured the ensuing ejecta cone, debris cloud, and tails at sub-arcsecond resolutions. We targeted the Didymos system over 11 nights from 26 September to 25 October 2022, and utilized both narrow and wide-field observations with and without adaptive optics, respectively. We took advantage of the spectral-spatial coupled measurements and produced both white-light images and spectral maps of the dust reflectance. We identified and characterized numerous dust features, such as the ejecta cone, spirals, wings, clumps, and tails. We found that the base of the Sunward edge of the wings, from 03 to 19 October, consistent with maximum grain sizes on the order of 0.05-0.2 mm, and that the earliest detected clumps have the highest velocities on the order of 10 m/s. We also see that three clumps in narrow-field mode (8x8'') exhibit redder colors and slower speeds, around 0.09 m/s, than the surrounding ejecta, likely indicating that the clump is comprised of larger, slower grains. We measured the properties of the primary tail, and resolved and measured the properties of the secondary tail earlier than any other published study, with first retrieval on 03 October. Both tails exhibit similarities in curvature and relative flux, however, the secondary tail appears thinner, which may be caused by lower energy ejecta and possibly a low energy formation mechanism such as secondary impacts. △ Less

Submitted 16 November, 2023; originally announced November 2023.

Comments: 22 pages, 14 figures, accepted for publication in the Planetary Science Journal

Waveform Modelling for the Laser Interferometer Space Antenna

Authors: LISA Consortium Waveform Working Group, Niayesh Afshordi, Sarp Akçay, Pau Amaro Seoane, Andrea Antonelli, Josu C. Aurrekoetxea, Leor Barack, Enrico Barausse, Robert Benkel, Laura Bernard, Sebastiano Bernuzzi, Emanuele Berti, Matteo Bonetti, Béatrice Bonga, Gabriele Bozzola, Richard Brito, Alessandra Buonanno, Alejandro Cárdenas-Avendaño, Marc Casals, David F. Chernoff, Alvin J. K. Chua, Katy Clough, Marta Colleoni, Mekhi Dhesi, Adrien Druart , et al. (121 additional authors not shown)

Abstract: LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmologic… ▽ More LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmological distances; from the beginnings of inspirals that will venture into the ground-based detectors' view to the death spiral of compact objects into massive black holes, and many sources in between. Central to realising LISA's discovery potential are waveform models, the theoretical and phenomenological predictions of the pattern of gravitational waves that these sources emit. This white paper is presented on behalf of the Waveform Working Group for the LISA Consortium. It provides a review of the current state of waveform models for LISA sources, and describes the significant challenges that must yet be overcome. △ Less

Submitted 20 December, 2023; v1 submitted 2 November, 2023; originally announced November 2023.

Comments: 239 pages, 11 figures, white paper from the LISA Consortium Waveform Working Group, invited for submission to Living Reviews in Relativity, updated with comments from community

Relativistic perturbation theory for black-hole boson clouds

Authors: Enrico Cannizzaro, Laura Sberna, Stephen R. Green, Stefan Hollands

Abstract: We develop a relativistic perturbation theory for scalar clouds around rotating black holes. We first introduce a relativistic product and corresponding orthogonality relation between modes, extending a recent result for gravitational perturbations. We then derive the analog of time-dependent perturbation theory in quantum mechanics, and apply it to calculate self-gravitational frequency shifts. T… ▽ More We develop a relativistic perturbation theory for scalar clouds around rotating black holes. We first introduce a relativistic product and corresponding orthogonality relation between modes, extending a recent result for gravitational perturbations. We then derive the analog of time-dependent perturbation theory in quantum mechanics, and apply it to calculate self-gravitational frequency shifts. This approach supersedes the non-relativistic "gravitational atom" approximation, brings close agreement with numerical relativity, and has practical applications for gravitational-wave astronomy. △ Less

Submitted 3 January, 2024; v1 submitted 18 September, 2023; originally announced September 2023.

Comments: 7+4 pages, 2+1 figures. v2: improved derivation of the main result, including two new appendices, and other minor changes. matches version accepted for publication in PRL

Project Aria: A New Tool for Egocentric Multi-Modal AI Research

Authors: Jakob Engel, Kiran Somasundaram, Michael Goesele, Albert Sun, Alexander Gamino, Andrew Turner, Arjang Talattof, Arnie Yuan, Bilal Souti, Brighid Meredith, Cheng Peng, Chris Sweeney, Cole Wilson, Dan Barnes, Daniel DeTone, David Caruso, Derek Valleroy, Dinesh Ginjupalli, Duncan Frost, Edward Miller, Elias Mueggler, Evgeniy Oleinik, Fan Zhang, Guruprasad Somasundaram, Gustavo Solaira , et al. (49 additional authors not shown)

Abstract: Egocentric, multi-modal data as available on future augmented reality (AR) devices provides unique challenges and opportunities for machine perception. These future devices will need to be all-day wearable in a socially acceptable form-factor to support always available, context-aware and personalized AI applications. Our team at Meta Reality Labs Research built the Aria device, an egocentric, mul… ▽ More Egocentric, multi-modal data as available on future augmented reality (AR) devices provides unique challenges and opportunities for machine perception. These future devices will need to be all-day wearable in a socially acceptable form-factor to support always available, context-aware and personalized AI applications. Our team at Meta Reality Labs Research built the Aria device, an egocentric, multi-modal data recording and streaming device with the goal to foster and accelerate research in this area. In this paper, we describe the Aria device hardware including its sensor configuration and the corresponding software tools that enable recording and processing of such data. △ Less

Submitted 1 October, 2023; v1 submitted 24 August, 2023; originally announced August 2023.

Physical modelling of near-Earth asteroid (23187) 2000 PN9 with ground-based optical and radar observations

Authors: L. Dover, S. C. Lowry, A. Rożek, B. Rozitis, S. L. Jackson, T. Zegmott, Yu. N. Krugly, I. N. Belskaya, A. Fitzsimmons, S. F. Green, C. Snodgrass, P. R. Weissman, M. Brozović, L. A. M. Benner, M. W. Busch, V. R. Ayvazian, V. Chiorny, R. Ya. Inasaridze, M. Krugov, S. Mykhailova, I. Reva, J. Hibbert

Abstract: We present a physical model and spin-state analysis of the potentially hazardous asteroid (23187) 2000 PN9. As part of a long-term campaign to make direct detections of the YORP effect, we collected optical lightcurves of the asteroid between 2006 and 2020. These observations were combined with planetary radar data to develop a detailed shape model which was used to search for YORP acceleration. W… ▽ More We present a physical model and spin-state analysis of the potentially hazardous asteroid (23187) 2000 PN9. As part of a long-term campaign to make direct detections of the YORP effect, we collected optical lightcurves of the asteroid between 2006 and 2020. These observations were combined with planetary radar data to develop a detailed shape model which was used to search for YORP acceleration. We report that 2000 PN9 is a relatively large top-shaped body with a sidereal rotation period of 2.53216$\pm$0.00015 h. Although we find no evidence for rotational acceleration, YORP torques smaller than $\sim$10$^{-8}$$\,\rm rad/day^{2}$ cannot be ruled out. It is likely that 2000 PN9 is a YORP-evolved object, and may be an example of YORP equilibrium or self limitation. △ Less

Submitted 18 August, 2023; originally announced August 2023.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society. Main paper: 15 pages, 12 figures. Appendix: 9 pages, 8 figures

Search for Eccentric Black Hole Coalescences during the Third Observing Run of LIGO and Virgo

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi , et al. (1750 additional authors not shown)

Abstract: Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effect… ▽ More Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass $M>70$ $M_\odot$) binaries covering eccentricities up to 0.3 at 15 Hzへるつ orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities $0 < e \leq 0.3$ at $0.33$ Gpc$^{-3}$ yr$^{-1}$ at 90\% confidence level. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Comments: 24 pages, 5 figures

Report number: LIGO-P2300080

Data-driven reactivity prediction of targeted covalent inhibitors using computed quantum features for drug discovery

Authors: Tom W. A. Montgomery, Peter Pogány, Alice Purdy, Mike Harris, Marek Kowalik, Alex Ferraro, Hikmatyar Hasan, Darren V. S. Green, Sam Genway

Abstract: We present an approach to combine novel molecular features with experimental data within a data-driven pipeline. The method is applied to the challenge of predicting the reactivity of a series of sulfonyl fluoride molecular fragments used for drug discovery of targeted covalent inhibitors. We demonstrate utility in predicting reactivity using features extracted from a workflow which employs quantu… ▽ More We present an approach to combine novel molecular features with experimental data within a data-driven pipeline. The method is applied to the challenge of predicting the reactivity of a series of sulfonyl fluoride molecular fragments used for drug discovery of targeted covalent inhibitors. We demonstrate utility in predicting reactivity using features extracted from a workflow which employs quantum embedding of the reactive warhead using density matrix embedding theory, followed by Hamiltonian simulation of the resulting fragment model from an initial reference state. These predictions are found to improve when studying both larger active spaces and longer evolution times. The calculated features form a `quantum fingerprint' which allows molecules to be clustered with regard to warhead properties. We identify that the quantum fingerprint is well suited to scalable calculation on future quantum computing hardware, and explore approaches to capture results on current quantum hardware using error mitigation and suppression techniques. We further discuss how this general framework may be applied to a wider range of challenges where the potential for future quantum utility exists. △ Less

Submitted 18 July, 2023; originally announced July 2023.

Flow Matching for Scalable Simulation-Based Inference

Authors: Maximilian Dax, Jonas Wildberger, Simon Buchholz, Stephen R. Green, Jakob H. Macke, Bernhard Schölkopf

Abstract: Neural posterior estimation methods based on discrete normalizing flows have become established tools for simulation-based inference (SBI), but scaling them to high-dimensional problems can be challenging. Building on recent advances in generative modeling, we here present flow matching posterior estimation (FMPE), a technique for SBI using continuous normalizing flows. Like diffusion models, and… ▽ More Neural posterior estimation methods based on discrete normalizing flows have become established tools for simulation-based inference (SBI), but scaling them to high-dimensional problems can be challenging. Building on recent advances in generative modeling, we here present flow matching posterior estimation (FMPE), a technique for SBI using continuous normalizing flows. Like diffusion models, and in contrast to discrete flows, flow matching allows for unconstrained architectures, providing enhanced flexibility for complex data modalities. Flow matching, therefore, enables exact density evaluation, fast training, and seamless scalability to large architectures--making it ideal for SBI. We show that FMPE achieves competitive performance on an established SBI benchmark, and then demonstrate its improved scalability on a challenging scientific problem: for gravitational-wave inference, FMPE outperforms methods based on comparable discrete flows, reducing training time by 30% with substantially improved accuracy. Our work underscores the potential of FMPE to enhance performance in challenging inference scenarios, thereby paving the way for more advanced applications to scientific problems. △ Less

Submitted 27 October, 2023; v1 submitted 26 May, 2023; originally announced May 2023.

Comments: NeurIPS 2023. Code available at https://github.com/dingo-gw/flow-matching-posterior-estimation

Acceleration of a Positron Bunch in a Hollow Channel Plasma

Authors: Spencer Gessner, Erik Adli, James M. Allen, Weiming An, Christine I. Clarke, Chris E. Clayton, Sebastien Corde, Antoine Doche, Joel Frederico, Selina Z. Green, Mark J. Hogan, Chan Joshi, Carl A. Lindstrom, Michael Litos, Kenneth A. Marsh, Warren B. Mori, Brendan O'Shea, Navid Vafaei-Najafabadi, Vitaly Yakimenko

Abstract: Plasmas are a compelling medium for particle acceleration owing to their natural ability to sustain electric fields that are orders of magnitude larger than those available in conventional radio-frequency accelerators. Plasmas are also unique amongst accelerator technologies in that they respond differently to beams of opposite charge. The asymmetric response of a plasma to highly-relativistic ele… ▽ More Plasmas are a compelling medium for particle acceleration owing to their natural ability to sustain electric fields that are orders of magnitude larger than those available in conventional radio-frequency accelerators. Plasmas are also unique amongst accelerator technologies in that they respond differently to beams of opposite charge. The asymmetric response of a plasma to highly-relativistic electron and positron beams arises from the fact that plasmas are composed of light, mobile electrons and heavy, stationary ions. Hollow channel plasma acceleration is a technique for symmetrizing the response of the plasma, such that it works equally well for high-energy electron and positron beams. In the experiment described here, we demonstrate the generation of a positron beam-driven wake in an extended, annular plasma channel, and acceleration of a second trailing witness positron bunch by the wake. The leading bunch excites the plasma wakefield and loses energy to the plasma, while the witness bunch experiences an accelerating field and gains energy, thus providing a proof-of-concept for hollow channel acceleration of positron beams. At a bunch separation of 330 um, the accelerating gradient is 70 MV/m, the transformer ratio is 0.55, and the energy transfer efficiency is 18% for a drive-to-witness beam charge ratio of 5:1. △ Less

Submitted 30 December, 2023; v1 submitted 4 April, 2023; originally announced April 2023.

Automated control and optimisation of laser driven ion acceleration

Authors: B. Loughran, M. J. V. Streeter, H. Ahmed, S. Astbury, M. Balcazar, M. Borghesi, N. Bourgeois, C. B. Curry, S. J. D. Dann, S. DiIorio, N. P. Dover, T. Dzelzanis, O. C. Ettlinger, M. Gauthier, L. Giuffrida, G. D. Glenn, S. H. Glenzer, J. S. Green, R. J. Gray, G. S. Hicks, C. Hyland, V. Istokskaia, M. King, D. Margarone, O. McCusker , et al. (10 additional authors not shown)

Abstract: The interaction of relativistically intense lasers with opaque targets represents a highly non-linear, multi-dimensional parameter space. This limits the utility of sequential 1D scanning of experimental parameters for the optimisation of secondary radiation, although to-date this has been the accepted methodology due to low data acquisition rates. High repetition-rate (HRR) lasers augmented by ma… ▽ More The interaction of relativistically intense lasers with opaque targets represents a highly non-linear, multi-dimensional parameter space. This limits the utility of sequential 1D scanning of experimental parameters for the optimisation of secondary radiation, although to-date this has been the accepted methodology due to low data acquisition rates. High repetition-rate (HRR) lasers augmented by machine learning present a valuable opportunity for efficient source optimisation. Here, an automated, HRR-compatible system produced high fidelity parameter scans, revealing the influence of laser intensity on target pre-heating and proton generation. A closed-loop Bayesian optimisation of maximum proton energy, through control of the laser wavefront and target position, produced proton beams with equivalent maximum energy to manually-optimized laser pulses but using only 60% of the laser energy. This demonstration of automated optimisation of laser-driven proton beams is a crucial step towards deeper physical insight and the construction of future radiation sources. △ Less

Submitted 1 March, 2023; originally announced March 2023.

Comments: 11 pages

Characterizing the nucleus of comet 162P/Siding Spring using ground-based photometry

Authors: Abbie Donaldson, Rosita Kokotanekova, Agata Rożek, Colin Snodgrass, Daniel Gardener, Simon F. Green, Nafiseh Masoumzadeh, James Robinson

Abstract: Comet 162P/Siding Spring is a large Jupiter-family comet with extensive archival lightcurve data. We report new r-band nucleus lightcurves for this comet, acquired in 2018, 2021 and 2022. With the addition of these lightcurves, the phase angles at which the nucleus has been observed range from $0.39^\circ$ to $16.33^\circ$. We absolutely-calibrate the comet lightcurves to r-band Pan-STARRS 1 magni… ▽ More Comet 162P/Siding Spring is a large Jupiter-family comet with extensive archival lightcurve data. We report new r-band nucleus lightcurves for this comet, acquired in 2018, 2021 and 2022. With the addition of these lightcurves, the phase angles at which the nucleus has been observed range from $0.39^\circ$ to $16.33^\circ$. We absolutely-calibrate the comet lightcurves to r-band Pan-STARRS 1 magnitudes, and use these lightcurves to create a convex shape model of the nucleus by convex lightcurve inversion. The best-fitting shape model for 162P has axis ratios $a/b = 1.56$ and $b/c = 2.33$, sidereal period $P = 32.864\pm0.001$ h, and a rotation pole oriented towards ecliptic longitude $λらむだ_E = 118^\circ \pm 26^\circ$ and latitude $βべーた_E=-50^\circ\pm21^\circ$. We constrain the possible nucleus elongation to lie within $1.4 < a/b < 2.0$ and discuss tentative evidence that 162P may have a bilobed structure. Using the shape model to correct the lightcurves for rotational effects, we derive a linear phase function with slope $βべーた=0.051\pm0.002$ mag deg$^{-1}$ and intercept $H_r(1,1,0) = 13.86 \pm 0.02$ for 162P. We find no evidence that the nucleus exhibited an opposition surge at phase angles down to 0.39$^\circ$. The challenges associated with modelling the shapes of comet nuclei from lightcurves are highlighted, and we comment on the extent to which we anticipate that LSST will alleviate these challenges in the coming decade. △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: 14 pages, 10 figures, accepted for publication in MNRAS

Open data from the third observing run of LIGO, Virgo, KAGRA and GEO

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, A. Al-Jodah, C. Alléné, A. Allocca , et al. (1719 additional authors not shown)

Abstract: The global network of gravitational-wave observatories now includes five detectors, namely LIGO Hanford, LIGO Livingston, Virgo, KAGRA, and GEO 600. These detectors collected data during their third observing run, O3, composed of three phases: O3a starting in April of 2019 and lasting six months, O3b starting in November of 2019 and lasting five months, and O3GK starting in April of 2020 and lasti… ▽ More The global network of gravitational-wave observatories now includes five detectors, namely LIGO Hanford, LIGO Livingston, Virgo, KAGRA, and GEO 600. These detectors collected data during their third observing run, O3, composed of three phases: O3a starting in April of 2019 and lasting six months, O3b starting in November of 2019 and lasting five months, and O3GK starting in April of 2020 and lasting 2 weeks. In this paper we describe these data and various other science products that can be freely accessed through the Gravitational Wave Open Science Center at https://gwosc.org. The main dataset, consisting of the gravitational-wave strain time series that contains the astrophysical signals, is released together with supporting data useful for their analysis and documentation, tutorials, as well as analysis software packages. △ Less

Submitted 7 February, 2023; originally announced February 2023.

Comments: 27 pages, 3 figures

Report number: LIGO-P2200316

Multi-Agent Dynamic Pricing in a Blockchain Protocol Using Gaussian Bandits

Authors: Alexis Asseman, Tomasz Kornuta, Anirudh Patel, Matt Deible, Sam Green

Abstract: The Graph Protocol indexes historical blockchain transaction data and makes it available for querying. As the protocol is decentralized, there are many independent Indexers that index and compete with each other for serving queries to the Consumers. One dimension along which Indexers compete is pricing. In this paper, we propose a bandit-based algorithm for maximization of Indexers' revenue via Co… ▽ More The Graph Protocol indexes historical blockchain transaction data and makes it available for querying. As the protocol is decentralized, there are many independent Indexers that index and compete with each other for serving queries to the Consumers. One dimension along which Indexers compete is pricing. In this paper, we propose a bandit-based algorithm for maximization of Indexers' revenue via Consumer budget discovery. We present the design and the considerations we had to make for a dynamic pricing algorithm being used by multiple agents simultaneously. We discuss the results achieved by our dynamic pricing bandits both in simulation and deployed into production on one of the Indexers operating on Ethereum. We have open-sourced both the simulation framework and tools we created, which other Indexers have since started to adapt into their own workflows. △ Less

Submitted 6 January, 2023; v1 submitted 13 December, 2022; originally announced December 2022.

Topological defect coarsening in quenched smectic-C films analyzed using artificial neural networks

Authors: Ravin A. Chowdhury, Adam A. S. Green, Cheol S. Park, Joseph E. Maclennan, Noel A. Clark

Abstract: Mechanically quenching a thin film of smectic-C liquid crystal results in the formation of a dense array of thousands of topological defects in the director field. The subsequent rapid coarsening of the film texture by the mutual annihilation of defects of opposite sign has been captured using high-speed, polarized light video microscopy. The temporal evolution of the texture has been characterize… ▽ More Mechanically quenching a thin film of smectic-C liquid crystal results in the formation of a dense array of thousands of topological defects in the director field. The subsequent rapid coarsening of the film texture by the mutual annihilation of defects of opposite sign has been captured using high-speed, polarized light video microscopy. The temporal evolution of the texture has been characterized using an object-detection convolutional neural network to determine the defect locations, and a binary classification network customized to evaluate the brush orientation dynamics around the defects in order to determine their topological signs. At early times following the quench, inherent limits on the spatial resolution result in undercounting of the defects and deviations from expected behavior. At intermediate to late times, the observed annihilation dynamics scale in agreement with theoretical predictions and simulations of the $2$D XY model. △ Less

Submitted 30 November, 2022; originally announced December 2022.

Comments: 11 pages, 10 figures

Adapting to noise distribution shifts in flow-based gravitational-wave inference

Authors: Jonas Wildberger, Maximilian Dax, Stephen R. Green, Jonathan Gair, Michael Pürrer, Jakob H. Macke, Alessandra Buonanno, Bernhard Schölkopf

Abstract: Deep learning techniques for gravitational-wave parameter estimation have emerged as a fast alternative to standard samplers $\unicode{x2013}$ producing results of comparable accuracy. These approaches (e.g., DINGO) enable amortized inference by training a normalizing flow to represent the Bayesian posterior conditional on observed data. By conditioning also on the noise power spectral density (PS… ▽ More Deep learning techniques for gravitational-wave parameter estimation have emerged as a fast alternative to standard samplers $\unicode{x2013}$ producing results of comparable accuracy. These approaches (e.g., DINGO) enable amortized inference by training a normalizing flow to represent the Bayesian posterior conditional on observed data. By conditioning also on the noise power spectral density (PSD) they can even account for changing detector characteristics. However, training such networks requires knowing in advance the distribution of PSDs expected to be observed, and therefore can only take place once all data to be analyzed have been gathered. Here, we develop a probabilistic model to forecast future PSDs, greatly increasing the temporal scope of DINGO networks. Using PSDs from the second LIGO-Virgo observing run (O2) $\unicode{x2013}$ plus just a single PSD from the beginning of the third (O3) $\unicode{x2013}$ we show that we can train a DINGO network to perform accurate inference throughout O3 (on 37 real events). We therefore expect this approach to be a key component to enable the use of deep learning techniques for low-latency analyses of gravitational waves. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Conserved currents for Kerr and orthogonality of quasinormal modes

Authors: Stephen R. Green, Stefan Hollands, Laura Sberna, Vahid Toomani, Peter Zimmerman

Abstract: We introduce a bilinear form for Weyl scalar perturbations of Kerr. The form is symmetric and conserved, and we show that, when combined with a suitable renormalization prescription involving complex r integration contours, quasinormal modes are orthogonal in the bilinear form for different (l, m, n). These properties are apparently not evident consequences of standard properties for the radial an… ▽ More We introduce a bilinear form for Weyl scalar perturbations of Kerr. The form is symmetric and conserved, and we show that, when combined with a suitable renormalization prescription involving complex r integration contours, quasinormal modes are orthogonal in the bilinear form for different (l, m, n). These properties are apparently not evident consequences of standard properties for the radial and angular solutions to the decoupled Teukolsky relations and rely on the Petrov type D character of Kerr and its t-$φふぁい$ reflection isometry. We show that quasinormal mode excitation coefficients are given precisely by the projection with respect to our bilinear form. These properties can make our bilinear form useful to set up a framework for nonlinear quasinormal mode coupling in Kerr. We also provide a general discussion on conserved local currents and their associated local symmetry operators for metric and Weyl perturbations, identifying a collection containing an increasing number of derivatives. △ Less

Submitted 19 February, 2024; v1 submitted 28 October, 2022; originally announced October 2022.

Comments: 14+11 pages, 2 figures. v3: fixed typos in Eq. (55) compared to the version accepted to PRD

Neural Importance Sampling for Rapid and Reliable Gravitational-Wave Inference

Authors: Maximilian Dax, Stephen R. Green, Jonathan Gair, Michael Pürrer, Jonas Wildberger, Jakob H. Macke, Alessandra Buonanno, Bernhard Schölkopf

Abstract: We combine amortized neural posterior estimation with importance sampling for fast and accurate gravitational-wave inference. We first generate a rapid proposal for the Bayesian posterior using neural networks, and then attach importance weights based on the underlying likelihood and prior. This provides (1) a corrected posterior free from network inaccuracies, (2) a performance diagnostic (the sa… ▽ More We combine amortized neural posterior estimation with importance sampling for fast and accurate gravitational-wave inference. We first generate a rapid proposal for the Bayesian posterior using neural networks, and then attach importance weights based on the underlying likelihood and prior. This provides (1) a corrected posterior free from network inaccuracies, (2) a performance diagnostic (the sample efficiency) for assessing the proposal and identifying failure cases, and (3) an unbiased estimate of the Bayesian evidence. By establishing this independent verification and correction mechanism we address some of the most frequent criticisms against deep learning for scientific inference. We carry out a large study analyzing 42 binary black hole mergers observed by LIGO and Virgo with the SEOBNRv4PHM and IMRPhenomXPHM waveform models. This shows a median sample efficiency of $\approx 10\%$ (two orders-of-magnitude better than standard samplers) as well as a ten-fold reduction in the statistical uncertainty in the log evidence. Given these advantages, we expect a significant impact on gravitational-wave inference, and for this approach to serve as a paradigm for harnessing deep learning methods in scientific applications. △ Less

Submitted 30 May, 2023; v1 submitted 11 October, 2022; originally announced October 2022.

Comments: 8+7 pages, 1+5 figures. [v2]: Minor updates to match published version, code available at https://github.com/dingo-gw/dingo

Report number: LIGO-P2200297

Journal ref: Phys. Rev. Lett. 130, 171403 (2023)

Reconstruction of interactions in the ProtoDUNE-SP detector with Pandora

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, Z. Ahmad, J. Ahmed, B. Aimard, F. Akbar, B. Ali-Mohammadzadeh, K. Allison, S. Alonso Monsalve, M. AlRashed, C. Alt, A. Alton, R. Alvarez, P. Amedo , et al. (1203 additional authors not shown)

Abstract: The Pandora Software Development Kit and algorithm libraries provide pattern-recognition logic essential to the reconstruction of particle interactions in liquid argon time projection chamber detectors. Pandora is the primary event reconstruction software used at ProtoDUNE-SP, a prototype for the Deep Underground Neutrino Experiment far detector. ProtoDUNE-SP, located at CERN, is exposed to a char… ▽ More The Pandora Software Development Kit and algorithm libraries provide pattern-recognition logic essential to the reconstruction of particle interactions in liquid argon time projection chamber detectors. Pandora is the primary event reconstruction software used at ProtoDUNE-SP, a prototype for the Deep Underground Neutrino Experiment far detector. ProtoDUNE-SP, located at CERN, is exposed to a charged-particle test beam. This paper gives an overview of the Pandora reconstruction algorithms and how they have been tailored for use at ProtoDUNE-SP. In complex events with numerous cosmic-ray and beam background particles, the simulated reconstruction and identification efficiency for triggered test-beam particles is above 80% for the majority of particle type and beam momentum combinations. Specifically, simulated 1 GeV/$c$ charged pions and protons are correctly reconstructed and identified with efficiencies of 86.1$\pm0.6$% and 84.1$\pm0.6$%, respectively. The efficiencies measured for test-beam data are shown to be within 5% of those predicted by the simulation. △ Less

Submitted 17 July, 2023; v1 submitted 29 June, 2022; originally announced June 2022.

Comments: 39 pages, 20 figures. Accepted version. Published version available in Eur. Phys. J. C 83, 618 (2023) https://doi.org/10.1140/epjc/s10052-023-11733-2

Report number: FERMILAB-PUB-22-488-AD-ESH-LBNF-ND-SCD, CERN-EP-DRAFT-MISC-2022-007

Journal ref: Eur. Phys. J. C 83, 618 (2023)

Physical properties of near-Earth asteroid (2102) Tantalus from multi-wavelength observations

Authors: Agata Rożek, Stephen C. Lowry, Benjamin Rozitis, Lord R. Dover, Patrick A. Taylor, Anne Virkki, Simon F. Green, Colin Snodgrass, Alan Fitzsimmons, Justyn Campbell-White, Sedighe Sajadian, Valerio Bozza, Martin J. Burgdorf, Martin Dominik, R. Figuera Jaimes, Tobias C. Hinse, Markus Hundertmark, Uffe G. Jørgensen, Penélope Longa-Peña, Markus Rabus, Sohrab Rahvar, Jesper Skottfelt, John Southworth

Abstract: Between 2010 and 2017 we have collected new optical and radar observations of the potentially hazardous asteroid (2102)~Tantalus from the ESO NTT and Danish telescopes at the La Silla Observatory and from the Arecibo planetary radar. The object appears to be nearly spherical, showing a low amplitude light-curve variation and limited large-scale features in the radar images. The spin-state is diffi… ▽ More Between 2010 and 2017 we have collected new optical and radar observations of the potentially hazardous asteroid (2102)~Tantalus from the ESO NTT and Danish telescopes at the La Silla Observatory and from the Arecibo planetary radar. The object appears to be nearly spherical, showing a low amplitude light-curve variation and limited large-scale features in the radar images. The spin-state is difficult to constrain with the available data; including a certain light-curve subset significantly changes the spin-state estimates, and the uncertainties on period determination are significant. Constraining any change in rotation rate was not possible, despite decades of observations. The convex lightcurve-inversion model, with rotational pole at $λらむだ=210{\pm}41$° and $βべーた=-30{\pm}35$°, is more flattened than the two models reconstructed by including radar observations: with prograde ($λらむだ=36{\pm}23$°, $βべーた=30{\pm}15$°), and with retrograde rotation mode ($λらむだ=180{\pm}24$°, $βべーた=-30{\pm}16$°). Using data from WISE we were able to determine that the prograde model produces the best agreement in size determination between radar and thermophysical modelling. Radar measurements indicate possible variation in surface properties, suggesting one side might have lower radar albedo and be rougher at centimetre-to-decimetre scale than the other. However, further observations are needed to confirm this. Thermophysical analysis indicates a surface covered in fine-grained regolith, consistent with radar albedo and polarisation ratio measurements. Finally, geophysical investigation of the spin-stability of Tantalus shows that it could be exceeding its critical spin-rate via cohesive forces. △ Less

Submitted 28 June, 2022; originally announced June 2022.

Comments: Main paper: 14 pages and 12 figure, Appendix: 13 pages and 12 figues

Neutron spectroscopy with a high-pressure nitrogen-filled spherical proportional counter

Authors: I. Giomataris, S. Green, I. Katsioulas, P. Knights, I. Manthos, T. Neep, K. Nikolopoulos, T. Papaevangelou, B. Phoenix, J. Sanders, R. Ward

Abstract: The spherical proportional counter is a large volume gaseous detector which finds application in several fields, including direct Dark Matter searches. When the detector is filled with nitrogen it becomes an effective neutron spectrometer thanks to the $^{14}$N(n,$\mathrmαあるふぁ$)$^{11}$B and $^{14}$N(n,p)$^{14}$C reactions. Nitrogen, however, is a challenging operating gas for proportional counters and… ▽ More The spherical proportional counter is a large volume gaseous detector which finds application in several fields, including direct Dark Matter searches. When the detector is filled with nitrogen it becomes an effective neutron spectrometer thanks to the $^{14}$N(n,$\mathrmαあるふぁ$)$^{11}$B and $^{14}$N(n,p)$^{14}$C reactions. Nitrogen, however, is a challenging operating gas for proportional counters and requires a high electric field strength to gas pressure ratio. Benefiting from the latest advances in spherical proportional counter instrumentation and simulation techniques, we report first neutron measurements at operating pressures of up to 1.8 bar. This achievement enhances the prospects of the spherical proportional counter to act as a neutron spectrometer appropriate for challenging environments, including underground laboratories, and industrial and medical settings. △ Less

Submitted 13 February, 2023; v1 submitted 9 June, 2022; originally announced June 2022.

Comments: 15 pages, 24 figures

Journal ref: Nuclear Inst. and Methods in Physics Research, A 1049 (2023) 168124

Recurrent Video Restoration Transformer with Guided Deformable Attention

Authors: Jingyun Liang, Yuchen Fan, Xiaoyu Xiang, Rakesh Ranjan, Eddy Ilg, Simon Green, Jiezhang Cao, Kai Zhang, Radu Timofte, Luc Van Gool

Abstract: Video restoration aims at restoring multiple high-quality frames from multiple low-quality frames. Existing video restoration methods generally fall into two extreme cases, i.e., they either restore all frames in parallel or restore the video frame by frame in a recurrent way, which would result in different merits and drawbacks. Typically, the former has the advantage of temporal information fusi… ▽ More Video restoration aims at restoring multiple high-quality frames from multiple low-quality frames. Existing video restoration methods generally fall into two extreme cases, i.e., they either restore all frames in parallel or restore the video frame by frame in a recurrent way, which would result in different merits and drawbacks. Typically, the former has the advantage of temporal information fusion. However, it suffers from large model size and intensive memory consumption; the latter has a relatively small model size as it shares parameters across frames; however, it lacks long-range dependency modeling ability and parallelizability. In this paper, we attempt to integrate the advantages of the two cases by proposing a recurrent video restoration transformer, namely RVRT. RVRT processes local neighboring frames in parallel within a globally recurrent framework which can achieve a good trade-off between model size, effectiveness, and efficiency. Specifically, RVRT divides the video into multiple clips and uses the previously inferred clip feature to estimate the subsequent clip feature. Within each clip, different frame features are jointly updated with implicit feature aggregation. Across different clips, the guided deformable attention is designed for clip-to-clip alignment, which predicts multiple relevant locations from the whole inferred clip and aggregates their features by the attention mechanism. Extensive experiments on video super-resolution, deblurring, and denoising show that the proposed RVRT achieves state-of-the-art performance on benchmark datasets with balanced model size, testing memory and runtime. △ Less

Submitted 12 November, 2022; v1 submitted 5 June, 2022; originally announced June 2022.

Comments: Accepted by NeurIPS 2022. Code: https://github.com/JingyunLiang/RVRT

The Effect of Aspect Changes on Near-Earth Asteroid Phase Curves

Authors: Samuel L. Jackson, Benjamin Rozitis, Lord R. Dover, Simon F. Green, Ulrich C. Kolb, Allison E. Andrews, Stephen C. Lowry

Abstract: Phase curves of asteroids are typically considered to depend solely on the scattering properties of airless particulate surfaces and the size of the object being studied. In this study, we demonstrate the additional dependence of phase curves on object shape, rotation pole orientation, and viewing geometry over an apparition. Variations in the phase curve of near-Earth asteroid (159402) 1999 AP10… ▽ More Phase curves of asteroids are typically considered to depend solely on the scattering properties of airless particulate surfaces and the size of the object being studied. In this study, we demonstrate the additional dependence of phase curves on object shape, rotation pole orientation, and viewing geometry over an apparition. Variations in the phase curve of near-Earth asteroid (159402) 1999 AP10 over its apparition from July 2020 - January 2021 are verified to be due to aspect changes over the apparition. This is achieved through shape modelling of the asteroid and simulation of the phase curve over the apparition. We present simulations of asteroid phase curves over a range of geometries to understand the potential magnitude of this aspect effect, and under which circumstances it can begin to dominate in the phase curves. This dependence on aspect may introduce significant additional uncertainty in the properties derived from phase curve data. We provide and demonstrate software code to estimate the aspect-related uncertainty in near-Earth asteroid phase curves through simulation and model fitting of a randomly generated sample of ellipsoidal asteroid models over the observed viewing geometry. We demonstrate how ignoring this effect may lead to misleading interpretations of the data and underestimation of uncertainties in further studies, such as those in the infrared that use phase curve derived parameters when fitting physical properties of an asteroid. △ Less

Submitted 12 April, 2022; originally announced April 2022.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 14 pages, 11 figures, 5 tables

Equitable Continuous Organizations with Self-Assessed Valuations

Authors: Howard Heaton, Sam Green

Abstract: Organizations are often unable to align the interests of all stakeholders with the financial success of the organization (e.g. due to regulation). However, continuous organizations (COs) introduce a paradigm shift. COs offer immediate liquidity, are permission-less and can align incentives. CO shares are issued continuously in the form of tokens via a smart contract on a blockchain. Token prices a… ▽ More Organizations are often unable to align the interests of all stakeholders with the financial success of the organization (e.g. due to regulation). However, continuous organizations (COs) introduce a paradigm shift. COs offer immediate liquidity, are permission-less and can align incentives. CO shares are issued continuously in the form of tokens via a smart contract on a blockchain. Token prices are designed to increase as more tokens are minted. When the share supply is low, near-zero prices make it advantageous to buy and hold tokens until interest in the CO increases, enabling a profitable sale. This attribute of COs, known as investment efficiency, is desirable. Yet, it can yield allocative inefficiency via the "holdout problem," i.e. latecomers may find a CO more valuable than early tokenholders, but be unable to attain the same token holdings due to inflated prices. With the aim of increasing overall equity, we introduce a voting mechanism into COs. We show this balances allocative and investment efficiency, and may dissuade speculative trading behaviors, thereby decreasing investment risk. △ Less

Submitted 20 March, 2022; originally announced March 2022.

Thermal emission from bow shocks II: 3D magnetohydrodynamic models of Zeta Ophiuchi

Authors: S. Green, J. Mackey, P. Kavanagh, T. J. Haworth, M. Moutzouri, V. V. Gvaramdaze

Abstract: The nearby, massive, runaway star Zeta Ophiuchi has a large bow shock detected in optical and infrared, and, uniquely among runaway O stars, diffuse X-ray emission is detected from the shocked stellar wind. Here we make the first detailed computational investigation of the bow shock of Zeta Ophiuchi, to test whether a simple model of the bow shock can explain the observed nebula, and to compare th… ▽ More The nearby, massive, runaway star Zeta Ophiuchi has a large bow shock detected in optical and infrared, and, uniquely among runaway O stars, diffuse X-ray emission is detected from the shocked stellar wind. Here we make the first detailed computational investigation of the bow shock of Zeta Ophiuchi, to test whether a simple model of the bow shock can explain the observed nebula, and to compare the detected X-ray emission with simulated emission maps. We re-analysed archival {\it Chandra} observations of the thermal diffuse X-ray emission from the shocked wind region of the bow shock, finding total unabsorbed X-ray flux (0.3-2 keV band) corresponding to a diffuse luminosity of $L_\mathrm{X}=2.33~(0.79-3.45)\times10^{29}$ergs$^{-1}$. 3D MHD simulations were used to model the interaction of the star's wind with a uniform ISM using a range of stellar and ISM parameters motivated by observational constraints. Synthetic infrared, Ha, soft X-ray, emission measure, and radio 6\,GHz emission maps were generated from three simulations, for comparison with relevant observations. Simulations where the space velocity of Zeta Ophiuchi has a significant radial velocity produce infrared emission maps with opening angle of the bow shock in better agreement with observations than for the case where motion is fully in the plane of the sky. The simulation with the highest pressure has the closest match, with flux level within a factor of 2 of the observational lower limit, and emission weighted temperature of $\log_{10}(T_\mathrm{A}/\mathrm{K})=6.4$, although the morphology of the diffuse emission appears somewhat different. Observed X-ray emission is a filled bubble brightest near the star whereas simulations predict brightening towards the contact discontinuity as density increases. △ Less

Submitted 30 May, 2022; v1 submitted 11 March, 2022; originally announced March 2022.

Comments: 19 pages, 19 figures

Journal ref: A&A 665, A35 (2022)

PION: Simulations of Wind-Blown Nebulae

Authors: Jonathan Mackey, Samuel Green, Maria Moutzouri, Thomas J. Haworth, Robert D. Kavanagh, Maggie Celeste, Robert Brose, Davit Zargaryan, Ciarán O'Rourke

Abstract: We present an overview of PION, an open-source software project for solving radiation-magnetohydrodynamics equations on a nested grid, aimed at modelling asymmetric nebulae around massive stars. A new implementation of hybrid OpenMP/MPI parallel algorithms is briefly introduced, and improved scaling is demonstrated compared with the current release version. Three-dimensional simulations of an expa… ▽ More We present an overview of PION, an open-source software project for solving radiation-magnetohydrodynamics equations on a nested grid, aimed at modelling asymmetric nebulae around massive stars. A new implementation of hybrid OpenMP/MPI parallel algorithms is briefly introduced, and improved scaling is demonstrated compared with the current release version. Three-dimensional simulations of an expanding nebula around a Wolf-Rayet star are then presented and analysed, similar to previous 2D simulations in the literature. The evolution of the emission measure of the gas and the X-ray surface brightness are calculated as a function of time, and some qualitative comparison with observations is made. △ Less

Submitted 9 March, 2022; originally announced March 2022.

Comments: 6 pages, 3 figures, conference proceedings, IAUS362: The predictive power of computational astrophysics as discovery tool (D. Bisikalo, T. Hanawa, C. Boily, J. Stone, eds.)

Calibration of BAS-TR image plate response to GeV gold ions

Authors: D. Doria, P. Martin, H. Ahmed, A. Alejo, M. Cerchez, S. Ferguson, J. Fernandez-Tobias, J. S. Green, D. Gwynne, F. Hanton, J. Jarrett, D. A. Maclellan, A. McIlvenny, P. McKenna, J. A. Ruiz, M. Swantusch, O. Willi, S. Zhai, M. Borghesi, S. Kar

Abstract: The response of the BAS-TR image plate (IP) was absolutely calibrated using CR-39 track detector for high linear energy transfer (LET) Au ions up to $\sim$1.6 GeV (8.2 MeV/nucleon), accelerated by high-power lasers. The calibration was carried out by employing a high-resolution Thomson parabola spectrometer, which allowed resolving Au ions with closely spaced ionization states up to 58$^+$. A resp… ▽ More The response of the BAS-TR image plate (IP) was absolutely calibrated using CR-39 track detector for high linear energy transfer (LET) Au ions up to $\sim$1.6 GeV (8.2 MeV/nucleon), accelerated by high-power lasers. The calibration was carried out by employing a high-resolution Thomson parabola spectrometer, which allowed resolving Au ions with closely spaced ionization states up to 58$^+$. A response function was obtained by fitting the photo-stimulated luminescence (PSL) per Au ion for different ion energies, which is broadly in agreement with that expected from ion stopping in the active layer of the IP. This calibration would allow quantifying the ion energy spectra for high energy Au ions, which is important for further investigation of the laser-based acceleration of heavy ion beams. △ Less

Submitted 21 February, 2022; originally announced February 2022.

Comments: 6 pages; 4 figures. Submitted to journal for review

Journal ref: Review of Scientific Instruments 93, 033304 (2022)

Using machine learning to parametrize postmerger signals from binary neutron stars

Authors: Tim Whittaker, William E. East, Stephen R. Green, Luis Lehner, Huan Yang

Abstract: There is growing interest in the detection and characterization of gravitational waves from postmerger oscillations of binary neutron stars. These signals contain information about the nature of the remnant and the high-density and out-of-equilibrium physics of the postmerger processes, which would complement any electromagnetic signal. However, the construction of binary neutron star postmerger w… ▽ More There is growing interest in the detection and characterization of gravitational waves from postmerger oscillations of binary neutron stars. These signals contain information about the nature of the remnant and the high-density and out-of-equilibrium physics of the postmerger processes, which would complement any electromagnetic signal. However, the construction of binary neutron star postmerger waveforms is much more complicated than for binary black holes: (i) there are theoretical uncertainties in the neutron-star equation of state and other aspects of the high-density physics, (ii) numerical simulations are expensive and available ones only cover a small fraction of the parameter space with limited numerical accuracy, and (iii) it is unclear how to parametrize the theoretical uncertainties and interpolate across parameter space. In this work, we describe the use of a machine-learning method called a conditional variational autoencoder (CVAE) to construct postmerger models for hyper/massive neutron star remnant signals based on numerical-relativity simulations. The CVAE provides a probabilistic model, which encodes uncertainties in the training data within a set of latent parameters. We estimate that training such a model will ultimately require $\sim 10^4$ waveforms. However, using synthetic training waveforms as a proof-of-principle, we show that the CVAE can be used as an accurate generative model and that it encodes the equation of state in a useful latent representation. △ Less

Submitted 17 January, 2022; originally announced January 2022.

Journal ref: Phys. Rev. D 105, 124021 (2022)

Nonlinear effects in the black hole ringdown: absorption-induced mode excitation

Authors: Laura Sberna, Pablo Bosch, William E. East, Stephen R. Green, Luis Lehner

Abstract: Gravitational-wave observations of black hole ringdowns are commonly used to characterize binary merger remnants and to test general relativity. These analyses assume linear black hole perturbation theory, in particular that the ringdown can be described in terms of quasinormal modes even for times approaching the merger. Here we investigate a nonlinear effect during the ringdown, namely how a mod… ▽ More Gravitational-wave observations of black hole ringdowns are commonly used to characterize binary merger remnants and to test general relativity. These analyses assume linear black hole perturbation theory, in particular that the ringdown can be described in terms of quasinormal modes even for times approaching the merger. Here we investigate a nonlinear effect during the ringdown, namely how a mode excited at early times can excite additional modes as it is absorbed by the black hole. This is a third-order secular effect: the change in the black-hole mass causes a shift in the mode spectrum, so that the original mode is projected onto the new ones. Using nonlinear simulations, we study the ringdown of a spherically-symmetric scalar field around an asymptotically anti-de Sitter black hole, and we find that this "absorption-induced mode excitation" (AIME) is the dominant nonlinear effect. We show that this effect takes place well within the nonadiabatic regime, so we can analytically estimate it using a sudden mass-change approximation. Adapting our estimation technique to asymptotically-flat Schwarzschild black holes, we expect AIME to play a role in the analysis and interpretation of current and future gravitational wave observations. △ Less

Submitted 21 December, 2021; originally announced December 2021.

Comments: 15+4 pages, 9+4 figures

Journal ref: Phys. Rev. D 105, 064046 (2022)

Group equivariant neural posterior estimation

Authors: Maximilian Dax, Stephen R. Green, Jonathan Gair, Michael Deistler, Bernhard Schölkopf, Jakob H. Macke

Abstract: Simulation-based inference with conditional neural density estimators is a powerful approach to solving inverse problems in science. However, these methods typically treat the underlying forward model as a black box, with no way to exploit geometric properties such as equivariances. Equivariances are common in scientific models, however integrating them directly into expressive inference networks… ▽ More Simulation-based inference with conditional neural density estimators is a powerful approach to solving inverse problems in science. However, these methods typically treat the underlying forward model as a black box, with no way to exploit geometric properties such as equivariances. Equivariances are common in scientific models, however integrating them directly into expressive inference networks (such as normalizing flows) is not straightforward. We here describe an alternative method to incorporate equivariances under joint transformations of parameters and data. Our method -- called group equivariant neural posterior estimation (GNPE) -- is based on self-consistently standardizing the "pose" of the data while estimating the posterior over parameters. It is architecture-independent, and applies both to exact and approximate equivariances. As a real-world application, we use GNPE for amortized inference of astrophysical binary black hole systems from gravitational-wave observations. We show that GNPE achieves state-of-the-art accuracy while reducing inference times by three orders of magnitude. △ Less

Submitted 30 May, 2023; v1 submitted 25 November, 2021; originally announced November 2021.

Comments: 13+11 pages, 5+8 figures. [v2]: Minor updates to match published version, code available at https://github.com/dingo-gw/dingo

Journal ref: ICLR 2022

SatGen -- II. Assessing the impact of a disc potential on subhalo populations

Authors: Sheridan B. Green, Frank C. van den Bosch, Fangzhou Jiang

Abstract: The demographics of dark matter substructure depend sensitively on the nature of dark matter. Optimally leveraging this probe requires accurate theoretical predictions regarding the abundance of subhaloes. These predictions are hampered by artificial disruption in numerical simulations, by large halo-to-halo variance, and by the fact that the results depend on the baryonic physics of galaxy format… ▽ More The demographics of dark matter substructure depend sensitively on the nature of dark matter. Optimally leveraging this probe requires accurate theoretical predictions regarding the abundance of subhaloes. These predictions are hampered by artificial disruption in numerical simulations, by large halo-to-halo variance, and by the fact that the results depend on the baryonic physics of galaxy formation. In particular, numerical simulations have shown that the formation of a central disc can drastically reduce the abundance of substructure compared to a dark matter-only simulation, which has been attributed to enhanced destruction of substructure due to disc shocking. We examine the impact of discs on substructure using the semi-analytical subhalo model SatGen, which accurately models the tidal evolution of substructure free of the numerical disruption that still hampers $N$-body simulations. Using a sample of 10,000 merger trees of Milky Way-like haloes, we study the demographics of subhaloes that are evolved under a range of composite halo-disc potentials with unprecedented statistical power. We find that the overall subhalo abundance is relatively insensitive to properties of the disc aside from its total mass. For a disc that contains $5\%$ of $M_\mathrm{vir}$, the mean subhalo abundance within $r_\mathrm{vir}$ is suppressed by ${\lesssim}10\%$ relative to the no-disc case, a difference that is dwarfed by halo-to-halo variance. For the same disc mass, the abundance of subhaloes within 50 kpc is reduced by ${\sim}30\%$. We argue that the disc mainly drives excess mass loss for subhaloes with small pericentric radii and that the impact of disc shocking is negligible. △ Less

Submitted 25 October, 2021; originally announced October 2021.

Comments: 13 pages, 8 figures. Accepted to MNRAS

Journal ref: MNRAS 509, 2624 (2022)

Emulating Sunyaev-Zeldovich Images of Galaxy Clusters using Auto-Encoders

Authors: Tibor Rothschild, Daisuke Nagai, Han Aung, Sheridan B. Green, Michelle Ntampaka, John ZuHone

Abstract: We develop a machine learning algorithm that generates high-resolution thermal Sunyaev-Zeldovich (SZ) maps of novel galaxy clusters given only halo mass and mass accretion rate. The algorithm uses a conditional variational autoencoder (CVAE) in the form of a convolutional neural network and is trained with SZ maps generated from the IllustrisTNG simulation. Our method can reproduce many of the det… ▽ More We develop a machine learning algorithm that generates high-resolution thermal Sunyaev-Zeldovich (SZ) maps of novel galaxy clusters given only halo mass and mass accretion rate. The algorithm uses a conditional variational autoencoder (CVAE) in the form of a convolutional neural network and is trained with SZ maps generated from the IllustrisTNG simulation. Our method can reproduce many of the details of galaxy clusters that analytical models usually lack, such as internal structure and aspherical distribution of gas created by mergers, while achieving the same computational feasibility, allowing us to generate mock SZ maps for over $10^5$ clusters in 30 seconds on a laptop. We show that the model is capable of generating novel clusters (i.e. not found in the training set) and that the model accurately reproduces the effects of mass and mass accretion rate on the SZ images, such as scatter, asymmetry, and concentration, in addition to modeling merging sub-clusters. This work demonstrates the viability of machine-learning--based methods for producing the number of realistic, high-resolution maps of galaxy clusters necessary to achieve statistical constraints from future SZ surveys. △ Less

Submitted 9 February, 2022; v1 submitted 5 October, 2021; originally announced October 2021.

Comments: 13 pages, 12 figures, 1 tables, accepted for publication in MNRAS

Trajectory Planning with Deep Reinforcement Learning in High-Level Action Spaces

Authors: Kyle R. Williams, Rachel Schlossman, Daniel Whitten, Joe Ingram, Srideep Musuvathy, Anirudh Patel, James Pagan, Kyle A. Williams, Sam Green, Anirban Mazumdar, Julie Parish

Abstract: This paper presents a technique for trajectory planning based on continuously parameterized high-level actions (motion primitives) of variable duration. This technique leverages deep reinforcement learning (Deep RL) to formulate a policy which is suitable for real-time implementation. There is no separation of motion primitive generation and trajectory planning: each individual short-horizon motio… ▽ More This paper presents a technique for trajectory planning based on continuously parameterized high-level actions (motion primitives) of variable duration. This technique leverages deep reinforcement learning (Deep RL) to formulate a policy which is suitable for real-time implementation. There is no separation of motion primitive generation and trajectory planning: each individual short-horizon motion is formed during the Deep RL training to achieve the full-horizon objective. Effectiveness of the technique is demonstrated numerically on a well-studied trajectory generation problem and a planning problem on a known obstacle-rich map. This paper also develops a new loss function term for policy-gradient-based Deep RL, which is analogous to an anti-windup mechanism in feedback control. We demonstrate the inclusion of this new term in the underlying optimization increases the average policy return in our numerical example. △ Less

Submitted 12 August, 2022; v1 submitted 30 September, 2021; originally announced October 2021.

Journal ref: IEEE Transactions on Aerospace and Electronic Systems, 59 (2023) 2513-2529

Detection of the YORP Effect on the contact-binary (68346) 2001 KZ66 from combined radar and optical observations

Authors: Tarik J. Zegmott, S. C. Lowry, A. Rożek, B. Rozitis, M. C. Nolan, E. S. Howell, S. F. Green, C. Snodgrass, A. Fitzsimmons, P. R. Weissman

Abstract: The YORP effect is a small thermal-radiation torque experienced by small asteroids, and is considered to be crucial in their physical and dynamical evolution. It is important to understand this effect by providing measurements of YORP for a range of asteroid types to facilitate the development of a theoretical framework. We are conducting a long-term observational study on a selection of near-Eart… ▽ More The YORP effect is a small thermal-radiation torque experienced by small asteroids, and is considered to be crucial in their physical and dynamical evolution. It is important to understand this effect by providing measurements of YORP for a range of asteroid types to facilitate the development of a theoretical framework. We are conducting a long-term observational study on a selection of near-Earth asteroids to support this. We focus here on (68346) 2001 KZ66, for which we obtained both optical and radar observations spanning a decade. This allowed us to perform a comprehensive analysis of the asteroid's rotational evolution. Furthermore, radar observations from the Arecibo Observatory enabled us to generate a detailed shape model. We determined that (68346) is a retrograde rotator with its pole near the southern ecliptic pole, within a $ 15^\circ$ radius of longitude $ 170^\circ$ and latitude $ -85^\circ$. By combining our radar-derived shape model with the optical light curves we developed a refined solution to fit all available data, which required a YORP strength of $ (8.43\pm0.69)\times10^{-8} \rm~rad ~day^{-2} $. (68346) has a distinct bifurcated shape comprising a large ellipsoidal component joined by a sharp neckline to a smaller non-ellipsoidal component. This object likely formed from either the gentle merging of a binary system, or from the deformation of a rubble pile due to YORP spin-up. The shape exists in a stable configuration close to its minimum in topographic variation, where regolith is unlikely to migrate from areas of higher potential. △ Less

Submitted 30 August, 2021; originally announced August 2021.

Comments: To be published in the Monthly Notices of the Royal Astronomical Society (MNRAS). 19 pages, 14 figures

New metric reconstruction scheme for gravitational self-force calculations

Authors: Vahid Toomani, Peter Zimmerman, Andrew Spiers, Stefan Hollands, Adam Pound, Stephen R. Green

Abstract: Inspirals of stellar-mass objects into massive black holes will be important sources for the space-based gravitational-wave detector LISA. Modelling these systems requires calculating the metric perturbation due to a point particle orbiting a Kerr black hole. Currently, the linear perturbation is obtained with a metric reconstruction procedure that puts it in a "no-string" radiation gauge which is… ▽ More Inspirals of stellar-mass objects into massive black holes will be important sources for the space-based gravitational-wave detector LISA. Modelling these systems requires calculating the metric perturbation due to a point particle orbiting a Kerr black hole. Currently, the linear perturbation is obtained with a metric reconstruction procedure that puts it in a "no-string" radiation gauge which is singular on a surface surrounding the central black hole. Calculating dynamical quantities in this gauge involves a subtle procedure of "gauge completion" as well as cancellations of very large numbers. The singularities in the gauge also lead to pathological field equations at second perturbative order. In this paper we re-analyze the point-particle problem in Kerr using the corrector-field reconstruction formalism of Green, Hollands, and Zimmerman (GHZ). We clarify the relationship between the GHZ formalism and previous reconstruction methods, showing that it provides a simple formula for the "gauge completion". We then use it to develop a new method of computing the metric in a more regular gauge: a Teukolsky puncture scheme. This scheme should ameliorate the problem of large cancellations, and by constructing the linear metric perturbation in a sufficiently regular gauge, it should provide a first step toward second-order self-force calculations in Kerr. Our methods are developed in generality in Kerr, but we illustrate some key ideas and demonstrate our puncture scheme in the simple setting of a static particle in Minkowski spacetime. △ Less

Submitted 10 December, 2021; v1 submitted 9 August, 2021; originally announced August 2021.

Comments: 73 pages in CQG format, 2 figures. v2 has an expanded introduction

Journal ref: Class. Quantum Grav. 39, 015019 (2022)

Gamma-radiation sky maps from compact binaries

Authors: Néstor Ortiz, Federico Carrasco, Stephen R. Green, Luis Lehner, Steven L. Liebling, John Ryan Westernacher-Schneider

Abstract: We study sky maps and light curves of gamma-ray emission from neutron stars in compact binaries, and in isolation. We briefly review some gamma-ray emission models, and reproduce sky maps from a standard isolated pulsar in the Separatrix Layer model. We consider isolated pulsars with several variations of a dipole magnetic field, including superpositions, and predict their gamma-ray emission. Our… ▽ More We study sky maps and light curves of gamma-ray emission from neutron stars in compact binaries, and in isolation. We briefly review some gamma-ray emission models, and reproduce sky maps from a standard isolated pulsar in the Separatrix Layer model. We consider isolated pulsars with several variations of a dipole magnetic field, including superpositions, and predict their gamma-ray emission. Our results provide new heuristics on what can and cannot be inferred about the magnetic field configuration of pulsars from high-energy observations. We find that typical double-peak light curves can be produced by pulsars with significant multipole structure beyond a single dipole. For binary systems, we also present a simple approximation that is useful for rapid explorations of binary magnetic field structure. Finally, we predict the gamma-ray emission pattern from a compact black hole-neutron star binary moments before merger by applying the Separatrix Layer model to data simulated in full general relativity; we find that face-on observers receive little emission, equatorial observers see one broad peak, and more generic observers typically see two peaks. △ Less

Submitted 17 February, 2022; v1 submitted 14 July, 2021; originally announced July 2021.

Comments: 29 pages, 19 figures. Matches published version

Journal ref: JCAP02 (2022) 027

Asteroid Photometry with PIRATE: Optimizations and Techniques for Small Aperture Telescopes

Authors: Samuel L. Jackson, Ulrich C. Kolb, Simon F. Green

Abstract: Small aperture telescopes provide the opportunity to conduct high frequency, targeted observations of near-Earth Asteroids that are not feasible with larger facilities due to highly competitive time allocation requirements. Observations of asteroids with these types of facilities often focus on rotational brightness variations rather than longer-term phase angle dependent variations (phase curves)… ▽ More Small aperture telescopes provide the opportunity to conduct high frequency, targeted observations of near-Earth Asteroids that are not feasible with larger facilities due to highly competitive time allocation requirements. Observations of asteroids with these types of facilities often focus on rotational brightness variations rather than longer-term phase angle dependent variations (phase curves) due to the difficulty of achieving high precision photometric calibration. We have developed an automated asteroid light curve extraction and calibration pipeline for images of moving objects from the 0.43 m Physics Innovations Robotic Telescope Explorer (PIRATE). This allows for the frequency and quality of observations required to construct asteroid phase curves. Optimisations in standard data reduction procedures are identified that may allow for similar small aperture facilities, constructed from commercially available/off-the-shelf components, to improve image and subsequent data quality. A demonstration of the hardware and software capabilities is expressed through observation statistics from a 10 month observing campaign, and through the photometric characterisation of near-Earth Asteroids 8014 (1990 MF) and 19764 (2000 NF5). △ Less

Submitted 9 July, 2021; originally announced July 2021.

Comments: 41 pages, 15 figures, 2 tables, 1 appendix. Accepted for publication in PASP

Journal ref: PASP, 133, 075003 (2021)

Neutron spectroscopy with N$_2$-filled high-pressure large-volume spherical proportional counters

Authors: I. Giomataris, S. Green, I. Katsioulas, P. Knights, I. Manthos, J. Matthews, T. Neep, K. Nikolopoulos, T. Papaevangelou, B. Phoenix, J. Sanders, R. Ward

Abstract: Precise in-situ measurements of the neutron flux in underground laboratories is crucial for direct dark matter searches, as neutron induced backgrounds can mimic the typical dark matter signal. The development of a novel neutron spectroscopy technique using Spherical Proportional Counters is investigated. The detector is operated with nitrogen and is sensitive to both fast and thermal neutrons thr… ▽ More Precise in-situ measurements of the neutron flux in underground laboratories is crucial for direct dark matter searches, as neutron induced backgrounds can mimic the typical dark matter signal. The development of a novel neutron spectroscopy technique using Spherical Proportional Counters is investigated. The detector is operated with nitrogen and is sensitive to both fast and thermal neutrons through the $^{14}$N(n, $αあるふぁ$)$^{11}$B and $^{14}$N(n, p)$^{14}$C reactions. This method holds potential to be a safe, inexpensive, effective, and reliable alternative to $^3$He-based detectors. Measurements of fast and thermal neutrons from an Am-Be source with a Spherical Proportional Counter operated at pressures up to 2 bar at Birmingham are discussed. △ Less

Submitted 6 July, 2021; originally announced July 2021.

Comments: 4 pages, 5 figures, International Conference on Technology and Instrumentation in Particle Physics (TIPP2021)

Real-time gravitational-wave science with neural posterior estimation

Authors: Maximilian Dax, Stephen R. Green, Jonathan Gair, Jakob H. Macke, Alessandra Buonanno, Bernhard Schölkopf

Abstract: We demonstrate unprecedented accuracy for rapid gravitational-wave parameter estimation with deep learning. Using neural networks as surrogates for Bayesian posterior distributions, we analyze eight gravitational-wave events from the first LIGO-Virgo Gravitational-Wave Transient Catalog and find very close quantitative agreement with standard inference codes, but with inference times reduced from… ▽ More We demonstrate unprecedented accuracy for rapid gravitational-wave parameter estimation with deep learning. Using neural networks as surrogates for Bayesian posterior distributions, we analyze eight gravitational-wave events from the first LIGO-Virgo Gravitational-Wave Transient Catalog and find very close quantitative agreement with standard inference codes, but with inference times reduced from O(day) to a minute per event. Our networks are trained using simulated data, including an estimate of the detector-noise characteristics near the event. This encodes the signal and noise models within millions of neural-network parameters, and enables inference for any observed data consistent with the training distribution, accounting for noise nonstationarity from event to event. Our algorithm -- called "DINGO" -- sets a new standard in fast-and-accurate inference of physical parameters of detected gravitational-wave events, which should enable real-time data analysis without sacrificing accuracy. △ Less

Submitted 30 May, 2023; v1 submitted 23 June, 2021; originally announced June 2021.

Comments: 7+12 pages, 4+11 figures. [v2]: Minor updates to match published version, code available at https://github.com/dingo-gw/dingo

Report number: LIGO-P2100223

Journal ref: Phys.Rev.Lett. 127, 241103 (2021)

(6478) Gault: Physical characterization of an active main-belt asteroid

Authors: Maxime Devogèle, Marin Ferrais, Emmanuel Jehin, Nicholas Moskovitz, Brian A. Skiff, Stephen E. Levine, Annika Gustafsson, Davide Farnocchia, Marco Micheli, Colin Snodgrass, Galin Borisov, Jean Manfroid, Youssef Moulane, Zouhair Benkhaldoun, Artem Burdanov, Francisco J. Pozuelos, Michael Gillon, Julien de Wit, Simon F. Green, Philippe Bendjoya, Jean-Pierre Rivet, Luy Abe, David Vernet, Colin Orion Chandler, Chadwick A. Trujillo

Abstract: In December 2018, the main-belt asteroid (6478)~Gault was reported to display activity. Gault is an asteroid belonging to the Phocaea dynamical family and was not previously known to be active, nor was any other member of the Phocaea family. In this work we present the results of photometric and spectroscopic observations that commenced soon after the discovery of activity. We obtained observation… ▽ More In December 2018, the main-belt asteroid (6478)~Gault was reported to display activity. Gault is an asteroid belonging to the Phocaea dynamical family and was not previously known to be active, nor was any other member of the Phocaea family. In this work we present the results of photometric and spectroscopic observations that commenced soon after the discovery of activity. We obtained observations over two apparitions to monitor its activity, rotation period, composition, and possible non-gravitational orbital evolution. We find that Gault has a rotation period of $P = 2.4929 \pm 0.0003$ hours with a lightcurve amplitude of $0.06$ magnitude. This short rotation period close to the spin barrier limit is consistent with Gault having a density no smaller than $ρろー= 1.85$~g/cm$^3$ and its activity being triggered by the YORP spin-up mechanism. Analysis of the Gault phase curve over phase angles ranging from $0.4^{\circ}$ to $23.6^{\circ}$ provides an absolute magnitude of $H = 14.81 \pm 0.04$, $G1=0.25 \pm 0.07$, and $G2= 0.38 \pm 0.04$. Model fits to the phase curve find the surface regolith grain size constrained between 100-500 $\rmμみゅー$m. Using relations between the phase curve and albedo we determine that the geometrical albedo of Gault is $p_{\rm v} = 0.26 \pm 0.05$ corresponding to an equivalent diameter of $D = 2.8^{+0.4}_{-0.2}$ km. Our spectroscopic observations are all consistent with an ordinary chondrite-like composition (S, or Q-type in the Bus-DeMeo taxonomic classification). A search through archival photographic plate surveys found previously unidentified detections of Gault dating back to 1957 and 1958. Only the latter had been digitized, which we measured to nearly double the observation arc of Gault. Finally, we did not find any signal of activity during the 2020 apparition or non-gravitational effects on its orbit. △ Less

Submitted 22 April, 2021; originally announced April 2021.

Comments: 15 pages, 11 figures; Accepted in MNRAS

PION: Simulating bow shocks and circumstellar nebulae

Authors: Jonathan Mackey, Samuel Green, Maria Moutzouri, Thomas J. Haworth, Robert D. Kavanagh, Davit Zargaryan, Maggie Celeste

Abstract: Expanding nebulae are produced by mass loss from stars, especially during late stages of evolution. Multi-dimensional simulation of these nebulae requires high resolution near the star and permits resolution that decreases with distance from the star, ideally with adaptive timesteps. We report the implementation and testing of static mesh-refinement in the radiation-magnetohydrodynamics code PION,… ▽ More Expanding nebulae are produced by mass loss from stars, especially during late stages of evolution. Multi-dimensional simulation of these nebulae requires high resolution near the star and permits resolution that decreases with distance from the star, ideally with adaptive timesteps. We report the implementation and testing of static mesh-refinement in the radiation-magnetohydrodynamics code PION, and document its performance for 2D and 3D calculations. The bow shock produced by a hot, magnetized, slowly rotating star as it moves through the magnetized ISM is simulated in 3D, highlighting differences compared with 2D calculations. Latitude-dependent, time-varying magnetized winds are modelled and compared with simulations of ring nebulae around blue supergiants from the literature. A 3D simulation of the expansion of a fast wind from a Wolf-Rayet star into the slow wind from a previous red supergiant phase of evolution is presented, with results compared with results in the literature and analytic theory. Finally the wind-wind collision from a binary star system is modelled with 3D MHD, and the results compared with previous 2D hydrodynamic calculations. A python library is provided for reading and plotting simulation snapshots, and the generation of synthetic infrared emission maps using torus is also demonstrated. It is shown that state-of-the-art 3D MHD simulations of wind-driven nebulae can be performed using PION with reasonable computational resources. The source code and user documentation is made available for the community under a BSD3 licence. △ Less

Submitted 12 March, 2021; originally announced March 2021.

Comments: 27 pages, accepted for publication in MNRAS

Experiment Simulation Configurations Approximating DUNE TDR

Authors: DUNE Collaboration, B. Abi, R. Acciarri, M. A. Acero, G. Adamov, D. Adams, M. Adinolfi, Z. Ahmad, J. Ahmed, T. Alion, S. Alonso Monsalve, C. Alt, J. Anderson, C. Andreopoulos, M. P. Andrews, F. Andrianala, S. Andringa, A. Ankowski, M. Antonova, S. Antusch, A. Aranda-Fernandez, A. Ariga, L. O. Arnold, M. A. Arroyave, J. Asaadi , et al. (949 additional authors not shown)

Abstract: The Deep Underground Neutrino Experiment (DUNE) is a next-generation long-baseline neutrino oscillation experiment consisting of a high-power, broadband neutrino beam, a highly capable near detector located on site at Fermilab, in Batavia, Illinois, and a massive liquid argon time projection chamber (LArTPC) far detector located at the 4850L of Sanford Underground Research Facility in Lead, South… ▽ More The Deep Underground Neutrino Experiment (DUNE) is a next-generation long-baseline neutrino oscillation experiment consisting of a high-power, broadband neutrino beam, a highly capable near detector located on site at Fermilab, in Batavia, Illinois, and a massive liquid argon time projection chamber (LArTPC) far detector located at the 4850L of Sanford Underground Research Facility in Lead, South Dakota. The long-baseline physics sensitivity calculations presented in the DUNE Physics TDR, and in a related physics paper, rely upon simulation of the neutrino beam line, simulation of neutrino interactions in the near and far detectors, fully automated event reconstruction and neutrino classification, and detailed implementation of systematic uncertainties. The purpose of this posting is to provide a simplified summary of the simulations that went into this analysis to the community, in order to facilitate phenomenological studies of long-baseline oscillation at DUNE. Simulated neutrino flux files and a GLoBES configuration describing the far detector reconstruction and selection performance are included as ancillary files to this posting. A simple analysis using these configurations in GLoBES produces sensitivity that is similar, but not identical, to the official DUNE sensitivity. DUNE welcomes those interested in performing phenomenological work as members of the collaboration, but also recognizes the benefit of making these configurations readily available to the wider community. △ Less

Submitted 18 March, 2021; v1 submitted 8 March, 2021; originally announced March 2021.

Comments: 15 pages, 6 figures, configurations in ancillary files, v2 corrects a typo

Report number: FERMILAB-FN-1125-ND

Neural 3D Video Synthesis from Multi-view Video

Authors: Tianye Li, Mira Slavcheva, Michael Zollhoefer, Simon Green, Christoph Lassner, Changil Kim, Tanner Schmidt, Steven Lovegrove, Michael Goesele, Richard Newcombe, Zhaoyang Lv

Abstract: We propose a novel approach for 3D video synthesis that is able to represent multi-view video recordings of a dynamic real-world scene in a compact, yet expressive representation that enables high-quality view synthesis and motion interpolation. Our approach takes the high quality and compactness of static neural radiance fields in a new direction: to a model-free, dynamic setting. At the core of… ▽ More We propose a novel approach for 3D video synthesis that is able to represent multi-view video recordings of a dynamic real-world scene in a compact, yet expressive representation that enables high-quality view synthesis and motion interpolation. Our approach takes the high quality and compactness of static neural radiance fields in a new direction: to a model-free, dynamic setting. At the core of our approach is a novel time-conditioned neural radiance field that represents scene dynamics using a set of compact latent codes. We are able to significantly boost the training speed and perceptual quality of the generated imagery by a novel hierarchical training scheme in combination with ray importance sampling. Our learned representation is highly compact and able to represent a 10 second 30 FPS multiview video recording by 18 cameras with a model size of only 28MB. We demonstrate that our method can render high-fidelity wide-angle novel views at over 1K resolution, even for complex and dynamic scenes. We perform an extensive qualitative and quantitative evaluation that shows that our approach outperforms the state of the art. Project website: https://neural-3d-video.github.io/. △ Less

Submitted 2 May, 2022; v1 submitted 3 March, 2021; originally announced March 2021.

Comments: Accepted as an oral presentation for CVPR 2022. Project website: https://neural-3d-video.github.io/

The tidal evolution of dark matter substructure -- II. The impact of artificial disruption on subhalo mass functions and radial profiles

Authors: Sheridan B. Green, Frank C. van den Bosch, Fangzhou Jiang

Abstract: Several recent studies have indicated that artificial subhalo disruption (the spontaneous, non-physical disintegration of a subhalo) remains prevalent in state-of-the-art dark matter-only cosmological simulations. In order to quantify the impact of disruption on the inferred subhalo demographics, we augment the semi-analytical SatGen dynamical subhalo evolution model with an improved treatment of… ▽ More Several recent studies have indicated that artificial subhalo disruption (the spontaneous, non-physical disintegration of a subhalo) remains prevalent in state-of-the-art dark matter-only cosmological simulations. In order to quantify the impact of disruption on the inferred subhalo demographics, we augment the semi-analytical SatGen dynamical subhalo evolution model with an improved treatment of tidal stripping that is calibrated using the DASH database of idealized high-resolution simulations of subhalo evolution, which are free from artificial disruption. We also develop a model of artificial disruption that reproduces the statistical properties of disruption in the Bolshoi simulation. Using this framework, we predict subhalo mass functions (SHMFs), number density profiles, and substructure mass fractions and study how these quantities are impacted by artificial disruption and mass resolution limits. We find that artificial disruption affects these quantities at the $10-20\%$ level, ameliorating previous concerns that it may suppress the SHMF by as much as a factor of two. We demonstrate that semi-analytical substructure modeling must include orbit integration in order to properly account for splashback haloes, which make up roughly half of the subhalo population. We show that the resolution limit of $N$-body simulations, rather than artificial disruption, is the primary cause of the radial bias in subhalo number density found in dark matter-only simulations. Hence, we conclude that the mass resolution remains the primary limitation of using such simulations to study subhaloes. Our model provides a fast, flexible, and accurate alternative to studying substructure statistics in the absence of both numerical resolution limits and artificial disruption. △ Less

Submitted 26 September, 2021; v1 submitted 1 March, 2021; originally announced March 2021.

Comments: 18 pages, 11 figures. Published in MNRAS. This arXiv version includes the changes described in the erratum

Journal ref: MNRAS 503, 4075 (2021, original manuscript) and MNRAS 508, 2944 (2021, erratum)

Seeing the Bigger Picture: The Rosetta Mission Amateur Observing Campaign and Lessons for the Future

Authors: Helen Usher, Colin Snodgrass, Simon F. Green, Andrew Norton, Paul Roche

Abstract: Amateur astronomers can make useful contributions to the study of comets. They add temporal coverage and multi-scale observations which can aid the study of fast-changing, and large-scale comet features. We document and review the amateur observing campaign set up to complement the Rosetta space mission, including the data submitted to date, and consider the campaign's effectiveness in the light o… ▽ More Amateur astronomers can make useful contributions to the study of comets. They add temporal coverage and multi-scale observations which can aid the study of fast-changing, and large-scale comet features. We document and review the amateur observing campaign set up to complement the Rosetta space mission, including the data submitted to date, and consider the campaign's effectiveness in the light of experience from previous comet amateur campaigns. We report the results of surveys of campaign participants, the amateur astronomy community, and schools who participated in a comet 46P observing campaign. We draw lessons for future campaigns which include the need for: clarity of objectives; recognising the wider impact campaigns can have on increasing science capital; clear, consistent, timely and tailored guidance; easy upload procedures with in-built quality control; and, regular communication, feedback and recognition. △ Less

Submitted 13 November, 2020; originally announced November 2020.

Comments: 32 pages, 8 figures To be published in The AAS Planetary Science Journal A'Hearn Symposium Focus Issue