Showing 1–50 of 105 results for author: Chang, P

Dual-comb correlation spectroscopy of thermal light

Authors: Eugene J. Tsao, Alexander J. Lind, Connor Fredrick, Ryan K. Cole, Peter Chang, Kristina F. Chang, Dahyeon Lee, Matthew Heyrich, Nazanin Hoghooghi, Franklyn Quinlan, Scott A. Diddams

Abstract: The detection of light of thermal origin is the principal means by which humanity has learned about our world and the cosmos. In optical astronomy, in particular, direct detection of thermal photons and the resolution of their spectra have enabled discoveries of the broadest scope and impact. Such measurements, however, do not capture the phase of the thermal fields--a parameter that has proven cr… ▽ More The detection of light of thermal origin is the principal means by which humanity has learned about our world and the cosmos. In optical astronomy, in particular, direct detection of thermal photons and the resolution of their spectra have enabled discoveries of the broadest scope and impact. Such measurements, however, do not capture the phase of the thermal fields--a parameter that has proven crucial to transformative techniques in radio astronomy such as synthetic aperture imaging. Over the last 25 years, tremendous progress has occurred in laser science, notably in the phase-sensitive, broad bandwidth, high resolution, and traceable spectroscopy enabled by the optical frequency comb. In this work, we directly connect the fields of frequency comb laser spectroscopy and passive optical sensing as applied to astronomy, remote sensing, and atmospheric science. We provide fundamental sensitivity analysis of dual-comb correlation spectroscopy (DCCS), whereby broadband thermal light is measured via interferometry with two optical frequency combs. We define and experimentally verify the sensitivity scaling of DCCS at black body temperatures relevant for astrophysical observations. Moreover, we provide comparison with direct detection techniques and more conventional laser heterodyne radiometry. Our work provides the foundation for future exploration of comb-based broadband synthetic aperture hyperspectral imaging across the infrared and optical spectrum. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: 54 pages, 7 figures

Galactic Structure From Binary Pulsar Accelerations: Beyond Smooth Models

Authors: Thomas Donlon II, Sukanya Chakrabarti, Lawrence M. Widrow, Michael T. Lam, Philip Chang, Alice C. Quillen

Abstract: We measure the line-of-sight accelerations of 26 binary pulsars due to the Milky Way's gravitational potential, and produce a 3-dimensional map of the acceleration field of the Galaxy. Acceleration measurements directly give us the change in the line-of-sight velocity at present day, without requiring any assumptions inherent to kinematic modeling. We measure the Oort limit ($ρ_0=0.062\pm0.017$ \m… ▽ More We measure the line-of-sight accelerations of 26 binary pulsars due to the Milky Way's gravitational potential, and produce a 3-dimensional map of the acceleration field of the Galaxy. Acceleration measurements directly give us the change in the line-of-sight velocity at present day, without requiring any assumptions inherent to kinematic modeling. We measure the Oort limit ($ρ_0=0.062\pm0.017$ \msun/pc$^3$) and the dark matter density in the midplane ($ρ_{0,\textrm{DM}}=-0.010\pm0.018$ \msun/pc$^3$); these values are similar to, but have smaller uncertainties than previous pulsar timing measurements of these quantities. Here, we provide for the first time, values for the Oort constants and the slope of the rotation curve from direct acceleration measurements. We find that $A=15.4\pm2.6$ km/s/kpc and $B=-13.1\pm2.6$ km/s/kpc (consistent with results from \textit{Gaia}), and the slope of the rotation curve near the Sun is $-2\pm5$ km/s/kpc. We show that the Galactic acceleration field is clearly asymmetric, but due to data limitations it is not yet clear which physical processes drive this asymmetry. We provide updated models of the Galactic potential that account for various sources of disequilibrium; these models are incompatible with commonly used kinematic potentials. This indicates that use of kinematically derived Galactic potentials in precision tests (e.g., in tests of general relativity with pulsar timing) may be subject to larger uncertainties than reported. The acceleration data indicates that the mass of the Galaxy within the Solar circle is $2.3 \times 10^{11}$ M$_\odot$, roughly twice as large as currently accepted models. Additionally, the residuals of the acceleration data compared to existing Galactic models have a dependence on radial position; this trend can be explained if the Sun has an additional acceleration away from the Galactic center. △ Less

Submitted 20 May, 2024; v1 submitted 28 January, 2024; originally announced January 2024.

LenSiam: Self-Supervised Learning on Strong Gravitational Lens Images

Authors: Po-Wen Chang, Kuan-Wei Huang, Joshua Fagin, James Hung-Hsu Chan, Joshua Yao-Yu Lin

Abstract: Self-supervised learning has been known for learning good representations from data without the need for annotated labels. We explore the simple siamese (SimSiam) architecture for representation learning on strong gravitational lens images. Commonly used image augmentations tend to change lens properties; for example, zoom-in would affect the Einstein radius. To create image pairs representing the… ▽ More Self-supervised learning has been known for learning good representations from data without the need for annotated labels. We explore the simple siamese (SimSiam) architecture for representation learning on strong gravitational lens images. Commonly used image augmentations tend to change lens properties; for example, zoom-in would affect the Einstein radius. To create image pairs representing the same underlying lens model, we introduce a lens augmentation method to preserve lens properties by fixing the lens model while varying the source galaxies. Our research demonstrates this lens augmentation works well with SimSiam for learning the lens image representation without labels, so we name it LenSiam. We also show that a pre-trained LenSiam model can benefit downstream tasks. We open-source our code and datasets at https://github.com/kuanweih/LenSiam . △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: 5 pages, 2 figures. Accepted by NeurIPS 2023 AI for Science Workshop

Envelope Ejection and the Transition to Homologous Expansion in Common-Envelope Events

Authors: Vinaya Valsan, Sarah V. Borges, Logan Prust, Philip Chang

Abstract: We conduct a long-timescale ($5000\,$d) 3-D simulation of a common-envelope event with a $2\,M_{\odot}$ red giant and a $1\,M_{\odot}$ main sequence companion, using the moving-mesh hydrodynamic solver MANGA. Starting with an orbital radius of $52\,R_{\odot}$, our binary shrinks to an orbital radius of $5\,R_{\odot}$ in $200\,$d. We show that over a timescale of about $1500\,$d, the envelope is co… ▽ More We conduct a long-timescale ($5000\,$d) 3-D simulation of a common-envelope event with a $2\,M_{\odot}$ red giant and a $1\,M_{\odot}$ main sequence companion, using the moving-mesh hydrodynamic solver MANGA. Starting with an orbital radius of $52\,R_{\odot}$, our binary shrinks to an orbital radius of $5\,R_{\odot}$ in $200\,$d. We show that over a timescale of about $1500\,$d, the envelope is completely ejected while $80$ per cent is ejected in about $400\,$d. The complete ejection of the envelope is solely powered by the orbital energy of the binary, without the need for late-time reheating from recombination or jets. Motivated by recent theoretical and observational results, we also find that the envelope enters a phase of homologous expansion about $550\,\rm d$ after the start of our simulation. We also run a simplified 1-D model to show that heating from the central binary in the envelope at late times does not influence the ejection. This homologous expansion of the envelope would likely simplify calculations of the observational implications such as light curves. △ Less

Submitted 27 September, 2023; originally announced September 2023.

A Joint Fermi-GBM and Swift-BAT Analysis of Gravitational-Wave Candidates from the Third Gravitational-wave Observing Run

Authors: C. Fletcher, J. Wood, R. Hamburg, P. Veres, C. M. Hui, E. Bissaldi, M. S. Briggs, E. Burns, W. H. Cleveland, M. M. Giles, A. Goldstein, B. A. Hristov, D. Kocevski, S. Lesage, B. Mailyan, C. Malacaria, S. Poolakkil, A. von Kienlin, C. A. Wilson-Hodge, The Fermi Gamma-ray Burst Monitor Team, M. Crnogorčević, J. DeLaunay, A. Tohuvavohu, R. Caputo, S. B. Cenko , et al. (1674 additional authors not shown)

Abstract: We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses,… ▽ More We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses, the Targeted Search and the Untargeted Search, we investigate whether there are any coincident GRBs associated with the GWs. We also search the Swift-BAT rate data around the GW times to determine whether a GRB counterpart is present. No counterparts are found. Using both the Fermi-GBM Targeted Search and the Swift-BAT search, we calculate flux upper limits and present joint upper limits on the gamma-ray luminosity of each GW. Given these limits, we constrain theoretical models for the emission of gamma-rays from binary black hole mergers. △ Less

Submitted 25 August, 2023; originally announced August 2023.

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

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

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

Submitted 7 August, 2023; originally announced August 2023.

Comments: 24 pages, 5 figures

Report number: LIGO-P2300080

Search for gravitational-lensing signatures in the full third observing run of the LIGO-Virgo network

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

Abstract: Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated… ▽ More Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated signals from strong lensing by 1) performing targeted searches for subthreshold signals, 2) calculating the degree of overlap amongst the intrinsic parameters and sky location of pairs of signals, 3) comparing the similarities of the spectrograms amongst pairs of signals, and 4) performing dual-signal Bayesian analysis that takes into account selection effects and astrophysical knowledge. We also search for distortions to the gravitational waveform caused by 1) frequency-independent phase shifts in strongly lensed images, and 2) frequency-dependent modulation of the amplitude and phase due to point masses. None of these searches yields significant evidence for lensing. Finally, we use the non-detection of gravitational-wave lensing to constrain the lensing rate based on the latest merger-rate estimates and the fraction of dark matter composed of compact objects. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 28 pages, 11 figures

Report number: LIGO-P2200031

Search for subsolar-mass black hole binaries in the second part of Advanced LIGO's and Advanced Virgo's third observing run

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

Abstract: We describe a search for gravitational waves from compact binaries with at least one component with mass 0.2 $M_\odot$ -- $1.0 M_\odot$ and mass ratio $q \geq 0.1$ in Advanced LIGO and Advanced Virgo data collected between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. No signals were detected. The most significant candidate has a false alarm rate of 0.2 $\mathrm{yr}^{-1}$. We estimate t… ▽ More We describe a search for gravitational waves from compact binaries with at least one component with mass 0.2 $M_\odot$ -- $1.0 M_\odot$ and mass ratio $q \geq 0.1$ in Advanced LIGO and Advanced Virgo data collected between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. No signals were detected. The most significant candidate has a false alarm rate of 0.2 $\mathrm{yr}^{-1}$. We estimate the sensitivity of our search over the entirety of Advanced LIGO's and Advanced Virgo's third observing run, and present the most stringent limits to date on the merger rate of binary black holes with at least one subsolar-mass component. We use the upper limits to constrain two fiducial scenarios that could produce subsolar-mass black holes: primordial black holes (PBH) and a model of dissipative dark matter. The PBH model uses recent prescriptions for the merger rate of PBH binaries that include a rate suppression factor to effectively account for PBH early binary disruptions. If the PBHs are monochromatically distributed, we can exclude a dark matter fraction in PBHs $f_\mathrm{PBH} \gtrsim 0.6$ (at 90% confidence) in the probed subsolar-mass range. However, if we allow for broad PBH mass distributions we are unable to rule out $f_\mathrm{PBH} = 1$. For the dissipative model, where the dark matter has chemistry that allows a small fraction to cool and collapse into black holes, we find an upper bound $f_{\mathrm{DBH}} < 10^{-5}$ on the fraction of atomic dark matter collapsed into black holes. △ Less

