We propose a method based on the process of extracting gravitational wave (GW) parameters from GW signals to simulate the binary neutron-star (BNS) merging events. We simulate 1000 GW standard sirens based on the observation of the Einstein Telescope (ET). Almost all the simulated GW data are in the redshift range of $[0,3]$. The role of the GW standard siren in the inference of the cosmological parameters is investigated. We find that the GW data can help improve the accuracy of cosmological parameters. Moreover, the degeneracy of cosmological parameters is broken by the GW data. The GW standard siren is helpful for the constraint of the cosmological parameters.

TianQin is a space-based gravitational-wave observatory scheduled to be launched in the 2030s. In this work, we make a preliminary forecast for the cosmological parameter estimation with the gravitational-wave standard siren observation from TianQin. We simulate the standard siren data of TianQin based on its 5-year observation after the completion of construction. In the simulation, three models for the population of massive black hole binary (MBHB), i.e., pop III, Q3nod, and Q3d, are considered to predict the event numbers of MBHB mergers. We find that: (i) among the three MBHB models, the Q3nod model can provide the tightest constraints on the cosmological parameters; (ii) TianQin's standard siren observation can effectively break the parameter degeneracies inherent in the cosmic microwave background observation; and (iii) the future standard siren observation from TianQin can significantly improve the cosmological parameter estimation under the current mainstream electromagnetic observations.

Fast radio bursts (FRBs) are a mysterious astrophysical phenomenon of bright pulses emitted at radio frequencies, and it is expected for them to be frequently detected in the future. The dispersion measures of FRBs are related to cosmological parameters, and thus FRBs have a potential to be developed into a new cosmological probe if their data can be largely accumulated in the future. In this work, we study the capability of future FRB data for improving cosmological parameter estimation in two dynamical dark energy models. We find that the simulated FRB data can break the parameter degeneracies inherent in the current cosmic microwave background (CMB) data. Therefore, the combination of the CMB and FRB data can significantly improve the constraints on the Hubble constant and dark energy parameters, compared to those using CMB or FRB alone. If ten thousand events of FRBs with known redshifts are detected in the future, they would behave better than the baryon acoustic oscillation (BAO) data in breaking the parameter degeneracies inherent in the CMB data. We also find that the combination of FRB and gravitational wave (GW) standard siren data provides an independent low-redshift probe to verify the results from CMB and BAO data. For the data combination of CMB, GW, and FRB, it is found that the main contribution to the constraints comes from the CMB and GW data, but the inclusion of FRB data still can evidently improve the constraint on the baryon density.

Taiji, a space-based gravitational-wave observatory, consists of three satellites forming an equilateral triangle with arm length of $3\times 10^6$ km, orbiting around the Sun. Taiji is able to observe the gravitational-wave standard siren events of massive black hole binary (MBHB) merger, which is helpful in probing the expansion of the universe. In this paper, we preliminarily forecast the capability of Taiji for improving cosmological parameter estimation with the gravitational-wave standard siren data. We simulate five-year standard siren data based on three fiducial cosmological models and three models of MBHB's formation and growth. It is found that the standard siren data from Taiji can effectively break the cosmological parameter degeneracies generated by the cosmic microwave background (CMB) anisotropies data, especially for dynamical dark energy models. The constraints on cosmological parameters are significantly improved by the data combination CMB+Taiji, compared to the CMB data alone. Compared to the current optical cosmological observations, Taiji can still provide help in improving the cosmological parameter estimation to some extent. In addition, we consider an ideal scenario to investigate the potential of Taiji on constraining cosmological parameters. We conclude that the standard sirens of MBHB from Taiji will become a powerful cosmological probe in the future.

