The LIGO/Virgo Scientific Collaboration (LVC) recently announced the detection of a compact object binary merger, GW190425, with a total mass of $3.4^{+0.3}_{-0.1}$ M$_{\odot}$, and individual component masses in the range of about 1.1 to 2.5 $M_{\odot}$. If the constituent compact objects are neutron stars, then the total mass is five standard deviations higher than the mean of $2.66\pm 0.12$ M$_{\odot}$ for Galactic binary neutron stars (BNSs). LVC suggests such massive BNS systems are born from a fast-merging channel, and therefore, their non-detection in the Galaxy to be due to a selection effect. However, we are unable to reconcile the inferred formation efficiency from the reported merger rate, $\mathcal{R}_{\rm GW190425}=460^{+1050}_{-390}$ yr$^{-1}$ Gpc$^{-3}$, with predictions from our own study for fast-merging BNS systems. Moreover, the comparable merger rates of GW190425, and GW170817 are possibly in tension with our results for two reasons: (i) more massive systems are expected to have a lower formation rate, and (ii) fast merging channels should constitute $\lesssim 10\%$ of the total BNS systems if case BB unstable mass transfer is permitted to take place as a formation pathway. We argue that to account for the high merger rate of GW190425 as a BNS system requires: (i) revisiting our understanding of NS formation in supernova explosions, or (ii): that more massive NSs need to be preferentially born with either very weak or very high magnetic fields so that they would be undetectable in the radio surveys. Perhaps massive NSs detected in NS-white dwarf binaries are our clues to the formation path of GW190425 systems.

The heaviest neutron stars and lightest black holes expected to be produced by stellar evolution leave the mass-range $2.2$ M$_{\odot}\lesssim m \lesssim 5$ M$_\odot$ largely unpopulated. Objects found in this so-called lower mass gap likely originate from a distinct astrophysical process. Such an object, with mass $2.6$ M$_\odot$ was recently detected in the binary merger GW190814 through gravitational waves by LIGO/Virgo. Here we show that black holes in the mass gap are naturally assembled through mergers and accretion in AGN disks, and can subsequently participate in additional mergers. We compute the properties of AGN-assisted mergers involving neutron stars and black holes, accounting for accretion. We find that mergers in which one of the objects is in the lower mass gap represent up to $4$% of AGN-assisted mergers detectable by LIGO/Virgo. The lighter object of GW190814, with mass $2.6$ M$_\odot$, could have grown in an AGN disk through accretion. We find that the unexpectedly high total mass of 3.4 M$_\odot$ observed in the neutron star merger GW190425 may also be due to accretion in an AGN disk.

HR 6819 was recently claimed to be a hierarchical triple system of a Be star in a wide orbit around an inner binary system of a black hole (BH) and a B III type star. We argue that this system is unlikely to be a hierarchical triple due to three reasons: (i) Given that this system is discovered in a magnitude limited Bright Star Catalog, the expected number of such systems in the Milky Way amounts to about $10^4$ while the estimate for the MW budget for such systems is between $10^2-10^3$ systems under generous assumptions. Such a large gap cannot be reconciled as it would otherwise likely overflow the MW budget for BHs; (ii) The dynamical stability of this system sets lower bounds on the orbital separation of the outer Be star, while it not being resolved by Gaia places an upper limit on its projected sky separation. We show that these two constraints would imply a narrow range for the outer orbit without resorting to geometrical fine-tuning; (iii) The triple system should have survived the stellar evolution prior to the formation of the BH in the inner binary. We perform numerical simulations starting with conservative initial conditions of this system and show that a small parameter space for BH progenitor star's mass loss, BH natal kicks, and initial orbital separation can reproduce HR 6819. Therefore, we propose this system is a chance superposition of a Be star with a binary.

The LIGO/Virgo Scientific Collaboration recently announced the detection of a compact object binary merger, GW190412, as the first asymmetric binary black hole (BBH) merger with mass ratio $q\approx0.25$. Other than the mass ratio, this BBH has shown to have a positive effective spin of around $\chi_{\rm eff}\approx0.28$. Assuming a field formation channel, associating this effective spin to either the primary or the secondary BH each has its implications: If the spin of the BBH comes form the primary BH, it has consequences for the efficiency of angular momentum transport in the formation of the BH. If, on the other hand, the spin is due to the secondary BH through tidal spin-up processes, one has to note that (i) such processes have very short delay-times, and (ii) subsequently, their local merger rate is determined by local star formation rate at assumed formation metallicity of the BBH. We show that the predicted merger rate density from this channel is $\lesssim 0.3~\rm Gpc^{-3} yr^{-1}$ and in tension with the rather high local merger rate of such systems which we estimate from this single event to be $\sim 1.7^{+2.5}_{-1.4}~\rm Gpc^{-3} yr^{-1}$ (90\% confidence interval, and assuming 50 days of observing time). Large natal kicks ($v\gtrsim 500\,{\rm km/s}$) would be required to get such BBHs with an in-plane spin component to account for the marginal detection of precession in GW190412. However, this would only exacerbate the tension as the estimated local merger rate would be further decreased. Similarly, the formation of such systems through the dynamical assembly is exceedingly rare, leaving this system a dilemma hard to account for with the currently accepted paradigms of BBH formation.

The LIGO/Virgo Collaboration (LVC) recently reported the detection of GW190814, a merger of a $23^{+1.0}_{-0.9}~M_{\odot}$ primary black hole (BH), and a $2.6^{+0.08}_{-0.08}~M_{\odot}$ secondary. The secondary's mass falls into the mass-gap regime, which refers to the scarcity of compact objects in the mass range of 2-5 $M_{\odot}$. The first clue to the formation of the GW190814 lies in the fact that the primary is a very massive BH. We suggest that the secondary was born as a neutron star (NS) where a significant amount of the supernova ejecta mass from its formation remained bound to the binary due to the presence of the massive BH companion. The bound mass forms a circumbinary accretion disk, and its accretion onto the NS created a mass-gap object. In this scenario, LIGO/Virgo will only detect mass-gap objects in binary mergers with an extreme mass ratio. We also predict a correlation between the mass of the secondary and the mass of the primary in such asymmetric mergers. Our model can be tested with future data from the LVC's third-observing run.

