Search

We find statistical evidence for a mismatch between the (global) spatial curvature parameter $K$ in the geodesic equation for incoming photons, and the corresponding parameter in the Friedmann equation that determines the time evolution of the background spacetime and its perturbations. The mismatch hereafter referred to as `curvature-slip' is especially evident when the SH0ES prior on the current expansion rate is assumed. This result is based on joint analyses of cosmic microwave background (CMB) observations with the PLANCK satellite (P18), first year results of the Dark Energy Survey (DES), Baryonic Oscillation (BAO) data, and - at a lower level of significance - also on Pantheon SNIa (SN) catalog. For example, the betting odds against the Null Hypothesis are greater than $10^7$:1, 1400:1 and 1000:1 when P18+SH0ES, P18+DES+SH0ES, and P18+BAO+SH0ES, respectively, are considered. Datasets involving SNIa weaken this curvature slip considerably. Notably, even when the SH0ES prior is not imposed the betting odds for the rejection of the Null Hypothesis are 70:1 and 160:1 in cases where P18+DES and P18+BAO are considered. When the SH0ES prior is imposed, global fit of the modified model (that allows for a nonvanishing `curvature slip') strongly outperforms that of $\Lambda$CDM as is manifested by significant Deviance Information Criterion (DIC) gains, ranging between 7 and 23, depending on the dataset combination considered. Even in comparison to K$\Lambda$CDM the proposed model results in significant, albeit smaller, DIC gains when SN data are excluded. Our finding could possibly be interpreted as an inherent inconsistency between the (idealized) maximally symmetric nature of the FRW metric, and the dynamical evolution of the GR-based homogeneous and isotropic $\Lambda$CDM model (abridged), Comment: Submitted. Comments are welcome

Spiral galaxies M31 and M33 are Fermi/LAT-detected gamma-ray sources. We model the broadband non-thermal (NT) emission of the central region of M31 (R < 5.5 kpc) and of the disk of M33 (R ~ 9 kpc). For either galaxy, we self-consistently model the broadband SED of the diffuse NT emission based on published radio and gamma-ray data. All relevant radiative processes involving relativistic and thermal electrons (synchrotron, Compton scattering, bremsstrahlung, and free-free emission and absorption), along with relativistic protons (neutral-pion decay following interaction with thermal protons), are considered, using exact emissivity formulae. We also use the Fermi/LAT validated gamma-ray emissivities for pulsars. We find that, in both sources, the radio emission is composed of primary and secondary electron synchrotron and thermal bremsstrahlung. The M33 gamma-ray emission appears to be mainly hadronic, similar to the Magellanic Clouds (Persic & Rephaeli 2022). In contrast, we find suggestions of a more complex situation in the central region of M31, whose emission could be a mix of pulsar emission and hadronic emission, with the latter possibly originating from both the disk and the vicinity of the nuclear black hole. The alternative modelling of the spectra of M31 and M33 is motivated by the different hydrogen distribution in the two galaxies: the hydrogen deficiency in the central region of M31 partially unveils emissions from the nuclear BH and the pulsar population in the bulge and inner disk. If this were to be the case in M33 as well, these emissions would be outshined by diffuse pionic emission originating within the flat central-peak gas distribution in M33., Comment: 19 pages, 4 figures, 2 appendices. References updated, typos corrected

Detailed modelling of the spectro-spatial distributions of energetic electrons and protons in galactic discs and haloes of starburst galaxies (SBGs) is needed in order to follow their interactions with the magnetized interstellar medium and radiation fields, determine their radiative yields, and for estimating their residual spectral densities in intergalactic environments. We have developed a semi-analytical approach for calculating the particle spectro-spatial distributions in the disc and halo based on a diffusion model for particle propagation from acceleration sites in the central SB and disc regions, including all their relevant interaction modes. Important overall normalization of our models is based on previous modelling of the Galactic disc (with the GALPROP code), scaled to the higher star-formations rate in NGC253, and on spatially resolved radio measurements of the central SB and disc. These provide the essential input for determining the particle distributions and their predicted radiative yields in the outer disc and inner halo for a range of values of the key parameters that affect diffusion rate and energy losses. Results of our work clearly indicate that quantitative description of non-thermal emission in SBGs has to be based on modelling of the particle distributions in the entire disc, not just the central SB region., Comment: 11 pages, 3 figures

The Magellanic Clouds, two dwarf galaxy companions to the Milky Way, are among the Fermi Large Area Telescope (LAT) brightest gamma-ray sources. Aiming at a comprehensive modeling of the non-thermal electromagnetic and neutrino emission in both Clouds, we self-consistently model the radio and gamma-ray spectral energy distribution from their disks based on recently published Murchison Widefield Array and Fermi/LAT data. All relevant radiative processes involving relativistic and thermal electrons (synchrotron, Compton scattering, and bremsstrahlung) and relativistic protons (neutral-pion decay following interaction with thermal protons) are considered, using exact emission formulae. Our joint spectral analyses indicate that radio emission in the Clouds has both primary and secondary electron synchrotron and thermal bremsstrahlung origin, whereas gamma rays originate mostly from neutral-pion decay with some contributions from relativistic bremsstrahlung and Compton scattering off starlight. The proton spectra in both galaxies are modeled as power laws in energy with similar spectral indices, ~2.4, and energy densities, ~1 eV/cm3. The predicted 0.1-10 GeV neutrino flux is too low for detection by current and upcoming experiments. Our analyses confirm earlier suggestions of a largely hadronic origin of the gamma-ray emission in both Magellanic Clouds., Comment: A&A, in press; 11 pages, 7 figures