Because the binary neutron star (BNS) merger GW170817 and its electromagnetic (EM) counterparts were successfully detected, there is a chance to explore the joint effect of the host galaxy and our Milky Way on the weak equivalence principle (WEP) test. In this paper, through using the Navarro-Frenk-White profile and the Herquist profile, we present a new analytic model to calculate the galactic gravitational potential, in which the possible locations of the source by the observed angle offset and the second SN (SN2) kick is accounted for. We show that the upper limit of $\Delta \gamma$ is $10^{-9}$ for the comparison of GW170817/GRB 170817A, and is $10^{-4}$ for the comparison of GW170817/AT2017gfo, which are tighter by one or two orders of magnitude than these of only by the measured Milky Way gravitational potential in literatures. We demonstrate that the WEP test could be strengthened by the contribution of the host galaxy NGC 4993 to the Shapiro time delay. Meanwhile, we also find the large natal kick can produce a maximum deviation of about $20\%$ on the results with a typical kick velocity $400\sim 500 km s^{-1}$. Finally, we study and analyse the impact from the halo mass of NGC 4993 with a typical 0.2 dex uncertainty, and find the upper limit of $\Delta \gamma$ of the maximum mass $10^{12.4}h^{-1} M_{\odot}$ is nearly two times tighter than that of the minimum mass $10^{12.0}h^{-1} M_{\odot}$. It suggests the constraint on WEP test would be also enhanced when the halo mass increases in the host galaxy NGC 4993.

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.

We present the first application of deep learning for gravitational wave parameter estimation of binary black hole mergers evolving on quasi-circular orbits with aligned or anti-aligned spins. We use root-leaf structured networks to ensure that common physical features are shared across all parameters. In order to cover a broad range of astrophysically motivated scenarios, we use a training dataset with over $10^7$ modeled waveforms to ensure local time- and scale-invariance. The trained models are applied to estimate the astrophysical parameters of the existing catalog of detected binary black hole mergers, and their corresponding black hole remnants, including the final spin and the gravitational wave quasi-normal frequencies. Using a deterministic neural network model, we are able to efficiently provide point-parameter estimation results, along with statistical errors caused by the noise spectrum uncertainty. We also introduce the first application of Bayesian neural networks for gravitational wave parameter estimation of real astrophysical events. These probabilistic models were trained with over $10^7$ modeled waveforms and using 1024 nodes (65,536 core processors) on the Theta supercomputer at Argonne Leadership Computing Facility to reduce the training stage to just thirty minutes. In inference mode, both the deterministic and Bayesian neural networks estimate the astrophysical parameters of binary black hole mergers within 2 milliseconds using a single Tesla V100 GPU. Both deterministic and Bayesian neural networks produce agreeing parameter estimation results, which are also consistent with Bayesian analyses used to characterize the catalog of binary black hole mergers observed by the advanced LIGO and Virgo detectors.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

Denoising of time domain data is a crucial task for many applications such as communication, translation, virtual assistants etc. For this task, a combination of a recurrent neural net (RNNs) with a Denoising Auto-Encoder (DAEs) has shown promising results. However, this combined model is challenged when operating with low signal-to-noise ratio (SNR) data embedded in non-Gaussian and non-stationary noise. To address this issue, we design a novel model, referred to as 'Enhanced Deep Recurrent Denoising Auto-Encoder' (EDRDAE), that incorporates a signal amplifier layer, and applies curriculum learning by first denoising high SNR signals, before gradually decreasing the SNR until the signals become noise dominated. We showcase the performance of EDRDAE using time-series data that describes gravitational waves embedded in very noisy backgrounds. In addition, we show that EDRDAE can accurately denoise signals whose topology is significantly more complex than those used for training, demonstrating that our model generalizes to new classes of gravitational waves that are beyond the scope of established denoising algorithms.

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.

A space-based very long baseline interferometry (VLBI) programme, named as the Cosmic Microscope, is proposed to involve dual VLBI telescopes in the space working together with giant ground-based telescopes (e.g., Square Kilometre Array, FAST, Arecibo) to image the low radio frequency Universe with the purpose of unraveling the compact structure of cosmic constituents including supermassive black holes and binaries, pulsars, astronomical masers and the underlying source, and exoplanets amongst others. The operational frequency bands are 30, 74, 330 and 1670 MHz, supporting broad science areas. The mission plans to launch two 30-m-diameter radio telescopes into 2,000 km x 90,000 km elliptical orbits. The two telescopes can work in flexibly diverse modes: (i) space-ground VLBI. The maximum space-ground baseline length is about 100,000 km; it provides a high-dynamic-range imaging capacity with unprecedented high resolutions at low frequencies (0.4 mas at 1.67 GHz and 20 mas at 30 MHz) enabling studies of exoplanets and supermassive black hole binaries (which emit nanoHz gravitational waves); (ii) space-space single-baseline VLBI. This unique baseline enables the detection of flaring hydroxyl masers, and more precise position measurement of pulsars and radio transients at milli-arcsecond level; (iii) single dish mode, where each telescope can be used to monitor transient bursts and rapidly trigger follow-up VLBI observations. The large space telescope will also contribute in measuring and constraining the total angular power spectrum from the Epoch of Reionization. In short, the Cosmic Microscope offers astronomers the opportunity to conduct novel, frontier science.

