Search

We investigate if the discrepancy between estimates of the total baryon mass fraction obtained from observations of the cosmic microwave background (CMB) and of galaxy groups/clusters persists when a large sample of groups is considered. To this purpose, 91 candidate X-ray groups/poor clusters at redshift 0.1 < z < 1 are selected from the COSMOS 2 deg^2 survey, based only on their X-ray luminosity and extent. This sample is complemented by 27 nearby clusters with a robust, analogous determination of the total and stellar mass inside R_500. The total sample of 118 groups and clusters with z < 1 spans a range in M_500 of ~10^13--10^15 M_sun. We find that the stellar mass fraction associated with galaxies at R_500 decreases with increasing total mass as (M_500)^-0.37 \pm 0.04, independent of redshift. Estimating the total gas mass fraction from a recently derived, high quality scaling relation, the total baryon mass fraction (f_500^stars+gas=f_500^stars+f_500^gas) is found to increase by ~ 25% when M_500 increases from =5 X 10^13 M_sun to = 7 X 10^14 M_sun. After consideration of a plausible contribution due to intra--cluster light (11--22% of the total stellar mass), and gas depletion through the hierarchical assembly process (10% of the gas mass), the estimated values of the total baryon mass fraction are still lower than the latest CMB measure of the same quantity (WMAP5), at a significance level of 3.3\sigma for groups of =5 X 10^13~M_sun. The discrepancy decreases towards higher total masses, such that it is 1\sigma at = 7 X 10^14~M_sun. We discuss this result in terms of non--gravitational processes such as feedback and filamentary heating., Comment: 14 pages, accepted for publication in ApJ. Revised discussion, added section on systematics

We present an update on the NZ-wide advances in the field of Radio Astronomy and Radio Engineering with a particular focus on contributions, not thus reported elsewhere, which hope to either directly or indirectly contribute to New Zealand's engagement with the international Square Kilometre Array (SKA) project. We discuss the status of the SKA project in New Zealand with particular reference to activities of the New Zealand Square Kilometre Array Research and Development Consortium., Comment: 6 pages, 2 figures, published in Proceedings of the 16th Electronics New Zealand Conference (ENZcon), held 18-20 Nov, 2009, Dunedin, New Zealand

The Devasthal Fast Optical Telescope (DFOT) is a 1.3 meter aperture optical telescope, recently installed at Devasthal, Nainital. We present here the first results using an \Ha filter with this telescope on a Wolf-Rayet dwarf galaxy Mrk 996. The instrumental response and the \Ha sensitivity obtained with the telescope are $(3.3 \pm 0.3)\times 10^{-15} (erg s^{-1} cm^{-2})/(counts s^{-1})$ and $7.5\times10^{-17} erg s^{-1} cm^{-2} arcsec^{-2}$ respectively. The \Ha flux and the equivalent width for Mrk~996 are estimated as $(132 \pm 37)\times 10^{-14} erg s^{-1} cm^{-2}$ and $\sim$96 $\mathring{A}$ respectively. The star formation rate is estimated as $0.4\pm0.1 M_\odot$yr$^{-1}$. Mrk~996 deviates from the radio-FIR correlation known for normal star forming galaxies with a deficiency in its radio continuum. The ionized gas as traced by \Ha emission is found in a disk shape which is misaligned with respect to the old stellar disk. This misalignment is indicative of a recent tidal interaction in the galaxy. We believe that galaxy-galaxy tidal interaction is the main cause of the WR phase in Mrk~996., Comment: 13 pages, 2 figures Accepted for publication in "Journal of Astrophysics and Astronomy"

Although dark matter is a central element of modern cosmology, the history of how it became accepted as part of the dominant paradigm is often ignored or condensed into a brief anecdotical account focused around the work of a few pioneering scientists. The aim of this review is to provide the reader with a broader historical perspective on the observational discoveries and the theoretical arguments that led the scientific community to adopt dark matter as an essential part of the standard cosmological model., 86 pages, 8 figures

The Massive and Distant Clusters of WISE Survey (MaDCoWS) provides a catalog of high-redshift ($0.7\lesssim z\lesssim 1.5$) infrared-selected galaxy clusters. However, the verification of the ionized intracluster medium, indicative of a collapsed and nearly virialized system, is made challenging by the high redshifts of the sample members. The main goal of this work is to test the capabilities of the Atacama Compact Array (ACA; also known as the Morita Array) Band 3 observations, centered at about 97.5 GHz, to provide robust validation of cluster detections via the thermal Sunyaev-Zeldovich (SZ) effect. Using a pilot sample that comprises ten MaDCoWS galaxy clusters, accessible to ACA and representative of the median sample richness, we infer the masses of the selected galaxy clusters and respective detection significance by means of a Bayesian analysis of the interferometric data. Our test of the "Verification with the ACA - Localization and Cluster Analysis" (VACA LoCA) program demonstrates that the ACA can robustly confirm the presence of the virialized intracluster medium in galaxy clusters previously identified in full-sky surveys. In particular, we obtain a significant detection of the SZ effect for seven out of the ten VACA LoCA clusters. We note that this result is independent of the assumed pressure profile. However, the limited angular dynamic range of the ACA in Band 3 alone, short observational integration times, and possible contamination from unresolved sources limit the detailed characterization of the cluster properties and the inference of the cluster masses within scales appropriate for the robust calibration of mass-richness scaling relations., 19 pages (including appendices), 14 figures, and 4 tables; accepted for publication in A&A

We study the spatial fluctuations in the hydrogen ionizing background in the epoch following reionization (z ~ 5--6). The rapid decrease with redshift in the photon mean free path (m.f.p.), combined with the clustering of increasingly rare ionizing sources, can result in a very inhomogenous ionizing background during this epoch. We systematically investigate the probability density functions (PDFs) and power spectra of ionizing flux, by varying several parameters such as the m.f.p., minimum halo mass capable of hosting stars, and halo duty cycle. In order to be versatile, we make use of analytic, semi-numeric and numeric approaches. Our models show that the ionizing background indeed has sizable fluctuations during this epoch sourced by the clustering of sources, with the PDFs being a factor of few wide at half of the maximum likelihood. The distributions also show marked asymmetries, with a high-value tail set by clustering on small scales, and a shorter low-value tail which is set by the mean free path. The power spectrum of the ionizing background is much more sensitive to source properties than the PDF and can be well-understood analytically with a framework similar to the halo model (usually used to describe dark matter clustering). Nevertheless, we find that Lya forest spectra are extremely insensitive to the details of the UVB, despite marked differences in the PDFs and power spectra of our various ionizing backgrounds. Assuming a uniform ionizing background only underestimates the value of the mean ionization rate inferred from the Lya forest by a few percent. Instead, analysis of the Lya forest is dominated by the uncertainties in the density field. Thus, our results justify the common assumption of a uniform ionizing background in Lya forest analysis even during this epoch., 11 pages, 11 figures, submitted to the MNRAS