This is a solicited whitepaper for the Snowmass 2021 community planning exercise. The paper focuses on measurements and science with the Cosmic Microwave Background (CMB). The CMB is foundational to our understanding of modern physics and continues to be a powerful tool driving our understanding of cosmology and particle physics. In this paper, we outline the broad and unique impact of CMB science for the High Energy Cosmic Frontier in the upcoming decade. We also describe the progression of ground-based CMB experiments, which shows that the community is prepared to develop the key capabilities and facilities needed to achieve these transformative CMB measurements., Comment: contribution to Snowmass 2021

The advent of new all-sky radio surveys such as the VLA Sky Survey (VLASS) and the Rapid ASKAP Continuum Survey (RACS), performed with the latest generation radio telescopes, is opening new possibilities on the classification and study of extragalactic $\gamma$-ray sources, specially the underrepresented ones like radio galaxies. In particular, the enhanced sensitivity (sub-mJy level) and resolution (a few arcsec) provides a better morphological and spectral classification. In this work, we present the reclassification of a Fermi-LAT source as a new FRII radio galaxy from the INTEGRAL sample found to emit at GeV energies. Through a broad-band spectral fitting from radio to $\gamma$-ray, we find that the commonly invoked jet contribution is not sufficient to account for the observed $\gamma$-ray flux. Our modeling suggests that the observed emission could mainly originate in the lobes (rather than in the radio core) by inverse Compton scattering of radio-emitting electrons off the ambient photon fields. In addition, we cross-correlated the latest generation radio surveys with a list of Fermi-LAT objects from the literature considered to be candidate misaligned AGN, finding four new radio galaxies with a double-lobed morphology. Additional four objects could be classified as such thanks to previous studies in the literature, for a total of nine new radio galaxies with GeV emission presented in this work. We foresee that further objects of this class might be found in the near future with the advent of the Square Kilometer Array (SKA), populating the GeV sky., Comment: Submitted to MNRAS

The cosmic microwave background (CMB) temperature, $T$, surely the most precisely measured cosmological parameter, has been inferred from {\it local} measurements of the blackbody spectrum to an exquisite precision of 1 part in $\sim 4700$. On the other hand, current precision allows inference of other basic cosmological parameters at the $\sim 1\%$ level from CMB power spectra, galaxy correlation and lensing, luminosity distance measurements of supernovae, as well as other cosmological probes. A basic consistency check of the standard cosmological model is an independent inference of $T$ at recombination. In this work we first use the recent Planck data, supplemented by either the first year data release of the dark energy survey (DES), baryon acoustic oscillations (BAO) data, and the Pantheon SNIa catalog, to extract $T$ at the $\sim 1\%$ precision level. We then explore correlations between $T$, the Hubble parameter, $H_{0}$, and the global spatial curvature parameter, $\Omega_{k}$. Our parameter estimation indicates that imposing the local constraint from the SH0ES experiment on $H_{0}$ results in significant statistical preference for departure at recombination from the locally inferred $T$. However, only moderate evidence is found in this analysis for tension between local and cosmological estimates of $T$, if the local constraint on $H_{0}$ is relaxed. All other dataset combinations that include the CMB with either BAO, SNIa, or both, disfavor the addition of a new free temperature parameter even in the presence of the local constraint on $H_{0}$. Analysis limited to the Planck dataset suggests the temperature at recombination was higher than expected at recombination at the $\gtrsim 95\%$ confidence level if space is globally flat., Comment: 12 pages, 3 figures. Abbreviated abstract. Accepted for publication in PRD

Recent analyses of the broad spectral energy distributions (SED) of extensive lobes of local radio-galaxies have confirmed the leptonic origin of their Fermi/LAT gamma-ray emission, significantly constraining the level of hadronic contribution. SED of distant (D > 125 Mpc) radio-galaxy lobes are currently limited to the radio and X-ray bands, hence give no information on the presence of non-thermal (NT) protons but are adequate to describe the properties of NT electrons. Modeling lobe radio and X-ray emission in 3C 98, Pictor A, DA 240, Cygnus A, 3C 326, and 3C 236, we fully determine the properties of intra-lobe NT electrons and estimate the level of the related gamma-ray emission from Compton scattering of the electrons off the superposed Cosmic Microwave Background, Extragalactic Background Light, and source-specific radiation fields., Comment: 8 pages, 4 figures; a few typos corrected; several missing references added

Radio and $\gamma$-ray measurements of large lobes of several radio galaxies provide adequate basis for determining whether emission in these widely separated spectral regions is largely by energetic electrons. This is very much of interest as there is of yet no unequivocal evidence for a significant energetic proton component to account for $\gamma$-ray emission by neutral pion decay. A quantitative assessment of the proton spectral distribution necessitates full accounting of the local and background radiation fields in the lobes; indeed, doing so in our recent analysis of the spectral energy distribution of the Fornax A lobes considerably weakened previous conclusions on the hadronic origin of the emission measured by the Fermi satellite. We present the results of similar analyses of the measured radio, X-ray and $\gamma$-ray emission from the lobes of Centaurus A, Centaurus B, and NGC 6251. The results indicate that the measured $\gamma$-ray emission from these lobes can be accounted for by Compton scattering of the radio-emitting electrons off the superposed radiation fields in the lobes; consequently, we set upper bounds on the energetic proton contents of the lobes., Comment: MNRAS, in press. 10 pages, 7 figures

