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

SciCode: A Research Coding Benchmark Curated by Scientists

Authors: Minyang Tian, Luyu Gao, Shizhuo Dylan Zhang, Xinan Chen, Cunwei Fan, Xuefei Guo, Roland Haas, Pan Ji, Kittithat Krongchon, Yao Li, Shengyan Liu, Di Luo, Yutao Ma, Hao Tong, Kha Trinh, Chenyu Tian, Zihan Wang, Bohao Wu, Yanyu Xiong, Shengzhu Yin, Minhui Zhu, Kilian Lieret, Yanxin Lu, Genglin Liu, Yufeng Du , et al. (5 additional authors not shown)

Abstract: Since language models (LMs) now outperform average humans on many challenging tasks, it has become increasingly difficult to develop challenging, high-quality, and realistic evaluations. We address this issue by examining LMs' capabilities to generate code for solving real scientific research problems. Incorporating input from scientists and AI researchers in 16 diverse natural science sub-fields,… ▽ More Since language models (LMs) now outperform average humans on many challenging tasks, it has become increasingly difficult to develop challenging, high-quality, and realistic evaluations. We address this issue by examining LMs' capabilities to generate code for solving real scientific research problems. Incorporating input from scientists and AI researchers in 16 diverse natural science sub-fields, including mathematics, physics, chemistry, biology, and materials science, we created a scientist-curated coding benchmark, SciCode. The problems in SciCode naturally factorize into multiple subproblems, each involving knowledge recall, reasoning, and code synthesis. In total, SciCode contains 338 subproblems decomposed from 80 challenging main problems. It offers optional descriptions specifying useful scientific background information and scientist-annotated gold-standard solutions and test cases for evaluation. Claude3.5-Sonnet, the best-performing model among those tested, can solve only 4.6% of the problems in the most realistic setting. We believe that SciCode demonstrates both contemporary LMs' progress towards becoming helpful scientific assistants and sheds light on the development and evaluation of scientific AI in the future. △ Less

Submitted 18 July, 2024; originally announced July 2024.

Comments: 25 pages, 9 figures, 7 tables

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

Well-forced graphs

Authors: Cheryl Grood, Ruth Haas, Bonnie Jacob, Erika King, Shahla Nasserasr

Abstract: A graph in which all minimal zero forcing sets are in fact minimum size is called ``well-forced." This paper characterizes well-forced trees and presents an algorithm for determining which trees are well-forced. Additionally, we characterize which vertices in a tree are contained in no minimal zero forcing set. A graph in which all minimal zero forcing sets are in fact minimum size is called ``well-forced." This paper characterizes well-forced trees and presents an algorithm for determining which trees are well-forced. Additionally, we characterize which vertices in a tree are contained in no minimal zero forcing set. △ Less

Submitted 21 December, 2023; originally announced December 2023.

MSC Class: 05C50

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

APACE: AlphaFold2 and advanced computing as a service for accelerated discovery in biophysics

Authors: Hyun Park, Parth Patel, Roland Haas, E. A. Huerta

Abstract: The prediction of protein 3D structure from amino acid sequence is a computational grand challenge in biophysics, and plays a key role in robust protein structure prediction algorithms, from drug discovery to genome interpretation. The advent of AI models, such as AlphaFold, is revolutionizing applications that depend on robust protein structure prediction algorithms. To maximize the impact, and e… ▽ More The prediction of protein 3D structure from amino acid sequence is a computational grand challenge in biophysics, and plays a key role in robust protein structure prediction algorithms, from drug discovery to genome interpretation. The advent of AI models, such as AlphaFold, is revolutionizing applications that depend on robust protein structure prediction algorithms. To maximize the impact, and ease the usability, of these novel AI tools we introduce APACE, AlphaFold2 and advanced computing as a service, a novel computational framework that effectively handles this AI model and its TB-size database to conduct accelerated protein structure prediction analyses in modern supercomputing environments. We deployed APACE in the Delta and Polaris supercomputers, and quantified its performance for accurate protein structure predictions using four exemplar proteins: 6AWO, 6OAN, 7MEZ, and 6D6U. Using up to 300 ensembles, distributed across 200 NVIDIA A100 GPUs, we found that APACE is up to two orders of magnitude faster than off-the-self AlphaFold2 implementations, reducing time-to-solution from weeks to minutes. This computational approach may be readily linked with robotics laboratories to automate and accelerate scientific discovery. △ Less

Submitted 1 July, 2024; v1 submitted 15 August, 2023; originally announced August 2023.

Comments: 7 pages, 4 figures, 2 tables

ACM Class: I.2

Journal ref: Proceedings of the National Academy of Sciences, 121, 27, (2024)

HPC-driven computational reproducibility in numerical relativity codes: A use case study with IllinoisGRMHD

Authors: Yufeng Luo, Qian Zhang, Roland Haas, Zachariah B. Etienne, Gabrielle Allen

Abstract: Reproducibility of results is a cornerstone of the scientific method. Scientific computing encounters two challenges when aiming for this goal. Firstly, reproducibility should not depend on details of the runtime environment, such as the compiler version or computing environment, so results are verifiable by third-parties. Secondly, different versions of software code executed in the same runtime… ▽ More Reproducibility of results is a cornerstone of the scientific method. Scientific computing encounters two challenges when aiming for this goal. Firstly, reproducibility should not depend on details of the runtime environment, such as the compiler version or computing environment, so results are verifiable by third-parties. Secondly, different versions of software code executed in the same runtime environment should produce consistent numerical results for physical quantities. In this manuscript, we test the feasibility of reproducing scientific results obtained using the IllinoisGRMHD code that is part of an open-source community software for simulation in relativistic astrophysics, the Einstein Toolkit. We verify that numerical results of simulating a single isolated neutron star with IllinoisGRMHD can be reproduced, and compare them to results reported by the code authors in 2015. We use two different supercomputers: Expanse at SDSC, and Stampede2 at TACC. By compiling the source code archived along with the paper on both Expanse and Stampede2, we find that IllinoisGRMHD reproduces results published in its announcement paper up to errors comparable to round-off level changes in initial data parameters. We also verify that a current version of IlliinoisGRMHD reproduces these results once we account for bug fixes which has occurred since the original publication △ Less

Submitted 8 December, 2023; v1 submitted 3 July, 2023; originally announced July 2023.

Comments: 23 pages, 6 figures, accepted to Classical and Quantum Gravity

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

Discovering Interpretable Directions in the Semantic Latent Space of Diffusion Models

