Search

The first two months of the third Advanced LIGO and Virgo observing run (2019 April-May) showed that distant gravitational-wave (GW) events can now be readily detected. Three candidate mergers containing neutron stars (NS) were reported in a span of 15 days, all likely located more than 100 Mpc away. However, distant events such as the three new NS mergers are likely to be coarsely localized, which highlights the importance of facilities and scheduling systems that enable deep observations over hundreds to thousands of square degrees to detect the electromagnetic counterparts. On 2019 May 10 02:59:39.292 UT the GW candidate S190510g was discovered and initially classified as a binary neutron star (BNS) merger with 98% probability. The GW event was localized within an area of 3462 deg2, later refined to 1166 deg2 (90%) at a distance of 227 92 Mpc. We triggered Target-of-Opportunity observations with the Dark Energy Camera (DECam), a wide-field optical imager mounted at the prime focus of the 4 m Blanco Telescope at Cerro Tololo Inter-American Observatory in Chile. This Letter describes our DECam observations and our real-time analysis results, focusing in particular on the design and implementation of the observing strategy. Within 24 hr of the merger time, we observed 65% of the total enclosed probability of the final skymap with an observing efficiency of 94%. We identified and publicly announced 13 candidate counterparts. S190510g was reclassified 1.7 days after the merger, after our observations were completed, with a "BNS merger" probability reduced from 98% to 42% in favor of a "terrestrial classification.

The first two months of the third Advanced LIGO and Virgo observing run (2019 April-May) showed that distant gravitational-wave (GW) events can now be readily detected. Three candidate mergers containing neutron stars (NS) were reported in a span of 15 days, all likely located more than 100 Mpc away. However, distant events such as the three new NS mergers are likely to be coarsely localized, which highlights the importance of facilities and scheduling systems that enable deep observations over hundreds to thousands of square degrees to detect the electromagnetic counterparts. On 2019 May 10 02:59:39.292 UT the GW candidate S190510g was discovered and initially classified as a binary neutron star (BNS) merger with 98% probability. The GW event was localized within an area of 3462 deg2, later refined to 1166 deg2 (90%) at a distance of 227 92 Mpc. We triggered Target-of-Opportunity observations with the Dark Energy Camera (DECam), a wide-field optical imager mounted at the prime focus of the 4 m Blanco Telescope at Cerro Tololo Inter-American Observatory in Chile. This Letter describes our DECam observations and our real-time analysis results, focusing in particular on the design and implementation of the observing strategy. Within 24 hr of the merger time, we observed 65% of the total enclosed probability of the final skymap with an observing efficiency of 94%. We identified and publicly announced 13 candidate counterparts. S190510g was reclassified 1.7 days after the merger, after our observations were completed, with a "BNS merger" probability reduced from 98% to 42% in favor of a "terrestrial classification.

The Galactic center is dominated by the gravity of a super-massive black hole (SMBH), Sagittarius A∗, and is suspected to contain a sizable population of binary stars. Such binaries form hierarchical triples with the SMBH, undergoing Eccentric Kozai-Lidov (EKL) evolution, which can lead to high-eccentricity excitations for the binary companions' mutual orbit. This effect can lead to stellar collisions or Roche-lobe crossings, as well as orbital shrinking due to tidal dissipation. In this work we investigate the dynamical and stellar evolution of such binary systems, especially with regards to the binaries' post-main-sequence evolution. We find that the majority of binaries (∼75%) is eventually separated into single stars, while the remaining binaries (∼25%) undergo phases of common-envelope evolution and/or stellar mergers. These objects can produce a number of different exotic outcomes, including rejuvenated stars, G2-like infrared-excess objects, stripped giant stars, Type Ia supernovae (SNe), cataclysmic variables, symbiotic binaries, or compact object binaries. We estimate that, within a sphere of 250 Mpc radius, about 7.5-15 SNe Ia per year should occur in galactic nuclei due to this mechanism, potentially detectable by the Zwicky Transient Facility and ASAS-SN. Likewise, we estimate that, within a sphere of 1 Gpc3 volume, about 10-20 compact object binaries form per year that could become gravitational wave sources. Based on results of EKL-driven compact object binary mergers in galactic nuclei by Hoang et al., this compact object binary formation rate translates to about 15-30 events per year that are detectable by Advanced LIGO.