We study how the cosmological constraints from growth data are improved by including the measurements of bias from Dark Energy Survey (DES). In particular, we utilize the biasing properties of the DES Luminous Red Galaxies (LRGs) and the growth data provided by the various galaxy surveys in order to constrain the growth index ($\gamma$) of the linear matter perturbations. Considering a constant growth index we can put tight constraints, up to $\sim 10\%$ accuracy, on $\gamma$. Specifically, using the priors of the Dark Energy Survey and implementing a joint likelihood procedure between theoretical expectations and data we find that the best fit value is in between $\gamma=0.64\pm 0.075$ and $0.65\pm 0.063$. On the other hand utilizing the Planck priors we obtain $\gamma=0.680\pm 0.089$ and $0.690\pm 0.071$. This shows a small but non-zero deviation from General Relativity ($\gamma_{\rm GR}\approx 6/11$), nevertheless the confidence level is in the range $\sim 1.3-2\sigma$. Moreover, we find that the estimated mass of the dark-matter halo in which LRGs survive lies in the interval $\sim 6.2 \times 10^{12} h^{-1} M_{\odot}$ and $1.2 \times 10^{13} h^{-1} M_{\odot}$, for the different bias models. Finally, allowing $\gamma$ to evolve with redshift [Taylor expansion: $\gamma(z)=\gamma_{0}+\gamma_{1}z/(1+z)$] we find that the $(\gamma_{0},\gamma_{1})$ parameter solution space accommodates the GR prediction at $\sim 1.7-2.9\sigma$ levels., Comment: 8 pages, 3 figures, discussion added, to appear in European Physical Journal C (EPJC)

The primary advantage of moderately superheated bubble chamber detectors is their simultaneous sensitivity to nuclear recoils from WIMP dark matter and insensitivity to electron recoil backgrounds. A comprehensive analysis of PICO gamma calibration data demonstrates for the first time that electron recoils in C$_3$F$_8$ scale in accordance with a new nucleation mechanism, rather than one driven by a hot-spike as previously supposed. Using this semi-empirical model, bubble chamber nucleation thresholds may be tuned to be sensitive to lower energy nuclear recoils while maintaining excellent electron recoil rejection. The PICO-40L detector will exploit this model to achieve thermodynamic thresholds as low as 2.8 keV while being dominated by single-scatter events from coherent elastic neutrino-nucleus scattering of solar neutrinos. In one year of operation, PICO-40L can improve existing leading limits from PICO on spin-dependent WIMP-proton coupling by nearly an order of magnitude for WIMP masses greater than 3 GeV c$^{-2}$ and will have the ability to surpass all existing non-xenon bounds on spin-independent WIMP-nucleon coupling for WIMP masses from 3 to 40 GeV c$^{-2}$., 18 pages, 12 figures

Cosmic inflation provides a mechanism for generating the early density perturbations that seeded the large-scale structures we see today. Primordial non-Gaussianity is among the most promising of few observational tests of physics at this epoch. At present non-Gaussianity is best constrained by the cosmic microwave background, but in the near term large-scale structure data may be competitive so long as the effects of primordial non-Gaussianity can be modeled through the non-linear process of structure formation. We discuss recent work modeling effects of a few types of primordial non-Gaussianity on the large-scale halo clustering and the halo mass function. More specifically, we compare analytic and N-body results for two variants of the curvaton model of inflation: (i) a "tauNL" scenario in which the curvaton and inflaton contribute equally to the primordial curvature perturbation and (ii) a "gNL" model where the usual quadratic fNL term in the potential cancels, but a large cubic term remains., Comment: 9 pages, 4 figures, Proceedings of the DPF-2011 Conference, Providence, RI, August 8-13, 2011

We present submillimetre observations of the z=3.653 quasar SDSS160705+533558 together with data in the optical and infrared. The object is unusually bright in the far-IR and submm with an IR luminosity of ~10^14 L_sun. We ascribe this luminosity to a combination of AGN and starburst emission, with the starburst forming stars at a rate of a few thousand solar masses per year. Submillimetre Array (SMA) imaging observations with a resolution ~1" show that the submm (850 micron) emission is extended on scales of 10--35kpc and is offset from the optical position by ~10 kpc. This morphology is dissimilar to that found in submm galaxies, which are generally un- or marginally resolved on arcsecond scales, or submm-luminous AGN where the AGN lies at the peak of the submm or molecular emission. The simplest explanation is that the object is in the early stages of a merger between a gas rich galaxy, which hosts the starburst, and a gas-poor AGN-host galaxy, which is responsible for the quasar emission. It is also possible that jet induced star formation might contribute to the unusual morphology., ApJ Letters, in press

We consider the problem of building inhomogeneous cosmological models in scalar-tensor theories of gravity. This starts by splitting the field equations of these theories into constraint and evolution equations, and then proceeds by identifying exact solutions to the constraints. We find exact, closed form expressions for geometries that correspond to the initial data for cosmological models containing regular arrays of point-like masses. These solutions extend similar methods that have recently been applied to Einstein's equations, and provides sufficient initial conditions to perform numerical integration of the evolution equations. We use our new solutions to study the effects of inhomogeneity in cosmologies governed by scalar-tensor theories of gravity, including the spatial inhomogeneity allowed in Newton's constant. Finally, we compare our solutions to their general relativistic counterparts, and investigate the effect of changing the coupling constant between the scalar and tensor degrees of freedom., Comment: 18 pages, 8 figures

