The spin distribution of binary black hole mergers contains key information concerning the formation channels of these objects, and the astrophysical environments where they form, evolve and coalesce. To quantify the suitability of deep learning to characterize the signal manifold of quasi-circular, spinning, non-precessing binary black hole mergers, we introduce a modified version of WaveNet trained with a novel optimization scheme that incorporates general relativistic constraints of the spin properties of astrophysical black holes. The neural network model is trained, validated and tested with 1.5 million $\ell=|m|=2$ waveforms generated within the regime of validity of NRHybSur3dq8, i.e., mass-ratios $q\leq8$ and individual black hole spins $ | s^z_{\{1,\,2\}} | \leq 0.8$. Using this neural network model, we quantify how accurately we can infer the astrophysical parameters of black hole mergers in the absence of noise. We do this by computing the overlap between waveforms in the testing data set and the corresponding signals whose mass-ratio and individual spins are predicted by our neural network. We find that the convergence of high performance computing and physics-inspired optimization algorithms enable an accurate reconstruction of the mass-ratio and individual spins of binary black hole mergers across the parameter space under consideration. This is a significant step towards an informed utilization of physics-inspired deep learning models to reconstruct the spin distribution of binary black hole mergers in realistic detection scenarios.

Significant investments to upgrade or construct large-scale scientific facilities demand commensurate investments in R&D to design algorithms and computing approaches to enable scientific and engineering breakthroughs in the big data era. The remarkable success of Artificial Intelligence (AI) algorithms to turn big-data challenges in industry and technology into transformational digital solutions that drive a multi-billion dollar industry, which play an ever increasing role shaping human social patterns, has promoted AI as the most sought after signal processing tool in big-data research. As AI continues to evolve into a computing tool endowed with statistical and mathematical rigor, and which encodes domain expertise to inform and inspire AI architectures and optimization algorithms, it has become apparent that single-GPU solutions for training, validation, and testing are no longer sufficient. This realization has been driving the confluence of AI and high performance computing (HPC) to reduce time-to-insight and to produce robust, reliable, trustworthy, and computationally efficient AI solutions. In this white paper, we present a summary of recent developments in this field, and discuss avenues to accelerate and streamline the use of HPC platforms to design accelerated AI algorithms.

A recent new measurement and re-analysis of past measurements suggested an improved quenching factor value and uncertainty for CsI[Na]. This implies a measurement of the COHERENT experiment of coherent elastic neutrino-nucleus scattering that is closer to the Standard Model prediction and has less uncertainty. We illustrate the impact of this improvement by revisiting fits to the Weinberg angle, neutrino magnetic moments, neutron rms and neutrino charge radii, neutrino non-standard interactions (in particular those relevant for LMA-Dark) and new scalar as well as vector bosons. Significant improvement is observed, particularly for those scenarios coherently affecting the electroweak SM process.

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.

The precise determination of the physical and dynamical parameters of the HIP 109951 triple star system (WDS J22161-0705AB) which is formed by the A, Ba, and Bb components are presented. The binary nature of component B was recently confirmed by studying the radial velocities. The analysis of the system follows Al-Wardat's complex method for analyzing CVBS which employs Kurucz (Atlas9) line blanketed plane-parallel atmospheres simultaneously with an analytic method for dynamical analysis (we used Docobo's method) to calculate the parameters of this triple system. The result of our study yielded the following parameters: $T^{A}_{\rm eff}=5836\pm 80$ K, $R^{A}=1.090\pm0.039R_\odot$, $\log$g$^{A}=4.45\pm0.06$, $\mathcal{M}^{A}=1.05\pm0.16\mathcal{M}_\odot$; $T^{Ba}_{\rm eff}=5115\pm 80$ K, $R^{Ba}=0.596\pm0.05R_\odot$, $\log$ g$^{Ba}=4.60\pm0.07$, $\mathcal{M}^{Ba}=0.83\pm0.16\mathcal{M}_\odot$, $T^{Bb}_{\rm eff}=4500\pm 80$ K, $R^{Bb}=0.490\pm0.06R_\odot$, $\log$ g$^{Bb}=4.65\pm0.07$, and $\mathcal{M}^{Bb}=0.67\pm0.16\mathcal{M}_\odot$ based on the revised \textit{Hipparcos} parallax. The orbital solution gave a total mass as $\mathcal{M}=2.59\mathcal{M}_\odot$ based on Gaia parallax and $\mathcal{M}=2.15\mathcal{M}_\odot$ based on the revised \textit{Hipparcos} parallax. The synthetic spectral energy distributions (SED) and synthetic stellar photometry of the entire system and individual components are given and compared with the available observational ones. Finally, the positions of the system components on the HR diagram and the evolutionary tracks are given and their formation and the evolution of the system are discussed.