The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4., Comment: Astro2020 Decadal Project Whitepaper. arXiv admin note: text overlap with arXiv:1808.07445

Starburst galaxies and star-forming active galactic nuclei (AGN) are among the candidate sources thought to contribute appreciably to the extragalactic gamma-ray and neutrino backgrounds. NGC 1068 is the brightest of the star-forming galaxies found to emit gamma rays from 0.1 to 50 GeV. Precise measurements of the high-energy spectrum are crucial to study the particle accelerators and probe the dominant emission mechanisms. We have carried out 125 hours of observations of NGC 1068 with the MAGIC telescopes in order to search for gamma-ray emission in the very high energy band. We did not detect significant gamma-ray emission, and set upper limits at 95\% confidence level to the gamma-ray flux above 200 GeV f<5.1x10^{-13} cm^{-2} s ^{-1} . This limit improves previous constraints by about an order of magnitude and allows us to put tight constraints on the theoretical models for the gamma-ray emission. By combining the MAGIC observations with the Fermi-LAT spectrum we limit the parameter space (spectral slope, maximum energy) of the cosmic ray protons predicted by hadronuclear models for the gamma-ray emission, while we find that a model postulating leptonic emission from a semi-relativistic jet is fully consistent with the limits. We provide predictions for IceCube detection of the neutrino signal foreseen in the hadronic scenario. We predict a maximal IceCube neutrino event rate of 0.07 yr^{-1}., Comment: 11 pages, 4 figures, accepted for publication in The Astrophysical Journal

Quantitative modeling of the spectro-spatial distributions of energetic electrons and protons in galactic halos is needed in order to determine their interactions with the local plasma and radiation fields, and also to estimate their residual spectral densities in intracluster and intergalactic environments. We develop a semi-analytic approach for calculating the particle distributions in the halo based on a detailed diffusion model for particle propagation from acceleration sites and interactions in the galactic disk. Important overall normalization of our models is based on results from detailed modeling in the Galactic disk with the GALPROP code. This provides the essential input for determining particle distributions in the outer disk, which are used as source terms for calculating the distributions in the extensive halo for a range of values of key parameters affecting energy losses and propagation mode. Our modeling approach is applied to the two edge-on star-forming galaxies NGC 4631 and NGC 4666, for which recent mapping of radio emission in the inner halo provides the required overall normalization. We predict the levels and spatial profiles of radio, X-ray, and gamma-ray emission in the halos of these galaxies. Our quantitative modeling enables us to estimate the total calorimetric efficiencies of electrons and protons in star-forming galaxies, and to predict their residual spectral distributions in the outer halo and intergalactic space., Comment: 12 pages, 4 figures

Current measurements of the spectral energy distribution in radio, X-and-gamma-ray provide a sufficiently wide basis for determining basic properties of energetic electrons and protons in the extended lobes of the radio galaxy Fornax A. Of particular interest is establishing observationally, for the first time, the level of contribution of energetic protons to the extended emission observed by the Fermi satellite. Two recent studies concluded that the observed gamma-ray emission is unlikely to result from Compton scattering of energetic electrons off the optical radiation field in the lobes, and therefore that the emission originates from decays of neutral pions produced in interactions of energetic protons with protons in the lobe plasma, implying an uncomfortably high proton energy density. However, our exact calculation of the emission by energetic electrons in the magnetized lobe plasma leads to the conclusion that all the observed emission can, in fact, be accounted for by energetic electrons scattering off the ambient optical radiation field, whose energy density (which, based on recent observations, is dominated by emission from the central galaxy NGC 1316) we calculate to be higher than previously estimated., Comment: Corrected typos; added reference; revised figure; result unchanged

The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes (SATs) and one large-aperture 6-m telescope (LAT), with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The SATs will target the largest angular scales observable from Chile, mapping ~10% of the sky to a white noise level of 2 $\mu$K-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma(r)=0.003$. The LAT will map ~40% of the sky at arcminute angular resolution to an expected white noise level of 6 $\mu$K-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the LSST sky region and partially with DESI. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources., Comment: This paper presents an overview of the Simons Observatory science goals, details about the instrument will be presented in a companion paper. The author contribution to this paper is available at https://simonsobservatory.org/publications.php (Abstract abridged) -- matching version published in JCAP

Nearby galaxy clusters can potentially induce sub-microkelvin polarization signals in the cosmic microwave background (CMB) at characteristic scales of a few arcminutes. We explore four such polarization signals induced in a rich nearby fiducial cluster and calculate the likelihood of their detection by a telescope project with capabilities such as those of the Simons Observatory (SO). In our feasibility analysis, we include instrumental noise, primordial CMB anisotropy, statistical thermal Sunyaev-Zeldovich (SZ) cluster signal, and point source confusion, assuming a few percent of the nominal telescope observation time of an SO-like project. Our analysis indicates that the thermal SZ intensity can be sensitively mapped in rich nearby clusters and that the kinematic SZ intensity can be measured with high statistical significance toward a fast moving nearby cluster. The detection of polarized SZ signals will be quite challenging but could still be feasible toward several very rich nearby clusters with very high SZ intensity. The polarized SZ signal from a sample of ~20 clusters can be statistically detected at S/N~3, if observed for several months., Comment: 11 pages, 4 figures

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA., Comment: Published on Journal of High Energy Astrophysics (Elsevier)