Submitted 26 January, 2024; v1 submitted 2 December, 2022; originally announced December 2022.

Comments: https://dcc.ligo.org/P2200139

Search for gravitational-wave transients associated with magnetar bursts in Advanced LIGO and Advanced Virgo data from the third observing run

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

Abstract: Gravitational waves are expected to be produced from neutron star oscillations associated with magnetar giant flares and short bursts. We present the results of a search for short-duration (milliseconds to seconds) and long-duration ($\sim$ 100 s) transient gravitational waves from 13 magnetar short bursts observed during Advanced LIGO, Advanced Virgo and KAGRA's third observation run. These 13 bu… ▽ More Gravitational waves are expected to be produced from neutron star oscillations associated with magnetar giant flares and short bursts. We present the results of a search for short-duration (milliseconds to seconds) and long-duration ($\sim$ 100 s) transient gravitational waves from 13 magnetar short bursts observed during Advanced LIGO, Advanced Virgo and KAGRA's third observation run. These 13 bursts come from two magnetars, SGR 1935$+$2154 and Swift J1818.0$-$1607. We also include three other electromagnetic burst events detected by Fermi GBM which were identified as likely coming from one or more magnetars, but they have no association with a known magnetar. No magnetar giant flares were detected during the analysis period. We find no evidence of gravitational waves associated with any of these 16 bursts. We place upper bounds on the root-sum-square of the integrated gravitational-wave strain that reach $2.2 \times 10^{-23}$ $/\sqrt{\text{Hz}}$ at 100 Hz for the short-duration search and $8.7 \times 10^{-23}$ $/\sqrt{\text{Hz}}$ at $450$ Hz for the long-duration search, given a detection efficiency of 50%. For a ringdown signal at 1590 Hz targeted by the short-duration search the limit is set to $1.8 \times 10^{-22}$ $/\sqrt{\text{Hz}}$. Using the estimated distance to each magnetar, we derive upper bounds on the emitted gravitational-wave energy of $3.2 \times 10^{43}$ erg ($7.3 \times 10^{43}$ erg) for SGR 1935$+$2154 and $8.2 \times 10^{42}$ erg ($2.8 \times 10^{43}$ erg) for Swift J1818.0$-$1607, for the short-duration (long-duration) search. Assuming isotropic emission of electromagnetic radiation of the burst fluences, we constrain the ratio of gravitational-wave energy to electromagnetic energy for bursts from SGR 1935$+$2154 with available fluence information. The lowest of these ratios is $3 \times 10^3$. △ Less

Submitted 19 October, 2022; originally announced October 2022.

Comments: 30 pages with appendices, 5 figures, 10 tables

Report number: LIGO-P2100387

A non-interacting Galactic black hole candidate in a binary system with a main-sequence star

Authors: Sukanya Chakrabarti, Joshua D. Simon, Peter A. Craig, Henrique Reggiani, Timothy D. Brandt, Puragra Guhathakurta, Paul A. Dalba, Evan N. Kirby, Philip Chang, Daniel R. Hey, Alessandro Savino, Marla Geha, Ian B. Thompson

