Electromagnetic Counterparts to LIGO-VIRGO Events: Expanded Very Large Array Observations
No abstract available
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
No abstract available
Modern astronomers enjoy access to all-sky images across a wide range of the electromagnetic spectrum from long-wavelength radio to high-energy gamma rays. The most prominent feature in many of these images is our own Galaxy, with different features revealed in each wave band. Gravitational waves (GWs) have recently been added to the astronomers’ toolkit as a nonelectromagnetic messenger. To date, all identified GW sources have been extra-Galactic and transient. However, the Milky Way hosts a population of ultracompact binaries (UCBs), which radiate persistent GWs in the milliHertz band that is not observable with today’s terrestrial gravitational-wave detectors. Space-based detectors such as the Laser Interferometer Space Antenna will measure this population and provide a census of their location, masses, and orbital properties. In this work, we will show how this data can be used to form a false-color image of the Galaxy that represents the intensity and frequency of the gravitational waves produced by the UCB population. Such images can be used to study the morphology of the Galaxy, identify interesting multimessenger sources through cross-matching, and for educational and outreach purposes.
With the advent of modern neutrino and gravitational wave (GW) detectors, the promise of multimessenger detections of the next galactic core-collapse supernova (CCSN) has become very real. Such detections will give insight into the CCSN mechanism and the structure of the progenitor star, and may resolve longstanding questions in fundamental physics. In order to properly interpret these detections, a thorough understanding of the landscape of possible CCSN events, and their multimessenger signals, is needed. We present detailed predictions of neutrino and GW signals from 1D simulations of stellar core collapse, spanning the landscape of core-collapse progenitors from 9 to 120 M ⊙ . In order to achieve explosions in 1D, we use the Supernova Turbulence In Reduced-dimensionality model, which includes the effects of turbulence and convection in 1D supernova simulations to mimic the 3D explosion mechanism. We study the GW emission from the 1D simulations using an astroseismology analysis of the protoneutron star. We find that the neutrino and GW signals are strongly correlated with the structure of the progenitor star and remnant compact object. Using these correlations, future detections of the first few seconds of neutrino and GW emission from a galactic CCSN may be able to provide constraints on stellar evolution independent of preexplosion imaging and the mass of the compact object remnant prior to fallback accretion.
Upcoming LIGO–Virgo–KAGRA (LVK) observing runs are expected to detect a variety of inspiralling gravitational-wave (GW) events that come from black hole and neutron star binary mergers. Detection of noninspiral GW sources is also anticipated. We report the discovery of a new class of noninspiral GW sources—the end states of massive stars—that can produce the brightest simulated stochastic GW burst signal in the LVK bands known to date, and could be detectable in LVK run A+. Some dying massive stars launch bipolar relativistic jets, which inflate a turbulent energetic bubble—cocoon—inside of the star. We simulate such a system using state-of-the-art 3D general relativistic magnetohydrodynamic simulations and show that these cocoons emit quasi-isotropic GW emission in the LVK band, ~10–100 Hz, over a characteristic jet activity timescale ~10–100 s. Our first-principles simulations show that jets exhibit a wobbling behavior, in which case cocoon-powered GWs might be detected already in LVK run A+, but it is more likely that these GWs will be detected by the third-generation GW detectors with an estimated rate of ~10 events yr –1 . The detection rate drops to ~1% of that value if all jets were to feature a traditional axisymmetric structure instead of a wobble. Accompanied by electromagnetic emission from the energetic core-collapse supernova and the cocoon, we predict that collapsars are powerful multimessenger events.
We present a comprehensive analysis of 653 optical candidate counterparts reported during the third gravitational-wave (GW) observing run. Our sample concentrates on candidates from the 15 events (published in GWTC-2, GWTC-3, or not retracted on GraceDB) that had a >1% chance of including a neutron star in order to assess their viability as true kilonovae. In particular, we leverage tools available in real time, including pre-merger detections and cross-matching with catalogs (i.e., point-source, variable-star, quasar and host-galaxy redshift data sets), to eliminate 65% of candidates in our sample. We further employ spectroscopic classifications, late-time detections, and light-curve behavior analyses and conclude that 66 candidates remain viable kilonovae. These candidates lack sufficient information to determine their classifications, and the majority would require luminosities greater than that of AT 2017gfo. Pre-merger detections in public photometric survey data and comparison of cataloged host-galaxy redshifts with the GW event distances are critical to incorporate into vetting procedures, as these tools eliminated >20% and >30% of candidates, respectively. We expect that such tools that leverage archival information will significantly reduce the strain on spectroscopic and photometric follow-up resources in future observing runs. Finally, we discuss the critical role prompt updates from GW astronomers to the EM community play in reducing the number of candidates requiring vetting.
Not Available
Numerical simulations of merging compact objects and their remnants form the theoretical foundation for gravitational wave and multimessenger astronomy. While Cartesian-coordinate-based adaptive mesh refinement is commonly used for simulations, spherical-like coordinates are more suitable for nearly spherical remnants and azimuthal flows due to lower numerical dissipation in the evolution of fluid angular momentum, as well as requiring fewer numbers of computational cells. However, the use of spherical coordinates to numerically solve hyperbolic partial differential equations can result in severe Courant-Friedrichs-Lewy (CFL) stability condition time step limitations, which can make simulations prohibitively expensive. This paper addresses this issue for the numerical solution of coupled spacetime and general relativistic magnetohydrodynamics evolutions by introducing a double fast Fourier transform (FFT) filter and implementing it within the fully message passing interface (mpi)-parallelized sphericalnr framework in the einstein toolkit. In conclusion, we demonstrate the effectiveness and robustness of the filtering algorithm by applying it to a number of challenging code tests, and show that it passes these tests effectively, demonstrating convergence while also increasing the time step significantly compared to unfiltered simulations.
Numerical relativity is central to the investigation of astrophysical sources in the dynamical and strong-field gravity regime, such as binary black hole and neutron star coalescences. Current challenges set by gravitational-wave and multimessenger astronomy call for highly performant and scalable codes on modern massively parallel architectures. We present GR-Athena++, a general-relativistic, high-order, vertex-centered solver that extends the oct-tree, adaptive mesh refinement capabilities of the astrophysical (radiation) magnetohydrodynamics code Athena++. To simulate dynamical spacetimes, GR-Athena++ uses the Z4c evolution scheme of numerical relativity coupled to the moving puncture gauge. We demonstrate stable and accurate binary black hole merger evolutions via extensive convergence testing, cross-code validation, and verification against state-of-the-art effective-one-body waveforms. GR-Athena++ leverages the task-based parallelism paradigm of Athena++ to achieve excellent scalability. We measure strong-scaling efficiencies above 95% for up to ~1.2 × 10 4 CPUs and excellent weak scaling is shown up to ~10 5 CPUs in a production binary black hole setup with adaptive mesh refinement. GR-Athena++ thus allows for the robust simulation of compact binary coalescences and offers a viable path toward numerical relativity at exascale.
Massive black hole mergers are perhaps the most energetic astronomical events, establishing their importance as gravitational wave sources for LISA, and also possibly leading to observable influences on their local environments. Advances in numerical relativity over the last five years have fueled the development of a rich physical understanding of general relativity's predictions for these events. Z will overview the understanding of these event emerging from numerical simulation studies. These simulations elucidate the pre-merger dynamics of the black hole binaries, the consequent gravitational waveform signatures ' and the resulting state, including its kick velocity, for the final black hole produced by the merger. Scenarios are now being considered for observing each of these aspects of the merger, involving both gravitational-wave and electromagnetic astronomy.
In this work, we consider the properties of compact stars in which quark matter has low- and high-density phases that are separated by a first-order phase transition. Thus, unlike the commonly considered case of a single phase transition from hadronic to quark matter, our models of hybrid stars contain sequential phase transitions from hadronic matter to low- and then to high-density quark matter phases. We extend our previous study of the parameter space of hybrid stars with a single phase transition to those with sequential phase transitions, taking into account the constraints on the mass and radius of neutron stars from the NICER experiment, the experimental inferences of the neutron skin thickness of the lead nucleus by the PREX-II experiment, and constraints on the tidal deformability from the gravitational-wave event GW170817. We determine the range of the masses for which both twin and triplet configurations, i.e., identical-mass stars with two and three different values of radii, arise.
Neutron star mergers have recently become a tool to study extreme gravity, nucleosynthesis, and the chemical composition of the Universe. To date, there has been one joint gravitational and electromagnetic observation of a binary neutron star merger, GW170817, as well as a solely gravitational observation, GW190425. In order to accurately identify and interpret electromagnetic signals of neutron star mergers, better models of the matter outflows generated by these mergers are required. We compare a series of ejecta models to see where they provide strong constraints on the amount of ejected mass expected from a system, and where systematic uncertainties in current models prevent us from reliably extracting information from observed events. We also examine 2396 neutron star equations of state compatible with GW170817 to see whether a given ejecta mass could be reasonably produced with a neutron star of said equation of state, and whether different ejecta models provide consistent predictions. We find that the difference between models is often comparable to or larger than the error generally assumed for these models, implying better constraints on the models are needed. In conclusion, we also note that the extrapolation of outflow models outside of their calibration window, while commonly needed to analyze gravitational wave events, is extremely unreliable and occasionally leads to completely unphysical results.
The LIGO HET Response (LIGHETR) project is an enterprise to follow up optical transients (OTs) discovered as gravitational-wave merger sources by the LIGO/Virgo collaboration (LVC). Early spectroscopy has the potential to constrain crucial parameters such as the aspect angle. The LIGHETR collaboration also includes the capacity to model the spectroscopic evolution of mergers to facilitate a real-time direct comparison of models with our data. The principal facility is the Hobby–Eberly Telescope. LIGHETR uses the massively replicated VIRUS array of spectrographs to search for associated OTs and obtain early blue spectra, and in a complementary role, the low-resolution LRS2 spectrograph is used to obtain spectra of viable candidates as well as a densely sampled series of spectra of true counterparts. Once an OT is identified, the anticipated cadence of spectra would match or considerably exceed anything achieved for GW170817 = AT2017gfo for which there were no spectra in the first 12 hr and thereafter only roughly once daily. We describe special HET-specific software written to facilitate the program and attempts to determine the flux limits to undetected sources. We also describe our campaign to follow up OT candidates during the third observational campaign of the LIGO and Virgo Scientific Collaborations. We obtained VIRUS spectroscopy of candidate galaxy hosts for five LVC gravitational-wave events and LRS2 spectra of one candidate for the OT associated with S190901ap. We identified that candidate, ZTF19abvionh = AT2019pip, as a possible Wolf–Rayet star in an otherwise unrecognized nearby dwarf galaxy.
Abstract Kilonovae are ultraviolet, optical, and infrared transients powered by the radioactive decay of heavy elements following a neutron star merger. Joint observations of kilonovae and gravitational waves can offer key constraints on the source of Galactic r -process enrichment, among other astrophysical topics. However, robust constraints on heavy element production require rapid kilonova detection (within ∼1 day of merger) as well as multiwavelength observations across multiple epochs. In this study, we quantify the ability of 13 wide-field-of-view instruments to detect kilonovae, leveraging a large grid of over 900 radiative transfer simulations with 54 viewing angles per simulation. We consider both current and upcoming instruments, collectively spanning the full kilonova spectrum. The Roman Space Telescope has the highest redshift reach of any instrument in the study, observing kilonovae out to z ∼ 1 within the first day post-merger. We demonstrate that BlackGEM, DECam, GOTO, the Vera C. Rubin Observatory’s LSST, ULTRASAT, VISTA, and WINTER can observe some kilonovae out to z ∼ 0.1 (∼475 Mpc), while DDOTI, MeerLICHT, PRIME, Swift/UVOT, and ZTF are confined to more nearby observations. Furthermore, we provide a framework to infer kilonova ejecta properties following nondetections and explore variation in detectability with these ejecta parameters.
ABSTRACT Non-axisymmetrical deformations of the crust on rapidly rotating neutron stars are one of the main targets of searches for continuous gravitational waves. The maximum ellipticity, or fractional difference in moments of inertia, which can be supported by deformations of the crust (known as ‘mountains’), provides an important upper limit on the strength of these continuous gravitational wave sources. We use the formalism developed by Gittins and Andersson, along with a deforming force that acts mainly in the transverse direction, to obtain a maximum ellipticity of 7.4 × 10−6. This is larger than the original results that Gittins and Andersson obtained but consistent with earlier calculations by Ushomirsky, Cutler, and Bildsten. This suggests that rotating neutron stars could be strong sources of continuous gravitational waves.
For velocity and magnetic fields, the turbulent pressure and, more generally, the squared fields such as the components of the turbulent stress tensor, play important roles in astrophysics. For both one and three dimensions, we derive the equations relating the energy spectra of the fields to the spectra of their squares. We solve the resulting integrals numerically and show that for turbulent energy spectra of Kolmogorov type, the spectral slope of the stress spectrum is also of Kolmogorov type. For shallower turbulence spectra, the slope of the stress spectrum quickly approaches that of white noise, regardless of how blue the spectrum of the field is. For fully helical fields, the stress spectrum is elevated by about a factor of two in the subinertial range, while that in the inertial range remains unchanged. We discuss possible implications for understanding the spectrum of primordial gravitational waves from causally generated magnetic fields during cosmological phase transitions in the early universe. We also discuss potential diagnostic applications to the interstellar medium, where polarization and scintillation measurements characterize the square of the magnetic field.
With ground-based interferometers detecting hundreds of gravitational-wave (GW) events, GW astronomy has continued to blossom. U.S. and European scientists are developing plans to construct third-generation ground-based interferometers. ESA has proceeded through the mission formulation phase of space-based interferometer in a lower-frequency band, 10 –4 –0.1 Hz. Despite all these exciting developments, there is still a missing frequency band, 0.1–10 Hz. This mid-frequency band is rich with interesting astrophysical events. Coalescence and merger of intermediate-mass black holes (IMBHs) will occur in this frequency band. Coalescing stellar-mass BHs will pass through this frequency band days before they reach the frequency band of LIGO and Virgo. Detection of such signals would enable a mid-frequency detector to issue an advance notice to the high-frequency GW detectors, as well as to optical, x-ray and γ-ray telescopes. We propose Lunar Accelerometer Network Gravitational Observatory (LANGO) to detect GWs in this frequency band. In the first phase (LANGO 1), we propose to deploy four ambient-temperature (250 K) accelerometers in the tetrahedral or in a square configuration on the hemisphere facing the Earth. After successful operation at 250 K, LANGO would be upgraded to a cryogenic (4 K) version with over two orders of magnitude increased sensitivity (LANGO 2), with coherent rejection of seismic noise implemented. LANGO is a full-tensor detector, capable of determining the source direction and wave polarization. Each test mass (TM) is suspended as a pendulum with resonance frequency ∼ 0.01 Hz, thus is only weakly coupled to the lunar surface horizontally. LANGO is designed to detect the relative motion of globally separated, nearly free, TMs by using the Moon as a large quiet platform. LANGO 1 and 2 accelerometers aim at sensitivities ⩽ 10 –11 m s –2 Hz –1/2 and ⩽ 10 –13 m s –2 Hz –1/2 in the horizontal axes over the frequency band of 1 mHz–10 Hz, which yield GW sensitivities 1.1 x 10 -21 Hz -1/2 and 3.8 x 10 -24 Hz -1/2 at 1 Hz, respectively. The LANGO accelerometers will be 10 3 –10 5 times more sensitive than Apollo seismometers. With such sensitivity, LANGO will also make great contribution to the advancement of lunar geophysics.
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 ℎ 0 are calculated. The best upper limits are obtained in the frequency range of 100 to 200 Hz and they are ~1.1 ×10 -25 at 95% confidence level. The minimum upper limit of 1.10 ×10 -25 is achieved at a frequency 111.5 Hz. In conclusion, 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.
During the inspiralling of a white dwarf (WD) into an intermediate-mass black hole (~10 2-5 M ⊙ ), both gravitational waves (GWs) and electromagnetic (EM) radiation are emitted. Once the eccentric orbit's pericenter radius approaches the tidal radius, the WD would be tidally stripped upon each pericenter passage. The accretion of this stripped mass would produce EM radiation. It is suspected that the recently discovered new types of transients, namely the quasiperiodic eruptions and the fast ultraluminous x-ray bursts, might originate from such systems. Modeling these flares requires a prediction of the amount of stripped mass from the WD and the details of the mass supply to the accretion disk. We run hydrodynamical simulations to study the orbital parameter dependence of the stripped mass. We find that our results match the analytical estimate that the stripped mass is proportional to z 5/2 , where z is the excess depth by which the WD overfills its instantaneous Roche lobe at the pericenter. The corresponding fallback rate of the stripped mass is calculated, which may be useful in interpreting the individual flaring light curve in candidate EM sources. We further calculate the long-term mass-loss evolution of a WD during its inspiral and the detectability of the GW and EM signals. The EM signal from the mass-loss stage can be easily detected: the limiting distance is ~320(M h /10 4 M ⊙ ) Mpc for the Einstein Probe. The GW signal, for space-borne detectors such as Laser Interferometer Space Antenna or TianQin, can be detected only within the Local Supercluster (~33 Mpc).