The gravitational wave source GW190412 is a binary black hole (BBH) merger with three unique properties: i) its mass ratio is about 0.28, the lowest found so far, ii) it has a relatively high positive effective spin parameter \chi_eff=0.25, and iii) it is observed to be precessing due to in-plane projected spin of the binary with an in-plane precession parameter \chi_p=0.3. The two main formation channels of BBH formation fail to account for GW190412: field formation scenarios cannot explain the observed precession unless by invoking large natal kicks, and dynamical assembly in dense stellar systems is inefficient in producing such low mass-ratio BBH mergers. Here, we investigate whether "double mergers" in wide hierarchical quadruple systems in the "3+1" configuration could explain the unique properties of GW190412. In this scenario, a compact object quadruple system experiences two mergers: first, two compact objects in the innermost orbit merge due to secular chaotic evolution. At a later time, the merged compact object coalesces with another compact object due to secular Lidov-Kozai oscillations. We find that our scenario is consistent with GW190412. In particular, we find a preferential projected spin around \chi_p=0.2. However, the likelihood of a double merger is small and the formation efficiency of these systems is uncertain. If GW190412 originated from a double merger in a 3+1 quadruple, we find a strong constraint that the first merger likely occurred between roughly equal-mass BHs in the innermost orbit, since the recoil velocity from unequal-mass BHs would otherwise have disrupted the system.

The circumgalactic medium (CGM) of nearby star-forming galaxies shows clear indications of OVI absorption accompanied by little to no detectable NV absorption. This unusual spectral signature, accompanied by highly non-uniform absorption from lower ionization state species, indicates that the CGM must be viewed as a dynamic, multiphase medium, such as occurs in the presence of turbulence. Motivated by previous isotropic turbulent simulations, we carry out chemodynamical simulations of stratified media in a Navarro-Frenk-White (NFW) gravitational potential with a total mass of $10^{12}$ solar masses and turbulence that decreases radially. The simulations assume a metallicity of 0.3 solar, a redshift zero metagalatic UV background, and they track ionizations, recombinations, and species-by-species radiative cooling using the MAIHEM package. We compare a suite of ionic column densities with the COS-Halos sample of low-redshift star-forming galaxies. Turbulence with an average one-dimensional velocity dispersion approximately 40 km/s, corresponding to an energy injection rate of approximately $10^{49}$ erg/yr, produces a CGM that matches many of the observed ionic column densities and ratios. In this simulation, the NVI to OVI ratio is suppressed from its equilibrium value due to a combination of radiative cooling and cooling from turbulent mixing. This level of turbulence is consistent with expectations from observations of better constrained, higher-mass systems, and could be sustained by energy input from supernovae, gas inflows, and dynamical friction from dark matter subhalos. We also conduct a higher resolution run which yields smaller-scale structures, but remains in agreement with observations.

The initial separation of massive star binaries sets the timescale over which their compact remnants merge through the emission of gravitational waves. We show that the delay time distribution (DTD) of binary neutron stars or black holes can be inferred from the stochastic gravitational wave background (SGWB). If the DTD of a population is long, most of the mergers take place at low redshifts and the background would be rather quiet compared to a scenario in which the DTD is short leading to few individual detections at low redshift but a rather loud background. We show that different DTDs predict a factor of 5 difference in the magnitude of the gravitational wave background energy density ($\Omega_{\rm GW}$) and have the dominant effect on $\Omega_{\rm GW}$ over other factors such as the mass function of the primary BH mass, $p(m_1)$, the maximum considered BH mass ($M_{\rm max}$), and the effective spin of the black hole ($\chi_{\rm eff}$). A non-detection of such a background can rule out the short DTD scenario. We show that SGWB searches can rule out the short DTD scenario for the BBHs within about four years of observing time at advanced LIGO design sensitivty for a local merger rate of 30 $\rm Gpc^{-3} yr^{-1}$ assuming $p(m_1)\propto m_1^{-1}$, and $M_{\rm max}=50 M_{\odot}$.

Binary black holes (BBHs) detected by gravitational wave (GW) observations could be broadly divided into two formation channels: those formed through field binary evolution and those assembled dynamically in dense stellar systems. Each of these formation channels, and their sub-channels, populate a distinct region in the effective spin-mass ($\chi_{\rm eff}-M$) plane. Depending on the branching ratio of different channels, an ensemble of BBHs could show a trend in this plane. Here we fit a mass-dependent distribution for $\chi_{\rm eff}$ to the GWTC-1 BBHs from the first and second observing runs of Advanced LIGO and Advanced Virgo. We find a negative correlation between mass and the mean effective spin ($\bar{\chi}_{\mathrm{eff}}$), and positive correlation with its dispersion ($\sigma_{\chi_\mathrm{eff}}$) at 75\% and 80\% confidence. This trend is robust against the choice of mass variable, but most pronounced when the mass variable is taken to be the chirp mass of the binary. The result is consistent with significant contributions from both dynamically assembled and field binaries in the GWTC-1 catalog. The upcoming LIGO O3a data release will critically test this interpretation.

Formation of binary black holes (BBHs) detected by gravitational-wave (GW) observations could be broadly divided into two categories: those formed through field binary evolution and those assembled dynamically in dense stellar systems. The branching ratio of the BBHs refers to the contribution of each channel. The dynamical assembly channel would predict a symmetric distribution in the effective spins of the BBHs while field formation predicts BBHs to have positive effective spins. By modeling these two populations based on their effective spin distribution we show that in the 10 BBHs detected by LIGO/Virgo the contribution of the dynamically assembled BBHs to be more than about 50\% with 90\% confidence. This result is based on the assumption that the field binaries are born with positive effective spins not restricted to have small values.