The Galactic Center (GC) is dominated by the gravity of a super-massive blackhole (SMBH), Sagittarius A$^*$, and is suspected to contain a sizablepopulation of binary stars. Such binaries form hierarchical triples with theSMBH, undergoing Eccentric Kozai-Lidov (EKL) evolution, which can lead to higheccentricity excitations for the binary companions' mutual orbit. This effectcan lead to stellar collisions or Roche-lobe crossings, as well as orbitalshrinking due to tidal dissipation. In this work we investigate the dynamicaland stellar evolution of such binary systems, especially with regards to thebinaries' post-main-sequence evolution. We find that the majority of binaries(~75%) is eventually separated into single stars, while the remaining binaries(~25%) undergo phases of common-envelope evolution and/or stellar mergers.These objects can produce a number of different exotic outcomes, includingrejuvenated stars, G2-like infrared-excess objects, stripped giant stars, TypeIa supernovae (SNe), cataclysmic variables (CVs), symbiotic binaries (SBs), orcompact object binaries. We estimate that, within a sphere of 250 Mpc radius,about 7.5 to 15 Type Ia SNe per year should occur in galactic nuclei due tothis mechanism, potentially detectable by ZTF and ASAS-SN. Likewise we estimatethat, within a sphere of 1 Gpc$^3$ volume, about 10 to 20 compact objectbinaries form per year that could become gravitational wave sources. Based onresults of EKL-driven compact object binary mergers in galactic nuclei by Hoangat al. (2018), this compact object binary formation rate translates to about 15to 30 events per year detectable by Advanced LIGO.

We report the discovery of a giant planet in event UKIRT-2017-BLG-001, detected by the United Kingdom Infrared Telescope (UKIRT) microlensing survey. The mass ratio between the planet and its host is =-q 1.50+0.170-0.14 0.17 3, about 1.5 times the Jupiter/Sun mass ratio. The event lies 0 .35 from the Galactic center and suffers from high extinction of AK=1.68. Therefore, it could be detected only by a near-infrared (NIR) survey. The field also suffers from large spatial differential extinction, which makes it difficult to estimate the source properties required to derive the angular Einstein radius. Nevertheless, we find evidence suggesting that the source is located in the far disk. If correct, this would be the first source star of a microlensing event to be identified as belonging to the far disk. We estimate the lens mass and distance using a Bayesian analysis to find that the planet's mass is 1.28+ M0.44 J 0.37 , and it orbits a-+ 0.81 +0.21-0.27 Ṁ 0.21 star at an instantaneous projected separation of 4.18+0.96-0.88au. The system is at a distance of-6.3+16-2.1kpc, and so likely resides in the Galactic bulge. In addition, we find a non-standard extinction curve in this field, in agreement with previous results toward high-extinction fields near the Galactic center.

We report the discovery of a giant planet in event UKIRT-2017-BLG-001, detected by the United Kingdom Infrared Telescope (UKIRT) microlensing survey. The mass ratio between the planet and its host is =-q 1.50+0.170-0.14 0.17 3, about 1.5 times the Jupiter/Sun mass ratio. The event lies 0 .35 from the Galactic center and suffers from high extinction of AK=1.68. Therefore, it could be detected only by a near-infrared (NIR) survey. The field also suffers from large spatial differential extinction, which makes it difficult to estimate the source properties required to derive the angular Einstein radius. Nevertheless, we find evidence suggesting that the source is located in the far disk. If correct, this would be the first source star of a microlensing event to be identified as belonging to the far disk. We estimate the lens mass and distance using a Bayesian analysis to find that the planet's mass is 1.28+ M0.44 J 0.37 , and it orbits a-+ 0.81 +0.21-0.27 Ṁ 0.21 star at an instantaneous projected separation of 4.18+0.96-0.88au. The system is at a distance of-6.3+16-2.1kpc, and so likely resides in the Galactic bulge. In addition, we find a non-standard extinction curve in this field, in agreement with previous results toward high-extinction fields near the Galactic center.