We present high-resolution 870-um ALMA continuum maps of 30 bright sub-millimeter sources in the UKIDSS UDS field. These sources are selected from deep, 1-square degrees 850-um maps from the SCUBA--2 Cosmology Legacy Survey, and are representative of the brightest sources in the field (median SCUBA2 flux S_850=8.7+/-0.4 mJy). We detect 52 sub-millimeter galaxies (SMGs) at >4-sigma significance in our 30 ALMA maps. In 61+/-17% of the ALMA maps the single-dish source comprises a blend of >=2 SMGs, where the secondary SMGs are Ultra--Luminous Infrared Galaxies (ULIRGs) with L_IR>10^12 Lo. The brightest SMG contributes on average 80+/-4% of the single-dish flux density, and in the ALMA maps containing >=2 SMGs the secondary SMG contributes 25+/-3% of the integrated ALMA flux. We construct source counts and show that multiplicity boosts the apparent single-dish cumulative counts by 20% at S_870>7.5mJy, and by 60% at S_870>12mJy. We combine our sample with previous ALMA studies of fainter SMGs and show that the counts are well-described by a double power-law with a break at 8.5+/-0.6mJy. The break corresponds to a luminosity of ~6x10^12Lsol or a star-formation rate of ~1000Mo/yr. For the typical sizes of these SMGs, which are resolved in our ALMA data with r=1.2+/-0.1kpc, this yields a limiting SFR density of ~100Msol/yr/kpc2. Finally, the number density of S_870>2mJy SMGs is 80+/-30 times higher than that derived from blank-field counts. An over-abundance of faint SMGs is inconsistent with line-of-sight projections dominating multiplicity in the brightest SMGs, and indicates that a significant proportion of these high-redshift ULIRGs must be physically associated., 14 pages, 6 figures. ApJ in press

We revisit the physical properties of global and local monopoles and discuss their implications in the dynamics of monopole networks. In particular, we review the Velocity-dependent One-Scale (VOS) model for global and local monopoles and propose physically motivated changes to its equations. We suggest a new form for the acceleration term of the evolution equation of the root-mean-squared velocity and show that, with this change, the VOS model is able to describe the results of radiation and matter era numerical simulations of global monopole networks with a single value of the acceleration parameter $k$, thus resolving the tension previously found in the literature. We also show that the fact that the energy of global monopoles is not localized within their cores affects their dynamics and, thus, the Hubble damping terms in the VOS equations. We study the ultra-relativistic linear scaling regime predicted by the VOS equations and demonstrate that it cannot be attained either on radiation or matter eras and, thus, cannot arise from the cosmological evolution of a global monopole network. We also briefly discuss the implications of our findings for the VOS model for local monopoles., Comment: 8 pages, 2 figures

We investigate the gravitational instability of galactic discs, treating stars and cold interstellar gas as two distinct components, and taking into account the phenomenology of turbulence in the interstellar medium (ISM), i.e. the Larson-type scaling relations observed in the molecular and atomic gas. Besides deriving general properties of such systems, we analyse a large sample of galaxies from The HI Nearby Galaxy Survey (THINGS), and show in detail how interstellar turbulence affects disc instability in star-forming spirals. We find that turbulence has a significant effect on both the inner and the outer regions of the disc. In particular, it drives the inner gas disc to a regime of transition between two instability phases and makes the outer disc more prone to star-dominated instabilities., MNRAS, in press. Moderate revision to match the accepted version

Gravitational lensing by large-scale structure significantly impacts observations of the cosmic microwave background (CMB): it smooths the acoustic peaks in temperature and $E$-mode polarization power spectra, correlating previously uncorrelated modes; and it converts $E$-mode polarization into $B$-mode polarization. The act of measuring and removing the effect of lensing from CMB maps, or delensing, has been well studied in the context of $B$ modes, but little attention has been given to the delensing of the temperature and $E$ modes. In this paper, we model the expected delensed $T$ and $E$ power spectra to all orders in the lensing potential, demonstrating the sharpening of the acoustic peaks and a significant reduction in lens-induced power spectrum covariances. We then perform cosmological forecasts, demonstrating that delensing will yield improved sensitivity to parameters with upcoming surveys. We highlight the breaking of the degeneracy between the effective number of neutrino species and primordial helium fraction as a concrete application. We also show that delensing increases cosmological information as long as the measured lensing reconstruction is included in the analysis. We conclude that with future data, delensing will be crucial not only for primordial $B$-mode science but for a range of other observables as well., Comment: To be submitted to JCAP. Comments welcome

The canonical acoustic dark energy model (cADE), which is based on a scalar field with a canonical kinetic term that rapidly converts potential to kinetic energy around matter radiation equality, alleviates the Hubble tension found in $\Lambda$CDM. We show that it successfully passes new consistency tests in the CMB damping tail provided by the ACT data, while being increasingly constrained and distinguished from alternate mechanisms by the improved CMB acoustic polarization data from Planck. The best fit cADE model to a suite of cosmological observations, including the SH0ES $H_0$ measurement, has $H_0=70.25$ compared with $68.23$ (km s$^{-1}$ Mpc$^{-1}$) in $\Lambda$CDM and a finite cADE component is preferred at the $2.8\sigma$ level. The ability to raise $H_0$ is now mainly constrained by the improved Planck acoustic polarization data, which also plays a crucial role in distinguishing cADE from the wider class of early dark energy models. ACT and Planck TE polarization data are currently mildly discrepant in normalization and drive correspondingly different preferences in parameters. Improved constraints on intermediate scale polarization approaching the cosmic variance limit will be an incisive test of the acoustic dynamics of these models and their alternatives., Comment: v2: Added more references; matched the published version. 10 pages, 8 figures

Detection of 21~cm emission of HI from the epoch of reionization, at redshifts z>6, is limited primarily by foreground emission. We investigate the signatures of wide-field measurements and an all-sky foreground model using the delay spectrum technique that maps the measurements to foreground object locations through signal delays between antenna pairs. We demonstrate interferometric measurements are inherently sensitive to all scales, including the largest angular scales, owing to the nature of wide-field measurements. These wide-field effects are generic to all observations but antenna shapes impact their amplitudes substantially. A dish-shaped antenna yields the most desirable features from a foreground contamination viewpoint, relative to a dipole or a phased array. Comparing data from recent Murchison Widefield Array observations, we demonstrate that the foreground signatures that have the largest impact on the HI signal arise from power received far away from the primary field of view. We identify diffuse emission near the horizon as a significant contributing factor, even on wide antenna spacings that usually represent structures on small scales. For signals entering through the primary field of view, compact emission dominates the foreground contamination. These two mechanisms imprint a characteristic "pitchfork" signature on the "foreground wedge" in Fourier delay space. Based on these results, we propose that selective down-weighting of data based on antenna spacing and time can mitigate foreground contamination substantially by a factor ~100 with negligible loss of sensitivity., Published in ApJ

We show that the term maintaining conservation of the charged vector current for the transitions "neutron proton" even for different masses of the neutron and proton (see T. Leitner et al., Phys. Rev. C {\bf 73}, 065502 (2006) and A. M. Ankowski, arXiv:1601.06169 [hep-ph]) is related to the first class current contribution but not to the second class one as has been pointed out by C. Giunti, arXiv: 1602.00215 [hep-ph]., Comment: 3 pages, no figures