The density profiles of dwarf galaxies are a highly varied set. If the dark matter is an ultra-light particle such as axions, then simulations predict a distinctive and unique profile. If the axion mass is large enough to fit the ultra-faint dwarf (UFD) satellites($m\gtrapprox 10^{-21}$ eV), then the models do not fit the density profile of Fornax and Sculptor and are ruled out by more than $3-\sigma$ confidence. If the axion mass is in the mass range that can fit mass profiles of Fornax and Sculptor dwarf spheroidals, then its extended profile implies enormous masses ($\approx10^{11}-10^{12}M_{\odot}$) for the UFDs. These large masses for the UFDS are ruled out by more than $3-\sigma$ confidence by dynamical friction arguments. The tension would increase further considering star formation histories and stellar masses of the UFDs. Unless future ultra-light dark matter (ULDM) simulations with baryonic feedback show a significant change in the density structure of the halos, the current data is incompatible with the ULDM scenario. Relaxing the slope constraint from classical dwarf galaxies would lead to excluding ULDM with mass less than $6\times10^{-22}$ eV.

If confirmed, a wide binary system of 70 $M_{\odot}$ black hole (BH) and an 8 $M_{\odot}$ main sequence star (LB-1) is observed to reside in the Milky Way (MW). We show that long term evolution of an 8 $M_{\odot}$ star around a BH with mass between 5-70 $M_{\odot}$ makes them visible as ultra-luminous X-ray (ULX) sources in the sky. Given the expected ULX phase lifetime ($\approx0.1$ Myr) and their lack of detection in the MW, we conclude that the frequency of an 8-20 $M_{\odot}$ star to be in binary around a stellar mass BH should be less than ($f<2\times10^{-3}$). This is in tension with Liu et al. (2019) claimed detection frequency of LB-1 like system around 8-20$M_{\odot}$ stars ($f\approx3\times10^{-2}$). Moreover, the 8 $M_{\odot}$ star is likely to end as a neutron star (NS) born with a very small kick from an electron capture supernova (ECSN), leaving behind a wide NS-BH binary. So far less than 1\% of all the detectable pulsars in the MW are mapped and there has been no detection of any pulsars in binary systems around BHs which sets an upper bound of about 100 possible pulsar-BH systems in the MW. We show if the NS is born from ECSN, a frequency upper limit of ($f=5\times10^{-4}$) for stars with masses $\approx 8-20~M_{\odot}$ in the MW to have a BH companion. The rate discrepancy will further increase as more pulsars are mapped in the MW, yet these searches would not be able to rule out the Liu et al. detection frequency if NSs are instead born in core collapse SNe with the commonly inferred high kick velocities.

We investigate secular evolution in hierarchical quadruple systems as a formation channel of mass-gap black holes (with masses of about $3-5$ $M_{\odot}$) in systems that will eventually lead to binary black hole mergers detectable by ground-based gravitational wave detectors (LIGO/Virgo). We show that in a 3+1 hierarchical system, two episodes of induced mergers would first cause two neutron stars to merge and form a mass-gap black hole, which will subsequently merge with another (more massive) black hole through a second induced merger. We demonstrate that such systems are stable to flybys, and their formation would predict a high mass ratio and eccentric merger of a mass-gap black hole with a more massive black hole companion. Such a formation channel may explain observed gravitational wave events such as the recently-discovered LIGO/Virgo events S190814bv and S190924h.

Using a Milky Way double neutron star (DNS) merger rate of 210 Myr$^{-1}$, as derived by the Laser Interferometer Gravitational-Wave Observatory (LIGO), we demonstrate that the Laser Interferometer Space Antenna (LISA) will detect on average 240 (330) DNSs within the Milky Way for a 4-year (8-year) mission with a signal-to-noise ratio greater than 7. Even adopting a more pessimistic rate of 42 Myr$^{-1}$, as derived by the population of Galactic DNSs, we find a significant detection of 46 (65) Milky Way DNSs. These DNSs can be leveraged to constrain formation scenarios. In particular, traditional NS-discovery methods using radio telescopes are unable to detect DNSs with $P_{\rm orb}$ $\lesssim$1 hour (merger times $\lesssim$10 Myr). If a fast-merging channel exists that forms DNSs at these short orbital periods, LISA affords, perhaps, the only opportunity to observationally characterize these systems; we show that toy models for possible formation scenarios leave unique imprints on DNS orbital eccentricities, which may be measured by LISA for values as small as $\sim$10$^{-2}$.

We can be biased against observing massive black holes to merge in the local universe as the bounds on the maximum black hole mass ($M_{BH}^{\rm max}$) depends on the assumptions regarding the metallicity evolution of the star forming gas across the cosmic time. We investigate the bounds on the metallicity evolution, mass distribution and delay times of the binary black hole sources based on the ten observed events by LIGO. We parametrize $M_{BH}^{\rm max}$ to be a function of metallicity which itself is modeled to evolve with redshift in either a modest or rapid fashion. Rapid metallicity evolution models predict a stringent bound of $M_{BH}^{\rm max}=44^{+9}_{-5}M_{\odot}$, while the bound on $M_{BH}^{\rm max}$ in the models with modest metallicity evolution is $M_{BH}^{\rm max}=52^{+16}_{-9} M_{\odot}$. Therefore, inferring $M_{BH}^{\rm max}$ from GW data depends on the assumed metal enrichment history of the universe that is not severely constrained at the moment.