We present an empirical approach for interpreting gravitational wave signals of binary black hole mergers under the assumption that the underlying black hole population is sourced by remnants of stellar evolution. Using the observed relationship between galaxy mass and stellar metallicity, we predict the black hole count as a function of galaxy stellar mass. We show, for example, that a galaxy like the Milky Way should host millions of ~30M⊙ black holes and dwarf satellite galaxies like Draco should host ~100 such remnants, with weak dependence on the assumed initial mass function and stellar evolution model. Most low-mass black holes (~10M⊙) typically reside within massive galaxies (M* ≃ 1011M⊙) while massive black holes (~50M⊙) typically reside within dwarf galaxies (M* ≃ 109M⊙) today. If roughly 1 per cent of black holes are involved in a binary black hole merger, then the reported merger rate densities from advanced Laser Interferometer Gravitational-Wave Observatory can be accommodated for a range of merger time-scales, and the detection of mergers with > 50M⊙ black holes should be expected within the next decade. Identifying the host galaxy population of the mergers provides a way to constrain both the binary neutron star or black hole formation efficiencies and the merger time-scale distributions; these events would be primarily localized in dwarf galaxies if the merger time-scale is short compared to the age of the Universe and in massive galaxies otherwise. As more mergers are detected, the prospect of identifying the host galaxy population, either directly through the detection of electromagnetic counterparts of binary neutron star mergers or indirectly through the anisotropy of the events, will become a realistic possibility.

We present an empirical approach for interpreting gravitational wave signals of binary black hole mergers under the assumption that the underlying black hole population is sourced by remnants of stellar evolution. Using the observed relationship between galaxy mass and stellar metallicity, we predict the black hole count as a function of galaxy stellar mass. We show, for example, that a galaxy like the Milky Way should host millions of ~30M⊙ black holes and dwarf satellite galaxies like Draco should host ~100 such remnants, with weak dependence on the assumed initial mass function and stellar evolution model. Most low-mass black holes (~10M⊙) typically reside within massive galaxies (M* ≃ 1011M⊙) while massive black holes (~50M⊙) typically reside within dwarf galaxies (M* ≃ 109M⊙) today. If roughly 1 per cent of black holes are involved in a binary black hole merger, then the reported merger rate densities from advanced Laser Interferometer Gravitational-Wave Observatory can be accommodated for a range of merger time-scales, and the detection of mergers with > 50M⊙ black holes should be expected within the next decade. Identifying the host galaxy population of the mergers provides a way to constrain both the binary neutron star or black hole formation efficiencies and the merger time-scale distributions; these events would be primarily localized in dwarf galaxies if the merger time-scale is short compared to the age of the Universe and in massive galaxies otherwise. As more mergers are detected, the prospect of identifying the host galaxy population, either directly through the detection of electromagnetic counterparts of binary neutron star mergers or indirectly through the anisotropy of the events, will become a realistic possibility.

We present new adaptive optics (AO) imaging and spectroscopic measurements of Galactic center source G1 from W. M. Keck Observatory. Our goal is to understand its nature and relationship to G2, which is the first example of a spatially resolved object interacting with a supermassive black hole (SMBH). Both objects have been monitored with AO for the past decade (2003-2014) and are comparatively close to the black hole (a min ∼ 200-300 au) on very eccentric orbits (e G1 ∼ 0.99; e G2 ∼ 0.96). While G2 has been tracked before and during periapsis passage (T 0 ∼ 2014.2), G1 has been followed since soon after emerging from periapsis (T 0 ∼ 2001.3). Our observations of G1 double the previously reported observational time baseline, which improves its orbital parameter determinations. G1's orbital trajectory appears to be in the same plane as that of G2 but with a significantly different argument of periapsis (Δω = 21 ± 4°). This suggests that G1 is an independent object and not part of a gas stream containing G2, as has been proposed. Furthermore, we show for the first time that (1) G1 is extended in the epochs closest to periapsis along the direction of orbital motion, and (2) it becomes significantly smaller over time (450 au in 2004 to less than 170 au in 2009). Based on these observations, G1 appears to be the second example of an object tidally interacting with an SMBH. G1's existence 14 yr after periapsis, along with its compactness in epochs further from the time of periapsis, suggest that this source is stellar in nature.