As energetic particles diffuse out of radio and star-forming galaxies (SFGs), their intracluster density builds up to a level that could account for a substantial part or all the emission from a radio halo. We calculate the particle time-dependent, spectro-spatial distributions from a solution of a diffusion equation with radio galaxies as sources of electrons, and SFGs as sources of both electrons and protons. Whereas strong radio galaxies are typically found in the cluster (e.g., Coma) core, the fraction of SFGs increases with distance from the cluster center. Scaling particle escape rates from their sources to the reasonably well determined Galactic rates, and for realistic gas density and magnetic field spatial profiles, we find that predicted spectra and spatial profiles of radio emission from primary and secondary electrons are roughly consistent with those deduced from current measurements of the Coma halo (after subtraction of emission from the relic Coma A). Nonthermal X-ray emission is predicted to be mostly by Compton scattering of electrons from radio galaxies off the CMB, whereas $\gamma$-ray emission is primarily from the decay of neutral pions produced in interactions of protons from SFGs with protons in intracluster gas., Comment: 6 pages; Rapid Communication, Physical Review D

We present results from {\gamma}-ray observations of the Coma cluster incorporating 6 years of Fermi-LAT data and the newly released {\emph{Pass 8}} event-level analysis. Our analysis of the region reveals low-significance residual structures within the virial radius of the cluster that are too faint for a detailed investigation with the current data. Using a likelihood approach that is free of assumptions on the spectral shape we derive upper limits on the {\gamma}-ray flux that is expected from energetic particle interactions in the cluster. We also consider a benchmark spatial and spectral template motivated by models in which the observed radio halo is mostly emission by secondary electrons. In this case, the median expected and observed upper limits for the flux above 100 MeV are $1.7\times10^{-9}\,\mathrm{ph\,cm^{-2}\,s^{-1}}$ and $5.2\times10^{-9}\,\mathrm{ph\,cm^{-2}\,s^{-1}}$ respectively (the latter corresponds to residual emission at the level of 1.8{\sigma}). These bounds are comparable to or higher than predicted levels of hadronic gamma-ray emission in cosmic-ray models with or without reacceleration of secondary electrons, although direct comparisons are sensitive to assumptions regarding the origin and propagation mode of cosmic rays and magnetic field properties. The minimal expected {\gamma}-ray flux from radio and star-forming galaxies within the Coma cluster is roughly an order of magnitude below the median sensitivity of our analysis., Comment: 9 pages, 4 figures, 2 tables, minor revised version accepted for publication; corresponding authors: S. Zimmer, J. Conrad, O. Reimer & Y. Rephaeli

Energy densities of relativistic electrons and protons in extended galactic and intracluster regions are commonly determined from spectral radio and (rarely) $\gamma$-ray measurements. The time-independent particle spectral density distributions are commonly assumed to have a power-law (PL) form over the relevant energy range. A theoretical relation between energy densities of electrons and protons is usually adopted, and energy equipartition is invoked to determine the mean magnetic field strength in the emitting region. We show that for typical conditions, in both star-forming and starburst galaxies, these estimates need to be scaled down substantially due to significant energy losses that (effectively) flatten the electron spectral density distribution, resulting in a much lower energy density than deduced when the distribution is assumed to have a PL form. The steady-state electron distribution in the nuclear regions of starburst galaxies is calculated by accounting for Coulomb, bremsstrahlung, Compton, and synchrotron losses; the corresponding emission spectra of the latter two processes are calculated and compared to the respective PL spectra. We also determine the proton steady-state distribution by taking into account Coulomb and pion production losses, and briefly discuss implications of our steady-state particle spectra for estimates of proton energy densities and magnetic fields., Comment: 8 pages, 4 figures

The mass fraction of hot gas in clusters is a basic quantity whose level and dependence on the cluster mass and redshift are intimately linked to all cluster X-ray and SZ measures. Modeling the evolution of the gas fraction is clearly a necessary ingredient in the description of the hierarchical growth of clusters through mergers of subclumps and mass accretion on the one hand, and the dispersal of gas from the cluster galaxies by tidal interactions, galactic winds, and ram pressure stripping on the other hand. A reasonably complete description of this evolution can only be given by very detailed hydrodynamical simulations, which are, however, resource-intensive, and difficult to implement in the mapping of parameter space. A much more practical approach is the use of semi-analytic modeling that can be easily implemented to explore a wide range of parameters. We present first results from a simple model that describes the build up of the gas mass fraction in clusters by following the overall impact of the above processes during the merger and accretion history of each cluster in the ensemble. Acceptable ranges for model parameters are deduced through comparison with results of X-ray observations. Basic implications of our work for modeling cluster statistical properties, and the use of these properties in joint cosmological data analyses, are discussed., Comment: 10 pages, 6 figures, updated to match MNRAS accepted version

