Showing 1–50 of 103 results for author: Martinez, G

General Relativistic effects and the NIR variability of Sgr A* II: A systematic approach to temporal asymmetry

Authors: Sebastiano D. von Fellenberg, Gunther Witzel, Michi Bauboeck, Hui-Hsuan Chung, Nicola Marchili, Greg Martinez, Matteo Sadun-Bordoni, Guillaume Bourdarot, Tuan Do, Antonia Drescher, Giovanni Fazio, Frank Eisenhauer, Reinhard Genzel, Stefan Gillessen, Joseph L. Hora, Felix Mang, Thomas Ott, Howard A. Smith, Eduardo Ros, Diogo C. Ribeiro, Felix Widmann, S. P. Willner, J. Anton Zensus

Abstract: A systematic study, based on the third-moment structure function, of Sgr A*'s variability finds an exponential rise time $τ_{1,\rm{obs}}=14.8^{+0.4}_{-1.5}~\mathrm{minutes}$ and decay time $τ_{2,\rm{obs}}=13.1^{+1.3}_{-1.4}~\mathrm{minutes}$. This symmetry of the flux-density variability is consistent with earlier work, and we interpret it as caused by the dominance of Doppler boosting, as opposed… ▽ More A systematic study, based on the third-moment structure function, of Sgr A*'s variability finds an exponential rise time $τ_{1,\rm{obs}}=14.8^{+0.4}_{-1.5}~\mathrm{minutes}$ and decay time $τ_{2,\rm{obs}}=13.1^{+1.3}_{-1.4}~\mathrm{minutes}$. This symmetry of the flux-density variability is consistent with earlier work, and we interpret it as caused by the dominance of Doppler boosting, as opposed to gravitational lensing, in Sgr~A*'s light curve. A relativistic, semi-physical model of Sgr~A* confirms an inclination angle $i<45$ degrees. The model also shows that the emission of the intrinsic radiative process can have some asymmetry even though the observed emission does not. The third-moment structure function, which is a measure of the skewness of the light-curve increments, may be a useful summary statistic in other contexts of astronomy because it senses only temporal asymmetry, i.e., it averages to zero for any temporally symmetric signal. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: Accepted for publication in A&A letters

A detailed study of the very-high-energy Crab pulsar emission with the LST-1

Authors: CTA-LST Project, :, K. Abe, S. Abe, A. Abhishek, F. Acero, A. Aguasca-Cabot, I. Agudo, N. Alvarez Crespo, L. A. Antonelli, C. Aramo, A. Arbet-Engels, C. Arcaro, M. Artero, K. Asano, P. Aubert, A. Baktash, A. Bamba, A. Baquero Larriva, L. Baroncelli, U. Barres de Almeida, J. A. Barrio, I. Batkovic, J. Baxter, J. Becerra González , et al. (272 additional authors not shown)

Abstract: Context: There are currently three pulsars firmly detected by imaging atmospheric Cherenkov telescopes (IACTs), two of them reaching TeV energies, challenging models of very-high-energy (VHE) emission in pulsars. More precise observations are needed to better characterize pulsar emission at these energies. The LST-1 is the prototype of the Large-Sized Telescope, that will be part of the Cherenkov… ▽ More Context: There are currently three pulsars firmly detected by imaging atmospheric Cherenkov telescopes (IACTs), two of them reaching TeV energies, challenging models of very-high-energy (VHE) emission in pulsars. More precise observations are needed to better characterize pulsar emission at these energies. The LST-1 is the prototype of the Large-Sized Telescope, that will be part of the Cherenkov Telescope Array Observatory (CTAO). Its improved performance over previous IACTs makes it well suited for studying pulsars. Aims: To study the Crab pulsar emission with the LST-1, improving and complementing the results from other telescopes. These observations can also be used to characterize the potential of the LST-1 to study other pulsars and detect new ones. Methods: We analyzed a total of $\sim$103 hours of gamma-ray observations of the Crab pulsar conducted with the LST-1 in the period from September 2020 to January 2023. The observations were carried out at zenith angles less than 50 degrees. A new analysis of the Fermi-LAT data was also performed, including $\sim$14 years of observations. Results: The Crab pulsar phaseogram, long-term light-curve, and phase-resolved spectra are reconstructed with the LST-1 from 20 GeV to 450 GeV for P1 and up to 700 GeV for P2. The pulsed emission is detected with a significance of 15.2$σ$. The two characteristic emission peaks of the Crab pulsar are clearly detected (>10$σ$), as well as the so-called bridge emission (5.7$σ$). We find that both peaks are well described by power laws, with spectral indices of $\sim$3.44 and $\sim$3.03 respectively. The joint analysis of Fermi-LAT and LST-1 data shows a good agreement between both instruments in the overlapping energy range. The detailed results obtained in the first observations of the Crab pulsar with LST-1 show the potential that CTAO will have to study this type of sources. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: Accepted by A&A

Orbital and Atmospheric Characterization of the 1RXS J034231.8+121622 System Using High-Resolution Spectroscopy Confirms That The Companion is a Low-Mass Star

Authors: Clarissa R. Do Ó, Ben Sappey, Quinn M. Konopacky, Jean-Baptiste Ruffio, Kelly K. O'Neil, Tuan Do, Gregory Martinez, Travis S. Barman, Jayke S. Nguyen, Jerry W. Xuan, Christopher A. Theissen, Sarah Blunt, William Thompson, Chih-Chun Hsu, Ashley Baker, Randall Bartos, Geoffrey A. Blake, Benjamin Calvin, Sylvain Cetre, Jacques-Robert Delorme, Greg Doppmann, Daniel Echeverri, Luke Finnerty, Michael P. Fitzgerald, Julie Inglis , et al. (11 additional authors not shown)

Abstract: The 1RXS J034231.8+121622 system consists of an M dwarf primary and a directly imaged low-mass stellar companion. We use high resolution spectroscopic data from Keck/KPIC to estimate the objects' atmospheric parameters and radial velocities (RVs). Using PHOENIX stellar models, we find that the primary has a temperature of 3460 $\pm$ 50 K a metallicity of 0.16 $\pm$ 0.04, while the secondary has a… ▽ More The 1RXS J034231.8+121622 system consists of an M dwarf primary and a directly imaged low-mass stellar companion. We use high resolution spectroscopic data from Keck/KPIC to estimate the objects' atmospheric parameters and radial velocities (RVs). Using PHOENIX stellar models, we find that the primary has a temperature of 3460 $\pm$ 50 K a metallicity of 0.16 $\pm$ 0.04, while the secondary has a temperature of 2510 $\pm$ 50 K and a metallicity of $0.13\substack{+0.12 \\ -0.11}$. Recent work suggests this system is associated with the Hyades, placing it an older age than previous estimates. Both metallicities agree with current $[Fe/H]$ Hyades measurements (0.11 -- 0.21). Using stellar evolutionary models, we obtain significantly higher masses for the objects, of 0.30 $\pm$ 0.15 $M_\odot$ and 0.08 $\pm$ 0.01 $M_\odot$ (84 $\pm$ 11 $M_{Jup}$) respectively. Using the RVs and a new astrometry point from Keck/NIRC2, we find that the system is likely an edge-on, moderately eccentric ($0.41\substack{+0.27 \\ -0.08}$) configuration. We also estimate the C/O ratio of both objects using custom grid models, obtaining 0.42 $\pm$ 0.10 (primary) and 0.55 $\pm$ 0.10 (companion). From these results, we confirm that this system most likely went through a binary star formation process in the Hyades. The significant changes in this system's parameters since its discovery highlight the importance of high resolution spectroscopy for both orbital and atmospheric characterization of directly imaged companions. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: 30 pages, 18 figures, accepted for publication in AJ

Dark Matter Line Searches with the Cherenkov Telescope Array

Authors: S. Abe, J. Abhir, A. Abhishek, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Alfaro, N. Alvarez-Crespo, R. Alves Batista, J. -P. Amans, E. Amato, G. Ambrosi, L. Angel, C. Aramo, C. Arcaro, T. T. H. Arnesen, L. Arrabito, K. Asano, Y. Ascasibar, J. Aschersleben, H. Ashkar , et al. (540 additional authors not shown)

Abstract: Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of sele… ▽ More Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of selected dwarf spheroidal galaxies. We find that current limits and detection prospects for dark matter masses above 300 GeV will be significantly improved, by up to an order of magnitude in the multi-TeV range. This demonstrates that CTA will set a new standard for gamma-ray astronomy also in this respect, as the world's largest and most sensitive high-energy gamma-ray observatory, in particular due to its exquisite energy resolution at TeV energies and the adopted observational strategy focussing on regions with large dark matter densities. Throughout our analysis, we use up-to-date instrument response functions, and we thoroughly model the effect of instrumental systematic uncertainties in our statistical treatment. We further present results for other potential signatures with sharp spectral features, e.g.~box-shaped spectra, that would likewise very clearly point to a particle dark matter origin. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: 43 pages JCAP style (excluding author list and references), 19 figures

Mars 2020 Perseverance rover studies of the Martian atmosphere over Jezero from pressure measurements

Authors: A. Sánchez-Lavega, T. del Rio-Gaztelurrutia, R. Hueso, M. de la Torre Juárez, G. M. Martínez, A. -M. Harri, M. Genzer, M. Hieta, J. Polkko, J. A. Rodríguez-Manfredi, M. T. Lemmon, J. Pla-García, D. Toledo, A. Vicente-Retortillo, Daniel Viúdez-Moreiras, A. Munguira, L. K. Tamppari, C. Newman, J. Gómez-Elvira, S. Guzewich, T. Bertrand, V. Apéstigue, I. Arruego, M. Wolff, D. Banfield , et al. (2 additional authors not shown)

Abstract: The pressure sensors on Mars rover Perseverance measure the pressure field in the Jezero crater on regular hourly basis starting in sol 15 after landing. The present study extends up to sol 460 encompassing the range of solar longitudes from Ls 13° - 241° (Martian Year (MY) 36). The data show the changing daily pressure cycle, the sol-to-sol seasonal evolution of the mean pressure field driven by… ▽ More The pressure sensors on Mars rover Perseverance measure the pressure field in the Jezero crater on regular hourly basis starting in sol 15 after landing. The present study extends up to sol 460 encompassing the range of solar longitudes from Ls 13° - 241° (Martian Year (MY) 36). The data show the changing daily pressure cycle, the sol-to-sol seasonal evolution of the mean pressure field driven by the CO2 sublimation and deposition cycle at the poles, the characterization of up to six components of the atmospheric tides and their relationship to dust content in the atmosphere. They also show the presence of wave disturbances with periods 2-5 sols, exploring their baroclinic nature, short period oscillations (mainly at night-time) in the range 8-24 minutes that we interpret as internal gravity waves, transient pressure drops with duration 1-150 s produced by vortices, and rapid turbulent fluctuations. We also analyze the effects on pressure measurements produced by a regional dust storm over Jezero at Ls 155°. △ Less

Submitted 23 January, 2024; originally announced January 2024.

Journal ref: Journal of Geophysical Research: Planets, 128, e2022JE007480 (2023)

Chasing Gravitational Waves with the Cherenkov Telescope Array

Authors: Jarred Gershon Green, Alessandro Carosi, Lara Nava, Barbara Patricelli, Fabian Schüssler, Monica Seglar-Arroyo, Cta Consortium, :, Kazuki Abe, Shotaro Abe, Atreya Acharyya, Remi Adam, Arnau Aguasca-Cabot, Ivan Agudo, Jorge Alfaro, Nuria Alvarez-Crespo, Rafael Alves Batista, Jean-Philippe Amans, Elena Amato, Filippo Ambrosino, Ekrem Oguzhan Angüner, Lucio Angelo Antonelli, Carla Aramo, Cornelia Arcaro, Luisa Arrabito , et al. (545 additional authors not shown)

Abstract: The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very… ▽ More The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very-high-energy (VHE, > 100GeV) photons which have yet to be detected in coincidence with a gravitational wave signal. The Cherenkov Telescope Array (CTA) is a next-generation VHE observatory which aims to be indispensable in this search, with an unparalleled sensitivity and ability to slew anywhere on the sky within a few tens of seconds. New observing modes and follow-up strategies are being developed for CTA to rapidly cover localization areas of gravitational wave events that are typically larger than the CTA field of view. This work will evaluate and provide estimations on the expected number of of gravitational wave events that will be observable with CTA, considering both on- and off-axis emission. In addition, we will present and discuss the prospects of potential follow-up strategies with CTA. △ Less

Submitted 5 February, 2024; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: Presented at the 38th International Cosmic Ray Conference (ICRC 2023), 2023 (arXiv:2309.08219)

Report number: CTA-ICRC/2023/30

Performance of the joint LST-1 and MAGIC observations evaluated with Crab Nebula data

Authors: H. Abe, K. Abe, S. Abe, V. A. Acciari, A. Aguasca-Cabot, I. Agudo, N. Alvarez Crespo, T. Aniello, S. Ansoldi, L. A. Antonelli, C. Aramo, A. Arbet-Engels, C. Arcaro, M. Artero, K. Asano, P. Aubert, D. Baack, A. Babić, A. Baktash, A. Bamba, A. Baquero Larriva, L. Baroncelli, U. Barres de Almeida, J. A. Barrio, I. Batković , et al. (344 additional authors not shown)

