We perform a series of controlled N-body simulations to study realizations of the recently discovered Antlia 2 galaxy in cold dark matter (CDM) and self-interacting dark matter (SIDM) scenarios. Our simulations contain six benchmark models, where we vary the initial halo concentration and the self-scattering cross section. We adopt well-motivated initial stellar and halo masses, and our fiducial orbit has a small pericenter. After evolving in the Milky Way's tidal field, the simulated galaxies experience significant mass loss and their stellar distributions expand accordingly. These tidal effects are more prominent if the initial halo concentration is lower and if the self-scattering cross section is larger. Our results show that Antlia 2-like galaxies could be realized in CDM if the halo concentration is low and the stellar distribution is diffuse at the infall time, while these conditions could be relaxed in SIDM. We also find all the simulated galaxies predict approximately the same stellar velocity dispersion after imposing selection criteria for stellar particles. This has important implications for testing dark matter models using tidally disturbed systems.

We propose two effective parameters that fully characterise galactic-scale structure formation at high redshifts ($z\gtrsim5$) for a variety of dark matter (DM) models that have a primordial cutoff in the matter power spectrum. Our description is within the recently proposed ETHOS framework and includes standard thermal Warm DM (WDM) and models with dark acoustic oscillations (DAOs). To define and explore this parameter space, we use high-redshift zoom-in simulations that cover a wide range of non-linear scales from those where DM should behave as CDM ($k\sim10\,h\,{\rm Mpc}^{-1}$), down to those characterised by the onset of galaxy formation ($k\sim500\,h\,{\rm Mpc}^{-1}$). We show that the two physically motivated parameters $h_{\rm peak}$ and $k_{\rm peak}$, the amplitude and scale of the first DAO peak, respectively, are sufficient to parametrize the linear matter power spectrum and classify the DM models as belonging to effective non-linear structure formation regions. These are defined by their relative departure from Cold DM ($k_{\rm peak}\rightarrow\infty$) and WDM ($h_{\rm peak}=0$) according to the non-linear matter power spectrum and halo mass function. We identify a region where the DAOs still leave a distinct signature from WDM down to $z=5$, while a large part of the DAO parameter space is shown to be degenerate with WDM. Our framework can then be used to seamlessly connect a broad class of particle DM models to their structure formation properties at high redshift without the need of additional $N$-body simulations.

We study evolution of self-interacting dark matter (SIDM) subhalos in the Milky Way (MW) tidal field. The interaction between the subhalos and the MW's tides lead to more diverse dark matter distribution in the inner region, compared to their cold dark matter counterparts. We test this scenario with two MW satellite galaxies, Draco and Fornax, opposite extremes in the inner dark matter content, and find that they can be accommodated within the SIDM model proposed to explain the diverse rotation curves of spiral galaxies in the field.

The nature of the dark matter can affect the collapse time of dark matter haloes, and can therefore be imprinted in observables such as the stellar population ages and star formation histories of dwarf galaxies. In this paper we use high resolution hydrodynamical simulations of Local Group-analogue (LG) volumes in cold dark matter (CDM), sterile neutrino warm dark matter (WDM) and self-interacting dark matter (SIDM) models with the EAGLE galaxy formation code to study how galaxy formation times change with dark matter model. We are able to identify the same haloes in different simulations, since they share the same initial density field phases. We show that the stellar mass varies systematically with resolution by over a factor of two, in a manner that depends on the final stellar mass. The evolution of the stellar populations in SIDM is largely identical to that of CDM, but in WDM early star formation is instead suppressed. The time at which LG haloes can begin to form stars through atomic cooling is delayed by $\sim$200~Myr in WDM models compared to CDM. 70~per~cent of WDM haloes of mass $>10^{8}M_{\odot}$ collapse early enough to form stars before $z=6$, compared to 90~per~cent of CDM and SIDM galaxies. It will be necessary to measure stellar ages of old populations to a precision of better than 100~Myr, and to address degeneracies with the redshift of reionization, in order to use these observables to distinguish between dark matter models.

We combine observations from the Atacama Large Millimeter/submillimeter Array (ALMA) and the NOrthern Extended Millimeter Array (NOEMA) to assess the redshift and to study the star formation conditions in AzTEC2: one of the brightest sub-millimeter galaxies (SMGs) in the COSMOS field ($S_{\rm 1.1mm}=10.5\pm1.4$mJy). Our high-resolution observations confirm that AzTEC2 splits into two components (namely AzTEC2-A and AzTEC2-B) for which we detect [C$\,$II] and $^{12}$CO(5$\to$4) line emission, implying a redshift of $4.626\pm0.001$ ($4.633\pm0.001$) for AzTEC2-A (AzTEC2-B) and ruling out previous associations with a galaxy at $z\sim1$. We use the $^{12}$CO(5$\to$4) line emission and adopt typical SMG-like gas excitation conditions to estimate the molecular gas mass, which is $M_{\rm gas}(\alpha_{\rm CO}/2.5)=2.1\pm0.4 \times10^{11}{\rm M}_\odot$ for AzTEC2-A, and a factor four lower for AzTEC2-B. With the infrared-derived star formation rate of AzTEC2-A ($1920\pm100 \,M_\odot{\rm \, yr}^{-1}$) and AzTEC2-B ($710\pm 35\,M_\odot{\rm \,yr}^{-1}$), they both will consume their current gas reservoir within $(30-200)\,$Myr. We find evidence of a rotation-dominated [C$\,$II] disk in AzTEC2-A, with a de-projected rotational velocity of $v_{\rm rot}(i=39^\circ)=660\pm130{\rm \,km\,s}^{-1}$, velocity dispersion $\lesssim100{\rm \,km\,s}^{-1}$, and dynamical mass of $M_{\rm dyn}(i=39^\circ)=2.6^{+1.2}_{-0.9}\times10^{11}\,M_\odot$. We propose that an elevated gas accretion rate from the cosmic web might be the main driver of the intense levels of star formation in AzTEC2-A, which might be further enhanced by gravitational torques induced by its minor companion (AzTEC2-B). These results strengthen the picture whereby the population of single-dish selected SMGs is rather heterogeneous, including a population of pairs of massive, highly-active galaxies in a pre-coalescence phase.

We report a ground-based detection of the [OI] 63-um line in a z=6.027 gravitationally lensed dusty star-forming galaxy (DSFG) G09.83808 using the APEX SEPIA 660 receiver, the first unambiguous detection of the [OI]63 line beyond redshift 3, and the first obtained from the ground. The [OI]63 line is robustly detected at 22$\pm$5 Jy km s$^{-1}$, corresponding to an intrinsic (de-lensed) luminosity of $(5.4\pm1.3)\times10^{9}$ L$_\odot$. With the [OI]63/[CII] luminosity ratio of 4, the [OI]63 line is the main coolant of the neutral gas in this galaxy, in agreement with model predictions. The high [OI]63 luminosity compensates for the pronounced [CII] deficit ([CII]/FIR$\simeq4\times10^{-4}$). Using photon-dominated region models, we derive a source-averaged gas density $n=10^{4.0}$ cm$^{-3}$, and far-UV field strength $G=10^4 G_0$, comparable to the z=2-4 DSFG population. If G09.83808 represents a typical high-redshift DSFG, the [OI]63 line from z=6 non-lensed DSFGs should be routinely detectable in ALMA Band 9 observations with $\sim$15 min on-source, opening a new window to study the properties of the earliest DSFGs.

Bose-Einstein Condensate Dark Matter (BECDM; also known as Fuzzy Dark Matter) is motivated by fundamental physics and has recently received significant attention as a serious alternative to the established Cold Dark Matter (CDM) model. We perform cosmological simulations of BECDM gravitationally coupled to baryons and investigate structure formation at high redshifts ($z \gtrsim 5$) for a boson mass $m=2.5\cdot 10^{-22}~{\rm eV}$, exploring the dynamical effects of its wavelike nature on the cosmic web and the formation of first galaxies. Our BECDM simulations are directly compared to CDM as well as to simulations where the dynamical quantum potential is ignored and only the initial suppression of the power spectrum is considered -- a Warm Dark Matter-like ("WDM") model often used as a proxy for BECDM. Our simulations confirm that "WDM" is a good approximation to BECDM on large cosmological scales even in the presence of the baryonic feedback. Similarities also exist on small scales, with primordial star formation happening both in isolated haloes and continuously along cosmic filaments; the latter effect is not present in CDM. Global star formation and metal enrichment in these first galaxies are delayed in BECDM/"WDM" compared to the CDM case: in BECDM/"WDM" first stars form at $z\sim 13$/$13.5$ while in CDM star formation starts at $z\sim 35$. The signature of BECDM interference, not present in "WDM", is seen in the evolved dark matter power spectrum: although the small scale structure is initially suppressed, power on kpc scales is added at lower redshifts. Our simulations lay the groundwork for realistic simulations of galaxy formation in BECDM.

We present a physical characterization of MMJ100026.36+021527.9 (a.k.a. ``MAMBO-9''), a dusty star-forming galaxy (DSFG) at $z=5.850\pm0.001$. This is the highest redshift unlensed DSFG (and fourth most distant overall) found to-date, and is the first source identified in a new 2mm blank-field map in the COSMOS field. Though identified in prior samples of DSFGs at 850$\mu$m-1.2mm with unknown redshift, the detection at 2mm prompted further follow-up as it indicated a much higher probability that the source was likely to sit at $z>4$. Deep observations from the Atacama Large Millimeter and submillimeter Array (ALMA) presented here confirm the redshift through the secure detection of $^{12}$CO($J\!=$6$\rightarrow$5) and p-H$_{2}$O(2$_{1,1}\!\rightarrow$2$_{0,2}$). MAMBO-9 is comprised of a pair of galaxies separated by 6kpc with corresponding star-formation rates of 590M$_\odot$yr$^{-1}$ and 220M$_\odot$yr$^{-1}$ total molecular hydrogen gas mass of (1.7$\pm$0.4)$\times10^{11}$M$_\odot$, dust mass of (1.3$\pm$0.3)$\times10^{9}$M$_\odot$ and stellar mass of (3.2$^{+1.0}_{-1.5}$)$\times10^{9}$M$_\odot$. The total halo mass, (3.3$\pm$0.8)$\times10^{12}$M$_\odot$, is predicted to exceed $>10^{15}$M$_\odot$ by $z=0$. The system is undergoing a merger-driven starburst which will increase the stellar mass of the system tenfold in $\tau_{\rm depl}=40-80$Myr, converting its large molecular gas reservoir (gas fraction of 96$^{+1}_{-2}$%) into stars. MAMBO-9 evaded firm spectroscopic identification for a decade, following a pattern that has emerged for some of the highest redshift DSFGs found. And yet, the systematic identification of unlensed DSFGs like MAMBO-9 is key to measuring the global contribution of obscured star-formation to the star-formation rate density at $z>4$, the formation of the first massive galaxies, and the formation of interstellar dust at early times ($<$1Gyr).

We present ALMA Band 6 (nu=233GHz, lambda=1.3mm) continuum observations towards 68 'normal' star-forming galaxies within two Coma-like progenitor structures at z=2.10 and 2.47, from which ISM masses are derived, providing the largest census of molecular gas mass in overdense environments at these redshifts. Our sample comprises galaxies with a stellar mass range of 1x10^9M_sun - 4x10^11M_sun with a mean M_*~6x10^10M_sun. Combining these measurements with multiwavelength observations and SED modeling, we characterize the gas mass fraction and the star formation efficiency, and infer the impact of the environment on galaxies' evolution. Most of our detected galaxies (~70%) have star formation efficiencies and gas fractions similar to those found for coeval field galaxies and in agreement with the field scaling relations. However, we do find that the proto-clusters contain an increased fraction of massive, gas-poor galaxies, with low gas fractions (f_gas<6-10%) and red rest-frame ultraviolet/optical colors typical of post-starburst and passive galaxies. The relatively high abundance of passive galaxies suggests an accelerated evolution of massive galaxies in proto-cluster environments. The large fraction of quenched galaxies in these overdense structures also implies that environmental quenching takes place during the early phases of cluster assembly, even before virialization. From our data, we derive a quenching efficiency of E_q~0.45 and an upper limit on the quenching timescale of T_q<1Gyr.

The rich information on (sub)millimeter dust continuum emission from distant galaxies in the public Atacama Large Millimeter/submillimeter Array (ALMA) archive is contained in thousands of inhomogeneous observations from individual PI-led programs. To increase the usability of these data for studies deepening our understanding of galaxy evolution, we have developed automated mining pipelines for the ALMA archive in the COSMOS field (A3COSMOS) that efficiently exploit the available information for large numbers of galaxies across cosmic time and keep the data products in sync with the increasing public ALMA archive: (a) a dedicated ALMA continuum imaging pipeline, (b) two complementary photometry pipelines for both blind source extraction and prior source fitting, (c) a counterpart association pipeline utilizing the multiwavelength data available (including quality assessment based on machine-learning techniques), (d) an assessment of potential (sub)millimeter line contribution to the measured ALMA continuum, and (e) extensive simulations to provide statistical corrections to biases and uncertainties in the ALMA continuum measurements. Application of these tools yields photometry catalogs with ~1000 (sub)millimeter detections (spurious fraction ~8%-12%) from over 1500 individual ALMA continuum images. Combined with ancillary photometric and redshift catalogs and the above quality assessments, we provide robust information on redshift, stellar mass, and star formation rate for ~700 galaxies at redshifts 0.5-6 in the COSMOS field (with undetermined selection function). The ALMA photometric measurements and galaxy properties are released publicly within our blind extraction, prior fitting, and galaxy property catalogs, plus the images. These products will be updated on a regular basis in the future.

