Showing 1–50 of 136 results for author: Haas, R

AsterX: a new open-source GPU-accelerated GRMHD code for dynamical spacetimes

Authors: Jay V. Kalinani, Liwei Ji, Lorenzo Ennoggi, Federico G. Lopez Armengol, Lucas Timotheo Sanches, Bing-Jyun Tsao, Steven R Brandt, Manuela Campanelli, Riccardo Ciolfi, Bruno Giacomazzo, Roland Haas, Erik Schnetter, Yosef Zlochower

Abstract: We present AsterX, a novel open-source, GPU-accelerated, fully general relativistic magnetohydrodynamic (GRMHD) code designed for dynamic spacetimes in 3D Cartesian coordinates, and tailored for exascale computing. We utilize block-structured adaptive mesh refinement (AMR) through CarpetX, the new driver for the Einstein Toolkit, which is built on AMReX, a software framework for massively parallel… ▽ More We present AsterX, a novel open-source, GPU-accelerated, fully general relativistic magnetohydrodynamic (GRMHD) code designed for dynamic spacetimes in 3D Cartesian coordinates, and tailored for exascale computing. We utilize block-structured adaptive mesh refinement (AMR) through CarpetX, the new driver for the Einstein Toolkit, which is built on AMReX, a software framework for massively parallel applications. AsterX employs the Valencia formulation for GRMHD, coupled with the 'Z4c' formalism for spacetime evolution, while incorporating high resolution shock capturing schemes to accurately handle the hydrodynamics. AsterX has undergone rigorous testing in both static and dynamic spacetime, demonstrating remarkable accuracy and agreement with other codes in literature. Benchmarking the code through scaling tests on OLCF's Frontier supercomputer, we demonstrate a weak scaling efficiency of about 67%-77% on 4096 nodes compared to an 8-node performance. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 41 pages, 12 figures, to be submitted to Classical and Quantum Gravity

General Relativistic Stability and Gravitational Wave Content of Rotating Triaxial Neutron Stars

Authors: Yufeng Luo, Antonios Tsokaros, Roland Haas, Koji Uryu

Abstract: Triaxial neutron stars can be sources of continuous gravitational radiation detectable by ground-based interferometers. The amplitude of the emitted gravitational wave can be greatly affected by the state of the hydrodynamical fluid flow inside the neutron star. In this work we examine the most triaxial models along two sequences of constant rest mass, confirming their dynamical stability. We also… ▽ More Triaxial neutron stars can be sources of continuous gravitational radiation detectable by ground-based interferometers. The amplitude of the emitted gravitational wave can be greatly affected by the state of the hydrodynamical fluid flow inside the neutron star. In this work we examine the most triaxial models along two sequences of constant rest mass, confirming their dynamical stability. We also study the response of a triaxial figure of quasiequilibrium under a variety of perturbations that lead to different fluid flows. Starting from the general relativistic compressible analog of the Newtonian Jacobi ellipsoid, we perform simulations of Dedekind-type flows. We find that in some cases the triaxial neutron star resembles a Riemann-S-type ellipsoid with minor rotation and gravitational wave emission as it evolves towards axisymmetry. The present results highlight the importance of understanding the fluid flow in the interior of a neutron star in terms of its gravitational wave content. △ Less

Submitted 27 December, 2023; originally announced December 2023.

Comments: 16 pages, 5 figs

Journal ref: Symmetry 2024, 16(3), 343

Waveform Modelling for the Laser Interferometer Space Antenna

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

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

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

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

Outflows from Short-Lived Neutron-Star Merger Remnants Can Produce a Blue Kilonova

Authors: Sanjana Curtis, Pablo Bosch, Philipp Mösta, David Radice, Sebastiano Bernuzzi, Albino Perego, Roland Haas, Erik Schnetter

Abstract: We present a 3D general-relativistic magnetohydrodynamic simulation of a short-lived neutron star remnant formed in the aftermath of a binary neutron star merger. The simulation uses an M1 neutrino transport scheme to track neutrino-matter interactions and is well-suited to studying the resulting nucleosynthesis and kilonova emission. We find that the ejecta in our simulations under-produce $r$-pr… ▽ More We present a 3D general-relativistic magnetohydrodynamic simulation of a short-lived neutron star remnant formed in the aftermath of a binary neutron star merger. The simulation uses an M1 neutrino transport scheme to track neutrino-matter interactions and is well-suited to studying the resulting nucleosynthesis and kilonova emission. We find that the ejecta in our simulations under-produce $r$-process abundances beyond the second $r$-process peak. For sufficiently long-lived remnants, these outflows \textit{alone} can produce blue kilonovae, including the blue kilonova component observed for AT2017gfo. △ Less

Submitted 12 May, 2023; originally announced May 2023.

Comments: 9 pages, 5 figures, submitted to ApJL

Theoretical and Experimental Constraints for the Equation of State of Dense and Hot Matter

Authors: Rajesh Kumar, Veronica Dexheimer, Johannes Jahan, Jorge Noronha, Jacquelyn Noronha-Hostler, Claudia Ratti, Nico Yunes, Angel Rodrigo Nava Acuna, Mark Alford, Mahmudul Hasan Anik, Debarati Chatterjee, Katerina Chatziioannou, Hsin-Yu Chen, Alexander Clevinger, Carlos Conde, Nikolas Cruz-Camacho, Travis Dore, Christian Drischler, Hannah Elfner, Reed Essick, David Friedenberg, Suprovo Ghosh, Joaquin Grefa, Roland Haas, Alexander Haber , et al. (35 additional authors not shown)

Abstract: This review aims at providing an extensive discussion of modern constraints relevant for dense and hot strongly interacting matter. It includes theoretical first-principle results from lattice and perturbative QCD, as well as chiral effective field theory results. From the experimental side, it includes heavy-ion collision and low-energy nuclear physics results, as well as observations from neutro… ▽ More This review aims at providing an extensive discussion of modern constraints relevant for dense and hot strongly interacting matter. It includes theoretical first-principle results from lattice and perturbative QCD, as well as chiral effective field theory results. From the experimental side, it includes heavy-ion collision and low-energy nuclear physics results, as well as observations from neutron stars and their mergers. The validity of different constraints, concerning specific conditions and ranges of applicability, is also provided. △ Less

Submitted 12 June, 2024; v1 submitted 29 March, 2023; originally announced March 2023.

Journal ref: Living Reviews in Relativity (2024) 27:3

A monitoring campaign (2013-2020) of ESA's Mars Express to study interplanetary plasma scintillation

Authors: P. Kummamuru, G. Molera Calvés, G. Cimò, S. V. Pogrebenko, T. M. Bocanegra-Bahamón, D. A. Duev, M. D. Md Said, J. Edwards, M. Ma, J. Quick, A. Neidhardt, P. de Vicente, R. Haas, J. Kallunki, 1 G. Maccaferri, G. Colucci, W. J. Yang, L. F. Hao, S. Weston, M. A. Kharinov, A. G. Mikhailov, T. Jung

Abstract: The radio signal transmitted by the Mars Express (MEX) spacecraft was observed regularly between the years 2013-2020 at X-band (8.42 GHz) using the European Very Long Baseline Interferometry (EVN) network and University of Tasmania's telescopes. We present a method to describe the solar wind parameters by quantifying the effects of plasma on our radio signal. In doing so, we identify all the uncom… ▽ More The radio signal transmitted by the Mars Express (MEX) spacecraft was observed regularly between the years 2013-2020 at X-band (8.42 GHz) using the European Very Long Baseline Interferometry (EVN) network and University of Tasmania's telescopes. We present a method to describe the solar wind parameters by quantifying the effects of plasma on our radio signal. In doing so, we identify all the uncompensated effects on the radio signal and see which coronal processes drive them. From a technical standpoint, quantifying the effect of the plasma on the radio signal helps phase referencing for precision spacecraft tracking. The phase fluctuation of the signal was determined for Mars' orbit for solar elongation angles from 0 - 180 deg. The calculated phase residuals allow determination of the phase power spectrum. The total electron content (TEC) of the solar plasma along the line of sight is calculated by removing effects from mechanical and ionospheric noises. The spectral index was determined as $-2.43 \pm 0.11$ which is in agreement with Kolomogorov's turbulence. The theoretical models are consistent with observations at lower solar elongations however at higher solar elongation ($>$160 deg) we see the observed values to be higher. This can be caused when the uplink and downlink signals are positively correlated as a result of passing through identical plasma sheets. △ Less

Submitted 20 February, 2023; originally announced February 2023.

Comments: The paper has 13 figures and one table. It has been accepted for publication in PASA and the article will receive its DOI in a week's time

GRaM-X: A new GPU-accelerated dynamical spacetime GRMHD code for Exascale computing with the Einstein Toolkit

Authors: Swapnil Shankar, Philipp Mösta, Steven R. Brandt, Roland Haas, Erik Schnetter, Yannick de Graaf

Abstract: We present GRaM-X (General Relativistic accelerated Magnetohydrodynamics on AMReX), a new GPU-accelerated dynamical-spacetime general relativistic magnetohydrodynamics (GRMHD) code which extends the GRMHD capability of Einstein Toolkit to GPU-based exascale systems. GRaM-X supports 3D adaptive mesh refinement (AMR) on GPUs via a new AMR driver for the Einstein Toolkit called CarpetX which in turn… ▽ More We present GRaM-X (General Relativistic accelerated Magnetohydrodynamics on AMReX), a new GPU-accelerated dynamical-spacetime general relativistic magnetohydrodynamics (GRMHD) code which extends the GRMHD capability of Einstein Toolkit to GPU-based exascale systems. GRaM-X supports 3D adaptive mesh refinement (AMR) on GPUs via a new AMR driver for the Einstein Toolkit called CarpetX which in turn leverages AMReX, an AMR library developed for use by the United States DOE's Exascale Computing Project (ECP). We use the Z4c formalism to evolve the equations of GR and the Valencia formulation to evolve the equations of GRMHD. GRaM-X supports both analytic as well as tabulated equations of state. We implement TVD and WENO reconstruction methods as well as the HLLE Riemann solver. We test the accuracy of the code using a range of tests on static spacetime, e.g. 1D MHD shocktubes, the 2D magnetic rotor and a cylindrical explosion, as well as on dynamical spacetimes, i.e. the oscillations of a 3D TOV star. We find excellent agreement with analytic results and results of other codes reported in literature. We also perform scaling tests and find that GRaM-X shows a weak scaling efficiency of $\sim 40-50\%$ on 2304 nodes (13824 NVIDIA V100 GPUs) with respect to single-node performance on OLCF's supercomputer Summit. △ Less

Submitted 21 November, 2022; v1 submitted 31 October, 2022; originally announced October 2022.

Comments: 22 pages, 8 figures, to be submitted to Classical and Quantum Gravity

General relativistic simulations of collapsing binary neutron star mergers with Monte-Carlo neutrino transport

Authors: Francois Foucart, Matthew D. Duez, Roland Haas, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheel, Elizabeth Spira-Savett

Abstract: Recent gravitational wave observations of neutron star-neutron star and neutron star-black hole binaries appear to indicate that massive neutron stars may not be too uncommon in merging systems. In this manuscript, we present a first set of evolution of massive neutron star binaries using Monte-Carlo radiation transport for the evolution of neutrinos. We study a range of systems, from nearly symme… ▽ More Recent gravitational wave observations of neutron star-neutron star and neutron star-black hole binaries appear to indicate that massive neutron stars may not be too uncommon in merging systems. In this manuscript, we present a first set of evolution of massive neutron star binaries using Monte-Carlo radiation transport for the evolution of neutrinos. We study a range of systems, from nearly symmetric binaries that collapse to a black hole before forming a disk or ejecting material, to more asymmetric binaries in which tidal disruption of the lower mass star leads to the production of more interesting post-merger remnants. For the latter type of systems, we additionally study the impact of viscosity on the properties of the outflows, and compare our results to two recent simulations of identical binaries performed with the WhiskyTHC code. We find agreement on the black hole properties, disk mass, and mass and velocity of the outflows within expected numerical uncertainties, and some minor but noticeable differences in the evolution of the electron fraction when using a subgrid viscosity model, with viscosity playing a more minor role in our simulations. The method used to account for r-process heating in the determination of the outflow properties appears to have a larger impact on our result than those differences between numerical codes. We also use the simulation with the most ejected material to verify that our newly implemented Lagrangian tracers provide a reasonable sampling of the matter outflows as they leave the computational grid. We note that, given the lack of production of hot outflows in these mergers, the main role of neutrinos in these systems is to set the composition of the post-merger remnant. One of the main potential use of our simulations is thus as improved initial conditions for longer evolutions of such remnants. △ Less

Submitted 12 April, 2023; v1 submitted 11 October, 2022; originally announced October 2022.

Comments: 21p, 10 figures, under review by PRD

GENESIS: Co-location of Geodetic Techniques in Space

Authors: Pacôme Delva, Zuheir Altamimi, Alejandro Blazquez, Mathis Blossfeld, Johannes Böhm, Pascal Bonnefond, Jean-Paul Boy, Sean Bruinsma, Grzegorz Bury, Miltiadis Chatzinikos, Alexandre Couhert, Clément Courde, Rolf Dach, Véronique Dehant, Simone Dell'Agnello, Gunnar Elgered, Werner Enderle, Pierre Exertier, Susanne Glaser, Rüdiger Haas, Wen Huang, Urs Hugentobler, Adrian Jäggi, Ozgur Karatekin, Frank G. Lemoine , et al. (17 additional authors not shown)