The expected sensitivity of cluster SZ number counts to neutrino mass in the sub-eV range is assessed. We find that from the ongoing {\it Planck}/SZ measurements the (total) neutrino mass can be determined at a (1-sigma) precision of 0.06 eV, if the mass is in the range 0.1-0.3 eV, and the survey detection limit is set at the 5-sigma significance level. The mass uncertainty is predicted to be lower by a factor ~2/3, if a similar survey is conducted by a cosmic-variance-limited experiment, a level comparable to that projected if CMB lensing extraction is accomplished with the same experiment. At present, the main uncertainty in modeling cluster statistical measures reflects the difficulty in determining the mass function at the high-mass end., Comment: 8 pages, Proceedings of the 13th Marcel Grossmann Meeting

A basic quantity in the characterization of relativistic particles is the proton-to-electron (p/e) energy density ratio. We derive a simple approximate expression suitable to estimate this quantity, U_p/U_e = (m_p/m_e)^(3-q)/2, valid when a nonthermal `gas' of these particles is electrically neutral and the particles' power-law spectral indices are equal -- e.g., at injection. This relation partners the well-known p/e number density ratio at 1 GeV, i.e. N_p/N_e = (m_p/m_e)^{(q-1)/2}., Comment: 4 pages; to be published in Proc. of MGM13 (13th Marcel Grossmann Meeting -- Stockholm July 1-7, 2012)

We show that the fast moving component of the "bullet cluster" (1E0657-56) can induce potentially resolvable redshift differences between multiply-lensed images of background galaxies. The moving cluster effect can be expressed as the scalar product of the lensing deflection angle with the tangential velocity of the mass components, and it is maximal for clusters colliding in the plane of the sky with velocities boosted by their mutual gravity. The bullet cluster is likely to be the best candidate for the first measurement of this effect due to the large collision velocity and because the lensing deflection and the cluster fields can be calculated in advance. We derive the deflection field using multiply-lensed background galaxies detected with the Hubble Space Telescope. The velocity field is modeled using self-consistent N-body/hydrodynamical simulations constrained by the observed X-ray and gravitational lensing features of this system. We predict that the triply-lensed images of systems "G" and "H" straddling the critical curve of the bullet component will show the largest frequency shifts up to ~0.5 km/sec. This is within the range of the Atacama Large Millimeter/sub-millimeter Array (ALMA) for molecular emission, and is near the resolution limit of the new generation high-throughput optical-IR spectrographs. A detection of this effect measures the tangential motion of the subclusters directly, thereby clarifying the tension with LCDM, which is inferred from gas motion less directly. This method may be extended to smaller redshift differences using the Ly-alpha forest towards QSOs lensed by more typical clusters of galaxies. More generally, the tangential component of the peculiar velocities of clusters derived by our method complements the radial component determined by the kinematic SZ effect, providing a full 3-dimensional description of velocities., Comment: 12 pages, 6 Figures and 2 Tables, accepted for publication in the Astrophysical Journal

The detections of high-energy gamma-ray emission from the nearby starburst galaxies M82 & NGC253, and other local group galaxies, broaden our knowledge of star-driven nonthermal processes and phenomena in non-AGN star-forming galaxies. We review basic aspects of the related processes and their modeling in starburst galaxies. Since these processes involve both energetic electrons and protons accelerated by SN shocks, their respective radiative yields can be used to explore the SN-particle-radiation connection. Specifically, the relation between SN activity, energetic particles, and their radiative yields, is assessed through respective measures of the particle energy density in several star-forming galaxies. The deduced energy densities range from O(0.1) eV/cm^3 in very quiet environments to O(100) eV/cm^3 in regions with very high star-formation rates., Comment: 17 pages, 5 figures, to be published in Astrophysics and Space Science Proceedings

In this article we focus on the astrophysical results and the related cosmological implications derived from recent microwave surveys, with emphasis to those coming from the Planck mission. We critically discuss the impact of systematics effects and the role of methods to separate the cosmic microwave background signal from the astrophysical emissions and each different astrophysical component from the others. We then review of the state of the art in diffuse emissions, extragalactic sources, cosmic infrared back- ground, and galaxy clusters, addressing the information they provide to our global view of the cosmic structure evolution and for some crucial physical parameters, as the neutrino mass. Finally, we present three different kinds of scientific perspectives for fundamental physics and cosmology offered by the analysis of on-going and future cosmic microwave background projects at different angular scales dedicated to anisotropies in total intensity and polarization and to absolute temperature., Comment: 38 pages, 9 figures. Accepted for publications on IJMPD. Based on talks presented at the Thirteenth Marcel Grossmann Meeting on General Relativity, Stockholm, July 2012 (Cosmic Background (CB) sessions). It will appear also on corresponding proceedings edited by Kjell Rosquist, Robert T. Jantzen, Remo Ruffini, World Scientific, Singapore, 2013

It is known that modeling uncertainties and astrophysical foregrounds can potentially introduce appreciable bias in the deduced values of cosmological parameters. While it is commonly assumed that these uncertainties will be accounted for to a sufficient level of precision, the level of bias has not been properly quantified in most cases of interest. We show that the requirement that the bias in derived values of cosmological parameters does not surpass nominal statistical error, translates into a maximal level of overall error $O(N^{-1/2})$ on $|\Delta P(k)|/P(k)$ and $|\Delta C_{l}|/C_{l}$, where $P(k)$, $C_{l}$, and $N$ are the matter power spectrum, angular power spectrum, and number of (independent Fourier) modes at a given scale $l$ or $k$ probed by the cosmological survey, respectively. This required level has important consequences on the precision with which cosmological parameters are hoped to be determined by future surveys: In virtually all ongoing and near future surveys $N$ typically falls in the range $10^{6}-10^{9}$, implying that the required overall theoretical modeling and numerical precision is already very high. Future redshifted-21-cm observations, projected to sample $\sim 10^{14}$ modes, will require knowledge of the matter power spectrum to a fantastic $10^{-7}$ precision level. We conclude that realizing the expected potential of future cosmological surveys, which aim at detecting $10^{6}-10^{14}$ modes, sets the formidable challenge of reducing the overall level of uncertainty to $10^{-3}-10^{-7}$., Comment: 19 pages, revised version accepted for publication in JCAP

Hot gas in filamentary structures induces CMB aniostropy through the SZ effect. Guided by results from N-body simulations, we model the morphology and gas properties of filamentary gas and determine the power spectrum of the anisotropy. Our treatment suggests that power levels can be an appreciable fraction of the cluster contribution at multipoles $\ell\lesssim 1500$. Its spatially irregular morphology and larger characteristic angular scales can help to distinguish this SZ signature from that of clusters. In addition to intrinsic interest in this most extended SZ signal as a probe of filaments, its impact on cosmological parameter estimation should also be assessed. We find that filament `noise' can potentially bias determination of $A_s$, $n_s$, and $w$ (the normalization of the primordial power spectrum, the scalar index, and the dark energy equation of state parameter, respectively) by more than the nominal statistical uncertainty in Planck SZ survey data. More generally, when inferred from future optimal cosmic-variance-limited CMB experiments, we find that virtually all parameters will be biased by more than the nominal statistical uncertainty estimated for these next generation CMB experiments., Comment: 14 pages, 2 figures, accepted for publication in JCAP