Authors: René Haas, Inbar Huberman-Spiegelglas, Rotem Mulayoff, Stella Graßhof, Sami S. Brandt, Tomer Michaeli

Abstract: Denoising Diffusion Models (DDMs) have emerged as a strong competitor to Generative Adversarial Networks (GANs). However, despite their widespread use in image synthesis and editing applications, their latent space is still not as well understood. Recently, a semantic latent space for DDMs, coined `$h$-space', was shown to facilitate semantic image editing in a way reminiscent of GANs. The $h$-spa… ▽ More Denoising Diffusion Models (DDMs) have emerged as a strong competitor to Generative Adversarial Networks (GANs). However, despite their widespread use in image synthesis and editing applications, their latent space is still not as well understood. Recently, a semantic latent space for DDMs, coined `$h$-space', was shown to facilitate semantic image editing in a way reminiscent of GANs. The $h$-space is comprised of the bottleneck activations in the DDM's denoiser across all timesteps of the diffusion process. In this paper, we explore the properties of h-space and propose several novel methods for finding meaningful semantic directions within it. We start by studying unsupervised methods for revealing interpretable semantic directions in pretrained DDMs. Specifically, we show that global latent directions emerge as the principal components in the latent space. Additionally, we provide a novel method for discovering image-specific semantic directions by spectral analysis of the Jacobian of the denoiser w.r.t. the latent code. Next, we extend the analysis by finding directions in a supervised fashion in unconditional DDMs. We demonstrate how such directions can be found by relying on either a labeled data set of real images or by annotating generated samples with a domain-specific attribute classifier. We further show how to semantically disentangle the found direction by simple linear projection. Our approaches are applicable without requiring any architectural modifications, text-based guidance, CLIP-based optimization, or model fine-tuning. △ Less

Submitted 29 May, 2024; v1 submitted 20 March, 2023; originally announced March 2023.

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

Experiments on Generalizability of BERTopic on Multi-Domain Short Text

Authors: Muriël de Groot, Mohammad Aliannejadi, Marcel R. Haas

Abstract: Topic modeling is widely used for analytically evaluating large collections of textual data. One of the most popular topic techniques is Latent Dirichlet Allocation (LDA), which is flexible and adaptive, but not optimal for e.g. short texts from various domains. We explore how the state-of-the-art BERTopic algorithm performs on short multi-domain text and find that it generalizes better than LDA i… ▽ More Topic modeling is widely used for analytically evaluating large collections of textual data. One of the most popular topic techniques is Latent Dirichlet Allocation (LDA), which is flexible and adaptive, but not optimal for e.g. short texts from various domains. We explore how the state-of-the-art BERTopic algorithm performs on short multi-domain text and find that it generalizes better than LDA in terms of topic coherence and diversity. We further analyze the performance of the HDBSCAN clustering algorithm utilized by BERTopic and find that it classifies a majority of the documents as outliers. This crucial, yet overseen problem excludes too many documents from further analysis. When we replace HDBSCAN with k-Means, we achieve similar performance, but without outliers. △ Less

Submitted 16 December, 2022; originally announced December 2022.

Comments: Accepted poster presentation at WiNLP 2022, as a part of EMNLP 2022, 2 pages

Controllable GAN Synthesis Using Non-Rigid Structure-from-Motion

Authors: René Haas, Stella Graßhof, Sami S. Brandt

Abstract: In this paper, we present an approach for combining non-rigid structure-from-motion (NRSfM) with deep generative models,and propose an efficient framework for discovering trajectories in the latent space of 2D GANs corresponding to changes in 3D geometry. Our approach uses recent advances in NRSfM and enables editing of the camera and non-rigid shape information associated with the latent codes wi… ▽ More In this paper, we present an approach for combining non-rigid structure-from-motion (NRSfM) with deep generative models,and propose an efficient framework for discovering trajectories in the latent space of 2D GANs corresponding to changes in 3D geometry. Our approach uses recent advances in NRSfM and enables editing of the camera and non-rigid shape information associated with the latent codes without needing to retrain the generator. This formulation provides an implicit dense 3D reconstruction as it enables the image synthesis of novel shapes from arbitrary view angles and non-rigid structure. The method is built upon a sparse backbone, where a neural regressor is first trained to regress parameters describing the cameras and sparse non-rigid structure directly from the latent codes. The latent trajectories associated with changes in the camera and structure parameters are then identified by estimating the local inverse of the regressor in the neighborhood of a given latent code. The experiments show that our approach provides a versatile, systematic way to model, analyze, and edit the geometry and non-rigid structures of faces. △ Less

Submitted 14 November, 2022; originally announced November 2022.

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

Numerical relativity higher order gravitational waveforms of eccentric, spinning, non-precessing binary black hole mergers

Authors: Abhishek V. Joshi, Shawn G. Rosofsky, Roland Haas, E. A. Huerta

Abstract: We use the open source, community-driven, numerical relativity software, the Einstein Toolkit to study the physics of eccentric, spinning, nonprecessing binary black hole mergers with mass-ratios $q=\{2, 4, 6\}$, individual dimensionless spin parameters $χ_{1z}=\pm0.6$, $χ_{2z}=\pm0.3$, that include higher order gravitational wave modes $\ell\leq4$, except for memory modes. Assuming stellar mass b… ▽ More We use the open source, community-driven, numerical relativity software, the Einstein Toolkit to study the physics of eccentric, spinning, nonprecessing binary black hole mergers with mass-ratios $q=\{2, 4, 6\}$, individual dimensionless spin parameters $χ_{1z}=\pm0.6$, $χ_{2z}=\pm0.3$, that include higher order gravitational wave modes $\ell\leq4$, except for memory modes. Assuming stellar mass binary black hole mergers that may be detectable by the advanced LIGO detectors, we find that including modes up to $\ell=4$ increases the signal-to-noise of compact binaries between $3.5\%$ to $35\%$, compared to signals that only include the $\ell=|m|=2$ mode. We use two waveform models, TEOBResumS and SEOBNRE, which incorporate spin and eccentricity corrections in the waveform dynamics, to quantify the orbital eccentricity of our numerical relativity catalog in a gauge-invariant manner through fitting factor calculations. Our findings indicate that the inclusion of higher order wave modes has a measurable effect in the recovery of moderately and highly eccentric black hole mergers, and thus it is essential to develop waveform models and signal processing tools that accurately describe the physics of these astrophysical sources. △ Less