Abstract: Improving and homogenizing time and space reference systems on Earth and, more directly, realizing the Terrestrial Reference Frame (TRF) with an accuracy of 1mm and a long-term stability of 0.1mm/year are relevant for many scientific and societal endeavors. The knowledge of the TRF is fundamental for Earth and navigation sciences. For instance, quantifying sea level change strongly depends on an a… ▽ More Improving and homogenizing time and space reference systems on Earth and, more directly, realizing the Terrestrial Reference Frame (TRF) with an accuracy of 1mm and a long-term stability of 0.1mm/year are relevant for many scientific and societal endeavors. The knowledge of the TRF is fundamental for Earth and navigation sciences. For instance, quantifying sea level change strongly depends on an accurate determination of the geocenter motion but also of the positions of continental and island reference stations, as well as the ground stations of tracking networks. Also, numerous applications in geophysics require absolute millimeter precision from the reference frame, as for example monitoring tectonic motion or crustal deformation for predicting natural hazards. The TRF accuracy to be achieved represents the consensus of various authorities which has enunciated geodesy requirements for Earth sciences. Today we are still far from these ambitious accuracy and stability goals for the realization of the TRF. However, a combination and co-location of all four space geodetic techniques on one satellite platform can significantly contribute to achieving these goals. This is the purpose of the GENESIS mission, proposed as a component of the FutureNAV program of the European Space Agency. The GENESIS platform will be a dynamic space geodetic observatory carrying all the geodetic instruments referenced to one another through carefully calibrated space ties. The co-location of the techniques in space will solve the inconsistencies and biases between the different geodetic techniques in order to reach the TRF accuracy and stability goals endorsed by the various international authorities and the scientific community. The purpose of this white paper is to review the state-of-the-art and explain the benefits of the GENESIS mission in Earth sciences, navigation sciences and metrology. △ Less

Submitted 30 September, 2022; originally announced September 2022.

Comments: 31 pages, 9 figures, submitted to Earth, Planets and Space (EPS)

Addition of tabulated equation of state and neutrino leakage support to IllinoisGRMHD

Authors: Leonardo R. Werneck, Zachariah B. Etienne, Ariadna Murguia-Berthier, Roland Haas, Federico Cipolletta, Scott C. Noble, Lorenzo Ennoggi, Federico G. Lopez Armengol, Bruno Giacomazzo, Thiago Assumpção, Joshua Faber, Tanmayee Gupte, Bernard J. Kelly, Julian H. Krolik

Abstract: We have added support for realistic, microphysical, finite-temperature equations of state (EOS) and neutrino physics via a leakage scheme to IllinoisGRMHD, an open-source GRMHD code for dynamical spacetimes in the Einstein Toolkit. These new features are provided by two new, NRPy+-based codes: NRPyEOS, which performs highly efficient EOS table lookups and interpolations, and NRPyLeakage, which imp… ▽ More We have added support for realistic, microphysical, finite-temperature equations of state (EOS) and neutrino physics via a leakage scheme to IllinoisGRMHD, an open-source GRMHD code for dynamical spacetimes in the Einstein Toolkit. These new features are provided by two new, NRPy+-based codes: NRPyEOS, which performs highly efficient EOS table lookups and interpolations, and NRPyLeakage, which implements a new, AMR-capable neutrino leakage scheme in the Einstein Toolkit. We have performed a series of strenuous validation tests that demonstrate the robustness of these new codes, particularly on the Cartesian AMR grids provided by Carpet. Furthermore, we show results from fully dynamical GRMHD simulations of single unmagnetized neutron stars, and magnetized binary neutron star mergers. This new version of IllinoisGRMHD, as well as NRPyEOS and NRPyLeakage, is pedagogically documented in Jupyter notebooks and fully open source. The codes will be proposed for inclusion in an upcoming version of the Einstein Toolkit. △ Less

Submitted 14 December, 2023; v1 submitted 30 August, 2022; originally announced August 2022.

Comments: 20 pages, 9 figures. v3 erratum: fixed typo in Eq. (21)

Journal ref: Phys. Rev. D 107, 044037 (2023)

Kepler K2 Campaign 9: II. First space-based discovery of an exoplanet using microlensing

Authors: D. Specht, R. Poleski, M. T. Penny, E. Kerins, I. McDonald, Chung-Uk Lee, A. Udalski, I. A. Bond, Y. Shvartzvald, Weicheng Zang, R. A. Street, D. W. Hogg, B. S. Gaudi, T. Barclay, G. Barentsen, S. B. Howell, F. Mullally, C. B. Henderson, S. T. Bryson, D. A. Caldwell, M. R. Haas, J. E. Van Cleve, K. Larson, K. McCalmont, C. Peterson , et al. (61 additional authors not shown)

Abstract: We present K2-2016-BLG-0005Lb, a densely sampled, planetary binary caustic-crossing microlensing event found from a blind search of data gathered from Campaign 9 of the Kepler K2 mission (K2C9). K2-2016-BLG-0005Lb is the first bound microlensing exoplanet discovered from space-based data. The event has caustic entry and exit points that are resolved in the K2C9 data, enabling the lens--source rela… ▽ More We present K2-2016-BLG-0005Lb, a densely sampled, planetary binary caustic-crossing microlensing event found from a blind search of data gathered from Campaign 9 of the Kepler K2 mission (K2C9). K2-2016-BLG-0005Lb is the first bound microlensing exoplanet discovered from space-based data. The event has caustic entry and exit points that are resolved in the K2C9 data, enabling the lens--source relative proper motion to be measured. We have fitted a binary microlens model to the Kepler data, and to simultaneous observations from multiple ground-based surveys. Whilst the ground-based data only sparsely sample the binary caustic, they provide a clear detection of parallax that allows us to break completely the microlensing mass--position--velocity degeneracy and measure the planet's mass directly. We find a host mass of $0.58\pm0.04 ~{\rm M}_\odot$ and a planetary mass of $1.1\pm0.1 ~{\rm M_J}$. The system lies at a distance of $5.2\pm0.2~$kpc from Earth towards the Galactic bulge, more than twice the distance of the previous most distant planet found by Kepler. The sky-projected separation of the planet from its host is found to be $4.2\pm0.3~$au which, for circular orbits, deprojects to a host separation $a = 4.4^{+1.9}_{-0.4}~$au and orbital period $P = 13^{+9}_{-2}~$yr. This makes K2-2016-BLG-0005Lb a close Jupiter analogue orbiting a low-mass host star. According to current planet formation models, this system is very close to the host mass threshold below which Jupiters are not expected to form. Upcoming space-based exoplanet microlensing surveys by NASA's Nancy Grace Roman Space Telescope and, possibly, ESA's Euclid mission, will provide demanding tests of current planet formation models. △ Less

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

Comments: 18 pages. Accepted for publication in MNRAS Published by Oxford University Press on behalf of the Royal Astronomical Society

Broad band flux-density monitoring of radio sources with the Onsala twin telescopes

Authors: E. Varenius, F. Maio, K. Le Bail, R. Haas

Abstract: Context and aims: The Onsala twin telescopes (OTT) are two 13 m telescopes located at the Onsala Space Observatory in Sweden. With dual linear polarized broad-band (3-14 GHz) receivers, they are part of the next generation Very Long Baseline Interferometry (VLBI) Global Observing System (VGOS) for geodesy and astrometry. In addition to purely geodetic data products, VGOS will regularly produce ful… ▽ More Context and aims: The Onsala twin telescopes (OTT) are two 13 m telescopes located at the Onsala Space Observatory in Sweden. With dual linear polarized broad-band (3-14 GHz) receivers, they are part of the next generation Very Long Baseline Interferometry (VLBI) Global Observing System (VGOS) for geodesy and astrometry. In addition to purely geodetic data products, VGOS will regularly produce full-polarisation images of hundreds of radio sources. These rich monitoring data will be valuable for both astronomy and geodesy. In this pilot study we aim to monitor 10 bright radio sources to search for flares or similar activity, and to verify the instrument calibration on long (months) and short (hours) time scales. Method: We observed and analysed 91 short (<30 min) sessions spanning 7 months. We monitored seven potentially variable radio sources (0059+581, 0552+398, 1144+402, 1156+295, 1617+229, 3C418, OJ287) and three reference calibrators (3C147, 3C286, 3C295). We used the Common Astronomy Software Applications (CASA) package to fringe-fit, bandpass-correct and scale the data to obtain flux densities in the four standard VGOS bands: 3.0-3.5 GHz (band 1), 5.2-5.7 GHz (band 2), 6.3-6.8 GHz (band 3), and 10.2-10.7 GHz (band 4). Results: We obtain simultaneous multi-frequency light curves for ten radio sources. A bright multi-frequency flare is observed in the radio source 0059+581. OJ287 and 1156+295 show significant long-term variability. Conclusions: After correcting for instrumental biases, we determine the empirical flux density uncertainty as $\sim$5%. Future refined analysis methods will allow regular monitoring of more and fainter sources. △ Less

Submitted 22 March, 2022; v1 submitted 16 March, 2022; originally announced March 2022.

Comments: 19 pages, 5 figures. Submitted to Experimental Astronomy

A New Moment-Based General-Relativistic Neutrino-Radiation Transport Code: Methods and First Applications to Neutron Star Mergers

Authors: David Radice, Sebastiano Bernuzzi, Albino Perego, Roland Haas

Abstract: We present a new moment-based energy-integrated neutrino transport code for neutron star merger simulations in general relativity. In the merger context, ours is the first code to include Doppler effects at all orders in $\upsilon/c$, retaining all nonlinear neutrino-matter coupling terms. The code is validated with a stringent series of tests. We show that the inclusion of full neutrino-matter co… ▽ More We present a new moment-based energy-integrated neutrino transport code for neutron star merger simulations in general relativity. In the merger context, ours is the first code to include Doppler effects at all orders in $\upsilon/c$, retaining all nonlinear neutrino-matter coupling terms. The code is validated with a stringent series of tests. We show that the inclusion of full neutrino-matter coupling terms is necessary to correctly capture the trapping of neutrinos in relativistically moving media, such as in differentially rotating merger remnants. We perform preliminary simulations proving the robustness of the scheme in simulating ab-initio mergers to black hole collapse and long-term neutron star remnants up to ${\sim}70\,$ms. The latter is the longest dynamical spacetime, 3D, general relativistic simulations with full neutrino transport to date. We compare results obtained at different resolutions and using two different closures for the moment scheme. We do not find evidences of significant out-of-thermodynamic equilibrium effects, such as bulk viscosity, on the postmerger dynamics or gravitational wave emission. Neutrino luminosities and average energies are in good agreement with theory expectations and previous simulations by other groups using similar schemes. We compare dynamical and early wind ejecta properties obtained with M1 and with our older neutrino treatment. We find that the M1 results have systematically larger proton fractions. However, the differences in the nucleosynthesis yields are modest. This work sets the basis for future detailed studies spanning a wider set of neutrino reactions, binaries and equations of state. △ Less

Submitted 9 May, 2023; v1 submitted 29 November, 2021; originally announced November 2021.

Comments: 24 pages, 27 figures. Correct typos in time integrator and spatial limiter equations

Journal ref: MNRAS 512:1499-1521 (2022)

Observing UT1-UTC with VGOS

Authors: Rüdiger Haas, Eskil Varenius, Saho Matsumoto, Matthias Schartner

Abstract: We present first results of UT1-UTC determinations using the VLBI Global Observing System (VGOS). During December 2019 through February 2020 a series of 1~hour long observing sessions were performed using the VGOS stations at Ishioka in Japan and the Onsala twin telescopes in Sweden. The data of this VGOS-B series were correlated, post-correlation processed, and analysed at the Onsala Space Observ… ▽ More We present first results of UT1-UTC determinations using the VLBI Global Observing System (VGOS). During December 2019 through February 2020 a series of 1~hour long observing sessions were performed using the VGOS stations at Ishioka in Japan and the Onsala twin telescopes in Sweden. The data of this VGOS-B series were correlated, post-correlation processed, and analysed at the Onsala Space Observatory. The derived UT1-UTC results were compared to corresponding results from standard legacy S/X Intensive sessions (INT1/INT2), as well to the final values of the International Earth Rotation and Reference Frame Service (IERS), provided in IERS Bulletin~B. The VGOS-B series achieve 3-4 times lower formal uncertainties for the UT1-UTC results than standard legacy S/X INT series. Furthermore, the root mean square (RMS) agreement with respect to the IERS Bulletin~B is 30-40 % better for the VGOS-B results than for the INT1/INT2 results. △ Less

Submitted 30 November, 2020; originally announced November 2020.

Comments: 12 pages, 3 figures, 2 tables

Short-baseline interferometry local-tie experiments at the Onsala Space Observatory

Authors: Eskil Varenius, Rüdiger Haas, Tobias Nilsson