The scale of ongoing and future electromagnetic surveys pose formidable challenges to classify astronomical objects. Pioneering efforts on this front include citizen science campaigns adopted by the Sloan Digital Sky Survey (SDSS). SDSS datasets have been recently used to train neural network models to classify galaxies in the Dark Energy Survey (DES) that overlap the footprint of both surveys. Herein, we demonstrate that knowledge from deep learning algorithms, pre-trained with real-object images, can be transferred to classify galaxies that overlap both SDSS and DES surveys, achieving state-of-the-art accuracy $\gtrsim99.6\%$. We demonstrate that this process can be completed within just eight minutes using distributed training. While this represents a significant step towards the classification of DES galaxies that overlap previous surveys, we need to initiate the characterization of unlabelled DES galaxies in new regions of parameter space. To accelerate this program, we use our neural network classifier to label over ten thousand unlabelled DES galaxies, which do not overlap previous surveys. Furthermore, we use our neural network model as a feature extractor for unsupervised clustering and find that unlabeled DES images can be grouped together in two distinct galaxy classes based on their morphology, which provides a heuristic check that the learning is successfully transferred to the classification of unlabelled DES images. We conclude by showing that these newly labeled datasets can be combined with unsupervised recursive training to create large-scale DES galaxy catalogs in preparation for the Large Synoptic Survey Telescope era.

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.

We perform magnetohydrodynamic simulations in full general relativity of disk accretion onto nonspinning black hole binaries with mass ratio 36:29. We survey different disk models which differ in their scale height, total size and magnetic field to quantify the robustness of previous simulations on the initial disk model. Scaling our simulations to LIGO GW150914 we find that such systems could explain possible gravitational wave and electromagnetic counterparts such as the Fermi GBM hard X-ray signal reported 0.4s after GW150915 ended. Scaling our simulations to supermassive binary black holes, we find that observable flow properties such as accretion rate periodicities, the emergence of jets throughout inspiral, merger and post-merger, disk temperatures, thermal frequencies, and the time-delay between merger and the boost in jet outflows that we reported in earlier studies display only modest dependence on the initial disk model we consider here.

We describe the construction and characterization of the 280 GHz bolometric focal plane units (FPUs) to be deployed on the second flight of the balloon-borne SPIDER instrument. These FPUs are vital to SPIDER's primary science goal of detecting or placing an upper limit on the amplitude of the primordial gravitational wave signature in the cosmic microwave background (CMB) by constraining the B-mode contamination in the CMB from Galactic dust emission. Each 280 GHz focal plane contains a 16 x 16 grid of corrugated silicon feedhorns coupled to an array of aluminum-manganese transition-edge sensor (TES) bolometers fabricated on 150 mm diameter substrates. In total, the three 280 GHz FPUs contain 1,530 polarization sensitive bolometers (765 spatial pixels) optimized for the low loading environment in flight and read out by time-division SQUID multiplexing. In this paper we describe the mechanical, thermal, and magnetic shielding architecture of the focal planes and present cryogenic measurements which characterize yield and the uniformity of several bolometer parameters. The assembled FPUs have high yields, with one array as high as 95% including defects from wiring and readout. We demonstrate high uniformity in device parameters, finding the median saturation power for each TES array to be ~3 pW at 300 mK with a less than 6% variation across each array at one standard deviation. These focal planes will be deployed alongside the 95 and 150 GHz telescopes in the SPIDER-2 instrument, slated to fly from McMurdo Station in Antarctica in December 2018.