The Jaynes-Cummings (JC) system, which describes the interaction between a cavity and a two-level atom, is one of the most important systems in quantum optics. We obtain analytic solutions for the one- and two-photon transport in a waveguide side-coupled to the JC system using input-output formalism in Fock space. With these results, we discuss the conditions under which the JC system functions as a photon switch for waveguide photons in both the strong and weak coupling regimes, Comment: 9 pages, 4 figures

We study stimulated emission from an excited two-level atom coupled to a waveguide containing an incident single-photon pulse. We show that the strong photon correlation, as induced by the atom, plays a very important role in stimulated emission. Additionally, the temporal duration of the incident photon pulse is shown to have a marked effect on stimulated emission and atomic lifetime., Comment: 6 pages, 3 figures

The common approach to compute the cosmic-ray distribution in an starburst galaxy or region is equivalent to assume that at any point within that environment, there is an accelerator inputing cosmic rays at a reduced rate. This rate should be compatible with the overall volume-average injection, given by the total number of accelerators that were active during the starburst age. These assumptions seem reasonable, especially under the supposition of an homogeneous and isotropic distribution of accelerators. However, in this approach the temporal evolution of the superposed spectrum is not explicitly derived; rather, it is essentially assumed ab-initio. Here, we test the validity of this approach by following the temporal evolution and spatial distribution of the superposed cosmic-ray spectrum and compare our results with those from theoretical models that treat the starburst region as a single source. In the calorimetric limit (with no cosmic-ray advection), homogeneity is reached (typically within 20%) across most of the starburst region. However, values of center-to-edge intensity ratios can amount to a factor of several. Differences between the common homogeneous assumption for the cosmic-ray distribution and our models are larger in the case of two-zone geometries, such as a central nucleus with a surrounding disc. We have also found that the decay of the cosmic-ray density following the duration of the starburst process is slow, and even approximately 1 Myr after the burst ends (for a gas density of 35 cm-3) it may still be within an order of magnitude of its peak value. Based on our simulations, it seems that the detection of a relatively hard spectrum up to the highest gamma-ray energies from nearby starburst galaxies favors a relatively small diffusion coefficient (i.e., long diffusion time) in the region where most of the emission originates., Comment: To appear in MNRAS

Statistical measures of galaxy clusters are sensitive to neutrino masses in the sub-eV range. We explore the possibility of using cluster number counts from the ongoing PLANCK/SZ and future cosmic-variance-limited surveys to constrain neutrino masses from CMB data alone. The precision with which the total neutrino mass can be determined from SZ number counts is limited mostly by uncertainties in the cluster mass function and intracluster gas evolution; these are explicitly accounted for in our analysis. We find that projected results from the PLANCK/SZ survey can be used to determine the total neutrino mass with a (1\sigma) uncertainty of 0.06 eV, assuming it is in the range 0.1-0.3 eV, and the survey detection limit is set at the 5\sigma significance level. Our results constitute a significant improvement on the limits expected from PLANCK/CMB lensing measurements, 0.15 eV. Based on expected results from future cosmic-variance-limited (CVL) SZ survey we predict a 1\sigma uncertainty of 0.04 eV, a level comparable to that expected when CMB lensing extraction is carried out with the same experiment. A few percent uncertainty in the mass function parameters could result in up to a factor \sim 2-3 degradation of our PLANCK and CVL forecasts. Our analysis shows that cluster number counts provide a viable complementary cosmological probe to CMB lensing constraints on the total neutrino mass., Comment: Replaced with a revised version to match the MNRAS accepted version. arXiv admin note: text overlap with arXiv:1009.4110