Abstract: We present results from observation, correlation and analysis of interferometric measurements between the three geodetic very long baseline interferometry (VLBI) stations at the Onsala Space Observatory. In total 25 sessions were observed in 2019 and 2020, most of them 24 hours long, all using X-band only. These involved the legacy VLBI station ONSALA60 and the Onsala twin telescopes, ONSA13NE and… ▽ More We present results from observation, correlation and analysis of interferometric measurements between the three geodetic very long baseline interferometry (VLBI) stations at the Onsala Space Observatory. In total 25 sessions were observed in 2019 and 2020, most of them 24 hours long, all using X-band only. These involved the legacy VLBI station ONSALA60 and the Onsala twin telescopes, ONSA13NE and ONSA13SW, two broadband stations for the next generation geodetic VLBI global observing system (VGOS). We used two analysis packages: nuSolve to pre-process the data and solve ambiguities, and ASCOT to solve for station positions, including modelling gravitational deformation of the radio telescopes and other significant effects. We obtained weighted root mean square postfit residuals for each session on the order of 10-15 ps using group delays and 2-5 ps using phase delays. The best performance was achieved on the (rather short) baseline between the VGOS stations. As the main result of this work we determined the coordinates of the Onsala twin telescopes in VTRF2020b with sub-millimeter precision. This new set of coordinates should be used from now on for scheduling, correlation, as a~priori for data analyses, and for comparison with classical local-tie techniques. Finally, we find that positions estimated from phase-delays are offset $\sim+3$ mm in the Up-component with respect to group-delays. Additional modelling of (elevation-dependent) effects may contribute to future understanding of this offset. △ Less

Submitted 12 March, 2021; v1 submitted 30 October, 2020; originally announced October 2020.

Comments: Submitted to Journal of Geodesy

The Occurrence of Rocky Habitable Zone Planets Around Solar-Like Stars from Kepler Data

Authors: Steve Bryson, Michelle Kunimoto, Ravi K. Kopparapu, Jeffrey L. Coughlin, William J. Borucki, David Koch, Victor Silva Aguirre, Christopher Allen, Geert Barentsen, Natalie. M. Batalha, Travis Berger, Alan Boss, Lars A. Buchhave, Christopher J. Burke, Douglas A. Caldwell, Jennifer R. Campbell, Joseph Catanzarite, Hema Chandrasekharan, William J. Chaplin, Jessie L. Christiansen, Jorgen Christensen-Dalsgaard, David R. Ciardi, Bruce D. Clarke, William D. Cochran, Jessie L. Dotson , et al. (57 additional authors not shown)

Abstract: We present occurrence rates for rocky planets in the habitable zones (HZ) of main-sequence dwarf stars based on the Kepler DR25 planet candidate catalog and Gaia-based stellar properties. We provide the first analysis in terms of star-dependent instellation flux, which allows us to track HZ planets. We define $η_\oplus$ as the HZ occurrence of planets with radius between 0.5 and 1.5 $R_\oplus$ orb… ▽ More We present occurrence rates for rocky planets in the habitable zones (HZ) of main-sequence dwarf stars based on the Kepler DR25 planet candidate catalog and Gaia-based stellar properties. We provide the first analysis in terms of star-dependent instellation flux, which allows us to track HZ planets. We define $η_\oplus$ as the HZ occurrence of planets with radius between 0.5 and 1.5 $R_\oplus$ orbiting stars with effective temperatures between 4800 K and 6300 K. We find that $η_\oplus$ for the conservative HZ is between $0.37^{+0.48}_{-0.21}$ (errors reflect 68\% credible intervals) and $0.60^{+0.90}_{-0.36}$ planets per star, while the optimistic HZ occurrence is between $0.58^{+0.73}_{-0.33}$ and $0.88^{+1.28}_{-0.51}$ planets per star. These bounds reflect two extreme assumptions about the extrapolation of completeness beyond orbital periods where DR25 completeness data are available. The large uncertainties are due to the small number of detected small HZ planets. We find similar occurrence rates using both a Poisson likelihood Bayesian analysis and Approximate Bayesian Computation. Our results are corrected for catalog completeness and reliability. Both completeness and the planet occurrence rate are dependent on stellar effective temperature. We also present occurrence rates for various stellar populations and planet size ranges. We estimate with $95\%$ confidence that, on average, the nearest HZ planet around G and K dwarfs is about 6 pc away, and there are about 4 HZ rocky planets around G and K dwarfs within 10 pc of the Sun. △ Less

Submitted 3 November, 2020; v1 submitted 28 October, 2020; originally announced October 2020.

Comments: To appear in The Astronomical Journal

Observation of eccentric binary black hole mergers with second and third generation gravitational wave detector networks

Authors: Zhuo Chen, E. A. Huerta, Joseph Adamo, Roland Haas, Eamonn O'Shea, Prayush Kumar, Chris Moore

Abstract: [Abridged] We introduce an improved version of the Eccentric, Non-spinning, Inspiral-Gaussian-process Merger Approximant (ENIGMA) waveform model. We find that this ready-to-use model can: (i) produce physically consistent signals when sampling over 1M samples chosen over the $m_{\{1,\,2\}}\in[5M_{\odot},\,50M_{\odot}]$ parameter space, and the entire range of binary inclination angles; (ii) produc… ▽ More [Abridged] We introduce an improved version of the Eccentric, Non-spinning, Inspiral-Gaussian-process Merger Approximant (ENIGMA) waveform model. We find that this ready-to-use model can: (i) produce physically consistent signals when sampling over 1M samples chosen over the $m_{\{1,\,2\}}\in[5M_{\odot},\,50M_{\odot}]$ parameter space, and the entire range of binary inclination angles; (ii) produce waveforms within 0.04 seconds from an initial gravitational wave frequency $f_{\textrm{GW}} =15\,\textrm{Hz}$ and at a sample rate of 8192 Hz; and (iii) reproduce the physics of quasi-circular mergers. We utilize ENIGMA to compute the expected signal-to-noise ratio (SNR) distributions of eccentric binary black hole mergers assuming the existence of second and third generation gravitational wave detector networks that include the twin LIGO detectors, Virgo, KAGRA, LIGO-India, a LIGO-type detector in Australia, Cosmic Explorer, and the Einstein Telescope. In the context of advanced LIGO-type detectors, we find that the SNR of eccentric mergers is always larger than quasi-circular mergers for systems with $e_0\leq0.4$ at $f_{\textrm{GW}} =10\,\textrm{Hz}$, even if the timespan of eccentric signals is just a third of quasi-circular systems with identical total mass and mass-ratio. For Cosmic Explorer-type detector networks, we find that eccentric mergers have similar SNRs than quasi-circular systems for $e_0\leq0.3$ at $f_{\textrm{GW}} =10\,\textrm{Hz}$. Systems with $e_0\sim0.5$ at $f_{\textrm{GW}} =10\,\textrm{Hz}$ have SNRs that range between 50%-90% of the SNR produced by quasi-circular mergers, even if these eccentric signals are just between a third to a tenth the length of quasi-circular systems. For Einstein Telescope-type detectors, we find that eccentric mergers have similar SNRs than quasi-circular systems for $e_0\leq0.4$ at $f_{\textrm{GW}} =5\,\textrm{Hz}$. △ Less

Submitted 13 April, 2021; v1 submitted 7 August, 2020; originally announced August 2020.

Comments: 25 pages, 19 figures, 1 appendix. Accepted to Phys. Rev. D

ACM Class: J.2

Journal ref: Phys. Rev. D 103, 084018 (2021)

A magnetar engine for short GRBs and kilonovae

Authors: Philipp Mösta, David Radice, Roland Haas, Erik Schnetter, Sebastiano Bernuzzi

Abstract: We investigate the influence of magnetic fields on the evolution of binary neutron-star (BNS) merger remnants via three-dimensional (3D) dynamical-spacetime general-relativistic (GR) magnetohydrodynamic (MHD) simulations. We evolve a postmerger remnant with an initial poloidal magnetic field, resolve the magnetoturbulence driven by shear flows, and include a microphysical finite-temperature equati… ▽ More We investigate the influence of magnetic fields on the evolution of binary neutron-star (BNS) merger remnants via three-dimensional (3D) dynamical-spacetime general-relativistic (GR) magnetohydrodynamic (MHD) simulations. We evolve a postmerger remnant with an initial poloidal magnetic field, resolve the magnetoturbulence driven by shear flows, and include a microphysical finite-temperature equation of state (EOS). A neutrino leakage scheme that captures the overall energetics and lepton number exchange is also included. We find that turbulence induced by the magnetorotational instability (MRI) in the hypermassive neutron star (HMNS) amplifies magnetic field to beyond magnetar-strength ($10^{15}\, \mathrm{G}$). The ultra-strong toroidal field is able to launch a relativistic jet from the HMNS. We also find a magnetized wind that ejects neutron-rich material with a rate of $\dot{M}_{\mathrm{ej}} \simeq 1 \times10^{-1}\, \mathrm{M_{\odot}\, s^{-1}}$. The total ejecta mass in our simulation is $5\times 10^{-3}\, \mathrm{M_{\odot}}$. This makes the ejecta from the HMNS an important component in BNS mergers and a promising source of $r$-process elements that can power a kilonova. The jet from the HMNS reaches a terminal Lorentz factor of $\sim 5$ in our highest-resolution simulation. The formation of this jet is aided by neutrino-cooling preventing the accretion disk from protruding into the polar region. As neutrino pair-annihilation and radiative processes in the jet (which were not included in the simulations) will boost the Lorentz factor in the jet further, our simulations demonstrate that magnetars formed in BNS mergers are a viable engine for short gamma-ray bursts (sGRBs). △ Less

Submitted 28 October, 2020; v1 submitted 12 March, 2020; originally announced March 2020.

Comments: ApJL published version

Journal ref: ApJL, 901, L37 (2020)

Enabling real-time multi-messenger astrophysics discoveries with deep learning

Authors: E. A. Huerta, Gabrielle Allen, Igor Andreoni, Javier M. Antelis, Etienne Bachelet, Bruce Berriman, Federica Bianco, Rahul Biswas, Matias Carrasco, Kyle Chard, Minsik Cho, Philip S. Cowperthwaite, Zachariah B. Etienne, Maya Fishbach, Francisco Förster, Daniel George, Tom Gibbs, Matthew Graham, William Gropp, Robert Gruendl, Anushri Gupta, Roland Haas, Sarah Habib, Elise Jennings, Margaret W. G. Johnson , et al. (35 additional authors not shown)

Abstract: Multi-messenger astrophysics is a fast-growing, interdisciplinary field that combines data, which vary in volume and speed of data processing, from many different instruments that probe the Universe using different cosmic messengers: electromagnetic waves, cosmic rays, gravitational waves and neutrinos. In this Expert Recommendation, we review the key challenges of real-time observations of gravit… ▽ More Multi-messenger astrophysics is a fast-growing, interdisciplinary field that combines data, which vary in volume and speed of data processing, from many different instruments that probe the Universe using different cosmic messengers: electromagnetic waves, cosmic rays, gravitational waves and neutrinos. In this Expert Recommendation, we review the key challenges of real-time observations of gravitational wave sources and their electromagnetic and astroparticle counterparts, and make a number of recommendations to maximize their potential for scientific discovery. These recommendations refer to the design of scalable and computationally efficient machine learning algorithms; the cyber-infrastructure to numerically simulate astrophysical sources, and to process and interpret multi-messenger astrophysics data; the management of gravitational wave detections to trigger real-time alerts for electromagnetic and astroparticle follow-ups; a vision to harness future developments of machine learning and cyber-infrastructure resources to cope with the big-data requirements; and the need to build a community of experts to realize the goals of multi-messenger astrophysics. △ Less

Submitted 26 November, 2019; originally announced November 2019.

Comments: Invited Expert Recommendation for Nature Reviews Physics. The art work produced by E. A. Huerta and Shawn Rosofsky for this article was used by Carl Conway to design the cover of the October 2019 issue of Nature Reviews Physics

Journal ref: Nature Reviews Physics volume 1, pages 600-608 (2019)

Deep Learning for Multi-Messenger Astrophysics: A Gateway for Discovery in the Big Data Era

Authors: Gabrielle Allen, Igor Andreoni, Etienne Bachelet, G. Bruce Berriman, Federica B. Bianco, Rahul Biswas, Matias Carrasco Kind, Kyle Chard, Minsik Cho, Philip S. Cowperthwaite, Zachariah B. Etienne, Daniel George, Tom Gibbs, Matthew Graham, William Gropp, Anushri Gupta, Roland Haas, E. A. Huerta, Elise Jennings, Daniel S. Katz, Asad Khan, Volodymyr Kindratenko, William T. C. Kramer, Xin Liu, Ashish Mahabal , et al. (23 additional authors not shown)

