The detection of binary black hole coalescences by LIGO/Virgo has aroused the interest in primordial black holes (PBHs), because they could be both the progenitors of these black holes and a compelling candidate of dark matter (DM). PBHs are formed soon after the enhanced scalar perturbations re-enter horizon during radiation dominated era, which would inevitably induce gravitational waves as well. Searching for such scalar induced gravitational waves (SIGWs) provides an elegant way to probe PBHs. We perform the first direct search for the signals of SIGWs accompanying the formation of PBHs in North American Nanohertz Observatory for Gravitational waves (NANOGrav) 11-year data set. No statistically significant detection has been made, and hence we place a stringent upper limit on the abundance of PBHs at $95\%$ confidence level. In particular, less than one part in a million of the total DM mass could come from PBHs in the mass range of $[2 \times 10^{-3}, 7\times 10^{-1}] \Msun$.

We analytically calculate the scalar induced gravitational waves (SIGWs) and find a log-dependent slope of SIGW in the infrared regions $(f<f_c)$, namely $n_{\mathrm{GW}}(f)=3-2/\ln(f_c/f)$, and $n_{\mathrm{GW}}(f)=2-2/\ln(f_c/f)$ near the peak if the power spectrum of scalar curvature perturbation is quite narrow, where $f_c$ is roughly the frequency at the peak of SIGW. Such a log-dependent slope can be taken as a new template for distinguishing SIGW from other sources.

In this letter we calculate the scalar induced gravitational waves (SIGWs) accompanying the formation of primordial black hole during the radiation dominated era in three different gauges, i.e. synchronous gauge, Newton gauge and uniform curvature gauge, and we find that the energy density spectra of SIGWs, $\ogw(k)$, are identical in these three different gauges.

We investigate the blazar contribution to the cumulative neutrino intensity assuming a generic relationship between neutrino and gamma-ray luminosities, $L_{\nu}\propto (L_{\rm ph})^{\gamma_{\rm lw}}$. Using the gamma-ray luminosity functions for blazars including flat spectrum radio quasars (FSRQs) and BL Lac objects, as well as the $Fermi$-LAT detection efficiency, we estimate contributions from {blazars resolved by $Fermi$-LAT as well as the unresolved counterpart.} Combining the existing upper limits from stacking analyses, the cumulative neutrino flux from all blazars ({including $Fermi$-LAT resolved and unresolved ones}) are constrained { in the range $0\lesssim\gamma_{\rm lw}\lesssim2.5$}. We also evaluate the effects of the redshift evolution and the effective local number densities for each class of FSRQs, BL Lacs, and all blazars, by which we place another type of constraints on the blazar contribution using the non-detection of high-energy neutrino multiplets. We demonstrate that these two upper limits are complementary, and that the joint consideration of the stacking and multiplet analyses { not only supports the argument that blazars are disfavored as the dominant sources of the 100-TeV neutrino background, but it extends this argument by including also $Fermi$-LAT-unresolved blazars as well, for a more generic luminosity correlation $L_{\nu}\propto (L_{\rm ph})^{\gamma_{\rm lw}}$}.

The possibility that primordial black holes (PBHs) represent all of the dark matter (DM) in the Universe and explain the coalescences of binary black holes detected by LIGO/Virgo has attracted a lot of attention. PBHs are generated by the enhancement of scalar perturbations which inevitably produce the induced gravitational waves (GWs). We calculate the induced GWs up to the third-order correction which not only enhances the amplitude of induced GWs, but also extends the cutoff frequency from $2k_*$ to $3k_*$. Such effects of the third-order correction lead to an around $10\%$ increase of the signal-to-noise ratio (SNR) for both LISA and pulsar timing array (PTA) observations, and significantly widen the mass range of PBHs in the stellar mass window accompanying detectable induced GWs for PTA observations including IPTA, FAST and SKA. On the other hand, the null detections of the induced GWs by LISA and PTA experiments will exclude the possibility that all of the DM is comprised of PBHs and the GW events detected by LIGO/Virgo are generated by PBHs.