We present new adaptive optics (AO) imaging and spectroscopic measurements of Galactic center source G1 from W. M. Keck Observatory. Our goal is to understand its nature and relationship to G2, which is the first example of a spatially resolved object interacting with a supermassive black hole (SMBH). Both objects have been monitored with AO for the past decade (2003-2014) and are comparatively close to the black hole (a min ∼ 200-300 au) on very eccentric orbits (e G1 ∼ 0.99; e G2 ∼ 0.96). While G2 has been tracked before and during periapsis passage (T 0 ∼ 2014.2), G1 has been followed since soon after emerging from periapsis (T 0 ∼ 2001.3). Our observations of G1 double the previously reported observational time baseline, which improves its orbital parameter determinations. G1's orbital trajectory appears to be in the same plane as that of G2 but with a significantly different argument of periapsis (Δω = 21 ± 4°). This suggests that G1 is an independent object and not part of a gas stream containing G2, as has been proposed. Furthermore, we show for the first time that (1) G1 is extended in the epochs closest to periapsis along the direction of orbital motion, and (2) it becomes significantly smaller over time (450 au in 2004 to less than 170 au in 2009). Based on these observations, G1 appears to be the second example of an object tidally interacting with an SMBH. G1's existence 14 yr after periapsis, along with its compactness in epochs further from the time of periapsis, suggest that this source is stellar in nature.

The coalescence of a binary neutron star pair is expected to produce gravitational waves (GW) and electromagnetic radiation, both of which may be detectable with currently available instruments. We describe a search for a predicted r-process optical transient from these mergers, dubbed the "kilonova" (KN), using griz broadband data from the Dark Energy Survey Supernova Program (DES-SN). Some models predict KNe to be redder, shorter-lived, and dimmer than supernovae (SNe), but the event rate of KNe is poorly constrained. We simulate KN and SN light curves with the Monte-Carlo simulation code SNANA to optimize selection requirements, determine search efficiency, and predict SN backgrounds. Our analysis of the first two seasons of DES-SN data results in 0 events, and is consistent with our prediction of 1.1 ± 0.2 background events based on simulations of SNe. From our prediction, there is a 33% chance of finding 0 events in the data. Assuming no underlying galaxy flux, our search sets 90% upper limits on the KN volumetric rate of 1.0 Gpc-3 yr-1 for the dimmest KN model we consider (peak i-band absolute magnitude mag) and 2.4 Gpc-3 yr-1 for the brightest ( mag). Accounting for anomalous subtraction artifacts on bright galaxies, these limits are ∼3 times higher. This analysis is the first untriggered optical KN search and informs selection requirements and strategies for future KN searches. Our upper limits on the KN rate are consistent with those measured by GW and gamma-ray burst searches.

The coalescence of a binary neutron star pair is expected to produce gravitational waves (GW) and electromagnetic radiation, both of which may be detectable with currently available instruments. We describe a search for a predicted r-process optical transient from these mergers, dubbed the "kilonova" (KN), using griz broadband data from the Dark Energy Survey Supernova Program (DES-SN). Some models predict KNe to be redder, shorter-lived, and dimmer than supernovae (SNe), but the event rate of KNe is poorly constrained. We simulate KN and SN light curves with the Monte-Carlo simulation code SNANA to optimize selection requirements, determine search efficiency, and predict SN backgrounds. Our analysis of the first two seasons of DES-SN data results in 0 events, and is consistent with our prediction of 1.1 ± 0.2 background events based on simulations of SNe. From our prediction, there is a 33% chance of finding 0 events in the data. Assuming no underlying galaxy flux, our search sets 90% upper limits on the KN volumetric rate of 1.0 Gpc-3 yr-1 for the dimmest KN model we consider (peak i-band absolute magnitude mag) and 2.4 Gpc-3 yr-1 for the brightest ( mag). Accounting for anomalous subtraction artifacts on bright galaxies, these limits are ∼3 times higher. This analysis is the first untriggered optical KN search and informs selection requirements and strategies for future KN searches. Our upper limits on the KN rate are consistent with those measured by GW and gamma-ray burst searches.