The Event Horizon Telescope (EHT), a global very long baseline interferometric array observing at a wavelength of 1.3 mm, detected the first image of the M87 supermassive black hole (SMBH). M87 is a giant elliptical galaxy at the center of Virgo cluster, which is expected to have formed through merging of cluster galaxies. Consequently M87* hosted mergers of black holes through dynamical friction and could have one or multiple binary companions with a low mass ratio at large separations. We show that a long-term monitoring of the M87 SMBH image over $\sim$1 year with absolute positional accuracy of 1$\approx\mu as$ could detect such binary companions and exclude a large parameter space in semi major axis ($a_0$) and mass ratio ($q$), which is currently not constrained. Moreover, the presence of the accretion disk around M87* excludes a binary companion with $a_0\approx$ of order a mili parsec, as otherwise the accretion disk would have been tidally disrupted.

In Paper I (Safarzadeh & Berger 2019) we studied the determination of the delay time distribution (DTD) of binary neutron stars (BNS) through scaling relations between halo/stellar mass and the star formation history (SFH) of galaxies hosting gravitational wave (GW) events in the local universe. Here we explore how a detailed reconstruction of the individual SFHs of BNS merger host galaxies can improve on the use of the scaling relations. We use galaxies from the Galaxy and Mass Assembly (GAMA) survey, which is mass complete at $M_*>10^9$ M$_\odot$ in the redshift range $0.05<z<0.08$. We use the reconstructed SFHs derived from the Prospector code, for two distinct sets of priors (favoring continuous and bursty SFHs), and convolve those with power law DTDs characterized by an index $\Gamma$ and a minimum delay time $t_{\rm min}$. We find that with this approach $\mathcal{O}(100)-\mathcal{O}(300)$ host galaxies are required to constrain the DTD parameters, with the number depending on the choice of SFH prior and on the parameters of the true DTD. We further show that using only the host galaxies of BNS mergers, as opposed to the full population of potential host galaxies in the relevant cosmic volume, leads to a minor bias in the recovered DTD parameters. The required host galaxy sample size is nearly an order of magnitude smaller relative to the approach of using scaling relations, and we expect such a host galaxy sample to be collected within a decade or two, prior to the advent of third-generation GW detectors.

The presence of primordial magnetic fields increases the minimum halo mass in which star formation is possible at high redshifts. Estimates of the dynamical mass of ultra-faint dwarf galaxies (UFDs) within their half-light radius constrain their virialized halo mass before their infall into the Milky Way. The inferred halo mass and formation redshift of the UFDs place upper bounds on the primordial comoving magnetic field, $B_0$. We derive an upper limit of $0.50\pm 0.086$ ($0.31\pm 0.04$) nG on $B_0~$ assuming the average formation redshift of the UFD host halos is $z_{\rm form}=$ 10 (20), respectively.

We investigate the ability of current and third-generation gravitational wave (GW) detectors to determine the delay time distribution (DTD) of binary neutron stars (BNS) through a direct measurement of the BNS merger rate as a function of redshift. We assume that the DTD follows a power law distribution with a slope $\Gamma$ and a minimum merger time $t_{\rm min}$, and also allow the overall BNS formation efficiency per unit stellar mass to vary. By convolving the DTD and mass efficiency with the cosmic star formation history, and then with the GW detector capabilities, we explore two relevant regimes. First, for the current generation of GW detectors, which are only sensitive to the local universe, but can lead to precise redshift determinations via the identification of electromagnetic counterparts and host galaxies, we show that the DTD parameters are strongly degenerate with the unknown mass efficiency and therefore cannot be determined uniquely. Second, for third-generation detectors such as Einstein Telescope (ET) and Cosmic Explorer (CE), which will detect BNS mergers at cosmological distances, but with a redshift uncertainty inherent to GW-only detections ($\delta(z)/z\approx 0.1z$), we show that the DTD and mass efficiency can be well-constrained to better than 10\% with a year of observations. This long-term approach to determining the DTD through a direct mapping of the BNS merger redshift distribution will be supplemented by more near term studies of the DTD through the properties of BNS merger host galaxies at $z\approx 0$ (Safarzadeh & Berger 2019).

The delay time distribution of (DTD) of binary neutron stars (BNS) remains poorly constrained, mainly by the small known population of Galactic binaries, the properties of short gamma-ray burst host galaxies, and inferences from $r$-process enrichment. In the new era of BNS merger detections through gravitational waves (GW), a new route to the DTD is the demographics of the host galaxies, localized through associated electromagnetic counterparts. This approach takes advantage of the correlation between star formation history (SFH) and galaxy mass, such that the convolution of the SFH and DTD impacts the BNS merger rate as a function of galaxy mass. Here we quantify this approach for a power law DTD governed by two parameters: the power law index ($\Gamma$) and a minimum delay time ($t_{\rm min}$). Under the reasonable assumption that EM counterparts are likely only detectable in the local universe, accessible by the current generation of GW detectors, we study how many host galaxies at $z\sim 0$ are required to constrain the DTD parameters. We find that the DTD is mainly imprinted in the statistics of massive galaxies (stellar mass of $M_*\gtrsim 10^{10.5}$ M$_\odot$, comparable to the host galaxy of GW170817). Taking account of relevant uncertainties we find that $\mathcal{O}(10^3)$ host galaxies are required to constrain the DTD; for a fixed value of $t_{\rm min}$, as done in previous analyses of the DTD, $\mathcal{O}(10^2)$ host galaxies will suffice. Such a sample might become available within the next two decades, prior to the advent of third-generation GW detectors.

We study the recently observed JellyFish galaxies (JFGs), which are found to have their gas content ram-pressure stripped away in galaxy clusters. These galaxies are observed to have an enhanced star formation rate of about 0.2 dex compared with a control sample of the same stellar mass in their disks. We model the increase in the star formation efficiency as a function of interacluster medium pressure and parametrize the cold gas content of the galaxies as a function of cluster-centric distance. We show that regarding the external pressure as a positive feedback results in agreement with the observed distribution of enhanced star formation in the JFGs if clouds are shielded from evaporation by magnetic fields. Our results predict that satellites with halo mass $< 10^{11}M_{\odot}$ moving with Mach numbers $\mathcal{M}\approx2$, and inclination angles below 60 degrees, are more likely to be detected as JFGs.