Shocks arising in galaxy mergers could accelerate cosmic-ray (CR) ions to TeV-PeV energies. While propagating in the intergalactic medium, these CRs can produce high-energy neutrinos, electron-positron pairs and gamma-rays. In the presence of intergalactic magnetic fields, the secondary pairs will radiate observable emissions through synchrotron radiation and inverse Compton scattering. In this paper, we demonstrate that these emissions can explain the radio and X-ray fluxes of merging galaxies such as NGC 660 and NGC 3256. Using our model in combination with the observations, we can constrain the gas mass, shock velocity, magnetic field and the CR spectral index $s$ of these systems. For NGC 660 a single-zone model with a spectral index $2.1\lesssim s\lesssim2.2$ is able to reproduce simultaneously the radio and X-ray observations, while a simple one-zone scenario with $s\sim2$ can describe the radio and a large fraction of X-ray observations of NGC 3256. Our work provides a useful approach for studying the dynamics and physical parameters of galaxy mergers, which can play an important part in future multi-messenger studies of similar and related extragalactic sources.

The merger of dark matter halos and the gaseous structures embedded in them, such as proto-galaxies, galaxies, and groups and clusters of galaxies, results in strong shocks that are capable of accelerating cosmic rays (CRs) to $\sim10~\rm PeV$. These shocks will produce high-energy neutrinos and $\gamma$-rays through inelastic $pp$ collisions with ambient gaseous environments. In this work, we study the contributions of these halo mergers to the diffuse neutrino flux measured in IceCube and to the non-blazar portion of the extragalactic $\gamma$-ray background measured by $Fermi$. In order to calculate them, we formulate the redshift dependence of the shock velocity, galactic radius, halo gas content and galactic/intergalactic magnetic fields over the dark matter halo distribution up to a redshift $z=10$. We find that high-redshift mergers contribute a significant amount of the cosmic-ray energy luminosity density, and the resulting neutrino spectra could explain a large part of the observed diffuse neutrino flux above 0.1 PeV up to several PeV. We also show that our model can somewhat alleviate tensions with the extragalactic $\gamma$-ray background. First, since a larger fraction of the CR energy luminosity density comes from high redshifts, the accompanying $\gamma$-rays are more strongly suppressed through $\gamma\gamma$ annihilations with the cosmic microwave background (CMB) and the extragalactic background light (EBL). Second, mildly radiative-cooled shocks may lead to a harder CR spectrum with spectral indices of $1.5\lesssim s\lesssim2.0$. Our study suggests that halo mergers, a fraction of which may also induce starbursts in the merged galaxies, can be promising neutrino emitters without violating the existing $Fermi$ $\gamma$-ray constraints on the non-blazar component of the extragalactic $\gamma$-ray background.

As one of the probes of universe, strong gravitational lensing systems allow us to compare different cosmological models and constrain vital cosmological parameters. This purpose can be reached from the dynamic and geometry properties of strong gravitational lensing systems, for instance, time-delay $\Delta\tau$ of images, the velocity dispersion $\sigma$ of the lensing galaxies and the combination of these two effects, $\Delta\tau/\sigma^2$. In this paper, in order to carry out one-on-one comparisons between $\Lambda$CDM universe and $R_h=ct$ universe, we use a sample containing 36 strong lensing systems with the measurement of velocity dispersion from the SLACS and LSD survey. Concerning the time-delay effect, 12 two-image lensing systems with $\Delta\tau$ are also used. In addition, Monte Carlo (MC) simulations are used to compare the efficiency of the three methods as mentioned above. From simulations, we estimate the number of lenses required to rule out one model at the $99.7\%$ confidence level. Comparing with constraints from $\Delta\tau$ and the velocity dispersion $\sigma$, we find that using $\Delta\tau/\sigma^2$ can improve the discrimination between cosmological models. Despite the independence tests of these methods reveal a correlation between $\Delta\tau/\sigma^2$ and $\sigma$, $\Delta\tau/\sigma^2$ could be considered as an improved method of $\sigma$ if more data samples are available.

NGC 6782 is an early-type barred spiral galaxy exhibiting a rich and complex morphology with multiple ring patterns. To provide a physical understanding of its structure and kinematical properties, two-dimensional hydrodynamical simulations have been carried out. Numerical calculations reveal that the striking features in NGC 6782 can be reproduced provided that the gas flow is governed by the gravitational potential associated with a slowly rotating strong bar. In particular, the response of the gaseous disk to the bar potential leads to the excitation of spiral density waves at the inner Lindblad resonance giving rise to the appearance of a nearly circular nuclear ring with a pair of dust lanes. For a sufficiently strong bar potential, the inner 4:1 spiral density waves are also excited. The interaction of the higher harmonic waves with the waves excited at the inner Lindblad resonance and confined by the outer Lindblad resonance results in the observed diamond-shaped (or pointy oval) inner ring structure. The overall gas morphology and kinematical features are both well reproduced by the model provided that the pattern speed of the bar is $\sim 25$ km s$^{-1}$ kpc$^{-1}$.