We present a strongly interacting quadruple system associated with the K2 target EPIC 220204960. The K2 target itself is a Kp = 12.7-mag star at Teff ≃ 6100 K, which we designate as 'B-N' (blue northerly image). The host of the quadruple system, however, is a Kp ≃ 17-mag star with a composite M-star spectrum, which we designate as 'R-S' (red southerly image). With a 3.2-arcsec separation and similar radial velocities and photometric distances, 'B-N' is likely physically associated with 'R-S', making this a quintuple system, but that is incidental to our main claim of a strongly interacting quadruple system in 'R-S'. The two binaries in 'R-S' have orbital periods of 13.27 and 14.41 d, respectively, and each has an inclination angle of ≳89°. From our analysis of radial-velocity (RV) measurements, and of the photometric light curve, we conclude that all four stars are very similar with masses close to 0.4M⊙. Both of the binaries exhibit significant eclipse-timing variations where those of the primary and secondary eclipses 'diverge' by 0.05 d over the course of the 80-d observations. Via a systematic set of numerical simulations of quadruple systems consisting of two interacting binaries, we conclude that the outer orbital period is very likely to be between 300 and 500 d. If sufficient time is devoted to RV studies of this faint target, the outer orbit should be measurable within a year.

In this paper searches for flaring astrophysical neutrino sources and sources with periodic emission with the IceCube neutrino telescope are presented. In contrast to time-integrated searches, where steady emission is assumed, the analyses presented here look for a time-dependent signal of neutrinos using the information from the neutrino arrival times to enhance the discovery potential. A search was performed for correlations between neutrino arrival times and directions, as well as neutrino emission following time-dependent light curves, sporadic emission, or periodicities of candidate sources. These include active galactic nuclei, soft γ-ray repeaters, supernova remnants hosting pulsars, microquasars, and X-ray binaries. The work presented here updates and extends previously published results to a longer period that covers 4 years of data from 2008 April 5 to 2012 May 16, including the first year of operation of the completed 86 string detector. The analyses did not find any significant time-dependent point sources of neutrinos, and the results were used to set upper limits on the neutrino flux from source candidates.

Many close stellar binaries are accompanied by a faraway star. The "eccentric Kozai-Lidov" (EKL) mechanism can cause dramatic inclination and eccentricity fluctuations, resulting in tidal tightening of inner binaries of triple stars. We run a large set of Monte Carlo simulations, including the secular evolution of the orbits, general relativistic precession, and tides, and we determine the semimajor axis, eccentricity, inclination, and spin-orbit angle distributions of the final configurations. We find that the efficiency of forming tight binaries (≲ 16 days) when taking the EKL mechanism into account is 21%, and about 4% of all simulated systems ended up in a merger event. These merger events can lead to the formation of blue stragglers. Furthermore, we find that the spin-orbit angle distribution of the inner binaries carries a signature of the initial setup of the system; thus, observations can be used to disentangle close binaries' birth configuration. The resulting inner and outer final orbits' period distributions and their estimated fraction suggest that secular dynamics may be a significant channel for the formation of close binaries in triples and even blue stragglers.

Many close stellar binaries are accompanied by a faraway star. The "eccentric Kozai-Lidov" (EKL) mechanism can cause dramatic inclination and eccentricity fluctuations, resulting in tidal tightening of inner binaries of triple stars. We run a large set of Monte Carlo simulations, including the secular evolution of the orbits, general relativistic precession, and tides, and we determine the semimajor axis, eccentricity, inclination, and spin-orbit angle distributions of the final configurations. We find that the efficiency of forming tight binaries (≲ 16 days) when taking the EKL mechanism into account is 21%, and about 4% of all simulated systems ended up in a merger event. These merger events can lead to the formation of blue stragglers. Furthermore, we find that the spin-orbit angle distribution of the inner binaries carries a signature of the initial setup of the system; thus, observations can be used to disentangle close binaries' birth configuration. The resulting inner and outer final orbits' period distributions and their estimated fraction suggest that secular dynamics may be a significant channel for the formation of close binaries in triples and even blue stragglers.

We report a series of observations of the interstellar scintillation (ISS) of the double pulsar J0737-3039 over the course of 18 months. As in earlier work, the basic phenomenon is the variation in the ISS caused by the changing transverse velocities of each pulsar, the ionized interstellar medium (IISM), and the Earth. The transverse velocity of the binary system can be determined both by very long baseline interferometry and timing observations. The orbital velocity and inclination is almost completely determined from timing observations, but the direction of the orbital angular momentum is not known. Since the Earth's velocity is known, and can be compared with the orbital velocity by its effect on the timescale of the ISS, we can determine the orientation Ω of the pulsar orbit with respect to equatorial coordinates (Ω = 65 ± 2°). We also resolve the ambiguity (i = 88.°7 or 91.°3) in the inclination of the orbit deduced from the measured Shapiro delay by our estimate i = 88.°1 ± 0.°5. This relies on the analysis of the ISS over both frequency and time, and provides a model for the location, anisotropy, turbulence level, and transverse phase gradient of the IISM. We find that the IISM can be well-modeled during each observation, typically of a few orbital periods, but its turbulence level and mean velocity vary significantly over the 18 months. © 2014. The American Astronomical Society. All rights reserved.