We model the history of Galactic r-process enrichment using high-redshift, high-resolution zoom cosmological simulations of a Milky Way (MW) type halo. We assume that all r-process sources are neutron star mergers (NSMs) with a power law delay time distribution. We model the time to mix pollutants at subgrid scales, which allows us to to better compute the properties of metal poor (MP) and carbon enhanced metal poor (CEMP) stars, along with statistics of their r-process enhanced subclasses. Our simulations underpredict the cumulative ratios of r-process enhanced MP and CEMP stars (MP-r, CEMP-r) over MP and CEMP stars by about one order of magnitude, even when the minimum coalescence time of the double neutron stars ($t_{\rm min}$) is set to 1 Myr. No r-process enhanced stars form if $t_{\rm min}=100$ Myr. Our results show that even when we adopt the r-process yield estimates observed in GW170817, NSMs by themselves can only explain the observed frequency of r-process enhanced stars if either the birth rate of double neutron stars per unit mass of stars is boosted to $\approx10^{-4} M_\odot^{-1}$ or the Europium yield of each NSM event is boosted to $\approx 10^{-4} M_{\odot}$.

The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network.

The recent aLIGO/aVirgo discovery of gravitational waves from the neutron star merger (NSM) GW170817 and the follow up kilonova observations have shown that NSMs produce copious amount of r-process material. However, it is difficult to reconcile the large natal kicks and long average merging times of Double Neutron Stars (DNSs), with the levels of r-process enrichment seen in ultra faint dwarf (UFD) galaxies such as Reticulum II and Tucana III. Assuming that such dwarf systems have lost a significant fraction of their stellar mass through tidal stripping, we conclude that contrary to most current models, it is the DNSs with rather large natal kicks but very short merging timescales that can enrich UFD-type galaxies. These binaries are either on highly eccentric orbits, or form with very short separations due to an additional mass-transfer between the first-born neutron star and a naked helium star, progenitor of the second-born neutron star. These DNSs are born with a frequency that agrees with the statistics of the r-process UFDs, and merge well within the virial radius of their host halos, therefore contributing significantly to their r-process enrichment.

Molecular clouds exhibit lognormal probability density functions (PDF) of mass densities, which are thought to arise as a consequence of isothermal, supersonic turbulence. Star formation is then widely assumed to occur in perturbations in which gravitational collapse is faster than the rate of change due to turbulent motions. Here we use direct numerical simulations to measure this rate as a function of density for a range of turbulent Mach numbers, and show it is faster at high densities than at low densities. Furthermore, we show that both the density PDF and rate of change arise naturally in a simple model of turbulence as a continuous Markov process. The one-dimensional Langevin equation that describes this evolution depends on only two parameters, yet it captures the full evolution seen in direct three-dimensional simulations. If it is modified to include gravity, the Langevin equation also reproduces the rate of material collapsing to high densities seen in turbulent simulations including self-gravity. When generalized to include both temperature and density, similar analyses are likely applicable throughout astrophysics.

Although it is assumed that magnetic fields in accretion disks are dragged from the interstellar medium, the idea is likely not applicable to primordial black holes (PBHs) formed in the early universe. Here we show that magnetic fields can be generated in initially unmagnetized accretion disks around PBHs through the Biermann battery mechanism, and therefore provide the small scale seeds of magnetic field in the universe. The radial temperature and vertical density profiles of these disks provide the necessary conditions for the battery to operate naturally. The generated seed fields have a toroidal structure with opposite sign in the upper and lower half of the disk. In the case of a thin accretion disk around a rotating PBH, the field generation rate increases with increasing PBH spin. At a fixed $r/r_{\rm isco}$, where $r$ is the radial distance from the PBH and $r_{\rm isco}$ is the radius of the innermost stable circular orbit, the battery scales as $M^{-9/4}$, where $M$ is the PBH's mass. The very weak dependency of the battery on accretion rate, makes this mechanism a viable candidate to provide seed fields in an initially unmagnetized accretion disk, following which the magnetorotational instability could take over.

The recent EDGES collaboration detection of an absorption signal at a central frequency of $\nu = 78 \pm 1$ MHz points to the presence of a significant Lyman-$\alpha$ background by a redshift of $z=18$. The timing of this signal constrains the dark matter particle mass ($m_\chi$) in the warm dark matter (WDM) cosmological model. WDM delays the formation of small-scale structures, and therefore a stringent lower limit can be placed on $m_\chi,$ based on the presence of a sufficiently strong Ly-$\alpha$ background due to star formation at $z=18$. Our results show that the coupling the spin temperature to the gas through Ly-$\alpha$ pumping requires a minimum mass of $m_\chi>3$ keV if atomic cooling halos dominate the star formation rate at $z=18,$ and $m_\chi>2$ keV if ${\rm\,H_2}$ cooling halos also form stars efficiently at this redshift. These limits match or exceed the most stringent limits cited to date in the literature, even in the face of the many uncertainties regarding star-formation at high redshift.

The smallest satellites of the Milky Way ceased forming stars during the epoch of reionization and thus provide archaeological access to galaxy formation at $z>6$. Numerical studies of these ultra-faint dwarf galaxies (UFDs) require expensive cosmological simulations with high mass resolution that are carried out down to $z=0$. However, if we are able to statistically identify UFD host progenitors at high redshifts \emph{with relatively high probabilities}, we can avoid this high computational cost. To find such candidates, we analyze the merger trees of Milky Way type halos from the high-resolution ${\it Caterpillar}$ suite of dark matter only simulations. Satellite UFD hosts at $z=0 $ are identified based on four different abundance matching (AM) techniques. All the halos at high redshifts are traced forward in time in order to compute the probability of surviving as satellite UFDs today. Our results show that selecting potential UFD progenitors based solely on their mass at z=12 (8) results in a 10\% (20\%) chance of obtaining a surviving UFD at $z=0$ in three of the AM techniques we adopted. We find that the progenitors of surviving satellite UFDs have lower virial ratios ($\eta$), and are preferentially located at large distances from the main MW progenitor, while they show no correlation with concentration parameter. Halos with favorable locations and virial ratios are $\approx 3$ times more likely to survive as satellite UFD candidates at $z=0.$

