Search

The existence of galaxies with a surface brightness $\mu$ lower than the night sky has been known since three decades. Yet, their formation mechanism and emergence within a $\rm\Lambda CDM$ universe has remained largely undetermined. For the first time, we investigated the origin of Low Surface Brightness (LSB) galaxies with M$_{\star}$$\sim$10$^{9.5-10}$M$_{\odot}$, which we are able to reproduce within hydrodynamical cosmological simulations from the NIHAO suite. The simulated and observed LSBs share similar properties, having large HI reservoir, extended star formation histories and effective radii, low S\'{e}rsic index and slowly rising rotation curves. The formation mechanism of these objects is explored: simulated LSBs form as a result of co-planar co-rotating mergers and aligned accretion of gas at early times, while perpendicular mergers and mis-aligned gas accretion result in higher $\mu$ galaxies by $z$=0. The larger the merger, the stronger the correlation between merger orbital configuration and final $\mu$. While the halo spin parameter is consistently high in simulated LSB galaxies, the impact of halo concentration, feedback-driven gas outflows and merger time only plays a minor-to-no role in determining $\mu$. Interestingly, the formation scenario of such `classical' LSBs differs from the one of less massive, M$_{\star}$$\sim$10$^{7-9}$M$_{\odot}$, Ultra-Diffuse Galaxies, the latter resulting from the effects of SNae driven gas outflows: a M$_{\star}$ of $\sim$10$^9$M$_{\odot}$ thus represents the transition regime between a feedback-dominated to an angular momentum-dominated formation scenario in the LSB realm. Observational predictions are offered regarding spatially resolved star formation rates through LSB discs: these, together with upcoming surveys, can be used to verify the proposed emergence scenario of LSB galaxies., Comment: MNRAS Accepted 2019 April 02; in original form 2019 January 24. 15 pages,17 figures. main results in Fig 12 and 14. Added appendix A

The GERmanium Detector Array (GERDA) is a low background experiment located at the Laboratori Nazionali del Gran Sasso in Italy, which searches for neutrinoless double beta decay of $^{76}$Ge into $^{76}$Se+2e$^-$. GERDA has been conceived in two phases. Phase II, which started in December 2015, features several novelties including 30 new Ge detectors. These were manufactured according to the Broad Energy Germanium (BEGe) detector design that has a better background discrimination capability and energy resolution compared to formerly widely-used types. Prior to their installation, the new BEGe detectors were mounted in vacuum cryostats and characterized in detail in the HADES underground laboratory in Belgium. This paper describes the properties and the overall performance of these detectors during operation in vacuum. The characterization campaign provided not only direct input for GERDA Phase II data collection and analyses, but also allowed to study detector phenomena, detector correlations as well as to test the strength of pulse shape simulation codes., Comment: 29 pages, 18 figures

State-of-the-art cosmological hydrodynamical simulations of galaxy formation have reached the point at which their outcomes result in galaxies with ever more realism. Still, the employed sub-grid models include several free parameters such as the density threshold, $n$, to localize the star-forming gas. In this work, we investigate the possibilities to utilize the observed clustered nature of star formation (SF) in order to refine SF prescriptions and constrain the density threshold parameter. To this end, we measure the clustering strength, correlation length and power-law index of the two-point correlation function of young ($\tau<50$ Myr) stellar particles and compare our results to observations from the HST Legacy Extragalactic UV Survey (LEGUS). Our simulations reveal a clear trend of larger clustering signal and power-law index and lower correlation length as the SF threshold increases with only mild dependence on galaxy properties such as stellar mass or specific star formation rate. In conclusion, we find that the observed clustering of SF is inconsistent with a low threshold for SF ($n<1$ cm$^{-3}$) and strongly favours a high value for the density threshold of SF ($n>10$ cm$^{-3}$), as for example employed in the NIHAO project., Comment: accepted by MNRAS, 7 pages, 5 figures, key figure: figure 5