Abstract: Aims. LST-1, the prototype of the Large-Sized Telescope for the upcoming Cherenkov Telescope Array Observatory, is concluding its commissioning in Observatorio del Roque de los Muchachos on the island of La Palma. The proximity of LST-1 (Large-Sized Telescope 1) to the two MAGIC (Major Atmospheric Gamma Imaging Cherenkov) telescopes permits observations of the same gamma-ray events with both syste… ▽ More Aims. LST-1, the prototype of the Large-Sized Telescope for the upcoming Cherenkov Telescope Array Observatory, is concluding its commissioning in Observatorio del Roque de los Muchachos on the island of La Palma. The proximity of LST-1 (Large-Sized Telescope 1) to the two MAGIC (Major Atmospheric Gamma Imaging Cherenkov) telescopes permits observations of the same gamma-ray events with both systems. Methods. We describe the joint LST-1+MAGIC analysis pipeline and use simultaneous Crab Nebula observations and Monte Carlo simulations to assess the performance of the three-telescope system. The addition of the LST-1 telescope allows the recovery of events in which one of the MAGIC images is too dim to survive analysis quality cuts. Results. Thanks to the resulting increase in the collection area and stronger background rejection, we find a significant improvement in sensitivity, allowing the detection of 30% weaker fluxes in the energy range between 200 GeV and 3 TeV. The spectrum of the Crab Nebula, reconstructed in the energy range ~60 GeV to ~10 TeV, is in agreement with previous measurements. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 680, A66 (2023)

Prospects for $γ$-ray observations of the Perseus galaxy cluster with the Cherenkov Telescope Array

Authors: The Cherenkov Telescope Array Consortium, :, K. Abe, S. Abe, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Alfaro, N. Alvarez-Crespo, R. Alves Batista, J. -P. Amans, E. Amato, E. O. Angüner, L. A. Antonelli, C. Aramo, M. Araya, C. Arcaro, L. Arrabito, K. Asano, Y. Ascasíbar, J. Aschersleben , et al. (542 additional authors not shown)

Abstract: Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster med… ▽ More Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster medium. We estimate the sensitivity of the Cherenkov Telescope Array (CTA) to detect diffuse gamma-ray emission from the Perseus galaxy cluster. We perform a detailed spatial and spectral modelling of the expected signal for the DM and the CRp components. For each, we compute the expected CTA sensitivity. The observing strategy of Perseus is also discussed. In the absence of a diffuse signal (non-detection), CTA should constrain the CRp to thermal energy ratio within the radius $R_{500}$ down to about $X_{500}<3\times 10^{-3}$, for a spatial CRp distribution that follows the thermal gas and a CRp spectral index $α_{\rm CRp}=2.3$. Under the optimistic assumption of a pure hadronic origin of the Perseus radio mini-halo and depending on the assumed magnetic field profile, CTA should measure $α_{\rm CRp}$ down to about $Δα_{\rm CRp}\simeq 0.1$ and the CRp spatial distribution with 10% precision. Regarding DM, CTA should improve the current ground-based gamma-ray DM limits from clusters observations on the velocity-averaged annihilation cross-section by a factor of up to $\sim 5$, depending on the modelling of DM halo substructure. In the case of decay of DM particles, CTA will explore a new region of the parameter space, reaching models with $τ_χ>10^{27}$s for DM masses above 1 TeV. These constraints will provide unprecedented sensitivity to the physics of both CRp acceleration and transport at cluster scale and to TeV DM particle models, especially in the decay scenario. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Comments: 93 pages (including author list, appendix and references), 143 figures. Submitted to JCAP

Expanding Mars Climate Modeling: Interpretable Machine Learning for Modeling MSL Relative Humidity

Authors: Nour Abdelmoneim, Dattaraj B. Dhuri, Dimitra Atri, Germán Martínez

Abstract: For the past several decades, numerous attempts have been made to model the climate of Mars with extensive studies focusing on the planet's dynamics and the understanding of its climate. While physical modeling and data assimilation approaches have made significant progress, uncertainties persist in comprehensively capturing and modeling the complexities of Martian climate. In this work, we propos… ▽ More For the past several decades, numerous attempts have been made to model the climate of Mars with extensive studies focusing on the planet's dynamics and the understanding of its climate. While physical modeling and data assimilation approaches have made significant progress, uncertainties persist in comprehensively capturing and modeling the complexities of Martian climate. In this work, we propose a novel approach to Martian climate modeling by leveraging machine learning techniques that have shown remarkable success in Earth climate modeling. Our study presents a deep neural network designed to accurately model relative humidity in Gale Crater, as measured by NASA's Mars Science Laboratory ``Curiosity'' rover. By utilizing simulated meteorological variables produced by the Mars Planetary Climate Model, a robust Global Circulation Model, our model accurately predicts relative humidity with a mean error of 3\% and an $R^2$ score of 0.92. Furthermore, we present an approach to predict quantile ranges of relative humidity, catering to applications that require a range of values. To address the challenge of interpretability associated with machine learning models, we utilize an interpretable model architecture and conduct an in-depth analysis of its internal mechanisms and decision making processes. We find that our neural network can effectively model relative humidity at Gale crater using a few meteorological variables, with the monthly mean surface H$_2$O layer, planetary boundary layer height, convective wind speed, and solar zenith angle being the primary contributors to the model predictions. In addition to providing a fast and efficient method to modeling climate variables on Mars, this modeling approach can also be used to expand on current datasets by filling spatial and temporal gaps in observations. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Comments: 13 pages, 7 figures

Near-Infrared Flux Distribution of Sgr A* from 2005-2022: Evidence for an Enhanced Accretion Episode in 2019

Authors: Grant C. Weldon, Tuan Do, Gunther Witzel, Andrea M. Ghez, Abhimat K. Gautam, Eric E. Becklin, Mark R. Morris, Gregory D. Martinez, Shoko Sakai, Jessica R. Lu, Keith Matthews, Matthew W. Hosek Jr., Zoë Haggard

Abstract: Sgr A* is the variable electromagnetic source associated with accretion onto the Galactic center supermassive black hole. While the near-infrared (NIR) variability of Sgr A* was shown to be consistent over two decades, unprecedented activity in 2019 challenges existing statistical models. We investigate the origin of this activity by re-calibrating and re-analyzing all of our Keck Observatory Sgr… ▽ More Sgr A* is the variable electromagnetic source associated with accretion onto the Galactic center supermassive black hole. While the near-infrared (NIR) variability of Sgr A* was shown to be consistent over two decades, unprecedented activity in 2019 challenges existing statistical models. We investigate the origin of this activity by re-calibrating and re-analyzing all of our Keck Observatory Sgr A* imaging observations from 2005-2022. We present light curves from 69 observation epochs using the NIRC2 imager at 2.12 $μ$m with laser guide star adaptive optics. These observations reveal that the mean luminosity of Sgr A* increased by a factor of $\sim$3 in 2019, and the 2019 light curves had higher variance than in all time periods we examined. We find that the 2020-2022 flux distribution is statistically consistent with the historical sample and model predictions, but with fewer bright measurements above 0.6 mJy at the $\sim$2$σ$ level. Since 2019, we have observed a maximum $K_s$ (2.2 $μ$m) flux of 0.9 mJy, compared to the highest pre-2019 flux of 2.0 mJy and highest 2019 flux of 5.6 mJy. Our results suggest that the 2019 activity was caused by a temporary accretion increase onto Sgr A*, possibly due to delayed accretion of tidally-stripped gas from the gaseous object G2 in 2014. We also examine faint Sgr A* fluxes over a long time baseline to search for a quasi-steady quiescent state. We find that Sgr A* displays flux variations over a factor of $\sim$500, with no evidence for a quiescent state in the NIR. △ Less

Submitted 18 August, 2023; originally announced August 2023.

Comments: Accepted to ApJ Letters, 26 pages, 15 figures

Observations of the Crab Nebula and Pulsar with the Large-Sized Telescope Prototype of the Cherenkov Telescope Array

Authors: CTA-LST Project, :, H. Abe, K. Abe, S. Abe, A. Aguasca-Cabot, I. Agudo, N. Alvarez Crespo, L. A. Antonelli, C. Aramo, A. Arbet-Engels, C. Arcaro, M. Artero, K. Asano, P. Aubert, A. Baktash, A. Bamba, A. Baquero Larriva, L. Baroncelli, U. Barres de Almeida, J. A. Barrio, I. Batkovic, J. Baxter, J. Becerra González, E. Bernardini , et al. (467 additional authors not shown)

Abstract: CTA (Cherenkov Telescope Array) is the next generation ground-based observatory for gamma-ray astronomy at very-high energies. The Large-Sized Telescope prototype (LST-1) is located at the Northern site of CTA, on the Canary Island of La Palma. LSTs are designed to provide optimal performance in the lowest part of the energy range covered by CTA, down to $\simeq 20$ GeV. LST-1 started performing a… ▽ More CTA (Cherenkov Telescope Array) is the next generation ground-based observatory for gamma-ray astronomy at very-high energies. The Large-Sized Telescope prototype (LST-1) is located at the Northern site of CTA, on the Canary Island of La Palma. LSTs are designed to provide optimal performance in the lowest part of the energy range covered by CTA, down to $\simeq 20$ GeV. LST-1 started performing astronomical observations in November 2019, during its commissioning phase, and it has been taking data since then. We present the first LST-1 observations of the Crab Nebula, the standard candle of very-high energy gamma-ray astronomy, and use them, together with simulations, to assess the basic performance parameters of the telescope. The data sample consists of around 36 hours of observations at low zenith angles collected between November 2020 and March 2022. LST-1 has reached the expected performance during its commissioning period - only a minor adjustment of the preexisting simulations was needed to match the telescope behavior. The energy threshold at trigger level is estimated to be around 20 GeV, rising to $\simeq 30$ GeV after data analysis. Performance parameters depend strongly on energy, and on the strength of the gamma-ray selection cuts in the analysis: angular resolution ranges from 0.12 to 0.40 degrees, and energy resolution from 15 to 50%. Flux sensitivity is around 1.1% of the Crab Nebula flux above 250 GeV for a 50-h observation (12% for 30 minutes). The spectral energy distribution (in the 0.03 - 30 TeV range) and the light curve obtained for the Crab Nebula agree with previous measurements, considering statistical and systematic uncertainties. A clear periodic signal is also detected from the pulsar at the center of the Nebula. △ Less

Submitted 19 July, 2023; v1 submitted 22 June, 2023; originally announced June 2023.

Comments: Accepted in ApJ. v3: updated author list and acknowledgements, fixed typos and other minor issues

Discovery and Timing of Millisecond Pulsars with the Arecibo 327 MHz Drift-Scan Survey

Authors: Evan F. Lewis, Timothy E. E. Olszanski, Julia S. Deneva, Paulo C. C. Freire, Maura A. McLaughlin, Kevin Stovall, Manjari Bagchi, Jose G. Martinez, Benetge B. P. Perera

Abstract: We present the discovery and timing solutions of four millisecond pulsars (MSPs) discovered in the Arecibo 327 MHz Drift-Scan Pulsar Survey. Three of these pulsars are in binary systems, consisting of a redback (PSR J2055+1545), a black widow (PSR J1630+3550), and a neutron star-white dwarf binary (PSR J2116+1345). The fourth MSP, PSR J2212+2450, is isolated. We present the multiyear timing soluti… ▽ More We present the discovery and timing solutions of four millisecond pulsars (MSPs) discovered in the Arecibo 327 MHz Drift-Scan Pulsar Survey. Three of these pulsars are in binary systems, consisting of a redback (PSR J2055+1545), a black widow (PSR J1630+3550), and a neutron star-white dwarf binary (PSR J2116+1345). The fourth MSP, PSR J2212+2450, is isolated. We present the multiyear timing solutions as well as polarization properties across a range of radio frequencies for each pulsar. We perform a multiwavelength search for emission from these systems and find an optical counterpart for PSR J2055+1545 in Gaia DR3, as well as a gamma-ray counterpart for PSR J2116+1345 with the Fermi-LAT telescope. Despite the close colocation of PSR J2055+1545 with a Fermi source, we are unable to detect gamma-ray pulsations, likely due to the large orbital variability of the system. This work presents the first two binaries found by this survey with orbital periods shorter than a day; we expect to find more in the 40% of the survey data that have yet to be searched. △ Less

Submitted 19 October, 2023; v1 submitted 16 June, 2023; originally announced June 2023.

Comments: 22 pages, 11 figures, 4 tables. Published in ApJ

The Orbital Eccentricities of Directly Imaged Companions Using Observable-Based Priors: Implications for Population-level Distributions

Authors: Clarissa R. Do Ó, Kelly K. O'Neil, Quinn M. Konopacky, Tuan Do, Gregory D. Martinez, Jean-Baptiste Ruffio, Andrea M. Ghez

Abstract: The eccentricity of a sub-stellar companion is an important tracer of its formation history. Directly imaged companions often present poorly constrained eccentricities. A recently developed prior framework for orbit fitting called ''observable-based priors'' has the advantage of improving biases in derived orbit parameters for objects with minimal phase coverage, which is the case for the majority… ▽ More The eccentricity of a sub-stellar companion is an important tracer of its formation history. Directly imaged companions often present poorly constrained eccentricities. A recently developed prior framework for orbit fitting called ''observable-based priors'' has the advantage of improving biases in derived orbit parameters for objects with minimal phase coverage, which is the case for the majority of directly imaged companions. We use observable-based priors to fit the orbits of 21 exoplanets and brown dwarfs in an effort to obtain the eccentricity distributions with minimized biases. We present the objects' individual posteriors compared to their previously derived distributions, showing in many cases a shift toward lower eccentricities. We analyze the companions' eccentricity distribution at a population level, and compare this to the distributions obtained with the traditional uniform priors. We fit a Beta distribution to our posteriors using observable-based priors, obtaining shape parameters $α= 1.09^{+0.30}_{-0.22}$ and $β= 1.42^{+0.33}_{-0.25}$. This represents an approximately flat distribution of eccentricities. The derived $α$ and $β$ parameters are consistent with the values obtained using uniform priors, though uniform priors lead to a tail at high eccentricities. We find that separating the population into high and low mass companions yields different distributions depending on the classification of intermediate mass objects. We also determine via simulation that the minimal orbit coverage needed to give meaningful posteriors under the assumptions made for directly imaged planets is $\approx$ 15% of the inferred period of the orbit. △ Less