Abstract: We describe the discovery of a solar neighborhood (d=468 pc) binary system with a main-sequence sunlike star and a massive non-interacting black hole candidate. The spectral energy distribution (SED) of the visible star is described by a single stellar model. We derive stellar parameters from a high signal-to-noise Magellan/MIKE spectrum, classifying the star as a main-sequence star with… ▽ More We describe the discovery of a solar neighborhood (d=468 pc) binary system with a main-sequence sunlike star and a massive non-interacting black hole candidate. The spectral energy distribution (SED) of the visible star is described by a single stellar model. We derive stellar parameters from a high signal-to-noise Magellan/MIKE spectrum, classifying the star as a main-sequence star with $T_{\rm eff} = 5972 \rm K$, $\log{g} = 4.54$, and $M = 0.91$ \msun. The spectrum shows no indication of a second luminous component. To determine the spectroscopic orbit of the binary, we measured radial velocities of this system with the Automated Planet Finder, Magellan, and Keck over four months. We show that the velocity data are consistent with the \textit{Gaia} astrometric orbit and provide independent evidence for a massive dark companion. From a combined fit of our spectroscopic data and the astrometry, we derive a companion mass of $11.39^{+1.51}_{-1.31}$\msun. We conclude that this binary system harbors a massive black hole on an eccentric $(e =0.46 \pm 0.02)$, $185.4 \pm 0.1$ d orbit. These conclusions are independent of \cite{ElBadry2022Disc}, who recently reported the discovery of the same system. A joint fit to all available data (including \cite{ElBadry2022Disc}'s) yields a comparable period solution, but a lower companion mass of $9.32^{+0.22}_{-0.21} M_{\odot}$. Radial velocity fits to all available data produce a unimodal solution for the period that is not possible with either data set alone. The combination of both data sets yields the most accurate orbit currently available. △ Less

Submitted 5 June, 2023; v1 submitted 10 October, 2022; originally announced October 2022.

Comments: 30 pages, 17 figures, accepted to AAS journals (various changes in response to referee comments)

Strong Gravitational Lensing Parameter Estimation with Vision Transformer

Authors: Kuan-Wei Huang, Geoff Chih-Fan Chen, Po-Wen Chang, Sheng-Chieh Lin, Chia-Jung Hsu, Vishal Thengane, Joshua Yao-Yu Lin

Abstract: Quantifying the parameters and corresponding uncertainties of hundreds of strongly lensed quasar systems holds the key to resolving one of the most important scientific questions: the Hubble constant ($H_{0}$) tension. The commonly used Markov chain Monte Carlo (MCMC) method has been too time-consuming to achieve this goal, yet recent work has shown that convolution neural networks (CNNs) can be a… ▽ More Quantifying the parameters and corresponding uncertainties of hundreds of strongly lensed quasar systems holds the key to resolving one of the most important scientific questions: the Hubble constant ($H_{0}$) tension. The commonly used Markov chain Monte Carlo (MCMC) method has been too time-consuming to achieve this goal, yet recent work has shown that convolution neural networks (CNNs) can be an alternative with seven orders of magnitude improvement in speed. With 31,200 simulated strongly lensed quasar images, we explore the usage of Vision Transformer (ViT) for simulated strong gravitational lensing for the first time. We show that ViT could reach competitive results compared with CNNs, and is specifically good at some lensing parameters, including the most important mass-related parameters such as the center of lens $θ_{1}$ and $θ_{2}$, the ellipticities $e_1$ and $e_2$, and the radial power-law slope $γ'$. With this promising preliminary result, we believe the ViT (or attention-based) network architecture can be an important tool for strong lensing science for the next generation of surveys. The open source of our code and data is in \url{https://github.com/kuanweih/strong_lensing_vit_resnet}. △ Less

Submitted 8 October, 2022; originally announced October 2022.

Comments: Accepted by ECCV 2022 AI for Space Workshop

High-energy neutrinos from choked-jet supernovae: Searches and implications

Authors: Po-Wen Chang, Bei Zhou, Kohta Murase, Marc Kamionkowski

Abstract: The origin of the high-energy astrophysical neutrinos discovered by IceCube remains largely unknown. Multimessenger studies have indicated that the majority of these neutrinos come from gamma-ray-dark sources. Choked-jet supernovae (cjSNe), which are supernovae powered by relativistic jets stalled in stellar materials, may lead to neutrino emission via photohadronic interactions while the coproduc… ▽ More The origin of the high-energy astrophysical neutrinos discovered by IceCube remains largely unknown. Multimessenger studies have indicated that the majority of these neutrinos come from gamma-ray-dark sources. Choked-jet supernovae (cjSNe), which are supernovae powered by relativistic jets stalled in stellar materials, may lead to neutrino emission via photohadronic interactions while the coproduced gamma rays are absorbed. In this paper, we perform an unbinned maximum-likelihood analysis to search for correlations between IceCube's ten-year muon-track events and our SN Ib/c sample, collected from publicly available catalogs. In addition to the conventional power-law models, we also consider the impacts of more realistic neutrino emission models for the first time, and we study the effects of the jet beaming factor in the analyses. Our results show no significant correlation. Even so, the conservative upper limits we set to the contribution of cjSNe to the diffuse astrophysical neutrino flux still allow SNe Ib/c to be the dominant source of astrophysical neutrinos observed by IceCube. We discuss implications to the cjSNe scenario from our results and the power of future neutrino and supernova observations. △ Less

Submitted 31 May, 2024; v1 submitted 6 October, 2022; originally announced October 2022.

Comments: 12 pages, 7 figures. Match the published version: Following IceCube's stacking analysis approach, we set new limits w.r.t. TS_obs (see Sec. IV and Fig. 4). As a result, our constraints have improved by a factor of ~2 compared to our previous results (see new results in Secs. V-C, V-D and Figs. 2, 6, and 7)

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

Model-based cross-correlation search for gravitational waves from the low-mass X-ray binary Scorpius X-1 in LIGO O3 data

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

Abstract: We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to bala… ▽ More We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to balance sensitivity with computing cost. The search covered a range of gravitational-wave frequencies from 25Hz to 1600Hz, as well as ranges in orbital speed, frequency and phase determined from observational constraints. No significant detection candidates were found, and upper limits were set as a function of frequency. The most stringent limits, between 100Hz and 200Hz, correspond to an amplitude h0 of about 1e-25 when marginalized isotropically over the unknown inclination angle of the neutron star's rotation axis, or less than 4e-26 assuming the optimal orientation. The sensitivity of this search is now probing amplitudes predicted by models of torque balance equilibrium. For the usual conservative model assuming accretion at the surface of the neutron star, our isotropically-marginalized upper limits are close to the predicted amplitude from about 70Hz to 100Hz; the limits assuming the neutron star spin is aligned with the most likely orbital angular momentum are below the conservative torque balance predictions from 40Hz to 200Hz. Assuming a broader range of accretion models, our direct limits on gravitational-wave amplitude delve into the relevant parameter space over a wide range of frequencies, to 500Hz or more. △ Less

Submitted 2 January, 2023; v1 submitted 6 September, 2022; originally announced September 2022.

Comments: 19 pages, Open Access Journal PDF

Report number: LIGO-P2100110-v13

Journal ref: The Astrophysical Journal Letters, 941, L30 (2022)

Towards Powerful Probes of Neutrino Self-Interactions in Supernovae

Authors: Po-Wen Chang, Ivan Esteban, John F. Beacom, Todd A. Thompson, Christopher M. Hirata

Abstract: Neutrinos remain mysterious. As an example, enhanced self-interactions ($ν$SI), which would have broad implications, are allowed. At the high neutrino densities within core-collapse supernovae, $ν$SI should be important, but robust observables have been lacking. We show that $ν$SI make neutrinos form a tightly coupled fluid that expands under relativistic hydrodynamics. The outflow becomes either… ▽ More Neutrinos remain mysterious. As an example, enhanced self-interactions ($ν$SI), which would have broad implications, are allowed. At the high neutrino densities within core-collapse supernovae, $ν$SI should be important, but robust observables have been lacking. We show that $ν$SI make neutrinos form a tightly coupled fluid that expands under relativistic hydrodynamics. The outflow becomes either a burst or a steady-state wind; which occurs here is uncertain. Though the diffusive environment where neutrinos are produced may make a wind more likely, further work is needed to determine when each case is realized. In the burst-outflow case, $ν$SI increase the duration of the neutrino signal, and even a simple analysis of SN 1987A data has powerful sensitivity. For the wind-outflow case, we outline several promising ideas that may lead to new observables. Combined, these results are important steps towards solving the 35-year-old puzzle of how $ν$SI affect supernovae. △ Less

Submitted 28 June, 2023; v1 submitted 24 June, 2022; originally announced June 2022.

Comments: 5+19 pages, 4+14 figures. Accepted for publication in Phys.Rev.Lett. Match the version to be published: minor revisions for clarity; moved discussion about supernova neutrino properties to Supplemental Material C; updated references

Journal ref: Phys. Rev. Lett. 131, 071002 (2023)

Search for continuous gravitational wave emission from the Milky Way center in O3 LIGO--Virgo data

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

Abstract: We present a directed search for continuous gravitational wave (CW) signals emitted by spinning neutron stars located in the inner parsecs of the Galactic Center (GC). Compelling evidence for the presence of a numerous population of neutron stars has been reported in the literature, turning this region into a very interesting place to look for CWs. In this search, data from the full O3 LIGO--Virgo… ▽ More We present a directed search for continuous gravitational wave (CW) signals emitted by spinning neutron stars located in the inner parsecs of the Galactic Center (GC). Compelling evidence for the presence of a numerous population of neutron stars has been reported in the literature, turning this region into a very interesting place to look for CWs. In this search, data from the full O3 LIGO--Virgo run in the detector frequency band $[10,2000]\rm~Hz$ have been used. No significant detection was found and 95$\%$ confidence level upper limits on the signal strain amplitude were computed, over the full search band, with the deepest limit of about $7.6\times 10^{-26}$ at $\simeq 142\rm~Hz$. These results are significantly more constraining than those reported in previous searches. We use these limits to put constraints on the fiducial neutron star ellipticity and r-mode amplitude. These limits can be also translated into constraints in the black hole mass -- boson mass plane for a hypothetical population of boson clouds around spinning black holes located in the GC. △ Less

Submitted 9 April, 2022; originally announced April 2022.

Comments: 25 pages, 5 figures

Snowmass2021 Cosmic Frontier White Paper: Observational Facilities to Study Dark Matter

Authors: Sukanya Chakrabarti, Alex Drlica-Wagner, Ting S. Li, Neelima Sehgal, Joshua D. Simon, Simon Birrer, Duncan A. Brown, Rebecca Bernstein, Alberto D. Bolatto, Philip Chang, Kyle Dawson, Paul Demorest, Daniel Grin, David L. Kaplan, Joseph Lazio, Jennifer Marshall, Eric J. Murphy, Scott Ransom, Brant E. Robertson, Rajeev Singh, Anže Slosar, Tommaso Treu, Yu-Dai Tsai, Benjamin F. Williams

Abstract: We present an overview of future observational facilities that will significantly enhance our understanding of the fundamental nature of dark matter. These facilities span a range of observational techniques including optical/near-infrared imaging and spectroscopy, measurements of the cosmic microwave background, pulsar timing, 21-cm observations of neutral hydrogen at high redshift, and the measu… ▽ More We present an overview of future observational facilities that will significantly enhance our understanding of the fundamental nature of dark matter. These facilities span a range of observational techniques including optical/near-infrared imaging and spectroscopy, measurements of the cosmic microwave background, pulsar timing, 21-cm observations of neutral hydrogen at high redshift, and the measurement of gravitational waves. Such facilities are a critical component of a multi-pronged experimental program to uncover the nature of dark matter, while often providing complementary measurements of dark energy, neutrino physics, and inflation. △ Less

Submitted 11 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021. 31 pages, 6 figures, 2 tables

First joint observation by the underground gravitational-wave detector, KAGRA, with GEO600

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

Abstract: We report the results of the first joint observation of the KAGRA detector with GEO600. KAGRA is a cryogenic and underground gravitational-wave detector consisting of a laser interferometer with three-kilometer arms, and located in Kamioka, Gifu, Japan. GEO600 is a British--German laser interferometer with 600 m arms, and located near Hannover, Germany. GEO600 and KAGRA performed a joint observing… ▽ More We report the results of the first joint observation of the KAGRA detector with GEO600. KAGRA is a cryogenic and underground gravitational-wave detector consisting of a laser interferometer with three-kilometer arms, and located in Kamioka, Gifu, Japan. GEO600 is a British--German laser interferometer with 600 m arms, and located near Hannover, Germany. GEO600 and KAGRA performed a joint observing run from April 7 to 20, 2020. We present the results of the joint analysis of the GEO--KAGRA data for transient gravitational-wave signals, including the coalescence of neutron-star binaries and generic unmodeled transients. We also perform dedicated searches for binary coalescence signals and generic transients associated with gamma-ray burst events observed during the joint run. No gravitational-wave events were identified. We evaluate the minimum detectable amplitude for various types of transient signals and the spacetime volume for which the network is sensitive to binary neutron-star coalescences. We also place lower limits on the distances to the gamma-ray bursts analysed based on the non-detection of an associated gravitational-wave signal for several signal models, including binary coalescences. These analyses demonstrate the feasibility and utility of KAGRA as a member of the global gravitational-wave detector network. △ Less

Submitted 19 August, 2022; v1 submitted 2 March, 2022; originally announced March 2022.

Comments: Matches with published version

Report number: LIGO-P2100286

Journal ref: Progress of Theoretical and Experimental Physics, Volume 2022, Issue 6, 063F01 (2022)

Constraining blazar heating with the 2<z<3 Lyman-$α$ forest

Authors: A. Lamberts, E. Puchwein, C. Pfrommer, P. Chang, M. Shalaby, A. Broderick, P. Tiede, G. Rudie

Abstract: The intergalactic medium (IGM) acts like a calorimeter recording energy injection by cosmic structure formation, shocks and photoheating from stars and active galactic nuclei. It was recently proposed that spatially inhomogeneous TeV-blazars could significantly heat up the underdense IGM, resulting in patches of both cold and warm IGM around $z\simeq2-3$. The goal of this study is to compare predi… ▽ More The intergalactic medium (IGM) acts like a calorimeter recording energy injection by cosmic structure formation, shocks and photoheating from stars and active galactic nuclei. It was recently proposed that spatially inhomogeneous TeV-blazars could significantly heat up the underdense IGM, resulting in patches of both cold and warm IGM around $z\simeq2-3$. The goal of this study is to compare predictions of different blazar heating models with recent observations of the IGM. We perform a set of cosmological simulations and carefully compute mock observables of the Lyman-$α ($Ly$α$) forest. We perform a detailed assessment of different systematic uncertainties which typically impact this type of observables and find that they are smaller than the differences between our models. We find that our inhomogeneous blazar heating model is in good agreement with the Ly$α$ line properties and the rescaled flux probability distribution function at high redshift ($2.5<z<3$) but that our blazar heating models are challenged by lower redshift data ($2<z<2.5$). Our results could be explained by HeII reionisation although state-of-the-art models fall short on providing enough heating to the low-density IGM, thus motivating further radiative transfer studies of inhomogeneous HeII reionisation. If blazars are indeed hosted by group-mass halos of $2\times10^{13}{M}_\odot$, a later onset of blazar heating in comparison to previous models would be favoured, which could bring our findings here in agreement with the evidence of blazar heating from local gamma-ray observations. △ Less

Submitted 2 February, 2022; v1 submitted 31 January, 2022; originally announced January 2022.

Comments: Resubmitted after minor revisions

All-sky search for continuous gravitational waves from isolated neutron stars using Advanced LIGO and Advanced Virgo O3 data

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

Abstract: We present results of an all-sky search for continuous gravitational waves which can be produced by spinning neutron stars with an asymmetry around their rotation axis, using data from the third observing run of the Advanced LIGO and Advanced Virgo detectors. Four different analysis methods are used to search in a gravitational-wave frequency band from 10 to 2048 Hz and a first frequency derivativ… ▽ More We present results of an all-sky search for continuous gravitational waves which can be produced by spinning neutron stars with an asymmetry around their rotation axis, using data from the third observing run of the Advanced LIGO and Advanced Virgo detectors. Four different analysis methods are used to search in a gravitational-wave frequency band from 10 to 2048 Hz and a first frequency derivative from $-10^{-8}$ to $10^{-9}$ Hz/s. No statistically-significant periodic gravitational-wave signal is observed by any of the four searches. As a result, upper limits on the gravitational-wave strain amplitude $h_0$ are calculated. The best upper limits are obtained in the frequency range of 100 to 200 Hz and they are ${\sim}1.1\times10^{-25}$ at 95\% confidence-level. The minimum upper limit of $1.10\times10^{-25}$ is achieved at a frequency 111.5 Hz. We also place constraints on the rates and abundances of nearby planetary- and asteroid-mass primordial black holes that could give rise to continuous gravitational-wave signals. △ Less

Submitted 3 January, 2022; originally announced January 2022.

Comments: 23 main text pages, 17 figures

Report number: LIGO-P2100367

Eclipse timing the Milky Way's gravitational potential

Authors: Sukanya Chakrabarti, Daniel J. Stevens, Jason Wright, Roman R. Rafikov, Philip Chang, Thomas Beatty, Daniel Huber

Abstract: We show that a small, but \textit{measurable} shift in the eclipse mid-point time of eclipsing binary (EBs) stars of $\sim$ 0.1 seconds over a decade baseline can be used to directly measure the Galactic acceleration of stars in the Milky Way at $\sim$ kpc distances from the Sun. We consider contributions to the period drift rate from dynamical mechanisms other than the Galaxy's gravitational fiel… ▽ More We show that a small, but \textit{measurable} shift in the eclipse mid-point time of eclipsing binary (EBs) stars of $\sim$ 0.1 seconds over a decade baseline can be used to directly measure the Galactic acceleration of stars in the Milky Way at $\sim$ kpc distances from the Sun. We consider contributions to the period drift rate from dynamical mechanisms other than the Galaxy's gravitational field, and show that the Galactic acceleration can be reliably measured using a sample of $\textit{Kepler}$ EBs with orbital and stellar parameters from the literature. Given the uncertainties on the formulation of tidal decay, our approach here is necessarily approximate, and the contribution from tidal decay is an upper limit assuming the stars are not tidally synchronized. We also use simple analytic relations to search for well-timed sources in the \textit{Kepler} field, and find $\sim$ 70 additional detached EBs with low eccentricities that have estimated timing precision better than 1 second. We illustrate the method with a prototypical, precisely timed EB using an archival \textit{Kepler} light curve and a modern synthetic \textit{HST} light curve (which provides a decade baseline). This novel method establishes a realistic possibility for obtaining fundamental Galactic parameters using eclipse timing to measure Galactic accelerations, along with other emerging new methods, including pulsar timing and extreme precision radial velocity observations. This acceleration signal grows quadratically with time. Therefore, given baselines established in the near-future for distant EBs, we can expect to measure the period drift in the future with space missions like \textit{JWST} and the \textit{Roman Space Telescope}. △ Less

Submitted 15 December, 2021; originally announced December 2021.

Comments: 7 pages, 3 figures, submitted to ApJ Letters

Tests of General Relativity with GWTC-3

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

Abstract: The ever-increasing number of detections of gravitational waves (GWs) from compact binaries by the Advanced LIGO and Advanced Virgo detectors allows us to perform ever-more sensitive tests of general relativity (GR) in the dynamical and strong-field regime of gravity. We perform a suite of tests of GR using the compact binary signals observed during the second half of the third observing run of th… ▽ More The ever-increasing number of detections of gravitational waves (GWs) from compact binaries by the Advanced LIGO and Advanced Virgo detectors allows us to perform ever-more sensitive tests of general relativity (GR) in the dynamical and strong-field regime of gravity. We perform a suite of tests of GR using the compact binary signals observed during the second half of the third observing run of those detectors. We restrict our analysis to the 15 confident signals that have false alarm rates $\leq 10^{-3}\, {\rm yr}^{-1}$. In addition to signals consistent with binary black hole (BH) mergers, the new events include GW200115_042309, a signal consistent with a neutron star--BH merger. We find the residual power, after subtracting the best fit waveform from the data for each event, to be consistent with the detector noise. Additionally, we find all the post-Newtonian deformation coefficients to be consistent with the predictions from GR, with an improvement by a factor of ~2 in the -1PN parameter. We also find that the spin-induced quadrupole moments of the binary BH constituents are consistent with those of Kerr BHs in GR. We find no evidence for dispersion of GWs, non-GR modes of polarization, or post-merger echoes in the events that were analyzed. We update the bound on the mass of the graviton, at 90% credibility, to $m_g \leq 1.27 \times 10^{-23} \mathrm{eV}/c^2$. The final mass and final spin as inferred from the pre-merger and post-merger parts of the waveform are consistent with each other. The studies of the properties of the remnant BHs, including deviations of the quasi-normal mode frequencies and damping times, show consistency with the predictions of GR. In addition to considering signals individually, we also combine results from the catalog of GW signals to calculate more precise population constraints. We find no evidence in support of physics beyond GR. △ Less

Submitted 13 December, 2021; originally announced December 2021.

Report number: LIGO-P2100275

All-sky search for gravitational wave emission from scalar boson clouds around spinning black holes in LIGO O3 data

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

Abstract: This paper describes the first all-sky search for long-duration, quasi-monochromatic gravitational-wave signals emitted by ultralight scalar boson clouds around spinning black holes using data from the third observing run of Advanced LIGO. We analyze the frequency range from 20~Hz to 610~Hz, over a small frequency derivative range around zero, and use multiple frequency resolutions to be robust to… ▽ More This paper describes the first all-sky search for long-duration, quasi-monochromatic gravitational-wave signals emitted by ultralight scalar boson clouds around spinning black holes using data from the third observing run of Advanced LIGO. We analyze the frequency range from 20~Hz to 610~Hz, over a small frequency derivative range around zero, and use multiple frequency resolutions to be robust towards possible signal frequency wanderings. Outliers from this search are followed up using two different methods, one more suitable for nearly monochromatic signals, and the other more robust towards frequency fluctuations. We do not find any evidence for such signals and set upper limits on the signal strain amplitude, the most stringent being $\approx10^{-25}$ at around 130~Hz. We interpret these upper limits as both an "exclusion region" in the boson mass/black hole mass plane and the maximum detectable distance for a given boson mass, based on an assumption of the age of the black hole/boson cloud system. △ Less

Submitted 9 May, 2022; v1 submitted 30 November, 2021; originally announced November 2021.

Comments: 28 pages, 16 figures

Report number: P2100343

Journal ref: Phys. Rev. D 105, 102001, 2022

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in the Second and Third LIGO-Virgo Observing Runs

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

Abstract: We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both… ▽ More We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both once and twice the rotation frequency (dual harmonic). No evidence of GWs was found so we present 95\% credible upper limits on the strain amplitudes $h_0$ for the single harmonic search along with limits on the pulsars' mass quadrupole moments $Q_{22}$ and ellipticities $\varepsilon$. Of the pulsars studied, 23 have strain amplitudes that are lower than the limits calculated from their electromagnetically measured spin-down rates. These pulsars include the millisecond pulsars J0437\textminus4715 and J0711\textminus6830 which have spin-down ratios of 0.87 and 0.57 respectively. For nine pulsars, their spin-down limits have been surpassed for the first time. For the Crab and Vela pulsars our limits are factors of $\sim 100$ and $\sim 20$ more constraining than their spin-down limits, respectively. For the dual harmonic searches, new limits are placed on the strain amplitudes $C_{21}$ and $C_{22}$. For 23 pulsars we also present limits on the emission amplitude assuming dipole radiation as predicted by Brans-Dicke theory. △ Less

Submitted 20 July, 2022; v1 submitted 25 November, 2021; originally announced November 2021.

Comments: 37 pages

Report number: LIGO-P2100049

Constraints on the cosmic expansion history from GWTC-3

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

Abstract: We use 47 gravitational-wave sources from the Third LIGO-Virgo-KAGRA Gravitational-Wave Transient Catalog (GWTC-3) to estimate the Hubble parameter $H(z)$, including its current value, the Hubble constant $H_0$. Each gravitational-wave (GW) signal provides the luminosity distance to the source and we estimate the corresponding redshift using two methods: the redshifted masses and a galaxy catalog.… ▽ More We use 47 gravitational-wave sources from the Third LIGO-Virgo-KAGRA Gravitational-Wave Transient Catalog (GWTC-3) to estimate the Hubble parameter $H(z)$, including its current value, the Hubble constant $H_0$. Each gravitational-wave (GW) signal provides the luminosity distance to the source and we estimate the corresponding redshift using two methods: the redshifted masses and a galaxy catalog. Using the binary black hole (BBH) redshifted masses, we simultaneously infer the source mass distribution and $H(z)$. The source mass distribution displays a peak around $34\, {\rm M_\odot}$, followed by a drop-off. Assuming this mass scale does not evolve with redshift results in a $H(z)$ measurement, yielding $H_0=68^{+12}_{-7} {\rm km\,s^{-1}\,Mpc^{-1}}$ ($68\%$ credible interval) when combined with the $H_0$ measurement from GW170817 and its electromagnetic counterpart. This represents an improvement of 17% with respect to the $H_0$ estimate from GWTC-1. The second method associates each GW event with its probable host galaxy in the catalog GLADE+, statistically marginalizing over the redshifts of each event's potential hosts. Assuming a fixed BBH population, we estimate a value of $H_0=68^{+8}_{-6} {\rm km\,s^{-1}\,Mpc^{-1}}$ with the galaxy catalog method, an improvement of 42% with respect to our GWTC-1 result and 20% with respect to recent $H_0$ studies using GWTC-2 events. However, we show that this result is strongly impacted by assumptions about the BBH source mass distribution; the only event which is not strongly impacted by such assumptions (and is thus informative about $H_0$) is the well-localized event GW190814. △ Less

Submitted 19 November, 2021; v1 submitted 5 November, 2021; originally announced November 2021.

Comments: Main paper: 30 pages, 15 figure, 7 tables

Report number: LIGO-P2100185-v6

The Effect of Impact Parameter on Tidal Disruption Events

Authors: Alexandra Spaulding, Philip Chang

Abstract: Stars that pass too close to a supermassive black hole are disrupted by the black hole's tidal gravity and leads to some debris being ejected while the remainder accretes into the black hole. To better study the physics of this debris, we use the moving mesh code MANGA to follow the evolution of the star from its initial encounter to its complete destruction. By varying the impact parameter ($β$)… ▽ More Stars that pass too close to a supermassive black hole are disrupted by the black hole's tidal gravity and leads to some debris being ejected while the remainder accretes into the black hole. To better study the physics of this debris, we use the moving mesh code MANGA to follow the evolution of the star from its initial encounter to its complete destruction. By varying the impact parameter ($β$) of the star, we studied the energy distribution of the remaining material and the fallback rate of the material into the black hole as a function of time. We show that the spread of energy in the debris and peak luminosity time ($t_{\rm peak}$) are both directly related to the impact parameter. In particular, we find a $β^{1/2}$ scaling for the energy spread for $β=2-10$ and a frozen evolution for $β\gtrsim 10$. We discuss implication of this scaling for the rise time of the light curve and broadness of the luminosity peak for these lower $β$'s. These relationships provide a possible means of inferring the impact parameters for observed TDEs. △ Less

Submitted 18 November, 2020; v1 submitted 14 October, 2020; originally announced October 2020.

Comments: 8 pages, 6 figures, submitted to MNRAS

A measurement of the Galactic plane mass density from binary pulsar accelerations

Authors: Sukanya Chakrabarti, Philip Chang, Michael T. Lam, Sarah J. Vigeland, Alice C. Quillen

Abstract: We use compiled high-precision pulsar timing measurements to directly measure the Galactic acceleration of binary pulsars relative to the Solar System barycenter. Given the vertical accelerations, we use the Poisson equation to derive the Oort limit, i.e., the total volume mass density in the Galactic mid-plane. Our best-fitting model gives an Oort limit of… ▽ More We use compiled high-precision pulsar timing measurements to directly measure the Galactic acceleration of binary pulsars relative to the Solar System barycenter. Given the vertical accelerations, we use the Poisson equation to derive the Oort limit, i.e., the total volume mass density in the Galactic mid-plane. Our best-fitting model gives an Oort limit of $0.08^{0.05}_{-0.02} M_{\odot}/\rm pc^{3}$, which is close to estimates from recent Jeans analyses. Given the accounting of the baryon budget from McKee et al. (2015), we obtain a local dark matter density of $-0.004^{0.05}_{-0.02}~M_{\odot}/\rm pc^{3}$, which is slightly below other modern estimates but consistent within the current uncertainties of our method. While this first measurement of the Oort limit (and other Galactic parameters) has error bars that are currently several times larger than kinematical estimates, they should improve in the future. We also constrain the oblateness of the potential, finding it consistent with that expected from the disk and inconsistent with a potential dominated by a spherical halo, as is appropriate for our sample which is within a $\sim$ kpc of the Sun. We find that the slope of the rotation curve is not constrained by current measurements of binary pulsar accelerations. We give a fitting function for the vertical acceleration $a_{z}$: $a_{z} = -α_{1}z$; $\log_{10} (α_{1}/{\rm Gyr}^{-2})=3.69^{0.19}_{-0.12}$. By analyzing interacting simulations of the Milky Way, we find that large asymmetric variations in $da_{z}/dz$ as a function of vertical height may be a signature of sub-structure. We end by discussing the power of combining constraints from pulsar timing and high-precision radial velocity (RV) measurements towards lines-of-sight near pulsars, to test theories of gravity and constrain dark matter sub-structure. △ Less

Submitted 1 January, 2021; v1 submitted 8 October, 2020; originally announced October 2020.

Comments: accepted to ApJ Letters, changes following referee comments

Towards a direct measure of the Galactic acceleration

Authors: Sukanya Chakrabarti, Jason Wright, Philip Chang, Alice Quillen, Peter Craig, Joey Territo, Elena D'Onghia, Kathryn V. Johnston, Robert J. De Rosa, Daniel Huber, Katherine L. Rhode, Eric Nielsen

Abstract: High precision spectrographs can enable not only the discovery of exoplanets, but can also provide a fundamental measurement in Galactic dynamics. Over about ten year baselines, the expected change in the line-of-sight velocity due to the Galaxy's gravitational field for stars at $\sim$ kpc scale distances above the Galactic mid-plane is $\sim$ few - 10 cm/s, and may be detectable by the current g… ▽ More High precision spectrographs can enable not only the discovery of exoplanets, but can also provide a fundamental measurement in Galactic dynamics. Over about ten year baselines, the expected change in the line-of-sight velocity due to the Galaxy's gravitational field for stars at $\sim$ kpc scale distances above the Galactic mid-plane is $\sim$ few - 10 cm/s, and may be detectable by the current generation of high precision spectrographs. Here, we provide theoretical expectations for this measurement based on both static models of the Milky Way and isolated Milky Way simulations, as well from controlled dynamical simulations of the Milky Way interacting with dwarf galaxies. We simulate a population synthesis model to analyze the contribution of planets and binaries to the Galactic acceleration signal. We find that while low-mass, long-period planetary companions are a contaminant to the Galactic acceleration signal, their contribution is very small. Our analysis of $\sim$ ten years of data from the LCES HIRES/Keck precision radial velocity (RV) survey shows that slopes of the RV curves of standard RV stars agree with expectations of the local Galactic acceleration near the Sun within the errors, and that the error in the slope scales inversely as the square root of the number of observations. Thus, we demonstrate that a survey of stars with low intrinsic stellar jitter at kpc distances above the Galactic mid-plane for realistic sample sizes can enable a direct determination of the dark matter density. △ Less

Submitted 22 September, 2020; v1 submitted 29 July, 2020; originally announced July 2020.

Comments: accepted to ApJ Letters

The effect of diffusive nuclear burning in neutron star envelopes on cooling in accreting systems

Authors: M. J. P Wijngaarden, Wynn C. G. Ho, Philip Chang, Dany Page, Rudy Wijnands, Laura S. Ootes, Andrew Cumming, Nathalie Degenaar, Mikhail Beznogov

Abstract: Valuable information about the neutron star interior can be obtained by comparing observations of thermal radiation from a cooling neutron star crust with theoretical models. Nuclear burning of lighter elements that diffuse to deeper layers of the envelope can alter the relation between surface and interior temperatures and can change the chemical composition over time. We calculate new temperatur… ▽ More Valuable information about the neutron star interior can be obtained by comparing observations of thermal radiation from a cooling neutron star crust with theoretical models. Nuclear burning of lighter elements that diffuse to deeper layers of the envelope can alter the relation between surface and interior temperatures and can change the chemical composition over time. We calculate new temperature relations and consider two effects of diffusive nuclear burning (DNB) for H-C envelopes. First, we consider the effect of a changing envelope composition and find that hydrogen is consumed on short timescales and our temperature evolution simulations correspond to those of a hydrogen-poor envelope within ~100 days. The transition from a hydrogen-rich to a hydrogen-poor envelope is potentially observable in accreting NS systems as an additional initial decline in surface temperature at early times after the outburst. Second, we find that DNB can produce a non-negligible heat flux, such that the total luminosity can be dominated by DNB in the envelope rather than heat from the deep interior. However, without continual accretion, heating by DNB in H-C envelopes is only relevant for <1-80 days after the end of an accretion outburst, as the amount of light elements is rapidly depleted. Comparison to crust cooling data shows that DNB does not remove the need for an additional shallow heating source. We conclude that solving the time-dependent equations of the burning region in the envelope self-consistently in thermal evolution models instead of using static temperature relations would be valuable in future cooling studies. △ Less

Submitted 20 March, 2020; originally announced March 2020.

Comments: 10 pages, 10 figures; Published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 493, Issue 4, April 2020, Pages 4936-4944

The growth of the longitudinal beam-plasma instability in the presence of an inhomogeneous background

Authors: Mohamad Shalaby, Avery E. Broderick, Philip Chang, Christoph Pfrommer, Ewald Puchwein, Astrid Lamberts

Abstract: We study the longitudinal stability of beam-plasma systems in the presence of a density inhomogeneity in the background plasma. Previous works have focused on the non-relativistic regime where hydrodynamical models are used to evolve pre-existing Langmuir waves within inhomogeneous background plasmas. Here, for the first time we study the problem with kinetic equations in a fully-relativistic way.… ▽ More We study the longitudinal stability of beam-plasma systems in the presence of a density inhomogeneity in the background plasma. Previous works have focused on the non-relativistic regime where hydrodynamical models are used to evolve pre-existing Langmuir waves within inhomogeneous background plasmas. Here, for the first time we study the problem with kinetic equations in a fully-relativistic way. We do not assume the existence of Langmuir waves, and we focus on the rate and the mechanism by which waves are excited in such systems from an initial perturbation. We derive the structure of the unstable modes and compute an analytical approximation for their growth rates. Our computation is limited to dilute and cold beams, and shows an excellent agreement with particle-in-cell simulations performed using the SHARP code. We show that, due to such an inhomogeneity, the virulent beam-plasma instabilities in the intergalactic medium are not suppressed but their counterparts in the solar wind can be suppressed as evidenced by propagating type-III solar radio bursts. △ Less

Submitted 5 March, 2020; originally announced March 2020.

Comments: 23 pages, 5 figures, accepted for publication in the Journal of Plasma Physics; comments are welcome

General Relativistic Hydrodynamics on a Moving-mesh I: Static Spacetimes

Authors: Philip Chang, Zachariah Etienne

Abstract: We present the first-ever moving-mesh general relativistic hydrodynamics solver for static spacetimes as implemented in the code, MANGA. Our implementation builds on the architectures of MANGA and the numerical relativity Python package NRPy+. We review the general algorithm to solve these equations and, in particular, detail the time stepping; Riemann solution across moving faces; conversion betw… ▽ More We present the first-ever moving-mesh general relativistic hydrodynamics solver for static spacetimes as implemented in the code, MANGA. Our implementation builds on the architectures of MANGA and the numerical relativity Python package NRPy+. We review the general algorithm to solve these equations and, in particular, detail the time stepping; Riemann solution across moving faces; conversion between primitive and conservative variables; validation and correction of hydrodynamic variables; and mapping of the metric to a Voronoi moving-mesh grid. We present test results for the numerical integration of an unmagnetized Tolman-Oppenheimer-Volkoff star for 24 dynamical times. We demonstrate that at a resolution of $10^6$ mesh generating points, the star is stable and its central density drifts downward by 2% over this timescale. At a lower resolution the central density drift increases in a manner consistent with the adopted second order spatial reconstruction scheme. These results agree well with the exact solutions, and we find the error behavior to be similar to Eulerian codes with second-order spatial reconstruction. We also demonstrate that the new code recovers the fundamental mode frequency for the same TOV star but with its initial pressure depleted by 10% △ Less

Submitted 28 May, 2020; v1 submitted 21 February, 2020; originally announced February 2020.

Comments: 10 pages, 4 figures, accepted version, update description of Riemann solver

Time Dependent Radiation Hydrodynamics on a Moving Mesh

Authors: Philip Chang, Shane W. Davis, Yan-Fei Jiang

Abstract: We describe the structure and implementation of a radiation hydrodynamic solver for MANGA, the moving-mesh hydrodynamics module of the large-scale parallel code, Charm N-body GrAvity solver (ChaNGa). We solve the equations of time dependent radiative transfer using a reduced speed of light approximation following the algorithm of Jiang et al (2014). By writing the radiative transfer equations as a… ▽ More We describe the structure and implementation of a radiation hydrodynamic solver for MANGA, the moving-mesh hydrodynamics module of the large-scale parallel code, Charm N-body GrAvity solver (ChaNGa). We solve the equations of time dependent radiative transfer using a reduced speed of light approximation following the algorithm of Jiang et al (2014). By writing the radiative transfer equations as a generalized conservation equation, we solve the transport part of these equations on an unstructured Voronoi mesh. We then solve the source part of the radiative transfer equations following Jiang et al (2014) using an implicit solver, and couple this to the hydrodynamic equations. The use of an implicit solver ensure reliable convergence and preserves the conservation properties of these equations even in situations where the source terms are stiff due to the small coupling timescales between radiation and matter. We present the results of a limited number of test cases (energy conservation, momentum conservation, dynamic diffusion, linear waves, crossing beams, and multiple shadows) to show convergence with analytic results and numerical stability. We also show that it produces qualitatively the correct results in the presence of multiple sources in the optically thin case. △ Less

Submitted 25 February, 2020; v1 submitted 19 February, 2020; originally announced February 2020.

Comments: 13 pages, 8 figures, accepted to MNRAS

Simulation of a Compact Object with Outflows Moving Through a Gaseous Background

Authors: Xinyu Li, Philip Chang, Yuri Levin, Christopher D. Matzner, Philip J. Armitage

Abstract: A compact object moving relative to surrounding gas accretes material and perturbs the density of gas in its vicinity. In the classical picture of Bondi-Hoyle-Lyttleton accretion, the perturbation takes the form of an overdense wake behind the object, which exerts a dynamical friction drag. We use hydrodynamic simulations to investigate how the accretion rate and strength of dynamical friction are… ▽ More A compact object moving relative to surrounding gas accretes material and perturbs the density of gas in its vicinity. In the classical picture of Bondi-Hoyle-Lyttleton accretion, the perturbation takes the form of an overdense wake behind the object, which exerts a dynamical friction drag. We use hydrodynamic simulations to investigate how the accretion rate and strength of dynamical friction are modified by the presence of outflow from the compact object. We show that the destruction of the wake by an outflow reduces dynamical friction, and reverses its sign when the outflow is strong enough, in good quantitative agreement with analytic calculations. For a strong isotropic outflow, the outcome on scales that we have simulated is a negative dynamical friction, i.e., net acceleration. For jet-like outflows driven by reprocessed accretion, both the rate of accretion and the magnitude of dynamical friction drop for more powerful jets. The accretion rate is strongly intermittent when the jet points to the same direction as the motion of the compact object. The dynamical effects of outflows may be important for the evolution of compact objects during the common envelope phase of binary systems, and for accreting compact objects and massive stars encountering AGN discs. △ Less

Submitted 30 March, 2020; v1 submitted 14 December, 2019; originally announced December 2019.

Comments: Accepted by MNRAS

Thermal Loop Effects on Large-Scale Curvature Perturbation in the Higgs Inflation

Authors: Po-Wen Chang, Cheng-Wei Chiang, Kin-Wang Ng

Abstract: It is known that the Higgs potential in the Standard Model can drive a successful inflation as long as the Higgs field couples non-minimally to gravity. It is then inevitable to take into account the loop corrections of the Standard Model particles to the Higgs potential in the Higgs inflation. In this paper, we discuss the one-loop corrections at finite temperature to the curvature perturbation g… ▽ More It is known that the Higgs potential in the Standard Model can drive a successful inflation as long as the Higgs field couples non-minimally to gravity. It is then inevitable to take into account the loop corrections of the Standard Model particles to the Higgs potential in the Higgs inflation. In this paper, we discuss the one-loop corrections at finite temperature to the curvature perturbation generated during the Higgs inflation. We find that the thermal loop effects can suppress the power of the curvature perturbation at large scales, thus resulting in a low quadrupole of the temperature anisotropy in the cosmic microwave background. △ Less

Submitted 29 April, 2020; v1 submitted 10 July, 2019; originally announced July 2019.

Comments: 24 pages, 3 figures, 2 tables, expanded discussion on RG evolution and thermal equilibrium condition, references added, matches published version

Journal ref: JHEP 2004 (2020) 163

Antlia2's role in driving the ripples in the outer gas disk of the Galaxy

Authors: Sukanya Chakrabarti, Philip Chang, Adrian Price-Whelan, Justin Read, Leo Blitz, Lars Hernquist

Abstract: We employ the earlier published proper motions of the newly discovered Antlia 2 dwarf galaxy derived from Gaia data to calculate its orbital distribution in the cosmologically recent past. Using these observationally motivated orbits, we calculate the effect of the Antlia 2 dwarf galaxy on the outer HI disk of the Milky Way, using both test particle and Smoothed Particle Hydrodynamics simulations.… ▽ More We employ the earlier published proper motions of the newly discovered Antlia 2 dwarf galaxy derived from Gaia data to calculate its orbital distribution in the cosmologically recent past. Using these observationally motivated orbits, we calculate the effect of the Antlia 2 dwarf galaxy on the outer HI disk of the Milky Way, using both test particle and Smoothed Particle Hydrodynamics simulations. We find that orbits with low pericenters, $\sim$ 10 kpc, produce disturbances that match the observed outer HI disk perturbations. We have independently recalculated the proper motion of the Antlia 2 dwarf from Gaia data and found a proper motion of $(μ_αcosδ, μ_δ) = (-0.068,0.032) \pm (0.023,-0.031)~\rm mas/yr$, which agrees with results from Torrealba et al. (2019) within the errors, but gives lower mean pericenters, e.g., $\sim$ 15 kpc for our fiducial model of the Milky Way. We also show that the Sagittarius dwarf galaxy interaction does not match the observed perturbations in the outer gas disk. Thus, Antlia 2 may be the driver of the observed large perturbations in the outer gas disk of the Galaxy. The current location of the Antlia 2 dwarf galaxy closely matches that predicted by an earlier dynamical analysis (Chakrabarti \& Blitz 2009) of the dwarf that drove ripples in the outer Galaxy, and, in particular, its orbit is nearly coplanar to the Galactic disk. If the Antlia 2 dwarf galaxy is responsible for the perturbations in the outer Galactic disk, it would have a specific range of proper motions that we predict here; this can be tested soon with Gaia DR-3 and Gaia DR-4 data. △ Less

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

Comments: accepted to ApJ

Common Envelope Evolution on a Moving Mesh

Authors: Logan J. Prust, Philip Chang

Abstract: We outline the methodology of simulating common envelope evolution (CEE) with the moving-mesh code MANGA. We extend MANGA to include multiple time-steps. This provides substantial speedups for problems with large dynamic range. We describe the implementation of realistic equations of state relevant in stellar structure and the generation of suitable initial conditions. We then carry out two exampl… ▽ More We outline the methodology of simulating common envelope evolution (CEE) with the moving-mesh code MANGA. We extend MANGA to include multiple time-steps. This provides substantial speedups for problems with large dynamic range. We describe the implementation of realistic equations of state relevant in stellar structure and the generation of suitable initial conditions. We then carry out two example simulations of a 2 $\rm{M}_{\odot}$ red giant with a 0.36 $\rm{M}_{\odot}$ core and a 1 $\rm{M}_{\odot}$ companion undergoing CEE for 240 days. In one simulation the red giant is set into corotation with the orbital motion and in the other it is non-rotating. We find that the separation between the companion and red giant core shrinks from 52 $\rm{R}_{\odot}$ to 3.6 $\rm{R}_{\odot}$ and 3.2 $\rm{R}_{\odot}$ respectively, ending with an eccentricity of 0.1. We also find that 66 and 63 per cent of the envelope mass is ejected. This is higher than in many previous works. Several reasons for this are discussed. These include our inclusion of recombination energy. Our simulations show that putting giants in corotation increases the fraction of mass ejected from the system and results in a larger final orbital separation. We conclude that the entire envelope of the red giant might be ejected during the plunge phase of CEE in this region of parameter space. △ Less

Submitted 19 April, 2019; originally announced April 2019.

Cooling of the Cassiopeia A neutron star and the effect of diffusive nuclear burning

Authors: Wynn C. G. Ho, M. J. P. Wijngaarden, Philip Chang, Craig O. Heinke, Dany Page, Mikhail Beznogov, Daniel J. Patnaude

Abstract: The study of how neutron stars cool over time can provide invaluable insights into fundamental physics such as the nuclear equation of state and superconductivity and superfluidity. A critical relation in neutron star cooling is the one between observed surface temperature and interior temperature. This relation is determined by the composition of the neutron star envelope and can be influenced by… ▽ More The study of how neutron stars cool over time can provide invaluable insights into fundamental physics such as the nuclear equation of state and superconductivity and superfluidity. A critical relation in neutron star cooling is the one between observed surface temperature and interior temperature. This relation is determined by the composition of the neutron star envelope and can be influenced by the process of diffusive nuclear burning (DNB). We calculate models of envelopes that include DNB and find that DNB can lead to a rapidly changing envelope composition which can be relevant for understanding the long-term cooling behavior of neutron stars. We also report on analysis of the latest temperature measurements of the young neutron star in the Cassiopeia A supernova remnant. The 13 Chandra observations over 18 years show that the neutron star's temperature is decreasing at a rate of 2-3 percent per decade, and this rapid cooling can be explained by the presence of a proton superconductor and neutron superfluid in the core of the star. △ Less

Submitted 16 April, 2019; originally announced April 2019.

Comments: 7 pages, 7 figures; to appear in the AIP Conference Proceedings of the Xiamen-CUSTIPEN Workshop on the EOS of Dense Neutron-Rich Matter in the Era of Gravitational Wave Astronomy (January 3-7, 2019, Xiamen, China)

Journal ref: AIP Conf. Proc. 2127, 020007 (2019)

Cyberinfrastructure Requirements to Enhance Multi-messenger Astrophysics

Authors: Philip Chang, Gabrielle Allen, Warren Anderson, Federica B. Bianco, Joshua S. Bloom, Patrick R. Brady, Adam Brazier, S. Bradley Cenko, Sean M. Couch, Tyce DeYoung, Ewa Deelman, Zachariah B Etienne, Ryan J. Foley, Derek B Fox, V. Zach Golkhou, Darren R Grant, Chad Hanna, Kelly Holley-Bockelmann, D. Andrew Howell, E. A. Huerta, Margaret W. G. Johnson, Mario Juric, David L. Kaplan, Daniel S. Katz, Azadeh Keivani , et al. (17 additional authors not shown)

Abstract: The identification of the electromagnetic counterpart of the gravitational wave event, GW170817, and discovery of neutrinos and gamma-rays from TXS 0506+056 heralded the new era of multi-messenger astrophysics. As the number of multi-messenger events rapidly grow over the next decade, the cyberinfrastructure requirements to handle the increase in data rates, data volume, need for event follow up,… ▽ More The identification of the electromagnetic counterpart of the gravitational wave event, GW170817, and discovery of neutrinos and gamma-rays from TXS 0506+056 heralded the new era of multi-messenger astrophysics. As the number of multi-messenger events rapidly grow over the next decade, the cyberinfrastructure requirements to handle the increase in data rates, data volume, need for event follow up, and analysis across the different messengers will also explosively grow. The cyberinfrastructure requirements to enhance multi-messenger astrophysics will both be a major challenge and opportunity for astronomers, physicists, computer scientists and cyberinfrastructure specialists. Here we outline some of these requirements and argue for a distributed cyberinfrastructure institute for multi-messenger astrophysics to meet these challenges. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 7 pages, astro2020 white paper

Diffusive nuclear burning in cooling simulations and application to new temperature data of the Cassiopeia A neutron star

Authors: M. J. P. Wijngaarden, Wynn C. G. Ho, Philip Chang, Craig O. Heinke, Dany Page, Mikhail Beznogov, Daniel J. Patnaude

Abstract: A critical relation in the study of neutron star cooling is the one between surface temperature and interior temperature. This relation is determined by the composition of the neutron star envelope and can be affected by the process of diffusive nuclear burning (DNB), which occurs when elements diffuse to depths where the density and temperature are sufficiently high to ignite nuclear burning. We… ▽ More A critical relation in the study of neutron star cooling is the one between surface temperature and interior temperature. This relation is determined by the composition of the neutron star envelope and can be affected by the process of diffusive nuclear burning (DNB), which occurs when elements diffuse to depths where the density and temperature are sufficiently high to ignite nuclear burning. We calculate models of H-He and He-C envelopes that include DNB and obtain analytic temperature relations that can be used in neutron star cooling simulations. We find that DNB can lead to a rapidly changing envelope composition and prevents the build-up of thermally stable hydrogen columns y$_H$ > 10$^{7}$ g cm$^{-2}$, while DNB can make helium envelopes more transparent to heat flux for surface temperatures $T_s$ > 2 $\times 10^6$ K. We perform neutron star cooling simulations in which we evolve temperature and envelope composition, with the latter due to DNB and accretion from the interstellar medium. We find that a time-dependent envelope composition can be relevant for understanding the long-term cooling behaviour of isolated neutron stars. We also report on the latest Chandra observations of the young neutron star in the Cassiopeia A supernova remnant; the resulting 13 temperature measurements over more than 18 years yield a ten-year cooling rate of $\approx$ 2%. Finally, we fit the observed cooling trend of the Cassiopeia A neutron star with a model that includes DNB in the envelope. △ Less

Submitted 4 January, 2019; originally announced January 2019.

Comments: 16 pages, 15 figures, 4 tables; accepted for publication in MNRAS

Journal ref: Mon. Not. R. Astron. Soc. 484, 974-988 (2019)

Relating the HI Gas Structure of Spiral Disks to Passing Satellites

Authors: Andrew Lipnicky, Sukanya Chakrabarti, Philip Chang

Abstract: We extend the work of Chang & Chakrabarti (2011) to find simple scaling relations between the density response of the gas disk of a spiral galaxy and the pericenter distance and mass ratio of a perturbing satellite. From the analysis of results from a test particle code, we obtained a simple scaling relation for the density response due to a single satellite interacting with a galactic disk, over… ▽ More We extend the work of Chang & Chakrabarti (2011) to find simple scaling relations between the density response of the gas disk of a spiral galaxy and the pericenter distance and mass ratio of a perturbing satellite. From the analysis of results from a test particle code, we obtained a simple scaling relation for the density response due to a single satellite interacting with a galactic disk, over a wide range of satellite masses and pericenter distances. We have also explored the effects of multiple satellites on the galactic disk, focusing on cases that are commonly found in cosmological simulations. Here, we use orbits for the satellites that are drawn from cosmological simulations. For these cases, we compare our approximate scaling relations to the density response generated by satellites, and find that for two satellite interactions, our scaling relations approximately recover the response of the galactic disk. We have also examined the observed HI data in the outskirts of several spiral galaxies from the THINGS sample and compared the observed perturbations to that of cosmological simulations and our own scaling relations. While small perturbations can be excited by satellites drawn from cosmological simulations, we find that large perturbations (such as those that are seen in some THINGS galaxies like M51) are not recovered by satellites drawn from cosmological simulations that are similar to Milky Way galaxies. △ Less

Submitted 28 August, 2018; originally announced August 2018.

Comments: 12 pages, 8 figures, accepted for publication in MNRAS

Missing Gamma-ray Halos and the Need for New Physics in the Gamma-ray Sky

Authors: Avery E. Broderick, Paul Tiede, Philip Chang, Astrid Lamberts, Christoph Pfrommer, Ewald Puchwein, Mohamad Shalaby, Maria Werhahn

Abstract: An intergalactic magnetic field stronger than $3\times10^{-13}$~G would explain the lack of a bright, extended degree-scale, GeV-energy inverse Compton component in the gamma-ray spectra of TeV-blazars. A robustly predicted consequence of the presence of such a field is the existence of degree-scale GeV-energy gamma-ray halos -- gamma-ray bow ties -- about TeV-bright active galactic nuclei, corres… ▽ More An intergalactic magnetic field stronger than $3\times10^{-13}$~G would explain the lack of a bright, extended degree-scale, GeV-energy inverse Compton component in the gamma-ray spectra of TeV-blazars. A robustly predicted consequence of the presence of such a field is the existence of degree-scale GeV-energy gamma-ray halos -- gamma-ray bow ties -- about TeV-bright active galactic nuclei, corresponding to more than half of all radio galaxies. However, the emitting regions of these halos are confined to and aligned with the direction of the relativistic jets associated with gamma-ray sources. Based on the orientation of radio jets, we align and stack corresponding degree-scale gamma-ray images of isolated Fanaroff-Riley class I and II objects and exclude the existence of these halos at overwhelming confidence, limiting the intergalactic field strength to $<10^{-15}$~G for large-scale fields and progressively larger in the diffusive regime when the correlation length of the field becomes small in comparison to 1 Mpc. When combined with prior limits on the strength of the intergalactic magnetic field, this excludes a purely magnetic explanation for the absence of halos. Thus, it requires the existence of novel physical processes that preempt the creation of halos, e.g., the presence of beam-plasma instabilities in the intergalactic medium or a drastic cutoff of the very high energy spectrum of these sources. △ Less

Submitted 26 October, 2018; v1 submitted 8 August, 2018; originally announced August 2018.

Comments: Accepted to ApJ, 20 pages, 15 figures

X-Ray Luminosity and Size Relationship of Supernova Remnants in the LMC

Authors: Po-Sheng Ou, You-Hua Chu, Pierre Maggi, Chuan-Jui Li, Un Pang Chang, Robert A. Gruendl

Abstract: The Large Magellanic Cloud (LMC) has $\sim$60 confirmed supernova remnants (SNRs). Because of the known distance, 50 kpc, the SNRs' angular sizes can be converted to linear sizes, and their X-ray observations can be used to assess X-ray luminosities ($L_X$). We have critically examined the LMC SNRs' sizes reported in the literature to determine the most plausible sizes. These sizes and the $L_X$ d… ▽ More The Large Magellanic Cloud (LMC) has $\sim$60 confirmed supernova remnants (SNRs). Because of the known distance, 50 kpc, the SNRs' angular sizes can be converted to linear sizes, and their X-ray observations can be used to assess X-ray luminosities ($L_X$). We have critically examined the LMC SNRs' sizes reported in the literature to determine the most plausible sizes. These sizes and the $L_X$ determined from \emph{XMM-Newton} observations are used to investigate their relationship in order to explore the environmental and evolutionary effects on the X-ray properties of SNRs. We find: (1) Small LMC SNRs, a few to 10 pc in size, are all of Type Ia with $L_X>10^{36}$ ergs s$^{-1}$. The scarcity of small core-collapse (CC) SNRs is a result of CCSNe exploding in the low-density interiors of interstellar bubbles blown by their massive progenitors during their main sequence phase. (2) Medium-sized (10-30 pc) CC SNRs show bifurcation in $L_X$, with the X-ray-bright SNRs either in an environment associated with molecular clouds or containing pulsars and pulsar wind nebulae and the X-ray-faint SNRs being located in low-density interstellar environments. (3) Large (size$>$30 pc) SNRs show a trend of $L_X$ fading with size, although the scatter is large. The observed relationship between $L_X$ and sizes can help constrain models of SNR evolution. △ Less

Submitted 4 August, 2018; originally announced August 2018.

Comments: Accepted for publication in the ApJ

Journal ref: The Astrophysical Journal, Volume 863, Number 2, Page 137, 2018

Growth of beam-plasma instabilities in the presence of background inhomogeneity

Authors: Mohamad Shalaby, Avery E. Broderick, Philip Chang, Christoph Pfrommer, Astrid Lamberts, Ewald Puchwein

Abstract: We explore how inhomogeneity in the background plasma number density alters the growth of electrostatic unstable wavemodes of beam plasma systems. This is particularly interesting for blazar-driven beam-plasma instabilities, which may be suppressed by inhomogeneities in the intergalactic medium as was recently claimed in the literature. Using high resolution Particle-In-Cell simulations with the S… ▽ More We explore how inhomogeneity in the background plasma number density alters the growth of electrostatic unstable wavemodes of beam plasma systems. This is particularly interesting for blazar-driven beam-plasma instabilities, which may be suppressed by inhomogeneities in the intergalactic medium as was recently claimed in the literature. Using high resolution Particle-In-Cell simulations with the SHARP code, we show that the growth of the instability is local, i.e., regions with almost homogeneous background density will support the growth of the Langmuir waves, even when they are separated by strongly inhomogeneous regions, resulting in an overall slower growth of the instability. We also show that if the background density is continuously varying, the growth rate of the instability is lower; though in all cases, the system remains within the linear regime longer and the instability is not extinguished. In all cases, the beam loses approximately the same fraction of its initial kinetic energy in comparison to the uniform case at non-linear saturation. Thus, inhomogeneities in the intergalactic medium are unlikely to suppress the growth of blazar-driven beam-plasma instabilities. △ Less

Submitted 13 April, 2018; originally announced April 2018.

Comments: 10 pages, 6 figures, Accepted by ApJ, comments welcome

The Accelerating Pace of Star Formation

Authors: Spencer Caldwell, Philip Chang

Abstract: We study the temporal and spatial distribution of star formation rates in four well-studied star-forming regions in local molecular clouds(MCs): Taurus, Perseus, $ρ$ Ophiuchi, and Orion A. Using published mass and age estimates for young stellar objects in each system, we show that the rate of star formation over the last 10 Myrs has been accelerating and is (roughly) consistent with a $t^2$ power… ▽ More We study the temporal and spatial distribution of star formation rates in four well-studied star-forming regions in local molecular clouds(MCs): Taurus, Perseus, $ρ$ Ophiuchi, and Orion A. Using published mass and age estimates for young stellar objects in each system, we show that the rate of star formation over the last 10 Myrs has been accelerating and is (roughly) consistent with a $t^2$ power law. This is in line with previous studies of the star formation history of molecular clouds and with recent theoretical studies. We further study the clustering of star formation in the Orion Nebula Cluster(ONC). We examine the distribution of young stellar objects as a function of their age by computing an effective half-light radius for these young stars subdivided into age bins. We show that the distribution of young stellar objects is broadly consistent with the star formation being entirely localized within the central region. We also find a slow radial expansion of the newly formed stars at a velocity of $v=0.17\,{\rm km\,s}^{-1}$, which is roughly the sound speed of the cold molecular gas. This strongly suggests the dense structures that form stars persist much longer than the local dynamical time. We argue that this structure is quasi-static in nature and is likely the result of the density profile approaching an attractor solution as suggested by recent analytic and numerical analysis. △ Less

Submitted 20 November, 2017; originally announced November 2017.

Comments: 7 pages, 4 figures, submitted to MNRAS

The Effects of Protostellar Jet Feedback on Turbulent Collapse

Authors: Daniel Murray, Shivam Goyal, Philip Chang

Abstract: We present results of hydrodynamic simulations of massive star forming regions with and without protostellar jets. We show that jets change the normalization of the stellar mass accretion rate, but do not strongly affect the dynamics of star formation. In particular, $M_*(t) \propto f^2 (t-t_*)^2$ where $f = 1 - f_{\rm jet}$ is the fraction of mass accreted onto the protostar, $f_{\rm jet}$ is the… ▽ More We present results of hydrodynamic simulations of massive star forming regions with and without protostellar jets. We show that jets change the normalization of the stellar mass accretion rate, but do not strongly affect the dynamics of star formation. In particular, $M_*(t) \propto f^2 (t-t_*)^2$ where $f = 1 - f_{\rm jet}$ is the fraction of mass accreted onto the protostar, $f_{\rm jet}$ is the fraction ejected by the jet, and $(t-t_*)^2$ is the time elapsed since the formation of the first star. The star formation efficiency is nonlinear in time. We find that jets have only a small effect (of order 25\%) on the accretion rate onto the protostellar disk (the "raw" accretion rate). We show that the small scale structure -- the radial density, velocity, and mass accretion profiles are very similar in the jet and no-jet cases. Finally, we show that the inclusion of jets does drive turbulence but only on small (parsec) scales. △ Less

Submitted 25 October, 2017; originally announced October 2017.

Comments: 15 pages, 14 figures

The Contribution Of Outer HI Disks To The Merging Binary Black Hole Population

Authors: Sukanya Chakrabarti, Philip Chang, Richard O'Shaughnessy, Alyson Brooks, Sijing Shen, Jillian Bellovary, Wojciech Gladysz, Chris Belczynski

Abstract: We investigate the contribution of outer HI disks to the observable population of merging black hole binaries. Like dwarf galaxies, the outer HI disks of spirals have low star formation rates and lower metallicities than the inner disks of spirals. Since low-metallicity star formation can produce more detectable compact binaries than typical star formation, the environments in the outskirts of spi… ▽ More We investigate the contribution of outer HI disks to the observable population of merging black hole binaries. Like dwarf galaxies, the outer HI disks of spirals have low star formation rates and lower metallicities than the inner disks of spirals. Since low-metallicity star formation can produce more detectable compact binaries than typical star formation, the environments in the outskirts of spiral galaxies may be conducive to producing a rich population of massive binary black holes. We consider here both detailed controlled simulations of spirals and cosmological simulations, as well as the current range of observed values for metallicity and star formation in outer disks. We find that outer HI disks contribute at least as much as dwarf galaxies do to the observed LIGO/Virgo detection rates. Identifying the host galaxies of merging massive black holes should provide constraints on cosmological parameters and insights into the formation channels of binary mergers. △ Less

Submitted 25 October, 2017; originally announced October 2017.

Comments: accepted to ApJL, 5 pages, 2 figures

GW170817: A Neutron Star Merger in a Mass-Transferring Triple System

Authors: Philip Chang, Norman Murray

Abstract: The light curve of GW170817 is surprisingly blue and bright. Assuming that the event is a binary neutron star merger, we argue that blueness and brightness of the light curve is the result of ejecta that contains an substantial amount of thermal energy. To achieve this, the ejecta must be reheated at a substantial distance (1 to 2000 solar radii) from the merger to avoid losing the energy to adiab… ▽ More The light curve of GW170817 is surprisingly blue and bright. Assuming that the event is a binary neutron star merger, we argue that blueness and brightness of the light curve is the result of ejecta that contains an substantial amount of thermal energy. To achieve this, the ejecta must be reheated at a substantial distance (1 to 2000 solar radii) from the merger to avoid losing the energy to adiabatic cooling. We show that this reheating can occur if the merger occurs in a hierarchical triple system where the outer star has evolved and filled its Roche lobe. The outer star feeds mass to the inner binary, forming a circumbinary disc, driving the inner binary to merge. Because the outer star fills its Roche lobe, a substantial fraction of the dynamical ejecta collides with the evolved star, reheating the ejecta in the process. We suggest that the process of mass transfer in heirarchical triples tends to form coplanar triple systems such as PSR J0337+1715, and may provide electromagnetic counterparts to binary black hole mergers. △ Less

Submitted 16 October, 2017; originally announced October 2017.

Comments: submitted to MNRAS Letters, 5 pages, 1 figure

A Moving Mesh Hydrodynamic Solver for ChaNGa

Authors: Philip Chang, James Wadsley, Thomas R. Quinn

Abstract: We describe the structure and implementation of a moving-mesh hydrodynamics solver in the large-scale parallel code, Charm N-body GrAvity solver (ChaNGa). While largely based on the algorithm described by Springel (2010) that is implemented in AREPO, our algorithm differs a few aspects. We describe our use of the Voronoi tessellation library, VORO++, to compute the Voronoi tessellation directly. W… ▽ More We describe the structure and implementation of a moving-mesh hydrodynamics solver in the large-scale parallel code, Charm N-body GrAvity solver (ChaNGa). While largely based on the algorithm described by Springel (2010) that is implemented in AREPO, our algorithm differs a few aspects. We describe our use of the Voronoi tessellation library, VORO++, to compute the Voronoi tessellation directly. We also incorporate some recent advances in gradient estimation and reconstruction that gives better accuracy in hydrodynamic solutions at minimal computational cost. We validate this module with a small battery of test problems against the smooth particle hydrodynamics solver included in ChaNGa. Finally, we study one example of a scientific problem involving the mergers of two main sequence stars and highlight the small quantitative differences between smooth particle and moving-mesh hydrodynamics. We close with a discussion of anticipated future improvements and advancements. △ Less

Submitted 17 July, 2017; originally announced July 2017.

Comments: 15 pages, 7 figures, accepted to MNRAS

Importance of resolving the spectral support of beam-plasma instabilities in simulations

Authors: Mohamad Shalaby, Avery E. Broderick, Philip Chang, Christoph Pfrommer, Astrid Lamberts, Ewald Puchwein

Abstract: Many astrophysical plasmas are prone to beam-plasma instabilities. For relativistic and dilute beams, the {\it spectral} support of the beam-plasma instabilities is narrow, i.e., the linearly unstable modes that grow with rates comparable to the maximum growth rate occupy a narrow range of wave numbers. This places stringent requirements on the box-sizes when simulating the evolution of the instab… ▽ More Many astrophysical plasmas are prone to beam-plasma instabilities. For relativistic and dilute beams, the {\it spectral} support of the beam-plasma instabilities is narrow, i.e., the linearly unstable modes that grow with rates comparable to the maximum growth rate occupy a narrow range of wave numbers. This places stringent requirements on the box-sizes when simulating the evolution of the instabilities. We identify the implied lower limits on the box size imposed by the longitudinal beam plasma instability, i.e., typically the most stringent condition required to correctly capture the linear evolution of the instabilities in multidimensional simulations. We find that sizes many orders of magnitude larger than the resonant wavelength are typically required. Using one-dimensional particle-in-cell simulations, we show that the failure to sufficiently resolve the spectral support of the longitudinal instability yields slower growth and lower levels of saturation, potentially leading to erroneous physical conclusion. △ Less

Submitted 4 October, 2017; v1 submitted 31 March, 2017; originally announced April 2017.

Comments: 7 pages, 9 figures, accepted by ApJ

SHARP: A Spatially Higher-order, Relativistic Particle-in-Cell Code

Authors: Mohamad Shalaby, Avery E. Broderick, Philip Chang, Christoph Pfrommer, Astrid Lamberts, Ewald Puchwein

Abstract: Numerical heating in particle-in-cell (PIC) codes currently precludes the accurate simulation of cold, relativistic plasma over long periods, severely limiting their applications in astrophysical environments. We present a spatially higher-order accurate relativistic PIC algorithm in one spatial dimension, which conserves charge and momentum exactly. We utilize the smoothness implied by the usage… ▽ More Numerical heating in particle-in-cell (PIC) codes currently precludes the accurate simulation of cold, relativistic plasma over long periods, severely limiting their applications in astrophysical environments. We present a spatially higher-order accurate relativistic PIC algorithm in one spatial dimension, which conserves charge and momentum exactly. We utilize the smoothness implied by the usage of higher-order interpolation functions to achieve a spatially higher-order accurate algorithm (up to fifth order). We validate our algorithm against several test problems -- thermal stability of stationary plasma, stability of linear plasma waves, and two-stream instability in the relativistic and non-relativistic regimes. Comparing our simulations to exact solutions of the dispersion relations, we demonstrate that SHARP can quantitatively reproduce important kinetic features of the linear regime. Our simulations have a superior ability to control energy non-conservation and avoid numerical heating in comparison to common second-order schemes. We provide a natural definition for convergence of a general PIC algorithm: the complement of physical modes captured by the simulation, i.e., those that lie above the Poisson noise, must grow commensurately with the resolution. This implies that it is necessary to simultaneously increase the number of particles per cell and decrease the cell size. We demonstrate that traditional ways for testing for convergence fail, leading to plateauing of the energy error. This new PIC code enables us to faithfully study the long-term evolution of plasma problems that require absolute control of the energy and momentum conservation. △ Less

Submitted 11 June, 2017; v1 submitted 15 February, 2017; originally announced February 2017.

Comments: 26 pages, 19 figures, discussion about performance is added, published in ApJ

Constraints on the Intergalactic Magnetic Field from Bow Ties in the Gamma-ray Sky

Authors: Paul Tiede, Avery E. Broderick, Mohamad Shalaby, Christoph Pfrommer, Ewald Puchwein, Philip Chang, Astrid Lamberts

Abstract: Pair creation on the cosmic infrared background and subsequent inverse-Compton scattering on the CMB potentially reprocesses the TeV emission of blazars into faint GeV halos with structures sensitive to intergalactic magnetic fields (IGMF). Previous work has shown that these halos are then highly-anisotropic and extended. If the coherence length of the IGMF is greater than the inverse-Compton cool… ▽ More Pair creation on the cosmic infrared background and subsequent inverse-Compton scattering on the CMB potentially reprocesses the TeV emission of blazars into faint GeV halos with structures sensitive to intergalactic magnetic fields (IGMF). Previous work has shown that these halos are then highly-anisotropic and extended. If the coherence length of the IGMF is greater than the inverse-Compton cooling length of the pairs, then the orientation of the gamma-ray halo will be correlated with the direction of the magnetic field which is unknown and expected to change for each source. In order to constructively add each source we then use angular power spectra which are insensitive to the jet orientation. By looking at known GeV blazars detected by Fermi, we exclude the existence of an IGMF with coherence lengths $>100$ Mpc at greater than $3.9σ$ with current-day strengths in the range $10^{-16}$ to $10^{-15}$ G, and at 2$σ$ from $10^{-17}$ to $10^{-14}$ G. This provides a direct measurement of the non-existence of gamma-ray halos, providing an important check on previous results. △ Less

Submitted 13 February, 2020; v1 submitted 8 February, 2017; originally announced February 2017.

Comments: 12 pages, 14 figures, 1 appendix. Accepted to ApJ