We report a series of observations of the interstellar scintillation (ISS) of the double pulsar J0737-3039 over the course of 18 months. As in earlier work, the basic phenomenon is the variation in the ISS caused by the changing transverse velocities of each pulsar, the ionized interstellar medium (IISM), and the Earth. The transverse velocity of the binary system can be determined both by very long baseline interferometry and timing observations. The orbital velocity and inclination is almost completely determined from timing observations, but the direction of the orbital angular momentum is not known. Since the Earth's velocity is known, and can be compared with the orbital velocity by its effect on the timescale of the ISS, we can determine the orientation Ω of the pulsar orbit with respect to equatorial coordinates (Ω = 65 ± 2°). We also resolve the ambiguity (i = 88.°7 or 91.°3) in the inclination of the orbit deduced from the measured Shapiro delay by our estimate i = 88.°1 ± 0.°5. This relies on the analysis of the ISS over both frequency and time, and provides a model for the location, anisotropy, turbulence level, and transverse phase gradient of the IISM. We find that the IISM can be well-modeled during each observation, typically of a few orbital periods, but its turbulence level and mean velocity vary significantly over the 18 months. © 2014. The American Astronomical Society. All rights reserved.

In recent years, the number of pulsars with secure mass measurements has increased to a level that allows us to probe the underlying neutron star (NS) mass distribution in detail. We critically review the radio pulsar mass measurements. For the first time, we are able to analyze a sizable population of NSs with a flexible modeling approach that can effectively accommodate a skewed underlying distribution and asymmetric measurement errors. We find that NSs that have evolved through different evolutionary paths reflect distinctive signatures through dissimilar distribution peak and mass cutoff values. NSs in double NS and NS-white dwarf (WD) systems show consistent respective peaks at 1.33 M and 1.55 M, suggesting significant mass accretion (Δm 0.22 M) has occurred during the spin-up phase. The width of the mass distribution implied by double NS systems is indicative of a tight initial mass function while the inferred mass range is significantly wider for NSs that have gone through recycling. We find a mass cutoff at 2.1 M for NSs with WD companions, which establishes a firm lower bound for the maximum NS mass. This rules out the majority of strange quark and soft equation of state models as viable configurations for NS matter. The lack of truncation close to the maximum mass cutoff along with the skewed nature of the inferred mass distribution both enforce the suggestion that the 2.1 M limit is set by evolutionary constraints rather than nuclear physics or general relativity, and the existence of rare supermassive NSs is possible. © 2013. The American Astronomical Society. All rights reserved..

In recent years, the number of pulsars with secure mass measurements has increased to a level that allows us to probe the underlying neutron star (NS) mass distribution in detail. We critically review the radio pulsar mass measurements. For the first time, we are able to analyze a sizable population of NSs with a flexible modeling approach that can effectively accommodate a skewed underlying distribution and asymmetric measurement errors. We find that NSs that have evolved through different evolutionary paths reflect distinctive signatures through dissimilar distribution peak and mass cutoff values. NSs in double NS and NS-white dwarf (WD) systems show consistent respective peaks at 1.33 M and 1.55 M, suggesting significant mass accretion (Δm 0.22 M) has occurred during the spin-up phase. The width of the mass distribution implied by double NS systems is indicative of a tight initial mass function while the inferred mass range is significantly wider for NSs that have gone through recycling. We find a mass cutoff at 2.1 M for NSs with WD companions, which establishes a firm lower bound for the maximum NS mass. This rules out the majority of strange quark and soft equation of state models as viable configurations for NS matter. The lack of truncation close to the maximum mass cutoff along with the skewed nature of the inferred mass distribution both enforce the suggestion that the 2.1 M limit is set by evolutionary constraints rather than nuclear physics or general relativity, and the existence of rare supermassive NSs is possible. © 2013. The American Astronomical Society. All rights reserved..