The composition of the astrophysical relativistic jets remains uncertain. By kinetic particle-in-cell simulations, we show that the baryon component in the jet, or the so-called baryon loading effect (BLE), heavily affects relativistic jets transport dynamics in the interstellar medium. On the one hand, with the BLE, relativistic jets can transport in a much longer distance, because jet electrons draw a significant amount of energy from jet baryons via the Buneman-induced electrostatic waves and the Weibel-mediated collisionless shock; on the other hand, the jet electron phase space distribution may transform from a bottom-wide-single-peak structure to a center-wide-multiple-peak one by increasing the BLE, which largely influences the observed jet morphology. Implications for related astrophysical studies are also discussed.

In this paper we introduce a software package, the Very Long Baseline Interferometry Network SIMulator (VNSIM), which provides an integrated platform assisting radio astronomers to design the Very Long Baseline Interferometry (VLBI) experiments and evaluate the network performance with a user-friendly interface. Though VNSIM is motivated to be designed for the East Asia VLBI Network, it can also be expandable to other VLBI networks and generic interferometers. The software package not only integrates the functionality of plotting $(u,v)$ coverage, scheduling the observation, and displaying the dirty and CLEAN images, but also further extends to add new features such as the sensitivity calculation of a certain network and multiple-satellite space VLBI simulations which are useful for future space VLBI mission. In addition, VNSIM provides flexible interactions on both command line and graphical user interface and offers friendly support for log report and database management. VNSIM also supports multiprocessing accelerations, enabling to handle large survey data. To facilitate the future development and update, each simulation function is encapsulated in different Python module allowing for independently invoking and testing. In order to verify the performance of VNSIM, we have carried out various simulations and compared the results with other simulation tools. All tests show good consistency.

We report a possible dipole anisotropy on acceleration scale $g_{\dag}$ with 147 rotationally supported galaxies in local Universe. It is found that a monopole and dipole correction for the radial acceleration relation can better describe the SPARC data set. The monopole term is negligible but the dipole magnitude is significant. It is also found that the dipole correction is mostly induced by the anisotropy on the acceleration scale. The magnitude of $\hat{g}_{\dag}$-dipole reaches up to $0.25\pm0.04$, and its direction is aligned to $(l,b) = (171.30^{\circ}\pm7.18^{\circ}, -15.41^{\circ}\pm4.87^{\circ})$, which is very close to the maximum anisotropy direction from the hemisphere comparison method. Furthermore, robust check shows that the dipole anisotropy couldn't be reproduced by isotropic mock data set. However, it is still premature to claim that the Universe is anisotropic due to the small data samples and uncertainty in the current observations.

We discuss the development of novel deep learning algorithms to enable real-time regression analysis for time series data. We showcase the application of this new method with a timely case study, and then discuss the applicability of this approach to tackle similar challenges across science domains.

The unprecedented volume and rate of transient events that will be discovered by the Large Synoptic Survey Telescope (LSST) demands that the astronomical community update its followup paradigm. Alert-brokers -- automated software system to sift through, characterize, annotate and prioritize events for followup -- will be critical tools for managing alert streams in the LSST era. The Arizona-NOAO Temporal Analysis and Response to Events System (ANTARES) is one such broker. In this work, we develop a machine learning pipeline to characterize and classify variable and transient sources only using the available multiband optical photometry. We describe three illustrative stages of the pipeline, serving the three goals of early, intermediate and retrospective classification of alerts. The first takes the form of variable vs transient categorization, the second, a multi-class typing of the combined variable and transient dataset, and the third, a purity-driven subtyping of a transient class. While several similar algorithms have proven themselves in simulations, we validate their performance on real observations for the first time. We quantitatively evaluate our pipeline on sparse, unevenly sampled, heteroskedastic data from various existing observational campaigns, and demonstrate very competitive classification performance. We describe our progress towards adapting the pipeline developed in this work into a real-time broker working on live alert streams from time-domain surveys.