Abstract: This report provides an overview of recent work that harnesses the Big Data Revolution and Large Scale Computing to address grand computational challenges in Multi-Messenger Astrophysics, with a particular emphasis on real-time discovery campaigns. Acknowledging the transdisciplinary nature of Multi-Messenger Astrophysics, this document has been prepared by members of the physics, astronomy, compu… ▽ More This report provides an overview of recent work that harnesses the Big Data Revolution and Large Scale Computing to address grand computational challenges in Multi-Messenger Astrophysics, with a particular emphasis on real-time discovery campaigns. Acknowledging the transdisciplinary nature of Multi-Messenger Astrophysics, this document has been prepared by members of the physics, astronomy, computer science, data science, software and cyberinfrastructure communities who attended the NSF-, DOE- and NVIDIA-funded "Deep Learning for Multi-Messenger Astrophysics: Real-time Discovery at Scale" workshop, hosted at the National Center for Supercomputing Applications, October 17-19, 2018. Highlights of this report include unanimous agreement that it is critical to accelerate the development and deployment of novel, signal-processing algorithms that use the synergy between artificial intelligence (AI) and high performance computing to maximize the potential for scientific discovery with Multi-Messenger Astrophysics. We discuss key aspects to realize this endeavor, namely (i) the design and exploitation of scalable and computationally efficient AI algorithms for Multi-Messenger Astrophysics; (ii) cyberinfrastructure requirements to numerically simulate astrophysical sources, and to process and interpret Multi-Messenger Astrophysics data; (iii) management of gravitational wave detections and triggers to enable electromagnetic and astro-particle follow-ups; (iv) a vision to harness future developments of machine and deep learning and cyberinfrastructure resources to cope with the scale of discovery in the Big Data Era; (v) and the need to build a community that brings domain experts together with data scientists on equal footing to maximize and accelerate discovery in the nascent field of Multi-Messenger Astrophysics. △ Less

Submitted 1 February, 2019; originally announced February 2019.

Comments: 15 pages, no figures. White paper based on the "Deep Learning for Multi-Messenger Astrophysics: Real-time Discovery at Scale" workshop, hosted at NCSA, October 17-19, 2018 http://www.ncsa.illinois.edu/Conferences/DeepLearningLSST/

Physics of eccentric binary black hole mergers: A numerical relativity perspective

Authors: E. A. Huerta, Roland Haas, Sarah Habib, Anushri Gupta, Adam Rebei, Vishnu Chavva, Daniel Johnson, Shawn Rosofsky, Erik Wessel, Bhanu Agarwal, Diyu Luo, Wei Ren

Abstract: Gravitational wave observations of eccentric binary black hole mergers will provide unequivocal evidence for the formation of these systems through dynamical assembly in dense stellar environments. The study of these astrophysically motivated sources is timely in view of electromagnetic observations, consistent with the existence of stellar mass black holes in the globular cluster M22 and in the G… ▽ More Gravitational wave observations of eccentric binary black hole mergers will provide unequivocal evidence for the formation of these systems through dynamical assembly in dense stellar environments. The study of these astrophysically motivated sources is timely in view of electromagnetic observations, consistent with the existence of stellar mass black holes in the globular cluster M22 and in the Galactic center, and the proven detection capabilities of ground-based gravitational wave detectors. In order to get insights into the physics of these objects in the dynamical, strong-field gravity regime, we present a catalog of 89 numerical relativity waveforms that describe binary systems of non-spinning black holes with mass-ratios $1\leq q \leq 10$, and initial eccentricities as high as $e_0=0.18$ fifteen cycles before merger. We use this catalog to quantify the loss of energy and angular momentum through gravitational radiation, and the astrophysical properties of the black hole remnant, including its final mass and spin, and recoil velocity. We discuss the implications of these results for gravitational wave source modeling, and the design of algorithms to search for and identify eccentric binary black hole mergers in realistic detection scenarios. △ Less

Submitted 5 September, 2019; v1 submitted 21 January, 2019; originally announced January 2019.

Comments: 11 pages, 5 figures, 2 appendices. A visualization of this numerical relativity waveform catalog is available at https://gravity.ncsa.illinois.edu/products/outreach/; v2: 13 pages, 5 figures, calculations for angular momentum emission and recoil velocities are now included, references added. Accepted to Phys. Rev. D

ACM Class: J.2

Journal ref: Phys. Rev. D 100, 064003 (2019)

Gravitational waveforms from SpEC simulations : neutron star-neutron star and low-mass black hole-neutron star binaries

Authors: Francois Foucart, Matthew D. Duez, Tanja Hinderer, Jesus Caro, Andrew R. Williamson, Michael Boyle, Alessandra Buonanno, Roland Haas, Daniel A. Hemberger, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheel

Abstract: Gravitational waveforms from numerical simulations are a critical tool to test and analytically calibrate the waveform models used to study the properties of merging compact objects. In this paper, we present a series of high-accuracy waveforms produced with the SpEC code for systems involving at least one neutron star. We provide for the first time waveforms with sub-radian accuracy over more tha… ▽ More Gravitational waveforms from numerical simulations are a critical tool to test and analytically calibrate the waveform models used to study the properties of merging compact objects. In this paper, we present a series of high-accuracy waveforms produced with the SpEC code for systems involving at least one neutron star. We provide for the first time waveforms with sub-radian accuracy over more than twenty cycles for low-mass black hole-neutron star binaries, including binaries with non-spinning objects, and binaries with rapidly spinning neutron stars that maximize the impact on the gravitational wave signal of the near-resonant growth of the fundamental excitation mode of the neutron star (f-mode). We also provide for the first time with SpEC a high-accuracy neutron star-neutron star waveform. These waveforms are made publicly available as part of the SxS catalogue. We compare our results to analytical waveform models currently implemented in data analysis pipelines. For most simulations, the models lie outside of the predicted numerical errors in the last few orbits before merger, but do not show systematic deviations from the numerical results: comparing different models appears to provide reasonable estimates of the modeling errors. The sole exception is the equal-mass simulation using a rapidly counter-rotating neutron star to maximize the impact of the excitation of the f-mode, for which all models perform poorly. This is however expected, as even the single model that takes f-mode excitation into account ignores the significant impact of the neutron star spin on the f-mode excitation frequency. △ Less

Submitted 17 December, 2018; originally announced December 2018.

Comments: 17p, 10 figures

Journal ref: Phys. Rev. D 99, 044008 (2019)

Supporting High-Performance and High-Throughput Computing for Experimental Science

Authors: E. A. Huerta, Roland Haas, Shantenu Jha, Mark Neubauer, Daniel S. Katz

Abstract: The advent of experimental science facilities-instruments and observatories, such as the Large Hadron Collider, the Laser Interferometer Gravitational Wave Observatory, and the upcoming Large Synoptic Survey Telescope-has brought about challenging, large-scale computational and data processing requirements. Traditionally, the computing infrastructure to support these facility's requirements were o… ▽ More The advent of experimental science facilities-instruments and observatories, such as the Large Hadron Collider, the Laser Interferometer Gravitational Wave Observatory, and the upcoming Large Synoptic Survey Telescope-has brought about challenging, large-scale computational and data processing requirements. Traditionally, the computing infrastructure to support these facility's requirements were organized into separate infrastructure that supported their high-throughput needs and those that supported their high-performance computing needs. We argue that to enable and accelerate scientific discovery at the scale and sophistication that is now needed, this separation between high-performance computing and high-throughput computing must be bridged and an integrated, unified infrastructure provided. In this paper, we discuss several case studies where such infrastructure has been implemented. These case studies span different science domains, software systems, and application requirements as well as levels of sustainability. A further aim of this paper is to provide a basis to determine the common characteristics and requirements of such infrastructure, as well as to begin a discussion of how best to support the computing requirements of existing and future experimental science facilities. △ Less

Submitted 8 February, 2019; v1 submitted 6 October, 2018; originally announced October 2018.

Comments: 13 pages, 7 figures. Accepted to Computing and Software for Big Science

MSC Class: 90C06; 68Q85

Journal ref: Comput Softw Big Sci (2019) 3: 5

Fusing numerical relativity and deep learning to detect higher-order multipole waveforms from eccentric binary black hole mergers

Authors: Adam Rebei, E. A. Huerta, Sibo Wang, Sarah Habib, Roland Haas, Daniel Johnson, Daniel George

Abstract: We determine the mass-ratio, eccentricity and binary inclination angles that maximize the contribution of the higher-order waveform multipoles $(\ell, \, |m|)= \{(2,\,2),\, (2,\,1),\, (3,\,3),\, (3,\,2), \, (3,\,1),\, (4,\,4),\, (4,\,3),\, (4,\,2),\,(4,\,1)\}$ for the gravitational wave detection of eccentric binary black hole mergers. We carry out this study using numerical relativity waveforms t… ▽ More We determine the mass-ratio, eccentricity and binary inclination angles that maximize the contribution of the higher-order waveform multipoles $(\ell, \, |m|)= \{(2,\,2),\, (2,\,1),\, (3,\,3),\, (3,\,2), \, (3,\,1),\, (4,\,4),\, (4,\,3),\, (4,\,2),\,(4,\,1)\}$ for the gravitational wave detection of eccentric binary black hole mergers. We carry out this study using numerical relativity waveforms that describe non-spinning black hole binaries with mass-ratios $1\leq q \leq 10$, and orbital eccentricities as high as $e_0=0.18$ fifteen cycles before merger. For stellar-mass, asymmetric mass-ratio, binary black hole mergers, and assuming LIGO's Zero Detuned High Power configuration, we find that in regions of parameter space where black hole mergers modeled with $\ell=|m|=2$ waveforms have vanishing signal-to-noise ratios, the inclusion of $(\ell, \, |m|)$ modes enables the observation of these sources with signal-to-noise ratios that range between 30\% to 45\% the signal-to-noise ratio of optimally oriented binary black hole mergers modeled with $\ell=|m|=2$ numerical relativity waveforms. Having determined the parameter space where $(\ell, \, |m|)$ modes are important for gravitational wave detection, we construct waveform signals that describe these astrophysically motivate scenarios, and demonstrate that these topologically complex signals can be detected and characterized in real LIGO noise with deep learning algorithms. △ Less

Submitted 15 August, 2019; v1 submitted 25 July, 2018; originally announced July 2018.

Comments: 17 pages, 11 figures, 1 appendix, accepted to Phys. Rev. D

ACM Class: J.2

Journal ref: Phys. Rev. D 100, 044025 (2019)

Kepler Data Validation I -- Architecture, Diagnostic Tests, and Data Products for Vetting Transiting Planet Candidates

Authors: Joseph D. Twicken, Joseph H. Catanzarite, Bruce D. Clarke, Forrest Girouard, Jon M. Jenkins, Todd C. Klaus, Jie Li, Sean D. McCauliff, Shawn E. Seader, Peter Tenenbaum, Bill Wohler, Stephen T. Bryson, Christopher J. Burke, Douglas A. Caldwell, Michael R. Haas, Christopher E. Henze, Dwight T. Sanderfer

Abstract: The Kepler Mission was designed to identify and characterize transiting planets in the Kepler Field of View and to determine their occurrence rates. Emphasis was placed on identification of Earth-size planets orbiting in the Habitable Zone of their host stars. Science data were acquired for a period of four years. Long-cadence data with 29.4 min sampling were obtained for ~200,000 individual stell… ▽ More The Kepler Mission was designed to identify and characterize transiting planets in the Kepler Field of View and to determine their occurrence rates. Emphasis was placed on identification of Earth-size planets orbiting in the Habitable Zone of their host stars. Science data were acquired for a period of four years. Long-cadence data with 29.4 min sampling were obtained for ~200,000 individual stellar targets in at least one observing quarter in the primary Kepler Mission. Light curves for target stars are extracted in the Kepler Science Data Processing Pipeline, and are searched for transiting planet signatures. A Threshold Crossing Event is generated in the transit search for targets where the transit detection threshold is exceeded and transit consistency checks are satisfied. These targets are subjected to further scrutiny in the Data Validation (DV) component of the Pipeline. Transiting planet candidates are characterized in DV, and light curves are searched for additional planets after transit signatures are modeled and removed. A suite of diagnostic tests is performed on all candidates to aid in discrimination between genuine transiting planets and instrumental or astrophysical false positives. Data products are generated per target and planet candidate to document and display transiting planet model fit and diagnostic test results. These products are exported to the Exoplanet Archive at the NASA Exoplanet Science Institute, and are available to the community. We describe the DV architecture and diagnostic tests, and provide a brief overview of the data products. Transiting planet modeling and the search for multiple planets on individual targets are described in a companion paper. The final revision of the Kepler Pipeline code base is available to the general public through GitHub. The Kepler Pipeline has also been modified to support the TESS Mission which will commence in 2018. △ Less

Submitted 24 April, 2018; v1 submitted 12 March, 2018; originally announced March 2018.

Comments: 84 pages, 26 figures, 2 tables. Accepted for publication in PASP on 9 March 2018. Published in PASP on 24 April 2018

The NANOGrav 11-year Data Set: Pulsar-timing Constraints On The Stochastic Gravitational-wave Background

Authors: Z. Arzoumanian, P. T. Baker, A. Brazier, S. Burke-Spolaor, S. J. Chamberlin, S. Chatterjee, B. Christy, J. M. Cordes, N. J. Cornish, F. Crawford, H. Thankful Cromartie, K. Crowter, M. DeCesar, P. B. Demorest, T. Dolch, J. A. Ellis, R. D. Ferdman, E. Ferrara, W. M. Folkner, E. Fonseca, N. Garver-Daniels, P. A. Gentile, R. Haas, J. S. Hazboun, E. A. Huerta , et al. (35 additional authors not shown)