Submitted 21 March, 2023; v1 submitted 4 October, 2022; originally announced October 2022.

Comments: 13 pages, 5 figures, updated manuscript with added waveform comparisons to TEOBResumS. accepted in Phys. Rev. D

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

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)

Isomorphisms and properties of TAR reconfiguration graphs for zero forcing and other $X$-set parameters

Authors: Novi H. Bong, Joshua Carlson, Bryan Curtis, Ruth Haas, Leslie Hogben

Abstract: An $X$-TAR (token addition/removal) reconfiguration graph has as its vertices sets that satisfy some property $X$, with an edge between two sets if one is obtained from the other by adding or removing one element. This paper considers the $X$-TAR graph for $X-$ sets of vertices of a base graph $G$ where the $X$-sets of $G$ must satisfy certain conditions. Dominating sets, power dominating sets, ze… ▽ More An $X$-TAR (token addition/removal) reconfiguration graph has as its vertices sets that satisfy some property $X$, with an edge between two sets if one is obtained from the other by adding or removing one element. This paper considers the $X$-TAR graph for $X-$ sets of vertices of a base graph $G$ where the $X$-sets of $G$ must satisfy certain conditions. Dominating sets, power dominating sets, zero forcing sets, and positive semidefinite zero forcing sets are all examples of $X$-sets. For graphs $G$ and $G'$ with no isolated vertices, it is shown that $G$ and $G'$ have isomorphic $X$-TAR reconfiguration graphs if and only if there is a relabeling of the vertices of $G'$ such that $G$ and $G'$ have exactly the same $X$-sets. The concept of an $X$-irrelevant vertex is introduced to facilitate analysis of $X$-TAR graph isomorphisms. Furthermore, results related to the connectedness of the zero forcing TAR graph are given. We present families of graphs that exceed known lower bounds for connectedness parameters. △ Less

Submitted 19 May, 2022; originally announced May 2022.

Tensor-based Emotion Editing in the StyleGAN Latent Space

Authors: René Haas, Stella Graßhof, Sami S. Brandt

Abstract: In this paper, we use a tensor model based on the Higher-Order Singular Value Decomposition (HOSVD) to discover semantic directions in Generative Adversarial Networks. This is achieved by first embedding a structured facial expression database into the latent space using the e4e encoder. Specifically, we discover directions in latent space corresponding to the six prototypical emotions: anger, dis… ▽ More In this paper, we use a tensor model based on the Higher-Order Singular Value Decomposition (HOSVD) to discover semantic directions in Generative Adversarial Networks. This is achieved by first embedding a structured facial expression database into the latent space using the e4e encoder. Specifically, we discover directions in latent space corresponding to the six prototypical emotions: anger, disgust, fear, happiness, sadness, and surprise, as well as a direction for yaw rotation. These latent space directions are employed to change the expression or yaw rotation of real face images. We compare our found directions to similar directions found by two other methods. The results show that the visual quality of the resultant edits are on par with State-of-the-Art. It can also be concluded that the tensor-based model is well suited for emotion and yaw editing, i.e., that the emotion or yaw rotation of a novel face image can be robustly changed without a significant effect on identity or other attributes in the images. △ Less

Submitted 12 May, 2022; originally announced May 2022.

Comments: Accepted for poster presentation at AI4CC @ CVPRW

Federated Learning Enables Big Data for Rare Cancer Boundary Detection

Authors: Sarthak Pati, Ujjwal Baid, Brandon Edwards, Micah Sheller, Shih-Han Wang, G Anthony Reina, Patrick Foley, Alexey Gruzdev, Deepthi Karkada, Christos Davatzikos, Chiharu Sako, Satyam Ghodasara, Michel Bilello, Suyash Mohan, Philipp Vollmuth, Gianluca Brugnara, Chandrakanth J Preetha, Felix Sahm, Klaus Maier-Hein, Maximilian Zenk, Martin Bendszus, Wolfgang Wick, Evan Calabrese, Jeffrey Rudie, Javier Villanueva-Meyer , et al. (254 additional authors not shown)

Abstract: Although machine learning (ML) has shown promise in numerous domains, there are concerns about generalizability to out-of-sample data. This is currently addressed by centrally sharing ample, and importantly diverse, data from multiple sites. However, such centralization is challenging to scale (or even not feasible) due to various limitations. Federated ML (FL) provides an alternative to train acc… ▽ More Although machine learning (ML) has shown promise in numerous domains, there are concerns about generalizability to out-of-sample data. This is currently addressed by centrally sharing ample, and importantly diverse, data from multiple sites. However, such centralization is challenging to scale (or even not feasible) due to various limitations. Federated ML (FL) provides an alternative to train accurate and generalizable ML models, by only sharing numerical model updates. Here we present findings from the largest FL study to-date, involving data from 71 healthcare institutions across 6 continents, to generate an automatic tumor boundary detector for the rare disease of glioblastoma, utilizing the largest dataset of such patients ever used in the literature (25,256 MRI scans from 6,314 patients). We demonstrate a 33% improvement over a publicly trained model to delineate the surgically targetable tumor, and 23% improvement over the tumor's entire extent. We anticipate our study to: 1) enable more studies in healthcare informed by large and diverse data, ensuring meaningful results for rare diseases and underrepresented populations, 2) facilitate further quantitative analyses for glioblastoma via performance optimization of our consensus model for eventual public release, and 3) demonstrate the effectiveness of FL at such scale and task complexity as a paradigm shift for multi-site collaborations, alleviating the need for data sharing. △ Less

Submitted 25 April, 2022; v1 submitted 22 April, 2022; originally announced April 2022.

Comments: federated learning, deep learning, convolutional neural network, segmentation, brain tumor, glioma, glioblastoma, FeTS, BraTS

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

Penning-trap mass measurements of the deuteron and the HD+ molecular ion

Authors: Sascha Rau, Fabian Heiße, Florian Köhler-Langes, Sangeetha Sasidharan, Raphael Haas, Dennis Renisch, Christoph E. Düllmann, Wolfgang Quint, Sven Sturm, Klaus Blaum