We propose an improved methodology to constrain spatial variations of the fine structure constant using clusters of galaxies. We use the {\it Planck} 2013 data to measure the thermal Sunyaev-Zeldovich effect at the location of 618 X-ray selected clusters. We then use a Monte Carlo Markov Chain algorithm to obtain the temperature of the Cosmic Microwave Background at the location of each galaxy cluster. When fitting three different phenomenological parameterizations allowing for monopole and dipole amplitudes in the value of the fine structure constant we improve the results of earlier analysis involving clusters and the CMB power spectrum, and we also found that the best-fit direction of a hypothetical dipole is compatible with the direction of other known anomalies. Although the constraining power of our current datasets do not allow us to test the indications of a fine-structure constant dipole obtained though high-resolution optical/UV spectroscopy, our results do highlight that clusters of galaxies will be a very powerful tool to probe fundamental physics at low redshift., 11 pages, 5 figures and 3 tables. Accepted for publication in Physical Review D

A significant fraction of the total photospheric light in nearby galaxy clusters is thought to be contained within the diffuse intracluster light (ICL), which extends 100s of kpc from cluster cores. The study of the ICL can reveal details of the evolutionary histories and processes occurring within galaxy clusters, however since it has a very low surface brightness it is often difficult to detect. We present here the first measurements of the ICL as a fraction of total cluster light at z \sim 1 using deep J-band (1.2 ��m) imaging from HAWK-I on the VLT. We investigate the ICL in 6 X-ray selected galaxy clusters at 0.8< z, 12 pages, 12 figures. Accepted for publication in MNRAS

We point out a possible generation mechanism of non-Gaussian bubbles in the sky due to bubble nucleation in the early universe. We consider a curvaton scenario for inflation and assume that the curvaton field phi, whose energy density is subdominant during inflation but which is responsible for the curvature perturbation of the universe, is coupled to another field sigma which undergoes false vacuum decay through quantum tunneling. For this model, we compute the skewness of the curvaton fluctuations due to its interaction with sigma during tunneling, that is, on the background of an instanton solution that describes false vacuum decay. We find that the resulting skewness of the curvaton can become large in the spacetime region inside the bubble. We then compute the corresponding skewness in the statistical distribution of the cosmic microwave background (CMB) temperature fluctuations. We find a non-vanishing skewness in a bubble-shaped region in the sky. It can be large enough to be detected in the near future, and if detected it will bring us invaluable information about the physics in the early universe., 6 pages, 6 figures

We explore the Minkowski functionals of weak lensing convergence map to distinguish between $f(R)$ gravity and the general relativity (GR). The mock weak lensing convergence maps are constructed with a set of high-resolution simulations assuming different gravity models. It is shown that the lensing MFs of $f(R)$ gravity can be considerably different from that of GR because of the environmentally dependent enhancement of structure formation. We also investigate the effect of lensing noise on our results, and find that it is likely to distinguish F5, F6 and GR gravity models with a galaxy survey of $\sim3000$ degree$^2$ and with a background source number density of $n_g=30~{\rm arcmin}^{-2}$, comparable to an upcoming survey dark energy survey (DES). We also find that the $f(R)$ signal can be partially degenerate with the effect of changing cosmology, but combined use of other observations, such as the cosmic microwave background (CMB) data, can help break this degeneracy., 7 pages, 3 figures, Accepted by Physical Review D

Light sub-GeV dark matter (DM) constitutes an underexplored target, beyond the optimized sensitivity of typical direct DM detection experiments. We comprehensively investigate hadrophilic light DM produced from cosmic-ray collisions with the atmosphere. The resulting relativistic DM, originating from meson decays, can be efficiently observed in variety of experiments, such as XENON1T. We include for the first time decays of $\eta$, $\eta^{\prime}$ and $K^+$ mesons, leading to improved limits for DM masses above few hundred MeV. We incorporate an exact treatment of the DM attenuation in Earth and demonstrate that nuclear form factor effects can significantly impact the resulting testable DM parameter space. Further, we establish projections for upcoming experiments, such as DARWIN, over a wide range of DM masses below the GeV scale., Comment: 7 pages, 4 figures; matches published version

In the Standard Model, neutrinos are massless, yet oscillation experiments show in fact they do have a small mass. Currently only the differences of the masses' squares are known, and an upper bound on the sum. However, upcoming surveys of the Universe's large-scale structure (LSS) can probe the neutrino mass by exposing how neutrinos modulate galaxy clustering. But these measurements are challenging: in looking at the clustering of galaxy pairs, the effect of neutrinos is degenerate with galaxy formation, the details of which are unknown. Marginalizing over them degrades the constraints. Here we show that using correlations of galaxy triplets---the 3-Point Correlation Function or its Fourier-space analog the bispectrum---can break the degeneracy between galaxy formation physics (known as biasing) and the neutrino mass. Specifically, we find a signature of neutrinos in the bispectrum's dipole moment (with respect to triangle opening angle) that is roughly orthogonal to the contribution of galaxy biases. This signature was missed in previous works by failing to account for how neutrinos alter mode-coupling between perturbations on different scales. Our proposed signature will contribute to upcoming LSS surveys' such as DESI making a robust detection of the neutrino mass. We estimate that it can offer several-$\sigma$ evidence for non-zero $m_{\nu}$ with DESI from the bispectrum alone, and that this is independent from information in the galaxy power spectrum., Comment: 9 pages including references, submitted

There is a growing interest of using the intrinsic alignment (IA) of galaxy images as a tool to extract cosmological information complimentary to galaxy clustering analysis. Recently, Okumura & Taruya derived useful formulas for the intrinsic ellipticity--ellipticity correlation, the gravitational shear--intrinsic ellipticity correlation, and the velocity--intrinsic ellipticity correlation functions based on the linear alignment (LA) model. In this paper, using large-volume $N$-body simulations, we measure these alignment statistics in real and redshift space and compare them to the LA and nonlinear alignment model predictions. We find that anisotropic features of baryon acoustic oscillations in the IA statistics can be accurately predicted by our models. The anisotropy due to redshift-space distortions (RSDs) is also well described in the large-scale limit. Our results indicate that one can extract the cosmological information encoded in the IA through the Alcock-Paczynski and RSD effects., 9 pages, 6 figures

A new class of inflation models within the context of G-inflation is proposed, in which the standard model Higgs boson can act as an inflaton thanks to Galileon-like non-linear derivative interaction. The generated primordial density perturbation is shown to be consistent with the present observational data. We also make a general discussion on potential-driven G-inflation models, and find a new consistency relation between the tensor-to-scalar ratio $r$ and the tensor spectral index $n_T$, $r = -32 \sqrt{6}n_T / 9$, which is crucial in discriminating the present models from standard inflation with a canonical kinetic term., 10 pages, 4 figures