Abstract: We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released $11$-year dataset from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no significant evidence for a GWB, we place constraints on a GWB from a population of supermassive black-hole binaries, cosmic strings, and a primordial GWB. For the first time, we find that… ▽ More We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released $11$-year dataset from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no significant evidence for a GWB, we place constraints on a GWB from a population of supermassive black-hole binaries, cosmic strings, and a primordial GWB. For the first time, we find that the GWB upper limits and detection statistics are sensitive to the Solar System ephemeris (SSE) model used, and that SSE errors can mimic a GWB signal. We developed an approach that bridges systematic SSE differences, producing the first PTA constraints that are robust against SSE uncertainties. We thus place a $95\%$ upper limit on the GW strain amplitude of $A_\mathrm{GWB}<1.45\times 10^{-15}$ at a frequency of $f=1$ yr$^{-1}$ for a fiducial $f^{-2/3}$ power-law spectrum, and with inter-pulsar correlations modeled. This is a factor of $\sim 2$ improvement over the NANOGrav $9$-year limit, calculated using the same procedure. Previous PTA upper limits on the GWB will need revision in light of SSE systematic uncertainties. We use our constraints to characterize the combined influence on the GWB of the stellar mass-density in galactic cores, the eccentricity of SMBH binaries, and SMBH--galactic-bulge scaling relationships. We constrain cosmic-string tension using recent simulations, yielding an SSE-marginalized $95\%$ upper limit on the cosmic string tension of $Gμ< 5.3\times 10^{-11}$---a factor of $\sim 2$ better than the published NANOGrav $9$-year constraints. Our SSE-marginalized $95\%$ upper limit on the energy density of a primordial GWB (for a radiation-dominated post-inflation Universe) is $Ω_\mathrm{GWB}(f)h^2<3.4\times10^{-10}$. △ Less

Submitted 7 June, 2018; v1 submitted 8 January, 2018; originally announced January 2018.

Comments: 21 pages, 12 figures, 9 tables. Published in The Astrophysical Journal. Please send any comments/questions to S. R. Taylor ([email protected])

Journal ref: The Astrophysical Journal, Volume 859, Number 1, 2018

R-process Nucleosynthesis from Three-Dimensional Magnetorotational Core-Collapse Supernovae

Authors: Philipp Mösta, Luke F. Roberts, Goni Halevi, Christian D. Ott, Jonas Lippuner, Roland Haas, Erik Schnetter

Abstract: We investigate r-process nucleosynthesis in three-dimensional (3D) general-relativistic magnetohydrodynamic simulations of rapidly rotating strongly magnetized core collapse. The simulations include a microphysical finite-temperature equation of state and a leakage scheme that captures the overall energetics and lepton number exchange due to postbounce neutrino emission and absorption. We track th… ▽ More We investigate r-process nucleosynthesis in three-dimensional (3D) general-relativistic magnetohydrodynamic simulations of rapidly rotating strongly magnetized core collapse. The simulations include a microphysical finite-temperature equation of state and a leakage scheme that captures the overall energetics and lepton number exchange due to postbounce neutrino emission and absorption. We track the composition of the ejected material using the nuclear reaction network SkyNet. Our results show that the 3D dynamics of magnetorotational core-collapse supernovae (CCSN) are important for their nucleosynthetic signature. We find that production of r-process material beyond the second peak is reduced by a factor of 100 when the magnetorotational jets produced by the rapidly rotating core undergo a kink instability. Our results indicate that 3D magnetorotationally powered CCSNe are a robust r-process source only if they are obtained by the collapse of cores with unrealistically large precollapse magnetic fields of order $10^{13}$G. Additionally, a comparison simulation that we restrict to axisymmetry, results in overly optimistic r-process production for lower magnetic field strengths. △ Less

Submitted 26 December, 2017; originally announced December 2017.

Comments: 10 pages, 9 figures, 2 tables. submitted to ApJ

Report number: YITP-17-129; LANL LA-UR-17-31278

The Progenitor Dependence of Three-Dimensional Core-Collapse Supernovae

Authors: C. D. Ott, L. F. Roberts, A. da Silva Schneider, J. M. Fedrow, R. Haas, E. Schnetter

Abstract: We present a first study of the progenitor star dependence of the three-dimensional (3D) neutrino mechanism of core-collapse supernovae. We employ full 3D general-relativistic multi-group neutrino radiation-hydrodynamics and simulate the post-bounce evolutions of progenitors with zero-age main sequence masses of $12$, $15$, $20$, $27$, and $40\,M_\odot$. All progenitors, with the exception of the… ▽ More We present a first study of the progenitor star dependence of the three-dimensional (3D) neutrino mechanism of core-collapse supernovae. We employ full 3D general-relativistic multi-group neutrino radiation-hydrodynamics and simulate the post-bounce evolutions of progenitors with zero-age main sequence masses of $12$, $15$, $20$, $27$, and $40\,M_\odot$. All progenitors, with the exception of the $12\,M_\odot$ star, experience shock runaway by the end of their simulations. In most cases, a strongly asymmetric explosion will result. We find three qualitatively distinct evolutions that suggest a complex dependence of explosion dynamics on progenitor density structure, neutrino heating, and 3D flow. (1) Progenitors with massive cores, shallow density profiles, and high post-core-bounce accretion rates experience very strong neutrino heating and neutrino-driven turbulent convection, leading to early shock runaway. Accretion continues at a high rate, likely leading to black hole formation. (2) Intermediate progenitors experience neutrino-driven, turbulence-aided explosions triggered by the arrival of density discontinuities at the shock. These occur typically at the silicon/silicon-oxygen shell boundary. (3) Progenitors with small cores and density profiles without strong discontinuities experience shock recession and develop the 3D standing-accretion shock instability (SASI). Shock runaway ensues late, once declining accretion rate, SASI, and neutrino-driven convection create favorable conditions. These differences in explosion times and dynamics result in a non-monotonic relationship between progenitor and compact remnant mass. △ Less

Submitted 4 December, 2017; originally announced December 2017.

Comments: 8 pages, 6 figures, submitted

Report number: YITP-17-122

Eccentric, nonspinning, inspiral, Gaussian-process merger approximant for the detection and characterization of eccentric binary black hole mergers

Authors: E. A. Huerta, C. J. Moore, Prayush Kumar, Daniel George, Alvin J. K. Chua, Roland Haas, Erik Wessel, Daniel Johnson, Derek Glennon, Adam Rebei, A. Miguel Holgado, Jonathan R. Gair, Harald P. Pfeiffer

Abstract: We present $\texttt{ENIGMA}$, a time domain, inspiral-merger-ringdown waveform model that describes non-spinning binary black holes systems that evolve on moderately eccentric orbits. The inspiral evolution is described using a consistent combination of post-Newtonian theory, self-force and black hole perturbation theory. Assuming eccentric binaries that circularize prior to coalescence, we smooth… ▽ More We present $\texttt{ENIGMA}$, a time domain, inspiral-merger-ringdown waveform model that describes non-spinning binary black holes systems that evolve on moderately eccentric orbits. The inspiral evolution is described using a consistent combination of post-Newtonian theory, self-force and black hole perturbation theory. Assuming eccentric binaries that circularize prior to coalescence, we smoothly match the eccentric inspiral with a stand-alone, quasi-circular merger, which is constructed using machine learning algorithms that are trained with quasi-circular numerical relativity waveforms. We show that $\texttt{ENIGMA}$ reproduces with excellent accuracy the dynamics of quasi-circular compact binaries. We validate $\texttt{ENIGMA}$ using a set of $\texttt{Einstein Toolkit}$ eccentric numerical relativity waveforms, which describe eccentric binary black hole mergers with mass-ratios between $1 \leq q \leq 5.5$, and eccentricities $e_0 \lesssim 0.2$ ten orbits before merger. We use this model to explore in detail the physics that can be extracted with moderately eccentric, non-spinning binary black hole mergers. We use $\texttt{ENIGMA}$ to show that GW150914, GW151226, GW170104, GW170814 and GW170608 can be effectively recovered with spinning, quasi-circular templates if the eccentricity of these events at a gravitational wave frequency of 10Hz satisfies $e_0\leq \{0.175,\, 0.125,\,0.175,\,0.175,\, 0.125\}$, respectively. We show that if these systems have eccentricities $e_0\sim 0.1$ at a gravitational wave frequency of 10Hz, they can be misclassified as quasi-circular binaries due to parameter space degeneracies between eccentricity and spin corrections. Using our catalog of eccentric numerical relativity simulations, we discuss the importance of including higher-order waveform multipoles in gravitational wave searches of eccentric binary black hole mergers. △ Less

Submitted 24 January, 2018; v1 submitted 16 November, 2017; originally announced November 2017.

Comments: 19 pages, 10 figures, 1 Appendix. v2: we use numerical relativity simulations to quantify the importance of including higher-order waveform multipoles for the detection of eccentric binary black hole mergers, references added. Accepted to Phys. Rev. D

ACM Class: J.2

Journal ref: Phys. Rev. D 97, 024031 (2018)

Probing the gravitational redshift with an Earth-orbiting satellite

Authors: D. A. Litvinov, V. N. Rudenko, A. V. Alakoz, U. Bach, N. Bartel, A. V. Belonenko, K. G. Belousov, M. Bietenholz, A. V. Biriukov, R. Carman, G. Cimó, C. Courde, D. Dirkx, D. A. Duev, A. I. Filetkin, G. Granato, L. I. Gurvits, A. V. Gusev, R. Haas, G. Herold, A. Kahlon, B. Z. Kanevsky, V. L. Kauts, G. D. Kopelyansky, A. V. Kovalenko , et al. (23 additional authors not shown)

Abstract: We present an approach to testing the gravitational redshift effect using the RadioAstron satellite. The experiment is based on a modification of the Gravity Probe A scheme of nonrelativistic Doppler compensation and benefits from the highly eccentric orbit and ultra-stable atomic hydrogen maser frequency standard of the RadioAstron satellite. Using the presented techniques we expect to reach an a… ▽ More We present an approach to testing the gravitational redshift effect using the RadioAstron satellite. The experiment is based on a modification of the Gravity Probe A scheme of nonrelativistic Doppler compensation and benefits from the highly eccentric orbit and ultra-stable atomic hydrogen maser frequency standard of the RadioAstron satellite. Using the presented techniques we expect to reach an accuracy of the gravitational redshift test of order $10^{-5}$, a magnitude better than that of Gravity Probe A. Data processing is ongoing, our preliminary results agree with the validity of the Einstein Equivalence Principle. △ Less

Submitted 27 October, 2017; originally announced October 2017.

Comments: Submitted to Physics Letters A

Evolution of the Magnetized, Neutrino-Cooled Accretion Disk in the Aftermath of a Black Hole Neutron Star Binary Merger

Authors: Fatemeh Hossein Nouri, Matthew D. Duez, Francois Foucart, M. Brett Deaton, Roland Haas, Milad Haddadi, Lawrence E. Kidder, Christian D. Ott, Harald P. Pfeiffer, Mark A. Scheel, Bela Szilagyi

Abstract: Black hole-torus systems from compact binary mergers are possible engines for gamma-ray bursts (GRBs). During the early evolution of the post-merger remnant, the state of the torus is determined by a combination of neutrino cooling and magnetically-driven heating processes, so realistic models must include both effects. In this paper, we study the post-merger evolution of a magnetized black hole-n… ▽ More Black hole-torus systems from compact binary mergers are possible engines for gamma-ray bursts (GRBs). During the early evolution of the post-merger remnant, the state of the torus is determined by a combination of neutrino cooling and magnetically-driven heating processes, so realistic models must include both effects. In this paper, we study the post-merger evolution of a magnetized black hole-neutron star binary system using the Spectral Einstein Code (SpEC) from an initial post-merger state provided by previous numerical relativity simulations. We use a finite-temperature nuclear equation of state and incorporate neutrino effects in a leakage approximation. To achieve the needed accuracy, we introduce improvements to SpEC's implementation of general-relativistic magnetohydrodynamics (MHD), including the use of cubed-sphere multipatch grids and an improved method for dealing with supersonic accretion flows where primitive variable recovery is difficult. We find that a seed magnetic field triggers a sustained source of heating, but its thermal effects are largely cancelled by the accretion and spreading of the torus from MHD-related angular momentum transport. The neutrino luminosity peaks at the start of the simulation, and then drops significantly over the first 20\,ms but in roughly the same way for magnetized and nonmagnetized disks. The heating rate and disk's luminosity decrease much more slowly thereafter. These features of the evolution are insensitive to grid structure and resolution, formulation of the MHD equations, and seed field strength, although turbulent effects are not fully converged △ Less

Submitted 20 October, 2017; originally announced October 2017.

Comments: 17 pages, 18 figures

Journal ref: Phys. Rev. D 97, 083014 (2018)

Planetary Candidates Observed by Kepler. VIII. A Fully Automated Catalog With Measured Completeness and Reliability Based on Data Release 25

Authors: Susan E. Thompson, Jeffrey L. Coughlin, Kelsey Hoffman, Fergal Mullally, Jessie L. Christiansen, Christopher J. Burke, Steve Bryson, Natalie Batalha, Michael R. Haas, Joseph Catanzarite, Jason F. Rowe, Geert Barentsen, Douglas A. Caldwell, Bruce D. Clarke, Jon M. Jenkins, Jie Li, David W. Latham, Jack J. Lissauer, Savita Mathur, Robert L. Morris, Shawn E. Seader, Jeffrey C. Smith, Todd C. Klaus, Joseph D. Twicken, Bill Wohler , et al. (36 additional authors not shown)