Abstract: The masses of the lightest atomic nuclei and the electron mass are interlinked and are crucial in a wide range of research fields, with their values affecting observables in atomic, molecular and neutrino physics as well as metrology. The most precise values for these fundamental parameters come from Penning-trap mass spectrometry, which achieves relative mass uncertainties in the range of… ▽ More The masses of the lightest atomic nuclei and the electron mass are interlinked and are crucial in a wide range of research fields, with their values affecting observables in atomic, molecular and neutrino physics as well as metrology. The most precise values for these fundamental parameters come from Penning-trap mass spectrometry, which achieves relative mass uncertainties in the range of $10^{-11}$. However, redundancy checks using data from different experiments reveal significant inconsistencies in the masses of the proton ($m_p$), the deuteron ($m_d$) and helion ($m_\text{he}$), amounting to $5$ standard deviations for the term $Δ=m_p+m_d-m_{\text{he}}$, which suggests that the uncertainty of these values may have been underestimated. Here we present results from absolute mass measurements of the deuteron and the ${HD}^+$ molecular ion against $^{12}C$ as a mass reference. Our value for the deuteron $m_d=2.013\,553\,212\,535 (17)$u supersedes the precision of the literature value by a factor of $2.4$ and deviates from this by $4.8$ standard deviations. With a relative uncertainty of $8$ parts per trillion (ppt) this is the most precise mass value measured directly in atomic mass units. Furthermore, the measurement of the ${HD}^+$ molecular ion, $m({HD}^+)=3.021\,378\,241\,561\,(61)$u, not only allows for a rigorous consistency check of our measurements of the masses of the deuteron (this work) and proton, but also establishes an additional link for the masses of tritium and helium-3 to the atomic mass unit. Combined with a recent measurement of the deuteron-to-proton mass ratio the uncertainty of the reference value of $m_p$ can be reduced by a factor of three. This is a post-peer-review, pre-copyedit version of an article published in Nature. The final authenticated version is available online at https://doi.org/10.1038/s41586-020-2628-7 △ Less

Submitted 11 March, 2022; originally announced March 2022.

Journal ref: Nature vol. 585, pages 43-47 (2020)

Hamilton Paths in Dominating Graphs of Trees and Cycles

Authors: Kira Adaricheva, Heather Smith Blake, Chassidy Bozeman, Nancy E. Clarke, Ruth Haas, Margaret-Ellen Messinger, Karen Seyffarth

Abstract: The dominating graph of a graph $H$ has as its vertices all dominating sets of $H$, with an edge between two dominating sets if one can be obtained from the other by the addition or deletion of a single vertex of $H$. In this paper we prove that the dominating graph of any tree has a Hamilton path. We also show how a result about binary strings leads to a proof that the dominating graph of a cycle… ▽ More The dominating graph of a graph $H$ has as its vertices all dominating sets of $H$, with an edge between two dominating sets if one can be obtained from the other by the addition or deletion of a single vertex of $H$. In this paper we prove that the dominating graph of any tree has a Hamilton path. We also show how a result about binary strings leads to a proof that the dominating graph of a cycle on $n$ vertices has a Hamilton path if and only if $n\not\equiv 0 \pmod 4$. △ Less

Submitted 4 February, 2022; v1 submitted 8 December, 2021; originally announced December 2021.

MSC Class: 05C45; 05C05; 05C38

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)

Tensor-based Subspace Factorization for StyleGAN

Authors: René Haas, Stella Graßhof, Sami Sebastian Brandt

Abstract: In this paper, we propose $τ$GAN a tensor-based method for modeling the latent space of generative models. The objective is to identify semantic directions in latent space. To this end, we propose to fit a multilinear tensor model on a structured facial expression database, which is initially embedded into latent space. We validate our approach on StyleGAN trained on FFHQ using BU-3DFE as a struct… ▽ More In this paper, we propose $τ$GAN a tensor-based method for modeling the latent space of generative models. The objective is to identify semantic directions in latent space. To this end, we propose to fit a multilinear tensor model on a structured facial expression database, which is initially embedded into latent space. We validate our approach on StyleGAN trained on FFHQ using BU-3DFE as a structured facial expression database. We show how the parameters of the multilinear tensor model can be approximated by Alternating Least Squares. Further, we introduce a tacked style-separated tensor model, defined as an ensemble of style-specific models to integrate our approach with the extended latent space of StyleGAN. We show that taking the individual styles of the extended latent space into account leads to higher model flexibility and lower reconstruction error. Finally, we do several experiments comparing our approach to former work on both GANs and multilinear models. Concretely, we analyze the expression subspace and find that the expression trajectories meet at an apathetic face that is consistent with earlier work. We also show that by changing the pose of a person, the generated image from our approach is closer to the ground truth than results from two competing approaches. △ Less

Submitted 8 November, 2021; originally announced November 2021.

Comments: Accepted for FG2021

Conference proceedings KI4Industry AI for SMEs -- The online congress for practical entry into AI for SMEs

Authors: Michael Arnemann, Per Olof Beckemeier, Thomas Bertram, Michael Eder, Maximilian Erschig, Matthias Feiner, Francisco Javier Fernandez Garcia, Frederic Foerster, Ruediger Haas, Martin Kipfmueller, Jan Kotschenreuther, Bernd Langer, Ivan Lozada Rodriguez, Thomas Meibert, Simon Ottenhaus, Stefan Paschek, Lars Pfotzer, Michael M. Roth, Tim Schanz, Philip Scherer, Janine Schwienke, Martin Simon, Robin Tenscher-Philipp

Abstract: The Institute of Materials and Processes, IMP, of the University of Applied Sciences in Karlsruhe, Germany in cooperation with VDI Verein Deutscher Ingenieure e.V, AEN Automotive Engineering Network and their cooperation partners present their competences of AI-based solution approaches in the production engineering field. The online congress KI 4 Industry on November 12 and 13, 2020, showed what… ▽ More The Institute of Materials and Processes, IMP, of the University of Applied Sciences in Karlsruhe, Germany in cooperation with VDI Verein Deutscher Ingenieure e.V, AEN Automotive Engineering Network and their cooperation partners present their competences of AI-based solution approaches in the production engineering field. The online congress KI 4 Industry on November 12 and 13, 2020, showed what opportunities the use of artificial intelligence offers for medium-sized manufacturing companies, SMEs, and where potential fields of application lie. The main purpose of KI 4 Industry is to increase the transfer of knowledge, research and technology from universities to small and medium-sized enterprises, to demystify the term AI and to encourage companies to use AI-based solutions in their own value chain or in their products. △ Less

Submitted 5 August, 2021; v1 submitted 14 June, 2021; originally announced June 2021.

Comments: Editors: Matthias Feiner and Manuel Schoellhorn, 72 pages, 48 figures, in German, Conference proceedings KI 4 Industry, 79 pages in total