Recent studies suggest that the quenching properties of galaxies are correlated over several mega-parsecs. The large-scale "galactic conformity" phenomenon around central galaxies has been regarded as a potential signature of "galaxy assembly bias" or "pre-heating", both of which interpret conformity as a result of direct environmental effects acting on galaxy formation. Building on the iHOD halo quenching framework developed in Zu & Mandelbaum (2015, 2016), we discover that our fiducial halo mass quenching model, without any galaxy assembly bias, can successfully explain the overall environmental dependence and the conformity of galaxy colours in SDSS, as measured by the mark correlation functions of galaxy colours and the red galaxy fractions around isolated primaries, respectively. Our fiducial iHOD halo quenching mock also correctly predicts the differences in the spatial clustering and galaxy-galaxy lensing signals between the more vs. less red galaxy subsamples, split by the red-sequence ridge-line at fixed stellar mass. Meanwhile, models that tie galaxy colours fully or partially to halo assembly bias have difficulties in matching all these observables simultaneously. Therefore, we demonstrate that the observed environmental dependence of galaxy colours can be naturally explained by the combination of 1) halo quenching and 2) the variation of halo mass function with environment --- an indirect environmental effect mediated by two separate physical processes., 19 pages, 10 figures, comments welcome

In single-field models of inflation the effect of a long mode with momentum q reduces to a diffeomorphism at zeroth and first order in q. This gives the well-known consistency relations for the n-point functions. At order q^2 the long mode has a physical effect on the short ones, since it induces curvature, and we expect that this effect is the same as being in a curved FRW universe. In this paper we verify this intuition in various examples of the three-point function, whose behaviour at order q^2 can be written in terms of the power spectrum in a curved universe. This gives a simple alternative understanding of the level of non-Gaussianity in single-field models. Non-Gaussianity is always parametrically enhanced when modes freeze at a physical scale k_{ph, f} shorter than H: f_{NL} \sim (k_{ph, f}/H)^2., Comment: 18 pages, 1 figure. v2: small changes, JCAP published version

Using the dynamical system approach, we describe the general dynamics of cosmological scalar fields in terms of critical points and heteroclinic lines. It is found that critical points describe the initial and final states of the scalar field dynamics, but that heteroclinic lines which give a more complete description of the evolution in between the critical points. In particular, the heteroclinic line that departs from the (saddle) critical point of perfect fluid-domination is the representative path in phase space of quintessence fields that may be viable dark energy candidates. We also discuss the attractor properties of the heteroclinic lines, and their importance for the description of thawing and freezing fields., Comment: Minor changes to the text and two new figures, main conclusions unchanged. 12 pages, 11 figures, uses RevTex

Galaxies accrete their mass by means of both smooth accretion from the cosmic web, and the mergers of smaller entities. We wish to quantify the respective role of these two modes of accretion, which could determine the morphological types of galaxies observed today. Multi-zoom cosmological simulations are used to estimate as a function of time the evolution of mass in bound systems, for dark matter as well as baryons. The baryonic contents of dark matter haloes are studied. Merger histories are followed as a function of external density, and the different ways in which mass is assembled in galaxies and the stellar component accumulated are quantified. We find that most galaxies assemble their mass through smooth accretion, and only the most massive galaxies also grow significantly through mergers. The mean fraction of mass assembled by accretion is 77 %, and by mergers 23 %. We present typical accretion histories of hundreds of galaxies: masses of the most massive galaxies increase monotonically in time, mainly through accretion, many intermediate-mass objects also experience mass-loss events such as tidal stripping and evaporation. However, our simulations suffer from the overcooling of massive galaxies caused by the neglect of active galaxy nuclei (AGN) feedback. The time by which half of the galay mass has assembled, both in dark matter and baryons, is a decreasing function of mass, which is compatible with the observations of a so-called downsizing. At every epoch in the universe, there are low-mass galaxies actively forming stars, while more massive galaxies form their stars over a shorter period of time within half the age of the universe., A&A Accepted, 19 pages, 17 figures

The cosmic infrared background (CIB) includes roughly half of the energy radiated by all galaxies at all wavelengths across cosmic time, as observed at the present epoch. The PACS Evolutionary Probe (PEP) survey is exploited here to study the CIB and its redshift differential, at 70, 100 and 160 micron, where the background peaks. Combining PACS observations of the GOODS-S, GOODS-N, Lockman Hole and COSMOS areas, we define number counts spanning over more than two orders of magnitude in flux: from ~1 mJy to few hundreds mJy. Stacking of 24 micron sources and P(D) statistics extend the analysis down to ~0.2 mJy. Taking advantage of the wealth of ancillary data in PEP fields, differential number counts and CIB are studied up to z=5. Based on these counts, we discuss the effects of confusion on PACS blank field observations and provide confusion limits for the three bands considered. The total CIB surface brightness emitted above PEP 3 sigma flux limits is 4.52 +/- 1.18, 8.35 +/- 0.95 and 9.49 +/- 0.59 [nW/m2/sr] at 70, 100, and 160 micron, respectively. These values correspond to 58 +/- 7% and 74 +/- 5% of the COBE/DIRBE CIB direct measurements at 100 and 160 micron. Employing the P(D) analysis, these fractions increase to ~65% and ~89%. More than half of the resolved CIB was emitted at redshift z1.0 sources reaches 50% at 160 micron. Most of the CIB resolved in the three PACS bands was emitted by galaxies with infrared luminosities in the range 1e11-1e12 L(sun)., Proposed for acceptance on A&A

We reassess the evidence that WMAP temperature maps contain a statistically significant "cold spot" by repeating the analysis using simple circular top-hat (disk) weights, as well as Gaussian weights of varying width. Contrary to previous results that used Spherical Mexican Hat Wavelets, we find no significant signal at any scale when we compare the coldest spot from our sky to ones from simulated Gaussian random, isotropic maps. We trace this apparent discrepancy to the fact that WMAP cold spot's temperature profile just happens to favor the particular profile given by the wavelet. Since randomly generated maps typically do not exhibit this coincidence, we conclude that the original cold spot significance originated at least partly due to a fortuitous choice of using a particular basis of weight functions. We also examine significance of a more general measure that returns the most significant result among several choices of the weighting function, angular scale of the spot, and the statistics applied, and again find a null result., 7 pages, 3 figures, minor additions