We investigate the stellar mass loss of gas rich galaxies falling into clusters due to the change in the gravitational potential caused by the ram pressure stripping of their gas. We model the satellites with exponential stellar and gas disk profiles, assume rapid ram pressure stripping, and follow the stellar orbits in the shocked potential. Due to the change of the potential, the stars move from circular orbits to elliptical orbits with apocenters that are often outside the tidal radius, causing those stars to be stripped. We explore the impact of the redshift of infall, gas fraction, satellite halo mass and cluster mass on this process. The puffing of the satellites makes them appear as ultra diffuse galaxies, and the stripped stars contribute to the intracluster light. Our results show these effects are most significant for less-massive satellites, which have larger gas fractions when they are accreted into clusters. The preferential destruction of low mass systems causes the red fraction of cluster galaxies to be smaller at lower masses, an observation that is otherwise difficult to explain.

Matching the number counts and redshift distribution of submillimeter galaxies (SMGs) without invoking modifications to the initial mass function (IMF) has proved challenging for semi-analytic models (SAMs) of galaxy formation. We adopt a previously developed SAM that is constrained to match the z=0 galaxy stellar mass function and makes various predictions that agree well with observational constraints; we do not recalibrate the SAM for this work. We implement three prescriptions to predict the submillimeter flux densities of the model galaxies; two depend solely on star formation rate, whereas the other also depends on the dust mass. By comparing the predictions of the models, we find that taking into account the dust mass, which affects the dust temperature and thus influences the far-infrared spectral energy distribution, is crucial for matching the number counts and redshift distribution of SMGs. Moreover, despite using a standard IMF, our model can match the observed SMG number counts and redshift distribution reasonably well, which contradicts the conclusions of some previous studies that a top-heavy IMF, in addition to taking into account the effect of dust mass, is needed to match these observations. Although we have not identified the key ingredient that is responsible for our model matching the observed SMG number counts and redshift distribution without IMF variation- which is challenging given the different prescriptions for physical processes employed in the SAMs of interest-our results demonstrate that in SAMs, IMF variation is degenerate with other physical processes, such as stellar feedback.

We study the impact of the neutron star binaries' (NSBs) natal kick distribution on the Galactic r-process enrichment. We model the growth of a Milky Way type halo based on N-body simulation results and its star formation history based on multi epoch abundance matching techniques. We consider the NSBs that merge well beyond the galaxy's effective radius ($>2\times R_\mathrm{eff}$) do not contribute to Galactic r-process enrichment. Assuming a power-law delay-time distribution (DTD) function ($\propto t^{-1}$) with $t_\mathrm{min}=30$ Myr for binaries' coalescence timescales, and an exponential profile for their natal kick distribution with an average value of 180 km s$^{-1}$, we show that up to $\sim$ 40% of all formed NSBs do not contribute to r-process enrichment by $z=0$, either because they merge far from the galaxy at a given redshift (up to $\sim$ 25%) or have not yet merged by today ($\sim$15%). Our result is largely insensitive to the details of the DTD function. Assuming a constant coalescence timescale of 100 Myr well approximates the adopted DTD with 30% of the NSBs not contributing to r-process enrichment. Our results, although rather dependent on the adopted natal kick distribution, represent a first step towards estimating the impact of natal kicks and DTD functions on r-process enrichment of galaxies that would need to be incorporated in the hydrodynamical simulations.

To explain the high observed abundances of r-process elements in local ultra- faint dwarf (UFD) galaxies, we perform cosmological zoom simulations that include r-process production from neutron star mergers (NSMs). We model star-formation stochastically and simulate two different halos with total masses $\approx 10^8 M_{\odot}$ at z = 6. We find that the final distribution of [Eu/H] vs. [Fe/H] is relatively insensitive to the energy by which the r-process material is ejected into the interstellar medium, but strongly sensitive to the environment in which the NSM event occurs. In one halo the NSM event takes place at the center of the stellar distribution, leading to high-levels of r-process enrichment such as seen in a local UFD, Reticulum II (Ret II). In a second halo, the NSM event takes place outside of the densest part of the galaxy, leading to a more extended r-process distribution. The subsequent star formation occurs in an interstellar medium with shallow levels of r-process enrichment which results in stars with low levels of [Eu/H] compared to Ret II stars even when the maximum possible r-process mass is assumed to be ejected. This suggests that the natal kicks of neutron stars may also play an important role in determining the r-process abundances in UFD galaxies, a topic that warrants further theoretical investigation.

We study the relationship between the UV continuum slope and infrared excess (IRX$\equiv L_{\rm IR}/L_{\rm FUV}$) predicted by performing dust radiative transfer on a suite of hydrodynamical simulations of galaxies. Our suite includes both isolated disk galaxies and mergers intended to be representative of galaxies at both $z \sim 0$ and $z \sim 2-3$. Our low-redshift isolated disks and mergers often populate a region around the the locally calibrated \citet[][M99]{M99} relation but move well above the relation during merger-induced starbursts. Our high-redshift simulated galaxies are blue and IR-luminous, which makes them lie above the M99 relation. The value of UV continuum slope strongly depends on the dust type used in the radiative transfer calculations: Milky Way-type dust leads to significantly more negative (bluer) slopes compared with Small Magellanic Cloud-type dust. The effect on $\beta$ due to variations in the dust composition with galaxy properties or redshift can dominate over other sources of $\beta$ variations and is the dominant model uncertainty. The dispersion in $\beta$ is anticorrelated with specific star formation rate and tends to be higher for the $z \sim 2-3$ simulations. In the actively star-forming $z \sim 2-3$ simulated galaxies, dust attenuation dominates the dispersion in $\beta$, whereas in the $z \sim 0$ simulations, the contributions of SFH variations and dust are similar. For low-SSFR systems at both redshifts, SFH variations dominate the dispersion. Finally, the simulated $z \sim 2-3$ isolated disks and mergers both occupy a region in the \irxbeta\ plane consistent with observed $z \sim 2-3$ dusty star-forming galaxies (DSFGs). Thus, contrary to some claims in the literature, the blue colors of high-z DSFGs do not imply that they are short-lived starbursts.