Submitted 18 July, 2023; v1 submitted 6 June, 2023; originally announced June 2023.

Comments: Published in The Astronomical Journal, Volume 166, Issue 2, id.48, 22 pp. Update includes updated references and ORCID profiles

Sensitivity of the Cherenkov Telescope Array to spectral signatures of hadronic PeVatrons with application to Galactic Supernova Remnants

Authors: The Cherenkov Telescope Array Consortium, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Aloisio, N. Álvarez Crespo, R. Alves Batista, L. Amati, E. Amato, G. Ambrosi, E. O. Angüner, C. Aramo, C. Arcaro, T. Armstrong, K. Asano, Y. Ascasibar, J. Aschersleben, M. Backes, A. Baktash, C. Balazs, M. Balbo , et al. (334 additional authors not shown)

Abstract: The local Cosmic Ray (CR) energy spectrum exhibits a spectral softening at energies around 3~PeV. Sources which are capable of accelerating hadrons to such energies are called hadronic PeVatrons. However, hadronic PeVatrons have not yet been firmly identified within the Galaxy. Several source classes, including Galactic Supernova Remnants (SNRs), have been proposed as PeVatron candidates. The pote… ▽ More The local Cosmic Ray (CR) energy spectrum exhibits a spectral softening at energies around 3~PeV. Sources which are capable of accelerating hadrons to such energies are called hadronic PeVatrons. However, hadronic PeVatrons have not yet been firmly identified within the Galaxy. Several source classes, including Galactic Supernova Remnants (SNRs), have been proposed as PeVatron candidates. The potential to search for hadronic PeVatrons with the Cherenkov Telescope Array (CTA) is assessed. The focus is on the usage of very high energy $γ$-ray spectral signatures for the identification of PeVatrons. Assuming that SNRs can accelerate CRs up to knee energies, the number of Galactic SNRs which can be identified as PeVatrons with CTA is estimated within a model for the evolution of SNRs. Additionally, the potential of a follow-up observation strategy under moonlight conditions for PeVatron searches is investigated. Statistical methods for the identification of PeVatrons are introduced, and realistic Monte--Carlo simulations of the response of the CTA observatory to the emission spectra from hadronic PeVatrons are performed. Based on simulations of a simplified model for the evolution for SNRs, the detection of a $γ$-ray signal from in average 9 Galactic PeVatron SNRs is expected to result from the scan of the Galactic plane with CTA after 10 hours of exposure. CTA is also shown to have excellent potential to confirm these sources as PeVatrons in deep observations with $\mathcal{O}(100)$ hours of exposure per source. △ Less

Submitted 27 March, 2023; originally announced March 2023.

Comments: 34 pages, 16 figures, Accepted for publication in Astroparticle Physics

The Swansong of the Galactic Center Source X7: An Extreme Example of Tidal Evolution near the Supermassive Black Hole

Authors: Anna Ciurlo, Randall D. Campbell, Mark R. Morris, Tuan Do, Andrea M. Ghez, Eric E. Becklin, Rory O. Bentley, Devin S. Chu, Abhimat K. Gautam, Yash A. Gursahani, Aurelien Hees, Kelly Kosmo O'Neil, Jessica R. Lu, Gregory D. Martinez, Smadar Naoz, Shoko Sakai, Rainer Schoedel

Abstract: We present two decades of new high-angular-resolution near-infrared data from the W. M. Keck Observatory that reveal extreme evolution in X7, an elongated dust and gas feature, presently located half an arcsecond from the Galactic Center supermassive black hole. With both spectro-imaging observations of Br-γ line-emission and Lp (3.8 μm) imaging data, we provide the first estimate of its orbital p… ▽ More We present two decades of new high-angular-resolution near-infrared data from the W. M. Keck Observatory that reveal extreme evolution in X7, an elongated dust and gas feature, presently located half an arcsecond from the Galactic Center supermassive black hole. With both spectro-imaging observations of Br-γ line-emission and Lp (3.8 μm) imaging data, we provide the first estimate of its orbital parameters and quantitative characterization of the evolution of its morphology and mass. We find that the leading edge of X7 appears to be on a mildly eccentric (e~0.3), relatively short-period (170 years) orbit and is headed towards periapse passage, estimated to occur in ~2036. Furthermore, our kinematic measurements rule out the earlier suggestion that X7 is associated with the stellar source S0-73 or with any other point source that has overlapped with X7 during our monitoring period. Over the course of our observations, X7 has (1) become more elongated, with a current length-to-width ratio of 9, (2) maintained a very consistent long-axis orientation (position angle of 50 deg), (3) inverted its radial velocity differential from tip to tail from -50 to +80 km/sec, and (4) sustained its total brightness (12.8 Lp magnitudes at the leading edge) and color temperature (425 K), which suggest a constant mass of ~50 MEarth. We present a simple model showing that these results are compatible with the expected effect of tidal forces exerted on it by the central black hole and we propose that X7 is the gas and dust recently ejected from a grazing collision in a binary system. △ Less

Submitted 16 January, 2023; originally announced January 2023.

Comments: accepted by ApJ, 20 pages, 16 figures

Multi-wavelength study of the galactic PeVatron candidate LHAASO J2108+5157

Authors: S. Abe, A. Aguasca-Cabot, I. Agudo, N. Alvarez Crespo, L. A. Antonelli, C. Aramo, A. Arbet-Engels, M. Artero, K. Asano, P. Aubert, A. Baktash, A. Bamba, A. Baquero Larriva, L. Baroncelli, U. Barres de Almeida, J. A. Barrio, I. Batkovic, J. Baxter, J. Becerra González, E. Bernardini, M. I. Bernardos, J. Bernete Medrano, A. Berti, P. Bhattacharjee, N. Biederbeck , et al. (245 additional authors not shown)

Abstract: LHAASO J2108+5157 is one of the few known unidentified Ultra-High-Energy (UHE) gamma-ray sources with no Very-High-Energy (VHE) counterpart, recently discovered by the LHAASO collaboration. We observed LHAASO J2108+5157 in the X-ray band with XMM-Newton in 2021 for a total of 3.8 hours and at TeV energies with the Large-Sized Telescope prototype (LST-1), yielding 49 hours of good quality data. In… ▽ More LHAASO J2108+5157 is one of the few known unidentified Ultra-High-Energy (UHE) gamma-ray sources with no Very-High-Energy (VHE) counterpart, recently discovered by the LHAASO collaboration. We observed LHAASO J2108+5157 in the X-ray band with XMM-Newton in 2021 for a total of 3.8 hours and at TeV energies with the Large-Sized Telescope prototype (LST-1), yielding 49 hours of good quality data. In addition, we analyzed 12 years of Fermi-LAT data, to better constrain emission of its High-Energy (HE) counterpart 4FGL J2108.0+5155. We found an excess (3.7 sigma) in the LST-1 data at energies E > 3 TeV. Further analysis in the whole LST-1 energy range assuming a point-like source, resulted in a hint (2.2 sigma) of hard emission which can be described with a single power law with photon index Gamma = 1.6 +- 0.2 between 0.3 - 100 TeV. We did not find any significant extended emission which could be related to a Supernova Remnant (SNR) or Pulsar Wind Nebula (PWN) in the XMM-Newton data, which puts strong constraints on possible synchrotron emission of relativistic electrons. The LST-1 and LHAASO observations can be explained as inverse Compton-dominated leptonic emission of relativistic electrons with a cutoff energy of $100^{+70}_{-30}$ TeV. The low magnetic field in the source imposed by the X-ray upper limits on synchrotron emission is compatible with a hypothesis of a PWN or a TeV halo. The lack of a pulsar in the neighborhood of the UHE source is a challenge to the PWN/TeV-halo scenario. The UHE gamma rays can also be explained as $π^0$ decay-dominated hadronic emission due to interaction of relativistic protons with one of the two known molecular clouds in the direction of the source. The hard spectrum in the LST-1 band is compatible with protons escaping a shock around a middle-aged SNR because of their high low-energy cut-off. △ Less

Submitted 16 March, 2023; v1 submitted 3 October, 2022; originally announced October 2022.

Comments: 26 pages, 9 figures, 5 tables; Astronomy and Astrophysics (accepted on 6 March 2023)

Journal ref: A&A 673, A75 (2023)

First measurements and upgrade plans of the MAGIC intensity interferometer

Authors: Juan Cortina, V. A. Acciari, A. Biland, E. Colombo, C. da Costa, C. Delgado, C. Diaz, M. Fiori, D. Fink, T. Hassan, I. Jimenez-Martinez, E. Lyard, M. Mariotti, G. Martinez, R. Mirzoyan, G. Naletto, M. Polo, N. Produit, J. J. Rodriguez, T. Schweizer, R. Walter, C. W. Wunderlich, L. Zampieri, the MAGIC, LST collaborations

Abstract: The two MAGIC 17-m diameter Imaging Atmospheric Cherenkov Telescopes have been equipped to work also as an intensity interferometer with a deadtime-free, 4-channel, GPU-based, real-time correlator. Operating with baselines between approx. 40 and 90 m the MAGIC interferometer is able to measure stellar diameters of 0.5-1 mas in the 400-440 nm wavelength range with a sensitivity roughly 10 times bet… ▽ More The two MAGIC 17-m diameter Imaging Atmospheric Cherenkov Telescopes have been equipped to work also as an intensity interferometer with a deadtime-free, 4-channel, GPU-based, real-time correlator. Operating with baselines between approx. 40 and 90 m the MAGIC interferometer is able to measure stellar diameters of 0.5-1 mas in the 400-440 nm wavelength range with a sensitivity roughly 10 times better than that achieved in the 1970s by the Narrabri Stellar Intensity Interferometer. Besides, active mirror control allows to split the primary mirrors into sub-mirrors. This allows to make simultaneous calibration measurements of the zero-baseline correlation or to simultaneously collect six baselines below 17 m with almost arbitrary orientation, corresponding to angular scales of approx. 1-50 mas. We plan to perform test observations adding the nearby Cherenkov Telescope Array (CTA) LST-1 23 m diameter telescope by next year. All three telescope pairs will be correlated simultaneously. Adding LST-1 is expected to increase the sensitivity by at least 1 mag and significantly improve the u-v plane coverage. If successful, the proposed correlator setup is scalable enough to be implemented to the full CTA arrays. △ Less

Submitted 29 September, 2022; originally announced September 2022.

Comments: 15 pages, 9 figures, 1 table, SPIE Astronomical Telescopes + Instrumentation 2022

Journal ref: Proc. SPIE 12183, Optical and Infrared Interferometry and Imaging VIII, 121830C (26 August 2022) 12183-11

Measuring the Orbits of the Arches and Quintuplet Clusters using HST and Gaia: Exploring Scenarios for Star Formation Near the Galactic Center

Authors: Matthew W. Hosek Jr., Tuan Do, Jessica R. Lu, Mark R. Morris, Andrea M. Ghez, Gregory D. Martinez, Jay Anderson

Abstract: We present new absolute proper motion measurements for the Arches and Quintuplet clusters, two young massive star clusters near the Galactic center. Using multi-epoch HST observations, we construct proper motion catalogs for the Arches ($\sim$35,000 stars) and Quintuplet ($\sim$40,000 stars) fields in ICRF coordinates established using stars in common with the Gaia EDR3 catalog. The bulk proper mo… ▽ More We present new absolute proper motion measurements for the Arches and Quintuplet clusters, two young massive star clusters near the Galactic center. Using multi-epoch HST observations, we construct proper motion catalogs for the Arches ($\sim$35,000 stars) and Quintuplet ($\sim$40,000 stars) fields in ICRF coordinates established using stars in common with the Gaia EDR3 catalog. The bulk proper motions of the clusters are measured to be ($μ_{α*}$, $μ_δ$) = (-0.80 $\pm$ 0.032, -1.89 $\pm$ 0.021) mas/yr for the Arches and ($μ_{α*}$, $μ_δ$) = (-0.96 $\pm$ 0.032, -2.29 $\pm$ 0.023) mas/yr for the Quintuplet, achieving $\sim$5x higher precision than past measurements. We place the first constraints on the properties of the cluster orbits that incorporate the uncertainty in their current line-of-sight distances. The clusters will not approach closer than $\sim$25 pc to SgrA*, making it unlikely that they will inspiral into the Nuclear Star Cluster within their lifetime. Further, the cluster orbits are not consistent with being circular; the average value of r$_{apo}$ / r$_{peri}$ is $\sim$1.9 (equivalent to eccentricity of $\sim$0.31) for both clusters. Lastly, we find that the clusters do not share a common orbit, challenging one proposed formation scenario in which the clusters formed from molecular clouds on the open stream orbit derived by Kruijssen et al. (2015). Meanwhile, our constraints on the birth location and velocity of the clusters offer mild support for a scenario in which the clusters formed via collisions between gas clouds on the x1 and x2 bar orbit families. △ Less

Submitted 17 August, 2022; originally announced August 2022.

Comments: Accepted for publication in ApJ. 38 pages, 25 figures. Proper motion catalogs included in ancillary materials

Calibration and performance of the readout system based on switched capacitor arrays for the Large-Sized Telescope of the Cherenkov Telescope Array

Authors: Seiya Nozaki, Kyosuke Awai, Aya Bamba, Juan Abel Barrio, Maria Isabel Bernardos, Oscar Blanch, Joan Boix, Franca Cassol, Yuki Choushi, Carlos Delgado, Carlos Diaz, Nadia Fouque, Lluis Freixas, Pawel Gliwny, Shunichi Gunji, Daniela Hadasch, Dirk Hoffmann, Julien Houles, Yusuke Inome, Yuki Iwamura, Léa Jouvin, Hideaki Katagiri, Kiomei Kawamura, Daniel Kerszberg, Yusuke Konno , et al. (37 additional authors not shown)