Many cosmological studies predict that early supermassive black holes (SMBHs) can only form in the most massive dark matter halos embedded within large scale structures marked by galaxy over-densities that may extend up to 10 physical Mpc. This scenario, however, has not been confirmed observationally, as the search for galaxy over-densities around high-z quasars has returned conflicting results. The field around the z=6.28 quasar SDSSJ1030+0524 (J1030) is unique for multi-band coverage and represents an excellent data legacy for studying the environment around a primordial SMBH. In this paper we present wide-area (25x25 arcmin), Y- and J-band imaging of the J1030 field obtained with the near infrared camera WIRCam at the Canada-France-Hawaii Telescope (CFHT). We built source catalogues in the Y- and J-band, and matched those with our photometric catalogue in the r, z, i bands presented in Morselli et al. (2014). We used these new infrared data together with H and K and Spitzer/IRAC data to refine our selection of Lyman Break Galaxies (LBGs), extending our selection criteria to galaxies in the range 25.2=1.3. We found a significant asymmetry in the distribution of the high-z galaxies in J1030, supporting the existence of a coherent large-scale structure around the quasar. We compared our results with those of Bowler et al. (2015), who adopted similar LBGs selection criteria, and estimated an over-density of galaxies in the field of delta = 2.4, which is significant at >4 sigma. The over-density value and its significance are higher than those found in Morselli et al. (2014), and we interpret this as evidence of an improved LBG selection., Comment: 13 pages, 8 figures, accepted for publication in A&A

We study the impact of ultraviolet background (UVB) radiation field and the local stellar radiation on the H_I column density distribution f(N_HI) of damped Ly-alpha systems (DLAs) and sub-DLAs at z=3 using cosmological smoothed particle hydrodynamics simulations. We find that, in the previous simulations with an optically thin approximation, the UVB was sinking into the H_I cloud too deeply, and therefore we underestimated the f(N_HI) at 19 < log(N_HI) < 21.2 compared to the observations. If the UVB is shut off in the high-density regions with n_gas > 6 x 10^{-3} cm^{-3}, then we reproduce the observed f(N_HI) at z=3 very well. We also investigate the effect of local stellar radiation by post-processing our simulation with a radiative transfer code, and find that the local stellar radiation does not change the f(N_HI) very much. Our results show that the shape of f(N_HI) is determined primarily by the UVB with a much weaker effect by the local stellar radiation and that the optically thin approximation often used in cosmological simulation is inadequate to properly treat the ionization structure of neutral gas in and out of DLAs. Our result also indicates that the DLA gas is closely related to the transition region from optically-thick neutral gas to optically-thin ionized gas within dark matter halos., 5 pages, 3 figures, uses emulateapj. ApJL in press. Corrected a plotting error in Fig. 3, which reduced the effect of local stellar radiation on f(N_HI)

The FourStar galaxy evolution survey (ZFOURGE) is a 45 night legacy program with the FourStar near-infrared camera on Magellan and one of the most sensitive surveys to date. ZFOURGE covers a total of $400\ \mathrm{arcmin}^2$ in cosmic fields CDFS, COSMOS and UDS, overlapping CANDELS. We present photometric catalogs comprising $>70,000$ galaxies, selected from ultradeep $K_s$-band detection images ($25.5-26.5$ AB mag, $5\sigma$, total), and $>80\%$ complete to $K_s\times15$., Comment: 39 pages, 35 figures, accepted for publication in ApJ, full data release available at http://zfourge.tamu.edu/

We identify the largest known lensed images of a single spiral galaxy, lying close to the centre of the distant cluster MACS J1149.5+2223 ($z=0.544$). These images cover a total area of $\simeq 150 \Box\arcsec$ and are magnified $\simeq 200$ times. Unusually, there is very little image distortion implying the central mass distribution is almost uniform over a wide area ($r\simeq200 kpc$) with a surface density equal to the critical density for lensing, corresponding to maximal lens magnification. Many fainter multiply-lensed galaxies are also uncovered by our model, outlining a very large tangential critical curve, of radius $r\simeq 170 kpc$, posing a potential challenge for the standard LCDM-Cosmology. Because of the uniform central mass distribution a particularly clean measurement of the mass of the brightest cluster galaxy is possible here, for which we infer stars contribute most of the mass within a limiting radius of $\simeq 30 kpc$, with a mass-to-light ratio of $M/L_{B}\simeq 4.5(M/L)_{\odot}$. This cluster with its uniform and central mass distribution acts analogously to a regular magnifying glass, converging light without distorting the images, resulting in the most powerful lens yet discovered for accessing the faint high-$z$ Universe., 5 pages, 5 figures, accepted to ApJ letters - minor changes made. High resolution figures available at ftp://wise-ftp.tau.ac.il/pub/adiz/macs1149.5+2223

The Weakly Interacting Massive Particles(WIMPs) have long been the favored CDM candidate in the standard $\Lambda$CDM model. However, owing to great improvement in the experimental sensitivity in the past decade, some parameter space of the SUSY-based WIMP model is ruled out. In addition, WIMP as the CDM particle is also at variance with other astrophysical observables at small scales. We consider a model that addresses both these issues. In the model, the WIMP decays into a massive particle and radiation. We study the background evolution and the first order perturbation theory (coupled Einstein-Boltzmann equations) for this model and show that the dynamics can be captured by a single parameter $r=m_L/q$, which is the ratio of the lighter mass and the comoving momentum of the decay particle. We incorporate the relevant equations in the existing Boltzmann code CLASS to compute the matter power spectra and CMB angular power spectra. The decaying WIMP model is akin to a non-thermal Warm Dark Matter(WDM) model and suppresses matter power at small scales, which could alleviate several issues that plague the CDM model. We compare the predictions of the model with CMB, galaxy clustering, and high-z HI data. Both these data sets yield $r\gtrsim 10^6$, which can be translated into the bounds on other parameters. In particular, we obtain the following lower bounds on the self-annihilation cross-section of WIMPs $\sigma$, and the lighter mass $m_L$: $\sigma\gtrsim 5\times 10^{-44}\,\rm cm^2$ and $m_L\gtrsim 2.4\,\rm keV$. The lower limit on $m_L$ is comparable to constraints on the mass of thermally-produced WDM particle. The limit on the self-annihilation cross-section greatly expands the available parameter space as compared to the stable WIMP scenario., Comment: 24 pages (main text 17 pages), 7 figures, 2 tables

We determine the scalar potential for a D3-brane in a warped conifold background subject to general ultraviolet perturbations. Incorporating the effects of imaginary anti-self-dual (IASD) fluxes and four-dimensional curvature at the nonlinear level, we compute the leading terms in the D3-brane potential. We then provide strong cross-checks of our results by reproducing them in the dual gauge theory. Finally, we observe that the D3-brane potential induced by nonperturbative effects on D7-branes can be represented by a ten-dimensional supergravity solution containing suitable IASD fluxes. Our method allows for the systematic inclusion of compactification effects and serves to constrain the D3-brane effective action in a large class of stabilized compactifications., 4 pages; v2: minor corrections

