We study the evolution of cosmological perturbations in a contracting universe. We aim to determine under which conditions density perturbations grow to form large inhomogeneities and collapse into black holes. Our method consists in solving the cosmological perturbation equations in complete generality for a hydrodynamical fluid. We then describe the evolution of the fluctuations over the different length scales of interest and as a function of the equation of state for the fluid, and we explore two different types of initial conditions: quantum vacuum and thermal fluctuations. We also derive a general requirement for black hole collapse on sub-Hubble scales, and we use the Press-Schechter formalism to describe the black hole formation probability. For a fluid with a small sound speed (e.g., dust), we find that both quantum and thermal initial fluctuations grow in a contracting universe, and the largest inhomogeneities that first collapse into black holes are of Hubble size and the collapse occurs well before reaching the Planck scale. For a radiation-dominated fluid, we find that no black hole can form before reaching the Planck scale. In the context of matter bounce cosmology, it thus appears that only models in which a radiation-dominated era begins early in the cosmological evolution are robust against the formation of black holes. Yet, the formation of black holes might be an interesting feature for other models. We comment on a number of possible alternative early universe scenarios that could take advantage of this feature.
We propose to use the flux variability of lensed quasar images induced by gravitational microlensing to measure the transverse peculiar velocity of lens galaxies over a wide range of redshift. Microlensing variability is caused by the motions of the observer, the lens galaxy (including the motion of the stars within the galaxy), and the source; hence, its frequency is directly related to the galaxy's transverse peculiar velocity. The idea is to count time-event rates (e.g., peak or caustic crossing rates) in the observed microlensing light curves of lensed quasars that can be compared with model predictions for different values of the transverse peculiar velocity. To compensate for the large time-scale of microlensing variability we propose to count and model the number of events in an ensemble of gravitational lenses. We develop the methodology to achieve this goal and apply it to an ensemble of 17 lensed quasar systems. In spite of the shortcomings of the available data, we have obtained tentative estimates of the peculiar velocity dispersion of lens galaxies at $z\sim 0.5$, $\sigma_{\rm pec}(0.53\pm0.18)\simeq(638\pm213)\sqrt{\langle m \rangle/0.3 M_\odot} \, \rm km\, s^{-1}$. Scaling at zero redshift we derive, $\sigma_{\rm pec}(0)\simeq(491\pm164) \sqrt{\langle m \rangle/0.3 M_\odot} \, \rm km\, s^{-1}$, consistent with peculiar motions of nearby galaxies and with recent $N$-body nonlinear reconstructions of the Local Universe based on $\Lambda$CDM. We analyze the different sources of uncertainty of the method and find that for the present ensemble of 17 lensed systems the error is dominated by Poissonian noise, but that for larger ensembles the impact of the uncertainty on the average stellar mass may be significant.
We perform a detailed analysis on a primordial gravitational-wave background amplified during a Kasner-like pre-inflationary phase allowing for general triaxial anisotropies. It is found that the predicted angular distribution map of gravitational-wave intensity on large scales exhibits topologically distinctive patterns according to the degree of the pre-inflationary anisotropy, thereby serving as a potential probe for the pre-inflationary early universe with future all-sky observations of gravitational waves. We also derive an observational limit on the amplitude of such anisotropic gravitational waves from the B-mode polarisation of the cosmic microwave background.
We update the ingredients of the Gaussian streaming model (GSM) for the redshift-space clustering of biased tracers using the techniques of Lagrangian perturbation theory, effective field theory (EFT) and a generalized Lagrangian bias expansion. After relating the GSM to the cumulant expansion, we present new results for the real-space correlation function, mean pairwise velocity and pairwise velocity dispersion including counter terms from EFT and bias terms through third order in the linear density, its leading derivatives and its shear up to second order. We discuss the connection to the Gaussian peaks formalism. We compare the ingredients of the GSM to a suite of large N-body simulations, and show the performance of the theory on the low order multipoles of the redshift-space correlation function and power spectrum. We highlight the importance of a general biasing scheme, which we find to be as important as higher-order corrections due to non-linear evolution for the halos we consider on the scales of interest to us.
We consider f(R,T) modified theory of gravity in which, in general, the gravitational Lagrangian is given by an arbitrary function of the Ricci scalar and the trace of the energy-momentum tensor. We indicate that in this type of the theory, the coupling energy-momentum tensor is~not conserved. However, we mainly focus on a particular model that matter is minimally coupled to the geometry in the metric formalism and wherein, its coupling energy--momentum tensor is also conserved. We obtain the corresponding Raychaudhuri dynamical equation that presents the evolution of the kinematic quantities. Then for the chosen model, we derive the behavior of the deceleration parameter, and show that the coupling term can lead to an acceleration phase after the matter dominated phase. On the other hand, the curvature of the universe corresponds with the deviation from parallelism in the geodesic motion. Thus, we also scrutinize the motion of the free test particles on their geodesics, and derive the geodesic deviation equation in this modified theory to study the accelerating universe within the spatially flat FLRW background. Actually, this equation gives the relative accelerations of adjacent particles as a measurable physical quantity, and provides an elegant tool to investigate the timelike and the null structures of spacetime geometries. Then, through the null deviation vector, we find the observer area-distance as a function of the redshift for the chosen model, and compare the results with the corresponding results obtained in the literature.
We consider f(R,T) modified theory of gravity, in which the gravitational Lagrangian is given by an arbitrary function of the Ricci scalar and the trace of the energy-momentum tensor of the matter, in order to investigate the dark matter effects on the galaxy scale. We obtain the metric components for a spherically symmetric and static spacetime in the vicinity of general relativity solutions. However, we concentrate on a specific model of the theory where the matter is minimally coupled to the geometry, and derive the metric components in the galactic halo. Then, we fix the components by the rotational velocities of galaxies for the model, and show that the mass corresponding to the interaction term (which appears in the Einstein modified field equation) leads to a flat rotation curve in the halo of galaxies. Also for the proposed model, the light deflection angle has been derived, and then, been drawn while using some observed data.
We consider the fact that noticing on the operational meaning of the physical concepts played an impetus role in the appearance of general relativity (GR). Thus, we have paid more attention to the operational definition of the gravitational coupling constant in this theory as a dimensional constant which is gained through an experiment. However, as all available experiments just provide the value of this constant locally, this coupling constant can operationally be meaningful only in a local area. Regarding this point, to obtain an extension of GR for the large scale, we replace it by a conformal invariant model and then, reduce this model to a theory for the cosmological scale via breaking down the conformal symmetry through singling out a specific conformal frame which is characterized by the large scale characteristics of the universe. Finally, we come to the same field equations that historically were proposed by Einstein for the cosmological scale (GR plus the cosmological constant) as the result of his endeavor for making GR consistent with the Mach principle. However, we declare that the obtained field equations in this alternative approach do~not carry the problem of the field equations proposed by Einstein for being consistent with Mach's principle (i.e., the existence of de Sitter solution), and can also be considered compatible with this principle in the Sciama view.
It has recently been shown that if the dark matter is in thermal equilibrium with a sector that is highly decoupled from the Standard Model, it can freeze-out with an acceptable relic abundance, even if the dark matter is as heavy as ~1-100 PeV. In such scenarios, both the dark and visible sectors are populated after inflation, but with independent temperatures. The lightest particle in the dark sector will be generically long-lived, and can come to dominate the energy density of the universe. Upon decaying, these particles can significantly reheat the visible sector, diluting the abundance of dark matter and thus allowing for dark matter particles that are much heavier than conventional WIMPs. In this paper, we present a systematic and pedagogical treatment of the cosmological history in this class of models, emphasizing the simplest scenarios in which a dark matter candidate annihilates into hidden sector particles which then decay into visible matter through the vector, Higgs, or lepton portals. In each case, we find ample parameter space in which very heavy dark matter particles can provide an acceptable thermal relic abundance. We also discuss possible extensions of models featuring these dynamics.
We, in the present paper, investigate the formation and evolution of primordial black hole (PBH) within the scenario of nonsingular bouncing cosmology. We first analyze the PBH formation during the phase of matter contraction, which is different from that in an expanding background, and then evaluate the PBH abundance at the end of the contracting phase. Our result shows that it is generally small unless the energy scale parameter associated with the bouncing phase is as high as the Planck scale, i.e., $|H_{-}|\gtrsim M_p$, or the sound speed parameter of cosmological perturbations is sufficiently small, which implies, $c_s \ll 1$. Afterwards, we study the subsequent evolution of generating PBHs during the bouncing phase. For the PBH growth ignoring the Hawking radiation, a relation upon model parameters of the bouncing phase $\Upsilon \geq c_s^2 \pi^2 H^2_{-}$ is expected to be satisfied, in case that PBHs would grow to infinity before the bouncing point. We also calculate the back-reaction of PBHs in order to theoretically constrain bounce cosmologies by considering the effects of both PBH growth and the associated Hawking radiation. The constraint is in accordance to the relation $\Upsilon \geq c_s^2 \pi^2 H^2_{-}$ for the bounce cosmology with a relatively low energy scale $H_{-}^2 \ll 10^9 c_s^5 M_p^2$.
We report a study on axially and reflection symmetric dissipative fluids, just after its departure from hydrostatic and thermal equilibrium, at the smallest time scale at which the first signs of dynamic evolution appear. Such a time scale is smaller than the thermal relaxation time, the thermal adjustment time and the hydrostatic time. It is obtained that the onset of non--equilibrium will critically depend on a single function directly related to the time derivative of the vorticity. Among all fluid variables (at the time scale under consideration), only the tetrad component of the anisotropic tensor in the subspace orthogonal to the four--velocity and the Killing vector of axial symmetry, shows signs of dynamic evolution. Also, the first step towards a dissipative regime begins with a non--vanishing time derivative of the heat flux component along the meridional direction. The magnetic part of the Weyl tensor vanishes (not so its time derivative), indicating that the emission of gravitational radiation will occur at later times. Finally, the decreasing of the effective inertial mass density, associated to thermal effects, is clearly illustrated.
Links to: arXiv, form interface, find, astro-ph, recent, 1609, contact, help (Access key information)
The effects of near neighbors on the galaxy-galaxy lensing signal are investigated using a suite of Monte Carlo simulations. The redshifts, luminosities, and relative coordinates for the simulated lenses were obtained from a set of galaxies with known spectroscopic redshifts and known luminosities. As expected, when all lenses are assigned a single, fixed redshift, the mean tangential shear is identically equal to the excess surface mass density, scaled by the critical surface mass density: $\gamma_T = \Delta\Sigma \times \Sigma_c^{-1}$. When the lenses are assigned their observed redshifts and $\Sigma_c$ is taken to be the critical surface mass density of the central lens, the relationship $\gamma_T = \Delta\Sigma \times \Sigma_c^{-1}$ is violated because $\gtrsim 90$% of the near neighbors are located at redshifts significantly different from the central lenses. For a given central lens, physically unrelated near neighbors give rise to a ratio of $\gamma_T$ to $\Delta\Sigma \times \Sigma_c^{-1}$ that spans a wide range of $\sim 0.5$ to $\sim 1.5$ at projected distances $r_p \sim 1$ Mpc. The magnitude and sense of the discrepancy between $\gamma_T$ and $\Delta\Sigma \times \Sigma_c^{-1}$ are functions of both $r_p$ and the velocity dispersions of the central lenses, $\sigma_v$. At large $r_p$, the difference between $\gamma_T$ and $\Delta\Sigma \times \Sigma_c^{-1}$ is, on average, much greater for low-$\sigma_v$ central lenses than it is for high-$\sigma_v$ central lenses.
We use cosmological hydrodynamical simulations to study the effect of screened modified gravity models on the mass estimates of galaxy clusters. In particular, we focus on two novel aspects: (i) we study modified gravity models in which baryons and dark matter are coupled with different strengths to the scalar field, and, (ii) we put the simulation results into the greater context of a general screened-modified gravity parametrization. We compare the mass of clusters inferred via lensing versus the mass inferred via kinematical measurements as a probe of violations of the equivalence principle at Mpc scales. We find that estimates of cluster masses via X-ray observations is mainly sensitive to the coupling between the scalar degree of freedom and baryons -- while the kinematical mass is mainly sensitive to the coupling to dark matter. Therefore, the relation between the two mass estimates is a probe of a possible non-universal coupling between the scalar field, the standard model fields, and dark matter. Finally, we use observational data of kinetic, thermal and lensing masses to place constraints on deviations from general relativity on cluster scales for a general parametrization of screened modified gravity theories which contains $f(R)$ and Symmetron models. We find that while the kinematic mass can be used to place competitive constraints, using thermal measurements is challenging as a potential non-thermal contribution is degenerate with the imprint of modified gravity.
We show that, if they exist, lepton number asymmetries ($L_\alpha$) of neutrino flavors should be distinguished from the ones ($L_i$) of mass eigenstates, since Big Bang Nucleosynthesis (BBN) bounds on the flavor eigenstates cannot be directly applied to the mass eigenstates. Similarly, Cosmic Microwave Background (CMB) constraints on mass eigenstates do not directly constrain flavor asymmetries. Due to the difference of mass and flavor eigenstates, the cosmological constraint on the asymmetries of neutrino flavors can be much stronger than conventional expectation, but not uniquely determined unless at least the asymmetry of the heaviest neutrino is well constrained. Cosmological constraint on $L_i$ for a specific case is presented as an illustration.
Improvements in the accuracy of shape measurements are essential to exploit the statistical power of planned imaging surveys that aim to constrain cosmological parameters using weak lensing by large-scale structure. Although a range of tests can be performed using the measurements, the performance of the algorithm can only be quantified using simulated images. This yields, however, only meaningful results if the simulated images resemble the real observations sufficiently well. In this paper we explore the sensitivity of the multiplicative bias to the input parameters of Euclid-like image simulations.We find that algorithms will need to account for the local density of sources. In particular the impact of galaxies below the detection limit warrants further study, because magnification changes their number density, resulting in correlations between the lensing signal and multiplicative bias. Although achieving sub-percent accuracy will require further study, we estimate that sufficient archival Hubble Space Telescope data are available to create realistic populations of galaxies.
We present cosmic microwave background (CMB) power spectra from recent numerical simulations of cosmic strings in the Abelian Higgs model and compare them to CMB power spectra measured by Planck. We obtain revised constraints on the cosmic string tension parameter $G\mu$. For example, in the $\Lambda$CDM model with the addition of strings and no primordial tensor perturbations, we find $G\mu < 2.0 \times 10^{-7}$ at 95% confidence, about 20% lower than the value obtained from previous simulations, which had 1/64 of the spatial volume. We investigate the source of the difference, showing that the main cause is an improved treatment of the string evolution across the radiation-matter transition. The increased computational volume also makes possible to simulate fully the physical equations of motion, in which the string cores shrink in comoving coordinates. This, and the larger dynamic range, changes the amplitude of the power spectra by only about 10%, demonstrating that field theory simulations of cosmic strings have now reached the required dynamic range for CMB calculations.
Treating dark matter at large scales as an effectively viscous fluid provides an improved framework for the calculation of the density and velocity power spectra compared to the standard assumption of an ideal pressureless fluid. We discuss how this framework can be made concrete through an appropriate coarse-graining procedure. We also review results that demonstrate that it improves the convergence of cosmological perturbation theory.
Dwarf spheroidals are low-luminosity satellite galaxies of the Milky Way highly dominated by dark matter. Therefore, they are prime targets to search for signals from dark matter annihilation using gamma-ray observations. We analyse about 7 years of PASS8 Fermi data for seven classical dwarf galaxies, including Draco, adopting both the widely used Navarro-Frenk-White (NFW) profile and observationally motivated axisymmetric density profiles. For four of the selected dwarfs (Sextans, Carina, Sculptor and Fornax) axisymmetric mass models suggest a cored density profile rather than the commonly adopted cusped profile. We found that upper limits on the annihilation cross section for some of these dwarfs are significantly higher than the ones achieved using an NFW profile. Therefore, upper limits in the literature obtained using cusped profiles like the NFW might have been overestimated. Our results eventually show that it is extremely important to use observationally motivated density profiles going beyond the usually adopted NFW in order to obtain accurate constraints on the dark matter annihilation cross section.
We illustrate how it is possible to calculate the quantum gravitational effects on the spectra of primordial scalar/tensor perturbations starting from the canonical, Wheeler-De Witt, approach to quantum cosmology. The composite matter-gravity system is analysed through a Born-Oppenheimer approach in which gravitation is associated with the heavy degrees of freedom and matter (here represented by a scalar field) with the light ones. Once the independent degrees of freedom are identified the system is canonically quantised. The differential equation governing the dynamics of the primordial spectra with its quantum-gravitational corrections is then obtained and is applied to diverse inflationary evolutions. Finally, the analytical results are compared to observations through a Monte Carlo Markov Chain technique and an estimate of the free parameters of our approach is finally presented and the results obtained are compared with previous ones.
Cosmological simulations suggest a strong correlation between high optical-depth Ly$\alpha$ absorbers, which arise from the intergalactic medium (IGM), and 3-D mass overdensities on scales of $10-30$ $h^{-1}$ comoving Mpc. By examining the absorption spectra of $\sim$ 80,000 QSO sight-lines over a volume of 0.1 Gpc$^3$ in the Sloan Digital Sky Survey III (SDSS-III), we have identified an extreme overdensity, BOSS1441, which contains a rare group of strong Ly$\alpha$ absorbers at $z=2.32\pm 0.02$. This absorber group is associated with six QSOs at the same redshift on a 30 comoving Mpc scale. Using Mayall/MOSAIC narrowband and broadband imaging, we detect Ly$\alpha$ emitters (LAEs) down to $0.7\times L_{\rm{Ly\alpha}}^*$, and reveal a large-scale structure of Ly$\alpha$ emitters (LAEs) in this field. Our follow-up Large Binocular Telescope (LBT) observations have spectroscopically confirmed 19 galaxies in the density peak. We show that BOSS1441 has an LAE overdensity of $10.8\pm 2.6$ on a 15 comoving Mpc scale which could collapse to a massive cluster with $M\gtrsim10^{15}$ M$_\odot$ at $z\sim0$. This overdensity is among the most massive large-scale structures at $z\sim2$ discovered to date.
Most types of supernovae (SNe) have yet to be connected with their progenitor stellar systems. Here, we re-analyze the ten-year 1998-2008 SN sample collected by the Lick Observatory Supernova Search (LOSS) in order to constrain the progenitors of SNe Ia and stripped-envelope SNe (SE SNe; i.e., SNe IIb, Ib, Ic, and broad-lined Ic). We matched the LOSS galaxy sample with spectroscopy from the Sloan Digital Sky Survey and measured SN rates as a function of galaxy stellar mass, specific star-formation rate (sSFR), and oxygen abundance (metallicity). We find significant correlations between the SN rates and all three galaxy properties. The SN Ia correlations are consistent with other measurements, as well as with our previous explanation of these measurements in the form of a combination of the SN Ia delay-time distribution and the correlation between galaxy mass and age. Intriguingly, we measure a deficiency in the SE SN rates, relative to the SN II rates, in galaxies with low stellar masses, high sSFR values, and low metallicities. Using well-known galaxy scaling relations, any correlation between the rates and one of the galaxy properties examined here can be expressed as a correlation with the other two. These redundant correlations preclude us from establishing causality, i.e., from ascertaining which of the galaxy properties (or their combination) is the physical driver for the difference between the SE SN and SN II rates.
In Paper I of this series, we showed that the ratio between stripped-envelope supernova (SE SN) and Type II SN rates reveals a significant SE SN deficiency in galaxies with stellar masses $\lesssim 10^{10}~{\rm M}_\odot$. Here, we test this result by splitting the volume-limited subsample of the Lick Observatory Supernova Search (LOSS) SN sample into low- and high-mass galaxies and comparing the relative rates of various SN types found in them. The LOSS volume-limited sample contains 180 SNe and SN impostors and is complete for SNe Ia out to 80 Mpc and core-collapse SNe out to 60 Mpc. All of these transients were recently reclassified by us in Shivvers et al. (in prep.) We find that the relative rates of some types of SNe differ between low- and high-mass galaxies: SNe Ib and Ic are underrepresented by a factor of ~3 in low-mass galaxies. These galaxies also contain the only examples of SN 1987A-like SNe in the sample and host ~9 times as many SN impostors. Normal SNe Ia are ~30% more common in low-mass galaxies, making these galaxies better sources for homogeneous SN Ia cosmology samples. The relative rates of SNe IIb are consistent in both low- and high-mass galaxies. The same is true for broad-line SNe Ic, though our sample includes only two such objects. The results presented here strongly disagree with a similar analysis from the Palomar Transient Factory, as presented by Arcavi et al. (2010), which we ascribe to the incompleteness of their SN sample.
Fluctuations of the surface brightness of cosmic X-ray background (CXB) carry unique information about faint and low luminosity source populations, which is inaccessible for conventional large-scale structure (LSS) studies based on resolved sources. We used Chandra data of the XBOOTES field ($\sim9\,\mathrm{deg^2}$) to conduct the most accurate measurement to date of the power spectrum of fluctuations of the unresolved CXB on the angular scales of $\sim3\,$arcsec $-$ $\sim17\,$arcmin. We find that at sub-arcmin angular scales, the power spectrum is consistent with the AGN shot noise, without much need for any significant contribution from their one-halo term. This is consistent with the theoretical expectation that low-luminosity AGN reside alone in their dark matter halos. However, at larger angular scales we detect a significant LSS signal above the AGN shot noise. Its power spectrum, obtained after subtracting the AGN shot noise, follows a power law with the slope of $-0.8\pm0.1$ and its amplitude is much larger than what can be plausibly explained by the two-halo term of AGN. We demonstrate that the detected LSS signal is produced by unresolved clusters and groups of galaxies. For the flux limit of the XBOOTES survey, their flux-weighted mean redshift equals $\left<z\right>\sim0.3$, and the mean temperature of their intracluster medium (ICM), $\left<T\right>\approx 1.4$ keV, corresponds to the mass of $M_{500} \sim 10^{13.5}\,\mathrm{M}_\odot$. The power spectrum of CXB fluctuations carries information about the redshift distribution of these objects and the spatial structure of their ICM on the linear scales of up to $\sim$Mpc, i.e. of the order of the virial radius.
NGC 1365 is a Seyfert 2 galaxy with a starburst ring in its nuclear region. In this work we look at the XMM Reflection Grating Spectrometer (RGS) data from four 2012-13, three 2007 and two 2004 observations of NGC 1365. We characterise the narrow-line emitting gas visible by XMM RGS and make comparisons between the 2012-13 spectra and those from 2004-07, already published. This source is usually absorbed within the soft X-ray band, with a typical neutral column density of >1.5 x 10$^{23}$ cm$^{-2}$, and only 1 observation of the 9 we investigate shows low enough absorption for the continuum to emerge in the soft X-rays. We stack all observations from 2004-07, and separately three of the four observations from 2012-13, analysing the less absorbed observation separately. We first model the spectra using gaussian profiles representing the narrow line emission. We fit physically motivated models to the 2012-13 stacked spectra, with collisionally ionised components representing the starburst emission and photoionised line emission models representing the AGN line emission. The collisional and photoionised emission line models are fitted together on top of a physical continuum and absorption model. The X-ray narrow emission line spectrum of NGC 1365 is well represented by a combination of two collisionally ionised and three photoionised phases of emitting gas, all with higher than solar nitrogen abundances. This physical model was fitted to the 2012-13 stacked spectrum, and yet also fits well to the 2004-07 stacked spectrum, without changing any characteristics of the emitting gas phases. Our 2004-07 results are consistent with previous emission line work using these data, with 5 additional emission lines detected in both this and the 2012-13 stacked spectra. We also estimate the distance of the X-ray line-emitting photoionised gas from the central source to be <300 pc.
Virial shocks at edges of cosmic-web structures are a clear prediction of
standard structure formation theories. We derive a criterion for the stability
of the post-shock gas and of the virial shock itself in spherical, filamentary
and planar infall geometries. When gas cooling is important, we find that
shocks become unstable, and gas flows uninterrupted towards the center of the
respective halo, filament or sheet. For filaments, we impose this criterion on
self-similar infall solutions. We find that instability is expected for
filament masses between $10^{11}-10^{13}M_\odot Mpc^{-1}.$ Using a simplified
toy model, we then show that these filaments will likely feed halos with
$10^{10}M_{\odot}\lesssim M_{halo}\lesssim 10^{13}M_{\odot}$ at redshift $z=3$,
as well as $10^{12}M_{\odot}\lesssim M_{halo}\lesssim 10^{15}M_{\odot}$ at
$z=0$.
The instability will affect the survivability of the filaments as they
penetrate gaseous halos in a non-trivial way. Additionally, smaller halos
accreting onto non-stable filaments will not be subject to ram-pressure inside
the filaments. The instreaming gas will continue towards the center, and stop
either once its angular momentum balances the gravitational attraction, or when
its density becomes so high that it becomes self-shielded to radiation.
Scale invariant fluctuations of metric are universal feature of quantum gravity in de Sitter spacetime. We construct an effective Lagrangian which summarizes their implications on local physics by integrating super-horizon metric fluctuations. It shows infrared quantum effects are local and render fundamental couplings time dependent. We impose Lorenz invariance on the effective Lagrangian as it is required by the principle of general covariance. We show that such a requirement leads to unique physical predictions by fixing the quantization ambiguities. We explain how the gauge parameter dependence of observables is canceled. In particular the relative evolution speed of the couplings are shown to be gauge invariant.
We carry out a study of the exterior of an axially and reflection symmetric source of gravitational radiation. The exterior of such a source is filled with a null fluid produced by the dissipative processes inherent to the emission of gravitational radiation, thereby representing a generalization of the Vaidya metric for axially and reflection symmetric spacetimes. The role of the vorticity, and its relationship with the presence of gravitational radiation is put in evidence. The spherically symmetric case (Vaidya) is, asymptotically, recovered within the context of the $1+3$ formalism.
We test whether halo age and galaxy age are correlated at fixed halo and galaxy mass. The formation histories, and thus ages, of dark matter halos correlate with their large-scale density $\rho$, an effect known as assembly bias. We test whether this correlation extends to galaxies by measuring the dependence of galaxy stellar age on $\rho$. To clarify the comparison between theory and observation, and to remove the strong environmental effects on satellites, we use galaxy group catalogs to identify central galaxies and measure their quenched fraction, $f_Q$, as a function of large-scale environment. Models that match halo age to central galaxy age predict a strong positive correlation between $f_Q$ and $\rho$. However, we show that the amplitude of this effect depends on the definition of halo age: assembly bias is significantly reduced when removing the effects of splashback halos---those halos that are central but have passed through a larger halo or experienced strong tidal encounters. Defining age using halo mass at its peak value rather than current mass removes these effects. In SDSS data, at M$_{\rm gal}\gtrsim 10^{10.0}$ M_sol/h$^2$, there is a $\sim 5\%$ increase in $f_Q$ from low to high densities, which is in agreement with predictions of dark matter halos using peak halo mass. At lower stellar mass there is little to no correlation of $f_Q$ with $\rho$. For these galaxies, age-matching is inconsistent with the data across the wide range the halo formation metrics that we tested. This implies that halo formation history has a small but statistically significant impact on quenching of star formation at high masses, while the quenching process in low-mass central galaxies is uncorrelated with halo formation history.
I will talk on our new theory on baryogenesis through type-II leptogenesis which is different from the well-known type-I leptogenesis. I will comment on the Jarlskog phases, $\delta_{\rm CKM}$ and $\delta_{\rm PMNS}$, in the CKM and PMNS matrices. In the type-II leptogenesis, the PMNS phase is used for Sakharov's condition on the global quantum number generation in the Universe. For this to be effective, the SU(2)$\times$U(1) gauge symmetry must be broken during the leptogenesis epoch.
Links to: arXiv, form interface, find, astro-ph, recent, 1609, contact, help (Access key information)
Turbulence is a key ingredient for the evolution of the intracluster medium, whose properties can be predicted with high resolution numerical simulations. We present initial results on the generation of solenoidal and compressive turbulence in the intracluster medium during the formation of a small-size cluster using highly resolved, non-radiative cosmological simulations, with a refined monitoring in time. In this first of a series of papers, we closely look at one simulated cluster whose formation was distinguished by a merger around $z \sim 0.3$. We separate laminar gas motions, turbulence and shocks with dedicated filtering strategies and distinguish the solenoidal and compressive components of the gas flows using Hodge-Helmholtz decomposition. Solenoidal turbulence dominates the dissipation of turbulent motions ($\sim 95\%$) in the central cluster volume at all epochs. The dissipation via compressive modes is found to be more important ($\sim 30 \%$ of the total) only at large radii ($\geq 0.5 ~r_{\rm vir}$) and close to merger events. We show that enstrophy (vorticity squared) is good proxy of solenoidal turbulence. All terms ruling the evolution of enstrophy (i.e. baroclinic, compressive, stretching and advective terms) are found to be significant, but in amounts that vary with time and location. Two important trends for the growth of enstrophy in our simulation are identified: first, enstrophy is continuously accreted into the cluster from the outside, and most of that accreted enstrophy is generated near the outer accretion shocks by baroclinic and compressive processes. Second, in the cluster interior vortex stretching is dominant, although the other terms also contribute substantially.
We present an improved photometric redshift estimator code, CuBAN$z$, that is publicly available at https://goo.gl/fpk90V}{https://goo.gl/fpk90V. It uses the back propagation neural network along with clustering of the training set, which makes it more efficient than existing neural network codes. In CuBAN$z$, the training set is divided into several self learning clusters with galaxies having similar photometric properties and spectroscopic redshifts within a given span. The clustering algorithm uses the color information (i.e. $u-g$, $g-r$ etc.) rather than the apparent magnitudes at various photometric bands as the photometric redshift is more sensitive to the flux differences between different bands rather than the actual values. Separate neural networks are trained for each cluster using all possible colors, magnitudes and uncertainties in the measurements. For a galaxy with unknown redshift, we identify the closest possible clusters having similar photometric properties and use those clusters to get the photometric redshifts using the particular networks that were trained using those cluster members. For galaxies that do not match with any training cluster, the photometric redshifts are obtained from a separate network that uses entire training set. This clustering method enables us to determine the redshifts more accurately. SDSS Stripe 82 catalog has been used here for the demonstration of the code. For the clustered sources with redshift range $z_{\rm spec}<0.7$, the residual error ($\langle (z_{{\rm spec}}-z_{{\rm phot}})^2\rangle^{1/2} $) in the training/testing phase is as low as 0.03 compared to the existing ANNz code that provides residual error on the same test data set of 0.05. Further, we provide a much better estimate of the uncertainty of the derived photometric redshift.
The standard paradigm of collisionless cold dark matter is in tension with measurements on large scales. In particular, the best fit values of the Hubble rate $H_0$ and the matter density perturbation $\sigma_8$ inferred from the cosmic microwave background seem inconsistent with the results from direct measurements. We show that both problems can be solved in a framework in which dark matter consists of two distinct components, a dominant component and a subdominant component. The primary component is cold and collisionless. The secondary component is also cold, but interacts strongly with dark radiation, which itself forms a tightly coupled fluid. The growth of density perturbations in the subdominant component is inhibited by dark acoustic oscillations due to its coupling to the dark radiation, solving the $\sigma_8$ problem, while the presence of tightly coupled dark radiation ameliorates the $H_0$ problem. The subdominant component of dark matter and dark radiation continue to remain in thermal equilibrium until late times, inhibiting the formation of a dark disk. We present an example of a simple model that naturally realizes this scenario in which both constituents of dark matter are thermal WIMPs. Our scenario can be tested by future stage-IV experiments designed to probe the CMB and large scale structure.
NoAM for "No Action Method" is a framework for reconstructing the past orbits of observed tracers of the large scale mass density field. It seeks exact solutions of the equations of motion (EoM), satisfying initial homogeneity and the final observed particle (tracer) positions. The solutions are found iteratively reaching a specified tolerance defined as the RMS of the distance between reconstructed and observed positions. Starting from a guess for the initial conditions, NoAM advances particles using standard N-body techniques for solving the EoM. Alternatively, the EoM can be replaced by any approximation such as Zel'dovich and second order perturbation theory (2LPT). NoAM is suitable for billions of particles and can easily handle non-regular volumes, redshift space, and other constraints. We implement NoAM to systematically compare Zel'dovich, 2LPT, and N-body dynamics over diverse configurations ranging from idealized high-res periodic simulation box to realistic galaxy mocks. Our findings are (i) Non-linear reconstructions with Zel'dovich, 2LPT, and full dynamics perform better than linear theory only for idealized catalogs in real space. For realistic catalogs, linear theory is the optimal choice for reconstructing velocity fields smoothed on scales > 5 Mpc/h. (ii) all non-linear back-in-time reconstructions tested here, produce comparable enhancement of the baryonic oscillation signal in the correlation function.
With the steadily improving sensitivity afforded by current and future galaxy surveys, a robust extraction of two-point correlation function measurements may become increasingly hampered by the presence of astrophysical foregrounds or observational systematics. The concept of mode projection has been introduced as a means to remove contaminants for which it is possible to construct a spatial map reflecting the expected signal contribution. Owing to its computational efficiency compared to minimum-variance methods, the sub-optimal pseudo-Cl (PCL) power spectrum estimator is a popular tool for the analysis of high-resolution data sets. Here, we integrate mode projection into the framework of PCL power spectrum estimation. In contrast to results obtained with optimal estimators, we show that the uncorrected projection of template maps leads to biased power spectra. Based on analytical calculations, we find exact closed-form expressions for the expectation value of the bias and demonstrate that they can be recast in a form that allows a numerically efficient evaluation, preserving the favorable O(l_max^3) time complexity of PCL estimator algorithms. Using simulated data sets, we assess the scaling of the bias with various analysis parameters and demonstrate that it can be reliably removed. We conclude that in combination with mode projection, PCL estimators allow for a fast and robust computation of power spectra in the presence of systematic effects -- properties in high demand for the analysis of ongoing and future large scale structure surveys.
The substructures of light bosonic (axion-like) dark matter may condense into compact Bose stars. We study collapses of the critical-mass stars caused by attractive self-interaction of the axion-like particles and find that these processes proceed in an unexpected universal way. First, nonlinear self-similar evolution (similar to "wave collapse" in plasma physics) forces the particles to fall into the star center. Second, collisions in the dense center create an outgoing stream of mildly relativistic particles which carries away an essential part of the star mass. The collapse stops when the star remnant is no longer able to support the self-similar infall feeding the collisions. We shortly discuss possible astrophysical and cosmological implications of these phenomena.
We present a measurement of two-dimensional (2D) redshift-space power spectrum for the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 11 CMASS galaxies in the North Galactic Cap (NGC) based on the method developed by Jing & Borner (2001). In this method, we first measure the 2D redshift-space correlation function for the CMASS galaxies, and obtain the 2D power spectrum based on Fourier Transform of the correlation function. The method is tested with an N-body mock galaxy catalog, which demonstrates that the method can yield an accurate and unbiased measurement of the redshift-space power spectrum given the input 2D correlation function is correct. Compared with previous measurements in literature that are usually based on direct Fourier Transform in redshift space, our method has the advantages that the window function and shot-noise are fully corrected. In fact, our 2D power spectrum, by its construction, can accurately reproduce the 2D correlation function, and in the meanwhile can reproduce, for example, the 2D power spectrum of Beutler et al. (2014) accurately if ours is convolved with the window function they provided. Thus, our measurement can facilitate a direct comparison with the theoretical predictions. With this accurate measurement of the 2D power spectrum, we then develop a method to measure the structure growth rate, by separating the anisotropic redshift-space power spectrum from the isotropic real-space power spectrum. We have also carefully corrected for the nonlinearities in the mapping from real space to redshift space, according to the theoretical model of Zhang et al. (2013). Finally, we obtain f(zeff)sigma_8(zeff)=0.438\pm0.037 at the effective redshift zeff=0.57, where f(zeff) is the linear growth rate at redshift zeff. The result is useful for constraining cosmological parameters. The measurements of 2D power spectrum will be released soon.
We present a new, semi-analytic framework for estimating the level of residuals present in CMB maps derived from multi-frequency Cosmic Microwave Background (CMB) data and forecasting their impact on cosmological parameters. The data are assumed to contain non-negligible signals of astrophysical and/or Galactic origin, which we clean using parametric component separation technique. We account for discrepancies between the foreground model assumed during the separation procedure and the true one, allowing for differences in scaling laws and/or their spatial variations. Our estimates and their uncertainties include both systematic and statistical effects and are averaged over the instrumental noise and CMB signal realizations. The framework can be further extended to account self-consistently for existing uncertainties in the foreground models. We demonstrate and validate the framework on simple study cases which aim at estimating the tensor-to-scalar ratio, r. The proposed approach is computationally efficient permitting an investigation of hundreds of set-ups and foreground models on a single CPU.
The brightest southern quasar above redshift $z=1$, HE 0515$-$4414, with its strong intervening metal absorption-line system at $z_{abs}=1.1508$, provides a unique opportunity to precisely measure or limit relative variations in the fine-structure constant ($\Delta\alpha/\alpha$). A variation of just $\sim$3 parts per million (ppm) would produce detectable velocity shifts between its many strong metal transitions. Using new and archival observations from the Ultraviolet and Visual Echelle Spectrograph (UVES) we obtain an extremely high signal-to-noise ratio spectrum (peaking at S/N $\approx250$ pix$^{-1}$). This provides the most precise measurement of $\Delta\alpha/\alpha$ from a single absorption system to date, $\Delta\alpha/\alpha=-1.42\pm0.55_{\rm stat}\pm0.65_{\rm sys}$ ppm, comparable with the precision from previous, large samples of $\sim$150 absorbers. The largest systematic error in all (but one) previous similar measurements, including the large samples, was long-range distortions in the wavelength calibration. These would add a $\sim$2 ppm systematic error to our measurement and up to $\sim$10 ppm to other measurements using Mg and Fe transitions. However, we corrected the UVES spectra using well-calibrated spectra of the same quasar from the High Accuracy Radial velocity Planet Searcher (HARPS), leaving a residual 0.59 ppm systematic uncertainty, the largest contribution to our total systematic error. A similar approach, using short observations on future, well-calibrated spectrographs to correct existing, high S/N spectra, would efficiently enable a large sample of reliable $\Delta\alpha/\alpha$ measurements. The high S/N UVES spectrum also provides insights into analysis difficulties, detector artifacts and systematic errors likely to arise from 25-40-m telescopes.
We have investigated a recently proposed halo-based model, Camelus, for predicting weak-lensing peak counts, and compared its results over a collection of 162 cosmologies with those from N-body simulations. While counts from both models agree for peaks with $\mathcal{S/N}>1$ (where $\mathcal{S/N}$ is the ratio of the peak height to the r.m.s. shape noise), we find $\approx 50\%$ fewer counts for peaks near $\mathcal{S/N}=0$ and significantly higher counts in the negative $\mathcal{S/N}$ tail. Adding shape noise reduces the differences to within $20\%$ for all cosmologies. We also found larger covariances that are more sensitive to cosmological parameters. As a result, credibility regions in the $\{\Omega_m, \sigma_8\}$ are $\approx 30\%$ larger. Even though the credible contours are commensurate, each model draws its predictive power from different types of peaks. Low peaks, especially those with $2<\mathcal{S/N}<3$, convey important cosmological information in N-body data, as shown in \cite{DietrichHartlap, Kratochvil2010}, but \textsc{Camelus} constrains cosmology almost exclusively from high significance peaks $(\mathcal{S/N}>3)$. Our results confirm the importance of using a cosmology-dependent covariance with at least a 14\% improvement in parameter constraints. We identified the covariance estimation as the main driver behind differences in inference, and suggest possible ways to make Camelus even more useful as a highly accurate peak count emulator.
CMB Stage-4 experiments will reduce the uncertainties on the gravitational lensing potential by an order of magnitude compared to current measurements, and will also produce a Sunyaev-Zel'dovich (SZ) cluster catalog containing $\sim10^{5}$ objects, two orders of magnitudes higher than what is currently available. In this paper we propose to combine these two observables and show that it is possible to calibrate the masses of the full Stage-4 cluster catalog internally owing to the high signal to noise measurement of the CMB lensing convergence field. We find that a CMB Stage-4 experiment will constrain the hydrostatic bias parameter to sub-percent accuracy. We also show constraints on a non parametric $Y-M$ relationship which could be used to study its evolution with mass and redshift. Finally we present a joint likelihood for thermal SZ (tSZ) flux and mass measurements, and show that it could lead to a $\sim5\sigma$ detection of the lower limit on the sum of the neutrino masses in the normal hierarchy ($\sum m_{\nu}=60 \textrm{meV}$) once combined with measurements of the primordial CMB and CMB lensing power spectra.
The detection of gravitational waves from the merger of binary black holes by the LIGO Collaboration has opened a new window to astrophysics. With the sensitivities of ground based detectors in the coming years we can only detect the local black hole binary mergers. The integrated merger rate can instead be probed by the gravitational-wave background, the incoherent superposition of the released energy in gravitational waves during binary-black-hole coalescence. Through that, the properties of the binary black holes can be studied. In this work we show that by measuring the energy density $\Omega_{GW}$ (in units of the cosmic critical density) of the gravitational-wave background, we can search for the rare $\sim 100 M_{\odot}$ massive black holes formed in the Universe. In addition, we can answer how often the least massive BHs of mass $> 3 M_{\odot}$ form. Finally, if there are multiple channels for the formation of binary black holes and if any of them predicts a narrow mass range for the black holes, then the gravitational-wave-background spectrum may have features that with the future Einstein Telescope can be detected.
When modelling inflaton fluctuations as a free quantum scalar field, the initial vacuum is conventionally imposed at the infinite past. This is called the Bunch-Davies (BD) vacuum. If however an asymptotically Minkowskian past does not exist, this becomes inconsistent. We derive corrections to the scalar spectral index $n_s$ and the tensor tilt $n_t$ descending from arbitrary mixed states or from explicit non-BD initial conditions. The former may stem from some pre-inflationary background and can redshift away whereas the latter are induced by a timelike hypersurface parametrising a physical cut-off. In both cases, we find that corrections scale in parts or fully as $\mathcal O(\epsilon)$ where $\epsilon$ is the first slow-roll parameter. The precise observational footprint is hence dependent on the model driving inflation. Further, we show how the inflationary consistency relation is altered. We thus provide an analytic handle on possible high scale or pre-inflationary physics.
The spatial distribution of the metals residing in the intra-cluster medium (ICM) of galaxy clusters records all the information on a cluster's nucleosynthesis and chemical enrichment history. We present measurements from deep Suzaku and Chandra observations of the cool-core galaxy cluster Abell 3112 out its virial radius (~1470 kpc). We find that the ratio of the observed supernova type Ia explosions to the total supernova explosions have a uniform distribution at a level of 12-16% out to the cluster's virial radius. The non-varying supernova enrichment suggests that the ICM was enriched by metals at an early stage before the cluster itself was formed. We also find that the 2D delayed detonations models CDDT produce significantly worse fits to the X-ray spectra compared to simple 1D W7 models. This may indicate that CDDT explosions are not a dominant process of enriching the ICM.
Cosmic backreaction refers to the general question of whether a homogeneous and isotropic cosmological model is able to predict the correct expansion dynamics of our inhomogeneous Universe. One aspect of this issue concerns the validity of the continuous approximation: does a system of point masses expand the same way as a fluid does? This article shows that it is not exactly the case in Newtonian gravity, although the associated corrections vanish in an infinite Universe. It turns out that Gauss's law is a key ingredient for such corrections to vanish. Backreaction therefore generically arises in alternative theories of gravitation, which threatens the trustworthiness of their cosmological tests. This phenomenon is illustrated with a toy-model of massive gravity.
We present the study of a sample of nine QSO fields, with damped-Ly-alpha (DLA) or sub-DLA systems at z~0.6, observed with the X-Shooter spectrograph at the Very Large Telescope. By suitably positioning the X-Shooter slit based on high spatial resolution images of HST/ACS we are able to detect absorbing galaxies in 7 out of 9 fields (~ 78\% success rate) at impact parameters from 10 to 30 kpc. In 5 out of 7 fields the absorbing galaxies are confirmed via detection of multiple emission lines at the redshift of DLAs where only 1 out of 5 also emits a faint continuum. In 2 out of these 5 fields we detect a second galaxy at the DLA redshift. Extinction corrected star formation rates (SFR) of these DLA-galaxies, estimated using their H-alpha fluxes, are in the range 0.3-6.7 M_sun yr^-1. The emission metallicities of these five DLA-galaxies are estimated to be from 0.2 to 0.9 Z_sun. Based on the Voigt profile fits to absorption lines we find the metallicity of the absorbing neutral gas to be in a range of 0.05--0.6 Z_sun. The two remaining DLA-galaxies are quiescent galaxies with SFR < 0.4 M_sun yr^-1 (3-sigma) presenting continuum emission but weak or no emission lines. Using X-Shooter spectrum we estimate i-band absolute magnitude of -19.5+/-0.2 for both these DLA-galaxies that indicates they are sub-L* galaxies. Comparing our results with that of other surveys in the literature we find a possible redshift evolution of the SFR of DLA-galaxies.
Observational evidence suggests that some very large supermassive black holes
(SMBHs) already existed less than 1 Gyr after the Big Bang. Explaining the
formation and growth of the 'seeds' of these SMBHs is quite challenging. We
explore the formation of such seeds in the direct collapse scenario. Using 3D
hydrodynamical simulations, we investigate the impact of turbulence and
rotation on the fragmentation behavior of collapsing primordial gas in the
presence of a strong UV radiation background, which keeps the gas hot.
Additionally, we explore different ways in which the collapsing gas may be able
to stay hot, and thus limit fragmentation. Using a one-zone model, we examine
the interplay between magnetic fields, turbulence, and a UV radiation
background.
Feedback processes from stars and black holes shape the interstellar medium
(ISM) out of which new generations of luminous objects form. To understand the
properties of these objects, e.g. the stellar initial mass function, it is
vital to have knowledge of the chemical and thermodynamical properties of the
feedback-regulated ISM. To better understand the chemo-thermal state and
fragmentation behavior of gas in high-redshift galaxies, we updated, improved,
and extended a photodissociation region code. Our computational code, PDR-Zz,
is described in detail. Using this code, a grid of models is run, covering a
sizable range in physical properties. This allows us to systematically explore
the overall impact of various feedback effects, both radiative and chemical, on
the chemical and thermal balance of the gas in different physical regimes.
We apply a novel spectral graph technique, that of locally-biased semi-supervised eigenvectors, to study the diversity of galaxies. This technique permits us to characterize empirically the natural variations in observed spectra data, and we illustrate how this approach can be used in an exploratory manner to highlight both large-scale global as well as small-scale local structure in Sloan Digital Sky Survey (SDSS) data. We use this method in a way that simultaneously takes into account the measurements of spectral lines as well as the continuum shape. Unlike Principal Component Analysis, this method does not assume that the Euclidean distance between galaxy spectra is a good global measure of similarity between all spectra, but instead it only assumes that local difference information between similar spectra is reliable. Moreover, unlike other nonlinear dimensionality methods, this method can be used to characterize very finely both small-scale local as well as large-scale global properties of realistic noisy data. The power of the method is demonstrated on the SDSS Main Galaxy Sample by illustrating that the derived embeddings of spectra carry an unprecedented amount of information. By using a straightforward global or unsupervised variant, we observe that the main features correlate strongly with star formation rate and that they clearly separate active galactic nuclei. Computed parameters of the method can be used to describe line strengths and their interdependencies. By using a locally-biased or semi-supervised variant, we are able to focus on typical variations around specific objects of astronomical interest. We present several examples illustrating that this approach can enable new discoveries in the data as well as a detailed understanding of very fine local structure that would otherwise be overwhelmed by large-scale noise and global trends in the data.
There exists a class of ultralight Dark Matter (DM) models which could form a Bose-Einstein condensate (BEC) in the early universe and behave as a single coherent wave instead of individual particles in galaxies. We show that a generic BEC DM halo intervening along the line of sight of a gravitational wave (GW) signal could induce an observable change in the speed of GW, with the effective refractive index depending only on the mass and self-interaction of the constituent DM particles and the GW frequency. Hence, we propose to use the deviation in the speed of GW as a new probe of the BEC DM parameter space. With a multi-messenger approach to GW astronomy and/or with extended sensitivity to lower GW frequencies, the entire BEC DM parameter space can be effectively probed by our new method in the near future.
Dark matter simulations can serve as a basis for creating galaxy histories
via the galaxy-dark matter connection. Here, one such model by Becker (2015) is
implemented with several variations on three different dark matter simulations.
Stellar mass and star formation rates are assigned to all simulation subhalos
at all times, using subhalo mass gain to determine stellar mass gain. The
observational properties of the resulting galaxy distributions are compared to
each other and observations for a range of redshifts from 0-2. Although many of
the galaxy distributions seem reasonable, there are noticeable differences as
simulations, subhalo mass gain definitions, or subhalo mass definitions are
altered, suggesting that the model should change as these properties are
varied. Agreement with observations may improve by including redshift
dependence in the added-by-hand random contribution to star formation rate.
There appears to be an excess of faint quiescent galaxies as well (perhaps due
in part to differing definitions of quiescence). The ensemble of galaxy
formation histories for these models tend to have more scatter around their
average histories (for a fixed final stellar mass) than the two more predictive
and elaborate semi-analytic models of Guo et al (2013) and Henriques et al
(2015), and require more basis fluctuations (using PCA) to capture 90 percent
of the scatter around their average histories.
The codes to plot model predictions (in some cases alongside observational
data) are publicly available to test other mock catalogues at
https://github.com/jdcphysics/validation/codes/vsuite . Information on how to
use these codes is in the appendix.
Links to: arXiv, form interface, find, astro-ph, recent, 1609, contact, help (Access key information)
In this paper we investigate how observational effects could possibly bias cosmological inferences from peculiar velocity measurements. Specifically, we look at how bulk flow measurements are compared with theoretical predictions. Usually bulk flow calculations try to approximate the flow that would occur in a sphere around the observer. Using the Horizon Run 2 simulation we show that the traditional methods for bulk flow estimation can overestimate the magnitude of the bulk flow for two reasons: when the survey geometry is not spherical (the data do not cover the whole sky), and when the observations undersample the velocity distributions. Our results may explain why several bulk flow measurements found bulk flow velocities that seem larger than those expected in standard {\Lambda}CDM cosmologies. We recommend a different approach when comparing bulk flows to cosmological models, in which the theoretical prediction for each bulk flow measurement is calculated specifically for the geometry and sampling rate of that survey. This means that bulk flow values will not be comparable between surveys, but instead they are comparable with cosmological models, which is the more important measure.
The recent analysis of low-redshift supernovae (SN) has increased the apparent tension between the value of $H_0$ estimated from low and high red-shift observations such as the cosmic microwave background (CMB) radiation. On the other hand other observations have provided strong evidence for the existence of a local underdensity extending up to a red-shift of about $0.07$. We compute with different methods the effects of this local void on the low-redshift luminosity distance using an exact solution of the Einstein's equations, linear perturbation theory and a low-redshift expansion. The correction is proportional to the volume averaged density contrast and to the comoving distance form the center and is able to completely resolve the apparent $H_0$ tension. The void does not affect the high red-shift luminosity distance because the volume averaged density contrast tends to zero asymptotically. Since all the Cepheids used for the luminosity distance calibration are inside this local void, not properly taking into account these effects leads to a miss-estimation of the value of the Hubble constant. The corrections to the luminosity distance are due to the monopole component of the peculiar velocity component associated with the void. Previous analysis missed this correction because the integral of the density field necessary to obtain the total peculiar velocity was performed up to red-shift $0.06$, which is not a sufficiently large scale to include the effects of the local void.
We present the first release of the "Galaxy Cluster Merger Catalog". This catalog provides an extensive suite of mock observations and related data for N-body and hydrodynamical simulations of galaxy cluster mergers. These mock observations consist of projections of a number of important observable quantities in several different wavebands, for the entire evolution of each simulation as well as along different lines of sight through the three-dimensional simulation domain. The web interface to the catalog consists of easily browseable images over epoch and projection direction, as well as download links for the raw data and a JS9 interface for interactive data exploration. All of the data products are provided in the standard FITS file format, in image and table form. Data is being stored on the yt Hub (this http URL), which allows for remote access and analysis using a Jupyter notebook server. Future versions of the catalog will include simulations from a number of research groups and a variety of research topics related to the study of interactions of galaxy clusters with each other and with their member galaxies. The catalog is located at this http URL
We demonstrate that Chromo-Natural Inflation can be made consistent with observational data if the SU(2) gauge symmetry is spontaneously broken. Working in the Stueckelberg limit, we show that isocurvature is negligible, and the resulting adiabatic fluctuations can match current observational constraints. Observable levels of chirally-polarized gravitational radiation ($r\sim 10^{-3}$) can be produced while the evolution of all background fields is sub-Planckian. The gravitational wave spectrum is amplified via linear mixing with the gauge field fluctuations, and its amplitude is not simply set by the Hubble rate during inflation. This allows observable gravitational waves to be produced for an inflationary energy scale below the GUT scale. The tilt of the resulting gravitational wave spectrum can be either blue or red.
A new frontier in the search for dark matter (DM) is based on the idea of detecting the decoherence caused by DM scattering against a mesoscopic superposition of normal matter. Such superpositions are uniquely sensitive to very small momentum transfers from new particles and forces, especially DM with a mass below 100 MeV. Here we investigate what sorts of dark sectors are inaccessible with existing methods but would induce noticeable decoherence in the next generation of matter interferometers. We show that very soft, but medium range (0.1 nm - 1 {\mu}m) elastic interactions between matter and DM are particularly suitable. We construct toy models for such interactions, discuss existing constraints, and delineate the expected sensitivity of forthcoming experiments. The first hints of DM in these devices would appear as small variations in the anomalous decoherence rate with a period of one sidereal day. This is a generic signature of interstellar sources of decoherence, clearly distinguishing it from terrestrial backgrounds. The OTIMA experiment under development in Vienna will begin to probe Earth-thermalizing DM once sidereal variations in the background decoherence rate are pushed below one part in a hundred for superposed 5-nm gold nanoparticles. The proposals by Bateman et al. and Geraci et al. could be similarly sensitive, although they would require at least a month of data taking. DM that is absorbed or elastically reflected by the Earth, and so avoids a greenhouse density enhancement, would not be detectable by those three experiment. On the other hand, aggressive proposals of the MAQRO collaboration and Pino et al. would immediately open up many orders of magnitude in DM mass, interaction range, and coupling strength, regardless of how DM behaves in bulk matter.
We consider thick domain walls in a de Sitter universe following paper by Basu and Vilenkin. However, we are interested not only in stationary solutions found therein, but also investigate the general case of domain wall evolution with time. When the wall thickness parameter, $\delta_0$, is smaller than $H^{-1}/\sqrt{2}$, where $H$ is the Hubble parameter in de Sitter space-time, then the stationary solutions exist, and initial field configurations tend with time to the stationary ones. However, there are no stationary solutions for $\delta_0 \geq H^{-1}/\sqrt{2}$. We have calculated numerically the rate of the wall expansion in this case and have found that the width of the wall grows exponentially fast for $\delta_0 \gg H^{-1}$. An explanation for the critical value $\delta_{0c} = H^{-1}/\sqrt{2}$ is also proposed.
We report the discovery of a new ultra-faint dwarf satellite companion of the Milky Way based on the early survey data from the Hyper Suprime-Cam Subaru Strategic Program. This new satellite, Virgo I, which is located in the constellation of Virgo, has been identified as a statistically significant (5.5 sigma) spatial overdensity of star-like objects with a well-defined main sequence and red giant branch in their color-magnitude diagram. The significance of this overdensity increases to 10.8 sigma when the relevant isochrone filter is adopted for the search. Based on the distribution of the stars around the likely main sequence turn-off at r ~ 24 mag, the distance to Virgo I is estimated as 87 kpc, and its most likely absolute magnitude calculated from a Monte Carlo analysis is M_V = -0.8 +/- 0.9 mag. This stellar system has an extended spatial distribution with a half-light radius of 38 +12/-11 pc, which clearly distinguishes it from a globular cluster with comparable luminosity. Thus, Virgo I is one of the faintest dwarf satellites known and is located beyond the reach of the Sloan Digital Sky Survey. This demonstrates the power of this survey program to identify very faint dwarf satellites. This discovery of VirgoI is based only on about 100 square degrees of data, thus a large number of faint dwarf satellites are likely to exist in the outer halo of the Milky Way.
QUBIC is an instrument aiming at measuring the B mode polarisation anisotropies at medium scales angular scales (30-200 multipoles). The search for the primordial CMB B-mode polarization signal is challenging, because of many difficulties: smallness of the expected signal, instrumental systematics that could possibly induce polarization leakage from the large E signal into B, brighter than anticipated polarized foregrounds (dust) reducing to zero the initial hope of finding sky regions clean enough to have a direct primordial B-modes observation. The QUBIC instrument is designed to address all aspects of this challenge with a novel kind of instrument, a Bolometric Interferometer, combining the background-limited sensitivity of Transition-Edge-Sensors and the control of systematics allowed by the observation of interference fringe patterns, while operating at two frequencies to disentangle polarized foregrounds from primordial B mode polarization. Its characteristics are described in details in this Technological Design Report.
Links to: arXiv, form interface, find, astro-ph, recent, 1609, contact, help (Access key information)
Semi-numeric methods have made it possible to efficiently model the epoch of reionisation (EoR). While most implementations involve a reduction to a simple three-parameter model, we introduce a new mass-dependent ionising efficiency parameter that folds in physical parameters that are constrained by the latest numerical simulations. This new parameterization enables the effective modeling of a broad range of host halo masses containing ionising sources, extending from the smallest Population III host halos with $M \sim 10^6 M_\odot$, which are often ignored, to the rarest cosmic peaks with $M \sim 10^{12} M_\odot$ during EoR. We compare the resulting ionising histories with a typical three-parameter model and also compare with the latest constraints from the Planck mission. Our model results in a optical depth due to Thomson scattering, $\tau_{\mathrm{e}}$ = 0.057, that is consistent with Planck. The largest difference in our model is shown in the resulting bubble size distributions which peak at lower characteristic sizes and are broadened. We also consider the uncertainties of the various physical parameters and comparing the resulting ionising histories broadly disfavors a small contribution from galaxies. As the smallest haloes cease a meaningful contribution to the ionising photon budget after $z = 10$, implying they play a role in determining the start of EoR and little else, an ionising photon escape fraction greater than 0.6 is also broadly disfavored.
Intensity mapping is a promising technique for surveying the large scale structure of our Universe from $z=0$ to $z \sim 150$, using the brightness temperature field of spectral lines to directly observe previously unexplored portions of out cosmic timeline. Examples of targeted lines include the $21\,\textrm{cm}$ hyperfine transition of neutral hydrogen, rotational lines of carbon monoxide, and fine structure lines of singly ionized carbon. Recent efforts have focused on detections of the power spectrum of spatial fluctuations, but have been hindered by systematics such as foreground contamination. This has motivated the decomposition of data into Fourier modes perpendicular and parallel to the line-of-sight, which has been shown to be a particularly powerful way to diagnose systematics. However, such a method is well-defined only in the limit of a narrow-field, flat-sky approximation. This limits the sensitivity of intensity mapping experiments, as it means that wide surveys must be separately analyzed as a patchwork of smaller fields. In this paper, we develop a framework for analyzing intensity mapping data in a spherical Fourier-Bessel basis, which incorporates curved sky effects without difficulty. We use our framework to generalize a number of techniques in intensity mapping data analysis from the flat sky to the curved sky. These include visibility-based estimators for the power spectrum, treatments of interloper lines, and the "foreground wedge" signature of spectrally smooth foregrounds.
Conventional Type Ia supernova (SN Ia) cosmology analyses currently use a simplistic linear regression of magnitude versus color and light curve shape, which does not model intrinsic SN Ia variations and host galaxy dust as physically distinct effects, resulting in low color-magnitude slopes. We construct a probabilistic generative model for the distribution of dusty extinguished absolute magnitudes and apparent colors as a convolution of the intrinsic SN Ia color-magnitude distribution and the host galaxy dust reddening-extinction distribution. If the intrinsic color-magnitude (M_B vs. B-V) slope beta_int differs from the host galaxy dust law R_B, this convolution results in a specific curve of mean extinguished absolute magnitude vs. apparent color. The derivative of this curve smoothly transitions from beta_int in the blue tail to R_B in the red tail of the apparent color distribution. The conventional linear fit approximates this effective curve at this transition near the average apparent color, resulting in an apparent slope beta_app between beta_int and R_B. We incorporate these effects into a hierarchical Bayesian statistical model for SN Ia light curve measurements, and analyze a dataset of SALT2 optical light curve fits of a compilation of 277 nearby SN Ia at z < 0.10. The conventional linear fit obtains beta_app = 3. Our model finds a beta_int = 2.2 +/- 0.3 and a distinct dust law of R_B = 3.7 +/- 0.3, consistent with the average for Milky Way dust, while correcting a systematic distance bias of ~0.10 mag in the tails of the apparent color distribution.
By monitoring a large number of stars in the Local Group galaxies such as M33 with an 8\,m-class telescope with time integration of $\sim 100\,$sec per shot, we can detect microlensing events by sub-lunar mass compact objects (SULCOs) such as primordial black holes (PBHs) and rogue (free-floating) dwarf planets. For one night observation, we would be able to detect $10^{3-4}$ microlensing events caused by SULCOs with a mass of $10^{-9}$ to $10^{-7}$ solar mass for sources with S/N$>5$ if SULCOs constitute all the dark matter components. Moreover, we expect $10^{1-2}$ events in which sources with S/N$>100$ are weakly amplified due to lensing by SULCOs with a mass range of $10^{-11}$ to $10^{-7}$ solar mass. The method would provide a stringent constraint on the abundance of SULCOs at the distance $0.1-100$ kpc from us.
The determination of the resolution of cosmological N-body simulations, i.e., the range of scales in which quantities measured in them represent accurately the continuum limit, is an important open question. We address it here using scale-free models, for which self-similarity provides a powerful tool to control resolution. Such models also provide a robust testing ground for the so-called stable clustering approximation, which gives simple predictions for them. Studying large N-body simulations of such models with different force smoothing, we find that these two issues are in fact very closely related: our conclusion is that resolution in the non-linear regime extends, in practice, down to the scale at which stable clustering breaks down. Physically the association of the two scales is in fact simple to understand: stable clustering fails to be a good approximation when there are strong interactions of structures (in particular merging) and it is precisely such non-linear processes which are sensitive to fluctuations at the smaller scales affected by discretisation. Resolution may be further degraded if the short distance gravitational smoothing scale is larger than the scale to which stable clustering can propagate. We examine in detail the very different conclusions of a study by Smith et al. (2003), and find that the strong deviations from stable clustering reported by this work are the results of over-optimistic assumptions about resolution of the measured power spectra, and the reliance on Fourier space analysis. We emphasize the much poorer resolution obtained with the power spectrum compared to the two point correlation function.
We present an analytical solution for the luminosity distance in spatially flat cosmology with pressureless matter and the cosmological constant. The complex analytical solution is made of a real part and a negligible imaginary part. The real part of the luminosity distance allows finding the two parameters $H_0$ and $\om$. A simple expression for the distance modulus for SNs of type Ia is reported in the framework of the minimax approximation.
The classical equations of motion for an axion with potential $V(\phi)=m_a^2f_a^2 [1-\cos (\phi/f_a)]$ possess quasi-stable, localized, oscillating solutions, which we refer to as "axion stars". We study, for the first time, collapse of axion stars numerically using the full non-linear Einstein equations of general relativity and the full non-perturbative cosine potential. We map regions on an "axion star stability diagram", parameterized by the initial ADM mass, $M_{\rm ADM}$, and axion decay constant, $f_a$. We identify three regions of the parameter space: i) long-lived oscillating axion star solutions, with a base frequency, $m_a$, modulated by self-interactions, ii) collapse to a BH and iii) complete dispersal due to gravitational cooling and interactions. We locate the boundaries of these three regions and an approximate "triple point" $(M_{\rm TP},f_{\rm TP})\sim (2.4 M_{pl}^2/m_a,0.3 M_{pl})$. For $f_a$ below the triple point BH formation proceeds during winding (in the complex $U(1)$ picture) of the axion field near the dispersal phase. This could prevent astrophysical BH formation from axion stars with $f_a\ll M_{pl}$. For larger $f_a\gtrsim f_{\rm TP}$, BH formation occurs through the stable branch and we estimate the mass ratio of the BH to the stable state at the phase boundary to be $\mathcal{O}(1)$ within numerical uncertainty. We discuss the observational relevance of our findings for axion stars as BH seeds, which are supermassive in the case of ultralight axions. For the QCD axion, the typical BH mass formed from axion star collapse is $M_{\rm BH}\sim 3.4 (f_a/0.6 M_{pl})^{1.2} M_\odot$.
In the hope of avoiding model dependence of the cosmological observables, phenomenological parametrizations of Cosmic Inflation have recently been proposed. Typically, they are expressed in terms of two parameters associated with an expansion of the inflationary quantities matching the belief that inflation is characterized by two numbers only, the tensor-to-scalar ratio and the scalar spectral index. We give different arguments and examples showing that these new approaches are either not generic or insufficient to make predictions at the accuracy level needed by the cosmological data. We conclude that disconnecting inflation from high energy physics and gravity might not be the most promising way to learn about the physics of the early Universe.
We study the structure of two-point correlators of the inflationary field fluctuations in order to improve the accuracy and efficiency of the existing spectral methods. We present a description motivated by the separation of the fast and slow evolving components of the spectrum. Our purpose is to rephrase all the relevant equations of motion in terms of slowly varying quantities. This is important in order to consider the contribution from high-frequency modes to the spectrum without affecting computational performance. The slow-roll approximation is not required to reproduce the main distinctive features in the power spectrum for each specific model of inflation.
We conduct an analysis of the Planck 2015 data that is complete in reionization observables from the large angle polarization $E$-mode spectrum in the redshift range $6 < z < 30$. Based on 5 principal components, all of which are constrained by the data, this single analysis can be used to infer constraints on any model for reionization in the same range; we develop an effective likelihood approach for applying these constraints to models. By allowing for an arbitrary ionization history, this technique tests the robustness of inferences on the total optical depth from the usual step-like transition assumption, which is important for the interpretation of many other cosmological parameters such as the dark energy and neutrino mass. The Planck 2015 data not only allow a high redshift $z>15$ component to the optical depth but prefer it at the $2\sigma$ level. This preference is associated with excess power in the multipole range $10 \lesssim \ell \lesssim 20$ and may indicate high redshift ionization sources or unaccounted for systematics and foregrounds in the 2015 data.
We present spatially-resolved two-dimensional stellar kinematics for the 41 most massive early-type galaxies (MK <~ -25.7 mag, stellar mass M* >~ 10^11.8 Msun) of the volume-limited (D < 108 Mpc) MASSIVE survey. For each galaxy, we obtain high-quality spectra in the wavelength range of 3650 to 5850 A from the 246-fiber Mitchell integral-field spectrograph (IFS) at McDonald Observatory, covering a 107" x 107" field of view (often reaching 2 to 3 effective radii). We measure the 2-D spatial distribution of each galaxy's angular momentum (lambda and fast or slow rotator status), velocity dispersion (sigma), and higher-order non-Gaussian velocity features (Gauss-Hermite moments h3 to h6). Our sample contains a high fraction (~80% ) of slow and non-rotators with lambda <~ 0.2. When combined with the lower-mass ETGs in the ATLAS3D survey, we find the fraction of slow-rotators to increase dramatically with galaxy mass, reaching ~50% at MK ~ -25.5 mag and ~90% at MK <~ -26 mag. All of our fast rotators show a clear anti-correlation between h3 and V/sigma, and the slope of the anti-correlation is steeper in more round galaxies. The radial profiles of sigma show a clear luminosity and environmental dependence: the 12 most luminous galaxies in our sample (MK <~ -26 mag) are all brightest cluster/group galaxies (except NGC 4874) and all have rising or nearly flat sigma profiles, whereas five of the seven "isolated" galaxies are all fainter than MK = -25.8 mag and have falling sigma. All of our galaxies have positive average h4; the most luminous galaxies have average h4 ~ 0.05 while less luminous galaxies have a range of values between 0 and 0.05. Most of our galaxies show positive radial gradients in h4, and those galaxies also tend to have rising sigma profiles. We discuss the implications for the relationship among dynamical mass, sigma, h4, and velocity anisotropy for these massive galaxies.
The formation of supermassive stars (SMSs) via rapid mass accretion and their direct collapse into black holes (BHs) is a promising pathway for sowing seeds of supermassive BHs in the early universe. We calculate the evolution of rapidly accreting SMSs by solving the stellar structure equations including nuclear burning as well as general relativistic (GR) effects up to the onset of the collapse. We find that such SMSs have less concentrated structure than fully-convective counterpart, which is often postulated for non-accreting ones. This effect stabilizes the stars against GR instability even above the classical upper mass limit $\gtrsim 10^5~M_\odot$ derived for the fully-convective stars. The accreting SMS begins to collapse at the higher mass with the higher accretion rate. The collapse occurs when the nuclear fuel is exhausted only for cases with $\dot M \lesssim 0.1~M_\odot~{\rm yr}^{-1}$. With $\dot{M} \simeq 0.3 - 1~M_\odot~{\rm yr}^{-1}$, the star becomes GR-unstable during the helium-burning stage at $M \simeq 2 - 3.5~\times 10^5~M_\odot$. In an extreme case with $10~M_\odot~{\rm yr}^{-1}$, the star does not collapse until the mass reaches $\simeq 8.0\times 10^5~M_\odot$, where it is still in the hydrogen-burning stage. We expect that BHs with roughly the same mass will be left behind after the collapse in all the cases.
We present photometry and time-series spectroscopy of the nearby type Ia supernova (SN Ia) SN 2015F over $-16$ days to $+80$ days relative to maximum light, obtained as part of the Public ESO Spectroscopic Survey of Transient Objects (PESSTO). SN 2015F is a slightly sub-luminous SN Ia with a decline rate of $\Delta m15(B)=1.35 \pm 0.03$ mag, placing it in the region between normal and SN 1991bg-like events. Our densely-sampled photometric data place tight constraints on the epoch of first light and form of the early-time light curve. The spectra exhibit photospheric C II $\lambda 6580$ absorption until $-4$ days, and high-velocity Ca II is particularly strong at $<-10$ days at expansion velocities of $\simeq$23000\kms. At early times, our spectral modelling with syn++ shows strong evidence for iron-peak elements (Fe II, Cr II, Ti II, and V II) expanding at velocities $>14000$ km s$^{-1}$, suggesting mixing in the outermost layers of the SN ejecta. Although unusual in SN Ia spectra, including V II in the modelling significantly improves the spectral fits. Intriguingly, we detect an absorption feature at $\sim$6800 \AA\ that persists until maximum light. Our favoured explanation for this line is photospheric Al II, which has never been claimed before in SNe Ia, although detached high-velocity C II material could also be responsible. In both cases the absorbing material seems to be confined to a relatively narrow region in velocity space. The nucleosynthesis of detectable amounts of Al II would argue against a low-metallicity white dwarf progenitor. We also show that this 6800 \AA\ feature is weakly present in other normal SN Ia events, and common in the SN 1991bg-like sub-class.
We investigate whether small perturbations can cause relaxation to quantum equilibrium over very long timescales. We consider in particular a two-dimensional harmonic oscillator, which can serve as a model of a field mode on expanding space. We assume an initial wave function with small perturbations to the ground state. We present evidence that the trajectories are highly confined so as to preclude relaxation to equilibrium even over very long timescales. Cosmological implications are briefly discussed.
We describe general features that might be observed in the line spectra of relic cosmological particles should quantum nonequilibrium be preserved in their statistics. According to our arguments, these features would represent a significant departure from those of a conventional origin. Among other features, we find a possible spectral broadening (for incident photons) that is proportional to the energy resolution of the recording telescope (and so could be orders of magnitude larger than any intrinsic broadening). Notably, for a range of possible initial conditions we find the possibility of spectral line `narrowing' whereby a telescope could observe a spectral line which is narrower than it should conventionally be able to resolve. We briefly discuss implications for the indirect search for dark matter.
We investigate a new method to search for keV-scale sterile neutrinos that could account for Dark Matter. Neutrinos trapped in our galaxy could be captured on stable $^{163}$Dy if their mass is greater than 2.83~keV. Two experimental realizations are studied, an integral counting of $^{163}$Ho atoms in dysprosium-rich ores and a real-time measurement of the emerging electron spectrum in a dysprosium-based detector. The capture rates are compared to the solar neutrino and radioactive backgrounds. An integral counting experiment using several kilograms of $^{163}$Dy could reach a sensitivity for the sterile-to-active mixing angle $\sin^2\theta_{e4}$ of $10^{-5}$ significantly exceeding current laboratory limits. Mixing angles as low as $\sin^2\theta_{e4} \sim 10^{-7}$ / $\rm m_{^{163}\rm Dy}\rm{(ton)}$ could possibly be explored with a real-time experiment.
The current status of the phenomenology of short-baseline neutrino oscillations induced by light sterile neutrinos in the framework of 3+1 mixing is reviewed.
We revisit the nonthermal gravitino production at the (p)reheating stage after inflation. Particular attention is paid to large field inflation models with a $\mathbb{Z}_2$ symmetry, for which the previous perturbative analysis is inapplicable; and inflation models with a stabilizer superfield, which have not been studied non-perturbatively. It is found that in single-superfield inflation models (without the stabilizer field), nonthermal production of the transverse gravitino can be cosmologically problematic while the abundance of the longitudinal gravitino is small enough. In multi-superfield inflation models (with the stabilizer field), production of the transverse gravitino is significantly suppressed, and it is cosmologically harmless. We also clarify the relation between the background field method used in the preheating context and the standard perturbative decay method to estimate the gravitino abundance.
Hot, Dust-Obscured Galaxies (Hot DOGs), selected from the WISE all sky infrared survey, host some of the most powerful Active Galactic Nuclei (AGN) known, and might represent an important stage in the evolution of galaxies. Most known Hot DOGs are at $z> 1.5$, due in part to a strong bias against identifying them at lower redshift related to the selection criteria. We present a new selection method that identifies 153 Hot DOG candidates at $z\sim 1$, where they are significantly brighter and easier to study. We validate this approach by measuring a redshift $z=1.009$, and an SED similar to higher redshift Hot DOGs for one of these objects, WISE J1036+0449 ($L_{\rm\,Bol}\simeq 8\times 10^{46}\rm\,erg\,s^{-1}$), using data from Keck/LRIS and NIRSPEC, SDSS, and CSO. We find evidence of a broadened component in MgII, which, if due to the gravitational potential of the supermassive black hole, would imply a black hole mass of $M_{\rm\,BH}\simeq 2 \times 10^8 M_{\odot}$, and an Eddington ratio of $\lambda_{\rm\,Edd}\simeq 2.7$. WISE J1036+0449 is the first Hot DOG detected by NuSTAR, and the observations show that the source is heavily obscured, with a column density of $N_{\rm\,H}\simeq(2-15)\times10^{23}\rm\,cm^{-2}$. The source has an intrinsic 2-10 keV luminosity of $\sim 6\times 10^{44}\rm\,erg\,s^{-1}$, a value significantly lower than that expected from the mid-infrared/X-ray correlation. We also find that the other Hot DOGs observed by X-ray facilities show a similar deficiency of X-ray flux. We discuss the origin of the X-ray weakness and the absorption properties of Hot DOGs. Hot DOGs at $z\lesssim1$ could be excellent laboratories to probe the characteristics of the accretion flow and of the X-ray emitting plasma at extreme values of the Eddington ratio.
Links to: arXiv, form interface, find, astro-ph, recent, 1609, contact, help (Access key information)