We investigate the relation between the turbulent Mach number (\mach) and the escape fraction of Lyman continuum photons ($f_{\rm esc}$) in high-redshift galaxies. Approximating the turbulence as isothermal and isotropic, we show that the increase in the variance in column densities from $\mathcal{M}=1$ to $\mathcal{M}=10$ causes $f_{\rm esc}$ to increase by $\approx 25$\%, and the increase from $\mathcal{M}=1$ to $\mathcal{M}=20$ causes $f_{\rm esc}$ to increases by $\approx 50$\% for a medium with opacity $\tau\approx1$. At a fixed Mach number, the correction factor for escape fraction relative to a constant column density case scales exponentially with the opacity in the cell, which has a large impact for simulated star forming regions. Furthermore, in simulations of isotropic turbulence with full atomic/ionic cooling and chemistry, the fraction of HI drops by a factor of $\approx 2.5$ at $\mathcal{M}\approx10$ even when the mean temperature is $\approx5\times10^3 K$. If turbulence is unresolved, these effects together enhance $f_{\rm esc}$ by a factor $>3$ at Mach numbers above 10. Such Mach numbers are common at high-redshifts where vigorous turbulence is driven by supernovae, gravitational instabilities, and merger activity, as shown both by numerical simulations and observations. These results, if implemented in the current hydrodynamical cosmological simulations to account for unresolved turbulence, can boost the theoretical predictions of the Lyman Continuum photon escape fraction and further constrain the sources of reionization.

Utilizing high-resolution cosmological hydrodynamic simulations we investigate various ultra-violet absorption lines in the circumgalactic medium of star forming galaxies at low redshift, in hopes of checking and alleviating the claimed observational conundrum of the ratio of NV to OVI absorbers, among others. We find a satisfactory agreement between simulations and extant observational data with respect to the ratios of the following four line pairs examined, NV/OVI, SiIV/OVI, NIII/OVI and NII/OVI. For the pairs involving nitrogen lines, we examine two cases of nitrogen abundance, one with constant N/O ratio and the other with varying N/O ratio, with the latter motivated by theoretical considerations of two different synthetic sources of nitrogen that is empirically verified independently. Along a separate vector, for all line pairs, we examine two cases of radiation field, one with the Haardt-Madau background radiation field and the other with an additional local radiation field sourced by hot gas in the host galaxy. In all cases, two-sample Kolmogorov-Smirnov tests indicate excellent agreements. We find that the apparent agreements between simulations and observations will be strongly tested, if the bulk of current upper limits of various line ratios are turned into actual detections. We show that an increase in observational sensitivity by 0.2 dex will already start to significantly constrain the models.

We perform a statistical analysis of strong gravitational lensing by quasar hosts of background galaxies, in the two competing models of dark matter halos of quasars, HOD and CS models. Utilizing the BolshoiP Simulation we demonstrate that strong gravitational lensing provides a potentially very powerful test of models of quasar hosting halos. For quasars at $z=0.5$, the lensing probability by quasars of background galaxies in the HOD model is higher than that of the CS model by two orders of magnitude or more for lensing image separations in the range of $\theta\sim 1.2-12~$arcsec. To observationally test this, we show that, as an example, at the depth of the CANDELS wide field survey and with a quasar sample of $1000$ at $z=0.5$, the two models can be differentiated at $3-4\sigma$ confidence level.

To explore the connection between the global physical properties of galaxies and their far-infrared (FIR) spectral energy distributions (SEDs), we study the variation in the FIR SEDs of a set of hydrodynamically simulated galaxies that are generated by performing dust radiative transfer in post-processing. Our sample includes both isolated and merging systems at various stages of the merging process and covers infrared (IR) luminosities and dust masses that are representative of both low- and high-redshift galaxies. We study the FIR SEDs using principle component analysis (PCA) and find that 97\% of the variance in the sample can be explained by two principle components (PCs). The first PC characterizes the wavelength of the peak of the FIR SED, and the second encodes the breadth of the SED. We find that the coefficients of both PCs can be predicted well using a double power law in terms of the IR luminosity and dust mass, which suggests that these two physical properties are the primary determinants of galaxies' FIR SED shapes. Incorporating galaxy sizes does not significantly improve our ability to predict the FIR SEDs. Our results suggest that the observed redshift evolution in the effective dust temperature at fixed IR luminosity is not driven by geometry: the SEDs of $z \sim 2-3$ ultraluminous IR galaxies (ULIRGs) are cooler than those of local ULIRGs not because the high-redshift galaxies are more extended but rather because they have higher dust masses at fixed IR luminosity. Finally, based on our simulations, we introduce a two-parameter set of SED templates that depend on both IR luminosity and dust mass.