The resonantly produced sterile neutrino ($N_1$) of the neutrino minimal standard model ($\nu$MSM) is a compelling dark matter candidate, especially in the reported possible detection of $N_1$ with mass $m_\mathrm{s}=7.1$~keV in X-ray decay. This particle will be accessible to the XRISM X-ray mission over the next 12 months. We revisit the physics behind $N_1$ and the uncertainty in its parameters. We compare predictions for the $m_\mathrm{s}=7.1$keV $N_1$ mixing angle, $\sin^{2}(2\theta)$, and half-mode mass, $M_\mathrm{hm}$, to existing X-ray observations and structure formation constraints. The strongest available constraints rule out $N_1$ as a dark matter candidate, and a more optimistic reading of the data prefers $\sin^{2}(2\theta)=5\times10^{-11}$ and $M_\mathrm{hm}=3.5\times10^{8}$$\mathrm{M}_{\odot}$. We highlight that the most promising upcoming opportunity for a detection is to find a line of velocity dispersion $\sim500$$\mathrm{kms^{-1}}$ in the Virgo cluster with XRISM, and then draw up a list of future objects of study to determine: (i) whether the line is from dark matter generally, and (ii) if from dark matter, whether that candidate is indeed $N_1$., Comment: 13 pages, 7 figures. To be submitted to MNRAS. Contact: lovell@hi.is

In this work I study the problem of E/B-mode separation with binned cosmic shear two-point correlation function data. Motivated by previous work on E/B-mode separation with shear two-point correlation functions and the practical considerations of data analysis, I consider E/B-mode estimators which are linear combinations of the binned shear correlation function data points. I demonstrate that these estimators mix E- and B-modes generally. I then show how to define estimators which minimize this E/B-mode mixing and give practical recipes for their construction and use. Using these optimal estimators, I demonstrate that the vector space composed of the binned shear correlation function data points can be decomposed into approximately ambiguous, E- and B-mode subspaces. With simple Fisher information estimates, I show that a non-trivial amount of information on typical cosmological parameters is contained in the ambiguous mode subspace computed in this formalism. Next, I give two examples which apply these practical estimators and recipes to generic problems in cosmic shear data analysis: data compression and spatially locating B-mode contamination. In particular, by using wavelet-like estimators with the shear correlation functions directly, one can pinpoint B-mode contamination to specific angular scales and extract information on its shape. Finally, I discuss how these estimators can be used as part of blinded or closed-box cosmic shear data analyses in order to assess and find B-mode contamination at high-precision while avoiding observer biases., Comment: 15 pages, 5 figures, 3 appendices, MNRAS submitted, comments welcome!

Many motivated extensions of the Standard Model predict the existence of cosmic strings. Gravitational waves originating from the dynamics of the resulting cosmic string network have the ability to probe many otherwise inaccessible properties of the early universe. In this study we show how the spectrum of gravitational waves from a cosmic string network can be used to test the equation of state of the early universe prior to Big Bang Nucleosynthesis (BBN). We also demonstrate that current and planned gravitational wave detectors such as LIGO, LISA, DECIGO/BBO, and ET/CE have the potential to detect signals of a non-standard pre-BBN equation of state and evolution of the early universe (e.g., early non-standard matter domination or kination domination) or new degrees of freedom active in the early universe beyond the sensitivity of terrestrial collider experiments and cosmic microwave background measurements., Journal version, 31 pages, 8 figures

QSO emission-line spectra are compared to predictions based on theoretical ionizing continua of accretion disks. Observed line intensities do not show the expected trend of higher ionization with higher accretion disk temperature as derived from the black hole mass and accretion rate. This suggests that, at least for accretion rates close to the Eddington limit, the inner disk does not reach temperatures as high as expected from standard disk theory. Modified radial temperature profiles, taking account of winds or advection in the inner disk, achieve better agreement with observation. This conclusion agrees with an earlier study of QSO continuum colors as a function of disk temperature. The emission lines of radio-detected and radio-undetected sources show different trends as a function of disk temperature., Comment: 9 pages, 8 figures, submitted to ApJ

Galaxy clusters are the largest virialized structures in the observable Universe. The knowledge of their properties provides many useful astrophysical and cosmological information. Our aim is to derive the luminosity and stellar mass profiles of the nearby galaxy clusters of the Omega-WINGS survey and to study the main scaling relations valid for such systems. We have merged the data of the WINGS and Omega-WINGS databases, sorted the sources according to the distance from the brightest cluster galaxy (BCG) and calculated the integrated luminosity profiles in the $B$ and $V$ bands, taking into account extinction, photometric and spatial completeness, K-correction and background contribution. Then, by exploiting the spectroscopic sample we derived the stellar mass profiles of the clusters. We got the luminosity profiles of 46 galaxy clusters, reaching $r_{200}$ in 30 cases, and the stellar mass profiles of 42 of our objects. We successfully fitted all the integrated luminosity growth profiles with one or two embedded S\'ersic components, deriving the main clusters parameters. Finally, we checked the main scaling relation among the clusters parameters in comparison with those obtained for a selected sample of early-type galaxies (ETGs) of the same clusters. We found that the nearby galaxy clusters are non-homologous structures like ETGs and exhibit a color-magnitude (CM) red-sequence relation very similar to that observed for galaxies in clusters. These properties are not expected in the current cluster formation scenarios. In particular the existence of a CM relation for clusters, shown here for the first time, suggests that the baryonic structures grow and evolve in a similar way at all scales., Comment: 45 pages, 41 figures

By 1917, V.M. Slipher had singlehandedly established a tendency for 'spiral nebulae' to be redshifted (21 out of 25 cases). From a modern perspective, it could seem surprising that the expansion of the universe was not announced at this point. Examination of Slipher's papers shows that he reached a more subtle conclusion: the identification of cosmological peculiar velocities, including the bulk motion of the Milky Way, leading to a beautiful argument in favour of nebulae as distant stellar systems. Nevertheless, Slipher's data actually contain evidence at >8sigma for a positive mean velocity, even after subtracting the dipole owing to the motion of the observer. In 1929, Hubble estimated distances for a sample of no greater depth, using redshifts due almost entirely to Slipher. Hubble's distances were flawed in two distinct ways: in addition to an incorrect absolute calibration, the largest distances were systematically under-estimated. Nevertheless, he claimed the detection of a linear distance-redshift relation. Statistically, the evidence for such a correlation is less strong than the simple evidence for a positive mean velocity in Hubble's sample. Comparison with modern data shows that a sample of more than twice Hubble's depth would generally be required in order to reveal clearly the global linear expansion in the face of the 'noise' from peculiar velocities. When the theoretical context of the time is examined, the role of the de Sitter model and its prediction of a linear distance-redshift relation looms large. A number of searches for this relation were performed prior to Hubble over the period 1924-1928, with a similar degree of success. All were based on the velocities measured by Slipher, whose work from a Century ago stands out both for the precision of his measurements and for the subtle clarity of the arguments he employed to draw correct conclusions from them., Refereed contribution to the proceedings of "Origins of the Expanding Universe: 1912-1932", M. J. Way & D. Hunter, eds., ASP Conf. Ser., Vol. 471 in press. A meeting held in September 2012 to mark the Centenary of Slipher's first measurement of the radial velocity of M31