Abstract: We present the Kepler Object of Interest (KOI) catalog of transiting exoplanets based on searching four years of Kepler time series photometry (Data Release 25, Q1-Q17). The catalog contains 8054 KOIs of which 4034 are planet candidates with periods between 0.25 and 632 days. Of these candidates, 219 are new and include two in multi-planet systems (KOI-82.06 and KOI-2926.05), and ten high-reliabil… ▽ More We present the Kepler Object of Interest (KOI) catalog of transiting exoplanets based on searching four years of Kepler time series photometry (Data Release 25, Q1-Q17). The catalog contains 8054 KOIs of which 4034 are planet candidates with periods between 0.25 and 632 days. Of these candidates, 219 are new and include two in multi-planet systems (KOI-82.06 and KOI-2926.05), and ten high-reliability, terrestrial-size, habitable zone candidates. This catalog was created using a tool called the Robovetter which automatically vets the DR25 Threshold Crossing Events (TCEs, Twicken et al. 2016). The Robovetter also vetted simulated data sets and measured how well it was able to separate TCEs caused by noise from those caused by low signal-to-noise transits. We discusses the Robovetter and the metrics it uses to sort TCEs. For orbital periods less than 100 days the Robovetter completeness (the fraction of simulated transits that are determined to be planet candidates) across all observed stars is greater than 85%. For the same period range, the catalog reliability (the fraction of candidates that are not due to instrumental or stellar noise) is greater than 98%. However, for low signal-to-noise candidates between 200 and 500 days around FGK dwarf stars, the Robovetter is 76.7% complete and the catalog is 50.5% reliable. The KOI catalog, the transit fits and all of the simulated data used to characterize this catalog are available at the NASA Exoplanet Archive. △ Less

Submitted 4 March, 2018; v1 submitted 18 October, 2017; originally announced October 2017.

Comments: 61 pages, 23 Figures, 9 Tables, Accepted to The Astrophysical Journal Supplement Series

Gravitational Waves from Binary Black Hole Mergers Inside of Stars

Authors: J. M. Fedrow, C. D. Ott, U. Sperhake, J. Blackman, R. Haas, C. Reisswig, A. De Felice

Abstract: We present results from a controlled numerical experiment investigating the effect of stellar density gas on the coalescence of binary black holes (BBHs) and the resulting gravitational waves (GWs). This investigation is motivated by the proposed stellar core fragmentation scenario for BBH formation and the associated possibility of an electromagnetic counterpart to a BBH GW event. We employ full… ▽ More We present results from a controlled numerical experiment investigating the effect of stellar density gas on the coalescence of binary black holes (BBHs) and the resulting gravitational waves (GWs). This investigation is motivated by the proposed stellar core fragmentation scenario for BBH formation and the associated possibility of an electromagnetic counterpart to a BBH GW event. We employ full numerical relativity coupled with general-relativistic hydrodynamics and set up a $30 + 30 M_\odot$ BBH (motivated by GW150914) inside gas with realistic stellar densities. Our results show that at densities $ρ\gtrsim 10^6 - 10^7 \, \mathrm{g \, cm}^{-3}$ dynamical friction between the BHs and gas changes the coalescence dynamics and the GW signal in an unmistakable way. We show that for GW150914, LIGO observations conclusively rule out BBH coalescence inside stellar gas of $ρ\gtrsim 10^7 \, \mathrm{g\,cm}^{-3}$. Typical densities in the collapsing cores of massive stars are in excess of this density. This excludes the fragmentation scenario for the formation of GW150914. △ Less

Submitted 18 September, 2017; v1 submitted 24 April, 2017; originally announced April 2017.

Comments: 8 pages, 5 figures, Accepted for publication in Phys. Rev. Lett

Report number: YITP-17-40

Time-Series Analysis of Broadband Photometry of Neptune from K2

Authors: Jason F. Rowe, Patrick Gaulme, Jack J. Lissauer, Mark S. Marley, Amy A. Simon, Heidi B. Hammel, Victor Silva Aguirre, Thomas Barclay, Othman Benomar, Patrick Boumier, Douglas A. Caldwell, Sarah L. Casewell, William J. Chaplin, Knicole D. Colon, Enrico Corsaro, G. R. Davies, Jonathan J. Fortney, Rafael A. Garcia, John E. Gizis, Michael R. Haas, Benoit Mosser, Francois-Xavier Schmider

Abstract: We report here on our search for excess power in photometry of Neptune collected by the K2 mission that may be due to intrinsic global oscillations of the planet Neptune. To conduct this search, we developed new methods to correct for instrumental effects such as intrapixel variability and gain variations. We then extracted and analyzed the time-series photometry of Neptune from 49 days of nearly… ▽ More We report here on our search for excess power in photometry of Neptune collected by the K2 mission that may be due to intrinsic global oscillations of the planet Neptune. To conduct this search, we developed new methods to correct for instrumental effects such as intrapixel variability and gain variations. We then extracted and analyzed the time-series photometry of Neptune from 49 days of nearly continuous broadband photometry of the planet. We find no evidence of global oscillations and place an upper limit of $\sim$5 ppm at 1000 \uhz\ for the detection of a coherent signal. With an observed cadence of 1-minute and point-to-point scatter less than 0.01\%, the photometric signal is dominated by reflected light from the Sun, which is in turn modulated by atmospheric variability of Neptune at the 2\% level. A change in flux is also observed due to the increasing distance between Neptune and the K2 spacecraft, and solar variability with convection-driven solar p modes present. △ Less

Submitted 8 February, 2017; originally announced February 2017.

Comments: 31 pages, 6 figure, accepted for publication in AJ

Complete waveform model for compact binaries on eccentric orbits

Authors: E. A. Huerta, Prayush Kumar, Bhanu Agarwal, Daniel George, Hsi-Yu Schive, Harald P. Pfeiffer, Roland Haas, Wei Ren, Tony Chu, Michael Boyle, Daniel A. Hemberger, Lawrence E. Kidder, Mark A. Scheel, Bela Szilagyi

Abstract: We present a time domain waveform model that describes the inspiral-merger-ringdown (IMR) of compact binary systems whose components are non-spinning, and which evolve on orbits with low to moderate eccentricity. The inspiral evolution is described using third order post-Newtonian equations both for the equations of motion of the binary, and its far-zone radiation field. This latter component also… ▽ More We present a time domain waveform model that describes the inspiral-merger-ringdown (IMR) of compact binary systems whose components are non-spinning, and which evolve on orbits with low to moderate eccentricity. The inspiral evolution is described using third order post-Newtonian equations both for the equations of motion of the binary, and its far-zone radiation field. This latter component also includes instantaneous, tails and tails-of-tails contributions, and a contribution due to non-linear memory. This framework reduces to the post-Newtonian approximant TaylorT4 at third post-Newtonian order in the zero eccentricity limit. To improve phase accuracy, we incorporate higher-order post-Newtonian corrections for the energy flux of quasi-circular binaries and gravitational self-force corrections to the binding energy of compact binaries. This enhanced inspiral evolution prescription is combined with an analytical prescription for the merger-ringdown evolution using a catalog of numerical relativity simulations. This IMR waveform model reproduces effective-one-body waveforms for systems with mass-ratios between 1 to 15 in the zero eccentricity limit. Using a set of eccentric numerical relativity simulations, not used during calibration, we show that our eccentric model accurately reproduces the features of eccentric compact binary coalescence throughout the merger. Using this model we show that the gravitational wave transients GW150914 and GW151226 can be effectively recovered with template banks of quasi-circular, spin-aligned waveforms if the eccentricity $e_0$ of these systems when they enter the aLIGO band at a gravitational wave frequency of 14 Hz satisfies $e_0^{\rm GW150914}\leq0.15$ and $e_0^{\rm GW151226}\leq0.1$. △ Less

Submitted 1 February, 2017; v1 submitted 19 September, 2016; originally announced September 2016.

Comments: 31 pages, 19 figures, 3 appendices. Submitted to Phys Rev D. v2: direct comparison to eccentric numerical relativity simulations, not used for calibration of the model, included. References added. Accepted to Phys Rev D

ACM Class: J.2

Journal ref: Phys. Rev. D 95, 024038 (2017)

Planetary Radio Interferometry and Doppler Experiment (PRIDE) technique: A test case of the Mars Express Phobos fly-by

Authors: Dmitry A. Duev, Sergei V. Pogrebenko, Giuseppe Cimò, Guifré Molera Calvés, Tatiana M. Bocanegra Bahamón, Leonid I. Gurvits, Mark M. Kettenis, Joseph Kania, Valeriu Tudose, Pascal Rosenblatt, Jean-Charles Marty, Valery Lainey, Pablo de Vicente, Jonathan Quick, Marisa Nickola, Alexander Neidhardt, Gerhard Kronschnabl, Christian Plötz, Rüdiger Haas, Michael Lindqvist, Andrea Orlati, Alexander V. Ipatov, Mikhail A. Kharinov, Andrey G. Mikhailov, Jim Lovell , et al. (11 additional authors not shown)

Abstract: The closest ever fly-by of the Martian moon Phobos, performed by the European Space Agency's Mars Express spacecraft, gives a unique opportunity to sharpen and test the Planetary Radio Interferometry and Doppler Experiments (PRIDE) technique in the interest of studying planet - satellite systems. The aim of this work is to demonstrate a technique of providing high precision positional and Doppler… ▽ More The closest ever fly-by of the Martian moon Phobos, performed by the European Space Agency's Mars Express spacecraft, gives a unique opportunity to sharpen and test the Planetary Radio Interferometry and Doppler Experiments (PRIDE) technique in the interest of studying planet - satellite systems. The aim of this work is to demonstrate a technique of providing high precision positional and Doppler measurements of planetary spacecraft using the Mars Express spacecraft. The technique will be used in the framework of Planetary Radio Interferometry and Doppler Experiments in various planetary missions, in particular in fly-by mode. We advanced a novel approach to spacecraft data processing using the techniques of Doppler and phase-referenced very long baseline interferometry spacecraft tracking. We achieved, on average, mHz precision (30 μm/s at a 10 seconds integration time) for radial three-way Doppler estimates and sub-nanoradian precision for lateral position measurements, which in a linear measure (at a distance of 1.4 AU) corresponds to ~50 m. △ Less

Submitted 19 June, 2016; originally announced June 2016.

Comments: 9 pages, 14 figures. Astronomy and Astrophysics, accepted on 2016/05/31

Journal ref: A&A 593, A34 (2016)

An improved analysis of GW150914 using a fully spin-precessing waveform model

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, M. R. Abernathy, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, A. Allocca, P. A. Altin, S. B. Anderson , et al. (951 additional authors not shown)

Abstract: This paper presents updated estimates of source parameters for GW150914, a binary black-hole coalescence event detected by the Laser Interferometer Gravitational-wave Observatory (LIGO) on September 14, 2015 [1]. Reference presented parameter estimation [2] of the source using a 13-dimensional, phenomenological precessing-spin model (precessing IMRPhenom) and a 11-dimensional nonprecessing effecti… ▽ More This paper presents updated estimates of source parameters for GW150914, a binary black-hole coalescence event detected by the Laser Interferometer Gravitational-wave Observatory (LIGO) on September 14, 2015 [1]. Reference presented parameter estimation [2] of the source using a 13-dimensional, phenomenological precessing-spin model (precessing IMRPhenom) and a 11-dimensional nonprecessing effective-one-body (EOB) model calibrated to numerical-relativity simulations, which forces spin alignment (nonprecessing EOBNR). Here we present new results that include a 15-dimensional precessing-spin waveform model (precessing EOBNR) developed within the EOB formalism. We find good agreement with the parameters estimated previously [2], and we quote updated component masses of $35^{+5}_{-3}\mathrm{M}_\odot$ and $30^{+3}_{-4}\mathrm{M}_\odot$ (where errors correspond to 90% symmetric credible intervals). We also present slightly tighter constraints on the dimensionless spin magnitudes of the two black holes, with a primary spin estimate $0.65$ and a secondary spin estimate $0.75$ at 90% probability. Reference [2] estimated the systematic parameter-extraction errors due to waveform-model uncertainty by combining the posterior probability densities of precessing IMRPhenom and nonprecessing EOBNR. Here we find that the two precessing-spin models are in closer agreement, suggesting that these systematic errors are smaller than previously quoted. △ Less

Submitted 3 June, 2016; originally announced June 2016.

Comments: 18 pages, 9 figures

Report number: LIGO-P1600048

Journal ref: Phys. Rev. X 6, 041014 (2016)

RadioAstron gravitational redshift experiment: status update

Authors: D. A. Litvinov, U. Bach, N. Bartel, K. G. Belousov, M. Bietenholz, A. V. Biriukov, G. Cimo, D. A. Duev, L. I. Gurvits, A. V. Gusev, R. Haas, V. L. Kauts, B. Z. Kanevsky, A. V. Kovalenko, G. Kronschnabl, V. V. Kulagin, M. Lindqvist, G. Molera Calves, A. Neidhardt, C. Ploetz, S. V. Pogrebenko, N. K. Porayko, V. N. Rudenko, A. I. Smirnov, K. V. Sokolovsky , et al. (3 additional authors not shown)