We have recently developed an extended merger-tree model that efficiently follows hierarchical evolution of galaxy clusters and provides a quantitative description of both their dark matter and gas properties. We employed this diagnostic tool to calculate the thermal SZ power spectrum and cluster number counts, accounting explicitly for uncertainties in the relevant statistical and intrinsic cluster properties, such as the halo mass function and the gas equation of state. Results of these calculations are compared with those obtained from a direct analytic treatment and from hydrodynamical simulations. We show that under certain assumptions on the gas mass fraction our results are consistent with the latest SPT measurement. Our approach can be particularly useful in predicting cluster number counts and their dependence on cluster and cosmological parameters., Comment: 8 pages, 7 figures, accepted for publication in MNRAS

The energy density of energetic protons, U_p, in several nearby starburst nuclei (SBNs) has been directly deduced from gamma-ray measurements of the radiative decay of neutral pions produced in interactions with ambient protons. Lack of sufficient sensitivity and spatial resolution makes this direct deduction unrealistic in the foreseeable future for even moderately distant SBNs. A more viable indirect method for determining U_p in star-forming galaxies is to use its theoretically based scaling to the energy density of energetic electrons, U_e, which can be directly deduced from radio synchrotron and possibly also nonthermal hard X-ray emission. In order to improve the quantitative basis and diagnostic power of this leptonic method we reformulate and clarify its main aspects. Doing so we obtain a basic expression for the ratio U_p/U_e in terms of the proton and electron masses and the power-law indices that characterize the particle spectral distributions in regions where the total particle energy density is at equipartition with that of the mean magnetic field. We also express the field strength and the particle energy density in the equipartition region in terms of the region's size, mean gas density, IR and radio fluxes, and distance from the observer, and determine values of U_p in a sample of nine nearby and local SBNs., Comment: Published in: A&A, 567, A101 (2014) - 7 pages, 1 figure

We calculate the one- and two-photon scattering matrices of a pair of collocated non-identical two-level atoms coupled to a waveguide. We show that by proper choice of a two-photon input, the background fluorescence by the atoms may be completely quenched, as a result of quantum interference, and that when the atoms' detuning is smaller than their linewidths, extremely narrow fluorescence features emerge. Furthermore, the system emits a two-photon bound state which can display spatial oscillations/quantum beats, and can be tuned from bunched to anti-bunched statistics as the total photon energy is varied.

We show how to calculate the first- and second-order statistics of the scattered fields for an arbitrary intensity coherent state light field interacting with a two-level system in a waveguide geometry. Specifically, we calculate the resonance fluorescence from the qubit, using input-output formalism. We derive the transmission and reflection coefficients, and illustrate the bunching and anti-bunching of light that is scattered in the forward and backward directions, respectively. Our results agree with previous calculations on one- and two-photon scattering as well as those that are based on the master equation approach., Comment: 8 pages, 3 figures and supplementary material (Mathematica code). This is the published version: typos are fixed, conclusion section is expanded, references are updated

Extensive surveys of galaxy clusters motivate us to assess the likelihood of cluster-cluster lensing (CCL), namely, gravitational-lensing of a background cluster by a foreground cluster. We briefly describe the characteristics of CCLs in optical, X-ray and SZ measurements, and calculate their predicted numbers for $\Lambda$CDM parameters and a viable range of cluster mass functions and their uncertainties. The predicted number of CCLs in the strong-lensing regime varies from several ($<10$) to as high as a few dozen, depending mainly on whether lensing triaxiality bias is accounted for, through the c-M relation. A much larger number is predicted when taking into account also CCL in the weak-lensing regime. In addition to few previously suggested CCLs, we report a detection of a possible CCL in A383, where background candidate high-$z$ structures are magnified, as seen in deep Subaru observations., Comment: 9 pages, 5 figures, submitted to MNRAS

We report statistical results for dark matter (DM) velocity anisotropy, \beta, from a sample of some 6000 cluster-size halos (at redshift zero) identified in a \Lambda CDM hydrodynamical adaptive mesh refinement simulation performed with the Enzo code. These include profiles of \beta\ in clusters with different masses, relaxation states, and at several redshifts, modeled both as spherical and triaxial DM configurations. Specifically, although we find a large scatter in the DM velocity anisotropy profiles of different halos (across elliptical shells extending to at least ~$1.5 r_{vir}$), universal patterns are found when these are averaged over halo mass, redshift, and relaxation stage. These are characterized by a very small velocity anisotropy at the halo center, increasing outward to about 0.27 and leveling off at about $0.2 r_{vir}$. Indirect measurements of the DM velocity anisotropy fall on the upper end of the theoretically expected range. Though measured indirectly, the estimations are derived by using two different surrogate measurements - X-ray and galaxy dynamics. Current estimates of the DM velocity anisotropy are based on very small cluster sample. Increasing this sample will allow testing theoretical predictions, including the speculation that the decay of DM particles results in a large velocity boost. We also find, in accord with previous works, that halos are triaxial and likely to be more prolate when unrelaxed, whereas relaxed halos are more likely to be oblate. Our analysis does not indicate that there is significant correlation (found in some previous studies) between the radial density slope, \gamma, and \beta\ at large radii, $0.3 r_{vir} < r < r_{vir}$., Comment: 12 pages, 17 figures, accepted to ApJ