We use cosmological hydrodynamical galaxy formation simulations from the NIHAO project to investigate the impact of the threshold for star formation on the response of the dark matter (DM) halo to baryonic processes. The fiducial NIHAO threshold, $n=10\, {\rm cm}^{-3}$, results in strong expansion of the DM halo in galaxies with stellar masses in the range $10^{7.5} < M_{star} < 10^{9.5} M_{\odot}$. We find that lower thresholds such as $n=0.1$ (as employed by the EAGLE/APOSTLE and Illustris/AURIGA projects) do not result in significant halo expansion at any mass scale. Halo expansion driven by supernova feedback requires significant fluctuations in the local gas fraction on sub-dynamical times (i.e., < 50 Myr at galaxy half-light radii), which are themselves caused by variability in the star formation rate. At one per cent of the virial radius, simulations with $n=10$ have gas fractions of $\simeq 0.2$ and variations of $\simeq 0.1$, while $n=0.1$ simulations have order of magnitude lower gas fractions and hence do not expand the halo. The observed DM circular velocities of nearby dwarf galaxies are inconsistent with CDM simulations with $n=0.1$ and $n=1$, but in reasonable agreement with $n=10$. Star formation rates are more variable for higher $n$, lower galaxy masses, and when star formation is measured on shorter time scales. For example, simulations with $n=10$ have up to 0.4 dex higher scatter in specific star formation rates than simulations with $n=0.1$. Thus observationally constraining the sub-grid model for star formation, and hence the nature of DM, should be possible in the near future., Comment: 18 pages, 13 figures, accepted to MNRAS

The NIHAO cosmological simulations form a collection of a hundred high-resolution galaxies. We used these simulations to test the impact of stellar feedback on the morphology of the HI distribution in galaxies. We ran a subsample of twenty of the galaxies with different parameterizations of stellar feedback, looking for differences in the HI spatial distribution and morphology. We found that different feedback models do leave a signature in HI, and can potentially be compared with current and future observations. These findings can help inform future modeling efforts in the parameterization of stellar feedback in cosmological simulations of galaxy formation and evolution., Comment: Conference proceedings presented at the First Sharjah International Conference on Particle Physics, Astrophysics, and Cosmology

We used high resolution dark matter only cosmological simulations to investigate the structural properties of Lambda Cold Dark Matter ($\Lambda$CDM) haloes over cosmic time. The haloes in our study range in mass from $\sim 10^{10}$ to $\sim 10^{12} \mathrm{M}_\odot$, and are resolved with $10^5$ to $10^7$ particles. We fit the spherically averaged density profiles of DM haloes with the three parameter Einasto function. For our sample of haloes, the Einasto shape parameter, $\alpha$, is uncorrelated with the concentration, $c$, at fixed halo mass, and at all redshifts. Previous reports of an anti-correlation are traced to fitting degeneracies, which our fits are less sensitive to due to our higher spatial resolution. However, for individual haloes the evolution in $\alpha$ and $c$ is anti-correlated: at redshift $z=7$, $\alpha \simeq 0.4$ and decreases with time, while $c\simeq 3$ and increases with time. The evolution in structure is primarily due to accretion of mass at larger radii. We suggest that $\alpha$ traces the evolutionary state of the halo, with dynamically young haloes having high $\alpha$ (closer to a top-hat: $\alpha^{-1}=0$), and dynamically relaxed haloes having low $\alpha$ (closer to isothermal: $\alpha=0$). Such an evolutionary dependence reconciles the increase of $\alpha$ vs peak height, $\nu$, with the dependence on the slope of the power spectrum of initial density fluctuations found by previous studies., Comment: accepted for publication in MNRAS