Dark matter (comprising a quarter of the Universe) is usually assumed to be due to one and only one weakly interacting particle which is neutral and absolutely stable. We consider the possibility that there are several coexisting dark-matter particles, and explore in some detail the generic case where there are two. We discuss how the second dark-matter particle may relax the severe constraints on the parameter space of the Minimal Supersymmetric Standard Model, as well as other verifiable predictions in both direct and indirect search experiments.

We performed a series of hydro-dynamic simulations to investigate the orbital migration of a Jovian planet embedded in a proto-stellar disk. In order to take into account of the effect of the disk's self gravity, we developed and adopted an \textbf{Antares} code which is based on a 2-D Godunov scheme to obtain the exact Reimann solution for isothermal or polytropic gas, with non-reflecting boundary conditions. Our simulations indicate that in the study of the runaway (type III) migration, it is important to carry out a fully self consistent treatment of the gravitational interaction between the disk and the embedded planet. Through a series of convergence tests, we show that adequate numerical resolution, especially within the planet's Roche lobe, critically determines the outcome of the simulations. We consider a variety of initial conditions and show that isolated, non eccentric protoplanet planets do not undergo type III migration. We attribute the difference between our and previous simulations to the contribution of a self consistent representation of the disk's self gravity. Nevertheless, type III migration cannot be completely suppressed and its onset requires finite amplitude perturbations such as that induced by planet-planet interaction. We determine the radial extent of type III migration as a function of the disk's self gravity.

The Seyfert galaxy NGC 1097 has an extended neutral hydrogen disk, a companion, a prominent bar and a luminous circumnuclear starburst "ring". Magnetic fields as revealed by nonthermal radio-continuum emissions correlate well with the optical barred spiral structure on large scales, have a gross enhancement overlapping with the optical/infrared "ring", and show a trailing swirl around and within the "ring". We propose a scenario of bar-excited long-trailing fast magnetohydrodynamic (MHD) density waves at the modified inner Lindblad resonance (mILR), physically identified with the outer rim of the "ring". These sustained outgoing long-waves are bounced back by the $Q_M$-barrier in the form of incoming short-trailing waves. The damping of these waves deposits a {\it negative} angular momentum into the magnetized circumnuclear gas disk. Thus, gas materials spiral inward, bring in frozen-in magnetic flux, and accumulate inside the mILR to create a circular zone of high density and magnetic flux vulnerable to massive star formation. Depending on the wave damping efficiency, this process may simultaneously sustain a net mass inflow across the "ring" and toward the nucleus. A wavelet analysis on a Hubble Space Telescope image of central NGC 1097 shows a distinct two-arm spiral structure extended down to the nucleus as a strong evidence for circumnuclear MHD density waves. We predict that magnetic-field observations with improved sensitivity and resolution would reveal a specific correspondence between circumnuclear optical and magnetic field spirals much as those known to exist on large scales in nearby spiral galaxies, including NGC 1097.

An improved linear stability theory of small-amplitude oscillations of a self-gravitating, infinitesimally thin gaseous disk of spiral galaxies has been developed. It was shown that in the differentially rotating disks for nonaxisymmetric perturbations Toomre's modified hydrodynamical simulations to test the validity of the modified local criterion.

The kinetic theory is used to study the evolution of the self-gravitating disk of planetesimals. The effects of frequent collisions between planetesimals are taken into account by using a Krook integral in the Boltzmann kinetic equation. It is shown that as a result of an aperiodic collision-dissipative instability of small gravity disturbances the disk is subdivided into numerous dense fragments. These can eventually condense into the planetary sequence.

