We study the global and radial stellar mass assembly of eight zoomed-in MW-sized galaxies produced in Hydrodynamics cosmological simulations. The disk-dominated galaxies (4) show a fast initial stellar mass growth in the innermost parts, driven mostly by in-situ SF, but since $z\sim2-1$ the SF enters in a long-term quenching phase. The outer regions follow this trend but more gentle as more external they are. As the result, the radial stellar mass growth is highly inside-out due to both the inside-out structural growth and inside-out SF quenching. The half-mass radius evolves fast; for instance, $R_{0.5}$($z=1$)$<0.5$$R_{0.5}$($z=0$). Two other runs resemble lenticular galaxies. One shows also a pronounced inside-out growth and the other one presents a nearly uniform radial mass assembly. The other two galaxies suffered late major mergers. Their normalized radial mass growth histories (MGHs) are nearly close among them but with periods of outside-in assembly during or after the mergers. For all the simulations, the archaeological radial MGHs calculated from the $z = 0$ stellar-particles age distribution are similar to the current MGHs, which evidences that the mass assembly by ex-situ stars and the radial mass transport do not change significantly their radial mass distributions. Our results agree qualitatively with observational inferences from the fossil record method applied to a survey of local galaxies and from look-back observations of progenitors of MW-sized galaxies. However, the inside-out growth mode is more pronounced and the $R_{0.5}$ growth is faster in simulations than in observational inferences.

We introduce a new set of simulations of Milky Way-sized galaxies using the AMR code ART + hydrodynamics in a $\Lambda$CDM cosmogony. The simulation series is named GARROTXA and follow the formation of a halo/galaxy from z~$=$~60 to z~$=$~0. The final virial mass of the system is $\sim$7.4$\times$10$^{11}$M$_{\odot}$. Our results are as follows: (a) contrary to many previous studies, the circular velocity curve shows no central peak and overall agrees with recent MW observations. (b) Other quantities, such as M$\_{*}$(6$\times$10$^{10}$M$_\odot$) and R$_d$ (2.56 kpc), fall well inside the observational MW range. (c) We measure the disk-to-total ratio kinematically and find that D/T=0.42. (d) The cold gas fraction and star formation rate (SFR) at z=0, on the other hand, fall short from the values estimated for the Milky Way. As a first scientific exploitation of the simulation series, we study the spatial distribution of the hot X-ray luminous gas. We have found that most of this X-ray emitting gas is in a halo-like distribution accounting for an important fraction but not all of the missing baryons. An important amount of hot gas is also present in filaments. In all our models there is not a massive disk-like hot gas distribution dominating the column density. Our analysis of hot gas mock observations reveals that the homogeneity assumption leads to an overestimation of the total mass by factors 3 to 5 or to an underestimation by factors $0.7-0.1$, depending on the used observational method. Finally, we confirm a clear correlation between the total hot gas mass and the dark matter halo mass of galactic systems.

Motivated by recent numerical simulations of molecular cloud (MC)evolution, in which the clouds engage in global gravitational contraction, and local collapse events culminate significantly earlier than the global collapse, we investigate the growth of density perturbations embedded in a collapsing background, to which we refer as an Inverse Hubble Flow (IHF). We use the standard procedure for the growth of perturbations in a universe that first expands (the usual Hubble Flow) and then recollapses (the IHF). We find that linear density perturbations immersed in an IHF grow faster than perturbations evolving in a static background (the standard Jeans analysis). A fundamental distinction between the two regimes is that, in the Jeans case, the time $\tau_\mathrm{nl}$ for a density fluctuation to become nonlinear increases without limit as its initial value approaches zero, while in the IHF case $\tau_\mathrm{nl} \le \tau_\mathrm{ff}$ always, where $\tau_\mathrm{ff}$ is the free-fall time of the background density. We suggest that this effect, although moderate, implies that small-scale density fluctuations embedded in globally-collapsing clouds must collapse earlier than their parent cloud, regardless of whether the initial amplitude of the fluctuations is moderate or strongly nonlinear, thus allowing the classical mechanism of Hoyle fragmentation to operate in multi-Jeans-mass MCs. More fundamentally, our results show that, contrary to the standard paradigm that fluctuations of all scales grow at the same rate in the linear regime, the hierarchical nesting of the fluctuations of different scales does affect their growth even in the linear stage.

Supermassive black hole (BH) binaries would comprise the strongest sources of gravitational waves (GW) once they reach <<1 pc separations, for both pulsar timing arrays (PTAs) and space based (SB) detectors. While BH binaries coalescences constitute a natural outcome of the cosmological standard model and galaxy mergers, their dynamical evolution is still poorly understood and therefore their abundances at different stages. We use a dynamical model for the decay of BH binaries coupled with a cosmological simulation and semi-empirical approaches to the occupation of haloes by galaxies and BHs, in order to follow the evolution of the properties distribution of galaxies hosting BH binaries candidates to decay due to GWs emission. Our models allow us to relax simplifying hypothesis about the binaries occupation in galaxies and their mass, as well as redshift evolution. Following previously proposed electromagnetic (EM) signatures of binaries in the subpc regime, that include spectral features and variability, we model possible distributions of such signatures and also set upper limits to their lifespan. We found a bimodal distribution of hosts properties, corresponding to BH binaries suitable to be detected by PTA and the ones detectable only from space missions, as eLISA. Although it has been discussed that the peak of eLISA sources may happen at high z, we show that there must be a population of such sources in the nearby Universe that might show detectable EM signatures, representing an important laboratory for multimessenger astrophysics. We found a weak dependence of galaxy host properties on the binaries occupation, that can be traced back to the BH origin. The combination of the host correlations reported here with the expected EM signal, may be helpful to verify the presence of nearby GW candidates, and to distinguish them from 'regular' intrinsic AGN variability.