We study the origin of bulge stars and their angular momentum (AM) evolution in 10 spiral galaxies with baryonic masses above $10^{10}$M$_\odot$ in the NIHAO galaxy formation simulations. The simulated galaxies are in good agreement with observations of the relation between specific AM and mass of the baryonic component and the stellar bulge-to-total ratio ($B/T$). We divide the star particles at $z=0$ into disc and bulge components using a hybrid photometric/kinematic decomposition method that identifies all central mass above an exponential disc profile as the `bulge'. By tracking the bulge star particles back in time, we find that on average 95\% of the bulge stars formed {\it in situ}, 3\% formed {\it ex situ} in satellites of the same halo, and only 2\% formed {\it ex situ} in external galaxies. The evolution of the AM distribution of the bulge stars paints an interesting picture: the higher the final $B/T$ ratio, the more the specific AM remains preserved during the bulge formation. In all cases, bulge stars migrate significantly towards the central region, reducing their average galactocentric radius by roughly a factor 2, independently of the final $B/T$ value. However, in the higher $B/T$ ($\gtrsim0.2$) objects, the velocity of the bulge stars increases and the AM of the bulge is almost conserved, whereas at lower $B/T$ values, the velocity of the bulge stars decreases and the AM of bulge reduces. The correlation between the evolution of the AM and $B/T$ suggests that bulge and disc formation are closely linked and cannot be treated as independent processes., Comment: 17 pages, 16 Figures, 1 table; accepted for publication in MNRAS

We analyze 88 independent high-resolution cosmological zoom-in simulations of disk galaxies in the NIHAO simulations suite to explore the connection between the atomic gas fraction and angular momentum of baryons throughout cosmic time. The study is motivated by the analytic model of \citet{obreschkow16}, which predicts a relation between the atomic gas fraction $f_{\rm atm}$ and the global atomic stability parameter $q \equiv j\sigma / (GM)$, where $M$ and $j$ are the mass and specific angular momentum of the galaxy (stars+cold gas) and $\sigma$ is the velocity dispersion of the atomic gas. We show that the simulated galaxies follow this relation from their formation ($z\simeq4$) to present within $\sim 0.5$ dex. To explain this behavior, we explore the evolution of the local Toomre stability and find that $90\%$--$100\%$ of the atomic gas in all simulated galaxies is stable at any time. In other words, throughout the entire epoch of peak star formation until today, the timescale for accretion is longer than the timescale to reach equilibrium, thus resulting in a quasi-static equilibrium of atomic gas at any time. Hence, the evolution of $f_{\rm atm}$ depends on the complex hierarchical growth history primarily via the evolution of $q$. An exception are galaxies subject to strong environmental effects., Comment: 12 pages, 7 figures; accepted to ApJ

We use 85 pairs of high resolution LCDM cosmological simulations from the NIHAO project to investigate why in dwarf galaxies neutral hydrogen (HI) linewidths measured at 50% of the peak flux W_50 /2 (from the hydrodynamical simulations) tend to underpredict the maximum circular velocity VmaxDMO (from the corresponding dark matter only simulations). There are two main contributing processes. 1) Lower mass galaxies are less rotationally supported. This is confirmed observationally from the skewness of linewidths in bins of HI mass in both ALFALFA and HIPASS observations. 2) The HI distributions are less extended (relative to the dark matter halo) in dwarf galaxies. Coupled to the lower baryon-to-halo ratio results in rotation curves that are still rising at the last measured data point, in agreement with observations from SPARC. Combining these two effects, in both simulations and observations lower mass galaxies have on average lower W_50 / W_20. Additionally, mass loss driven by supernovae and projection effects (dwarf galaxies are in general not thin disks) further reduce the linewidths. The implied HI linewidth velocity function from NIHAO is in good agreement with observations in the nearby Universe of dwarf galaxies: 10 < W_50 /2 < 80 km/s. The dark matter only slope of -2.9 is reduced to -1.0 in the hydro simulations. Future radio observations of unbiased samples with higher spatial resolution will enable stricter tests of the simulations, and thus of the LCDM model., Comment: 21 pages, 17 figures, accepted to MNRAS

The stellar populations in the inner kiloparsecs of the Milky Way (MW) show complex kinematical and chemical structures. The origin and evolution of these structures is still under debate. Here we study the central region of a fully cosmological hydrodynamical simulation of a disk galaxy that reproduces key properties of the inner kiloparsecs of the MW: it has a boxy morphology and shows an overall rotation and dispersion profile in agreement with observations. We use a clustering algorithm on stellar kinematics to identify a number of discrete kinematic components: a high- and low-spin disk, a stellar halo and two bulge components; one fast rotating and one slow-rotating. We focus on the two bulge components and show that the slow rotating one is spherically symmetric while the fast rotating component shows a boxy/peanut morphology. Although the two bulge components are kinematically discrete populations at present-day, they are both mostly formed over similar time scales, from disk material. We find that stellar particles with lower initial birth angular momentum (most likely thick disc stars) end up in the slow-rotating low-spin bulge, while stars with higher birth angular momentum (most likely thin disc stars) are found in the high-spin bulge. This has the important consequence that a bulge population with a spheroidal morphology does not necessarily indicate a merger origin. In fact, we do find that only $\sim2.3$\% of the stars in the bulge components are ex-situ stars brought in by accreted dwarf galaxies early on. We identify these ex-situ stars as the oldest and most metal-poor stars on highly radial orbits with large vertical excursions from the disk., Comment: accepted by ApJ, 14 pages (plus 3 pages appendix), 13 + 8 figures

The similarity between the distributions of spins for galaxies ($\lambda_{\rm g}$) and for dark-matter haloes ($\lambda_{\rm h}$), indicated both by simulations and observations, is naively interpreted as a one-to-one correlation between the spins of a galaxy and its host halo. This is used to predict galaxy sizes in semi-analytic models via $R_{\rm e}\simeq\lambda_{\rm h} R_{\rm v}$, with $R_{\rm e}$ the half-mass radius of the galaxy and $R_{\rm v}$ the halo radius. Utilizing two different suites of zoom-in cosmological simulations, we find that $\lambda_{\rm g}$ and $\lambda_{\rm h}$ are in fact only barely correlated, especially at $z\geq 1$. A general smearing of this correlation is expected based on the different spin histories, where the more recently accreted baryons through streams gain and then lose significant angular momentum compared to the gradually accumulated dark matter. Expecting the spins of baryons and dark matter to be correlated at accretion into $R_{\rm v}$, the null correlation at the end reflects an anti-correlation between $\lambda_{\rm g}/\lambda_{\rm h}$ and $\lambda_{\rm h}$, which can partly arise from mergers and a compact star-forming phase that many galaxies undergo. On the other hand, the halo and galaxy spin vectors tend to be aligned, with a median $\cos\theta=0.6$-0.7 between galaxy and halo, consistent with instreaming within a preferred plane. The galaxy spin is better correlated with the spin of the inner halo, but this largely reflects the effect of the baryons on the halo. Following the null spin correlation, $\lambda_{\rm h}$ is not a useful proxy for $R_{\rm e}$. While our simulations reproduce a general relation of the sort $R_{\rm e}=AR_{\rm vir}$, in agreement with observational estimates, the relation becomes tighter with $A=0.02(c/10)^{-0.7}$, where $c$ is the halo concentration, which in turn introduces a dependence on mass and redshift., Comment: 20 pages, 16 figures, submitted to MNRAS

We use 25 simulated galaxies from the NIHAO project to define and characterize a variety of kinematic stellar structures: thin and thick discs, large scale single discs, classical and pseudo bulges, spheroids, inner discs, and stellar haloes. These structures have masses, spins, shapes and rotational support in good agreement with theoretical expectations and observational data. Above a dark matter halo mass of $2.5\times10^{\rm~11}M_{\rm\odot}$, all galaxies have a classical bulge and 70\% have a thin and thick disc. The kinematic (thin) discs follow a power-law relation between angular momentum and stellar mass $J_{\rm *}=3.4M_{\rm *}^{\rm1.26\pm0.06}$, in very good agreement with the prediction based on the empirical stellar-to-halo mass relation in the same mass range, and show a strong correlation between maximum `observed' rotation velocity and dark matter halo circular velocity $v_{\rm c}=6.4v_{\rm max}^{0.64\pm0.04}$. Tracing back in time these structures' progenitors, we find all to lose a fraction $1-f_j$ of their maximum angular momentum. Thin discs are significantly better at retaining their high-redshift spins ($f_j\sim0.70$) than thick ones ($f_j\sim0.40$). Stellar haloes have their progenitor baryons assembled the latest ($z_{\rm~1/2}\sim1.1$) and over the longest timescales ($\tau\sim6.2$~Gyr), and have the smallest fraction of stars born in-situ ($f_{\rm in-situ}=0.35\pm0.14$). All other structures have $1.5\lesssim z_{\rm1/2}\lesssim3$, $\tau=4\pm2$~Gyr and $f_{\rm in-situ}\gtrsim0.9$., Comment: Accepted by MNRAS. First application of the methods described in arXiv:1804.05576