We present results from strong-lens modelling of 10,000 SDSS clusters, to establish the universal distribution of Einstein radii. Detailed lensing analyses have shown that the inner mass distribution of clusters can be accurately modelled by assuming light traces mass, successfully uncovering large numbers of multiple-images. Approximate critical curves and the effective Einstein radius of each cluster can therefore be readily calculated, from the distribution of member galaxies and scaled by their luminosities. We use a subsample of 10 well-studied clusters covered by both SDSS and HST to calibrate and test this method, and show that an accurate determination of the Einstein radius and mass can be achieved by this approach "blindly", in an automated way, and without requiring multiple images as input. We present the results of the first 10,000 clusters analysed in the range $0.1=0.73^{+0.02}_{-0.03}$, $\sigma=0.316^{+0.004}_{-0.002}$, and with higher abundance of large $\theta_{e}$ clusters than predicted by $\Lambda$CDM. We visually inspect each of the clusters with $\theta_{e}>40 \arcsec$ ($z_{s}=2$) and find that $\sim20%$ are boosted by various projection effects detailed here, remaining with $\sim40$ real giant-lens candidates, with a maximum of $\theta_{e}=69\pm12 \arcsec$ ($z_{s}=2$) for the most massive candidate, in agreement with semi-analytic calculations. The results of this work should be verified further when an extended calibration sample is available., Comment: 18 pages, 15 figures, 1 table; V2 accepted to MNRAS, includes a significant revision, in particular a new discussion of the results

Little is known about the origin and basic properties of magnetic fields in clusters of galaxies. High conductivity in magnetized interstellar plasma suggests that galactic magnetic fields are (at least partly) ejected into intracluster (IC) space by the same processes that enrich IC gas with metals. We explore the dispersal of galactic fields by hydrodynamical simulations with our new {\em Enzo-Galcon} code, which is capable of tracking a large number galaxies during cluster assembly, and modeling the processes that disperse their interstellar media. Doing so we are able to describe the evolution of the mean strength of the field and its profile across the cluster. With the known density profile of dispersed gas and an estimated range of coherence scales, we predict the spatial distribution of Faraday rotation measure and find it to be consistent with observational data.

Clusters of galaxies are uniquely important cosmological probes of the evolution of the large scale structure, whose diagnostic power depends quite significantly on the ability to reliably determine their masses. Clusters are typically modeled as spherical systems whose intracluster gas is in strict hydrostatic equilibrium (i.e., the equilibrium gas pressure is provided entirely by thermal pressure), with the gravitational field dominated by dark matter, assumptions that are only rough approximations. In fact, numerical simulations indicate that galaxy clusters are typically triaxial, rather than spherical, and that turbulent gas motions (induced during hierarchical merger events) provide an appreciable pressure component. Extending our previous work, we present results of a joint analysis of X-ray, weak and strong lensing measurements of Abell 1689. The quality of the data allows us to determine both the triaxial shape of the cluster and the level of non-thermal pressure that is required if the intracluster gas is in hydrostatic equilibrium. We find that the dark matter axis ratios are 1.24 +/- 0.13 and 2.02 +/- 0.01 on the plane of the sky and along the line of sight, respectively, and that about 20% of the pressure is non-thermal. Our treatment demonstrates that the dynamical properties of clusters can be determined in a (mostly) bias-free way, enhancing the use of clusters as more precise cosmological probes., Comment: Accepted for publication in MNRAS

Adopting the convection-diffusion model for energetic electron and proton propagation, and accounting for all the relevant hadronic and leptonic processes, the steady-state energy distributions of these particles in the starburst galaxies M82 and NGC253 can be determined with a detailed numerical treatment. The electron distribution is directly normalized by the measured synchrotron radio emission from the central starburst region; a commonly expected theoretical relation is then used to normalize the proton spectrum in this region, and a radial profile is assumed for the magnetic field. The resulting radiative yields of electrons and protons are calculated: the predicted >100MeV and >100GeV fluxes are in agreement with the corresponding quantities measured with the orbiting Fermi telescope and the ground-based VERITAS and HESS Cherenkov telescopes. The cosmic-ray energy densities in central regions of starburst galaxies, as inferred from the radio and gamma-ray measurements of (respectively) non-thermal synchrotron and neutral-pion-decay emission, are U=O(100) eV/cm3, i.e. at least an order of magnitude larger than near the Galactic center and in other non-very-actively star-forming galaxies. These very different energy density levels reflect a similar disparity in the respective supernova rates in the two environments. A L(gamma) ~ SFR^(1.4) relationship is then predicted, in agreement with preliminary observational evidence., Comment: Invited talk at SciNeGHE2010 (8th Wotkshop on Science with the New Generation of High Energy Gamma-ray Experiments): Gamma-ray Astrophysics in the Multimessenger Context (Trieste, Sept.8-10, 2010)

We model the formation and evolution of galaxy clusters in the framework of an extended dark matter halo merger-tree algorithm that includes baryons and incorporates basic physical considerations. Our modified treatment is employed to calculate the probability density functions of the halo concentration parameter, intracluster gas temperature, and the integrated Comptonization parameter for different cluster masses and observation redshifts. Scaling relations between cluster mass and these observables are deduced that are somewhat different than previous results. Modeling uncertainties in the predicted probability density functions are estimated. Our treatment and the insight gained from the results presented in this paper can simplify the comparison of theoretical predictions with results from ongoing and future cluster surveys., Comment: 11 pages, 9 figures, submitted to MNRAS