Abstract: The Cherenkov Telescope Array (CTA) is the next-generation ground-based very-high-energy gamma-ray observatory. The Large-Sized Telescope (LST) of CTA is designed to detect gamma rays between 20 GeV and a few TeV with a 23-meter diameter mirror. We have developed the focal plane camera of the first LST, which has 1855 photomultiplier tubes (PMTs) and the readout system which samples a PMT waveform… ▽ More The Cherenkov Telescope Array (CTA) is the next-generation ground-based very-high-energy gamma-ray observatory. The Large-Sized Telescope (LST) of CTA is designed to detect gamma rays between 20 GeV and a few TeV with a 23-meter diameter mirror. We have developed the focal plane camera of the first LST, which has 1855 photomultiplier tubes (PMTs) and the readout system which samples a PMT waveform at GHz with switched capacitor arrays, Domino Ring Sampler ver4 (DRS4). To measure the precise pulse charge and arrival time of Cherenkov signals, we developed a method to calibrate the output voltage of DRS4 and the sampling time interval, as well as an analysis method to correct the spike noise of DRS4. Since the first LST was inaugurated in 2018, we have performed the commissioning tests and calibrated the camera. We characterised the camera in terms of the charge pedestal under various conditions of the night sky background, the charge resolution of each pixel, the charge uniformity of the whole camera, and the time resolutions with a test pulse and calibration laser. △ Less

Submitted 13 March, 2022; originally announced March 2022.

Comments: 14 pages, 13 figures, Proceedings of the SPIE Astronomical Telescopes + Instrumentation, 114470H (2020)

Kinematics and Feedback in H II regions in the Dwarf Starburst Galaxy IC 10

Authors: Maren Cosens, Shelley A. Wright, Norman Murray, Lee Armus, Karin Sandstrom, Tuan Do, Kirsten Larson, Gregory Martinez, Sanchit Sabhlok, Andrey Vayner, James Wiley

Abstract: We present a survey of the central region of the nearest starburst galaxy, IC 10, using the W. M. Keck Observatory Keck Cosmic Web Imager (KCWI) at high spectral and spatial resolution. We map the central starburst of IC 10 to sample the kinematic and ionization properties of the individual star-forming regions. Using the low spectral resolution mode of KCWI we map the oxygen abundance and with th… ▽ More We present a survey of the central region of the nearest starburst galaxy, IC 10, using the W. M. Keck Observatory Keck Cosmic Web Imager (KCWI) at high spectral and spatial resolution. We map the central starburst of IC 10 to sample the kinematic and ionization properties of the individual star-forming regions. Using the low spectral resolution mode of KCWI we map the oxygen abundance and with the high spectral resolution mode we identify 46 individual H II regions. These H II regions have an average radius of 4.0 pc, star formation rate $\sim1.3\times10^{-4}$ M$_\odot$ yr$^{-1}$, and velocity dispersion $\sim$16 km s$^{-1}$. None of the H II regions appear to be virialized ($\rm α_{vir}>>1$), and, on average, they show evidence of ongoing expansion. IC 10's H II regions are offset from the star forming region size-luminosity scaling relationships, as well as Larson's Law that relates size and velocity dispersion. We investigate the balance of inward and outward pressure, $\rm P_{in}$ and $\rm P_{out}$, finding $\rm P_{out}>P_{in}$ in 89% of H II regions, indicating feedback driven expansion even in these low mass H II regions. We find warm gas pressure ($\rm P_{gas}$) provides the dominant contribution to the outward pressure ($\rm P_{out}$). This counteracts the inward pressure which is dominated by turbulence in the surrounding gas rather than self-gravity. Five H II regions show evidence of outflows which are most likely supported by either stellar winds (2 regions) or champagne flows (3 regions). These observations provide new insights into the state of the star-forming regions in IC 10 and negative feedback from low mass clusters. △ Less

Submitted 8 February, 2022; originally announced February 2022.

Comments: 57 pages, 35 figures, 9 tables, accepted for publication in The Astrophysical Journal

Commissioning of the camera of the first Large Size Telescope of the Cherenkov Telescope Array

Authors: T. Saito, C. Delgado, O. Blanch, M. Artero, J. A. Barrio, F. Cassol, C. Diaz, D. Hadasch, D. Hoffmann, J. Houles, Y. Inome, M. Iori, L. Jouvin, D. Kerszberg, Y. Kobayashi, H. Kubo, G. Martinez, D. Mazin, E. Moretti, T. Nakamori, S. Nozaki, T. Oka, A. Okumura, M. Palatiello, M. Polo , et al. (8 additional authors not shown)

Abstract: The first Large Size Telescope (LST-1) of the Cherenkov Telescope Array has been operational since October 2018 at La Palma, Spain. We report on the results obtained during the camera commissioning. The noise level of the readout is determined as a 0.2 p.e. level. The gain of PMTs are well equalized within 2\% variation, using the calibration flash system. The effect of the night sky background on… ▽ More The first Large Size Telescope (LST-1) of the Cherenkov Telescope Array has been operational since October 2018 at La Palma, Spain. We report on the results obtained during the camera commissioning. The noise level of the readout is determined as a 0.2 p.e. level. The gain of PMTs are well equalized within 2\% variation, using the calibration flash system. The effect of the night sky background on the signal readout noise as well as the PMT gain estimation are also well evaluated. Trigger thresholds are optimized for the lowest possible gamma-ray energy threshold and the trigger distribution synchronization has been achieved within 1~ns precision. Automatic rate control realizes the stable observation with 1.5\% rate variation over 3 hours. The performance of the novel DAQ system demonstrates a less than 10\% dead time for 15 kHz trigger rate even with sophisticated online data correction. △ Less

Submitted 4 August, 2021; originally announced August 2021.

Comments: 10 pages, 9 figures, Proceedings of the 37th International Cosmic Ray Conference (ICRC 2021), Berlin, Germany

Report number: PoS-ICRC2021-718

Simple and statistically sound recommendations for analysing physical theories

Authors: Shehu S. AbdusSalam, Fruzsina J. Agocs, Benjamin C. Allanach, Peter Athron, Csaba Balázs, Emanuele Bagnaschi, Philip Bechtle, Oliver Buchmueller, Ankit Beniwal, Jihyun Bhom, Sanjay Bloor, Torsten Bringmann, Andy Buckley, Anja Butter, José Eliel Camargo-Molina, Marcin Chrzaszcz, Jan Conrad, Jonathan M. Cornell, Matthias Danninger, Jorge de Blas, Albert De Roeck, Klaus Desch, Matthew Dolan, Herbert Dreiner, Otto Eberhardt , et al. (50 additional authors not shown)

Abstract: Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by mul… ▽ More Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by multiple experiments, and random or grid sampling of model parameters. Whilst these methods are easy to apply, they exhibit pathologies even in low-dimensional parameter spaces, and quickly become problematic to use and interpret in higher dimensions. In this article we give clear guidance for going beyond these procedures, suggesting where possible simple methods for performing statistically sound inference, and recommendations of readily-available software tools and standards that can assist in doing so. Our aim is to provide any physicists lacking comprehensive statistical training with recommendations for reaching correct scientific conclusions, with only a modest increase in analysis burden. Our examples can be reproduced with the code publicly available at https://doi.org/10.5281/zenodo.4322283. △ Less

Submitted 11 April, 2022; v1 submitted 17 December, 2020; originally announced December 2020.

Comments: 15 pages, 4 figures. extended discussions. closely matches version accepted for publication

Report number: PSI-PR-20-23, BONN-TH-2020-11, CP3-20-59, KCL-PH-TH/2020-75, P3H-20-080, TTP20-044, TUM-HEP-1310/20, IFT-UAM/CSIC-20-180, TTK-20-47, CERN-TH-2020-215, FTPI-MINN-20-36, UMN-TH-4005/20, HU-EP-20/37, DESY 20-222, ADP-20-33/T1143, Imperial/TP/2020/RT/04, UCI-TR-2020-19, gambit-review-2020

Journal ref: Rep. Prog. Phys. 85 052201 (2022)

Distribution and habitability of (meta)stable brines on present-day Mars

Authors: Edgard G. Rivera-Valentín, Vincent F. Chevrier, Alejandro Soto, Germán Martínez

Abstract: Special Regions on Mars are defined as environments able to host liquid water that meets certain temperature and water activity requirements that allow known terrestrial organisms to replicate, and therefore could be habitable. Such regions would be a concern for planetary protection policies owing to the potential for forward contamination (biological contamination from Earth). Pure liquid water… ▽ More Special Regions on Mars are defined as environments able to host liquid water that meets certain temperature and water activity requirements that allow known terrestrial organisms to replicate, and therefore could be habitable. Such regions would be a concern for planetary protection policies owing to the potential for forward contamination (biological contamination from Earth). Pure liquid water is unstable on the Martian surface, but brines may be present. Experimental work has shown that brines persist beyond their predicted stability region, leading to metastable liquids. Here we show that (meta)stable brines can form and persist from the equator to high latitudes on the surface of Mars for a few percent of the year for up to six consecutive hours, a broader range than previously thought. However, only the lowest eutectic solutions can form, leading to brines with temperatures of less than 225 K. Our results indicate that (meta)stable brines on the Martian surface and shallow subsurface (a few centimeters deep) are not habitable because their water activities and temperatures fall outside the known tolerances for terrestrial life. Furthermore, (meta)stable brines do not meet the Special Regions requirements, reducing the risk for forward contamination and easing threats related to the exploration of the Martian surface. △ Less

Submitted 30 November, 2020; originally announced December 2020.

Journal ref: Rivera-Valentin, E. G., Chevrier, V. F., Soto, A., Martinez, G. (2020) Distribution and habitability of (meta)stable brines on present-day Mars. Nat Astron 4, 756-761

Rapid Variability of Sgr A* across the Electromagnetic Spectrum

Authors: G. Witzel, G. Martinez, S. P. Willner, E. E. Becklin, 4 H. Boyce, T. Do, A. Eckart, G. G. Fazio, A. Ghez, M. A. Gurwell, D. Haggard, R. Herrero-Illana, J. L. Hora, Z. Li, J. Liu, N. Marchili, Mark R. Morris, Howard A. Smith, M. Subroweit, J. A. Zensus

Abstract: Sagittarius A* (Sgr A*) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We have analyzed a comprehensive submillimeter (including new observations simultaneous with NIR monitoring), NIR, and 2-8 keV dataset. Submillimeter variations tend to lag those in the NIR by $\sim$30 minutes. An approximate Bayesian computation (ABC)… ▽ More Sagittarius A* (Sgr A*) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We have analyzed a comprehensive submillimeter (including new observations simultaneous with NIR monitoring), NIR, and 2-8 keV dataset. Submillimeter variations tend to lag those in the NIR by $\sim$30 minutes. An approximate Bayesian computation (ABC) fit to the X-ray first-order structure function shows significantly less power at short timescales in the X-rays than in the NIR. Less X-ray variability at short timescales combined with the observed NIR-X-ray correlations means the variability can be described as the result of two strictly correlated stochastic processes, the X-ray process being the low-pass-filtered version of the NIR process. The NIR--X-ray linkage suggests a simple radiative model: a compact, self-absorbed synchrotron sphere with high-frequency cutoff close to NIR frequencies plus a synchrotron self-Compton scattering component at higher frequencies. This model, with parameters fit to the submillimeter, NIR, and X-ray structure functions, reproduces the observed flux densities at all wavelengths, the statistical properties of all light curves, and the time lags between bands. The fit also gives reasonable values for physical parameters such as magnetic flux density $B\approx13$ G, source size $L \approx2.2R_{S}$, and high-energy electron density $n_{e}\approx4\times10^{7}$ cm$^{-3}$. An animation illustrates typical light curves, and we make public the parameter chain of our Bayesian analysis, the model implementation, and the visualization code. △ Less

Submitted 7 June, 2021; v1 submitted 18 November, 2020; originally announced November 2020.

Comments: Accepted by ApJS, June 3, 2021. The journal version of figure 21 is animated. The animation can also be found here: https://doi.org/10.17617/1.kctx3s25. This version (version 2) has been revised according to the referee's suggestions and includes optimized figures and some changes and corrections of the text; no scientific conclusions have changed

Searching for MHz Gravitational Waves from Harmonic Sources

Authors: Jeronimo G. C. Martinez, Brittany Kamai

Abstract: A MHz gravitational wave search for harmonic sources was conducted using a 704-hr dataset obtained from the Holometer, a pair of 40-meter power recycled Michelson interferometers. Our search was designed to look for cosmic string loops and eccentric black hole binaries in an entirely unexplored frequency range from 1 to 25 MHz. The measured cross-spectral density between both interferometers was u… ▽ More A MHz gravitational wave search for harmonic sources was conducted using a 704-hr dataset obtained from the Holometer, a pair of 40-meter power recycled Michelson interferometers. Our search was designed to look for cosmic string loops and eccentric black hole binaries in an entirely unexplored frequency range from 1 to 25 MHz. The measured cross-spectral density between both interferometers was used to perform four different searches. First, we search to identify any fundamental frequencies bins that have excess power above 5$σ$. Second, we reduce the per-bin threshold on any individual frequency bin by employing that a fundamental frequency and its harmonics all collectively lie above a threshold. We vary the number of harmonics searched over from $n= 4$ up to $n=23$. Third, we perform an agnostic approach to identify harmonic candidates that may have a single contaminated frequency bin or follow a power-law dependence. Lastly, we expand on the agnostic approach for individual candidates and search for a potential underlying population of harmonic sources. Each method was tested on the interferometer dataset, as well as a dark noise, photon shot-noise-limited, and simulated Gaussian-noise datasets. We conclude that these four different search methods did not find any candidate frequencies that would be consistent with harmonic sources. This work presents a new way of searching for gravitational wave candidates, which allowed us to survey a previously unexplored frequency range. △ Less