The Outer Scutum-Centaurus arm (OSC) is the most distant molecular spiral arm known in the Milky Way. The OSC may be the very distant end of the well-known Scutum-Centaurus arm, which stretches from the end of the Galactic bar to the outer Galaxy. At this distance the OSC is seen in the first Galactic quadrant. The population of star formation tracers in the OSC remains largely uncharacterized. Extragalactic studies show a strong correlation between molecular gas and star formation, and carbon monoxide (CO) emission was recently discovered in the OSC. Here we use the Arizona Radio Observatory (ARO) 12-m telescope to observe the $^{12}$CO J = 1-0 and $^{13}$CO J = 1-0 transitions toward 78 HII region candidates chosen from the WISE Catalog of Galactic HII Regions. These targets are spatially coincident with the Galactic longitude-latitude ($\ell, b$) OSC locus as defined by HI emission. We detect CO emission in $\sim 80$% of our targets. In total, we detect 117 $^{12}$CO and 40 $^{13}$CO emission lines. About 2/3 of our targets have at least one emission line originating beyond the Solar orbit. Most of the detections beyond the Solar orbit are associated with the Outer Arm, but there are 17 $^{12}$CO emission lines and 8 $^{13}$CO emission lines with LSR velocities that are consistent with the velocities of the OSC. There is no apparent difference between the physical properties (e.g., molecular column density) of these OSC molecular clouds and non--OSC molecular clouds within our sample.

We aim to present a generalized Bayesian inference method for constraining interiors of super Earths and sub-Neptunes. Our methodology succeeds in quantifying the degeneracy and correlation of structural parameters for high dimensional parameter spaces. Specifically, we identify what constraints can be placed on composition and thickness of core, mantle, ice, ocean, and atmospheric layers given observations of mass, radius, and bulk refractory abundance constraints (Fe, Mg, Si) from observations of the host star's photospheric composition. We employed a full probabilistic Bayesian inference analysis that formally accounts for observational and model uncertainties. Using a Markov chain Monte Carlo technique, we computed joint and marginal posterior probability distributions for all structural parameters of interest. We included state-of-the-art structural models based on self-consistent thermodynamics of core, mantle, high-pressure ice, and liquid water. Furthermore, we tested and compared two different atmospheric models that are tailored for modeling thick and thin atmospheres, respectively. First, we validate our method against Neptune. Second, we apply it to synthetic exoplanets of fixed mass and determine the effect on interior structure and composition when (1) radius, (2) atmospheric model, (3) data uncertainties, (4) semi-major axes, (5) atmospheric composition (i.e., a priori assumption of enriched envelopes versus pure H/He envelopes), and (6) prior distributions are varied. Our main conclusions are: [...]

We present an inversion method based on Bayesian analysis to constrain the interior structure of terrestrial exoplanets, in the form of chemical composition of the mantle and core size. Specifically, we identify what parts of the interior structure of terrestrial exoplanets can be determined from observations of mass, radius, and stellar elemental abundances. We perform a full probabilistic inverse analysis to formally account for observational and model uncertainties and obtain confidence regions of interior structure models. This enables us to characterize how model variability depends on data and associated uncertainties. We test our method on terrestrial solar system planets and find that our model predictions are consistent with independent estimates. Furthermore, we apply our method to synthetic exoplanets up to 10 Earth masses and up to 1.7 Earth radii as well as to exoplanet Kepler-36b. Importantly, the inversion strategy proposed here provides a framework for understanding the level of precision required to characterize the interior of exoplanets. Our main conclusions are: (1) observations of mass and radius are sufficient to constrain core size; (2) stellar elemental abundances (Fe, Si, Mg) are key constraints to reduce degeneracy in interior structure models and to constrain mantle composition; (3) the inherent degeneracy in determining interior structure from mass and radius observations does not only depend on measurement accuracies but also on the actual size and density of the exoplanet. We argue that precise observations of stellar elemental abundances are central in order to place constraints on planetary bulk composition and to reduce model degeneracy. [...]