Linear theory is used to determine the stability of the self-gravitating, rapidly (and nonuniformly) rotating, two-dimensional, and collisional particulate disk against small-amplitude gravity perturbations. A gas-kinetic theory approach is used by exploring the combined system of the Boltzmann and the Poisson equations. The effects of physical collisions between particles are taken into account by using in the Boltzmann kinetic equation a Krook model integral of collisions modified to allow collisions to be inelastic. It is shown that as a direct result of the classical Jeans instability and a secular dissipative-type instability of small-amplitude gravity disturbances (e.g. those produced by a spontaneous perturbation and/or a companion system) the disk is subdivided into numerous irregular ringlets, with size and spacing of the order of the Jeans length (= (4-6) mean epicyclic radius). The present research is aimed above all at explaining the origin of various structures in highly flattened, rapidly rotating systems of mutually gravitating particles. In particular, it is suggested that forthcoming Cassini spacecraft high-resolution images may reveal this kind of hyperfine of the order of 100 m structure in the main rings A, B, and C of the Saturnian ring system.

The reaction of collective oscillations excited in the interaction between aperiodically growing Jeans-type gravity perturbations and stars of a rapidly rotating disk of flat galaxies is considered. An equation is derived which describes the change in the main body of equilibrium distribution of stars in the framework of the nonresonant weakly nonlinear theory. Certain applications of the theory to the problem of relaxation of the Milky Way at radii where two-body relaxation is not effective are explored. The theory, as applied to the Solar neighborhood, accounts for observed features, such as the shape for the velocity ellipsoid of stars and the increase in star random velocities with age.

The time evolution of initially balanced, rapidly rotating models for an isolated disk of highly flattened galaxies of stars is calculated. The method of direct integration of the Newtonian equations of motion of stars over a time span of many galactic rotations is applied. Use of modern concurrent supercomputers has enabled us to make long simulation runs using a sufficiently large number of particles N=30,000. One of the goals of the present simulation is to test the validities of a modified local criterion for stability of Jeans-type gravity perturbations (e.g. those produced by a barlike structure, a spontaneous perturbation and/or a companion galaxy) in a self-gravitating, infinitesimally thin and collisionless disk. In addition to the local criterion we are interested in how model particles diffuse in velocity. This is of considerable interest in the kinetic theory of stellar disks. Certain astronomical implications of the simulations to actual disk-shaped (i.e. rapidly rotating) galaxies are explored. The weakly nonlinear, or quasi-linear kinetic theory of the Jeans instability in disk galaxies of stars is described as well.

Recent NICMOS observations of the central regions of nearby barred galaxies, NGC5383, NGC1530 and NGC1667 by Regan, Sheth and Mulchaey (1999), all show, in the unsharp mask map, a pair of distinct trailing spirals which can be followed all the way to the center. This result is in apparent conflict with the density wave theory. In this report, we argue that this difficulty may be removed if we take the self-gravitation of the disc into consideration, or the IILR does not exist in these galaxies.

The stability of a rotating fluid disk to the formation of spiral arms is studied in the tightwinding approximation in the linear regime. The dispersion relation for spirals that was derived by Bertin et al. is shown to contain a new, acoustic instability beyond the Lindblad resonances that depends only on pressure and rotation. In this regime, pressure and gravity exchange roles as drivers and inhibitors of spiral wave structures. Other instabilities that are enhanced by pressure are also found in the general dispersion relation by including higher order terms in the small parameter 1/kr for wavenumber k and radius r. These instabilities are present even for large values of Toomre's parameter Q. Unstable growth rates are determined in four cases: a self-gravitating disk with a flat rotation curve, a self-gravitating disk with solid body rotation, a non-self-gravitating disk with solid body rotation, and a non-self-gravitating disk with Keplerian rotation. The most important application appears to be as a source of spiral structure, possibly leading to accretion in non-self-gravitating disks, such as some galactic nuclear disks, disks around black holes, and proto-planetary disks. All of these examples have short orbital times so the unstable growth time can be small.

Observations with the Hubble Space Telescope reveal an irregular network of dust spiral arms in the nuclear region of the interacting disk galaxy NGC 2207. The spirals extend from ~50 pc to ~300 pc in galactocentric radius, with a projected width of ~20 pc. Radiative transfer calculations determine the gas properties of the spirals and the inner disk, and imply a factor of ~4 local gas compression in the spirals. The gas is not strongly self-gravitating, nor is there a nuclear bar, so the spirals could not have formed by the usual mechanisms applied to main galaxy disks. Instead, they may result from acoustic instabilities that amplify at small galactic radii. Such instabilities may promote gas accretion into the nucleus.