DataVault: A Data Storage Infrastructure for the Einstein Toolkit

Authors: Yufeng Luo, Roland Haas, Qian Zhang, Gabrielle Allen

Abstract: Data sharing is essential in the numerical simulations research. We introduce a data repository, DataVault, that is designed for data sharing, search and analysis. A comparative study of existing repositories is performed to analyze features that are critical to a data repository. We describe the architecture, workflow, and deployment of DataVault, and provide three use-case scenarios for differen… ▽ More Data sharing is essential in the numerical simulations research. We introduce a data repository, DataVault, that is designed for data sharing, search and analysis. A comparative study of existing repositories is performed to analyze features that are critical to a data repository. We describe the architecture, workflow, and deployment of DataVault, and provide three use-case scenarios for different communities to facilitate the use and application of DataVault. Potential features are proposed and we outline the future development for these features. △ Less

Submitted 15 February, 2021; v1 submitted 11 December, 2020; originally announced December 2020.

Comments: 17 pages, 3 figures, 2 tables

Discriminating Between Similar Nordic Languages

Authors: René Haas, Leon Derczynski

Abstract: Automatic language identification is a challenging problem. Discriminating between closely related languages is especially difficult. This paper presents a machine learning approach for automatic language identification for the Nordic languages, which often suffer miscategorisation by existing state-of-the-art tools. Concretely we will focus on discrimination between six Nordic languages: Danish,… ▽ More Automatic language identification is a challenging problem. Discriminating between closely related languages is especially difficult. This paper presents a machine learning approach for automatic language identification for the Nordic languages, which often suffer miscategorisation by existing state-of-the-art tools. Concretely we will focus on discrimination between six Nordic languages: Danish, Swedish, Norwegian (Nynorsk), Norwegian (Bokmål), Faroese and Icelandic. △ Less

Submitted 23 March, 2023; v1 submitted 11 December, 2020; originally announced December 2020.

Comments: Proceedings of the Eighth Workshop on NLP for Similar Languages, Varieties and Dialects

Deep Learning with Quantized Neural Networks for Gravitational Wave Forecasting of Eccentric Compact Binary Coalescence

Authors: Wei Wei, E. A. Huerta, Mengshen Yun, Nicholas Loutrel, Md Arif Shaikh, Prayush Kumar, Roland Haas, Volodymyr Kindratenko

Abstract: We present the first application of deep learning forecasting for binary neutron stars, neutron star - black hole systems, and binary black hole mergers that span an eccentricity range e <= 0.9. We train neural networks that describe these astrophysical populations, and then test their performance by injecting simulated eccentric signals in advanced LIGO noise available at the \texttt{Gravitationa… ▽ More We present the first application of deep learning forecasting for binary neutron stars, neutron star - black hole systems, and binary black hole mergers that span an eccentricity range e <= 0.9. We train neural networks that describe these astrophysical populations, and then test their performance by injecting simulated eccentric signals in advanced LIGO noise available at the \texttt{Gravitational Wave Open Science Center} to: 1) quantify how fast neural networks identify these signals before the binary components merge; 2) quantify how accurately neural networks estimate the time to merger once gravitational waves are identified; and 3) estimate the time-dependent sky localization of these events from early detection to merger. Our findings show that deep learning can identify eccentric signals from a few seconds (for binary black holes) up to tens of seconds (for binary neutron stars) prior to merger. A quantized version of our neural networks achieves 4x reduction in model size, and up to 2.5x inference speed up. These novel algorithms may be used to facilitate time-sensitive multi-messenger astrophysics observations of compact binaries in dense stellar environments. △ Less

Submitted 17 October, 2021; v1 submitted 7 December, 2020; originally announced December 2020.

Comments: 14 pages, 7 figures

ACM Class: I.2; I.5

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

Initial Data and Eccentricity Reduction Toolkit for Binary Black Hole Numerical Relativity Waveforms

Authors: Sarah Habib, Antoni Ramos-Buades, E. A. Huerta, Sascha Husa, Roland Haas, Zachariah Etienne

Abstract: The production of numerical relativity waveforms that describe quasicircular binary black hole mergers requires high-quality initial data, and an algorithm to iteratively reduce residual eccentricity. To date, these tools remain closed source, or in commercial software that prevents their use in high performance computing platforms. To address these limitations, and to ensure that the broader nume… ▽ More The production of numerical relativity waveforms that describe quasicircular binary black hole mergers requires high-quality initial data, and an algorithm to iteratively reduce residual eccentricity. To date, these tools remain closed source, or in commercial software that prevents their use in high performance computing platforms. To address these limitations, and to ensure that the broader numerical relativity community has access to these tools, herein we provide all the required elements to produce high-quality numerical relativity simulations in supercomputer platforms, namely: open source parameter files to numerical simulate spinning black hole binaries with asymmetric mass-ratios; open source $\texttt{Python}$ tools to produce high-quality initial data for numerical relativity simulations of spinning black hole binaries on quasi-circular orbits; open source $\texttt{Python}$ tools for eccentricity reduction, both as stand-alone software and deployed in the $\texttt{Einstein Toolkit}$'s software infrastructure. This open source toolkit fills in a critical void in the literature at a time when numerical relativity has an ever increasing role in the study and interpretation of gravitational wave sources. As part of our community building efforts, and to streamline and accelerate the use of these resources, we provide tutorials that describe, step by step, how to obtain and use these open source numerical relativity tools. △ Less

Submitted 17 November, 2020; originally announced November 2020.

Comments: 11 pages, 1 figure, 1 table. Open source software available at: https://github.com/NCSAGravity/eccred

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

Source term method for binary neutron stars initial data

Authors: Bing-Jyun Tsao, Roland Haas, Antonios Tsokaros

Abstract: The initial condition problem for a binary neutron star system requires a Poisson equation solver for the velocity potential with a Neumann-like boundary condition on the surface of the star. Difficulties that arise in this boundary value problem are: a) the boundary is not known a-priori, but constitutes part of the solution of the problem; b) various terms become singular at the boundary. In thi… ▽ More The initial condition problem for a binary neutron star system requires a Poisson equation solver for the velocity potential with a Neumann-like boundary condition on the surface of the star. Difficulties that arise in this boundary value problem are: a) the boundary is not known a-priori, but constitutes part of the solution of the problem; b) various terms become singular at the boundary. In this work, we present a new method to solve the fluid Poisson equation for irrotational/spinning binary neutron stars. The advantage of the new method is that it does not require complex fluid surface fitted coordinates and it can be implemented in a Cartesian grid, which is a standard choice in numerical relativity calculations. This is accomplished by employing the source term method proposed by Towers, where the boundary condition is treated as a jump condition and is incorporated as additional source terms in the Poisson equation, which is then solved iteratively. The issue of singular terms caused by vanishing density on the surface is resolved with an additional separation that shifts the computation boundary to the interior of the star. We present two-dimensional tests to show the convergence of the source term method, and we further apply this solver to a realistic three-dimensional binary neutron star problem. By comparing our solution with the one coming from the initial data solver COCAL, we demonstrate agreement to approximately $1\%$. Our method can be used in other problems with non-smooth solutions like in magnetized neutron stars. △ Less