Gravitational wave astronomy is a rapidly growing field of modern astrophysics, with observations being made frequently by the LIGO detectors. Gravitational wave signals are often extremely weak and the data from the detectors, such as LIGO, is contaminated with non-Gaussian and non-stationary noise, often containing transient disturbances which can obscure real signals. Traditional denoising methods, such as principal component analysis and dictionary learning, are not optimal for dealing with this non-Gaussian noise, especially for low signal-to-noise ratio gravitational wave signals. Furthermore, these methods are computationally expensive on large datasets. To overcome these issues, we apply state-of-the-art signal processing techniques, based on recent groundbreaking advancements in deep learning, to denoise gravitational wave signals embedded either in Gaussian noise or in real LIGO noise. We introduce SMTDAE, a Staired Multi-Timestep Denoising Autoencoder, based on sequence-to-sequence bi-directional Long-Short-Term-Memory recurrent neural networks. We demonstrate the advantages of using our unsupervised deep learning approach and show that, after training only using simulated Gaussian noise, SMTDAE achieves superior recovery performance for gravitational wave signals embedded in real non-Gaussian LIGO noise.

The main purpose of this work is to distinguish various holographic type dark energy (DE) models, including the $\Lambda$HDE, HDE, NADE and RDE model, by using various diagnostic tools. The first diagnostic tool is the Statefinder hierarchy, in which the evolution of Statefinder hierarchy parmeter $S^{(1)}_3(z)$ and $S^{(1)}_4(z)$ are studied. The second is composite null diagnostic (CND), in which the trajectories of $\{S^{(1)}_3, \epsilon\}$ and $\{S^{(1)}_4, \epsilon\}$ are investigated, where $\epsilon$ is the fractional growth parameter. The last is $w-w'$ analysis, where $w$ is the equation of state for DE and the prime denotes derivative with respect to $ln a$. In the analysis we consider two cases: varying current fractional DE density $\Omega_{de0}$ and varying DE model parameter $C$. We find that: (1) Both the Statefinder hierarchy and the CND have qualitative impact on $\Lambda$HDE, but only have quantitative impact on HDE. (2) $S_4^{(1)}$ can lead to larger differences than $S_3^{(1)}$, while the CND pair has a stronger ability to distinguish different models than the Statefinder hierarchy. (3) For the case of varying $C$, the $\{w, w'\}$ pair has qualitative impact on $\Lambda$HDE; for the case of varying $\Omega_{de0}$, the $\{w, w'\}$ pair only has quantitative impact; these results are different from the cases of HDE, RDE and NADE, in which the $\{w, w'\}$ pair only has quantitative impact on these models. In conclusion, compared with HDE, RDE and NADE, the $\Lambda$HDE model can be easily distinguished by using these diagnostic tools.

It is well known that the Milgrom's MOND (modified Newtonian dynamics) explains well the mass discrepancy problem in galaxy rotation curves. The MOND predicts a universal acceleration scale below which the Newtonian dynamics is invalid yet. The universal acceleration scale we got from the SPARC dataset is $g_{\dag}=1.02\times10^{-10} \rm m~s^{-2}$. Milgrom suggested that the acceleration scale may be a fingerprint of cosmology on local dynamics and related with the Hubble constant $g_{\dag}\sim cH_0$. In this paper, we use the hemisphere comparison method with the SPARC dataset to investigate the spatial anisotropy on the acceleration scale. We find that the hemisphere of the maximum acceleration scale is in the direction $(l,b) = ({175.5^\circ}^{+6^\circ}_{-10^\circ}, {-6.5^\circ}^{+8^\circ}_{-3^\circ})$ with $g_{\dag,max}=1.10\times10^{-10} \rm m~s^{-2}$, while the hemisphere of the minimum acceleration scale is in the opposite direction $(l,b) = ({355.5^\circ}^{+6^\circ}_{-10^\circ}, {6.5^\circ}^{+3^\circ}_{-8^\circ})$ with $g_{\dag,min}=0.76\times10^{-10} \rm m~s^{-2}$. The maximum anisotropy level reaches up to $0.37\pm0.04$. Robust tests present that such a level of anisotropy can't be reproduced by a statistically isotropic data. In addition, we show that the spatial anisotropy on the acceleration scale has little correlation with the non-uniform distribution of the SPARC data points in sky. We also find that the maximum anisotropy direction is close with other cosmological preferred directions, especially the direction of the "Australia dipole" for the fine structure constant.