Submitted 12 October, 2020; originally announced October 2020.

Comments: 14 pages, 5 figures

Journal ref: Classical and Quantum Gravity, Volume 37, Number 20, Pages 205006, Year 2020

Sensitivity of the Cherenkov Telescope Array for probing cosmology and fundamental physics with gamma-ray propagation

Authors: The Cherenkov Telescope Array Consortium, :, H. Abdalla, H. Abe, F. Acero, A. Acharyya, R. Adam, I. Agudo, A. Aguirre-Santaella, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves B, L. Amati, E. Amato, G. Ambrosi, E. O. Angüner, A. Araudo, T. Armstrong, F. Arqueros, L. Arrabito, K. Asano, Y. Ascasíbar, M. Ashley , et al. (474 additional authors not shown)

Abstract: The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $γ$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $γ$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nucle… ▽ More The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $γ$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $γ$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nuclei (AGN) and of their relativistic jets. Observations of AGN with CTA will enable a measurement of $γ$-ray absorption on the extragalactic background light with a statistical uncertainty below 15% up to a redshift $z=2$ and to constrain or detect $γ$-ray halos up to intergalactic-magnetic-field strengths of at least 0.3pG. Extragalactic observations with CTA also show promising potential to probe physics beyond the Standard Model. The best limits on Lorentz invariance violation from $γ$-ray astronomy will be improved by a factor of at least two to three. CTA will also probe the parameter space in which axion-like particles could constitute a significant fraction, if not all, of dark matter. We conclude on the synergies between CTA and other upcoming facilities that will foster the growth of $γ$-ray cosmology. △ Less

Submitted 26 February, 2021; v1 submitted 3 October, 2020; originally announced October 2020.

Comments: 71 pages (including affiliations and references), 13 figures, 6 tables. Accepted in JCAP; matches published version. Corresponding authors: Jonathan Biteau, Julien Lefaucheur, Humberto Martinez-Huerta, Manuel Meyer, Santiago Pita, Ievgen Vovk

Journal ref: JCAP 02 (2021) 048

Revealing the Formation of the Milky Way Nuclear Star Cluster via Chemo-Dynamical Modeling

Authors: Tuan Do, Gregory David Martinez, Wolfgang Kerzendorf, Anja Feldmeier-Krause, Manuel Arca Sedda, Nadine Neumayer, Alessia Gualandris

Abstract: The Milky Way nuclear star cluster (MW NSC) has been used as a template to understand the origin and evolution of galactic nuclei and the interaction of nuclear star clusters with supermassive black holes. It is the only nuclear star cluster with a supermassive black hole where we can resolve individual stars to measure their kinematics and metal abundance to reconstruct its formation history. Her… ▽ More The Milky Way nuclear star cluster (MW NSC) has been used as a template to understand the origin and evolution of galactic nuclei and the interaction of nuclear star clusters with supermassive black holes. It is the only nuclear star cluster with a supermassive black hole where we can resolve individual stars to measure their kinematics and metal abundance to reconstruct its formation history. Here, we present results of the first chemo-dynamical model of the inner 1 pc of the MW NSC using metallicity and radial velocity data from the KMOS spectrograph on the Very Large Telescope. We find evidence for two kinematically and chemically distinct components in this region. The majority of the stars belong to a previously known super-solar metallicity component with a rotation axis perpendicular to the Galactic plane. However, we identify a new kinematically distinct sub-solar metallicity component which contains about 7\% of the stars and appears to be rotating faster than the main component with a rotation axis that may be misaligned. This second component may be evidence for an infalling star cluster or remnants of a dwarf galaxy, merging with the MW NSC. These measurements show that the combination of chemical abundances with kinematics is a promising method to directly study the MW NSC's origin and evolution. △ Less

Submitted 4 September, 2020; originally announced September 2020.

Comments: 6 pages, 5 figures, accepted to ApJ Letters

Young stellar cluster dilution near supermassive black holes: the impact of Vector Resonant Relaxation on neighbour separation

Authors: Juan Giral Martínez, Jean-Baptiste Fouvry, Christophe Pichon