An important aspect of quenching star formation is the removal of the cold interstellar medium (ISM; non-ionised gas and dust) from a galaxy. In addition, dust grains can be destroyed in a hot or turbulent medium. The adopted timescale of dust removal usually relies on uncertain theoretical estimates. It is tricky to track the dust removal, because usually dust is constantly replenished by consecutive generations of stars. Our objective is to measure observationally the timescale of dust removal. We here explore an approach to select galaxies which do have detectable amounts of dust and cold ISM but exhibit a low current dust production rate. Any decrease of the dust and gas content as a function of the age of such galaxies therefore must be attributed to processes governing the ISM removal. We used a sample of galaxies detected by Herschel in the far-infrared with visually assigned early-type morphology or spirals with red colours. We also obtained JCMT/SCUBA-2 observations for five of them. We discovered an exponential decline of the dust-to-stellar mass ratio with age, which we interpret as an evolutionary trend of dust removal from these galaxies. For the first time we directly measure the dust removal timescale in such galaxies to be tau=(2.5+-0.4) Gyr (the corresponding half-life time is (1.75+-0.25) Gyr). This quantity may be used in models in which it must be assumed a priori and cannot be derived. Any process which removes dust in these galaxies, such as dust grain destruction, cannot happen on shorter timescales. The timescale is comparable to the quenching timescales found in simulations for galaxies with similar stellar masses. The dust is likely of internal, not external origin. It was either formed in the past directly by supernovae, or from seeds produced by SNe and with grain growth in the ISM contributing substantially to the dust mass accumulation.

In hierarchical models of structure formation, the first galaxies form in low-mass dark matter potential wells, probing the behavior of dark matter on kiloparsec (kpc) scales. Even though these objects are below the detection threshold of current telescopes, future missions will open an observational window into this emergent world. In this Letter we investigate how the first galaxies are assembled in a `fuzzy' dark matter (FDM) cosmology where dark matter is an ultralight $\sim 10^{-22}$~eV boson and the primordial stars are expected to form along dense dark matter filaments. Using a first-of-its-kind cosmological hydrodynamical simulation, we explore the interplay between baryonic physics and unique wavelike features inherent to FDM. In our simulation, the dark matter filaments show coherent interference patterns on the boson de Broglie scale and develop cylindrical soliton-like cores which are unstable under gravity and collapse into kpc-scale spherical solitons. Features of the dark matter distribution are largely unaffected by the baryonic feedback. On the contrary, the distributions of gas and stars, which do form along the entire filament, exhibit central cores imprinted by dark matter -- a smoking gun signature of FDM.

Milky Way (MW) satellites reside within dark matter (DM) subhalos with a broad distribution of circular velocity profiles. This diversity is enhanced with the inclusion of ultra-faint satellites, which seemingly have very high DM densities, albeit with large systematic uncertainties. We argue that if confirmed, this large diversity in the MW satellite population poses a serious test for the structure formation theory with possible implications for the DM nature. For the Cold Dark Matter model, the diversity might be a signature of the combined effects of subhalo tidal disruption by the MW disk and strong supernova feedback. For models with a dwarf-scale cutoff in the power spectrum, the diversity is a consequence of the lower abundance of dwarf-scale halos. This diversity is most challenging for Self-Interacting Dark Matter (SIDM) models with cross sections $\sigma/m_\chi\gtrsim1~$cm$^2$g$^{-1}$ where subhalos have too low densities to explain the ultra-faint galaxies. We propose a novel solution to explain the diversity of MW satellites based on the gravothermal collapse of SIDM haloes. This solution requires a velocity-dependent cross section that predicts a bimodal distribution of cuspy dense (collapsed) subhaloes consistent with the ultra-faint satellites, and cored lower density subhaloes consistent with the brighter satellites.

The development of methods and algorithms to solve the $N$-body problem for classical, collisionless, non-relativistic particles has made it possible to follow the growth and evolution of cosmic dark matter structures over most of the Universe's history. In the best studied case $-$ the cold dark matter or CDM model $-$ the dark matter is assumed to consist of elementary particles that had negligible thermal velocities at early times. Progress over the past three decades has led to a nearly complete description of the assembly, structure and spatial distribution of dark matter haloes, and their substructure in this model, over almost the entire mass range of astronomical objects. On scales of galaxies and above, predictions from this standard CDM model have been shown to provide a remarkably good match to a wide variety of astronomical data over a large range of epochs, from the temperature structure of the cosmic background radiation to the large-scale distribution of galaxies. The frontier in this field has shifted to the relatively unexplored subgalactic scales, the domain of the central regions of massive haloes, and that of low-mass haloes and subhaloes, where potentially fundamental questions remain. Answering them may require: (i) the effect of known but uncertain baryonic processes (involving gas and stars), and/or (ii) alternative models with new dark matter physics. Here we present a review of the field, focusing on our current understanding of dark matter structure from $N$-body simulations and on the challenges ahead.

We perform a series of cosmological hydrodynamic simulations to investigate the effects of non-gravitational dark matter (DM) interactions on the intergalactic medium (IGM). In particular, we use the ETHOS framework (Cyr-Racine et al. 2016; Vogelsberger et al. 2016) to compare statistics of the Lyman-$\alpha$ forest in cold dark matter (CDM) with an alternative model in which the DM couples strongly with a relativistic species in the early universe. These models are characterised by a cutoff in the linear power spectrum, followed by a series of 'dark acoustic oscillations' (DAOs) on sub-dwarf scales. While the primordial cutoff delays the formation of the first galaxies, structure builds-up more rapidly in the interacting DM model compared to CDM. We show that although DAOs are quickly washed away in the non-linear clustering of DM at $z\lesssim10$, their signature can be imprinted prominently in the Lyman-$\alpha$ flux power spectrum at $z>5$. On scales larger than the cutoff ($k\sim0.08$ s/km for the specific model considered here), the relative difference to CDM is reminiscent of a warm dark matter (WDM) model with a similar initial cutoff; however, the redshift evolution on smaller scales is distinctly different. The appearance and disappearance of DAOs in the Lyman-$\alpha$ flux spectrum provides a powerful way to distinguish interacting DM models from WDM and, indeed, variations in the thermal history of the IGM that may also induce a small-scale cutoff.

We present deep continuum observations at a wavelength of 2mm centered on the COSMOS field using the Goddard IRAM Superconducting Millimeter Observer (GISMO) at the IRAM 30m-telescope. These data constitute the widest deep 2mm survey to-date, reaching a uniform $\sigma\sim0.23$mJy beam$^{-1}$ sensitivity over $\sim250$arcmin$^2$ at $\sim\,24''$ resolution. We detect four sources at high significance (S/N$\,\geq\,$4.4) with an expected number of false detection of 0.09 sources, and five sources at $4.4>\,$S/N$\,\geq\,3.7$ with an expected number of false detection of 1.65 sources. Combined with deep GISMO observations in GOODS-N, we constrain the 2mm number counts over one decade in flux density. These measurements agree with most galaxy evolution models tested here, except those with large population of dusty star-forming galaxies at $z>7$. Five GISMO sources have counterparts in (sub-)millimeter catalogs available in COSMOS. Their redshifts suggest that all but one lie above $z\sim3$. These four high-redshift ($z>3$) galaxies have $\tilde{z}=3.9$, SFRs $\sim\,$400 - 1200M$_{\odot}\,$yr$^{-1}$ and $M_{\rm dust}\sim10^{9.5}\,$M$_{\odot}$. They provide a relatively complete selection ($\sim66\%$) of the most luminous ($L_{\rm IR}>10^{12.6}\,$L$_{\odot}$) and highest redshift ($z>3$) galaxies detected within our survey area by AzTEC at 1.1mm. We thus conclude that 2mm surveys favor the selection of massive, vigorously star-forming, high-redshift galaxies. This is corroborated by GISMO-C4, a source with a low false detection probability ($\sim\,6.2\%$), for which the absence of a (sub-)millimeter counterpart supports a high redshift origin ($z\gtrsim3$).

Self-interacting dark matter provides a promising alternative for the cold dark matter paradigm to solve potential small-scale galaxy formation problems. Nearly all self-interacting dark matter simulations so far have considered only elastic collisions. Here we present simulations of a galactic halo within a generic inelastic model using a novel numerical implementation in the Arepo code to study arbitrary multi-state inelastic dark matter scenarios. For this model we find that inelastic self-interactions can: (i) create larger subhalo density cores compared to elastic models for the same cross section normalisation; (ii) lower the abundance of satellites without the need for a power spectrum cutoff; (iii) reduce the total halo mass by about 10%; (iv) inject the energy equivalent of O(100) million Type II supernovae in galactic haloes through level de-excitation; (v) avoid the gravothermal catastrophe due to removal of particles from halo centers. We conclude that a ~5 times larger elastic cross section is required to achieve the same central density reduction as the inelastic model. This implies that well-established constraints on self-interacting cross sections have to be revised if inelastic collisions are the dominant mode. In this case significantly smaller cross sections can achieve the same core density reduction thereby increasing the parameter space of allowed models considerably.

We present the science case for mapping several thousand galaxy (proto)clusters at z=1-10 with a large aperture single dish sub-mm facility, producing a high-redshift counterpart to local large surveys of rich clusters like the well-studied Abell catalogue. Principal goals of a large survey of distant clusters are the evolution of galaxy clusters over cosmic time and the impact of environment on the evolution and formation of galaxies. To make a big leap forward in this emerging research field, the community would benefit from a large-format, wide-band, direct-detection spectrometer (e.g., based on MKID technology), covering a wide field of ~1 square degree and a frequency coverage from 70 to 700 GHz.

Two decades of effort have been poured into both single-dish and interferometric millimeter-wave surveys of the sky to infer the volume density of dusty star-forming galaxies (DSFGs, with SFR>100M$_\odot$ yr$^{-1}$) over cosmic time. Though obscured galaxies dominate cosmic star-formation near its peak at $z\sim2$, the contribution of such heavily obscured galaxies to cosmic star-formation is unknown beyond $z\sim2.5$ in contrast to the well-studied population of Lyman-break galaxies (LBGs) studied through deep, space- and ground-based pencil beam surveys in the near-infrared. Unlocking the volume density of DSFGs beyond $z>3$, particularly within the first 1 Gyr after the Big Bang is critical to resolving key open questions about early Universe galaxy formation: (1) What is the integrated star-formation rate density of the Universe in the first few Gyr and how is it distributed among low-mass galaxies (e.g. Lyman-break galaxies) and high-mass galaxies (e.g. DSFGs and quasar host galaxies)? (2) How and where do the first massive galaxies assemble? (3) What can the most extreme DSFGs teach us about the mechanisms of dust production (e.g. supernovae, AGB stars, grain growth in the ISM) <1 Gyr after the Big Bang? We summarize the types of observations needed in the next decade to address these questions.

It is well established that the central deficit of dark matter (DM) observed in many dwarf galaxies disagrees with the cuspy DM haloes predicted in the collision-less and cold DM (CDM) model. Plausible solutions to this problem are based on an effective energy deposition into the central halo with an origin that is either based on baryonic physics (e.g. supernova-driven gas blowouts; SNF) or on new DM physics (e.g. self-interacting DM, SIDM). We argue that the fundamental difference between the two is whether the process is impulsive or adiabatic, and explore novel ways to distinguishing them by looking at the response of stellar orbits. We perform idealised simulations of tracers embedded in a $1.48\times 10^{10}{\rm M}_{\odot}$ spherical halo, and model the creation of a $\sim1$kpc DM matter core in SIDM with $\sigma_T/m_{\chi} = 2\,{\rm cm}^2{\rm g}^{-1}$ and through SNF by impulsively removing a mass equivalent to $\mathcal{O}(10\%)$ of the halo's potential energy. Choosing idealised initial orbital configurations for the tracers, we find that radial actions are conserved (changed) in the SIDM (impulsive) case. The adiabaticity of the SIDM case prevents tracers from changing their orbital family during core formation whereas SNF separates tracers of initially the same family to a variety of orbits. We show that these key features remain in a cosmological halo, albeit for a few dynamical timescales. The number density and velocity dispersion profile of a Plummer sphere with $r_{1/2} = 500$ pc change only marginally under adiabatic core formation, whereas SNF causes a substantial expansion of the sphere driving it out of Jeans equilibrium. Our results point towards promising ways of differentiating adiabatic from impulsive core formation.

In this work, we discuss the relation between the strength of the self-interaction of dark matter particles and the predicted properties of the inner density distributions of dark matter haloes. We present the results of $N$-body simulations for 28 haloes performed with the same initial conditions for cold dark matter and for self-interacting dark matter, for a range of cross-sections. We provide a simple phenomenological description of these results and compare with the semi-analytical model typically used in the literature. Using these results, we then predict how the inner dark matter surface density and core radius should depend on the self-interaction cross-section for observed haloes.

The "sloshing" of the cold gas in the cores of relaxed clusters of galaxies is a widespread phenomenon, evidenced by the presence of spiral-shaped "cold fronts" in X-ray observations of these systems. In simulations, these flows of cold gas readily form by interactions of the cluster core with small subclusters, due to a separation of the cold gas from the dark matter (DM), due to their markedly different collisionalities. In this work, we use numerical simulations to investigate the effects of increasing the DM collisionality on sloshing cold fronts in a cool-core cluster. For clusters in isolation, the formation of a flat DM core via self-interactions results in modest adiabatic expansion and cooling of the core gas. In merger simulations, cold fronts form in the same manner as in previous simulations, but the flattened potential in the core region enables the gas to expand to larger radii in the initial stages. Upon infall, the subcluster's DM mass decreases via collisions, reducing its influence on the core. Thus, the sloshing gas moves slower, inhibiting the growth of fluid instabilities relative to simulations where the DM cross section is zero. This also inhibits turbulent mixing and the increase in entropy that would otherwise result. For values of the cross section $\sigma/m > 1$, subclusters do not survive as self-gravitating structures for more than two core passages. Additionally, separations between the peaks in the X-ray emissivity and thermal Sunyaev-Zeldovich effect signals during sloshing may place constraints on DM self-interactions.

We use high-resolution hydrodynamical simulation to test the difference of halo properties in cold dark matter (CDM) and a self-interacting dark matter (SIDM) scenario with a constant cross-section of $\sigma^\text{T}/m_{\chi}=1\;\text{cm}^{2}\text{g}^{-1}$. We find that the interplay between dark matter self-interaction and baryonic physics induces a complex evolution of the halo properties, which depends on the halo mass and morphological type, as well as on the halo mass accretion history. While high mass haloes, selected as analogues of early-type galaxies, show cored profiles in the SIDM run, systems of intermediate mass and with a significant disk component can develop a profile that is similar or cuspier than in CDM. The final properties of SIDM haloes - measured at z=0.2 - correlate with the halo concentration and formation time, suggesting that the differences between different systems are due to the fact that we are observing the impact self-interaction. We also search for signatures of self-interacting dark matter in the lensing signal of the main haloes and hints of potential differences in the distribution of Einstein radii, which suggests that future wide-field survey might be able to distinguish between CDM and SIDM models on this basis. Finally, we find that the subhalo abundances are not altered in the adopted SIDM model with respect to CDM.

We contrast predictions for the high-redshift galaxy population and reionization history between cold dark matter (CDM) and an alternative self-interacting dark matter model based on the recently developed ETHOS framework that alleviates the small-scale CDM challenges within the Local Group. We perform the highest resolution hydrodynamical cosmological simulations (a 36~Mpc$^3$ volume with gas cell mass of $\sim10^5\mathrm{M}_{\odot}$ and minimum gas softening of $\sim180$~pc) within ETHOS to date -- plus a CDM counterpart -- to quantify the abundance of galaxies at high redshift and their impact on reionization. We find that ETHOS predicts galaxies with higher ultraviolet (UV) luminosities than their CDM counterparts and a faster build-up of the faint end of the UV luminosity function. These effects, however, make the optical depth to reionization less sensitive to the power spectrum cut-off: the ETHOS model differs from the CDM $\tau$ value by only 10 per cent and is consistent with Planck limits if the effective escape fraction of UV photons is 0.1-0.5. We conclude that current observations of high-redshift luminosity functions cannot differentiate between ETHOS and CDM models, but deep JWST surveys of strongly-lensed, inherently faint galaxies have the potential to test non-CDM models that offer attractive solutions to CDM's Local Group problems.

A cutoff in the linear matter power spectrum at dwarf galaxy scales has been shown to affect the abundance, formation mechanism and age of dwarf haloes and their galaxies at high and low redshift. We use hydrodynamical simulations of galaxy formation within the ETHOS framework in a benchmark model that has such a cutoff, and that has been shown to be an alternative to the cold dark matter (CDM) model that alleviates its dwarf-scale challenges. We show how galaxies in this model form differently to CDM on a halo-by-halo basis, at redshifts $z\ge6$. We show that ETHOS haloes at the half-mode mass scale form with 50~per~cent less mass than their CDM counterparts due to their later formation times, yet they retain more of their gas reservoir due to the different behaviour of gas and dark matter during the monolithic collapse of the first haloes in models with a galactic-scale cutoff. As a result, galaxies in ETHOS haloes near the cutoff scale grow rapidly between $z=10-6$ and by $z=6$ end up having very similar stellar masses, higher gas fractions and higher star formation rates relative to their CDM counterparts. We highlight these differences by making predictions for how the number of galaxies with old stellar populations is suppressed in ETHOS for both $z=6$ galaxies and for gas-poor Local Group fossil galaxies. Interestingly, we find an age gradient in ETHOS between galaxies that form in high and low density environments.

Observations of the baryon to dark matter fraction in galaxies through cosmic time are a fundamental test for galaxy formation models. Recent observational studies have suggested that some disk galaxies at $z>1$ host declining rotation curves, in contrast with observations of low redshift disk galaxies where stellar or HI rotation curves flatten at large radii. We present an observational counterexample, a galaxy named DSFG850.95 at $z=1.555$ (4.1 Gyr after the big bang) that hosts a flat rotation curve between radii of $\sim$6--14 kpc (1.2--2.8 disk scale lengths) and has a dark matter fraction of $0.44\pm0.08$ at the H-band half light radius, similar to the Milky Way. We create position-velocity and position-dispersion diagrams using Keck/MOSFIRE spectroscopic observations of H$\alpha$ and [NII] emission features, which reveal a flat rotation velocity of $V_{\rm flat}=285\pm12$ km/s and an ionized gas velocity dispersion of $\sigma_{0}=48\pm4$ km/s. This galaxy has a rotation-dominated velocity field with $V_{\rm flat}/\sigma_{0}\sim6$. Ground-based H-band imaging reveals a disk with S\'ersic index of $1.29\pm0.03$, an edge-on inclination angle of $87\pm2^{\circ}$, and an H-band half light radius of $8.4\pm0.1$ kpc. Our results point to DSFG850.95 being a massive, rotationally-supported disk galaxy with a high dark-matter-to-baryon fraction in the outer galaxy, similar to disk galaxies at low redshift.

We carry out a blind search of 3mm continuum sources using the ALMA Science Archive to derive the first galaxy number counts at this wavelength. The analyzed data are drawn from observations towards three extragalactic legacy fields: COSMOS, CDF-S, and the UDS comprising more than 130 individual ALMA Band 3 pointings and an effective survey area of ~200 sq. arcmin with a continuum sensitivity that allows for the direct detection of unlensed Dusty Star-Forming Galaxies (DSFGs) dust emission beyond the epoch of reionization. We present a catalog of 16 sources detected at >5sigma with flux densities S_3mm = 60 - 600 uJy from which number counts are derived. These number counts are then used to place constraints on the volume density of DSFGs with an empirical backward evolution model. Our measured 3mm number counts indicate that the contribution of DSFGs to the cosmic star formation rate density at z >~ 4 is non-negligible. This is contrary to the generally adopted assumption of a sharply decreasing contribution of obscured galaxies at z > 4 as inferred by optical and near-infrared surveys. This work demonstrates the power of ALMA 3mm observations which can reach outstanding continuum sensitivities during typical spectral line science programs. Further constraints on 3mm-selected galaxies will be essential to refine models of galaxy formation and evolution as well as models of early Universe dust production mechanisms.

Strong gravitational lensing has been identified as a promising astrophysical probe to study the particle nature of dark matter. In this paper we present a detailed study of the power spectrum of the projected mass density (convergence) field of substructure in a Milky Way-sized halo. This power spectrum has been suggested as a key observable that can be extracted from strongly lensed images and yield important clues about the matter distribution within the lens galaxy. We use two different $N$-body simulations from the ETHOS framework: one with cold dark matter and another with self-interacting dark matter and a cutoff in the initial power spectrum. Despite earlier works that identified $ k \gtrsim 100$ kpc$^{-1}$ as the most promising scales to learn about the particle nature of dark matter we find that even at lower wavenumbers - which are actually within reach of observations in the near future - we can gain important information about dark matter. Comparing the amplitude and slope of the power spectrum on scales $0.1 \lesssim k/$kpc$^{-1} \lesssim 10$ from lenses at different redshifts can help us distinguish between cold dark matter and other exotic dark matter scenarios that alter the abundance and central densities of subhalos. Furthermore, by considering the contribution of different mass bins to the power spectrum we find that subhalos in the mass range $10^7 - 10^8$ M$_{\odot}$ are on average the largest contributors to the power spectrum signal on scales $2 \lesssim k/$kpc$^{-1} \lesssim 15$, despite the numerous subhalos with masses $> 10^8$ M$_{\odot}$ in a typical lens galaxy. Finally, by comparing the power spectra obtained from the subhalo catalogs to those from the particle data in the simulation snapshots we find that the seemingly-too-simple halo model is in fact a fairly good approximation to the much more complex array of substructure in the lens.

We devise an optimal method to measure the temporal power spectrum of the lensing and intrinsic fluctuations of multiply-imaged strongly lensed quasar light curves, along with the associated time delays. The method is based on a Monte-Carlo Markov Chain (MCMC) sampling of a putative gaussian likelihood, and accurately recovers the input properties of simulated light curves, as well as the "Time Delay Challenge". We apply this method to constrain the dimensionless cosmological (non-linear) matter power spectrum on milliparsec scales (comparable to the size of the solar system), to $\Delta_{\rm NL}^2< 4 \times 10^7$ at $k_{\rm NL} \sim 10^3 {\rm pc}^{-1}$. Using a semi-analytic nonlinear clustering model which is calibrated to simulations, the corresponding constraint on the primordial (linear) scalar power spectrum is ${\cal P}_{\cal R} < 3 \times 10^{-9}$ at $k_{\rm L}\sim$ 3 pc$^{-1}$. This is the strongest constraint on primordial power spectrum at these scales, and is within an order of magnitude from the standard $\Lambda$CDM prediction. We also report measurements of temporal spectra for intrinsic variabilities of quasar light curves, which can be used to constrain the size of the emitting region in accretion disks. Future cadenced optical imaging surveys, such as LSST, should increase the number of observed strongly lensed quasars by 3 orders of magnitude and significantly improve these measurements, even though improvements in modelling quasar accretion and stellar microlensing are necessary.

Deep, pencil-beam surveys from ALMA at 1.1-1.3mm have uncovered an apparent absence of high-redshift dusty galaxies, with existing redshift distributions peaking around $z\sim1.5-2.5$. This has led to a perceived dearth of dusty systems at $z>4$, and the conclusion, according to some models, that the early Universe was relatively dust-poor. In this paper, we extend the backward evolution galaxy model described by Casey et al. (2018) to the ALMA regime (in depth and area) and determine that the measured number counts and redshift distributions from ALMA deep field surveys are fully consistent with constraints of the infrared luminosity function (IRLF) at $z<2.5$ determined by single-dish submillimeter and millimeter surveys conducted on much larger angular scales ($\sim1-10$deg$^{2}$). We find that measured 1.1-1.3mm number counts are most constraining for the measurement of the faint-end slope of the IRLF at $z<2.5$ instead of the prevalence of dusty galaxies at $z>4$. Recent studies have suggested that UV-selected galaxies at $z>4$ may be particularly dust-poor, but we find their millimeter-wave emission cannot rule out consistency with the Calzetti dust attenuation law even by assuming relatively typical, cold-dust ($T_{\rm dust}\approx30\,$K) SEDs. Our models suggest that the design of ALMA deep fields requires substantial revision to constrain the prevalence of $z>4$ early Universe obscured starbursts. The most promising avenue for detection and characterization of such early dusty galaxies will come from future ALMA 2mm blank field surveys covering a few hundred arcmin$^{2}$ and the combination of existing and future dual-purpose 3mm datasets.

We employ isolated N-body simulations to study the response of self-interacting dark matter (SIDM) halos in the presence of the baryonic potentials. Dark matter self-interactions lead to kinematic thermalization in the inner halo, resulting in a tight correlation between the dark matter and baryon distributions. A deep baryonic potential shortens the phase of SIDM core expansion and triggers core contraction. This effect can be further enhanced by a large self-scattering cross section. We find the final SIDM density profile is sensitive to the baryonic concentration and the strength of dark matter self-interactions. Assuming a spherical initial halo, we also study evolution of the SIDM halo shape together with the density profile. The halo shape at later epochs deviates from spherical symmetry due to the influence of the non-spherical disc potential, and its significance depends on the baryonic contribution to the total gravitational potential, relative to the dark matter one. In addition, we construct a multi-component model for the Milky Way, including an SIDM halo, a stellar disc and a bulge, and show it is consistent with observations from stellar kinematics and streams.

Though half of cosmic starlight is absorbed by dust and reradiated at long wavelengths (3$\mu$m-3mm), constraints on the infrared through millimeter galaxy luminosity function (the `IRLF') are poor in comparison to the rest-frame ultraviolet and optical galaxy luminosity function, particularly at z>2.5. Here we present a backward evolution model for interpreting number counts, redshift distributions, and cross-band flux density correlations in the infrared and submillimeter sky, from 70$\mu$m-2mm, using a model for the IRLF out to the epoch of reionization. Mock submillimeter maps are generated by injecting sources according to the prescribed IRLF and flux densities drawn from model spectral energy distributions that mirror the distribution of SEDs observed in $0<z<5$ dusty star-forming galaxies (DSFGs). We explore two extreme hypothetical case-studies: a dust-poor early Universe model, where DSFGs contribute negligibly ($<$10%) to the integrated star-formation rate density at $z>4$, and an alternate dust-rich early Universe model, where DSFGs dominate $\sim$90% of $z>4$ star-formation. We find that current submm/mm datasets do not clearly rule out either of these extreme models. We suggest that future surveys at 2mm will be crucial to measuring the IRLF beyond $z\sim4$. The model framework developed in this paper serves as a unique tool for the interpretation of multiwavelength IR/submm extragalactic datasets and will enable more refined constraints on the IRLF than can be made from direct measurements of individual galaxies' integrated dust emission.