A statistical analysis of stacked Compton$-y$ maps of quasar hosts with a median redshift of $1.5$ using Millennium Simulation is performed to address two issues, one on the feedback energy from quasars and the other on testing dark matter halo models for quasar hosts. On the first, we find that, at the resolution of FWHM=$10$ arcmin obtained by Planck data, the observed thermal Sunyaev-Zeldovich (tSZ) effect can be entirely accounted for and explained by the thermal energy of halos sourced by gravitational collapse of halos, without a need to invoke additional, large energy sources, such as quasar or stellar feedback. Allowing for uncertainties of dust temperature in the calibration of observed Comton$-y$ maps, the maximum additional feedback energy is $\sim 25\%$ of that previously suggested. Second, we show that, with FWHM=$1$ arcmin beam, tSZ measurements will provide a potentially powerful test of quasar-hosting dark matter halo models, limited only by possible observational systematic uncertainties, not by statistical ones, even in the presence of possible quasar feedback.

A new model for quasar-hosting dark matter halos, meeting two physical conditions, is put forth. First, significant interactions are taken into consideration to trigger quasar activities. Second, satellites in very massive halos at low redshift are removed from consideration, due to their deficiency of cold gas. We analyze the {\em Millennium Simulation} to find halos that meet these two conditions and simultaneously match two-point auto-correlation functions of quasars and cross-correlation functions between quasars and galaxies at $z=0.5-3.2$. %The found halos have some distinct properties worth noting. The masses of found quasar hosts decrease with decreasing redshift, with the mass thresholds being $[(2-5)\times 10^{12}, (2-5)\times 10^{11}, (1-3)\times 10^{11}]\msun$ for median luminosities of $\sim[10^{46}, 10^{46}, 10^{45}]$erg/s at $z=(3.2, 1.4, 0.53)$, respectively, an order of magnitude lower than those inferred based on halo occupation distribution modeling. In this model quasar hosts are primarily massive central halos at $z\ge 2-3$ but increasingly dominated by lower mass satellite halos experiencing major interactions towards lower redshift. But below $z=1$ satellite halos in groups more massive than $\sim 2\times 10^{13}\msun$ do not host quasars. Whether for central or satellite halos, imposing the condition of significant interactions substantially boosts the clustering strength compared to the total population with the same mass cut. The inferred lifetimes of quasars at $z=0.5-3.2$ of $3-30$Myr are in agreement with observations. Quasars at $z\sim 2$ would be hosted by halos of mass $\sim 5\times 10^{11}\msun$ in this model, compared to $\sim 3\times 10^{12}\msun$ previously thought, which would help reconcile with the observed, otherwise puzzling high covering fractions for Lyman limit systems around quasars.

We present a new technique for overcoming confusion noise in deep far-infrared \Herschel space telescope images making use of prior information from shorter $\lambda<2$\micron wavelengths. For the deepest images obtained by \Herschels, the flux limit due to source confusion is about a factor of three brighter than the flux limit due to instrumental noise and (smooth) sky background. We have investigated the possibility of de-confusing simulated \Herschel PACS-160\micron images by using strong Bayesian priors on the positions and weak priors on the flux of sources. We find the blended sources and group them together and simultaneously fit their fluxes. We derive the posterior probability distribution function of fluxes subject to these priors through Monte Carlo Markov Chain (MCMC) sampling by fitting the image. Assuming we can predict FIR flux of sources based on ultraviolet-optical part of their SEDs to within an order of magnitude, the simulations show that we can obtain reliable fluxes and uncertainties at least a factor of three fainter than the confusion noise limit of $3\sigma_{c} $=2.7 mJy in our simulated PACS-160 image. This technique could in principle be used to mitigate the effects of source confusion in any situation where one has prior information of positions and plausible fluxes of blended sources. For \Herschel, application of this technique will improve our ability to constrain the dust content in normal galaxies at high redshift.

We compare the predictions of three independently developed semi-analytic galaxy formation models that are being used to aid in the interpretation of results from the CANDELS survey. These models are each applied to the same set of halo merger trees extracted from the "Bolshoi" simulation and are carefully tuned to match the local galaxy stellar mass function using the powerful method of Bayesian Inference coupled with MCMC or by hand. The comparisons reveal that in spite of the significantly different parameterizations for star formation and feedback processes, the three models yield qualitatively similar predictions for the assembly histories of galaxy stellar mass and star formation over cosmic time. We show that the SAMs generally require strong outflows to suppress star formation in low-mass halos to match the present day stellar mass function. However, all of the models considered produce predictions for the star formation rates and metallicities of low-mass galaxies that are inconsistent with existing data and diverge between the models. We suggest that large differences in the metallicity relations and small differences in the stellar mass assembly histories of model galaxies stem from different assumptions for the outflow mass-loading factor. Importantly, while more accurate observational measurements for stellar mass, SFR and metallicity of galaxies at 1<z<5 will discriminate between models, the discrepancies between the models and existing data of these observables have already revealed challenging problems in understanding star formation and its feedback in galaxy formation. The three sets of models are being used to construct catalogs of mock galaxies on light cones that have the same geometry as the CANDELS survey, which should be particularly useful for quantifying the biases and uncertainties on measurements and inferences from the real observations. -ABRIDGED

The rate of ionization by cosmic rays in interstellar gas directly associated with gamma-ray emitting supernova remnants is for the first time calculated to be several orders of magnitude larger than the Galactic average. Analysis of ionization-induced chemistry yields the first quantitative prediction of the astrophysical H_{2}^{+} emission line spectrum, which should be detectable together with H_3^+ lines. The predicted coincident observation of those emission lines and gamma-rays will help prove that supernova remnants are sources of cosmic rays.

We study the process of spinning up black holes by accretion from slim disks in a wide range of accretion rates. We show that for super-Eddington accretion rates and low values of the viscosity parameter alpha (\lesssim 0.01) the limiting value of the dimensionless spin parameter a_* can reach values higher than a_*=0.9978 inferred by Thorne (1974) in his seminal study. For \mdot=10\mdot_{Edd} and alpha=0.01 spin equilibrium is reached at a_*=0.9994. We show that the equilibrium spin value depends strongly on the assumed value of alpha. We also prove that for high accretion rates the impact of captured radiation on spin evolution is negligible.