We present high-resolution N-body/hydrodynamics simulations of dwarf galaxies formed in isolated CDM halos with the same virial mass, Mv~2.5x10^10 Msun at z=0, in order to (1) study the mass assembly histories (MAHs) of the halo, stars, and gas components, and (2) explore the effects of the halo MAHs on the stellar/baryonic assembly of the simulated dwarfs and on their z~0 properties. Overall, the simulated dwarfs are roughly consistent with observations. Our main results are: a) The stellar-to-halo mass ratio is ~0.01 and remains roughly constant since z~1 (the stellar MAHs follow closely the halo MAHs), with a smaller value at higher z's for those halos that assemble their mass later. b) The evolution of the galaxy gas fraction, fg, is episodic and higher, most of the time, than the stellar fraction. When fg decreases (increases), the gas fraction in the halo typically increases (decreases), showing that the SN driven outflows play an important role in regulating the gas fractions -and hence the SFR- of the dwarfs. However, in most cases, an important fraction of the gas escapes the virial radius, Rv; at z=4 the total baryon fraction inside Rv is 1.5-2 times smaller than the universal one, while at z=0 is 2-6 times smaller, with the earlier assembled halos ejecting more gas. c) The SF histories are episodic with changes in the SFRs of factors 2-10 on average. d) Although the dwarfs formed in late assembled halos show more extended SF histories, their z~0 SFRs are still below the ones measured for local isolated dwarfs. e) The effects of baryons on Mv are such that at almost any time Mv is 10-20% smaller than the corresponding Mv obtained in pure N-body simulations. Our results suggest that rather than increasing the strength of the SN-driven outflows, processes that reduce the SF efficiency even more will help to solve the potential issues faced by the CDM-based simulations of dwarfs.

We present an idealized, semi-empirical model for the evolution of gravitationally contracting molecular clouds (MCs) and their star formation rate (SFR) and efficiency (SFE). The model assumes that the instantaneous SFR is given by the mass above a certain density threshold divided by its free-fall time. The instantaneous number of massive stars is computed assuming a Kroupa IMF. These stars feed back on the cloud through ionizing radiation, eroding it. The main controlling parameter of the evolution turns out to be the maximum cloud mass, $\Mmax$. This allows us to compare various properties of the model clouds against their observational counterparts. A giant molecular cloud (GMC) model ($\Mmax \sim 10^5 \Msun$) adheres very well to the evolutionary scenario recently inferred by Kawamura et al. (2009) for GMCs in the Large Magellanic Cloud. A model cloud with $\Mmax \approx 2000 \Msun$ evolves in the Kennicutt-Schmidt diagram first passing through the locus of typical low- to-intermediate mass star-forming clouds, and then moving towards the locus of high-mass star-forming ones over the course of $\sim 10$ Myr. Also, the stellar age histograms for this cloud a few Myr before its destruction agree very well with those observed in the $\rho$-Oph stellar association, whose parent cloud has a similar mass, and imply that the SFR of the clouds increases with time. Our model thus agrees well with various observed properties of star-forming MCs, suggesting that the scenario of gravitationally collapsing MCs, with their SFR regulated by stellar feedback, is entirely feasible and in agreement with key observed properties of molecular clouds.

We study the effect that the dark matter background (DMB) has on the gravitational energy content and, in general, on the star formation efficiency of a molecular cloud (MC). We first analyze the effect that a dark matter halo, described by the Navarro et al. (1996) density profile, has on the energy budget of a spherical, homogeneous, cloud located at different distances from the halo center. We found that MCs located in the innermost regions of a massive galaxy can feel a contraction force greater than their self-gravity due to the incorporation of the potential of the galaxy's dark matter halo. We also calculated analytically the gravitational perturbation that a MC produces over a uniform DMB (uniform at the scales of a MC) and how this perturbation will affect the evolution of the MC itself. The study shows that the star formation in a MC will be considerably enhanced if the cloud is located in a dense and low velocity dark matter environment. We confirm our results by measuring the star formation efficiency in numerical simulations of the formation and evolution of MCs within different DMBs. Our study indicates that there are situations where the dark matter's gravitational contribution to the evolution of the molecular clouds should not be neglected.

(Abridged) By means of high-resolution cosmological simulations in the context of the LCDM scenario, the specific star formation rate (SSFR=SFR/Ms, Ms is the stellar mass)--Ms and stellar mass fraction (Fs=Ms/Mh, Mh is the halo mass)--Ms relations of low-mass galaxies (2.5< Mh/10^10 Msun <50 at redshift z=0) at different epochs are predicted. The Hydrodynamics ART code was used and some variations of the sub-grid parameters were explored. Most of simulated galaxies, specially those with the highest resolutions, have significant disk components and their structural and dynamical properties are in reasonable agreement with observations of sub-M* field galaxies. However, the SSFRs are 5-10 times smaller than the averages of several (compiled and homogenized here) observational determinations for field blue/star-forming galaxies at z<0.3 (at low masses, most of observed field galaxies are actually blue/star-forming). This inconsistency seems to remain even at z~1.5 though less drastic. The Fs of simulated galaxies increases with Mh as semi-empirical inferences show, but in absolute values the former are ~5-10 times larger than the latter at z=0; this difference increases probably to larger factors at z~1-1.5. The inconsistencies reported here imply that simulated low-mass galaxies (0.2<Ms/10^9 Msun <30 at z=0) assembled their stellar masses much earlier than observations suggest. This confirms the predictions previously found by means of LCDM-based models of disk galaxy formation and evolution for isolated low-mass galaxies (Firmani & Avila-Reese 2010), and highlight that our implementation of astrophysics into simulations and models are still lacking vital ingredients.