Since their discovery, submillimetre-selected galaxies (SMGs) have revolutionized the field of galaxy formation and evolution. From the hundreds of square degrees mapped at submillimetre wavelengths, only a handful of sources have been confirmed to lie at z>5 and only two at z>6. All of these SMGs are rare examples of extreme starburst galaxies with star formation rates (SFRs) of >1000 M_sun/yr and therefore are not representative of the general population of dusty star-forming galaxies. Consequently, our understanding of the nature of these sources, at the earliest epochs, is still incomplete. Here we report the spectroscopic identification of a gravitationally amplified (mu = 9.3 +/- 1.0) dusty star-forming galaxy at z=6.027. After correcting for gravitational lensing, we derive an intrinsic less-extreme SFR of 380 +/- 50 M_sun/yr for this source and find that its gas and dust properties are similar to those measured for local Ultra Luminous Infrared Galaxies (ULIRGs), extending the local trends to a poorly explored territory in the early Universe. The star-formation efficiency of this galaxy is similar to those measured in its local analogues, despite a ~12 Gyr difference in cosmic time.

We present a multi-wavelength analysis of galaxies selected at 450 and 850um from the deepest SCUBA-2 observations in the Extended Groth Strip (EGS) field, which have an average depth of sigma_450=1.9 and sigma_850=0.46 mJy/beam over ~70 sq. arcmin. The final sample comprises 95 sources: 56 (59%) are detected at both wavelengths, 31 (33%) are detected only at 850um, and 8 (8%) are detected only at 450um. We identify counterparts for 75% of the whole sample. The redshift distributions of the 450 and 850um samples peak at different redshifts with median values of z=1.66 +/- 0.18 and z=2.30 +/- 0.20, respectively. However, the two populations have similar IR luminosities, SFRs, and stellar masses, with mean values of (1.5 +/- 0.2) x 10^12 L_sun, (150 +/- 20) M_sun/yr, and (9.0 +/- 0.6) x 10^10 M_sun, respectively. This places most of our sources (>85%) on the high-mass end of the `main-sequence' of star-forming galaxies. Exploring the IR excess vs UV-slope (IRX-beta) relation we find that the most luminous galaxies are consistent with the Meurer law, while the less luminous galaxies lie below this relation. Using the results of a two-dimensional modelling of the HST H_160-band imaging, we derive a median Sersic index of n=1.4 +0.3 -0.1 and a median half-light radius of R_1/2 = 4.8 +/ 0.4 kpc. Based on a visual-like classification in the same band, we find that the dominant component for most of the galaxies at all redshifts is a disk-like structure, although there is a transition from irregular disks to disks with a spheroidal component at z~1.4, which morphologically supports the scenario of SMGs as progenitors of massive elliptical galaxies.

As a binary pulsar moves through a wind of dark matter particles, the resulting dynamical friction modifies the binary's orbit. We study this effect for the double disk dark matter (DDDM) scenario, where a fraction of the dark matter is dissipative and settles into a thin disk. For binaries within the dark disk, this effect is enhanced due to the higher dark matter density and lower velocity dispersion of the dark disk, and due to its co-rotation with the baryonic disk.We estimate the effect and compare it with observations for two different limits in the Knudsen number ($Kn$). First, in the case where DDDM is effectively collisionless within the characteristic scale of the binary ($Kn\gg1$) and ignoring the possible interaction between the pair of dark matter wakes. Second, in the fully collisional case ($Kn\ll1$), where a fluid description can be adopted and the interaction of the pair of wakes is taken into account. We find that the change in the orbital period is of the same order of magnitude in both limits. A comparison with observations reveals good prospects to probe currently allowed DDDM models with timing data from binary pulsars in the near future. We finally comment on the possibility of extending the analysis to the intermediate (rarefied gas) case with $Kn\sim1$.

We perform dark-matter-only simulations of 28 relaxed massive cluster-sized haloes for Cold Dark Matter (CDM) and Self-Interacting Dark Matter (SIDM) models, to study structural differences between the models at large radii, where the impact of baryonic physics is expected to be very limited. We find that the distributions for the radial profiles of the density, ellipsoidal axis ratios, and velocity anisotropies ($\beta$) of the haloes differ considerably between the models (at the $\sim1\sigma$ level), even at $\gtrsim10\%$ of the virial radius, if the self-scattering cross section is $\sigma/m_\chi=1$ cm$^2$ gr$^{-1}$. Direct comparison with observationally inferred density profiles disfavours SIDM for $\sigma/m_\chi=1$ cm$^2$ gr$^{-1}$, but in an intermediate radial range ($\sim3\%$ of the virial radius), where the impact of baryonic physics is uncertain. At this level of the cross section, we find a narrower $\beta$ distribution in SIDM, clearly skewed towards isotropic orbits, with no SIDM (90\% of CDM) haloes having $\beta>0.12$ at $7\%$ of the virial radius. We estimate that with an observational sample of $\sim30$ ($\sim10^{15}$ M$_\odot$) relaxed clusters, $\beta$ can potentially be used to put competitive constraints on SIDM, once observational uncertainties improve by a factor of a few. We study the suppression of the memory of halo assembly history in SIDM clusters. For $\sigma/m_\chi=1$ cm$^2$ gr$^{-1}$, we find that this happens only in the central halo regions ($\sim1/4$ of the scale radius of the halo), and only for haloes that assembled their mass within this region earlier than a formation redshift $z_f\sim2$. Otherwise, the memory of assembly remains and is reflected in ways similar to CDM, albeit with weaker trends.

We formulate an effective theory of structure formation (ETHOS) that enables cosmological structure formation to be computed in almost any microphysical model of dark matter physics. This framework maps the detailed microphysical theories of particle dark matter interactions into the physical effective parameters that shape the linear matter power spectrum and the self-interaction transfer cross section of non-relativistic dark matter. These are the input to structure formation simulations, which follow the evolution of the cosmological and galactic dark matter distributions. Models with similar effective parameters in ETHOS but with different dark particle physics would nevertheless result in similar dark matter distributions. We present a general method to map an ultraviolet complete or effective field theory of low energy dark matter physics into parameters that affect the linear matter power spectrum and carry out this mapping for several representative particle models. We further propose a simple but useful choice for characterizing the dark matter self-interaction transfer cross section that parametrizes self-scattering in structure formation simulations. Taken together, these effective parameters in ETHOS allow the classification of dark matter theories according to their structure formation properties rather than their intrinsic particle properties, paving the way for future simulations to span the space of viable dark matter physics relevant for structure formation.

We present a theoretical analysis of some unexplored aspects of relaxed Bose-Einstein condensate dark matter (BECDM) haloes. This type of ultralight bosonic scalar field dark matter is a viable alternative to the standard cold dark matter (CDM) paradigm, as it makes the same large-scale predictions as CDM and potentially overcomes CDM's small-scale problems via a galaxy-scale de Broglie wavelength. We simulate BECDM halo formation through mergers, evolved under the Schr\"odinger-Poisson equations. The formed haloes consist of a soliton core supported against gravitational collapse by the quantum pressure tensor and an asymptotic $r^{-3}$ NFW-like profile. We find a fundamental relation of the core=to-halo mass with the dimensionless invariant $\Xi \equiv \lvert E \rvert/M^3/(Gm/\hbar)^2$ or $M_{\rm c}/M \simeq 2.6 \Xi^{1/3}$, linking the soliton to global halo properties. For $r \geq 3.5 \,r_{\rm c}$ core radii, we find equipartition between potential, classical kinetic, and quantum gradient energies. The haloes also exhibit a conspicuous turbulent behavior driven by the continuous reconnection of vortex lines due to wave interference. We analyse the turbulence 1D velocity power spectrum and find a $k^{-1.1}$ power-law. This suggests the vorticity in BECDM haloes is homogeneous, similar to thermally-driven counterflow BEC systems from condensed matter physics, in contrast to a $k^{-5/3}$ Kolmogorov power-law seen in mechanically-driven quantum systems. The mode where the power spectrum peaks is approximately the soliton width, implying the soliton-sized granules carry most of the turbulent energy in BECDM haloes.

We explore for the first time the effect of self-interacting dark matter (SIDM) on the dark matter (DM) and baryonic distribution in massive galaxies formed in hydrodynamical cosmological simulations, including explicit baryonic physics treatment. A novel implementation of Super-Massive Black Hole (SMBH) formation and evolution is used, as in Tremmel et al.(2015, 2016), allowing to explicitly follow SMBH dynamics at the center of galaxies. A high SIDM constant cross-section is chosen, $\sigma$=10 $\rm cm^2/gr$, to amplify differences from CDM models. Milky Way-like galaxies form a shallower DM density profile in SIDM than they do in CDM, with differences already at 20 kpc scales. This demonstrates that even for the most massive spirals the effect of SIDM dominates over the adiabatic contraction due to baryons. Strikingly, the dynamics of SMBHs differs in the SIDM and reference CDM case. SMBHs in massive spirals have sunk to the centre of their host galaxy in both the SIDM and CDM run, while in less massive galaxies about 80$\%$ of the SMBH population is off-centered in the SIDM case, as opposed to the CDM case in which $\sim$90$\%$ of SMBHs have reached their host's centre. SMBHs are found as far as $\sim$9 kpc away from the centre of their host SIDM galaxy. This difference is due to the increased dynamical friction timescale caused by the lower DM density in SIDM galaxies compared to CDM, resulting in 'core stalling'. This pilot work highlights the importance of simulating in a full hydrodynamical context different DM models combined to SMBH physics to study their influence on galaxy formation.

Galactic rotation curves are a fundamental constraint for any cosmological model. We use controlled N-body simulations of galaxies to study the gravitational effect of baryons in a scenario with collisionless cold dark matter (CDM) versus one with a self-interacting dark matter (SIDM) component. In particular, we examine the inner profiles of the rotation curves in the velocity range Vmax = 30-250 km/s, whose diversity has been found to be greater than predicted by the Lambda-CDM scenario. We find that the scatter in the observed rotation curves exceeds that predicted by dark matter only mass-concentration relations in either the CDM nor SIDM models. Allowing for realistic baryonic content and spatial distributions, however, helps create a large variety of rotation curve shapes; which is in better agreement with observations in the case of self-interactions due to the characteristic cored profiles being more accommodating to the slowly rising rotation curves than CDM. We find individual fits to model two of the most remarkable outliers of similar Vmax, UGC 5721 and IC 2574, the former a cusp-like rotation curve and the latter a seemingly 8 kpc cored profile. This diversity in SIDM arises as permutations of overly concentrated halos with compact baryonic distributions versus underdense halos with extended baryonic disks. The SIDM solution is promising and its feasibility ultimately depends on the sampling of the halo mass-concentration relation and its interplay with the baryonic profiles, emphasising the need for a better understanding of the frequency of extreme outliers present in current observational samples.

We directly detect dust emission in an optically-detected, multiply-imaged galaxy lensed by the Frontier Fields cluster MACSJ0717.5+3745. We detect two images of the same galaxy at 1.1mm with the AzTEC camera on the Large Millimeter Telescope leaving no ambiguity in the counterpart identification. This galaxy, MACS071_Az9, is at z>4 and the strong lensing model (mu=7.5) allows us to calculate an intrinsic IR luminosity of 9.7e10 Lsun and an obscured star formation rate of 14.6 +/- 4.5 Msun/yr. The unobscured star formation rate from the UV is only 4.1 +/- 0.3 Msun/yr which means the total star formation rate (18.7 +/- 4.5 Msun/yr) is dominated (75-80%) by the obscured component. With an intrinsic stellar mass of only 6.9e9Msun, MACS0717_Az9 is one of only a handful of z>4 galaxies at these lower masses that is detected in dust emission. This galaxy lies close to the estimated star formation sequence at this epoch. However, it does not lie on the dust obscuration relation (IRX-beta) for local starburst galaxies and is instead consistent with the Small Magellanic Cloud (SMC) attenuation law. This remarkable lower mass galaxy showing signs of both low metallicity and high dust content may challenge our picture of dust production in the early Universe.

The gravitationally-lensed galaxy A1689-zD1 is one of the most distant spectroscopically confirmed sources ($z=7.5$). It is the earliest known galaxy where the interstellar medium (ISM) has been detected; dust emission was detected with the Atacama Large Millimetre Array (ALMA). A1689-zD1 is also unusual among high-redshift dust emitters as it is a sub-L* galaxy and is therefore a good prospect for the detection of gaseous ISM in a more typical galaxy at this redshift. We observed A1689-zD1 with ALMA in bands 6 and 7 and with the Green Bank Telescope (GBT) in band $Q$. To study the structure of A1689-zD1, we map the mm thermal dust emission and find two spatial components with sizes about $0.4-1.7$\,kpc (lensing-corrected). The rough spatial morphology is similar to what is observed in the near-infrared with {\it HST} and points to a perturbed dynamical state, perhaps indicative of a major merger or a disc in early formation. The ALMA photometry is used to constrain the far-infrared spectral energy distribution, yielding a dust temperature ($T_{\rm dust} \sim 35$--$45$\,K for $\beta = 1.5-2$). We do not detect the CO(3-2) line in the GBT data with a 95\% upper limit of 0.3\,mJy observed. We find a slight excess emission in ALMA band~6 at 220.9\,GHz. If this excess is real, it is likely due to emission from the [CII] 158.8\,$\mu$m line at $z_{\rm [CII]} = 7.603$. The stringent upper limits on the [CII]/$L_{\rm FIR}$ luminosity ratio suggest a [CII] deficit similar to several bright quasars and massive starbursts.

The isotropic gamma-ray background arises from the contribution of unresolved sources, including members of confirmed source classes and proposed gamma-ray emitters such as the radiation induced by dark matter annihilation and decay. Clues about the properties of the contributing sources are imprinted in the anisotropy characteristics of the gamma-ray background. We use 81 months of Pass 7 Reprocessed data from the Fermi Large Area Telescope to perform a measurement of the anisotropy angular power spectrum of the gamma-ray background. We analyze energies between 0.5 and 500 GeV, extending the range considered in the previous measurement based on 22 months of data. We also compute, for the first time, the cross-correlation angular power spectrum between different energy bins. We find that the derived angular spectra are compatible with being Poissonian, i.e. constant in multipole. Moreover, the energy dependence of the anisotropy suggests that the signal is due to two populations of sources, contributing, respectively, below and above 2 GeV. Finally, using data from state-of-the-art numerical simulations to model the dark matter distribution, we constrain the contribution from dark matter annihilation and decay in Galactic and extragalactic structures to the measured anisotropy. These constraints are competitive with those that can be derived from the average intensity of the isotropic gamma-ray background. Data are available at https://www-glast.stanford.edu/pub_data/552.

We present deep observations at 450 um and 850 um in the Extended Groth Strip field taken with the SCUBA-2 camera mounted on the James Clerk Maxwell Telescope as part of the deep SCUBA-2 Cosmology Legacy Survey (S2CLS), achieving a central instrumental depth of $\sigma_{450}=1.2$ mJy/beam and $\sigma_{850}=0.2$ mJy/beam. We detect 57 sources at 450 um and 90 at 850 um with S/N > 3.5 over ~70 sq. arcmin. From these detections we derive the number counts at flux densities $S_{450}>4.0$ mJy and $S_{850}>0.9$ mJy, which represent the deepest number counts at these wavelengths derived using directly extracted sources from only blank-field observations with a single-dish telescope. Our measurements smoothly connect the gap between previous shallower blank-field single-dish observations and deep interferometric ALMA results. We estimate the contribution of our SCUBA-2 detected galaxies to the cosmic infrared background (CIB), as well as the contribution of 24 um-selected galaxies through a stacking technique, which add a total of $0.26\pm0.03$ and $0.07\pm0.01$ MJy/sr, at 450 um and 850 um, respectively. These surface brightnesses correspond to $60\pm20$ and $50\pm20$ per cent of the total CIB measurements, where the errors are dominated by those of the total CIB. Using the photometric redshifts of the 24 um-selected sample and the redshift distributions of the submillimetre galaxies, we find that the redshift distribution of the recovered CIB is different at each wavelength, with a peak at $z\sim1$ for 450 um and at $z\sim2$ for 850um, consistent with previous observations and theoretical models.

We investigate the effects of self-interacting dark matter (SIDM) on the tidal stripping and evaporation of satellite galaxies in a Milky Way-like host. We use a suite of five zoom-in, dark-matter-only simulations, two with velocity-independent SIDM cross sections, two with velocity-dependent SIDM cross sections, and one cold dark matter simulation for comparison. After carefully assigning stellar mass to satellites at infall, we find that stars are stripped at a higher rate in SIDM than in CDM. In contrast, the total bound dark matter mass loss rate is minimally affected, with subhalo evaporation having negligible effects on satellites for viable SIDM models. Centrally located stars in SIDM haloes disperse out to larger radii as cores grow. Consequently, the half-light radius of satellites increases, stars become more vulnerable to tidal stripping, and the stellar mass function is suppressed. We find that the ratio of core radius to tidal radius accurately predicts the relative strength of enhanced SIDM stellar stripping. Velocity-independent SIDM models show a modest increase in the stellar stripping effect with satellite mass, whereas velocity-dependent SIDM models show a large increase in this effect towards lower masses, making observations of ultra-faint dwarfs prime targets for distinguishing between and constraining SIDM models. Due to small cores in the largest satellites of velocity-dependent SIDM, no identifiable imprint is left on the all-sky properties of the stellar halo. While our results focus on SIDM, the main physical mechanism of enhanced tidal stripping of stars apply similarly to satellites with cores formed via other means.

Dark matter annihilation or decay could have a significant impact on the ionisation and thermal history of the universe. In this paper, we study the potential contribution of dark matter annihilation (s-wave- or p-wave-dominated) or decay to cosmic reionisation, via the production of electrons, positrons and photons. We map out the possible perturbations to the ionisation and thermal histories of the universe due to dark matter processes, over a broad range of velocity-averaged annihilation cross-sections/decay lifetimes and dark matter masses. We have employed recent numerical studies of the efficiency with which annihilation/decay products induce heating and ionization in the intergalactic medium, and in this work extended them down to a redshift of $1+z = 4$ for two different reionisation scenarios. We also improve on earlier studies by using the results of detailed structure formation models of dark matter haloes and subhaloes that are consistent with up-to-date $N$-body simulations, with estimates on the uncertainties that originate from the smallest scales. We find that for dark matter models that are consistent with experimental constraints, a contribution of more than 10% to the ionisation fraction at reionisation is disallowed for all annihilation scenarios. Such a contribution is possible only for decays into electron/positron pairs, for light dark matter with mass $m_\chi \lesssim $ 100 MeV, and a decay lifetime $\tau_\chi \sim 10^{24} - 10^{25}$ s.

We present a catalogue of nearly 3,000 submillimetre sources detected at 850um over ~5 square degrees surveyed as part of the James Clerk Maxwell Telescope (JCMT) SCUBA-2 Cosmology Legacy Survey (S2CLS). This is the largest survey of its kind at 850um, probing a meaningful cosmic volume at the peak of star formation activity and increasing the sample size of submillimetre galaxies selected at 850um by an order of magnitude. We describe the wide 850um survey component of S2CLS, which covers the key extragalactic survey fields: UKIDSS-UDS, COSMOS, Akari-NEP, Extended Groth Strip, Lockman Hole North, SSA22 and GOODS-North. The average 1-sigma depth of S2CLS is 1.2 mJy/beam, approaching the SCUBA-2 850um confusion limit, which we determine to be ~0.8 mJy/beam. We measure the single dish 850um number counts to unprecedented accuracy, reducing the Poisson errors on the differential counts to approximately 4% at S_850~3mJy. With several independent fields, we investigate field-to-field variance, finding that the number counts on 0.5-1 degree scales are generally within 50% of the S2CLS mean for S_850>3mJy, with scatter consistent with the Poisson and estimated cosmic variance uncertainties, although there is a marginal (2-sigma) density enhancement in the GOODS-North field. The observed number counts are in reasonable agreement with recent phenomenological and semi-analytic models. Finally, the large solid angle of S2CLS allows us to measure the bright-end counts: at S_850>10mJy there are approximately ten sources per square degree, and we detect the distinctive up-turn in the number counts indicative of the detection of local sources of 850um emission and strongly lensed high-redshift galaxies. Here we describe the data collection and reduction procedures and present calibrated maps and a catalogue of sources; these are made publicly available.

Epsilon Eridani is a nearby, young Sun-like star that hosts a ring of cool debris analogous to the solar system's Edgeworth-Kuiper belt. Early observations at (sub-)mm wavelengths gave tentative evidence of the presence of inhomogeneities in the ring, which have been ascribed to the effect of a putative low eccentricity planet, orbiting close to the ring. The existence of these structures have been recently challenged by high resolution interferometric millimeter observations. Here we present the deepest single-dish image of Epsilon Eridani at millimeter wavelengths, obtained with the Large Millimeter Telescope Alfonso Serrano (LMT). The main goal of these LMT observations is to confirm (or refute) the presence of non-axisymmetric structure in the disk. The dusty ring is detected for the first time along its full projected elliptical shape. The radial extent of the ring is not spatially resolved and shows no evidence, to within the uncertainties, of dust density enhancements. Additional features of the 1.1 mm map are: (i) the presence of significant flux in the gap between the ring and the star, probably providing the first exo-solar evidence of Poynting-Robertson drag, (ii) an unambiguous detection of emission at the stellar position with a flux significantly above that expected from Epsilon Eridani's photosphere, and (iii) the identification of numerous unresolved sources which could correspond to background dusty star-forming galaxies.

We present the first simulations within an effective theory of structure formation (ETHOS), which includes the effect of interactions between dark matter and dark radiation on the linear initial power spectrum and dark matter self-interactions during non-linear structure formation. We simulate a Milky Way-like halo in four different dark matter models and the cold dark matter case. Our highest resolution simulation has a particle mass of $2.8\times 10^4\,{\rm M}_\odot$ and a softening length of $72.4\,{\rm pc}$. We demonstrate that all alternative models have only a negligible impact on large scale structure formation. On galactic scales, however, the models significantly affect the structure and abundance of subhaloes due to the combined effects of small scale primordial damping in the power spectrum and late time self-interactions. We derive an analytic mapping from the primordial damping scale in the power spectrum to the cutoff scale in the halo mass function and the kinetic decoupling temperature. We demonstrate that certain models within this extended effective framework that can alleviate the too-big-to-fail and missing satellite problems simultaneously, and possibly the core-cusp problem. The primordial power spectrum cutoff of our models naturally creates a diversity in the circular velocity profiles, which is larger than that found for cold dark matter simulations. We show that the parameter space of models can be constrained by contrasting model predictions to astrophysical observations. For example, some models may be challenged by the missing satellite problem if baryonic processes were to be included and even over-solve the too-big-to-fail problem; thus ruling them out.

We present 8.5 arcsec resolution 1.1mm continuum imaging and CO spectroscopic redshift measurements of eight extremely bright submillimetre galaxies identified from the Planck and Herschel surveys, taken with the Large Millimeter Telescope's AzTEC and Redshift Search Receiver instruments. We compiled a candidate list of high redshift galaxies by cross-correlating the Planck Surveyor mission's highest frequency channel (857 GHz, FWHM = 4.5 arcmin) with the archival Herschel Spectral and Photometric Imaging Receiver (SPIRE) imaging data, and requiring the presence of a unique, single Herschel counterpart within the 150 arcsec search radius of the Planck source positions with 350 micron flux density larger than 100 mJy, excluding known blazars and foreground galaxies. All eight candidate objects observed are detected in 1.1mm continuum by AzTEC bolometer camera, and at least one CO line is detected in all cases with a spectroscopic redshift between 1.3 < z(CO) < 3.3. Their infrared spectral energy distributions mapped using the Herschel and AzTEC photometry are consistent with cold dust emission with characteristic temperature between $T_d$ = 43 K and 84 K. With apparent infrared luminosity of up to L(IR) = $3\times10^{14} \mu^{-1} L_\odot$, they are some of the most luminous galaxies ever found (with yet unknown gravitational magnification factor $\mu$). The analysis of their spectral energy distributions (SEDs) suggests that star formation is powering the bulk of their extremely large IR luminosities. Derived molecular gas masses of $M_{H2}=(0.6-7.8)\times 10^{11} M_\odot$ (for $\mu$~10) also make them some of the most gas-rich high redshift galaxies ever detected.

We have recently introduced a novel statistical measure of dark matter clustering in phase space, the particle phase space average density ($P^2SAD$). In a two-paper series, we studied the structure of $P^2SAD$ in the Milky-Way-size Aquarius haloes, constructed a physically motivated model to describe it, and illustrated its potential as a powerful tool to predict signals sensitive to the nanostructure of dark matter haloes. In this letter, we report a remarkable universality of the clustering of dark matter in phase space as measured by $P^2SAD$ within the subhaloes of host haloes across different environments covering a range from dwarf-size to cluster-size haloes ($10^{10}-10^{15}$ M$_\odot$). Simulations show that the universality of $P^2SAD$ holds for more than 7 orders of magnitude, over a 2D phase space, covering over 3 orders of magnitude in distance/velocity, with a simple functional form that can be described by our model. Invoking the universality of $P^2SAD$, we can accurately predict the non-linear power spectrum of dark matter at small scales all the way down to the decoupling mass limit of cold dark matter particles. As an application, we compute the subhalo boost to the annihilation of dark matter in a wide range of host halo masses.

We present multiwavelength identifications for the counterparts of 1088 submillimeter sources detected at 850$\mu$m in the SCUBA-2 Cosmology Legacy Survey study of the UKIDSS-UDS field. By utilising an ALMA pilot study on a subset of our bright SCUBA-2 sample as a training set, along with the deep optical-near-infrared data available in this field, we develop a novel technique, Optical-IR Triple Color (OIRTC), using $z-K$, $K-[3.6]$, $[3.6]-[4.5]$ colors to select the candidate submillimeter galaxy (SMG) counterparts. By combining radio identification and the OIRTC technique, we find counterpart candidates for 80% of the Class = 1 $\geq4\,\sigma$ SCUBA-2 sample, defined as those that are covered by both radio and OIR imaging and the base sample for our scientific analyses. Based on the ALMA training set, we expect the accuracy of these identifications to be $82\pm20$%, with a completeness of $69\pm16$%, essentially as accurate as the traditional $p$-value technique but with higher completeness. We find that the fraction of SCUBA-2 sources having candidate counterparts is lower for fainter 850$\mu$m sources, and we argue that for follow-up observations sensitive to SMGs with $S_{850}\gtrsim 1$ mJy across the whole ALMA beam, the fraction with multiple counterparts is likely to be $>40$% for SCUBA-2 sources at $S_{850} \gtrsim 4$ mJy. We find that the photometric redshift distribution for the SMGs is well fit by a lognormal distribution, with a median redshift of $z=2.3\pm0.1$. After accounting for the sources without any radio and/or OIRTC counterpart, we estimate the median redshift to be $z=2.6\pm0.1$ for SMGs with $S_{850} >1$ mJy. We also use this new large sample to study the clustering of SMGs and the the far-infrared properties of the unidentified submillimeter sources by stacking their Herschel SPIRE far-infrared emission.

We explore the co-evolution of the specific angular momentum of dark matter haloes and the cold baryons that comprise the galaxies within. We study over two thousand central galaxies within the reference cosmological hydrodynamical simulation of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) project. We employ a methodology within which the evolutionary history of a system is specified by the time-evolving properties of the Lagrangian particles that define it at z=0. We find a strong correlation between the evolution of the specific angular momentum of today's stars (cold gas) and that of the inner (whole) dark matter halo they are associated with. This link is particularly strong for the stars formed before the epoch of maximum expansion and subsequent collapse of the central dark matter halo (turnaround). Spheroids are typically assembled primarily from stars formed prior to turnaround, and are therefore destined to suffer a net loss of angular momentum associated with the strong merging activity during the assembly of the inner dark matter halo. Stellar discs retain their specific angular momentum since they are comprised of stars formed mainly after turnaround, from gas that mostly preserves the high specific angular momentum it acquired by tidal torques during the linear growth of the halo. Since the specific angular momentum loss of the stars is tied to the galaxy's morphology today, it may be possible to use our results to predict, statistically, the assembly history of a halo given the morphology of the galaxy it hosts.

An understanding of the mass build-up in galaxies over time necessitates tracing the evolution of cold gas (molecular and atomic) in galaxies. To that end, we have conducted a pilot study called CO Observations with the LMT of the Blind Ultra-Deep H I Environment Survey (COOL BUDHIES). We have observed 23 galaxies in and around the two clusters Abell 2192 (z = 0.188) and Abell 963 (z = 0.206), where 12 are cluster members and 11 are slightly in the foreground or background, using about 28 total hours on the Redshift Search Receiver (RSR) on the Large Millimeter Telescope (LMT) to measure the $^{12}$CO J = 1 --> 0 emission line and obtain molecular gas masses. These new observations provide a unique opportunity to probe both the molecular and atomic components of galaxies as a function of environment beyond the local Universe. For our sample of 23 galaxies, nine have reliable detections (S/N$\geq$3.6) of the $^{12}$CO line, and another six have marginal detections (2.0 < S/N < 3.6). For the remaining eight targets we can place upper limits on molecular gas masses roughly between $10^9$ and $10^{10} M_\odot$. Comparing our results to other studies of molecular gas, we find that our sample is significantly more abundant in molecular gas overall, when compared to the stellar and the atomic gas component, and our median molecular gas fraction lies about $1\sigma$ above the upper limits of proposed redshift evolution in earlier studies. We discuss possible reasons for this discrepancy, with the most likely conclusion being target selection and Eddington bias.

Recent observations with the GISMO (Goddard-IRAM Superconducting 2 Millimeter Observer) 2 mm camera revealed a detection 8 arcsec away from the lensed galaxy MACS1149-JD1 at z=9.6. Within the 17.5 arcsec FWHM GISMO beam, this detection is consistent with the position of the high-redshift galaxy and therefore, if confirmed, this object could be claimed to be the youngest galaxy producing significant quantities of dust. We present higher resolution (8.5 arcsec) observations of this system taken with the AzTEC 1.1 mm camera mounted on the Large Millimeter Telescope Alfonso Serrano. Dust continuum emission at the position of MACS1149-JD1 is not detected with an r.m.s. of 0.17 mJy/beam. However, we find a detection ~ 11 arcsec away from MACS1149-JD1, still within the GISMO beam which is consistent with an association to the GISMO source. Combining the AzTEC and GISMO photometry, together with Herschel ancillary data, we derive a z_phot= 0.7-1.6 for the dusty galaxy. We conclude therefore that the GISMO and AzTEC detections are not associated with MACS1149-JD1. From the non-detection of MACS1149-JD1 we derive the following (3 \sigma) upper limits corrected for gravitational lensing magnification and for cosmic microwave background effects: dust mass < 1.6 x 10^7 M_sun, IR luminosity < 8 x 10^10 L_sun, star formation rate < 14 M_sun/yr, and UV attenuation < 2.7 mag. These limits are comparable to those derived for other high-redshift galaxies from deep Atacama Large Millimeter/submillimeter Array (ALMA) observations.

Measuring redshifted CO line emission is an unambiguous method for obtaining an accurate redshift and total cold gas content of optically faint, dusty starburst systems. Here, we report the first successful spectroscopic redshift determination of AzTEC J095942.9+022938 ("COSMOS AzTEC-1"), the brightest 1.1mm continuum source found in the AzTEC/JCMT survey (Scott et al. 2008), through a clear detection of the redshifted CO (4-3) and CO (5-4) lines using the Redshift Search Receiver on the Large Millimeter Telescope. The CO redshift of $z=4.3420\pm0.0004$ is confirmed by the detection of the redshifted 158 micron [C II] line using the Submillimeter Array. The new redshift and Herschel photometry yield $L_{FIR}=(1.1\pm0.1)\times 10^{13} L_\odot$ and $SFR = 1300\, M_\odot$ yr$^{-1}$. Its molecular gas mass derived using the ULIRG conversion factor is $1.4\pm0.2 \times 10^{11} M_\odot$ while the total ISM mass derived from the 1.1mm dust continuum is $3.7\pm0.7 \times 10^{11} M_\odot$ assuming dust temperature of 35 K. Our dynamical mass analysis suggests that the compact gas disk ($r\approx 1.1$ kpc, inferred from dust continuum and SED analysis) has to be nearly face-on, providing a natural explanation for the uncommonly bright, compact stellar light seen by the HST. The [C II] line luminosity $L_{[C~II]} = 7.8\pm1.1 \times 10^9 L_\odot$ is remarkably high, but it is only 0.04 per cent of the total IR luminosity. AzTEC COSMOS-1 and other high redshift sources with a spatially resolved size extend the tight trend seen between [C II]/FIR ratio and $\Sigma_{FIR}$ among IR-bright galaxies reported by Diaz-Santos et al. (2013) by more than an order of magnitude, supporting the explanation that the higher intensity of the IR radiation field is responsible for the "[C II] deficiency" seen among luminous starburst galaxies.

We present Early Science observations with the Large Millimeter Telescope, AzTEC 1.1 mm continuum images and wide bandwidth spectra (73-111 GHz) acquired with the Redshift Search Receiver, towards four bright lensed submillimetre galaxies identified through the Herschel Lensing Survey-snapshot and the SCUBA-2 Cluster Snapshot Survey. This pilot project studies the star formation history and the physical properties of the molecular gas and dust content of the highest redshift galaxies identified through the benefits of gravitational magnification. We robustly detect dust continuum emission for the full sample and CO emission lines for three of the targets. We find that one source shows spectroscopic multiplicity and is a blend of three galaxies at different redshifts (z=2.040, 3.252 and 4.680), reminiscent of previous high-resolution imaging follow-up of unlensed submillimetre galaxies, but with a completely different search method, that confirm recent theoretical predictions of physically unassociated blended galaxies. Identifying the detected lines as 12CO (J_up=2-5) we derive spectroscopic redshifts, molecular gas masses, and dust masses from the continuum emission. The mean H_2 gas mass of the full sample is (2.0 +- 0.2) x 10^11 M_sun/\mu, and the mean dust mass is (2.0+-0.2) x 10^9 M_sun/\mu, where \mu=2-5 is the expected lens amplification. Using these independent estimations we infer a gas-to-dust ratio of \delta_GDR=55-75, in agreement with other measurements of submillimetre galaxies. Our magnified high-luminosity galaxies fall on the same locus as other high-redshift submillimetre galaxies, extending the L'_CO - L_FIR correlation observed for local luminous and ultraluminous infrared galaxies to higher FIR and CO luminosities.

We present the first cosmological simulations of dwarf galaxies, which include dark matter self-interactions and baryons. We study two dwarf galaxies within cold dark matter, and four different elastic self-interacting scenarios with constant and velocity-dependent cross sections, motivated by a new force in the hidden dark matter sector. Our highest resolution simulation has a baryonic mass resolution of $1.8\times 10^2\,{\rm M}_\odot$ and a gravitational softening length of $34\,{\rm pc}$ at $z=0$. In this first study we focus on the regime of mostly isolated dwarf galaxies with halo masses $\sim10^{10}\,{\rm M}_\odot$ where dark matter dynamically dominates even at sub-kpc scales. We find that while the global properties of galaxies of this scale are minimally affected by allowed self-interactions, their internal structures change significantly if the cross section is large enough within the inner sub-kpc region. In these dark-matter-dominated systems, self-scattering ties the shape of the stellar distribution to that of the dark matter distribution. In particular, we find that the stellar core radius is closely related to the dark matter core radius generated by self-interactions. Dark matter collisions lead to dwarf galaxies with larger stellar cores and smaller stellar central densities compared to the cold dark matter case. The central metallicity within $1\,{\rm kpc}$ is also larger by up to $\sim 15\%$ in the former case. We conclude that the mass distribution, and characteristics of the central stars in dwarf galaxies can potentially be used to probe the self-interacting nature of dark matter.

Using simulations we demonstrate that some of the published redshift distributions of submillimetre galaxies (SMGs) at different wavelengths, which were previously reported to be statistically different, are consistent with a parent distribution of the same population of galaxies. The redshift distributions which peak at z_med=2.9, 2.6, 2.2, 2.2, and 2.0 for galaxies selected at 2 and 1.1 mm, and 870, 850, and 450 um respectively, can be derived from a single parent redshift distribution, in contrast with previous studies. The differences can be explained through wavelength selection, depth of the surveys, and to a lesser degree, angular resolution. The main differences are attributed to the temperature of the spectral energy distributions, as shorter-wavelength maps select a hotter population of galaxies. Using the same parent distribution and taking into account lensing bias we can also reproduce the redshift distribution of 1.4 mm-selected ultra-bright galaxies, which peaks at z_med=3.4. However, the redshift distribution of 450 um-selected galaxies in the deepest surveys, which peaks at z_med=1.4, cannot be reproduced from the same parent population with just these selection effects. In order to explain this distribution we have to add another population of galaxies or include different selection biases.

The IceCube Neutrino Observatory has observed highly energetic neutrinos in excess of the expected atmospheric neutrino background. It is intriguing to consider the possibility that such events are probing physics beyond the standard model. In this context, $\mathcal{O}$(PeV) dark matter particles decaying to neutrinos have been considered while dark matter annihilation has been dismissed invoking the unitarity bound as a limiting factor. However, the latter claim was done ignoring the contribution from dark matter substructure, which for PeV Cold Dark Matter would extend down to a free streaming mass of $\mathcal{O}$($10^{-18}$M$_\odot$). Since the unitarity bound is less stringent at low velocities, ($\sigma_{\rm ann}$v)$\leq4\pi/m_\chi^2v$, then, it is possible that these cold and dense subhalos would contribute dominantly to a dark-matter-induced neutrino flux and easily account for the events observed by IceCube. A Sommerfeld-enhanced dark matter model can naturally support such scenario. Interestingly, the spatial distribution of the events shows features that would be expected in a dark matter interpretation. Although not conclusive, 9 of the 37 events appear to be clustered around a region near the Galactic Center while 6 others spatially coincide, within the reported angular errors, with 5 of 26 Milky Way satellites. However, a simple estimate of the probability of the latter occurring by chance is $\sim35\%$. More events are needed to statistically test this hypothesis. PeV dark matter particles are massive enough that their abundance as standard thermal relics would overclose the Universe. This issue can be solved in alternative scenarios, for instance if the decay of new massive unstable particles generates significant entropy reheating the Universe to a slightly lower temperature than the freeze-out temperature, $T_{\rm RH} \lesssim T_{\rm f}\sim4\times10^4$~GeV.

We present a novel perspective on the clustering of dark matter in phase space by defining the particle phase space average density ($P^2SAD$) as a two-dimensional extension of the two-point correlation function averaged within a certain volume in phase space. This statistics is a sensitive measure of small scale (sub-)structure of dark matter haloes. By analysing the structure of $P^2SAD$ in Milky-Way-size haloes using the Aquarius simulations, we find it to be nearly universal at small scales, i.e. small separations in phase space, where substructures dominate. This remarkable universality occurs across time and in regions of substantially different ambient densities (by nearly four orders of magnitude), with typical variations in $P^2SAD$ of a factor of a few. The maximum variations occur in regions where substructures have been strongly disrupted. The universality is also preserved across haloes of similar mass but diverse mass accretion histories and subhalo distributions. The universality is also broken at large scales, where the smooth dark matter distribution in the halo dominates. Although at small scales the structure of $P^2SAD$ is roughly described by a subhalo model, we argue that the simulation data is better fitted by a family of superellipse contours. This functional shape is inspired by a model that extends the stable clustering hypothesis into phase space. In a companion paper, we refine this model and show its advantages as a method to obtain predictions for non-gravitational signatures of dark matter.

We present a model for the structure of the particle phase space average density ($P^2SAD$) in galactic haloes, introduced recently as a novel measure of the clustering of dark matter. Our model is based on the stable clustering hypothesis in phase space, the spherical collapse model, and tidal disruption of substructures, which is calibrated against the Aquarius simulations. Using this model, we can predict the behaviour of $P^2SAD$ in the numerically unresolved regime, down to the decoupling mass limit of generic WIMP models. This prediction can be used to estimate signals sensitive to the small scale structure of dark matter. For example, the dark matter annihilation rate can be estimated for arbitrary velocity-dependent cross sections in a convenient way using a limit of $P^2SAD$ to zero separation in physical space. We illustrate our method by computing the global and local subhalo annihilation boost to that of the smooth dark matter distribution in a Milky-Way-size halo. Two cases are considered, one where the cross section is velocity independent and one that approximates Sommerfeld-enhanced models. We find that the global boost is $\sim10-30$, which is at the low end of current estimates (weakening expectations of large extragalactic signals), while the boost at the solar radius is below the percent level. We make our code to compute $P^2SAD$ publicly available, which can be used to estimate various observables that probe the nanostructure of dark matter haloes.

We demonstrate that self-interacting dark matter models with interactions mediated by light particles can have significant deviations in the matter power-spectrum and detailed structure of galactic halos when compared to a standard cold dark matter scenario. While these deviations can take the form of suppression of small scale structure that are in some ways similar to that of warm dark matter, the self-interacting models have a much wider range of possible phenomenology. A long-range force in the dark matter can introduce multiple scales to the initial power spectrum, in the form of dark acoustic oscillations and an exponential cut-off in the power spectrum. Using simulations we show that the impact of these scales can remain observationally relevant up to the present day. Furthermore, the self-interaction can continue to modify the small-scale structure of the dark matter halos, reducing their central densities and creating a dark matter core. The resulting phenomenology is unique to this type of models.

The mass aggregation and merger histories of present-day distinct haloes selected from the cosmological Millennium Simulations I and II are mapped into stellar mass aggregation and galaxy merger histories of central galaxies by using empirical stellar-to-halo and stellar-to-gas mass relations. The growth of bulges driven by the galaxy mergers/interactions is calculated using dynamical prescriptions. The predicted bulge demographics at redshift z~0 is consistent with observations (Zavala+2012). Here we present the evolution of the morphological mix (traced by the bulge-to-total mass ratio, B/T) as a function of mass up to z=3. This mix remains qualitatively the same up to z~1: B/T<0.1 galaxies dominate at low masses, 0.1<B/T<0.45 at intermediate masses, and B/T>0.45 at large masses. At z>1, the fractions of disc-dominated and bulgeless galaxies increase strongly, and by z~2 the era of pure disc galaxies is reached. Bulge-dominated galaxies acquire such a morphology, and most of their mass, following a downsizing trend. Since our results are consistent with most of the recent observational studies of the morphological mix at different redshifts, a LCDM-based scenario of merger-driven bulge assembly does not seem to face critical issues. However, if the stellar-to-halo mass relation changes too little with redshift, then some tensions with observations appear.

We combine the detailed Star Formation Histories of the Fornax and Sculptor dwarf Spheroidals with the Mass Assembly History of their dark matter (DM) halo progenitors to estimate if the energy deposited by Supernova type II (SNeII) is sufficient to create a substantial DM core. Assuming the efficiency of energy injection of the SNeII into DM particles is $\epsilon_{\rm gc}=0.05$, we find that a single early episode, $z \gtrsim z_{\rm infall}$, that combines the energy of all SNeII due to explode over 0.5 Gyr, is sufficient to create a core of several hundred parsecs in both Sculptor and Fornax. Therefore, our results suggest that it is energetically plausible to form cores in Cold Dark Matter (CDM) halos via early episodic gas outflows triggered by SNeII. Furthermore, based on CDM merger rates and phase-space density considerations, we argue that the probability of a subsequent complete regeneration of the cusp is small for a substantial fraction of dwarf-size haloes.

We study the effect of baryons on the abundance of structures and substructures in a Lambda-CDM cosmology, using a pair of high resolution cosmological simulations from the GIMIC project. Both simulations use identical initial conditions, but while one contains only dark matter, the other also includes baryons. We find that gas pressure, reionisation, supernova feedback, stripping, and truncated accretion systematically reduce the total mass and the abundance of structures below ~10^12 solar masses compared to the pure dark matter simulation. Taking this into account and adopting an appropriate detection threshold lowers the abundance of observed galaxies with maximum circular velocities below 100 km/s, significantly reducing the reported discrepancy between Lambda-CDM and the measured HI velocity function of the ALFALFA survey. We also show that the stellar-to-total mass ratios of galaxies with stellar masses of ~10^5 - 10^7 solar masses inferred from abundance matching of the (sub)halo mass function to the observed galaxy mass function increase by a factor of ~2. In addition, we find that an important fraction of low-mass subhaloes are completely devoid of stars. Accounting for the presence of dark subhaloes below 10^10 solar masses further reduces the abundance of observable objects, and leads to an additional increase in the inferred stellar-to-total mass ratio by factors of 2 - 10 for galaxies in haloes of 10^9 - 10^10 solar masses. This largely reconciles the abundance matching results with the kinematics of individual dwarf galaxies in Lambda-CDM. We propose approximate corrections to the masses of objects derived from pure dark matter calculations to account for baryonic effects.

The detailed origin of the diffuse gamma-ray background is still unknown. However, the contribution of unresolved sources is expected to induce small-scale anisotropies in this emission, which may provide a way to identify and constrain the properties of its contributors. Recent studies have predicted the contributions to the angular power spectrum (APS) from extragalactic and galactic dark matter (DM) annihilation or decay. The Fermi-LAT collaboration reported detection of angular power with a significance larger than $3\sigma$ in the energy range from 1 GeV to 10 GeV on 22 months of data [Ackermann et al. 2012]. For these preliminary results the already published Fermi-LAT APS measurements [Ackermann et al. 2012] are compared to the accurate predictions for DM anisotropies from state-of-the-art cosmological simulations as presented in [Fornasa et al. 2013] to derive constraints on different DM candidates.

Self-Interacting Dark Matter is an attractive alternative to the Cold Dark Matter paradigm only if it is able to substantially reduce the central densities of dwarf-size haloes while keeping the densities and shapes of cluster-size haloes within current constraints. Given the seemingly stringent nature of the latter, it was thought for nearly a decade that SIDM would be viable only if the cross section for self-scattering was strongly velocity-dependent. However, it has recently been suggested that a constant cross section per unit mass of sigma_T/m~0.1cm^2/g is sufficient to accomplish the desired effect. We explicitly investigate this claim using high resolution cosmological simulations of a Milky-Way size halo and find that, similarly to the Cold Dark Matter case, such cross section produces a population of massive subhaloes that is inconsistent with the kinematics of the classical dwarf spheroidals, in particular with the inferred slopes of the mass profiles of Fornax and Sculptor. This problem is resolved if sigma_T/m~1cm^2/g at the dwarf spheroidal scales. Since this value is likely inconsistent with the halo shapes of several clusters, our results leave only a small window open for a velocity-independent Self-Interacting Dark Matter model to work as a distinct alternative to Cold Dark Matter.

Self-interacting dark matter offers an interesting alternative to collisionless dark matter because of its ability to preserve the large-scale success of the cold dark matter model, while seemingly solving its challenges on small scales. We present here the first study of the expected dark matter detection signal taking into account different self-scattering models. We demonstrate that models with constant and velocity dependent cross sections, which are consistent with observational constraints, lead to distinct signatures in the velocity distribution, because non-thermalised features found in the cold dark matter distribution are thermalised through particle scattering. Depending on the model, self-interaction can lead to a 10% reduction of the recoil rates at high energies, corresponding to a minimum speed that can cause recoil larger than 300 km/s, compared to the cold dark matter case. At lower energies these differences are smaller than 5% for all models. The amplitude of the annual modulation signal can increase by up to 25%, and the day of maximum amplitude can shift by about two weeks with respect to the cold dark matter expectation. Furthermore, the exact day of phase reversal of the modulation signal can also differ by about a week between the different models. In general, models with velocity dependent cross sections peaking at the typical velocities of dwarf galaxies lead only to minor changes in the detection signals, whereas allowed constant cross section models lead to significant changes. We conclude that different self-interacting dark matter scenarios might be distinguished from each other through the details of direct detection signals. Furthermore, detailed constraints on the intrinsic properties of dark matter based on null detections, should take into account the possibility of self-scattering and the resulting effects on the detector signal.

We use the Millennium I and II cosmological simulations to revisit the impact of mergers in the growth of bulges in central galaxies in the LCDM scenario. We seed galaxies within the growing CDM haloes using semi-empirical relations to assign stellar and gaseous masses, and an analytic treatment to estimate the transfer of stellar mass to the bulge of the remnant after a galaxy merger. We find that this model roughly reproduces the observed correlation between the bulge-to-total (B/T) mass ratio and stellar mass in present-day central galaxies as well as their observed demographics, although low-mass B/T<0.1 (bulgeless) galaxies might be scarce relative to the observed abundance. In our merger-driven scenario, bulges have a composite population made of (i) stars acquired from infalling satellites, (ii) stars transferred from the primary disc due to merger-induced perturbations, and (iii) newly formed stars in starbursts triggered by mergers. We find that (i) and (ii) are the main channels of mass assembly, with the first being dominant for massive galaxies, creating large bulges with different stellar populations than those of the inner discs, while the second is dominant for intermediate/low-mass galaxies creating small bulges with similar stellar populations to the inner discs. We associate the dominion of the first (second) channel to classical (pseudo) bulges, and compare the predicted fractions to observations. We remark that our treatment does not include other mechanisms of bulge growth such as intrinsic secular disc instabilities or misaligned gas accretion. We find that the evolution of the stellar and gaseous contents of the satellite as it moves towards the central galaxy is a key ingredient in setting the morphology of the remnant, and that a good match to the observed bulge demographics occurs when this evolution proceeds closely to that of the central galaxy.

The existence of an extended hot gaseous corona surrounding clusters, groups and massive galaxies is well established by observational evidence and predicted by current theories of galaxy formation. When a small galaxy collides with a larger one, their coronae are the first to interact, producing disturbances that remove gas from the smaller system and settle it into the corona of the larger one. For a Milky-Way-size galaxy merging into a low-mass group, ram pressure stripping and the Kelvin-Helmholtz instability are the most relevant of these disturbances. We argue that the turbulence generated by the latter mixes the material of both coronae in the wake of the orbiting satellite creating a "warm phase" mixture with a cooling time a factor of several shorter than that of the ambient intragroup gas. We reach this conclusion using analytic estimates, as well as adiabatic and dissipative high resolution numerical simulations of a spherical corona subject to the ablation process of a constant velocity wind with uniform density and temperature. Although this is a preliminary analysis, our results are promising and we speculate that the mixture could potentially trigger in situ star formation and/or be accreted into the central galaxy as a cold gas flow resulting in a new mode of star formation in galaxy groups and clusters.

The Fermi-LAT collaboration has recently reported the detection of angular power above the photon noise level in the diffuse gamma-ray background between 1 and 50 GeV. Such signal can be used to constrain a possible contribution from Dark-Matter-induced photons. We estimate the intensity and features of the angular power spectrum (APS) of this potential Dark Matter (DM) signal, for both decaying and annihilating DM candidates, by constructing template all-sky gamma-ray maps for the emission produced in the galactic halo and its substructures, as well as in extragalactic (sub)halos. The DM distribution is given by state-of-the-art N-body simulations of cosmic structure formation, namely Millennium-II for extragalactic (sub)halos, and Aquarius for the galactic halo and its subhalos. We use a hybrid method of extrapolation to account for (sub)structures that are below the resolution limit of the simulations, allowing us to estimate the total emission all the way down to the minimal self-bound halo mass. We describe in detail the features appearing in the APS of our template maps and we estimate the effect of various uncertainties such as the value of the minimal halo mass, the fraction of substructures hosted in a halo and the shape of the DM density profile. Our results indicate that the fluctuation APS of the DM-induced emission is of the same order as the Fermi-LAT APS, suggesting that one can constrain this hypothetical emission from the comparison with the measured anisotropy. We also quantify the uncertainties affecting our results, finding "theoretical error bands" spanning more than two orders of magnitude and dominated (for a given particle physics model) by our lack of knowledge of the abundance of low-mass (sub)halos.

We present N-body simulations of a new class of self-interacting dark matter models, which do not violate any astrophysical constraints due to a non-power-law velocity dependence of the transfer cross section which is motivated by a Yukawa-like new gauge boson interaction. Specifically, we focus on the formation of a Milky Way-like dark matter halo taken from the Aquarius project and re-simulate it for a couple of representative cases in the allowed parameter space of this new model. We find that for these cases, the main halo only develops a small core (~1 kpc) followed by a density profile identical to that of the standard cold dark matter scenario outside of that radius. Neither the subhalo mass function nor the radial number density of subhaloes are altered in these models but there is a significant change in the inner density structure of subhaloes resulting in the formation of a large density core. As a consequence, the inner circular velocity profiles of the most massive subhaloes differ significantly from the cold dark matter predictions and we demonstrate that they are compatible with the observational data of the brightest Milky Way dSphs in such a velocity-dependent self-interacting dark matter scenario. Specifically, and contrary to the cold dark matter case, there are no subhaloes that are more concentrated than what is inferred from the kinematics of the Milky Way dSphs. We conclude that these models offer an interesting alternative to the cold dark matter model that can reduce the recently reported tension between the brightest Milky Way satellites and the dense subhaloes found in cold dark matter simulations.

For the first time, the Fermi-LAT measured the angular power spectrum (APS) of anisotropies in the diffuse gamma-ray background. The data is found to be broadly compatible with a model with contributions from the point sources in the 1-year catalog, the galactic diffuse background, and the extragalactic isotropic emission; however deviations are present at both large and small angular scales. In this study, we complement the model with a contribution from Dark Matter (DM) whose distribution is modeled exploiting the results of the most recent N-body simulations, considering the contribution of extragalactic halos and subhalos (from Millennium-II) and of galactic substructures (from Aquarius). With the use of the Fermi Science Tools, these simulations serve as templates to produce mock gamma-ray count maps for DM gamma-ray emission, both in the case of an annihilating and a decaying DM candidate. The APS will then be computed and compared with the Fermi-LAT results to derive constraints on the DM particle physics properties. The possible systematic due to an imperfect model of the galactic foreground is also studied and taken into account properly. The present paper reports on the status of the project.

(Abridged) The extragalactic background light (EBL) observed at multiple wavelengths is a promising tool to probe the nature of dark matter since it might contain a significant contribution from gamma-rays produced promptly by dark matter annihilation. Additionally, the electrons and positrons produced in the annihilation give energy to the CMB photons to populate the EBL with X-rays and gamma-rays. We here create full-sky maps of the radiation from both of these contributions using the high-resolution Millennium-II simulation. We use upper limits on the contributions of unknown sources to the EBL to constrain the intrinsic properties of dark matter using a model-independent approach that can be employed as a template to test different particle physics models (including those with a Sommerfeld enhancement). These upper limits are based on observations spanning eight orders of magnitude in energy (from soft X-rays measured by CHANDRA to gamma-rays measured by Fermi), and on expectations for the contributions from blazars and star forming galaxies. To exemplify this approach, we analyze a set of benchmark Sommerfeld-enhanced models that give the correct dark matter abundance, satisfy CMB constraints, and fit the cosmic ray spectra measured by PAMELA and Fermi without any contribution from local subhalos. We find that these models are in conflict with the EBL constraints unless the contribution of unresolved subhalos is small and the annihilation signal dominates the EBL. We conclude that provided the collisionless cold dark matter paradigm is accurate, even for conservative estimates of the contribution from unresolved subhalos and astrophysical backgrounds, the EBL is at least as sensitive a probe of these types of scenarios as the CMB. Our results disfavor an explanation of the positron excess measured by PAMELA based only on dark matter annihilation in the smooth Galactic halo.

We calculate how the relic density of dark matter particles is altered when their annihilation is enhanced by the Sommerfeld mechanism due to a Yukawa interaction between the annihilating particles. Maintaining a dark matter abundance consistent with current observational bounds requires the normalization of the s-wave annihilation cross section to be decreased compared to a model without enhancement. The level of suppression depends on the specific parameters of the particle model, with the kinetic decoupling temperature having the most effect. We find that the cross section can be reduced by as much as an order of magnitude for extreme cases. We also compute the mu-type distortion of the CMB energy spectrum caused by energy injection from such Sommerfeld-enhanced annihilation. Our results indicate that in the vicinity of resonances, associated with bound states, distortions can be large enough to be excluded by the upper limit |mu|<9.0x10^(-5) found by the COBE/FIRAS experiment.

If dark matter is composed of neutralinos, one of the most exciting prospects for its detection lies in observations of the gamma-ray radiation created in pair annihilations between neutralinos, a process that may contribute significantly to the extragalactic gamma-ray background (EGB) radiation. We here use the high-resolution Millennium-II simulation of cosmic structure formation to produce the first full-sky maps of the expected radiation coming from extragalactic dark matter structures. Our map making procedure takes into account the total gamma-ray luminosity from all haloes and their subhaloes, and includes corrections for unresolved components of the emission as well as an extrapolation to the damping scale limit of neutralinos. Our analysis also includes a proper normalization of the signal according to a specific supersymmetric model based on minimal supergravity. The new simulated maps allow a study of the angular power spectrum of the gamma-ray background from dark matter annihilation, which has distinctive features associated with the nature of the annihilation process and may be detectable in forthcoming observations by the recently launched FERMI satellite. Our results are in broad agreement with analytic models for the gamma-ray background, but they also include higher-order correlations not readily accessible in analytic calculations and, in addition, provide detailed spectral information for each pixel. In particular, we find that difference maps at different energies can reveal cosmic large-scale structure at low and intermediate redshifts. If the intrinsic emission spectrum is characterized by an emission peak, cosmological tomography with gamma ray annihilation radiation is in principle possible.

If dark matter is composed of neutralinos, the gamma-ray radiation produced in their annihilation offers an attractive possibility for dark matter detection. This process may contribute significantly to the extragalactic gamma-ray background (EGB) radiation, which is being measured by the FERMI satellite with unprecedented sensitivity. Using the high-resolution Millennium-II simulation of cosmic structure formation we have produced the first full-sky maps of the expected contribution of dark matter annihilation to the EGB radiation. Our maps include a proper normalization of the signal according to a specific supersymmetric model based on minimal supergravity. The new simulated maps allow a study of the angular power spectrum of the gamma-ray background from dark matter annihilation, which has distinctive features associated with the nature of the annihilation process. Our results are in broad agreement with analytic models for the gamma-ray background, but they also include higher-order correlations not readily accessible in analytic calculations and, in addition, provide detailed spectral information for each pixel. In particular, we find that color maps combining different energies can reveal the cosmic large-scale structure at low and intermediate redshifts.

Using constrained simulations of the local Universe for generic cold dark matter and for 1keV warm dark matter, we investigate the difference in the abundance of dark matter halos in the local environment. We find that the mass function within 20 Mpc/h of the Local Group is ~2 times larger than the universal mass function in the 10^9-10^13 M_odot/h mass range. Imposing the field of view of the on-going HI blind survey ALFALFA in our simulations, we predict that the velocity function in the Virgo-direction region exceeds the universal velocity function by a factor of 3. Furthermore, employing a scheme to translate the halo velocity function into a galaxy velocity function, we compare the simulation results with a sample of galaxies from the early catalog release of ALFALFA. We find that our simulations are able to reproduce the velocity function in the 80-300 km/s velocity range, having a value ~10 times larger than the universal velocity function in the Virgo-direction region. In the low velocity regime, 35-80 km/s, the warm dark matter simulation reproduces the observed flattening of the velocity function. On the contrary, the simulation with cold dark matter predicts a steep rise in the velocity function towards lower velocities; for V_max=35 km/s, it forecasts ~10 times more sources than the ones observed. If confirmed by the complete ALFALFA survey, our results indicate a potential problem for the cold dark matter paradigm or for the conventional assumptions about energetic feedback in dwarf galaxies.

We explore how the slopes and scatters of the scaling relations of disk galaxies (Vm-L[-M], R-L[-M], and Vm-R) do change when moving from B to K bands and to stellar and baryonic quantities. For our compiled sample of 76 normal, non-interacting high and low surface brightness galaxies, we find some changes, which evidence evolution effects, mainly related to gas infall and star formation (SF). We also explore correlations among the (B-K) color, stellar mass fraction fs, mass M (luminosity L), and surface density (SB), as well as correlations among the residuals of the scaling relations. Some of our findings are: (i) the scale length Rb is a third parameter in the baryonic TF relation and the residuals of this relation follow a trend (slope ~-0.15) with the residuals of the Rb-Mb relation; for the stellar and K band cases, R is not anymore a third parameter and the mentioned trend disappears; (ii) among the TFRs, the B-band TFR is the most scattered; in this case, the color is a third parameter; (iii) the LSB galaxies break some observed trends, which suggest a threshold in the gas surface density Sg, below which the SF becomes independent of the gas infall rate and Sg. Our results are interpreted and discussed in the light of LCDM-based models of galaxy evolution. The models explain not only the baryonic scaling relations, but also most of the processes responsible for the observed changes in the slopes, scatters, and correlations among the residuals when changing to stellar and luminous quantities. The baryon fraction is required to be smaller than 0.05 on average. We detect some potential difficulties for the models: the observed color-M and surface density-M correlations are steeper, and the intrinsic scatter in the baryonic TFR is smaller than those predicted. [abridged]

We present 3D hydrodynamical simulations of a precessing jet propagating inside a supernova remnant (SNR) shell, particularly applied to the W50-SS433 system in a search for the origin of its peculiar elongated morphology. Several runs were carried out with different values for the mass loss rate of the jet, the initial radius of the SNR, and the opening angle of the precession cone. We found that our models successfully reproduce the scale and morphology of W50 when the opening angle of the jets is set to 10$\degr$ or if this angle linearly varies with time. For these models, more realistic runs were made considering that the remnant is expanding into an interstellar medium (ISM) with an exponential density profile (as HI observations suggest). Taking into account all these ingredients, the large scale morphology of the W50-SS 433 system, including the asymmetry between the lobes (formed by the jet-SNR interaction), is well reproduced.

We derive an expression for the entropy of a dark matter halo described by a Navarro-Frenk-White model with a core. The comparison of this entropy with the one of dark matter at the freeze-out era allows us to constraint the parameter space in mSUGRA models. Moreover, joining these constraints with the ones obtained from the usual abundance criteria and demanding both criteria to be consistent with the 2$\sigma$ bounds for the abundance of dark matter: $0.112\leq\Omega_{DM}h^2\leq0.122$, we are able to clearly discriminate validity regions among the values of $\tan \beta$, one of the parameters of the mSUGRA model. We found that for the explored regions of the parameter space, small values of tan$\beta$ are not favored; only for tan$\beta\simeq50$ are both criteria significantly consistent. In the region where both criteria are consistent we also found a lower bound for the neutralino mass, $m_{\chi}\geq 141$ GeV.

We investigate the evolution of angular momentum in simulations of galaxy formation in a cold dark matter universe. We analyse two model galaxies generated in the N-body/hydrodynamic simulations of Okamoto et al. Starting from identical initial conditions, but using different assumptions for the baryonic physics, one of the simulations produced a bulge-dominated galaxy and the other one a disc-dominated galaxy. The main difference is the treatment of star formation and feedback, both of which were designed to be more efficient in the disc-dominated object. We find that the specific angular momentum of the disc-dominated galaxy tracks the evolution of the angular momentum of the dark matter halo very closely: the angular momentum grows as predicted by linear theory until the epoch of maximum expansion and remains constant thereafter. By contrast, the evolution of the angular momentum of the bulge-dominated galaxy resembles that of the central, most bound halo material: it also grows at first according to linear theory, but 90% of it is rapidly lost as pre-galactic fragments, into which gas had cooled efficiently, merge, transferring their orbital angular momentum to the outer halo by tidal effects. The disc-dominated galaxy avoids this fate because the strong feedback reheats the gas which accumulates in an extended hot reservoir and only begins to cool once the merging activity has subsided. Our analysis lends strong support to the classical theory of disc formation whereby tidally torqued gas is accreted into the centre of the halo conserving its angular momentum.

We derive an expression for the entropy of a present dark matter halo described by a Navarro-Frenk-White modified model with a central core. The comparison of this entropy with the one of the halo at the freeze-out era allows us to obtain an expression for the relic abundance of neutralinos, which in turn is used to constrain the parameter space in mSUGRA models, when used with the WMAP observations. Moreover, by joning these results with the ones obtained from the usual abundance criteria, we are able to clearly discriminate validity regions among tan beta values of the mSUGRA model, by demanding both criteria to be consistent with the 2 sigma bounds of the WMAP observations for the relic density: 0.112 < Omega h^2 < 0.122. We found that for sign mu positive, small values of tan beta are not favored; only for tan beta ~ 50 are both criteria significantly consistent. The use of both criteria also allows us to put a lower bound on the neutralino mass of > 151 GeV.

Motivated by the possible conflict between the Navarro-Frenk-White(NFW) model predictions for the dark matter contents of galactic systems and its correlation with baryonic surface density, we will explore an alternative paradigm for the description of dark matter halos. Such an alternative emerges from Tsallis' non-extensive thermodynamics applied to self-gravitating systems and leads to the so-called ``stellar polytrope'' (SP) model. We consider that this could be a better approach to real structures rather than the isothermal model, given the fact that the first one takes into account the non-extensivity of energy and entropy present in these type of systems characterized by long-range interactions. We compare a halo based on the Navarro-Frenk-White (NFW) and one which follows the SP description. Analyzing the dark matter contents estimated by means of global physical parameters of galactic disks, obtained from a sample of actual galaxies, with the ones of the unobserved dark matter halos, we conclude that the SP model is favored over the NFW model in such a comparison.

We explore the use of two criteria to constraint the allowed parameter space in mSUGRA models. Both criteria are based in the calculation of the present density of neutralinos as dark matter in the Universe. The first one is the usual ``abundance'' criterion which is used to calculate the relic density after the ``freeze-out'' era. To compute the relic density we used the numerical public code micrOMEGAs. The second criterion applies the microcanonical definition of entropy to a weakly interacting and self-gravitating gas evaluating then the change in the entropy per particle of this gas between the ``freeze-out'' era and present day virialized structures. An ``entropy-consistency'' criterion emerges by comparing theoretical and empirical estimates of this entropy. The main objective of our work is to determine for which regions of the parameter space in the mSUGRA model are both criteria consistent with the 2$\sigma$ bounds according to WMAP for the relic density: $0.0945<\Omega_{CDM}h^2<0.1287$. As a first result, we found that for $A_0=0$, sgn$\mu=+$, small values of tan$\beta$ are not favored; only for tan$\beta\simeq50$ are both criteria significantly consistent.

We study a dark matter halo model from two points of view: the ``stellar polytrope'' (SP) model coming from Tsallis' thermodynamics, and the one coming from the Navarro-Frenk-White (NFW) paradigm. We make an appropriate comparison between both halo models and analyzing the relations between the global physical parameters of observed galactic disks, coming from a sample of actual galaxies, with the ones of the unobserved dark matter halos, we conclude that the SP model is favored over the NFW model in such a comparison.

In this work we study the gravitational influence of the material extending from Uranus orbit to the Kuiper belt and beyond on objects moving within these regions. We conclude that a density distribution given by $\rho(r)=\frac{1}{r}$ (for $r\geq 20 UA$) generates a constant acceleration towards the Sun on those objects, which, with the proper amount of mass, accounts for the blue shift detected on the Pioneers space crafts. We also discuss the effect of this gravitational pull on Neptune, and comment on the possible origin of such a matter distribution.

We present an alternative for the description of galactic halos based on Tsallis' non--extensive entropy formalism; on this scheme, halos are stellar polytropes characterized by three parameters, the central density, $\rho_c$, the central velocity dispersion, $\sigma_c$ and the polytropic index, $n$. To evaluate these parameters we take the Navarro-Frenk-White paradigm as a comparative model and make the following assumptions: both halo models must have the same virial mass, the same total energy and the same maximal velocity. These three conditions fix all the parameters for a given stellar polytrope allowing us to compare both halo models. The halos studied have virial masses on the range $10^{12}-10^{15} M_\odot$, and it was found after the analysis that they are described, at all scales, by almost the same polytropic index, $n\approx 4.8$, implying an empirical estimation of Tsallis non--extensive parameter for this type of dynamical systems: $q\approx 1.3$.

For a compiled sample of disk galaxies with available photometry (B and K bands), velocity line-widths and HI integral fluxes, several parameters which trace the luminous, baryonic and dark matter contents were inferred. We investigated how these parameters do vary with different galaxy properties, and confronted the results with predictions of galaxy evolutionary models in the context of the LCDM cosmogony. The ratio of disk-to-total maximum circular velocity, vd/vt, depends mainly on the central disk surface density Sig_d (or surface brightness, SB), increasing roughly as Sig_d^0.15. While a fraction of high SB galaxies have a vd/vt ratio corresponding to the maximum disk solution, the low SB are completely dark matter dominated. The trend is similar for the models, although they have slightly smaller vd/vt ratios than observations. An analysis of residuals of the vd/vt-Sig_d relation shows that vd/vt tends to decrease as the galaxy is redder, more luminous (massive), and of earlier type. The models allow us to explain the physics of these results, which imply a connexion between halo structure and luminous properties. The estimated dynamical-to-baryon mass and dynamical mass-to-light ratios at a given radius, for observations and models, decrease with the SB and do not correlate with the galaxy scale, contrary to what has been reported in previous works, based on the analysis of rotation curve shapes. We discuss this difference and state the importance to solve the controversy on whether the dark and luminous contents in disk galaxies depend on SB or luminosity. The broad agreement between models and observations favors the LCDM scenario. However, the excess of dark matter inside the optical region of disk galaxies remains as the main difficulty. (abridged)

We discuss the shape and decomposition of rotation curves (RCs) of galaxies formed within growing cold dark matter halos. The outer RC shape correlates mainly with the surface brightness (SB), the luminous mass fraction, fd, and the bulge fraction. In order the shapes of RC depend significantly on luminosity, fd should be a strong function of mass (feedback). For the preferred values of fd (\lesssim 0.03), the high SB models can be maximum disks only when the halos have a shallow core. The low SB models are sub-maximum disks. The residuals of the baryonic Tully-Fisher (TF) and disk mass-radius relations show a clear anti-correlation among them, but when one passes to the TF and luminosity-radius relations, the anti-correlation almost disappears. Therefore, the observed lack of correlation among the residuals of the last two relations should not be interpreted as an evidence of sub-maximal disks.