We present the first results of applying Gaussian Mixture Models in the stellar kinematic space of normalized angular momentum and binding energy on NIHAO high resolution galaxies to separate the stars into multiple components. We exemplify this method using a simulated Milky Way analogue, whose stellar component hosts: thin and thick discs, classical and pseudo bulges, and a stellar halo. The properties of these stellar structures are in good agreement with observational expectations in terms of sizes, shapes and rotational support. Interestingly, the two kinematic discs show surface mass density profiles more centrally concentrated than exponentials, while the bulges and the stellar halo are purely exponential. We trace back in time the Lagrangian mass of each component separately to study their formation history. Between z~3 and the end of halo virialization, z~1.3, all components lose a fraction of their angular momentum. The classical bulge loses the most (~95%) and the thin disc the least (~60%). Both bulges formed their stars in-situ at high redshift, while the thin disc formed ~98% in-situ, but with a constant SFR~1.5M$_{\rm\odot}$/yr$^{\rm-1}$ over the last ~11 Gyr. Accreted stars (6% of total stellar mass) are mainly incorporated to the thick disc or the stellar halo, which formed ex-situ 8% and 45% of their respective masses. Our analysis pipeline is freely available at https://github.com/aobr/gsf., Comment: MNRAS accepted. First results described in arXiv:1804.06635