Submitted 20 June, 2021; v1 submitted 17 October, 2020; originally announced October 2020.

Comments: 24 pages, 15 figures

Reconfiguration graphs of zero forcing sets

Authors: Jesse Geneson, Ruth Haas, Leslie Hogben

Abstract: This paper begins the study of reconfiguration of zero forcing sets, and more specifically, the zero forcing graph. Given a base graph $G$, its zero forcing graph, $\mathscr{Z}(G)$, is the graph whose vertices are the minimum zero forcing sets of $G$ with an edge between vertices $B$ and $B'$ of $\mathscr{Z}(G)$ if and only if $B$ can be obtained from $B'$ by changing a single vertex of $G$. It is… ▽ More This paper begins the study of reconfiguration of zero forcing sets, and more specifically, the zero forcing graph. Given a base graph $G$, its zero forcing graph, $\mathscr{Z}(G)$, is the graph whose vertices are the minimum zero forcing sets of $G$ with an edge between vertices $B$ and $B'$ of $\mathscr{Z}(G)$ if and only if $B$ can be obtained from $B'$ by changing a single vertex of $G$. It is shown that the zero forcing graph of a forest is connected, but that many zero forcing graphs are disconnected. We characterize the base graphs whose zero forcing graphs are either a path or the complete graph, and show that the star cannot be a zero forcing graph. We show that computing $\mathscr{Z}(G)$ takes $2^{Θ(n)}$ operations in the worst case for a graph $G$ of order $n$. △ Less

Submitted 1 September, 2020; originally announced September 2020.

MSC Class: 68R10; 05C50; 05C57

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)

Solving the Bethe-Salpeter equation on massively parallel architectures

Authors: Xiao Zhang, Sebastian Achilles, Jan Winkelmann, Roland Haas, André Schleife, Edoardo Di Napoli

Abstract: The last ten years have witnessed fast spreading of massively parallel computing clusters, from leading supercomputing facilities down to the average university computing center. Many companies in the private sector have undergone a similar evolution. In this scenario, the seamless integration of software and middleware libraries is a key ingredient to ensure portability of scientific codes and gu… ▽ More The last ten years have witnessed fast spreading of massively parallel computing clusters, from leading supercomputing facilities down to the average university computing center. Many companies in the private sector have undergone a similar evolution. In this scenario, the seamless integration of software and middleware libraries is a key ingredient to ensure portability of scientific codes and guarantees them an extended lifetime. In this work, we describe the integration of the ChASE library, a modern parallel eigensolver, into an existing legacy code for the first-principles computation of optical properties of materials via solution of the Bethe-Salpeter equation for the optical polarization function. Our numerical tests show that, as a result of integrating ChASE and parallelizing the reading routine, the code experiences a remarkable speedup and greatly improved scaling behavior on both multi- and many-core architectures. We demonstrate that such a modernized BSE code will, by fully exploiting parallel computing architectures and file systems, enable domain scientists to accurately study complex material systems that were not accessible before. △ Less

Submitted 15 June, 2020; originally announced June 2020.

Comments: 17 Pages plus 7 pages of supplemental information, 6 figures and 3 tables. To be submitted to Computer Physics Communications

Multipolar Effective-One-Body Waveforms for Precessing Binary Black Holes: Construction and Validation

Authors: Serguei Ossokine, Alessandra Buonanno, Sylvain Marsat, Roberto Cotesta, Stanislav Babak, Tim Dietrich, Roland Haas, Ian Hinder, Harald P. Pfeiffer, Michael Pürrer, Charles J. Woodford, Michael Boyle, Lawrence E. Kidder, Mark A. Scheel, Béla Szilágyi

Abstract: As gravitational-wave detectors become more sensitive, we will access a greater variety of signals emitted by compact binary systems, shedding light on their astrophysical origin and environment. A key physical effect that can distinguish among formation scenarios is the misalignment of the spins with the orbital angular momentum, causing the spins and the binary's orbital plane to precess. To acc… ▽ More As gravitational-wave detectors become more sensitive, we will access a greater variety of signals emitted by compact binary systems, shedding light on their astrophysical origin and environment. A key physical effect that can distinguish among formation scenarios is the misalignment of the spins with the orbital angular momentum, causing the spins and the binary's orbital plane to precess. To accurately model such systems, it is crucial to include multipoles beyond the dominant quadrupole. Here, we develop the first multipolar precessing waveform model in the effective-one-body (EOB) formalism for the inspiral, merger and ringdown (IMR) of binary black holes: SEOBNRv4PHM. In the nonprecessing limit, the model reduces to SEOBNRv4HM, which was calibrated to numerical-relativity (NR) simulations, and waveforms from perturbation theory. We validate SEOBNRv4PHM by comparing it to the public catalog of 1405 precessing NR waveforms of the Simulating eXtreme Spacetimes (SXS) collaboration, and also to new 118 precessing NR waveforms, which span mass ratios 1-4 and spins up to 0.9. We stress that SEOBNRv4PHM is not calibrated to NR simulations in the precessing sector. We compute the unfaithfulness against the 1523 SXS precessing NR waveforms, and find that, for $94\%$ ($57\%$) of the cases, the maximum value, in the total mass range $20-200 M_\odot$, is below $3\%$ ($1\%$). Those numbers become $83\%$ ($20\%$) when using the IMR, multipolar, precessing phenomenological model IMRPhenomPv3HM. We investigate the impact of such unfaithfulness values with two parameter-estimation studies on synthetic signals. We also compute the unfaithfulness between those waveform models and identify in which part of the parameter space they differ the most. We validate them also against the multipolar, precessing NR surrogate model NRSur7dq4, and find that the SEOBNRv4PHM model outperforms IMRPhenomPv3HM. △ Less

Submitted 20 April, 2020; originally announced April 2020.