Abstract: We investigate the rate of orbital orientation dilution of young stellar clusters in the vicinity of supermassive black holes. Within the framework of vector resonant relaxation, we predict the time evolution of the two-point correlation function of the stellar orbital plane orientations as a function of their initial angular separation and diversity in orbital parameters (semi-major axis, eccentr… ▽ More We investigate the rate of orbital orientation dilution of young stellar clusters in the vicinity of supermassive black holes. Within the framework of vector resonant relaxation, we predict the time evolution of the two-point correlation function of the stellar orbital plane orientations as a function of their initial angular separation and diversity in orbital parameters (semi-major axis, eccentricity). As expected, the larger the spread in initial orientations and orbital parameters, the more efficient the dilution of a given set of co-eval stars, with a characteristic timescale set up by the coherence time of the background potential fluctuations. A Markovian prescription which matches numerical simulations allows us to efficiently probe the underlying kinematic properties of the unresolved nucleus when requesting consistency with a given dilution efficiency, imposed by the observed stellar disc within the one arcsecond of Sgr A*. As a proof of concept, we compute maps of constant dilution times as a function of the semi major axis cusp index and fraction of intermediate mass black holes in the old background stellar cluster. This computation suggests that vector resonant relaxation should prove useful in this context since it impacts orientations on timescales comparable to the stars' age. △ Less

Submitted 5 August, 2020; originally announced August 2020.

Comments: 20 pages, 16 figures, submitted to MNRAS

Sensitivity of the Cherenkov Telescope Array to a dark matter signal from the Galactic centre

Authors: The Cherenkov Telescope Array Consortium, :, A. Acharyya, R. Adam, C. Adams, I. Agudo, A. Aguirre-Santaella, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves Batista, L. Amati, G. Ambrosi, E. O. Angüner, L. A. Antonelli, C. Aramo, A. Araudo, T. Armstrong, F. Arqueros, K. Asano, Y. Ascasíbar, M. Ashley, C. Balazs, O. Ballester , et al. (427 additional authors not shown)

Abstract: We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models giv… ▽ More We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models given a standard cuspy profile of the dark matter density distribution. Importantly, even for a cored profile, the projected sensitivity of CTA will be sufficient to probe various well-motivated models of thermally produced dark matter at the TeV scale. This is due to CTA's unprecedented sensitivity, angular and energy resolutions, and the planned observational strategy. The survey of the inner Galaxy will cover a much larger region than corresponding previous observational campaigns with imaging atmospheric Cherenkov telescopes. CTA will map with unprecedented precision the large-scale diffuse emission in high-energy gamma rays, constituting a background for dark matter searches for which we adopt state-of-the-art models based on current data. Throughout our analysis, we use up-to-date event reconstruction Monte Carlo tools developed by the CTA consortium, and pay special attention to quantifying the level of instrumental systematic uncertainties, as well as background template systematic errors, required to probe thermally produced dark matter at these energies. "Full likelihood tables complementing our analysis are provided here [ https://doi.org/10.5281/zenodo.4057987 ]" △ Less

Submitted 30 January, 2021; v1 submitted 31 July, 2020; originally announced July 2020.

Comments: 68 pages (including references) and 26 figures; text identical to the version published in JCAP

Journal ref: JCAP01(2021)057

Measuring Mars Atmospheric Winds From Orbit

Authors: Scott Guzewich J. B. Abshire. M. M. Baker, J. M. Battalio, T. Bertrand, A. J. Brown, A. Colaprete, A. M. Cook, D. R. Cremons, M. M. Crismani, A. I. Dave, M. Day, M. -C. Desjean, M. Elrod, L. K. Fenton, J. Fisher, L. L. Gordley, P. O. Hayne, N. G. Heavens, J. L. Hollingsworth, D. Jha, V. Jha, M. A. Kahre, A. SJ. Khayat, A. M. Kling, S. R. Lewis, B. T. Marshall , et al. (16 additional authors not shown)

Abstract: Wind is the process that connects Mars' climate system. Measurements of Mars atmospheric winds from orbit would dramatically advance our understanding of Mars and help prepare for human exploration of the Red Planet. Multiple instrument candidates are in development and will be ready for flight in the next decade. We urge the Decadal Survey to make these measurements a priority for 2023-2032. Wind is the process that connects Mars' climate system. Measurements of Mars atmospheric winds from orbit would dramatically advance our understanding of Mars and help prepare for human exploration of the Red Planet. Multiple instrument candidates are in development and will be ready for flight in the next decade. We urge the Decadal Survey to make these measurements a priority for 2023-2032. △ Less

Submitted 10 July, 2020; originally announced July 2020.

Comments: A White Paper submitted to the Planetary Science and Astrobiology Decadal Survey 2023-2032

A population of dust-enshrouded objects orbiting the Galactic black hole

Authors: Anna Ciurlo, Randall D. Campbell, Mark R. Morris, Tuan Do, Andrea M. Ghez, Aurelien Hees, Breann N. Sitarski, Kelly Kosmo O'Neil, Devin S. Chu, Gregory D. Martinez, Smadar Naoz, Alexander P. Stephan

Abstract: The central 0.1 parsecs of the Milky Way host a supermassive black hole identified with the position of the radio and infrared source Sagittarius A*, a cluster of young, massive stars (the S stars) and various gaseous features. Recently, two unusual objects have been found to be closely orbiting Sagittarius A*: the so-called G sources, G1 and G2. These objects are unresolved (having a size of the… ▽ More The central 0.1 parsecs of the Milky Way host a supermassive black hole identified with the position of the radio and infrared source Sagittarius A*, a cluster of young, massive stars (the S stars) and various gaseous features. Recently, two unusual objects have been found to be closely orbiting Sagittarius A*: the so-called G sources, G1 and G2. These objects are unresolved (having a size of the order of 100 astronomical units, except at periapse, where the tidal interaction with the black hole stretches them along the orbit) and they show both thermal dust emission and line emission from ionized gas. G1 and G2 have generated attention because they appear to be tidally interacting with the supermassive Galactic black hole, possibly enhancing its accretion activity. No broad consensus has yet been reached concerning their nature: the G objects show the characteristics of gas and dust clouds but display the dynamical properties of stellar-mass objects. Here we report observations of four additional G objects, all lying within 0.04 parsecs of the black hole and forming a class that is probably unique to this environment. The widely varying orbits derived for the six G objects demonstrate that they were commonly but separately formed. △ Less

Submitted 22 January, 2020; originally announced January 2020.

Journal ref: Nature volume 577, pages337-340(2020)

Unprecedented variability of Sgr A* in NIR

Authors: Tuan Do, Gunther Witzel, Abhimat K. Gautam, Zhuo Chen, Andrea M. Ghez, Mark R. Morris, Eric E. Becklin, Anna Ciurlo, Matthew Hosek Jr., Gregory D. Martinez, Keith Matthews, Shoko Sakai, Rainer Schödel

Abstract: The electromagnetic counterpart to the Galactic center supermassive black hole, Sgr A*, has been observed in the near-infrared for over 20 years and is known to be highly variable. We report new Keck Telescope observations showing that Sgr A* reached much brighter flux levels in 2019 than ever measured at near-infrared wavelengths. In the K$^\prime$ band, Sgr A* reached flux levels of $\sim6$ mJy,… ▽ More The electromagnetic counterpart to the Galactic center supermassive black hole, Sgr A*, has been observed in the near-infrared for over 20 years and is known to be highly variable. We report new Keck Telescope observations showing that Sgr A* reached much brighter flux levels in 2019 than ever measured at near-infrared wavelengths. In the K$^\prime$ band, Sgr A* reached flux levels of $\sim6$ mJy, twice the level of the previously observed peak flux from $>13,000$ measurements over 130 nights with the VLT and Keck Telescopes. We also observe a factor of 75 change in flux over a 2-hour time span with no obvious color changes between 1.6 $μ$m and 2.1 $μ$m. The distribution of flux variations observed this year is also significantly different than the historical distribution. Using the most comprehensive statistical model published, the probability of a single night exhibiting peak flux levels observed this year, given historical Keck observations, is less than $0.3\%$. The probability to observe the flux levels similar to all 4 nights of data in 2019 is less than $0.05\%$. This increase in brightness and variability may indicate a period of heightened activity from Sgr A* or a change in its accretion state. It may also indicate that the current model is not sufficient to model Sgr A* at high flux levels and should be updated. Potential physical origins of Sgr A*'s unprecedented brightness may be from changes in the accretion-flow as a result of the star S0-2's closest passage to the black hole in 2018 or from a delayed reaction to the approach of the dusty object G2 in 2014. Additional multi-wavelength observations will be necessary to both monitor Sgr A* for potential state changes and to constrain the physical processes responsible for its current variability. △ Less

Submitted 5 August, 2019; originally announced August 2019.

Comments: 9 pages, 4 figures, accepted to ApJ Letters

Relativistic redshift of the star S0-2 orbiting the Galactic center supermassive black hole

Authors: Tuan Do, Aurelien Hees, Andrea Ghez, Gregory D. Martinez, Devin S. Chu, Siyao Jia, Shoko Sakai, Jessica R. Lu, Abhimat K. Gautam, Kelly Kosmo O'Neil, Eric E. Becklin, Mark R. Morris, Keith Matthews, Shogo Nishiyama, Randy Campbell, Samantha Chappell, Zhuo Chen, Anna Ciurlo, Arezu Dehghanfar, Eulalia Gallego-Cano, Wolfgang E. Kerzendorf, James E. Lyke, Smadar Naoz, Hiromi Saida, Rainer Schödel , et al. (4 additional authors not shown)

Abstract: General Relativity predicts that a star passing close to a supermassive black hole should exhibit a relativistic redshift. We test this using observations of the Galactic center star S0-2. We combine existing spectroscopic and astrometric measurements from 1995-2017, which cover S0-2's 16-year orbit, with measurements in 2018 March to September which cover three events during its closest approach… ▽ More General Relativity predicts that a star passing close to a supermassive black hole should exhibit a relativistic redshift. We test this using observations of the Galactic center star S0-2. We combine existing spectroscopic and astrometric measurements from 1995-2017, which cover S0-2's 16-year orbit, with measurements in 2018 March to September which cover three events during its closest approach to the black hole. We detect the combination of special relativistic- and gravitational-redshift, quantified using a redshift parameter, $Υ$. Our result, $Υ=0.88 \pm 0.17$, is consistent with General Relativity ($Υ=1$) and excludes a Newtonian model ($Υ=0$ ) with a statistical significance of 5 $σ$. △ Less

Submitted 24 July, 2019; originally announced July 2019.

Comments: 78 pages, 21 figures, accepted to Science, astrometry and radial velocity measurements, and posterior chains available here, http://www.astro.ucla.edu/~ghezgroup/gc/gr/data.zip

The Discovery of Six Recycled Pulsars from the Arecibo 327-MHz Drift-Scan Pulsar Survey

Authors: J. G. Martinez, P. Gentile, P. C. C. Freire, K. Stovall, J. S. Deneva, G. Desvignes, F. A. Jenet, M. A. McLaughlin, M. Bagchi, T. Devine

Abstract: Recycled pulsars are old ($\gtrsim10^{8}$ yr) neutron stars that are descendants from close, interacting stellar systems. In order to understand their evolution and population, we must find and study the largest number possible of recycled pulsars in a way that is as unbiased as possible. In this work, we present the discovery and timing solutions of five recycled pulsars in binary systems (PSRs J… ▽ More Recycled pulsars are old ($\gtrsim10^{8}$ yr) neutron stars that are descendants from close, interacting stellar systems. In order to understand their evolution and population, we must find and study the largest number possible of recycled pulsars in a way that is as unbiased as possible. In this work, we present the discovery and timing solutions of five recycled pulsars in binary systems (PSRs J0509$+$0856, J0709$+$0458, J0732$+$2314, J0824$+$0028, J2204$+$2700) and one isolated millisecond pulsar (PSR J0154$+$1833). These were found in data from the Arecibo 327-MHz Drift-Scan Pulsar Survey (AO327). All these pulsars have a low dispersion measure (DM) ($\lesssim 45 \, \rm{pc}\, cm^{-3}$), and have a DM-determined distance of $\lesssim$ 3 kpc. Their timing solutions, have data spans ranging from 1 to $\sim$ 7 years, include precise estimates of their spin and astrometric parameters, and for the binaries, precise estimates of their Keplerian binary parameters. Their orbital periods range from about 4 to 815 days and the minimum companion masses (assuming a pulsar mass of 1.4 $\rm{M_{\odot}}$) range from $\sim$ 0.06--1.11 $\rm{M_{\odot}}$. For two of the binaries we detect post-Keplerian parameters; in the case of PSR~J0709$+$0458 we measure the component masses but with a low precision, in the not too distant future the measurement of the rate of advance of periastron and the Shapiro delay will allow very precise mass measurements for this system. Like several other systems found in the AO327 data, PSRs J0509$+$0854, J0709$+$0458 and J0732$+$2314 are now part of the NANOGrav timing array for gravitational wave detection. △ Less

Submitted 12 June, 2019; originally announced June 2019.

An adaptive scheduling tool to optimize measurements to reach a scientific objective: methodology and application to the measurements of stellar orbits in the Galactic Center

Authors: A. Hees, A. Dehghanfar, T. Do, A. M. Ghez, G. D. Martinez, R. Campbell, J. R. Lu

Abstract: In various fields of physics and astronomy, access to experimental facilities or to telescopes is becoming more and more competitive and limited. It becomes therefore important to optimize the type of measurements and their scheduling to reach a given scientific objective and to increase the chances of success of a scientific project. In this communication, extending the work of Ford (2008) and of… ▽ More In various fields of physics and astronomy, access to experimental facilities or to telescopes is becoming more and more competitive and limited. It becomes therefore important to optimize the type of measurements and their scheduling to reach a given scientific objective and to increase the chances of success of a scientific project. In this communication, extending the work of Ford (2008) and of Loredo et al. (2012), we present an efficient adaptive scheduling tool aimed at prioritzing measurements in order to reach a scientific goal. The algorithm, based on the Fisher matrix, can be applied to a wide class of measurements. We present this algorithm in detail and discuss some practicalities such as systematic errors or measurements losses due to contigencies (such as weather, experimental failure, ...). As an illustration, we consider measurements of the short-period star S0-2 in our Galactic Center. We show that the radial velocity measurements at the two turning points of the radial velocity curve are more powerful for detecting the gravitational redshift than measurements at the maximal relativistic signal. We also explicitly present the methodology that was used to plan measurements in order to detect the relativistic redshift considering systematics and possible measurements losses. For the future, we identify the astrometric turning points to be highly sensitive to the relativistic advance of the periastron. Finally, we also identify measurements particularly sensitive to the distance to our Galactic Center: the radial velocities around periastron and the astrometric measurements just before closest approach and at the maximal right ascension astrometric turning point. △ Less

Submitted 7 June, 2019; originally announced June 2019.

Comments: 12 pages, 6 figures, submitted

Journal ref: ApJ 880, 87, 2019

Monte Carlo studies for the optimisation of the Cherenkov Telescope Array layout

Authors: A. Acharyya, I. Agudo, E. O. Angüner, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves Batista, J. -P. Amans, L. Amati, E. Amato, G. Ambrosi, L. A. Antonelli, C. Aramo, T. Armstrong, F. Arqueros, L. Arrabito, K. Asano, H. Ashkar, C. Balazs, M. Balbo, B. Balmaverde, P. Barai, A. Barbano, M. Barkov , et al. (445 additional authors not shown)

Abstract: The Cherenkov Telescope Array (CTA) is the major next-generation observatory for ground-based very-high-energy gamma-ray astronomy. It will improve the sensitivity of current ground-based instruments by a factor of five to twenty, depending on the energy, greatly improving both their angular and energy resolutions over four decades in energy (from 20 GeV to 300 TeV). This achievement will be possi… ▽ More The Cherenkov Telescope Array (CTA) is the major next-generation observatory for ground-based very-high-energy gamma-ray astronomy. It will improve the sensitivity of current ground-based instruments by a factor of five to twenty, depending on the energy, greatly improving both their angular and energy resolutions over four decades in energy (from 20 GeV to 300 TeV). This achievement will be possible by using tens of imaging Cherenkov telescopes of three successive sizes. They will be arranged into two arrays, one per hemisphere, located on the La Palma island (Spain) and in Paranal (Chile). We present here the optimised and final telescope arrays for both CTA sites, as well as their foreseen performance, resulting from the analysis of three different large-scale Monte Carlo productions. △ Less

Submitted 2 April, 2019; originally announced April 2019.

Comments: 48 pages, 16 figures, accepted for publication in Astroparticle Physics

Astro2020 Science White Paper: The fundamentals of outflows from evolved stars

Authors: E. De Beck, M. L. Boyer, V. Bujarrabal, L. Decin, J. P. Fonfría, M. Groenewegen, S. Höfner, O. Jones, T. Kaminski, M. Maercker, P. Marigo, M. Matsuura, M. Meixner, G. Quintana Lacaci Martínez, P. Scicluna, R. Szczerba, L. Velilla Prieto, W. Vlemmings, M. Wiedner

Abstract: Models of the chemical evolution of the interstellar medium, galaxies, and the Universe rely on our understanding of the amounts and chemical composition of the material returned by stars and supernovae. Stellar yields are obtained from stellar-evolution models, which currently lack predictive prescriptions of stellar mass loss, although it significantly affects stellar lifetimes, nucleosynthesis,… ▽ More Models of the chemical evolution of the interstellar medium, galaxies, and the Universe rely on our understanding of the amounts and chemical composition of the material returned by stars and supernovae. Stellar yields are obtained from stellar-evolution models, which currently lack predictive prescriptions of stellar mass loss, although it significantly affects stellar lifetimes, nucleosynthesis, and chemical ejecta. Galaxy properties are derived from observations of the integrated light of bright member stars. Stars in the late stages of their evolution are among the infrared-brightest objects in galaxies. An unrealistic treatment of the mass-loss process introduces significant uncertainties in galaxy properties derived from their integrated light. We describe current efforts and future needs and opportunities to characterize AGB outflows: driving mechanisms, outflow rates, underlying fundamental physical and chemical processes such as dust grain formation, and dependency of these on metallicity. △ Less

Submitted 28 March, 2019; originally announced March 2019.

Comments: Submitted to Astro2020 decadal survey; 8 pages, 2 figures

The Green Bank North Celestial Cap Pulsar Survey. IV: Four New Timing Solutions

Authors: R. J. Aloisi, A. Cruz, L. Daniels, N. Meyers, R. Roekle, A. Schuett, J. K. Swiggum, M. E. DeCesar, D. L. Kaplan, R. S. Lynch, K. Stovall, Lina Levin, A. M. Archibald, S. Banaszak, C. M. Biwer, J. Boyles, P. Chawla, L. P. Dartez, B. Cui, D. F. Day, A. J. Ford, J. Flanigan, E. Fonseca, J. W. T. Hessels, J. Hinojosa , et al. (18 additional authors not shown)

Abstract: We present timing solutions for four pulsars discovered in the Green Bank Northern Celestial Cap (GBNCC) survey. All four pulsars are isolated with spin periods between 0.26$\,$s and 1.84$\,$s. PSR J0038$-$2501 has a 0.26$\,$s period and a period derivative of ${7.6} \times {10}^{-19}\,{\rm s\,s}^{-1}$, which is unusually low for isolated pulsars with similar periods. This low period derivative ma… ▽ More We present timing solutions for four pulsars discovered in the Green Bank Northern Celestial Cap (GBNCC) survey. All four pulsars are isolated with spin periods between 0.26$\,$s and 1.84$\,$s. PSR J0038$-$2501 has a 0.26$\,$s period and a period derivative of ${7.6} \times {10}^{-19}\,{\rm s\,s}^{-1}$, which is unusually low for isolated pulsars with similar periods. This low period derivative may be simply an extreme value for an isolated pulsar or it could indicate an unusual evolution path for PSR J0038$-$2501, such as a disrupted recycled pulsar (DRP) from a binary system or an orphaned central compact object (CCO). Correcting the observed spin-down rate for the Shklovskii effect suggests that this pulsar may have an unusually low space velocity, which is consistent with expectations for DRPs. There is no X-ray emission detected from PSR J0038$-$2501 in an archival swift observation, which suggests that it is not a young orphaned CCO. The high dispersion measure of PSR J1949+3426 suggests a distance of 12.3$\,$kpc. This distance indicates that PSR J1949+3426 is among the most distant 7% of Galactic field pulsars, and is one of the most luminous pulsars. △ Less

Submitted 8 March, 2019; originally announced March 2019.

Comments: 7 pages, 5 figures

An Adaptive Optics Survey of Stellar Variability at the Galactic Center

Authors: Abhimat Krishna Gautam., Tuan Do, Andrea M. Ghez, Mark R. Morris, Gregory D. Martinez, Matthew W. Hosek Jr., Jessica R. Lu, Shoko Sakai, Gunther Witzel, Siyao Jia, Eric E. Becklin, Keith Matthews

Abstract: We present a $\approx 11.5$ year adaptive optics (AO) study of stellar variability and search for eclipsing binaries in the central $\sim 0.4$ pc ($\sim 10''$) of the Milky Way nuclear star cluster. We measure the photometry of 563 stars using the Keck II NIRC2 imager ($K'$-band, $λ_0 = 2.124 \text{ } μ\text{m}$). We achieve a photometric uncertainty floor of $Δm_{K'} \sim 0.03$ ($\approx 3\%$), c… ▽ More We present a $\approx 11.5$ year adaptive optics (AO) study of stellar variability and search for eclipsing binaries in the central $\sim 0.4$ pc ($\sim 10''$) of the Milky Way nuclear star cluster. We measure the photometry of 563 stars using the Keck II NIRC2 imager ($K'$-band, $λ_0 = 2.124 \text{ } μ\text{m}$). We achieve a photometric uncertainty floor of $Δm_{K'} \sim 0.03$ ($\approx 3\%$), comparable to the highest precision achieved in other AO studies. Approximately half of our sample ($50 \pm 2 \%$) shows variability. $52 \pm 5\%$ of known early-type young stars and $43 \pm 4 \%$ of known late-type giants are variable. These variability fractions are higher than those of other young, massive star populations or late-type giants in globular clusters, and can be largely explained by two factors. First, our experiment time baseline is sensitive to long-term intrinsic stellar variability. Second, the proper motion of stars behind spatial inhomogeneities in the foreground extinction screen can lead to variability. We recover the two known Galactic center eclipsing binary systems: IRS 16SW and S4-258 (E60). We constrain the Galactic center eclipsing binary fraction of known early-type stars to be at least $2.4 \pm 1.7\%$. We find no evidence of an eclipsing binary among the young S-stars nor among the young stellar disk members. These results are consistent with the local OB eclipsing binary fraction. We identify a new periodic variable, S2-36, with a 39.43 day period. Further observations are necessary to determine the nature of this source. △ Less

Submitted 12 November, 2018; originally announced November 2018.

Comments: 69 pages, 28 figures, 12 tables. Accepted for publication in The Astrophysical Journal

Journal ref: The Astrophysical Journal, 871:103 (32pp), 2019 January 20

Constraining the potential liquid water environment at Gale crater, Mars

Authors: Edgard G. Rivera-Valentín, Raina V. Gough, Vincent F. Chevrier, Katherine M. Primm, German M. Martínez, Margaret Tolbert

Abstract: The Mars Science Laboratory (MSL) Rover Environmental Monitoring Station (REMS) has now made continuous in-situ meteorological measurements for several martian years at Gale crater, Mars. Of importance in the search for liquid formation are REMS' measurements of ground temperature and in-air measurements of temperature and relative humidity, which is with respect to ice. Such data can constrain th… ▽ More The Mars Science Laboratory (MSL) Rover Environmental Monitoring Station (REMS) has now made continuous in-situ meteorological measurements for several martian years at Gale crater, Mars. Of importance in the search for liquid formation are REMS' measurements of ground temperature and in-air measurements of temperature and relative humidity, which is with respect to ice. Such data can constrain the surface and subsurface stability of brines. Here we use updated calibrations to REMS data and consistent relative humidity comparisons (i.e., w.r.t. liquid vs w.r.t. ice) to investigate the potential formation of surface and subsurface liquids throughout MSL's traverse. We specifically study the potential for the deliquescence of calcium perchlorate. Our data analysis suggests that surface brine formation is not favored within the first 1648 sols as there are only two times (sols 1232 and 1311) when humidity-temperature conditions were within error consistent with a liquid phase. On the other hand, modeling of the subsurface environment would support brine production in the shallow subsurface. Indeed, we find that the shallow subsurface for terrains with low thermal inertia ($Γ\lesssim 300$ J m$^{-2}$ K$^{-1}$ s$^{-1/2}$) may be occasionally favorable to brine formation through deliquescence. Terrains with $Γ\lesssim175$ J m$^{-2}$ K$^{-1}$ s$^{-1/2}$ and albedos of $\gtrsim0.25$ are the most apt to subsurface brine formation. Should brines form, they would occur around Ls 100$^{\circ}$. Their predicted properties would not meet the Special nor Uncertain Region requirements, as such they would not be potential habitable environments to life as we know it. △ Less

Submitted 2 November, 2018; originally announced November 2018.

Comments: 25 pages, 4 figures

Journal ref: Rivera-Valentin, E. G., et al. (2018) Constraining the potential liquid water environment at Gale crater, Mars. Journal of Geophysical Research: Planets, 123, 1156-1167

Size-Luminosity Scaling Relations of Local and Distant Star Forming Regions

Authors: Maren Cosens, Shelley A. Wright, Etsuko Mieda, Norman Murray, Lee Armus, Tuan Do, James E. Larkin, Kirsten Larson, Gregory Martinez, Gregory Walth, Andrey Vayner

Abstract: We investigate star forming scaling relations using Bayesian inference on a comprehensive data sample of low- (z<0.1) and high-redshift (1<z<5) star forming regions. This full data set spans a wide range of host galaxy stellar mass ($M_{*} \sim10^6-10^{11} M_{\odot}$) and clump star formation rates (SFR $ \sim10^{-5}-10^2 M_\odot yr^{-1}$). We fit the power-law relationship between the size (r… ▽ More We investigate star forming scaling relations using Bayesian inference on a comprehensive data sample of low- (z<0.1) and high-redshift (1<z<5) star forming regions. This full data set spans a wide range of host galaxy stellar mass ($M_{*} \sim10^6-10^{11} M_{\odot}$) and clump star formation rates (SFR $ \sim10^{-5}-10^2 M_\odot yr^{-1}$). We fit the power-law relationship between the size (r$_{Hα}$) and luminosity (L$_{Hα}$) of the star forming clumps using the Bayesian statistical modeling tool Stan that makes use of Markov Chain Monte Carlo (MCMC) sampling techniques. Trends in the scaling relationship are explored for the full sample and subsets based on redshift and selection effects between samples. In our investigation we find no evidence of redshift evolution of the size-luminosity scaling relationship, nor a difference in slope between lensed and unlensed data. There is evidence of a break in the scaling relationship between high and low star formation rate surface density ($Σ_{SFR}$) clumps. The size-luminosity power law fit results are L$_{Hα}\sim$ r$_{Hα}^{2.8}$ and L$_{Hα}\sim$ r$_{Hα}^{1.7}$ for low and high $Σ_{SFR}$ clumps, respectively. We present a model where star forming clumps form at locations of gravitational instability and produce an ionized region represented by the Strömgren radius. A radius smaller than the scale height of the disk results in a scaling relationship of $L \propto r^3$ (high $Σ_{SFR}$ clumps), and a scaling of $L \propto r^2$ (low $Σ_{SFR}$ clumps) if the radius is larger than the disk scale height. △ Less

Submitted 24 October, 2018; originally announced October 2018.

Comments: 33 pages, 16 figures, 8 tables. Accepted for publication in The Astrophysical Journal

Improving Orbit Estimates for Incomplete Orbits with a New Approach to Priors -- with Applications from Black Holes to Planets

Authors: K. Kosmo O'Neil, G. D. Martinez, A. Hees, A. M. Ghez, T. Do, G. Witzel, Q. Konopacky, E. E. Becklin, D. S. Chu, J. Lu, K. Matthews, S. Sakai

Abstract: We propose a new approach to Bayesian prior probability distributions (priors) that can improve orbital solutions for low-phase-coverage orbits, where data cover less than approximately 40% of an orbit. In instances of low phase coverage such as with stellar orbits in the Galactic center or with directly-imaged exoplanets, data have low constraining power and thus priors can bias parameter estimat… ▽ More We propose a new approach to Bayesian prior probability distributions (priors) that can improve orbital solutions for low-phase-coverage orbits, where data cover less than approximately 40% of an orbit. In instances of low phase coverage such as with stellar orbits in the Galactic center or with directly-imaged exoplanets, data have low constraining power and thus priors can bias parameter estimates and produce under-estimated confidence intervals. Uniform priors, which are commonly assumed in orbit fitting, are notorious for this. We propose a new observable-based prior paradigm that is based on uniformity in observables. We compare performance of this observable-based prior and of commonly assumed uniform priors using Galactic center and directly-imaged exoplanet (HR 8799) data. The observable-based prior can reduce biases in model parameters by a factor of two and helps avoid under-estimation of confidence intervals for simulations with less than about 40% phase coverage. Above this threshold, orbital solutions for objects with sufficient phase coverage such as S0-2, a short-period star at the Galactic center with full phase coverage, are consistent with previously published results. Below this threshold, the observable-based prior limits prior influence in regions of prior dominance and increases data influence. Using the observable-based prior, HR 8799 orbital analyses favor lower eccentricity orbits and provide stronger evidence that the four planets have a consistent inclination around 30 degrees to within 1-sigma. This analysis also allows for the possibility of coplanarity. We present metrics to quantify improvements in orbital estimates with different priors so that observable-based prior frameworks can be tested and implemented for other low-phase-coverage orbits. △ Less

Submitted 10 June, 2019; v1 submitted 14 September, 2018; originally announced September 2018.

Comments: Published in AJ. 23 pages, 14 figures. Monte Carlo chains are available in the published article, or are available upon request

Journal ref: AJ, 158, 4 (2019)

PSR J2234+0611: A new laboratory for stellar evolution

Authors: K. Stovall, P. C. C. Freire, J. Antoniadis, M. Bagchi, J. S. Deneva, N. Garver-Daniels, J. G. Martinez, M. A. McLaughlin, Z. Arzoumanian, H. Blumer, P. R. Brook, H. T. Cromartie, P. B. Demorest, M. E. Decesar, T. Dolch, J. A. Ellis, R. D. Ferdman, E. C. Ferrara, E. Fonseca, P. A. Gentile, M. L. Jones, M. T. Lam, D. R. Lorimer, R. S. Lynch, C. Ng , et al. (8 additional authors not shown)

Abstract: We report timing results for PSR J2234+0611, a 3.6-ms pulsar in a 32-day, eccentric (e = 0.13) orbit with a helium white dwarf companion discovered as part of the Arecibo Observatory 327 MHz drift scan survey. The precise timing and the eccentric nature of the orbit allow precise measurements of an unusual number of parameters: a) a precise proper motion of 27.10(3) mas/yr and a parallax of 1.05(4… ▽ More We report timing results for PSR J2234+0611, a 3.6-ms pulsar in a 32-day, eccentric (e = 0.13) orbit with a helium white dwarf companion discovered as part of the Arecibo Observatory 327 MHz drift scan survey. The precise timing and the eccentric nature of the orbit allow precise measurements of an unusual number of parameters: a) a precise proper motion of 27.10(3) mas/yr and a parallax of 1.05(4) mas resulting in a pulsar distance of 0.95(4) kpc; this allows a precise estimate of the transverse velocity, 123(5) km/s. Together with previously published spectroscopic measurements of the systemic radial velocity, this allows a full 3-D determination of the system's velocity; b) precise measurements of the rate of advance of periastron, which after subtraction of the contribution of the proper motion yields a total system mass of $1.6518^{+0.0033}_{-0.0035}$ solar masses; c) a Shapiro delay measurement, h_3 = $82 \pm 14$ ns despite the orbital inclination not being near 90 deg; combined with the measurement of the total mass, this yields a pulsar mass of $1.353^{+0.014}_{-0.017}$ solar masses and a companion mass of $0.298^{+0.015}_{-0.012}$ solar masses; d) we measure precisely the secular variation of the projected semi-major axis and detect significant annual orbital parallax; together these allow a determination of the full 3-D orbital geometry, including an unambiguous orbital inclination (i = $138.7^{+2.5}_{-2.2}$ deg) and a position angle for the line of nodes (Omega = $44^{+5}_{-4}$ deg). We discuss the component masses to investigate hypotheses previously advanced to explain the origin of eccentric MSPs. The unprecedented determination of the full 3-D position, motion and orbital orientation of the system, plus the precisely measured pulsar and WD mass and the latter's optical detection make this system an unique test of our understanding of white dwarfs and their atmospheres. △ Less

Submitted 26 November, 2018; v1 submitted 13 September, 2018; originally announced September 2018.

Comments: 16 pages, 5 figures, accepted by ApJ

Combined collider constraints on neutralinos and charginos

Authors: The GAMBIT Collaboration, Peter Athron, Csaba Balázs, Andy Buckley, Jonathan M. Cornell, Matthias Danninger, Ben Farmer, Andrew Fowlie, Tomás E. Gonzalo, Julia Harz, Paul Jackson, Rose Kudzman-Blais, Anders Kvellestad, Gregory D. Martinez, Andreas Petridis, Are Raklev, Christopher Rogan, Pat Scott, Abhishek Sharma, Martin White, Yang Zhang

Abstract: Searches for supersymmetric electroweakinos have entered a crucial phase, as the integrated luminosity of the Large Hadron Collider is now high enough to compensate for their weak production cross-sections. Working in a framework where the neutralinos and charginos are the only light sparticles in the Minimal Supersymmetric Standard Model, we use gambit to perform a detailed likelihood analysis of… ▽ More Searches for supersymmetric electroweakinos have entered a crucial phase, as the integrated luminosity of the Large Hadron Collider is now high enough to compensate for their weak production cross-sections. Working in a framework where the neutralinos and charginos are the only light sparticles in the Minimal Supersymmetric Standard Model, we use gambit to perform a detailed likelihood analysis of the electroweakino sector. We focus on the impacts of recent ATLAS and CMS searches with 36 fb$^{-1}$ of 13 TeV proton-proton collision data. We also include constraints from LEP and invisible decays of the $Z$ and Higgs bosons. Under the background-only hypothesis, we show that current LHC searches do not robustly exclude any range of neutralino or chargino masses. However, a pattern of excesses in several LHC analyses points towards a possible signal, with neutralino masses of $(m_{\tildeχ_1^0}, m_{\tildeχ_2^0}, m_{\tildeχ_3^0}, m_{\tildeχ_4^0})$ = (8-155, 103-260, 130-473, 219-502) GeV and chargino masses of $(m_{\tildeχ_1^{\pm}}, m_{\tildeχ_2^{\pm}})$ = (104-259, 224-507) GeV at the 95% confidence level. The lightest neutralino is mostly bino, with a possible modest Higgsino or wino component. We find that this excess has a combined local significance of $3.3σ$, subject to a number of cautions. If one includes LHC searches for charginos and neutralinos conducted with 8 TeV proton-proton collision data, the local significance is lowered to 2.9$σ$. We briefly consider the implications for dark matter, finding that the correct relic density can be obtained through the Higgs-funnel and $Z$-funnel mechanisms, even assuming that all other sparticles are decoupled. All samples, gambit input files and best-fit models from this study are available on Zenodo. △ Less

Submitted 20 June, 2019; v1 submitted 6 September, 2018; originally announced September 2018.

Comments: 38 pages, 16 figures, v3 is the version accepted by EPJC

Report number: CoEPP-MN-18-7, gambit-physics-2018

Journal ref: Eur.Phys.J. C79 (2019) no.5, 395

Global analyses of Higgs portal singlet dark matter models using GAMBIT

Authors: The GAMBIT Collaboration, Peter Athron, Csaba Balázs, Ankit Beniwal, Sanjay Bloor, José Eliel Camargo-Molina, Jonathan M. Cornell, Ben Farmer, Andrew Fowlie, Tomás E. Gonzalo, Felix Kahlhoefer, Anders Kvellestad, Gregory D. Martinez, Pat Scott, Aaron C. Vincent, Sebastian Wild, Martin White, Anthony G. Williams

Abstract: We present global analyses of effective Higgs portal dark matter models in the frequentist and Bayesian statistical frameworks. Complementing earlier studies of the scalar Higgs portal, we use GAMBIT to determine the preferred mass and coupling ranges for models with vector, Majorana and Dirac fermion dark matter. We also assess the relative plausibility of all four models using Bayesian model com… ▽ More We present global analyses of effective Higgs portal dark matter models in the frequentist and Bayesian statistical frameworks. Complementing earlier studies of the scalar Higgs portal, we use GAMBIT to determine the preferred mass and coupling ranges for models with vector, Majorana and Dirac fermion dark matter. We also assess the relative plausibility of all four models using Bayesian model comparison. Our analysis includes up-to-date likelihood functions for the dark matter relic density, invisible Higgs decays, and direct and indirect searches for weakly-interacting dark matter including the latest XENON1T data. We also account for important uncertainties arising from the local density and velocity distribution of dark matter, nuclear matrix elements relevant to direct detection, and Standard Model masses and couplings. In all Higgs portal models, we find parameter regions that can explain all of dark matter and give a good fit to all data. The case of vector dark matter requires the most tuning and is therefore slightly disfavoured from a Bayesian point of view. In the case of fermionic dark matter, we find a strong preference for including a CP-violating phase that allows suppression of constraints from direct detection experiments, with odds in favour of CP violation of the order of 100:1. Finally, we present DDCalc 2.0.0, a tool for calculating direct detection observables and likelihoods for arbitrary non-relativistic effective operators. △ Less

Submitted 31 July, 2020; v1 submitted 30 August, 2018; originally announced August 2018.

Comments: 30 pages, 12 figures. v2 matches the published version

Report number: ADP-18-22/T1070, COEPP-MN-18-6, DESY-18-141, TTK-18-34, gambit-physics-2018

Journal ref: Eur.Phys.J. C79 (2019) no.1, 38

Variability Timescale and Spectral Index of Sgr A* in the Near Infrared: Approximate Bayesian Computation Analysis of the Variability of the Closest Supermassive Black Hole

Authors: G. Witzel, G. Martinez, J. Hora, S. P. Willner, M. R. Morris, C. Gammie, E. E. Becklin, M. L. N. Ashby, F. Baganoff, S. Carey, T. Do, G. G. Fazio, A. Ghez, W. J. Glaccum, D. Haggard, R. Herrero-Illana, J. Ingalls, R. Narayan, H. A. Smith

Abstract: Sagittarius A* (Sgr A*) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We present an analysis of the most comprehensive NIR variability dataset of Sgr A* to date: eight 24-hour epochs of continuous monitoring of Sgr A* at 4.5 $μ$m with the IRAC instrument on the Spitzer Space Telescope, 93 epochs of 2.18 $μ$m data from Na… ▽ More Sagittarius A* (Sgr A*) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We present an analysis of the most comprehensive NIR variability dataset of Sgr A* to date: eight 24-hour epochs of continuous monitoring of Sgr A* at 4.5 $μ$m with the IRAC instrument on the Spitzer Space Telescope, 93 epochs of 2.18 $μ$m data from Naos Conica at the Very Large Telescope, and 30 epochs of 2.12 $μ$m data from the NIRC2 camera at the Keck Observatory, in total 94,929 measurements. A new approximate Bayesian computation method for fitting the first-order structure function extracts information beyond current Fast Fourier Transformation (FFT) methods of power spectral density (PSD) estimation. With a combined fit of the data of all three observatories, the characteristic coherence timescale of Sgr A* is $τ_{b} = 243^{+82}_{-57}$ minutes ($90\%$ credible interval). The PSD has no detectable features on timescales down to 8.5 minutes ($95\%$ credible level), which is the ISCO orbital frequency for a dimensionless spin parameter $a = 0.92$. One light curve measured simultaneously at 2.12 and 4.5 $μ$m during a low flux-density phase gave a spectral index $α_s = 1.6 \pm 0.1$ ($F_ν\propto ν^{-α_s}$). This value implies that the Sgr A* NIR color becomes bluer during higher flux-density phases. The probability densities of flux densities of the combined datasets are best fit by log-normal distributions. Based on these distributions, the Sgr A* spectral energy distribution is consistent with synchrotron radiation from a non-thermal electron population from below 20 GHz through the NIR. △ Less

Submitted 5 September, 2018; v1 submitted 1 June, 2018; originally announced June 2018.

Comments: Accepted for publication in ApJ on May 30, 2018. A machine readable version of the light curve data is included in the journal's online publication. Version 2 includes proof corrections

The Green Bank North Celestial Cap Pulsar Survey III: 45 New Pulsar Timing Solutions

Authors: Ryan S. Lynch, Joseph K. Swiggum, Vlad I. Kondratiev, David L. Kaplan, Kevin Stovall, Emmanuel Fonseca, Mallory S. E. Roberts, Lina Levin, Megan E. DeCesar, Bingyi Cui, S. Bradley Cenko, Pradip Gatkine, Anne M. Archibald, Shawn Banaszak, Christopher M. Biwer, Jason Boyles, Pragya Chawla, Louis P. Dartez, David Day, Anthony J. Ford, Joseph Flanigan, Jason W. T. Hessels, Jesus Hinojosa, Fredrick A. Jenet, Chen Karako-Argaman , et al. (15 additional authors not shown)

Abstract: We provide timing solutions for 45 radio pulsars discovered by the Robert C. Byrd Green Bank Telescope. These pulsars were found in the Green Bank North Celestial Cap pulsar survey, an all-GBT-sky survey being carried out at a frequency of 350 MHz. We include pulsar timing data from the Green Bank Telescope and Low Frequency Array. Our sample includes five fully recycled millisecond pulsars (MSPs,… ▽ More We provide timing solutions for 45 radio pulsars discovered by the Robert C. Byrd Green Bank Telescope. These pulsars were found in the Green Bank North Celestial Cap pulsar survey, an all-GBT-sky survey being carried out at a frequency of 350 MHz. We include pulsar timing data from the Green Bank Telescope and Low Frequency Array. Our sample includes five fully recycled millisecond pulsars (MSPs, three of which are in a binary system), a new relativistic double neutron star system, an intermediate mass binary pulsar, a mode-changing pulsar, a 138-ms pulsar with a very low magnetic field, and several nulling pulsars. We have measured two post-Keplerian parameters and thus the masses of both objects in the double neutron star system. We also report a tentative companion mass measurement via Shapiro delay in a binary MSP. Two of the MSPs can be timed with high precision and have been included in pulsar timing arrays being used to search for low-frequency gravitational waves, while a third MSP is a member of the black widow class of binaries. Proper motion is measurable in five pulsars and we provide an estimate of their space velocity. We report on an optical counterpart to a new black widow system and provide constraints on the optical counterparts to other binary MSPs. We also present a preliminary analysis of nulling pulsars in our sample. These results demonstrate the scientific return of long timing campaigns on pulsars of all types. △ Less

Submitted 13 May, 2018; originally announced May 2018.

Comments: Accepted for publication in The Astrophysical Journal; 30 pages, 9 figures, 8 tables

Journal ref: Astrophysical Journal, 2018, 859, 93

The Green Bank Northern Celestial Cap Pulsar Survey II: The Discovery and Timing of Ten Pulsars

Authors: A. M. Kawash, M. A. McLaughlin, D. L. Kaplan, M. E. DeCesar, L. Levin, D. R. Lorimer, R. S. Lynch, K. Stovall, J. K. Swiggum, E. Fonseca, A. M. Archibald, S. Banaszak, C. M. Biwer, J. Boyles, B. Cui, L. P. Dartez, D. Day, S. Ernst, A. J. Ford, J. Flanigan, S. A. Heatherly, J. W. T. Hessels, J. Hinojosa, F. A. Jenet, C. Karako-Argaman , et al. (19 additional authors not shown)

Abstract: We present timing solutions for ten pulsars discovered in 350 MHz searches with the Green Bank Telescope. Nine of these were discovered in the Green Bank Northern Celestial Cap survey and one was discovered by students in the Pulsar Search Collaboratory program in analysis of drift-scan data. Following discovery and confirmation with the Green Bank Telescope, timing has yielded phase-connected sol… ▽ More We present timing solutions for ten pulsars discovered in 350 MHz searches with the Green Bank Telescope. Nine of these were discovered in the Green Bank Northern Celestial Cap survey and one was discovered by students in the Pulsar Search Collaboratory program in analysis of drift-scan data. Following discovery and confirmation with the Green Bank Telescope, timing has yielded phase-connected solutions with high precision measurements of rotational and astrometric parameters. Eight of the pulsars are slow and isolated, including PSR J0930$-$2301, a pulsar with nulling fraction lower limit of $\sim$30\% and nulling timescale of seconds to minutes. This pulsar also shows evidence of mode changing. The remaining two pulsars have undergone recycling, accreting material from binary companions, resulting in higher spin frequencies. PSR J0557$-$2948 is an isolated, 44 \rm{ms} pulsar that has been partially recycled and is likely a former member of a binary system which was disrupted by a second supernova. The paucity of such so-called `disrupted binary pulsars' (DRPs) compared to double neutron star (DNS) binaries can be used to test current evolutionary scenarios, especially the kicks imparted on the neutron stars in the second supernova. There is some evidence that DRPs have larger space velocities, which could explain their small numbers. PSR J1806+2819 is a 15 \rm{ms} pulsar in a 44 day orbit with a low mass white dwarf companion. We did not detect the companion in archival optical data, indicating that it must be older than 1200 Myr. △ Less

Submitted 9 March, 2018; originally announced March 2018.

Comments: 9 pages, 5 figures

Pulsar J1411+2551: A Low Mass New Double Neutron Star System

Authors: J. G. Martinez, K. Stovall, P. C. C. Freire, J. S. Deneva, T. M. Tauris, A. Ridolfi, N. Wex, F. A. Jenet, M. A. McLaughlin, M. Bagchi

Abstract: In this work, we report the discovery and characterization of PSR J1411+2551, a new binary pulsar discovered in the Arecibo 327 MHz Drift Pulsar Survey. Our timing observations of the radio pulsar in the system span a period of about 2.5 years. This timing campaign allowed a precise measurement of its spin period (62.4 ms) and its derivative (9.6 $\pm$ 0.7) $\times 10^{-20}\, \rm s\, s^{-1}$; from… ▽ More In this work, we report the discovery and characterization of PSR J1411+2551, a new binary pulsar discovered in the Arecibo 327 MHz Drift Pulsar Survey. Our timing observations of the radio pulsar in the system span a period of about 2.5 years. This timing campaign allowed a precise measurement of its spin period (62.4 ms) and its derivative (9.6 $\pm$ 0.7) $\times 10^{-20}\, \rm s\, s^{-1}$; from these, we derive a characteristic age of $\sim 10\,$Gyr and a surface magnetic field strength of 2.5 $\times 10^{9}$ G. These numbers indicate that this pulsar was mildly recycled by accretion of matter from the progenitor of the companion star. The system has an eccentric ($e\, = \, 0.17$) 2.61 day orbit. This eccentricity allows a highly significant measurement of the rate of advance of periastron, $\dotω = 0.07686 \pm 0.00046 ^{\circ}~{\rm yr}^{-1}$. Assuming general relativity accurately models the orbital motion, this implies a total system mass M = $2.538 \pm 0.022 M_{\odot}$. The minimum companion mass is $0.92\, M_{\odot}$ and the maximum pulsar mass is $1.62\, M_{\odot}$. The large companion mass and the orbital eccentricity suggest that PSR J1411+2551 is a double neutron star system; the lightest known to date including the DNS merger GW 170817. Furthermore, the relatively low orbital eccentricity and small proper motion limits suggest that the second supernova had a relatively small associated kick; this and the low system mass suggest that it was an ultra-stripped supernova. △ Less

Submitted 7 December, 2017; v1 submitted 27 November, 2017; originally announced November 2017.

Comments: Accepted for publication in APJ letters

Retrieval of Water Vapor Column Abundance and Aerosol Properties from ChemCam Passive Sky Spectroscopy

Authors: Timothy H. McConnochie, Michael D. Smith, Michael J. Wolff, Steve Bender, Mark Lemmon, Roger C. Wiens, Sylvestre Maurice, Olivier Gasnault, Jeremie Lasue, Pierre-Yves Meslin, Ari-Matti Harri, Maria Genzer, Osku Kemppinen, Germán M. Martínez, Lauren DeFlores, Diana Blaney, Jeffrey R. Johnson, James F. Bell III

Abstract: We derive water vapor column abundances and aerosol properties from Mars Science Laboratory (MSL) ChemCam passive mode observations of scattered sky light. Each ChemCam passive sky observation acquires spectra at two different elevation angles. We fit these spectra with a discrete-ordinates multiple scattering radiative transfer model, using the correlated-k approximation for gas absorption band… ▽ More We derive water vapor column abundances and aerosol properties from Mars Science Laboratory (MSL) ChemCam passive mode observations of scattered sky light. Each ChemCam passive sky observation acquires spectra at two different elevation angles. We fit these spectra with a discrete-ordinates multiple scattering radiative transfer model, using the correlated-k approximation for gas absorption bands. The retrieval proceeds by first fitting the continuum of the ratio of the two elevation angles to solve for aerosol properties, and then fitting the continuum-removed ratio to solve for gas abundances. The final step of the retrieval makes use of the observed CO2 absorptions and the known CO2 abundance to correct the retrieved water vapor abundance for the effects of the vertical distribution of scattering aerosols and to derive an aerosol scale height parameter. The ChemCam-retrieved water abundances show, with only a few exceptions, the same seasonal behavior and the same timing of seasonal minima and maxima as the TES, CRISM, and REMS-H data sets that we compare them to. However ChemCam-retrieved water abundances are generally lower than zonal and regional scale from-orbit water vapor data, while at the same time being significantly larger than pre-dawn REMS-H abundances. Pending further analysis of REMS-H volume mixing ratio uncertainties, the differences between ChemCam and REMS-H pre-dawn mixing ratios appear to be much too large to be explained by large scale circulations and thus they tend to support the hypothesis of substantial diurnal interactions of water vapor with the surface. Our preliminary aerosol results, meanwhile, show the expected seasonal pattern in dust particle size but also indicate a surprising inter-annual increase in water-ice cloud opacities. △ Less

Submitted 29 October, 2017; v1 submitted 24 October, 2017; originally announced October 2017.

Comments: 64 pages with embedded figures; this is the accepted version of the manuscript; the meta-data version of the abstract has been shorted to meet arXiv rules