We study the impact of the host on dwarf galaxy properties using four new Milky Way-like, ultra high-resolution simulations, ($N_{\rm part} >10^7$) from the NIHAO project. We split our sample in satellite ($R2.5 R_{\rm 200}$) galaxies. Simulated galaxies from all three groups are in excellent agreement with Local Group dwarf galaxies in terms of: stellar mass-velocity dispersion, stellar mass-metallicity relation, star formation histories, and stellar mass functions. Satellites and nearby galaxies show lower velocity dispersions and gas fractions compared to field galaxies. While field galaxies follow global abundance matching relations, satellites and nearby galaxies deviate from them, showing lower dark matter masses for given stellar mass. The reason for this deficit in dark matter mass is substantial mass loss experienced by satellites and $\sim80$\% of the nearby galaxies, while orbiting inside $R_{200}$ at earlier times. However, both satellites and nearby objects fall back onto the relation for field galaxies if we use the maximum of their virial mass instead of the present-day value. This allows us to provide estimates for the peak masses of observed Local Group galaxies. Finally, using radial velocities, distances, and the velocity dispersion-stellar mass relation from our simulations, we derive a metric to distinguish between galaxies harassed by the central object and unaffected ones. Applying this metric to observed objects we find that even far away dwarf galaxies like Eri II ($D\approx$ 370 kpc) have a strong probability ($\approx 83$\%) of having been affected by the Milky Way in the past. This naturally explains the lack of gas and recent star formation seen in Eri II., Comment: accepted by MNRAS, 31 pages, 17 pages main text, 14 pages appendix with tables of dwarf galaxy properties and convergence tests, 17 main figures

The GERDA experiment searches for the lepton number violating neutrinoless double beta decay of $^{76}$Ge ($^{76}$Ge $\rightarrow$ $^{76}$Se + 2e$^-$) operating bare Ge diodes with an enriched $^{76}$Ge fraction in liquid argon. The exposure for BEGe-type detectors is increased threefold with respect to our previous data release. The BEGe detectors feature an excellent background suppression from the analysis of the time profile of the detector signals. In the analysis window a background level of $1.0_{-0.4}^{+0.6}\cdot10^{-3}$ cts/(keV$\cdot$kg$\cdot$yr) has been achieved; if normalized to the energy resolution this is the lowest ever achieved in any 0$\nu\beta\beta$ experiment. No signal is observed and a new 90 \% C.L. lower limit for the half-life of $8.0\cdot10^{25}$ yr is placed when combining with our previous data. The median expected sensitivity assuming no signal is $5.8\cdot10^{25}$ yr., Comment: 5 pages, 2 figures