Comments: 24 pages, 18 figures. Abstract abridged

Journal ref: Phys. Rev. D 102, 044055 (2020)

Alpha spectrometric characterization of thin $^{233}$U sources for $^{229\text{(m)}}$Th production

Authors: Raphael Haas, Michelle Hufnagel, Roman Abrosimov, Christoph E. Düllmann, Dominik Krupp, Christoph Mokry, Dennis Renisch, Jörg Runke, Ulrich W. Scherer

Abstract: Four different techniques were applied for the production of $^{233}$U alpha recoil ion sources, providing $^{229}$Th ions. They were compared with respect to a minimum energy spread of the $^{229}$Th recoil ions, using the emitted alpha particles as an indicator. The techniques of Molecular Plating, Drop-on-Demand inkjet printing, chelation from dilute nitric acid solution on chemically functiona… ▽ More Four different techniques were applied for the production of $^{233}$U alpha recoil ion sources, providing $^{229}$Th ions. They were compared with respect to a minimum energy spread of the $^{229}$Th recoil ions, using the emitted alpha particles as an indicator. The techniques of Molecular Plating, Drop-on-Demand inkjet printing, chelation from dilute nitric acid solution on chemically functionalized silicon surfaces, and self-adsorption on passivated titanium surfaces were used. All fabricated sources were characterized by using alpha spectrometry, radiographic imaging, and scanning electron microscopy. A direct validation for the estimated recoil ion rate was obtained by collecting $^{228}$Th recoil ions from $^{232}$U recoil ion sources prepared by self-adsorption and Molecular Plating. The chelation and the self-adsorption based approaches appear most promising for the preparation of recoil ion sources delivering monochromatic recoil ions. △ Less

Submitted 16 November, 2020; v1 submitted 6 April, 2020; originally announced April 2020.

Comments: 20 pages, 14 figures, 11 tables, ahead of print at Radiochim. Acta (2020)

Journal ref: Radiochimica Acta (2020) ahead of print

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)

Development of a recoil ion source providing slow Th ions including $^{229(m)}$Th in a broad charge state distribution

Authors: Raphael Haas, Tom Kieck, Dmitry Budker, Christoph E. Düllmann, Karin Groot-Berning, Wenbing Li, Dennis Renisch, Ferdinand Schmidt-Kaler, Felix Stopp, Anna Viatkina

Abstract: Ions of the isomer $^{229m}$Th are a topic of high interest for the construction of a "nuclear clock" and in the field of fundamental physics for testing symmetries of nature. They can be efficiently captured in Paul traps which are ideal for performing high precision quantum logic spectroscopy. Trapping and identification of long-lived $^{232}$Th$^{+}$ ions from a laser ablation source was alread… ▽ More Ions of the isomer $^{229m}$Th are a topic of high interest for the construction of a "nuclear clock" and in the field of fundamental physics for testing symmetries of nature. They can be efficiently captured in Paul traps which are ideal for performing high precision quantum logic spectroscopy. Trapping and identification of long-lived $^{232}$Th$^{+}$ ions from a laser ablation source was already demonstrated by the TACTICa collaboration on Trapping And Cooling of Thorium Ions with Calcium. The $^{229m}$Th is most easily accessible as $α$-decay daughter of the decay of $^{233}$U. We report on the development of a source for slow Th ions, including $^{229(m)}$Th for the TACTICa experiment. The $^{229(m)}$Th source is currently under construction and comprises a $^{233}$U monolayer, from which $^{229(m)}$Th ions recoil. These are decelerated in an electric field. Conservation of the full initial charge state distribution of the $^{229(m)}$Th recoil ions is one of the unique features of this source. We present ion-flight simulations for our adopted layout and give a final design. This source will provide Th ions in their original charge state at energies suitable for capture in a linear Paul trap for spectroscopy investigations. △ Less

Submitted 12 February, 2020; v1 submitted 26 November, 2019; originally announced November 2019.

Comments: 6 pages, 3 figures, PLATAN19 conference proceeding published in Hyperfine Interact 2020

Journal ref: Hyperfine Interact 241, 25 (2020)

Unequal Mass Binary Neutron Star Simulations with Neutrino Transport: Ejecta and Neutrino Emission

Authors: Trevor Vincent, Francois Foucart, Matthew D. Duez, Roland Haas, Lawrence E. Kidder, Harald P. Pfeiffer, Mark A. Scheel

Abstract: We present twelve new simulations of unequal mass neutron star mergers. The simulations were preformed with the SpEC code, and utilize nuclear-theory based equations of state and a two-moment gray neutrino transport scheme with an improved energy estimate based on evolving the number density. We model the neutron stars with the SFHo, LS220 and DD2 equations of state (EOS) and we study the neutrino… ▽ More We present twelve new simulations of unequal mass neutron star mergers. The simulations were preformed with the SpEC code, and utilize nuclear-theory based equations of state and a two-moment gray neutrino transport scheme with an improved energy estimate based on evolving the number density. We model the neutron stars with the SFHo, LS220 and DD2 equations of state (EOS) and we study the neutrino and matter emission of all twelve models to search for robust trends between binary parameters and emission characteristics. We find that the total mass of the dynamical ejecta exceeds $0.01M_\odot$ only for SFHo with weak dependence on the mass-ratio across all models. We find that the ejecta have a broad electron fraction ($Y_e$) distribution ($\approx 0.06-0.48$), with mean $0.2$. $Y_e$ increases with neutrino irradiation over time, but decreases with increasing binary asymmetry. We also find that the models have ejecta with a broad asymptotic velocity distribution ($\approx 0.05-0.7c$). The average velocity lies in the range $0.2c - 0.3c$ and decreases with binary asymmetry. Furthermore, we find that disk mass increases with binary asymmetry and stiffness of the EOS. The $Y_e$ of the disk increases with softness of the EOS. The strongest neutrino emission occurs for the models with soft EOS. For (anti) electron neutrinos we find no significant dependence of the magnitude or angular distribution or neutrino luminosity with mass-ratio. The heavier neutrino species have a luminosity dependence on mass-ratio but an angular distribution which does not change with mass-ratio. △ Less

Submitted 1 August, 2019; originally announced August 2019.

Comments: 23 pages, 15 figures

MSC Class: 83-08

Journal ref: Phys. Rev. D 101, 044053 (2020)

Cographs

Authors: Robert Haas