About 0.4s after the Laser Interferometer Gravitational-Wave Observatory (LIGO) detected a transient gravitational-wave (GW) signal GW150914, the Fermi Gamma-ray Burst Monitor (GBM) also found a weak electromagnetic transient (GBM transient 150914). Time and location coincidences favor a possible association between GW150904 and GBM transient 150914. Under this possible association, we adopt Fermi's electromagnetic (EM) localization and derive constraints on possible violations of the Weak Equivalence Principle (WEP) from the observations of two events. Our calculations are based on four comparisons: (1)The first is the comparison of the initial GWs detected at the two LIGO sites. From the different polarizations of these initial GWs, we obtain a limit on any difference in the parametrized post-Newtonian (PPN) parameter $\Delta\gamma\lesssim 10^{-10}$. (2) The second is a comparison of GWs and possible EM waves. Using a traditional super-Eddington accretion model for GBM transient 150914, we again obtain an upper limit $\Delta\gamma\lesssim 10^{-10}$. Compared with previous results for photons and neutrinos, our limits are five orders of magnitude stronger than those from PeV neutrinos in blazar flares, and seven orders stronger than those from MeV neutrinos in SN1987A. (3) The third is a comparison of GWs with different frequencies in the range [35 Hz, 250 Hz]. (4) The fourth is a comparison of EM waves with different energies in the range [1 keV, 10 MeV]. These last two comparisons lead to an even stronger limit, $\Delta\gamma\lesssim 10^{-8}$. Our results highlight the potential of multi-messenger signals exploiting different emission channels to strengthen existing tests of the WEP.

The response of a detector to gravitational wave is a function of frequency. When the time a photon moving around in the Fabry-Perot cavities is the same order of the period of a gravitational wave, the phase-difference due to the gravitational wave should be an integral along the path. We present a formula description for detector response to gravitational wave with varied frequencies. The LIGO data for GW150914 and GW 151226 are reexamined in this framework. For GW150924, the traveling time of a photon in the LIGO detector is just a bit larger than a half period of the highest frequency of gravitational wave and the similar result is obtained with LIGO and Virgo collaborations. However, we are not always so luck. In the case of GW151226, the time of a photon traveling in the detector is larger than the period of the highest frequency of gravitational wave and the announced signal cannot match well the template with the initial black hole masses 14.2M$_\odot$ and 7.5M$_\odot$.

The complex structures of warm and hot dense matter are essential to understand the behaviors of materials in high energy density physics processes and provide new features of matter constitutions. Here, along a new unified first-principle determined Hugoniot curve of iron from normal condensed condition up to 1 Gbar, the novel structures characterized by the ionic clusters and separated "electron bubbles" are revolutionarily unraveled using newly developed quantum Langevin molecular dynamics (QLMD). Subsistence of complex clusters, with bonds formed by inner shell electrons of neighbor ions, can persist in the time length of 50 femto-seconds dynamically with quantum flowing bubbles, which are produced by the interplay of Fermi electron degeneracy, the ionic coupling and the dynamical nature. With the inclusion of those complicated features in QLMD, the present data could serve as a first-principle benchmark in a wide range of temperatures and densities.

We present stellar parameters of 19 K-type giants and their abundances of 13 chemical elements (Al, Ba, Ca, Fe, K, Mg, Mn, Na, Ni, Sc, Si, Ti and V), selected from two moving groups, covering the metallicity range of -0.6 < [Fe/H] < 0.2, based on high resolution spectra. Most elemental abundances show similar trends with previous studies except for Al, Na and Ba, which are affected by evolution seriously. The abundance ratios of [Na/Mg] increase smoothly with higher [Mg/H] and [Al/Mg] decrease slightly with increasing [Mg/H]. [Mg/Ba] show distinction between these two moving groups which is mainly induced by chemical evolution and partly by kinematic effects. The inhomogeneous metallicity of each star from the moving groups demonstrate that these stars have different chemical origins before they were kinematically aggregated and favor the dynamical resonant theory.

The absolute branching ratio of the beta-delayed gamma-ray emission of N18 was measured, providing absolute normalization of the previous work by Olness et al. [J.W. Olness et al., Nuc. Phys. A373, 13 (1982)] who measured the relative branching ratios for the individual gamma-rays. We find the total absolute branching ratio for beta-delayed gamma emission of N18 to be 76.7 +/- 7.2(stat) +/- 5.8(norm)%. A combination of other results suggests a value consistent with our result, but smaller than that calculated by Millener as quoted in Olness et al.