We provide and describe a library of diffuse stellar radiation fields in spiral galaxies derived using calculations of the transfer of stellar radiation from the main morphological components - disks, thin disks, and bulges - through the dusty interstellar medium. These radiation fields are self-consistent with the solutions for the integrated panchromatic spectral energy distributions (SEDs) previously presented using the same model. Because of this, observables calculated from the radiation fields, such as gamma-ray or radio emission, can be self-consistently combined with the solutions for the ultraviolet/optical/submillimeter SEDs, thus expanding the range of applicability of the radiation transfer model to a broader range of wavelengths and physical quantities. We also give analytic solutions for radiation fields in optically thin stellar disks, in stellar disks with one dust disk and in stellar disks with two dust disks. The analytic solutions for the direct light are exact and can be used as benchmarks. The analytic solutions with scattering are only approximate, becoming exact only in the extreme optically thick limit. We find strongly contrasting solutions for the spatial distribution of the radiation fields for disks, thin disks and bulges. For bulges we find a strong dependence of the radiation fields on Sersic index., accepted for publication in MNRAS; revised version contains full analytic solutions

We argue that calculating vacuum energy requires quantum field theory whose axioms are adapted to curved spacetime. In this context, we suggest that non-zero vacuum energy is connected to dynamical breaking of electroweak symmetry. The observed meV scale can be understood in terms of electroweak physics via a naive estimate. The scenario requires all particle masses to have a dynamical origin. Any Higgs particle has to be a composite, and the origin of vacuum energy might be probed at the LHC., 6 pages, written for the Gravity Research Foundation essay competition

LUX-ZEPLIN (LZ) is a dark matter detector expected to obtain world-leading sensitivity to weakly interacting massive particles (WIMPs) interacting via nuclear recoils with a ~7-tonne xenon target mass. This manuscript presents sensitivity projections to several low-energy signals of the complementary electron recoil signal type: 1) an effective neutrino magnetic moment and 2) an effective neutrino millicharge, both for pp-chain solar neutrinos, 3) an axion flux generated by the Sun, 4) axion-like particles forming the galactic dark matter, 5) hidden photons, 6) mirror dark matter, and 7) leptophilic dark matter. World-leading sensitivities are expected in each case, a result of the large 5.6t 1000d exposure and low expected rate of electron recoil backgrounds in the $, Comment: v2 updates exclusion sensitivities from single-sided to two-sided

The recent discovery of diffuse dwarf galaxies that are deficient in dark matter appears to challenge the current paradigm of structure formation in our Universe. We describe the numerical experiments to determine if the so-called dark matter deficient galaxies (DMDGs) could be produced when two gas-rich, dwarf-sized galaxies collide with a high relative velocity of $\sim 300\,{\rm kms^{-1}}$. Using idealized high-resolution simulations with both mesh-based and particle-based gravito-hydrodynamics codes, we find that DMDGs can form as high-velocity galaxy collisions separate dark matter from the warm disk gas which subsequently is compressed by shock and tidal interaction to form stars. Then using a large simulated universe IllustrisTNG, we discover a number of high-velocity galaxy collision events in which DMDGs are expected to form. However, we did not find evidence that these types of collisions actually produced DMDGs in the TNG100-1 run. We argue that the resolution of the numerical experiment is critical to realize the "collision-induced" DMDG formation scenario. Our results demonstrate one of many routes in which galaxies could form with unconventional dark matter fractions., Comment: 18 pages including 13 figures. Accepted for publication in ApJ. For associated movies and information, see https://ejshin0.github.io/page/category/visualizations.html and https://physics.snu.ac.kr/cosmo/index/research.html#dmdg

We estimate the Milky Way (MW) halo properties using satellite kinematic data including the latest measurements from Gaia DR2. With a simulation-based 6D phase-space distribution function (DF) of satellite kinematics, we can infer halo properties efficiently and without bias, and handle the selection function and measurement errors rigorously in the Bayesian framework. Applying our DF from the EAGLE simulation to 28 satellites, we obtain an MW halo mass of $M=1.23_{-0.18}^{+0.21}\times 10^{12} M_\odot$ and a concentration of $c=9.4_{ -2.1}^{ +2.8}$ with the prior based on the $M$-$c$ relation. The inferred mass profile is consistent with previous measurements but with better precision and reliability due to the improved methodology and data. Potential improvement is illustrated by combining satellite data and stellar rotation curves. Using our EAGLE DF and best-fit MW potential, we provide much more precise estimates of kinematics for those satellites with uncertain measurements. Compared to the EAGLE DF, which matches the observed satellite kinematics very well, the DF from the semi-analytical model based on the dark-matter-only simulation Millennium II (SAM-MII) over-represents satellites with small radii and velocities. We attribute this difference to less disruption of satellites with small pericenter distances in the SAM-MII simulation. By varying the disruption rate of such satellites in this simulation, we estimate a $\sim 5\%$ scatter in the inferred MW halo mass among hydrodynamics-based simulations., Comment: Accepted for publication in ApJ; data tables available at https://github.com/syrte/mw_sats_kin

A new cosmological scenario is proposed in which a light scalaron of $f (R)$ gravity plays the role of dark matter. In this scenario, the scalaron initially resides at the minimum of its effective potential while the electroweak symmetry is unbroken. At the beginning of the electroweak crossover, the evolving expectation value of the Higgs field triggers the evolution of the scalaron due to interaction between these fields. After the electroweak crossover, the oscillating scalaron can represent cold dark matter. Its current energy density depends on a single free parameter, the scalaron mass $m$, and the value $m \simeq 4 \times 10^{-3}\, \text{eV}$ is required to explain the observed dark-matter abundance. Larger mass values would be required in scenarios where the scalaron is excited before the electroweak crossover., 13 pages, references added, discussion extended, version accepted for publication in Phys. Lett. B