Abstract: Cographs--defined most simply as complete graphs with colored lines--both dualize and generalize ordinary graphs, and promise a comparably wide range of applications. This article introduces them by examples, catalogues, and elementary properties. Any finite cograph may be realized in several ways, including inner products, polynomials, geometrically, or by "fat intersections." Particular classes… ▽ More Cographs--defined most simply as complete graphs with colored lines--both dualize and generalize ordinary graphs, and promise a comparably wide range of applications. This article introduces them by examples, catalogues, and elementary properties. Any finite cograph may be realized in several ways, including inner products, polynomials, geometrically, or by "fat intersections." Particular classes then considered include sum cographs (points in Z or Zn; the line C(P,Q) joining points P and Q defined C(P,Q) = P+Q); difference cographs (C(P,Q) = |P-Q|; and intersection cographs (points are sets; C(P,Q) = P intersect Q). Intersection cographs, especially, promise many applications; described here are some to aesthetics. Point-line cographs turn out equivalent to linear spaces. Finally solved here is an interesting group-theoretic problem arising from group cographs (points in a group; C(P,Q) = {PQ,QP}). △ Less

Submitted 28 May, 2019; originally announced May 2019.

Catching, trapping and in-situ-identification of thorium ions inside Coulomb crystals of $^{40}$Ca$^+$ ions

Authors: Felix Stopp, Karin Groot-Berning, Georg Jacob, Dmitry Budker, Raphael Haas, Dennis Renisch, Jörg Runke, Petra Thörle-Pospiech, Christoph E. Düllmann, Ferdinand Schmidt-Kaler

Abstract: Thorium ions exhibit unique nuclear properties with high relevance for testing symmetries of nature, and Paul traps feature an ideal experimental platform for performing high precision quantum logic spectroscopy. Loading of stable or long-lived isotopes is well-established and relies on ionization from an atomic beam. A different approach allows trapping short-lived isotopes available as alpha-dec… ▽ More Thorium ions exhibit unique nuclear properties with high relevance for testing symmetries of nature, and Paul traps feature an ideal experimental platform for performing high precision quantum logic spectroscopy. Loading of stable or long-lived isotopes is well-established and relies on ionization from an atomic beam. A different approach allows trapping short-lived isotopes available as alpha-decay daughters, which recoil from a thin sample of the precursor nuclide. A prominent example is the short-lived $^{229\text{m}}$Th, populated in a decay of long-lived $^{233}$U. Here, ions are provided by an external source and are decelerated to be available for trapping. Such setups offer the option to trap various isotopes and charge states of thorium. Investigating this complex procedure, we demonstrate the observation of single $^{232}$Th$^+$ ions trapped, embedded into and sympathetically cooled via Coulomb interactions by co-trapped $^{40}$Ca$^+$ ions. Furthermore, we discuss different options for a non-destructive identification of the sympathetically cooled thorium ions in the trap, and describe in detail our chosen experimental method, identifying mass and charge of thorium ions from the positions of calcium ions, as their fluorescence is imaged on a CCD camera. These findings are verified by means of a time-of-flight signal when extracting ions of different mass-to-charge ratio from the Paul trap and steering them into a detector. △ Less

Submitted 23 April, 2019; originally announced April 2019.

Journal ref: Hyperfine Interact (2019) 240: 33

The concept of laser-based conversion electron Mössbauer spectroscopy for a precise energy determination of $^{229m}$Th

Authors: Lars C. von der Wense, Benedict Seiferle, Christian Schneider, Justin Jeet, Ines Amersdorffer, Nicolas Arlt, Florian Zacherl, Raphael Haas, Dennis Renisch, Patrick Mosel, Philip Mosel, Milutin Kovacev, Uwe Morgner, Christoph E. Düllmann, Eric R. Hudson, Peter G. Thirolf

Abstract: $^{229}$Th is the only nucleus currently under investigation for the development of a nuclear optical clock (NOC) of ultra-high accuracy. The insufficient knowledge of the first nuclear excitation energy of $^{229}… ▽ More $^{229}$Th is the only nucleus currently under investigation for the development of a nuclear optical clock (NOC) of ultra-high accuracy. The insufficient knowledge of the first nuclear excitation energy of $^{229}$Th has so far hindered direct nuclear laser spectroscopy of thorium ions and thus the development of a NOC. Here, a nuclear laser excitation scheme is detailed, which makes use of thorium atoms instead of ions. This concept, besides potentially leading to the first nuclear laser spectroscopy, would determine the isomeric energy to 40 $μ$eV resolution, corresponding to 10 GHz, which is a $10^4$ times improvement compared to the current best energy constraint. This would determine the nuclear isomeric energy to a sufficient accuracy to allow for nuclear laser spectroscopy of individual thorium ions in a Paul trap and thus the development of a single-ion nuclear optical clock. △ Less

Submitted 2 April, 2019; originally announced April 2019.

Journal ref: Hyperfine Interact (2019) 240: 23

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/

Time-symmetry, symplecticity and stability of Euler-Maclaurin and Lanczos-Dyche integration

Authors: Charalampos M. Markakis, Michael F. O'Boyle, Derek Glennon, Khoa Tran, Pablo Brubeck, Roland Haas, Hsi-Yu Schive, Kōji Uryū

Abstract: Numerical evolution of time-dependent differential equations via explicit Runge-Kutta or Taylor methods typically fails to preserve symmetries of a system. It is known that there exists no numerical integration method that in general preserves both the energy and the symplectic structure of a Hamiltonian system. One is thus normally forced to make a choice. Nevertheless, a symmetric integration fo… ▽ More Numerical evolution of time-dependent differential equations via explicit Runge-Kutta or Taylor methods typically fails to preserve symmetries of a system. It is known that there exists no numerical integration method that in general preserves both the energy and the symplectic structure of a Hamiltonian system. One is thus normally forced to make a choice. Nevertheless, a symmetric integration formula, obtained by Lanczos-Dyche via two-point Taylor expansion (or Hermite interpolation), is shown here to preserve both energy as well as symplectic structure for linear systems. This formula shares similarities with the Euler-Maclaurin formula, but is superconvergent rather than asymptotically convergent. For partial differential equations, the resulting evolution methods are unconditionally stable, i.e, not subject to a Courant-Friedrichs-Lewy limit. Although generally implicit, these methods become explicit for linear systems. △ Less

Submitted 28 January, 2019; originally announced January 2019.

Comments: Preprint submitted to J. Comp. Phys

MSC Class: 65P10 (Primary) 65L20; 65M20 (Secondary)