Abstract: A test of a cornerstone of general relativity, the gravitational redshift effect, is currently being conducted with the RadioAstron spacecraft, which is on a highly eccentric orbit around Earth. Using ground radio telescopes to record the spacecraft signal, synchronized to its ultra-stable on-board H-maser, we can probe the varying flow of time on board with unprecedented accuracy. The observation… ▽ More A test of a cornerstone of general relativity, the gravitational redshift effect, is currently being conducted with the RadioAstron spacecraft, which is on a highly eccentric orbit around Earth. Using ground radio telescopes to record the spacecraft signal, synchronized to its ultra-stable on-board H-maser, we can probe the varying flow of time on board with unprecedented accuracy. The observations performed so far, currently being analyzed, have already allowed us to measure the effect with a relative accuracy of $4\times10^{-4}$. We expect to reach $2.5\times10^{-5}$ with additional observations in 2016, an improvement of almost a magnitude over the 40-year old result of the GP-A mission. △ Less

Submitted 19 May, 2016; originally announced May 2016.

Comments: Submitted to the 14th Marcel Grossmann Meeting Proceedings

Measuring Transit Signal Recovery in the Kepler Pipeline. III. Completeness of the Q1-Q17 DR24 Planet Candidate Catalogue, with Important Caveats for Occurrence Rate Calculations

Authors: Jessie L. Christiansen, Bruce D. Clarke, Christopher J. Burke, Jon M. Jenkins, Stephen T. Bryson, Jeffrey L. Coughlin, Fergal Mullally, Susan E. Thompson, Joseph D. Twicken, Natalie M. Batalha, Michael R. Haas, Joseph Catanzarite, Jennifer R. Campbell, AKM Kamal Uddin, Khadeejah Zamudio, Jeffrey C. Smith, Christopher E. Henze

Abstract: With each new version of the Kepler pipeline and resulting planet candidate catalogue, an updated measurement of the underlying planet population can only be recovered with an corresponding measurement of the Kepler pipeline detection efficiency. Here, we present measurements of the sensitivity of the pipeline (version 9.2) used to generate the Q1-Q17 DR24 planet candidate catalog (Coughlin et al.… ▽ More With each new version of the Kepler pipeline and resulting planet candidate catalogue, an updated measurement of the underlying planet population can only be recovered with an corresponding measurement of the Kepler pipeline detection efficiency. Here, we present measurements of the sensitivity of the pipeline (version 9.2) used to generate the Q1-Q17 DR24 planet candidate catalog (Coughlin et al. 2016). We measure this by injecting simulated transiting planets into the pixel-level data of 159,013 targets across the entire Kepler focal plane, and examining the recovery rate. Unlike previous versions of the Kepler pipeline, we find a strong period dependence in the measured detection efficiency, with longer (>40 day) periods having a significantly lower detectability than shorter periods, introduced in part by an incorrectly implemented veto. Consequently, the sensitivity of the 9.2 pipeline cannot be cast as a simple one-dimensional function of the signal strength of the candidate planet signal as was possible for previous versions of the pipeline. We report on the implications for occurrence rate calculations based on the Q1-Q17 DR24 planet candidate catalog and offer important caveats and recommendations for performing such calculations. As before, we make available the entire table of injected planet parameters and whether they were recovered by the pipeline, enabling readers to derive the pipeline detection sensitivity in the planet and/or stellar parameter space of their choice. △ Less

Submitted 18 May, 2016; originally announced May 2016.

Comments: 8 pages, 5 figures, full electronic version of Table 1 available at the NASA Exoplanet Archive; accepted by ApJ May 2nd, 2016

False positive probabilties for all Kepler Objects of Interest: 1284 newly validated planets and 428 likely false positives

Authors: Timothy D. Morton, Stephen T. Bryson, Jeffrey L. Coughlin, Jason F. Rowe, Ganesh Ravichandran, Erik A. Petigura, Michael R. Haas, Natalie M. Batalha

Abstract: We present astrophysical false positive probability calculations for every Kepler Object of Interest (KOI)---the first large-scale demonstration of a fully automated transiting planet validation procedure. Out of 7056 KOIs, we determine that 1935 have probabilities <1% to be astrophysical false positives, and thus may be considered validated planets. 1284 of these have not yet been validated or co… ▽ More We present astrophysical false positive probability calculations for every Kepler Object of Interest (KOI)---the first large-scale demonstration of a fully automated transiting planet validation procedure. Out of 7056 KOIs, we determine that 1935 have probabilities <1% to be astrophysical false positives, and thus may be considered validated planets. 1284 of these have not yet been validated or confirmed by other methods. In addition, we identify 428 KOIs likely to be false positives that have not yet been identified as such, though some of these may be a result of unidentified transit timing variations. A side product of these calculations is full stellar property posterior samplings for every host star, modeled as single, binary, and triple systems. These calculations use 'vespa', a publicly available Python package able to be easily applied to any transiting exoplanet candidate. △ Less

Submitted 9 May, 2016; originally announced May 2016.

Comments: 20 pages, 8 figures. Published in ApJ. Instructions to reproduce results can be found at https://github.com/timothydmorton/koi-fpp

General Relativistic Three-Dimensional Multi-Group Neutrino Radiation-Hydrodynamics Simulations of Core-Collapse Supernovae

Authors: Luke F. Roberts, Christian D. Ott, Roland Haas, Evan P. O'Connor, Peter Diener, Erik Schnetter

Abstract: We report on a set of long-term general-relativistic three-dimensional (3D) multi-group (energy-dependent) neutrino-radiation hydrodynamics simulations of core-collapse supernovae. We employ a full 3D two-moment scheme with the local M1 closure, three neutrino species, and 12 energy groups per species. With this, we follow the post-core-bounce evolution of the core of a nonrotating $27$-$M_\odot$… ▽ More We report on a set of long-term general-relativistic three-dimensional (3D) multi-group (energy-dependent) neutrino-radiation hydrodynamics simulations of core-collapse supernovae. We employ a full 3D two-moment scheme with the local M1 closure, three neutrino species, and 12 energy groups per species. With this, we follow the post-core-bounce evolution of the core of a nonrotating $27$-$M_\odot$ progenitor in full unconstrained 3D and in octant symmetry for $\gtrsim$$ 380\,\mathrm{ms}$. We find the development of an asymmetric runaway explosion in our unconstrained simulation. We test the resolution dependence of our results and, in agreement with previous work, find that low resolution artificially aids explosion and leads to an earlier runaway expansion of the shock. At low resolution, the octant and full 3D dynamics are qualitatively very similar, but at high resolution, only the full 3D simulation exhibits the onset of explosion. △ Less

Submitted 6 April, 2017; v1 submitted 26 April, 2016; originally announced April 2016.

Comments: Accepted to ApJ

Report number: YITP-16-54

Detection of Potential Transit Signals in 17 Quarters of Kepler Data: Results of the Final Kepler Mission Transiting Planet Search (DR25)

Authors: Joseph D. Twicken, Jon M. Jenkins, Shawn E. Seader, Peter Tenenbaum, Jeffrey C. Smith, Lee S. Brownston, Christopher J. Burke, Joseph H. Catanzarite, Bruce D. Clarke, Miles T. Cote, Forrest R. Girouard, Todd C. Klaus, Jie Li, Sean D. McCauliff, Robert L. Morris, Bill Wohler, Jennifer R. Campbell, Akm Kamal Uddin, Khadeejah A. Zamudio, Anima Sabale, Steven T. Bryson, Douglas A. Caldwell, Jessie L. Christiansen, Jeffrey L. Coughlin, Michael R. Haas , et al. (3 additional authors not shown)

Abstract: We present results of the final Kepler Data Processing Pipeline search for transiting planet signals in the full 17-quarter primary mission data set. The search includes a total of 198,709 stellar targets, of which 112,046 were observed in all 17 quarters and 86,663 in fewer than 17 quarters. We report on 17,230 targets for which at least one transit signature is identified that meets the specifie… ▽ More We present results of the final Kepler Data Processing Pipeline search for transiting planet signals in the full 17-quarter primary mission data set. The search includes a total of 198,709 stellar targets, of which 112,046 were observed in all 17 quarters and 86,663 in fewer than 17 quarters. We report on 17,230 targets for which at least one transit signature is identified that meets the specified detection criteria: periodicity, minimum of three observed transit events, detection statistic (i.e., signal-to-noise ratio) in excess of the search threshold, and passing grade on three statistical transit consistency tests. Light curves for which a transit signal is identified are iteratively searched for additional signatures after a limb-darkened transiting planet model is fitted to the data and transit events are removed. The search for additional planets adds 16,802 transit signals for a total of 34,032; this far exceeds the number of transit signatures identified in prior pipeline runs. There was a strategic emphasis on completeness over reliability for the final Kepler transit search. A comparison of the transit signals against a set of 3402 well-established, high-quality Kepler Objects of Interest yields a recovery rate of 99.8%. The high recovery rate must be weighed against a large number of false-alarm detections. We examine characteristics of the planet population implied by the transiting planet model fits with an emphasis on detections that would represent small planets orbiting in the habitable zone of their host stars. △ Less

Submitted 16 November, 2016; v1 submitted 20 April, 2016; originally announced April 2016.

Comments: 60 pages, 20 figures, 3 tables. Submitted to AAS Journals on 4/19/16. Resubmitted following referee report on 6/23/16. Accepted for publication in AJ on 6/30/16. Published in AJ on 11/11/16

Identifying False Alarms in the Kepler Planet Candidate Catalog

Authors: F. Mullally, Jeffery L. Coughlin, Susan E. Thompson, Jessie Christiansen, Christopher Burke, Bruce D. Clarke, Michael R. Haas

Abstract: We present a new automated method to identify instrumental features masquerading as small, long period planets in the \kepler\ planet candidate catalog. These systematics, mistakenly identified as planet transits, can have a strong impact on occurrence rate calculations because they cluster in a region of parameter space where Kepler's sensitivity to planets is poor. We compare individual transit-… ▽ More We present a new automated method to identify instrumental features masquerading as small, long period planets in the \kepler\ planet candidate catalog. These systematics, mistakenly identified as planet transits, can have a strong impact on occurrence rate calculations because they cluster in a region of parameter space where Kepler's sensitivity to planets is poor. We compare individual transit-like events to a variety of models of real transits and systematic events, and use a Bayesian Information Criterion to evaluate the likelihood that each event is real. We describe our technique and test its performance on simulated data. Results from this technique are incorporated in the \kepler\ Q1-17 DR24 planet candidate catalog of \citet{Coughlin15}. △ Less

Submitted 9 February, 2016; originally announced February 2016.

Comments: 6 pages. Accepted for publication in PASP

Campaign 9 of the $K2$ Mission: Observational Parameters, Scientific Drivers, and Community Involvement for a Simultaneous Space- and Ground-based Microlensing Survey

Authors: Calen B. Henderson, Radosław Poleski, Matthew Penny, Rachel A. Street, David P. Bennett, David W. Hogg, B. Scott Gaudi, W. Zhu, T. Barclay, G. Barentsen, S. B. Howell, F. Mullally, A. Udalski, M. K. Szymański, J. Skowron, P. Mróz, S. Kozłowski, Ł. Wyrzykowski, P. Pietrukowicz, I. Soszyński, K. Ulaczyk, M. Pawlak, T. Sumi, F. Abe, Y. Asakura , et al. (96 additional authors not shown)

Abstract: $K2$'s Campaign 9 ($K2$C9) will conduct a $\sim$3.7 deg$^{2}$ survey toward the Galactic bulge from 7/April through 1/July of 2016 that will leverage the spatial separation between $K2$ and the Earth to facilitate measurement of the microlens parallax $π_{\rm E}$ for $\gtrsim… ▽ More $K2$'s Campaign 9 ($K2$C9) will conduct a $\sim$3.7 deg$^{2}$ survey toward the Galactic bulge from 7/April through 1/July of 2016 that will leverage the spatial separation between $K2$ and the Earth to facilitate measurement of the microlens parallax $π_{\rm E}$ for $\gtrsim$127 microlensing events. These will include several that are planetary in nature as well as many short-timescale microlensing events, which are potentially indicative of free-floating planets (FFPs). These satellite parallax measurements will in turn allow for the direct measurement of the masses of and distances to the lensing systems. In this white paper we provide an overview of the $K2$C9 space- and ground-based microlensing survey. Specifically, we detail the demographic questions that can be addressed by this program, including the frequency of FFPs and the Galactic distribution of exoplanets, the observational parameters of $K2$C9, and the array of resources dedicated to concurrent observations. Finally, we outline the avenues through which the larger community can become involved, and generally encourage participation in $K2$C9, which constitutes an important pathfinding mission and community exercise in anticipation of $WFIRST$. △ Less

Submitted 7 March, 2016; v1 submitted 30 December, 2015; originally announced December 2015.

Comments: 19 pages, 11 figures, 3 tables; submitted to PASP

Planetary Candidates Observed by Kepler. VII. The First Fully Uniform Catalog Based on The Entire 48 Month Dataset (Q1-Q17 DR24)

Authors: Jeffrey L. Coughlin, F. Mullally, Susan E. Thompson, Jason F. Rowe, Christopher J. Burke, David W. Latham, Natalie M. Batalha, Aviv Ofir, Billy L. Quarles, Christopher E. Henze, Angie Wolfgang, Douglas A. Caldwell, Stephen T. Bryson, Avi Shporer, Joseph Catanzarite, Rachel Akeson, Thomas Barclay, William J. Borucki, Tabetha S. Boyajian, Jennifer R. Campbell, Jessie L. Christiansen, Forrest R. Girouard, Michael R. Haas, Steve B. Howell, Daniel Huber , et al. (10 additional authors not shown)

Abstract: We present the seventh Kepler planet candidate catalog, which is the first to be based on the entire, uniformly processed, 48 month Kepler dataset. This is the first fully automated catalog, employing robotic vetting procedures to uniformly evaluate every periodic signal detected by the Q1-Q17 Data Release 24 (DR24) Kepler pipeline. While we prioritize uniform vetting over the absolute correctness… ▽ More We present the seventh Kepler planet candidate catalog, which is the first to be based on the entire, uniformly processed, 48 month Kepler dataset. This is the first fully automated catalog, employing robotic vetting procedures to uniformly evaluate every periodic signal detected by the Q1-Q17 Data Release 24 (DR24) Kepler pipeline. While we prioritize uniform vetting over the absolute correctness of individual objects, we find that our robotic vetting is overall comparable to, and in most cases is superior to, the human vetting procedures employed by past catalogs. This catalog is the first to utilize artificial transit injection to evaluate the performance of our vetting procedures and quantify potential biases, which are essential for accurate computation of planetary occurrence rates. With respect to the cumulative Kepler Object of Interest (KOI) catalog, we designate 1,478 new KOIs, of which 402 are dispositioned as planet candidates (PCs). Also, 237 KOIs dispositioned as false positives (FPs) in previous Kepler catalogs have their disposition changed to PC and 118 PCs have their disposition changed to FP. This brings the total number of known KOIs to 8,826 and PCs to 4,696. We compare the Q1-Q17 DR24 KOI catalog to previous KOI catalogs, as well as ancillary Kepler catalogs, finding good agreement between them. We highlight new PCs that are both potentially rocky and potentially in the habitable zone of their host stars, many of which orbit solar-type stars. This work represents significant progress in accurately determining the fraction of Earth-size planets in the habitable zone of Sun-like stars. The full catalog is publicly available at the NASA Exoplanet Archive. △ Less

Submitted 18 February, 2016; v1 submitted 18 December, 2015; originally announced December 2015.

Comments: Accepted to the Astrophysical Journal Supplement Series. 30 pages, 9 figures, 7 tables. We make the DR24 robovetter decision code publicly available at http://github.com/JeffLCoughlin/robovetter, with input and output examples provided using the same data as contained in the full paper's tables

The K2 Ecliptic Plane Input Catalog (EPIC) and Stellar Classifications of 138,600 Targets in Campaigns 1-8

Authors: Daniel Huber, Stephen T. Bryson, Michael R. Haas, Thomas Barclay, Geert Barentsen, Steve B. Howell, Sanjib Sharma, Dennis Stello, Susan E. Thompson

Abstract: The K2 Mission uses the Kepler spacecraft to obtain high-precision photometry over ~80 day campaigns in the ecliptic plane. The Ecliptic Plane Input Catalog (EPIC) provides coordinates, photometry and kinematics based on a federation of all-sky catalogs to support target selection and target management for the K2 mission. We describe the construction of the EPIC, as well as modifications and short… ▽ More The K2 Mission uses the Kepler spacecraft to obtain high-precision photometry over ~80 day campaigns in the ecliptic plane. The Ecliptic Plane Input Catalog (EPIC) provides coordinates, photometry and kinematics based on a federation of all-sky catalogs to support target selection and target management for the K2 mission. We describe the construction of the EPIC, as well as modifications and shortcomings of the catalog. Kepler magnitudes (Kp) are shown to be accurate to ~0.1 mag for the Kepler field, and the EPIC is typically complete to Kp~17 (Kp~19 for campaigns covered by SDSS). We furthermore classify 138,600 targets in Campaigns 1-8 (~88% of the full target sample) using colors, proper motions, spectroscopy, parallaxes, and galactic population synthesis models, with typical uncertainties for G-type stars of ~3% in Teff, ~0.3 dex in log(g), ~40% in radius, ~10% in mass, and ~40% in distance. Our results show that stars targeted by K2 are dominated by K-M dwarfs (~41% of all selected targets), F-G dwarfs (~36%) and K giants (~21%), consistent with key K2 science programs to search for transiting exoplanets and galactic archeology studies using oscillating red giants. However, we find a significant variation of the fraction of cool dwarfs with galactic latitude, indicating a target selection bias due to interstellar reddening and the increased contamination by giant stars near the galactic plane. We discuss possible systematic errors in the derived stellar properties, and differences to published classifications for K2 exoplanet host stars. The EPIC is hosted at the Mikulski Archive for Space Telescopes (MAST): http://archive.stsci.edu/k2/epic/search.php. △ Less

Submitted 19 May, 2016; v1 submitted 8 December, 2015; originally announced December 2015.

Comments: 19 pages, 12 figures, 5 tables, accepted for publication in ApJS. An electronic version of Table 5 is available as an ancillary file (sidebar on the right), and source codes are available at https://github.com/danxhuber/k2epic and https://github.com/danxhuber/galclassify v3: minor text changes and updated uncertainties in Table 5; v4: minor text changes to match published version

A large scale dynamo and magnetoturbulence in rapidly rotating core-collapse supernovae

Authors: Philipp Mösta, Christian D. Ott, David Radice, Luke F. Roberts, Erik Schnetter, Roland Haas

Abstract: Magnetohydrodynamic (MHD) turbulence is of key importance in many high-energy astrophysical systems, where MHD instabilities can amplify local magnetic field over very short time scales. Specifically, the magnetorotational instability (MRI) and dynamo action have been suggested as a mechanism to grow magnetar-strength magnetic field ($\ge10^{15} G$) and magnetorotationally power the explosion of a… ▽ More Magnetohydrodynamic (MHD) turbulence is of key importance in many high-energy astrophysical systems, where MHD instabilities can amplify local magnetic field over very short time scales. Specifically, the magnetorotational instability (MRI) and dynamo action have been suggested as a mechanism to grow magnetar-strength magnetic field ($\ge10^{15} G$) and magnetorotationally power the explosion of a rotating massive star. Such stars are progenitor candidates for type Ic-bl hypernova explosions and make up all supernovae connected to long gamma-ray bursts (GRBs). The MRI has been studied with local high-resolution shearing box simulations in 3D or with global 2D simulations, but it is an open question whether MRI-driven turbulence can result in the creation of a large-scale ordered and dynamically relevant field. Here we report results from global 3D general-relativistic magnetohydrodynamic (GRMHD) turbulence simulations and show that MRI-driven MHD turbulence in rapidly rotating protoneutron stars produces an inverse cascade of energy. We find a large-scale ordered toroidal field that is consistent with the formation of bipolar magnetorotationally driven outflows. Our results demonstrate that rapidly rotating massive stars are plausible progenitors for both type Ic-bl supernovae and long GRBs, present a viable formation scenario for magnetars, and may account for potentially magnetar-powered superluminous supernovae. △ Less

Submitted 28 March, 2018; v1 submitted 2 December, 2015; originally announced December 2015.

Comments: Author's accepted version of Nature letter DOI:10.1038/nature15755, published online, Nature, Nov 30, 2015

Journal ref: Volume 528, pp. 376-379 (2015)

Low mass binary neutron star mergers : gravitational waves and neutrino emission

Authors: Francois Foucart, Roland Haas, Matthew D. Duez, Evan O'Connor, Christian D. Ott, Luke Roberts, Lawrence E. Kidder, Jonas Lippuner, Harald P. Pfeiffer, Mark A. Scheel

Abstract: Neutron star mergers are among the most promising sources of gravitational waves for advanced ground-based detectors. These mergers are also expected to power bright electromagnetic signals, in the form of short gamma-ray bursts, infrared/optical transients, and radio emission. Simulations of these mergers with fully general relativistic codes are critical to understand the merger and post-merger… ▽ More Neutron star mergers are among the most promising sources of gravitational waves for advanced ground-based detectors. These mergers are also expected to power bright electromagnetic signals, in the form of short gamma-ray bursts, infrared/optical transients, and radio emission. Simulations of these mergers with fully general relativistic codes are critical to understand the merger and post-merger gravitational wave signals and their neutrinos and electromagnetic counterparts. In this paper, we employ the SpEC code to simulate the merger of low-mass neutron star binaries (two $1.2M_\odot$ neutron stars) for a set of three nuclear-theory based, finite temperature equations of state. We show that the frequency peaks of the post-merger gravitational wave signal are in good agreement with predictions obtained from simulations using a simpler treatment of gravity. We find, however, that only the fundamental mode of the remnant is excited for long periods of time: emission at the secondary peaks is damped on a millisecond timescale in the simulated binaries. For such low-mass systems, the remnant is a massive neutron star which, depending on the equation of state, is either permanently stable or long-lived. We observe strong excitations of l=2, m=2 modes, both in the massive neutron star and in the form of hot, shocked tidal arms in the surrounding accretion torus. We estimate the neutrino emission of the remnant using a neutrino leakage scheme and, in one case, compare these results with a gray two-moment neutrino transport scheme. We confirm the complex geometry of the neutrino emission, also observed in previous simulations with neutrino leakage, and show explicitly the presence of important differences in the neutrino luminosity, disk composition, and outflow properties between the neutrino leakage and transport schemes. △ Less

Submitted 22 January, 2016; v1 submitted 21 October, 2015; originally announced October 2015.

Comments: Accepted by PRD; 23 pages; 24 figures; 4 tables

Journal ref: Phys. Rev. D 93, 044019 (2016)

Neutrino-Driven Convection in Core-Collapse Supernovae: High-Resolution Simulations

Authors: David Radice, Christian D. Ott, Ernazar Abdikamalov, Sean M. Couch, Roland Haas, Erik Schnetter

Abstract: We present results from high-resolution semi-global simulations of neutrino-driven convection in core-collapse supernovae. We employ an idealized setup with parametrized neutrino heating/cooling and nuclear dissociation at the shock front. We study the internal dynamics of neutrino-driven convection and its role in re-distributing energy and momentum through the gain region. We find that even if b… ▽ More We present results from high-resolution semi-global simulations of neutrino-driven convection in core-collapse supernovae. We employ an idealized setup with parametrized neutrino heating/cooling and nuclear dissociation at the shock front. We study the internal dynamics of neutrino-driven convection and its role in re-distributing energy and momentum through the gain region. We find that even if buoyant plumes are able to locally transfer heat up to the shock, convection is not able to create a net positive energy flux and overcome the downwards transport of energy from the accretion flow. Turbulent convection does, however, provide a significant effective pressure support to the accretion flow as it favors the accumulation of energy, mass and momentum in the gain region. We derive an approximate equation that is able to explain and predict the shock evolution in terms of integrals of quantities such as the turbulent pressure in the gain region or the effects of non-radial motion of the fluid. We use this relation as a way to quantify the role of turbulence in the dynamics of the accretion shock. Finally, we investigate the effects of grid resolution, which we change by a factor 20 between the lowest and highest resolution. Our results show that the shallow slopes of the turbulent kinetic energy spectra reported in previous studies are a numerical artefact. Kolmogorov scaling is progressively recovered as the resolution is increased. △ Less

Submitted 22 March, 2016; v1 submitted 16 October, 2015; originally announced October 2015.

Comments: 19 pages, 14 figures. Version accepted on ApJ. An animated version of Fig. 3 is available at: https://youtu.be/yLub83WP3WA

Journal ref: The Astrophysical Journal, 820, 76 (2016)

A Machine Learning Technique to Identify Transit Shaped Signals

Authors: Susan E. Thompson, Fergal Mullally, Jeff Coughlin, Jessie L. Christiansen, Christopher E. Henze, Michael R. Haas, Christopher J. Burke

Abstract: We describe a new metric that uses machine learning to determine if a periodic signal found in a photometric time series appears to be shaped like the signature of a transiting exoplanet. This metric uses dimensionality reduction and k-nearest neighbors to determine whether a given signal is sufficiently similar to known transits in the same data set. This metric is being used by the Kepler Robove… ▽ More We describe a new metric that uses machine learning to determine if a periodic signal found in a photometric time series appears to be shaped like the signature of a transiting exoplanet. This metric uses dimensionality reduction and k-nearest neighbors to determine whether a given signal is sufficiently similar to known transits in the same data set. This metric is being used by the Kepler Robovetter to determine which signals should be part of the Q1-Q17 DR24 catalog of planetary candidates. The Kepler Mission reports roughly 20,000 potential transiting signals with each run of its pipeline, yet only a few thousand appear sufficiently transit shaped to be part of the catalog. The other signals tend to be variable stars and instrumental noise. With this metric we are able to remove more than 90% of the non-transiting signals while retaining more than 99% of the known planet candidates. When tested with injected transits, less than 1% are lost. This metric will enable the Kepler mission and future missions looking for transiting planets to rapidly and consistently find the best planetary candidates for follow-up and cataloging. △ Less

Submitted 31 August, 2015; originally announced September 2015.

Comments: Accepted for Publication to The Astrophysical Journal