LCDM cosmological models with Early Dark Energy (EDE) have been proposed to resolve tensions between the Hubble constant H0 = 100h km/s/Mpc measured locally, giving h ~ 0.73, and H0 deduced from Planck cosmic microwave background (CMB) and other early universe measurements plus LCDM, giving h ~ 0.67. EDE models do this by adding a scalar field that temporarily adds dark energy equal to about 10% of the cosmological energy density at the end of the radiation-dominated era at redshift z ~ 3500. Here we compare linear and nonlinear predictions of a Planck-normalized LCDM model including EDE giving h = 0.728 with those of standard Planck-normalized LCDM with h = 0.678. We find that nonlinear evolution reduces the differences between power spectra of fluctuations at low redshifts. As a result, at z = 0 the halo mass functions on galactic scales are nearly the same, with differences only 1-2%. However, the differences dramatically increase at high redshifts. The EDE model predicts 50% more massive clusters at z = 1 and twice more galaxy-mass halos at z = 4. Even greater increases in abundances of galaxy-mass halos at higher redshifts may make it easier to reionize the universe with EDE. Predicted galaxy abundances and clustering will soon be tested by JWST observations. Positions of baryonic acoustic oscillations (BAOs) and correlation functions differ by about 2% between the models -- an effect that is not washed out by nonlinearities. Both standard LCDM and the EDE model studied here agree well with presently available acoustic-scale observations, but DESI and Euclid measurements will provide stringent new tests.

We report the bivariate HI- and H$_2$-stellar mass distributions of local galaxies in addition of an inventory of galaxy mass functions, MFs, for HI, H$_2$, cold gas, and baryonic mass, separately into early- and late-type galaxies. The MFs are determined using the HI and H$_2$ conditional distributions and the galaxy stellar mass function, GSMF. For the conditional distributions we use the compilation presented in Calette et al. 2018. For determining the GSMF from $M_{\ast}\sim3\times10^{7}$ to $3\times10^{12}$ $M_{\odot}$, we combine two spectroscopic samples from the SDSS at the redshift range $0.0033<z<0.2$. We find that the low-mass end slope of the GSMF, after correcting from surface brightness incompleteness, is $\alpha\approx-1.4$, consistent with previous determinations. The obtained HI MFs agree with radio blind surveys. Similarly, the H$_2$ MFs are consistent with CO follow-up optically-selected samples. We estimate the impact of systematics due to mass-to-light ratios and find that our MFs are robust against systematic errors. We deconvolve our MFs from random errors to obtain the intrinsic MFs. Using the MFs, we calculate cosmic density parameters of all the baryonic components. Baryons locked inside galaxies represent 5.4% of the universal baryon content, while $\sim96$% of the HI and H$_2$ mass inside galaxies reside in late-type morphologies. Our results imply cosmic depletion times of H$_2$ and total neutral H in late-type galaxies of $\sim 1.3$ and 7.2 Gyr, respectively, which shows that late type galaxies are on average inefficient in converting H$_2$ into stars and in transforming HI gas into H$_2$. Our results provide a fully self-consistent empirical description of galaxy demographics in terms of the bivariate gas--stellar mass distribution and their projections, the MFs. This description is ideal to compare and/or to constrain galaxy formation models.

We use the stellar kinematics for $2458$ galaxies from the MaNGA survey to explore dynamical scaling relations between the stellar mass $M_{\star}$ and the total velocity parameter at the effective radius, $R_e$, defined as $S_{K}^{2}=KV_{R_e}^{2}+\sigma_{\star_e}^{2}$, which combines rotation velocity $V_{R_e}$, and velocity dispersion $\sigma_{\star_e}$. We confirm that spheroidal and spiral galaxies follow the same $M_{\star}-S_{0.5}$ scaling relation with lower scatter than the $M_{\star}-V_{R_e}$ and $M_{\star}-\sigma_{\star_e}$ ones. We also explore a more general Universal Fundamental Plane described by the equation $log(\Upsilon_{e}) = log (S_{0.5}^{2}) - log (I_{e}) - log (R_{e}) + C$, which in addition to kinematics, $S_{0.5}$, and effective radius, $R_e$, it includes surface brightness, $I_e$, and dynamical mass-to-light ratio, $\Upsilon_e$. We use sophisticated Schwarzschild dynamical models for a sub-sample of 300 galaxies from the CALIFA survey to calibrate the so called Universal Fundamental Plane. That calibration allows us to propose both: (i) a parametrization to estimate the difficult-to-measure dynamical mass-to-light ratio at the effective radius; and (ii) a new dynamical mass proxy consistent with dynamical models within $0.09\ dex$. We reproduce the relation between the dynamical mass and the stellar mass in the inner regions of galaxies. We use the estimated dynamical mass-to-light ratio from our analysis, $\Upsilon_{e}^{fit}$, to explore the Universal Fundamental Plane with the MaNGA data set. We find that all classes of galaxies, from spheroids to disks, follow this Universal Fundamental Plane with a scatter significantly smaller $(0.05\ dex)$ than the one reported for the $M_{\star}-S_{0.5}$ relation $(0.1\ dex)$, the Fundamental Plane $(\sim 0.09\ dex)$ and comparable with Tully-Fisher studies $(\sim 0.05\ dex)$, but for a wider range of galaxy types.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

We perform a "fossil record" analysis for ~800 low-redshift spiral galaxies, using STARLIGHT applied to integral field spectroscopic observations from the SDSS-IV MaNGA survey to obtain fully spatially-resolved high-resolution star formation histories (SFHs). From the SFHs, we are able to build maps indicating the present-day distribution of stellar populations of different ages in each galaxy. We find small negative mean age gradients in most spiral galaxies, especially at high stellar mass, which reflects the formation times of stellar populations at different galactocentric radii. We show that the youngest (<10^{8.5} years) populations exhibit significantly more extended distributions than the oldest (>10^{9.5} years), again with a strong dependence on stellar mass. By interpreting the radial profiles of "time slices" as indicative of the size of the galaxy at the time those populations had formed, we are able to trace the simultaneous growth in mass and size of the spiral galaxies over the last 10 Gyr. Despite finding that the evolution of the measured light-weighted radius is consistent with inside-out growth in the majority of spiral galaxies, the evolution of an equivalent mass-weighted radius has changed little over the same time period. Since radial migration effects are likely to be small, we conclude that the growth of disks in spiral galaxies has occurred predominantly through an inside-out mode (with the effect greatest in high-mass galaxies), but this has not had anywhere near as much impact on the distribution of mass within spiral galaxies.

We study the spatially resolved properties of 333 elliptical galaxies with the MaNGA/SDSS-IV survey. The aim is to understand the fundamental processes of formation and quenching of elliptical galaxies. We used the DESI Legacy Imaging Surveys for accurate morphological classification. Based on integrated spectroscopic properties and colors, we classified galaxies into classical "red and dead", recently quenched, and blue star-forming ellipticals (CLEs, RQEs, and BSFs corresponding to 75%, 10%, and 4% of the sample, respectively). We inferred their stellar age and stellar metallicity gradients out to 1.5 effective radius, and reconstructed their global and radial histories of mass growth and star formation. We find the mass- and light-weighted age gradients of CLEs are nearly flat or mildly negative, with small differences between both ages. The respective metallicity gradients are negative, being flatter as less massive are the CLEs. The more massive CLEs assembled earlier and quenched faster than the less massive ones. The CLEs show a weak inside-out growth and a clear inside-out quenching. At masses < $10^{11}$ $M_{\bigodot}$, the age and $Z$ gradients of the RQEs and BSFs are flatter than those of the CLEs but with larger scatters. They show very weak inside-out growth and quenching, being the quenching slow and even not completed at $z\sim$ 0 for the BSFs. Instead, the massive RQEs show an outside-in quenching, and positive gradients in the light-weighted age and stellar metallicities. Our results are consistent with a scenario where the inner parts of CLEs formed by an early and coeval dissipative collapse with a consequent burst of star formation and further quenching, whereas the outer parts continued their assembly likely by dry mergers. We also discuss some evolutionary scenarios for the RQE and BSF galaxies that would agree with their generic results presented here.

We present the HI-MaNGA programme of HI follow-up for the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. MaNGA, which is part of the Fourth phase of the Sloan Digital Sky Surveys (SDSS-IV), is in the process of obtaining integral field unit (IFU) spectroscopy for a sample of ~10,000 nearby galaxies. We give an overview of the HI 21cm radio follow-up observing plans and progress and present data for the first 331 galaxies observed in the 2016 observing season at the Robert C. Bryd Green Bank Telescope (GBT). We also provide a cross match of the current MaNGA (DR15) sample with publicly available HI data from the Arecibo Legacy Fast Arecibo L-band Feed Array (ALFALFA) survey. The addition of HI data to the MaNGA data set will strengthen the survey's ability to address several of its key science goals that relate to the gas content of galaxies, while also increasing the legacy of this survey for all extragalactic science.

We investigate the dependence of the galaxy-halo connection and galaxy density field in modified gravity models using the $N-$body simulations for $f(R)$ and nDGP models at $z=0$. Because of the screening mechanisms employed by these models, chameleon and Vainshtein, halos are clustered differently in the non-linear regime of structure formation. We quantify their deviations in the galaxy density field from the standard $\Lambda$CDM model under different environments. We populate galaxies in halos via the (Sub)Halo Abundance Matching. Our main results are: 1) The galaxy-halo connection {\it strongly} depends on the gravity model; a maximum variation of $\sim40\%$ is observed between Halo Occupational Distribution (HOD) parameters; 2) $f(R)$ gravity models predict an excess of galaxies in low density environments of $\sim10\%$ but predict a deficit of $\sim10\%$ at high density environments for $|f_{R0}| = 10^{-4}$ and $10^{-6}$ while $|f_{R0}| = 10^{-5}$ predicts more high density structures; nDGP models are consistent with $\Lambda$CDM; 3) Different gravity models predict different dependences of the galaxy luminosity function (GLF) with the environment, especially in void-like regions we find differences around $\sim10\%$ for the $f(R)$ models while nDPG models remain closer to $\Lambda$CDM for low-luminosity galaxies but there is a deficit of $\sim11\%$ for high-luminosity galaxies in all environments. We conclude that the dependence of the GLF with environment might provide a test to distinguish between gravity models and their screening mechanisms from the $\Lambda$CDM. We provide HOD parameters for the gravity models analyzed in this paper.

Crucial progress in our understanding of star formation and feedback will depend on the ability to obtain spatially resolved spectroscopic observations of \ion{H}{ii} regions, from which reliable instantaneous measurements of their physical conditions can be obtained. Acquiring these datasets across full galactic systems will prove crucial for obtaining population samples that enable us to understand the time evolution of similar regions, and the variability of conditions among coeval regions. Separating the spatial and temporal dependencies in such way for different physical processes involved in star formation and the injection of feedback is crucial to overcome the inherit degeneracies associated with observing instantaneous snapshots of a dynamic ISM at any given time. Emission line diagnostics are at the core of measuring the physical condition in \ion{H}{ii} regions (e.g. dynamics, SFR, chemical abundances, dust extinction, ionization and excitation, etc.). These measurements require high spatial resolution, contiguous coverage across full galactic systems, and sensitivities significantly deeper than past efforts. The spatial scale required to resolve the \ion{H}{ii} regions of a few pc is only attainable in the Local Group where very large sky coverage is necessary.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

By means of post-processed cosmological Hydrodynamics simulations we explore the ability of the fossil record method to recover the stellar mass, age gradients, and global/radial star formation and mass growth histories of galaxies observed with an optical Integral Field Unity. We use two simulations of representative Milky Way-sized galaxies, one disk dominated and another bulge-dominated. We generate data cubes emulating MaNGA/SDSS-IV instrumental and observational conditions, and analyze them using the Pipe3D pipeline. For optimal MaNGA instrumental/observational setups, the global masses, ages, and dust extinctions of both galaxies are well recovered within the uncertainties, while for the radial distributions, the trend is to get slightly flatter age gradients and weaker inside-out growth modes than the true ones. The main bias is in recovering younger stellar populations in the inner, older regions, and slightly older stellar populations in the outer, younger regions, as well as lower global stellar masses. These trends are enhanced as lower is the spatial sampling/resolution (less number of fibers), more inclined and dust-attenuated is the observed galaxy, lower is the signal-to-noise ratio, and worse is the seeing. The recovered properties and histories are more affected for the disk-dominated galaxy with a prominent inside-out growth mode than for the bulge-dominated galaxy assembled earlier and more coherently. We make publicly available the code and the mock IFU datacubes used in this paper.

We present the results of the archaeological analysis of the stellar populations of a sample of ~4,000 galaxies observed by the SDSS-IV-MaNGA survey using Pipe3D. Based on this analysis we extract a sample of ~150,000 SFRs and stellar masses that mimic a single cosmological survey covering the redshift range between z~0 to z~7. We confirm that the Star-Forming Main Sequence holds as a tight relation in this range of redshifts, and evolves strongly in both the zero-point and slope. This evolution is different for the population of local star-forming (SFGs) and retired (RGs) galaxies, with the latter presenting a stronger evolution in the zero-point and a weaker evolution in the slope. The fraction of RGs decreases rapidly with z, particularly for those classified as RGs at z~0. Contrary to previous studies we detect RGs well above z>1, although not all of them are progenitors of local RGs. Finally, adopting the required corrections to make the survey complete in mass in a limited volume, we recover the cosmic star-formation rate (SFR), stellar mass density, and average specific SFR histories of the Universe in this wide range of look-back times, with a remarkable agreement with the values reported by various cosmological surveys. We demonstrate that the progenitors of local RGs were more actively forming stars in the past, contributing to most of the cosmic SFR density at z>0.5, and to most of the cosmic stellar mass density at any redshift. They suffer a general quenching in the SFR at z~0.35. Below this redshift the progenitors of local SFGs dominate the SFR density of the Universe.

We compile and carrefully homogenize local galaxy samples with available information on stellar, $\rm HI$ and/or $\rm H_{2}$ masses, and morphology. After processing the information on upper limits in the case of non gas detections, we determine the $\rm HI$- and $\rm H_{2}$-to-stellar mass relations and their $1\sigma$ scatter for both late- and early-type galaxies. The obtained relations are fitted to single or double power laws. Late-type galaxies are significantly gas richer than early-type ones, specially at high masses. The respective $\rm H_{2}$-to-$\rm HI$ mass ratios as a function of $M_{\ast}$ are discussed. Further, we constrain the full mass-dependent distribution functions of the $\rm HI$- and $\rm H_{2}$-to-stellar mass ratios. We find that they can be described by a Schechter function for late types and a (broken) Schechter + uniform function for early types. By using the observed galaxy stellar mass function and the volume-complete late-to-early-type galaxy ratio as a function of $M_{\ast}$, these empirical distribution functions are mapped into $\rm HI$ and $\rm H_{2}$ mass functions. The obtained mass functions are consistent with those inferred from large surveys. The empirical gas-to-stellar mass relations and their distributions for local late- and early-type galaxies presented here can be used to constrain models and simulations of galaxy evolution.

Understanding the formation and evolution of early-type, spheroid-dominated galaxies is an open question within the context of the hierarchical clustering scenario, particularly, in low-density environments. Our goal is to study the main structural, dynamical, and stellar population properties and assembly histories of field spheroid-dominated galaxies formed in a LCDM scenario to assess to what extend they are consistent with observations. We selected spheroid-dominated systems from a LCDM simulation that includes star formation, chemical evolution and Supernova feedback. A sample of 18 field systems with Mstar <= 6x10^10 Msun that are dominated by the spheroid component. For this sample we estimate the fundamental relations of ellipticals and then compared with current observations. The simulated spheroid galaxies have sizes in good agreement with observations. The bulges follow a Sersic law with Sersic indexes that correlate with the bulge-to-total mass ratios. The structural-dynamical properties of the simulated galaxies are consistent with observed Faber-Jackson, Fundamental Plane, and Tully-Fisher relations. However, the simulated galaxies are bluer and with higher star formation rates than observed isolated early-type galaxies. The archaeological mass growth histories show a slightly delayed formation and more prominent inside-out growth mode than observational inferences based on the fossil record method. The main structural and dynamical properties of the simulated spheroid-dominated galaxies are consistent with observations. This is remarkable since none of them has been tuned to be reproduced. However, the simulated galaxies are blue and star-forming, and with later stellar mass growth histories as compared to observational inferences. This is mainly due to the persistence of extended discs in the simulations. Abridged

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.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

We present here the characterization of the main properties of a sample of 98 AGN host galaxies, both type-II and type-I, in comparison with those of about 2700 non-active galaxies observed by the MaNGA survey. We found that AGN hosts are morphologically early-type or early-spirals. For a given morphology AGN hosts are, in average, more massive, more compact, more central peaked and rather pressurethan rotational-supported systems. We confirm previous results indicating that AGN hosts are located in the intermediate/transition region between star-forming and non-star-forming galaxies (i.e., the so-called green valley), both in the ColorMagnitude and the star formation main sequence diagrams. Taking into account their relative distribution in terms of the stellar metallicity and oxygen gas abundance and a rough estimation of their molecular gas content, we consider that these galaxies are in the process of halting/quenching the star formation, in an actual transition between both groups. The analysis of the radial distributions of the starformation rate, specific star-formation rate, and molecular gas density shows that the quenching happens from inside-out involving both a decrease of the efficiency of the star formation and a deficit of molecular gas. All the intermediate data-products used to derive the results of our analysis are distributed in a database including the spatial distribution and average properties of the stellar populations and ionized gas, published as a Sloan Digital Sky Survey Value Added Catalog being part of the 14th Data Release: http://www.sdss.org/dr14/manga/manga-data/manga-pipe3d-value-added-catalog/

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.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

We present new determinations of the stellar-to-halo mass relation (SHMR) at $z=0-10$ that match the evolution of the galaxy stellar mass function, the SFR$-M_*$ relation,and the cosmic star formation rate. We utilize a compilation of 40 observational studies from the literature and correct them for potential biases. Using our robust determinations of halo mass assembly and the SHMR, we infer star formation histories, merger rates, and structural properties for average galaxies, combining star-forming and quenched galaxies. Our main findings: (1) The halo mass $M_{50}$ above which 50\% of galaxies are quenched coincides with sSFR/sMAR$\sim1$, where sMAR is the specific halo mass accretion rate. (2) $M_{50}$ increases with redshift, presumably due to cold streams being more efficient at high redshift while virial shocks and AGN feedback become more relevant at lower redshifts. (3) The ratio sSFR/sMAR has a peak value, which occurs around $M_{\rm vir}\sim2\times10^{11}M_{\odot}$. (4) The stellar mass density within 1 kpc, $\Sigma_1$, is a good indicator of the galactic global sSFR. (5) Galaxies are statistically quenched after they reach a maximum in $\Sigma_1$, consistent with theoretical expectations of the gas compaction model; this maximum depends on redshift. (6) In-situ star formation is responsible for most galactic stellar mass growth, especially for lower-mass galaxies. (7) Galaxies grow inside out. The marked change in the slope of the size--mass relation when galaxies became quenched, from $d\log R_{\rm eff}/d\log M_*\sim0.35$ to $\sim2.5$, could be the result of dry minor mergers.

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 dark matter spike induced by the adiabatic growth of a massive black hole in a cuspy environment, may explain the thermal dark matter density required to fit the cut-off in the HESSJ1745-290 gamma-ray spectra as TeV dark matter signal with a background component. The spike extension appears comparable with the HESS angular resolution.

We analyze the possibility that the HESS gamma-ray source at the Galactic Center could be explained as the secondary flux produced by annihilation of TeV Dark Matter (TeVDM) particles with locally enhanced density, in a region spatially compatible with the HESS observations themselves. We study the inner 100 pc considering (i) the extrapolation of several density profiles from state-of-the-art N-body + Hydrodynamics simulations of Milky Way-like galaxies, (ii) the DM spike induced by the black hole, and (iii) the DM particles scattering off by bulge stars. We show that in some cases the DM spike may provide the enhancement in the flux required to explain the cut-off in the HESS J1745-290 gamma-ray spectra as TeVDM. In other cases, it may helps to describe the spatial tail reported by HESS II at angular scales < 0.54 degrees towards Sgr A.

By means of the fossil record method implemented through Pipe3D, we reconstruct the global and radial stellar mass growth histories (MGHs) of an unprecedentedly large sample of galaxies, ranging from dwarf to giant objects, from the "Mapping Nearby Galaxies at the Apache Point Observatory" survey. We confirm that the main driver of the global MGHs is mass, with more massive galaxies assembling their masses earlier (downsizing), though for a given mass, the global MGHs segregate by color, specific star formation rate (sSFR), and morphological type. From the inferred radial mean MGHs, we find that at the late evolutionary stages (or for fractions of assembled mass larger than ~ 80%), the innermost regions formed stars on average earlier than the outermost ones (inside-out). At earlier epochs, when the age resolution of the method becomes poor, the mass assembly seems to be spatially homogeneous or even in the outside-in mode, specially for the red/quiescent/early-type galaxies. The innermost MGHs are in general more regular (less scatter around the mean) than the outermost ones. For dwarf and low-mass galaxies, we do not find evidence of an outside-in formation mode; instead their radial MGHs are very diverse most of the time, with periods of outside- in and inside-out modes (or strong radial migration), suggesting this an episodic SF history. Blue/star-forming/late-type galaxies present on average a significantly more pronounced inside-out formation mode than red/quiescent/early-type galaxies, independently of mass. We discuss our results in the light of the processes of galaxy formation, quenching, and radial migration. We discuss also on the uncertainties and biases of the fossil record method and how they could affect our results.

We introduce a new set of eight Milky Way-sized cosmological simulations performed using the AMR code ART + Hydrodynamics in a LCDM cosmology. The set of zoom-in simulations covers present-day virial masses in the 0.83-1.56 x 10^12 msun range and is carried out with our simple but effective deterministic star formation (SF) and ``explosive' stellar feedback prescriptions. The work is focused on showing the goodness of the simulated set of ``field' Milky Way-sized galaxies. Our results are as follows. (a) The circular velocity curves of our simulated galaxies are nearly flat. (b) Runs ending with a significant disk component, for their stellar masses, have V_max, radius, SF rate, gas fraction, and specific angular momentum values consistent with observations of late-type galaxies. (C) The two most spheroid-dominated galaxies formed in halos with late active merger histories, but other run that ends also as spheroid-dominated, never had major mergers. (d) Our simulations are consistent with the empirical stellar-to-halo mass correlation, and those that end as disk-dominated, evolve mostly along the low-mass branch of this correlation. (e) Moreover, since the last 6.5-10 Gyr, the baryonic/stellar and halo mass growth histories are proportional. (f) Within Rvir ~ 25-50% of the baryons are missed. (g) The z ~ 0 gas velocity dispersion profiles, sigma_z(r), are nearly flat and can be mostly explained by the kinetic energy injected by stars. (h) The average values of sigma_z increase at higher redshifts, following roughly the shape of the SF history.

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.

By means of N-body+Hydrodynamic zoom-in simulations we study the evolution of the inner dark matter and stellar mass distributions of central dwarf galaxies formed in halos of virial masses Mv=2-3x10^10 Msun at z=0, both in a WDM and CDM cosmology. The half-mode mass in the WDM power spectrum of our simulations is Mf= 2x 10^10 Msun. In the dark matter only simulations halo density profiles are well described by the NFW parametric fit in both cosmologies, though the WDM halos have concentrations lower by factors 1.5--2.0 than their CDM counterparts. In the hydrodynamic simulations, the effects of baryons significantly flatten the inner density, velocity dispersion, and pseudo phase-space density profiles of the WDM halos but not of the CDM ones. The density slope measured at ~ 0.02xRv, alpha, becomes shallow in periods of 2 to 5 Gyr in the WDM runs. We explore whether this flattening process correlates with the global SF, Ms/Mv ratio, gas outflow, and internal specific angular momentum histories. We do not find any clear trends but when alpha is shallower than -0.5, Ms/Mv is always between 0.25 and 1%. We conclude that the main reason of the formation of the shallow core is the presence of strong gas mass fluctuations inside the inner halo, which are a consequence of the feedback driven by a very bursty and sustained SF history in shallow gravitational potentials. Our WDM halos, which assemble late and are less concentrated than the CDM ones, obey these conditions. There are also (rare) CDM systems with extended mass assembly histories that obey these conditions and form indeed shallow cores. The dynamical heating and expansion processes, behind the DM core flattening, apply also to the stars in a such a way that the stellar age and metallicity gradients of the dwarfs are softened, their stellar half-mass radii strongly grow with time, and their central surface densities decrease.

We have studied a sample of 89 very isolated, elliptical galaxies at z < 0.08 and compared their properties with elliptical galaxies located in a high-density environment such as the Coma supercluster. Our aim is to probe the role of environment on the morphological transformation and quenching of elliptical galaxies as a function of mass. In addition, we elucidate the nature of a particular set of blue and star-forming isolated ellipticals identified here. We study physical properties of ellipticals such as color, specific star formation rate, galaxy size, and stellar age, as a function of stellar mass and environment based on SDSS data. We analyze the blue star-forming isolated ellipticals in more detail, through photometric characterization using GALFIT, and infer their star formation history using STARLIGHT. Among the isolated ellipticals ~ 20% are blue, 8% are star forming, and ~ 10% are recently quenched, while among the Coma ellipticals ~ 8% are blue and just <= 1% are star forming or recently quenched. There are four isolated galaxies (~ 4.5%) that are blue and star forming at the same time. These galaxies, with masses between 7 x 10^9 and 2 x 10^10 h-2 M_sun, are also the youngest galaxies with light-weighted stellar ages <= 1 Gyr and exhibit bluer colors toward the galaxy center. Around 30-60% of their present-day luminosity, but only < 5% of their present-day mass, is due to star formation in the last 1 Gyr. The processes of morphological transformation and quenching seem to be in general independent of environment since most of elliptical galaxies are 'red and dead', although the transition to the red sequence should be faster for isolated ellipticals. In some cases, the isolated environment seems to propitiate the rejuvenation of ellipticals by recent (< 1 Gyr) cold gas accretion.

We present zoom-in N-body + Hydrodynamic simulations of dwarf central galaxies formed in Warm Dark Matter (WDM) halos with masses at present-day of $2-4\times 10^{10}$ \msun. Two different cases are considered, the first one when halo masses are close to the corresponding half-mode filtering scale \Mhm\ (\mwdm =1.2 keV), and the second when they are 20 to 30 times the corresponding \Mhm\ (\mwdm = 3.0 keV). The WDM simulations are compared with the respective Cold Dark Matter (CDM) simulations. The dwarfs formed in halos of masses (20-30)\Mhm have roughly similar properties and evolution than their CDM counterparts; on the contrary, those formed in halos of masses around \Mhm, are systematically different from their CDM counterparts. As compared to the CDM dwarfs, they assemble the dark and stellar masses later, having mass-weighted stellar ages 1.4--4.8 Gyr younger; their circular velocity profiles are shallower, with maximal velocities 20--60% lower; their stellar distributions are much less centrally concentrated and with larger effective radii, by factors 1.3--3. The WDM dwarfs at the filtering scale (\mwdm =1.2 keV) have disk-like structures, and end in most cases with higher gas fractions and lower stellar-to-total mass ratios than their CDM counterparts. The late halo assembly, low halo concentrations, and the absence of satellites of the former with respect to the latter, are at the basis of the differences.

We derive the stellar-to-halo mass relations, SHMR, of local blue and red central galaxies separately, as well as the fraction of halos hosting blue/red central galaxies. We find that: 1) the SHMR of central galaxies is segregated by color, with blue centrals having a SHMR above the one of red centrals; at logMh~12, the Ms/Mh ratio of the blue centrals is ~0.05, which is ~1.7 times larger than the value of red centrals. 2) The intrinsic scatters of the SHMRs of red and blue centrals are ~0.14 and ~0.11dex, respectively. The intrinsic scatter of the average SHMR of all central galaxies changes from ~0.20dex to ~0.14dex in the 11.3<logMh<15 range. 3) The fraction of halos hosting blue centrals at Mh=1E11Msun is 87%, but at 2x1E12Msun decays to ~20%, approaching to a few per cents at higher masses. The characteristic mass at which this fraction is the same for blue and red galaxies is Mh~7x1E11Msun. Our results suggest that the SHMR of central galaxies at large masses is shaped by halo mass quenching (likely through shock virial heating and AGN feedback), but group richness also plays an important role: central galaxies living in less dense environments quenched their star formation later or did not quench it yet. At low masses, processes that delay star formation without invoking too strong supernova-driven outflows could explain the high Ms/Mh ratios of blue centrals as compared to those of the scarce red centrals.

We present an overview of a new integral field spectroscopic survey called MaNGA (Mapping Nearby Galaxies at Apache Point Observatory), one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) that began on 2014 July 1. MaNGA will investigate the internal kinematic structure and composition of gas and stars in an unprecedented sample of 10,000 nearby galaxies. We summarize essential characteristics of the instrument and survey design in the context of MaNGA's key science goals and present prototype observations to demonstrate MaNGA's scientific potential. MaNGA employs dithered observations with 17 fiber-bundle integral field units that vary in diameter from 12" (19 fibers) to 32" (127 fibers). Two dual-channel spectrographs provide simultaneous wavelength coverage over 3600-10300 A at R~2000. With a typical integration time of 3 hr, MaNGA reaches a target r-band signal-to-noise ratio of 4-8 (per A, per 2" fiber) at 23 AB mag per sq. arcsec, which is typical for the outskirts of MaNGA galaxies. Targets are selected with stellar mass greater than 1e9 Msun using SDSS-I redshifts and i-band luminosity to achieve uniform radial coverage in terms of the effective radius, an approximately flat distribution in stellar mass, and a sample spanning a wide range of environments. Analysis of our prototype observations demonstrates MaNGA's ability to probe gas ionization, shed light on recent star formation and quenching, enable dynamical modeling, decompose constituent components, and map the composition of stellar populations. MaNGA's spatially resolved spectra will enable an unprecedented study of the astrophysics of nearby galaxies in the coming 6 yr.

We analyse the dark, gas, and stellar mass assembly histories of low-mass halos (Mvir ~ 10^10.3 - 10^12.3 M_sun) identified at redshift z = 0 in cosmological numerical simulations. Our results indicate that for halos in a given present-day mass bin, the gas-to-baryon fraction inside the virial radius does not evolve significantly with time, ranging from ~0.8 for smaller halos to ~0.5 for the largest ones. Most of the baryons are located actually not in the galaxies but in the intrahalo gas; for the more massive halos, the intrahalo gas-to-galaxy mass ratio is approximately the same at all redshifts, z, but for the least massive halos, it strongly increases with z. The intrahalo gas in the former halos gets hotter with time, being dominant at z = 0, while in the latter halos, it is mostly cold at all epochs. The multiphase ISM and thermal feedback models in our simulations work in the direction of delaying the stellar mass growth of low-mass galaxies.

We perform an exhaustive comparison among central galaxies from SDSS catalogs in different local environments at 0.01<=z<=0.08. The central galaxies are separated into two categories: group centrals (host halos containing satellites) and field centrals (host halos without satellites). From the latter, we select other two subsamples: isolated centrals and bright field centrals, both with the same magnitude limit. The stellar mass (Ms) distributions of the field and group central galaxies are different, which explains why in general the field central galaxies are mainly located in the blue cloud/star forming regions, whereas the group central galaxies are strongly biased to the red sequence/passive regions. The isolated centrals occupy the same regions as the bright field centrals since both populations have similar Ms distributions. At parity of Ms, the color and specific star formation rate (sSFR) distributions of the samples are similar, specially between field and group centrals. Furthermore, we find that the stellar-to-halo mass (Ms-Mh) relation of isolated galaxies does not depend on the color, sSFR and morphological type. For systems without satellites, the Ms-Mh relation steepens at high halo masses compared to group centrals, which is a consequence of assuming a one-to-one relation between group total stellar mass and halo mass. Under the same assumption, the scatter around the Ms-Mh relation of centrals with satellites increases with halo mass. Our results suggest that the mass growth of central galaxies is mostly driven by the halo mass, with environment and mergers playing a secondary role.

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.

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.

In order to attain a statistical description of the evolution of cosmic density fluctuations in agreement with results from the numerical simulations, we introduce a probability conditional formalism (CF) based on an inventory of isolated overdense regions in a density random field. This formalism is a useful tool for describing at the same time the mass function (MF) of dark haloes, their mass aggregation histories (MAHs) and merging rates (MRs). The CF focuses on virialized regions in a self-consistent way rather than in mass elements, and it offers an economical description for a variety of random fields. Within the framework of the CF, we confirm that, for a Gaussian field, it is not possible to reproduce at the same time the MF, MAH, and MR of haloes, both for a constant and moving barrier. Then, we develop an inductive method for constraining the cumulative conditional probability from a given halo MF description, and thus, using the CF, we calculate the halo MAHs and MRs. By applying this method to the MF measured in numerical simulations by Tinker et al. 2008, we find that a reasonable solution, justified by a mass conservation argument, is obtained if ones introduce a rescaling -increment by ~30% - of the virial mass used in simulations and a (slight) deviation from Gaussianity. Thus, both the MAH and MR obtained by a Monte Carlo merger tree agree now with the predictions of numerical simulations. We discuss on the necessity of rescaling the virial mass in simulations when comparing with analytical approaches on the ground of the matter not accounted as part of the halos and the halo mass limit due to numerical. Our analysis supports the presence of a diffuse dark matter component that is not taken into account in the measured halo MFs inasmuch as it is not part of the collapsed structures.

The mass assembly of a whole population of sub-Milky Way galaxies is studied by means of hydrodynamical simulations within the $\Lambda$-CDM cosmology. Our results show that while dark halos assemble hierarchically, in stellar mass this trend is inverted in the sense that the smaller the galaxy, the later is its stellar mass assembly on average. Our star formation and supernovae feedback implementation in a multi-phase interstellar medium seems to play a key role on this process. However, the obtained downsizing trend is not yet as strong as observations show.

Cosmological hydrodynamical simulations are studied in order to analyse generic trends for the stellar, baryonic and halo mass assembly of low-mass galaxies (M_* < 3 x 10^10 M_sun) as a function of their present halo mass, in the context of the Lambda-CDM scenario and common subgrid physics schemes. We obtain that smaller galaxies exhibit higher specific star formation rates and higher gas fractions. Although these trends are in rough agreement with observations, the absolute values of these quantities tend to be lower than observed ones since z~2. The simulated galaxy stellar mass fraction increases with halo mass, consistently with semi-empirical inferences. However, the predicted correlation between them shows negligible variations up to high z, while these inferences seem to indicate some evolution. The hot gas mass in z=0 halos is higher than the central galaxy mass by a factor of ~1-1.5 and this factor increases up to ~5-7 at z~2 for the smallest galaxies. The stellar, baryonic and halo evolutionary tracks of simulated galaxies show that smaller galaxies tend to delay their baryonic and stellar mass assembly with respect to the halo one. The Supernova feedback treatment included in this model plays a key role on this behaviour albeit the trend is still weaker than the one inferred from observations. At z>2, the overall properties of simulated galaxies are not in large disagreement with those derived from observations.

Several occupational distributions for satellite galaxies more massive than ms~4E7 Msun around MW-sized hosts are presented and used to predict their internal dynamics. For this, a galaxy group mock catalog is constructed on the basis of (sub)halo-to-stellar mass relations constrained with the observed Galaxy Stellar Mass Function decomposed into centrals and satellites, and the 2-point correlation functions. For the MW-sized hosts we find that: 6.6% have 2 satellites in the SMC-LMC mass range; their probabilities to have 1, 2 or 3 satellites equal or larger than the LMC, the SMC or Sgr, respectively, are 0. 26, 0. 14, and 0. 14; those hosts with 3 satellites >= Sgr (as the MW) are among the most common cases, though the most and 2nd most massive satellites in these systems are SMALLER than the LMC and SMC by 0.7 and 0.8 dex, respectively. The cumulative satellite mass function is quite broad, the particular configuration of the MW being of low frequency but NOT an outlier (it is within the upper 1sig yet). For the MW satellite configuration, we constrain the halo mass to be Mh>1.4E12 Msun (1sig level). We find that the abundance of massive subhalos should agree with the abundance of massive satellites in all MW-sized hosts (there is not a massive satellite missing problem for LCDM). However, the vmax of the subhalos of satellites with ms<1E8 Msun is systematically larger than the vmax inferred from current observational studies of the MW bright dwarfs; at difference of previous works, this conclusion is based on an analysis of the overall population of MW-sized galaxies. Some pieces of evidence suggest that the issue could refer only to satellite dwarfs but not to central dwarfs; then, environmental processes associated to dwarfs inside host halos combined with SN-driven core expansion could be at the basis of the lowering of vmax.[Abridged]

We infer the local stellar-to-halo/subhalo mass relations (MRs) for central and satellite galaxies separately. We constraint this relations by using several combinations of observational data, consisting of the total galaxy stellar mass function, its decomposition into centrals and satellites, and the projected two-point correlation functions measured in different stellar mass bins. The differences among the resulting MRs are within the model-fit uncertainties (which are very small, smaller than the intrinsic scatter between galaxy and halo mass), no matter what combination of data are used. We also constrain the values of the intrinsic scatter around the central-halo (CH) and satellite-subhalo (SS) MRs assuming them to be constant: sigma_c= 0.168+-0.051 dex and sigma_s = 0.172+-0.057 dex, respectively. The CH and SS MRs are actually different, in particular when we take the subhalo mass at the present-day epoch instead of at their accretion time. We calculate several population statistics. We find that the central galaxy is not on average within the mass distribution of the most-massive satellite, even for cluster-sized halos, i.e., centrals are not a mere realization of the high-end of the satellite mass function; however for > 3x10^13 msun halos, ~15% of centrals could be. We also find that the probabilities of MW-sized halos of having N Magellanic-Clouds (MCs)-sized satellites agree well with observational measures; for a halo mass of 2x10^12 msun, the probability to have 2 MCs is 5.4%, but if we exclude those systems with satellites larger than the MCs, then the probability decreases to <2.2%.

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.

We extend the abundance matching technique (AMT) to infer the satellite-subhalo and central-halo mass relations (MRs) of galaxies, as well as the corresponding satellite conditional mass functions (CMFs). We use the observed galaxy stellar mass function (GSMF) decomposed into centrals and satellites and the LCDM halo/subhalo mass functions as inputs. We explore the effects of defining the subhalo mass at the time of accretion (m_acc) vs. at the time of observation (m_obs). We test the standard assumption that centrals and satellites follow the same MRs, showing that this assumption leads to predictions in disagreement with observations, specially for m_obs. Instead, when the satellite-subhalo MRs are constrained following our AMT, they are always different from the central-halo MR: the smaller the stellar mass (Ms), the less massive is the subhalo of satellites as compared to the halo of centrals of the same Ms. On average, for Ms<2x10^11Msol, the dark mass of satellites decreased by 60-65% with respect to their masses at accretion time. The resulting MRs for both definitions of subhalo mass yield satellite CMFs in agreement with observations. Also, when these MRs are used in a HOD model, the predicted correlation functions agree with observations. We show that the use of m_obs leads to less uncertain MRs than m_acc, and discuss implications of the obtained satellite-subhalo MR. For example, we show that the tension between abundance and dynamics of MW satellites in LCDM gives if the slope of the GSMF faint-end slope upturns to -1.6.

An unified picture of stellar and halo mass build-up as a function of mass is presented. Inferred stellar-dark halo mass relations of galaxies, Ms-Mh, out to z=4 together with average LCDM halo mass aggregation histories (MAHs) are used for inferring average Ms growth histories, the Galaxian Hybrid Evolutionary Tracks (GHETs). The more massive the galaxy, the earlier transited in average from an active regime of Ms growth to a passive one: log(Mtran/Msun)=10.30+0.55z ("population downsizing"), where Mtran is the typical transition stellar mass. This result agrees with independent observational determinations based on the evolution of the galaxy stellar mass function decomposition into blue and red galaxies. The specific star formation rate, SSFR, predicted from the derivative of the GHET is consistent with direct measures of the SSFR for galaxies at different z's. The average GHETs of galaxies smaller than Mtran at z=0 (Ms~10^10.3 Msun) did not reach the quiescent regime, and for them, the lower the mass, the faster the later Ms growth rate ("downsizing in SSFR"). The GHETs allow to predict the transition rate in number density of active to passive population; the predicted values agree with direct estimates of growth rate in number density for the (massive) red population up to z~1. We show that LCDM-based models of disk galaxy evolution are able to reproduce the low-mass side of the Ms-Mh relation at z~0, but at higher z's disagree strongly with the GHETs: models do not reproduce the downsizing in SSFR and the high SSFR of low mass galaxies. (Abridged)

By means of the abundance matching technique, we infer the local stellar and baryonic mass-halo mass (Ms-Mh and Mb-Mh) relation for central blue and red galaxies separately in the mass range Ms~10^8.5-10^12.0 Msun. The observational inputs are the SDSS central blue and red Galaxy Stellar Mass Functions reported in Yang et al. 2009, and the measured local gas mass-Ms relations for blue and red galaxies. For the Halo Mass Function associated to central blue galaxies, the distinct LCDM one is used and set up to exclude: (i) the observed group/cluster mass function (blue galaxies are rare as centers of groups/clusters), and (ii) halos with a central major merger at resdshifts z<0.8 (dry late major mergers destroy the disks of blue galaxies). For red galaxies, we take the complement of this function to the total. The obtained mean Ms-Mh and Mb-Mh relations of central blue and red galaxies do not differ significantly from the respective relations for all central galaxies. For Mh>10^11.5 Msun, the Ms's of red galaxies tend to be higher than those of blue ones for a given Mh, the difference not being larger than 1.7. For Mh<10^11.5 Msun, this trend is inverted. For blue (red) galaxies: (a) the maximum value of fs=Ms/Mh is 0.021^{+0.016}_{-0.009} (0.034{+0.026}_{-0.015}) and it is attained atlog(Mh/Msun)~12.0 (log(Mh/Msun)~11.9); (b) fs\propto Mh (fs\propto Mh^3) at the low-mass end while at the high-mass end, fs\propto Mh^-0.4 (fs\propto Mh^-0.6). The baryon mass fractions, fb=Mb/Mh, of blue and red galaxies reach maximum values of fb=0.028^{+0.018}_{-0.011} and fb=0.034^{+0.025}_{-0.014}, respectively. For Mh<10^11.3 Msun, a much steeper dependence of fb on Mh is obtained for the red galaxies than for the blue ones. We discuss on the differences found in the fs-Mh and fb-Mh relations between blue and red galaxies in the light of of semi-empirical galaxy models.

(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.

The emerging empirical picture of galaxy stellar mass (Ms) assembly shows that galaxy population buildup proceeds from top to down in Ms. By connecting galaxies to LCDM halos and their histories, individual (average) Ms growth tracks can be inferred. These tracks show that massive galaxies assembled their Ms the earlier the more massive the halo, and that less massive galaxies are yet actively growing in Ms, the more active the less massive is the halo. The predicted star formation rates as a function of mass and the downsizing of the typical mass that separate active galaxies from the passive ones agree with direct observational determinations. This implies that the LCDM scenario is consistent with these observations. The challenge is now to understand the baryonic physics that drives the significant and systematical shift of the stellar mass assembly of galaxies from the mass assembly of their corresponding halos (from halo upsizing to galaxy downsizing).

We present numerical simulations aimed at exploring the effects of varying the sub-grid physics parameters on the evolution and the properties of the galaxy formed in a low-mass dark matter halo (~7 times 10^10 Msun/h at redshift z=0). The simulations are run within a cosmological setting with a nominal resolution of 218 pc comoving and are stopped at z = 0.43. In all of our simulations, an extended old/intermediate-age stellar halo and a more compact younger stellar disk are formed. We found that a non negligible fraction of the halo stars are formed in situ in a spheroidal distribution. Changes in the sub-grid physics parameters affect significantly and in a complex way the evolution and properties of the galaxy: (i) Lower threshold densities nsf produce larger stellar effective radii Re, less peaked circular velocity curves V_c(R), and greater amounts of low-density and hot gas in the disk mid-plane; (ii) When stellar feedback is modeled by temporarily switching off radiative cooling in the star forming regions, Re increases (by a factor of ~ 2 in our particular model), the circular velocity curve becomes flatter, and a complex multi-phase gaseous disk structure develops; (iii) A more efficient local conversion of gas mass to stars, measured by a stellar particle mass distribution biased toward larger values, increases the strength of the feedback energy injection -driving outflows and inducing burstier SF histories; iv) If feedback is too strong, gas loss by galactic outflows -which are easier to produce in low-mass galaxies- interrupts SF, whose history becomes episodic. The simulations exhibit two important shortcomings: the baryon fractions are higher, and the specific SF rates are much smaller, than observationally inferred values for redshifts ~ 0.4-1.

The observations show that less massive the galaxies are, the higher on average is their specific star formation rate (SSFR = SFR/Ms, Ms is the stellar mass). Such a trend, called the 'SSFR downsizing' (SSFR-DS) phenomenon, is seen for local and high-z (back to z~1-2) galaxy samples. We use observational data related only to disc galaxies and explore the average SSFR change with z for different masses. For Ms in the range ~10^9.5-10^10.5 Msun, the SSFR increases with (1+z) to a power that barely depends on Ms, and at all z's smaller galaxies have ever higher SSFRs. The latter strongly disagree with the LCDM hierarchical mass accretion rates. By means of self-consistent models of disc galaxy evolution inside growing LCDM halos, the effects that disc feedback-driven outflows and gas re-accretion have on the galaxy SSFR histories are explored. The parameters of the outflow and re-accretion schemes are tuned to reproduce the z~0 Mh-Ms relation inferred from observations. In the case of outflows only, the SSFR of individual model galaxies is roughly proportional to (1+z)^2.2 for all the masses with a normalization factor that depends on mass as Ms^0.1, i.e more massive galaxies have slightly larger SSFRs, contrary to the observed strong SSFR-DS trend. For the re-accretion cases, the dependence on z remains approximately the same as without re-infall, but the correlation on Ms even increases for most of the reasonable values of the model parameters. The comparison of models and observations in the SSFR-Ms plane at z~0 shows the divergent trend in SSFR as lower are the masses (upsizing vs downsizing). We explain why the models show the reported trends, and conclude that the SSFR-DS for low-mass galaxies poses a challenge for LCDM-based disc galaxy evolution models. (Abridged)

We are proposing to conduct a multicolor, synoptic infrared (IR) imaging survey of the Northern sky with a new, dedicated 6.5-meter telescope at San Pedro M\'artir (SPM) Observatory. This initiative is being developed in partnership with astronomy institutions in Mexico and the University of California. The 4-year, dedicated survey, planned to begin in 2017, will reach more than 100 times deeper than 2MASS. The Synoptic All-Sky Infrared (SASIR) Survey will reveal the missing sample of faint red dwarf stars in the local solar neighborhood, and the unprecedented sensitivity over such a wide field will result in the discovery of thousands of z ~ 7 quasars (and reaching to z > 10), allowing detailed study (in concert with JWST and Giant Segmented Mirror Telescopes) of the timing and the origin(s) of reionization. As a time-domain survey, SASIR will reveal the dynamic infrared universe, opening new phase space for discovery. Synoptic observations of over 10^6 supernovae and variable stars will provide better distance measures than optical studies alone. SASIR also provides significant synergy with other major Astro2010 facilities, improving the overall scientific return of community investments. Compared to optical-only measurements, IR colors vastly improve photometric redshifts to z ~ 4, enhancing dark energy and dark matter surveys based on weak lensing and baryon oscillations. The wide field and ToO capabilities will enable a connection of the gravitational wave and neutrino universe - with events otherwise poorly localized on the sky - to transient electromagnetic phenomena.

By means of galaxy evolutionary models, we explore the direct consequences of the LCDM cosmogony on the size evolution of galactic discs, avoiding intentionally the introduction of intermediate (uncertain) astrophysical processes. Based on the shape of the rotation curves and guided by a simplicity criterion, we adopt an average galaxy mass baryon fraction of 0.03. In order to study general behaviors, only models with the average initial conditions are analyzed. The stellar and B-band effective radii, R* and RB, of individual galaxies grow significantly with time (inside-out disc formation) with laws that are weakly dependent on mass, M*,or luminosity, LB. However, the change of R* with z at fixed M* is slow; for z<2.5, R*(M*=const) ~ (1+z)^-0.4. On the other hand, the change of RB with z at a fixed LB is strong and resembles the RB decreasing law of the individual models; roughly RB(LB=const) ~ (1+z)^-0.85 for z<0.75, and ~(1+z)^-1.1 for z>0.75. We find also that at z=0, R* ~ M*^0.38 and RB ~ LB^0.40, remaining the slopes of these relations almost the same up to z ~ 3. Our predictions are in reasonable agreement with observational inferences on the typical radius change with z of late-type galaxies more luminous (massive) than high values imposed by the selection effects. The models seem also to be consistent, within the large scatter, with the RB and LB values obtained from non complete samples of sub-L* late-type galaxies with available rest-frame photometric information at different z's. The properties and evolution of the LCDM haloes seem to be the main drivers of galaxy disc size evolution. Nevertheless, the models reveal a potential difficulty in explaining the observed steepening of the RB-LB relation with respect to the R*-M* one, an effect related to the well established color-magnitude relation.

For a sample of long GRBs with known redshift, we study the distribution of the evolutionary tracks on the rest-frame luminosity-peak energy Liso-Ep' diagram. We are interested in exploring the extension of the `Yonetoku' correlation to any phase of the prompt light curve, and in verifying how the high-signal prompt duration time, Tf, in the rest frame correlates with the residuals of such correlation (Firmani et al. 2006). For our purpose, we analyse separately two samples of time-resolved spectra corresponding to 32 GRBs with peak fluxes >1.8 phot cm^-2 s^-1 from the Swift-BAT detector, and 7 bright GRBs from the CGRO-BATSE detector previously processed by Kaneko et al. (2006). After constructing the Liso-Ep' diagram, we discuss the relevance of selection effects, finding that they could affect significantly the correlation. However, we find that these effects are much less significant in the Liso x Tf-Ep' diagram, where the intrinsic scatter reduces significantly. We apply further corrections for reducing the intrinsic scatter even more. For the sub-samples of GRBs (7 from Swift and 5 from CGRO) with measured jet break time, we analyse the effects of correcting Liso by jet collimation. We find that (i) the scatter around the correlation is reduced, and (ii) this scatter is dominated by the internal scatter of the individual evolutionary tracks. These results suggest that the time, integrated `Amati' and `Ghirlanda' correlations are consequences of the time resolved features, not of selection effects, and therefore call for a physical origin. We finally remark the relevance of looking inside the nature of the evolutionary tracks.

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 a revaluation of the optical morphology for 549 galaxies from the Catalog of Isolated Galaxies (Karachentseva 1973) that are available in the SDSS (DR6). Both the high resolution and high dynamic range of the SDSS images and our semi-automatic image processing scheme, allow for a major quality and uniform morphological analysis. The processing scheme includes (1) sky-subtracted, cleaned, and logarithmic scaled g-band images; (2) filtered-enhanced versions of the images in (1); and (3) the RGB images available in the SDSS database. Special attention was paid for distinguishing between E, S0, and Sa candidates through an additional analysis of (4) the surface brightness, \epsilon, PA, and A_4/B_4 Fourier series expansion profiles. The median type in the sample corresponds to Sbc, with 65% of the sample being of this type or later. A scarce population of E (3.5%) and S0 (5%) galaxies, and a non-negligible 14% fraction of SaSab spirals are identified. We calculate also the gri absolute magnitudes corrected by Galactic and internal extinctions and present the g-i color distribution and the color-magnitude diagram. Among the spirals, we find fractions of strong and suspected bars of 65.8% and of 33.3% of rings. For the E galaxies we find a richness of distinct substructure in their isophotal shape and also of morphological distortions. At least 78% of the ellipticals show some kind of distortion, suggesting that these galaxies suffered late dry mergers. Finally, we calculate the CAS (concentration, asymmetry and clumpiness) structural parameters in the ugriz bands. We analyze the loci of our isolated galaxies in different projections of the CAS volume and discuss some trends of these parameters with the color band, as well as with the type and the galaxy color. [abridged]

Long Gamma-Ray Bursts (GRBs) are the brightest electromagnetic explosions in the Universe, associated to the death of massive stars. As such, GRBs are potential tracers of the evolution of the cosmic massive star formation, metallicity, and Initial Mass Function. GRBs also proved to be appealing cosmological distance indicators. This opens a unique opportunity to constrain the cosmic expansion history up to redshifts 5-6. A brief review on both subjects is presented here.

We report a series of high-resolution cosmological N-body simulations designed to explore the formation and properties of dark matter halos with masses close to the damping scale of the primordial power spectrum of density fluctuations. We further investigate the effect that the addition of a random component, v_rms, into the particle velocity field has on the structure of halos. We adopted as a fiducial model the Lambda Warm Dark Matter cosmology with a non-thermal sterile neutrino mass of 0.5 keV. The filtering mass corresponds then to M_f = 2.6x10^12 M_sun/h. Halos of masses close to M_f were simulated with several million of particles. The results show that, on one hand, the inner density slope of these halos (at radii <~0.02 the virial radius Rvir) is systematically steeper than the one corresponding to the NFW fit or to the CDM counterpart. On the other hand, the overall density profile (radii larger than 0.02Rvir) is less curved and less concentrated than the NFW fit, with an outer slope shallower than -3. For simulations with v_rms, the inner halo density profiles flatten significantly at radii smaller than 2-3 kpc/h (<~0.010-0.015Rvir). A constant density core is not detected in our simulations, with the exception of one halo for which the flat core radius is ~1 kpc/h. Nevertheless, if ``cored'' density profiles are used to fit the halo profiles, the inferred core radii are ~0.1-0.8 kpc/h, in rough agreement with theoretical predictions based on phase-space constrains, and on dynamical models of warm gravitational collapse. A reduction of v_rms by a factor of 3 produces a modest decrease in core radii, less than a factor of 1.5. We discuss the extension of our results into several contexts, for example, to the structure of the cold DM micro-halos at the damping scale of this model.

We study the spectral and energetics properties of 47 long-duration gamma-ray bursts (GRBs) with known redshift, all of them detected by the Swift satellite. Due to the narrow energy range (15-150 keV) of the Swift-BAT detector, the spectral fitting is reliable only for fitting models with 2 or 3 parameters. As high uncertainty and correlation among the errors is expected, a careful analysis of the errors is necessary. We fit both the power law (PL, 2 parameters) and cut--off power law (CPL, 3 parameters) models to the time-integrated spectra of the 47 bursts, and present the corresponding parameters, their uncertainties, and the correlations among the uncertainties. The CPL model is reliable only for 29 bursts for which we estimate the nuf_nu peak energy Epk. For these GRBs, we calculate the energy fluence and the rest- frame isotropic-equivalent radiated energy, Eiso, as well as the propagated uncertainties and correlations among them. We explore the distribution of our homogeneous sample of GRBs on the rest-frame diagram E'pk vs Eiso. We confirm a significant correlation between these two quantities (the "Amati" relation) and we verify that, within the uncertainty limits, no outliers are present. We also fit the spectra to a Band model with the high energy power law index frozen to -2.3, obtaining a rather good agreement with the "Amati" relation of non-Swift GRBs.

A release of multicolor broad band (BVRI) photometry for a subsample of 44 isolated spirals drawn from the Catalogue of Isolated Galaxies (CIG) is presented. Total magnitudes and colors at various circular apertures, as well as some global structural/morphological parameters are estimated. Morphology is reevaluated through optical and sharp/filtered R band images, (B-I) color index maps, and archive near-IR JHK images from the Two-Micron Survey. The CAS structural parameters (Concentration, Asymmetry, and Clumpiness) were calculated from the images in each one of the bands. The fraction of galaxies with well identified optical/near-IR bars (SB) is 63%, while a 17% more shows evidence of weak or suspected bars (SAB). The sample average value of the maximum bar ellipticity is 0.4. Half of the galaxies in the sample shows rings. We identify two candidates for isolated galaxies with disturbed morphology. The structural CAS parameters change with the observed band, and the tendencies they follow with the morphological type and global color are more evident in the redder bands. In any band, the major difference between our isolated spirals and a sample of interacting spirals is revealed in the A-S plane. A deep and uniformly observed sample of isolated galaxies is intended for various purposes including (i) comparative studies of environmental effects, (ii) confronting model predictions of galaxy evolution and (iii) evaluating the change of galaxy properties with redshift.

The early X-ray light curve of Gamma Ray Bursts (GRBs) is complex, and shows a typical steep-flat-steep behaviour. The time T_a at which the flat (plateau) part ends may bear some important physical information, especially if it plays the same role of the so called jet break time t_jet. To this aim, stimulated by the recent analysis of Willingale et al., we have assembled a sample of GRBs of known redshifts, spectral parameters of the prompt emission, and T_a. By using T_a as a jet angle indicator, and then estimating the collimation corrected prompt energetics, we find a correlation between the latter quantity and the peak energy of the prompt emission. However, this correlation has a large dispersion, similar to the dispersion of the Amati correlation and it is not parallel to the Ghirlanda correlation. Furthermore, we show that the correlation itself results mainly from the dependence of the jet opening angle on the isotropic prompt energy, with the time T_a playing no role, contrary to what we find for the jet break time t_jet. We also find that for the bursts in our sample T_a weakly correlates with E_iso of the prompt emission, but that this correlation disappears when considering all bursts of known redshift and T_a. There is no correlation between T_a and the isotropic energy of the plateau phase.

Recently, a tight correlation among three quantities that characterize the prompt emission of long Gamma-Ray Bursts (GRBs) with known redshift z, was discovered (Firmani et al. 2006). We use this correlation to construct the Hubble diagram (HD) with a sample of 19 GRBs in the broad range of z=0.17-4.5, and carry out a full statistical analysis to constrain cosmological parameters (CPs). To optimally solve the problem of circularity, a Bayesian approach is applied. The main result is that the concordance LambdaCDM cosmology is fully consistent with the GRB data at the level of several tests. If we assume the Lambda cosmology, then we find Om_M=0.31^{+0.09}_{-0.08} and Om_Lambda=0.80^{+0.20}_{-0.30}$ (1sigma); the flat-geometry case is within 1sigma. Assuming flatness, we find Om_M=0.29^{+0.08}_{-0.06}, and fixing Om_M=0.28, we obtain a dark energy equation of state parameter w=-1.07^{+0.25}_{-0.38}, i.e. the ambdaCDM model (w=-1) is within 1sigma. Given the low number of usable GRBs we cannot yet constrain well the possible evolution of w=w(z). However, the case w(z)=-1 (LambdaCDM) is consistent at the 68.3% CL with GRBs. It is shown also how a broad range of z's in the used sample improves the determination of CPs from the HD, which is the case of GRBs as distance indicators.

Tight constraints on cosmological parameters can be obtained with standard candles spanning a range of redshifts as large as possible. We propose to treat SN Ia and long Gamma-Ray Bursts (GRBs) as a single class of candles. Taking advantage of the recent release of the Supernova Legacy Survey and {\it the recent finding of a tight correlation among the energetics and other prompt gamma-ray emission properties of GRBs}, we are able to standardize the luminosities/energetics of both classes of objects. In this way we can jointly use GRB and SNIa as cosmological probes to constrain Omega_m and Omega_L and the Dark Energy equation of state parameters through the same Bayesian method that we have, so far, applied to GRBs alone. Despite the large disparity in number (115 SNIa versus 19 GRBs) we show that the constraints on Omega_m and Omega_L are greatly improved by the inclusion of GRBs. More importantly, the result of the combined sample is in excellent agreement with the Lambda-CDM concordance cosmological model and does not require an evolving equation of state for the Dark Energy.

We show by means of a high-resolution N-body simulation how the mass assembly histories of galaxy-size cold dark matter (CDM) halos depend on environment. Halos in high density environments form earlier and a higher fraction of their mass is assembled in major mergers, compared to low density environments. The distribution of the present-day specific mass aggregation rate is strongly dependent on environment. While in low density environments only ~20% of the halos are not accreting mass at the present epoch, this fraction rises to ~80% at high densities. At z=1 the median of the specific aggregation rate is ~4 times larger than at z=0 and almost independent on environment. All the dependences on environment found here are critically enhanced by local processes associated to subhalos because the fraction of subhalos increases as the environment gets denser. The distribution of the halo specific mass aggregation rate as well as its dependence on environment resemble the relations for the specific star formation rate distribution of galaxies. An analogue of the morphology-density relation is also present at the level of CDM halos, being driven by the halo major merging history. Nevertheless, baryonic processes are necessary in order to explain further details and the evolution of the star formation rate-, color- and morphology-environment relations.

The old dream of integrating into one the study of micro and macrocosmos is now a reality. Cosmology, astrophysics, and particle physics intersect in a scenario (but still not a theory) of cosmic structure formation and evolution called Lambda Cold Dark Matter (LCDM) model. This scenario emerged mainly to explain the origin of galaxies. In these lecture notes, I first present a review of the main galaxy properties, highlighting the questions that any theory of galaxy formation should explain. Then, the cosmological framework and the main aspects of primordial perturbation generation and evolution are pedagogically detached. Next, I focus on the ``dark side'' of galaxy formation, presenting a review on LCDM halo assembling and properties, and on the main candidates for non-baryonic dark matter. It is shown how the nature of elemental particles can influence on the features of galaxies and their systems. Finally, the complex processes of baryon dissipation inside the non-linearly evolving CDM halos, formation of disks and spheroids, and transformation of gas into stars are briefly described, remarking on the possibility of a few driving factors and parameters able to explain the main body of galaxy properties. A summary and a discussion of some of the issues and open problems of the LCDM paradigm are given in the final part of these notes.

Possible detection of signatures of structure formation at the end of the 'dark age' epoch (z~40-20) is examined. We discuss the spectral-spatial fluctuations in the CMBR temperature produced by elastic resonant scattering of CMBR photons on HD molecules located in protostructures moving with peculiar velocity. Detailed chemical kinematic evolution of HD molecules in the expanding homogeneous medium is calculated. Then, the HD abundances are linked to protostructures at their maximum expansion, whose properties are estimated by using the top-hat spherical approach and the LambdaCDM cosmology. We find that the optical depths in the HD three lowest pure rotational lines for high-peak protohaloes at their maximum expansion are much higher than those in LiH molecule. The corresponding spectral-spatial fluctuation amplitudes however are probably too weak as to be detected by current and forthcoming millimeter-telescope facilities. We extend our estimates of spectral-spatial fluctuations to gas clouds inside collapsed CDM haloes by using results from a crude model of HD production in these clouds. The fluctuations for the highest-peak CDM haloes at redshifts ~20-30 could be detected in the future. Observations will be important to test model predictions of early structure formation in the universe.

We report the discovery of a correlation among three prompt emission properties of GRBs. These are the isotropic peak luminosity L_iso, the peak energy of the time-integrated prompt emission spectrum E_pk, and the ``high signal" timescale T_0.45, previously used to characterize the variability behavior of bursts. In the rest frame of the source the found correlation reads L_iso\propto E_pk^1.62 T_0.45^-0.49. We find other strong correlations, but at the cost of increasing the number of variables, involving the variability and the isotropic energy of the prompt emission. With respect to the previous tight correlations found in GRBs the newly found correlation does not require any information from the afterglow phase of the bursts, nor any model-dependent assumption. In the popular scenario in which we are receiving beamed radiation originating in a fireball pointing at us, the found correlation preserves its form in the comoving frame. This helps to explain the small scatter of the correlation, and underlines the role of the local brightness (i.e. the brightness of the visible fraction of the fireball surface). This correlation has been found for 19 objects, and it is hard to establish if any selection bias affects it. Its connection with the prompt local brightness is promising, but a solid physical understanding is still to be found. Despite all that, we find that some properties of the correlation, which we discuss, support its true existence, and this has important implications for the GRB physics. Furthermore, it is possible to use such correlation as an accurate redshift estimator, and its tightness will allow us to use it as a tool to constrain the cosmological parameters (abridged)

We present a photometric and structural analysis of 42 E/S0 galaxies in (E/S0 + S) pairs observed in the BVRI color bands. We empirically determine the effects of interactions on their morphology, structure and stellar populations as seen from the light concentration (C), asymmetry (A), and clumpiness (S) parameters. We further compare these values to a control sample of 67 mostly isolated, non-interacting E/S0 galaxies. The paired E/S0 galaxies occupy a more scattered loci in CAS space than non-interacting E/S0's, and the structural effects of interactions on E/S0's are minor, in contrast to disk galaxies involved in interactions. This suggests that observational methods for recognizing interactions at high z, such the CAS methodology, would hardly detect E/S0's involved in interactions (related to early phases of the so called `dry-mergers'). We however find statistical differences in A when comparing isolated and interacting E/S0s. In the mean, paired E/S0 galaxies have A values 2.96+-0.72 times larger than the ones of non-interacting E/S0's. For the subset of presumably strongly interacting E/S0's, A and S can be several times larger than the typical values of the isolated E/S0's. We show that the asymmetries are consistent with several internal and external morphological distortions. We conclude that the interacting E/S0s in pairs should be dense, gas poor galaxies in systems spaning a wide range of interaction stages, with typical merging timescales >~ 0.1-0.5 Gyr. We use the observed phenomenology of these galaxies to predict the approximate loci of `dry pre-mergers' in the CAS space.(Abridged)

High-resolution LCDM cosmological N-body simulations are used to study the properties of galaxy-size dark halos in different environments (cluster, void, and "field"). Halos in clusters and their surroundings have a median spin parameter ~1.3 times lower, and tend to be more spherical and to have less aligned internal angular momentum than halos in voids and the field. For halos in clusters the concentration parameters decrease on average with mass with a slope of ~0.1; for halos in voids these concentrations do not change with mass. For masses <5 10^11 M_sh^-1, halos in clusters are on average ~30-40% more concentrated and have ~2 times higher central densities than halos in voids. When comparing only parent halos, the differences are less pronounced but they are still significant. The Vmax-and Vrms-mass relations are shallower and more scattered for halos in clusters than in voids, and for a given Vmax or Vrms, the mass is smaller at z=1 than at z=0 in all the environments. At z=1, the differences in the halo properties with environment almost dissapear, suggesting this that the differences were stablished mainly after z~1. The halos in clusters undergo more dramatic changes than those in the field or the voids. The differences with environment are owing to (i) the dependence of halo formation time on environment, and (ii) local effects as tidal stripping and the tumultuos histories that halos suffer in high-density regions. We calculate seminumerical models of disk galaxy evolution in halos with the properties found for the different environments. For a given disk mass, the galaxy disks have higher surface density, larger Vd,max and secular bulge-to-disk ratio, lower gas fraction, and are redder as one goes from cluster to void environments, in rough agreement with observations. (abridged)

[Abridged] Very recently, relations between the peak energy of Gamma-Ray Burst burst spectra, the isotropic-equivalent energy of the burst, and the radiated energy of the burst have been found. In a way that is exactly analogous to the way in which the relation between the peak luminosity and the rate of decline of the light curve of Type Ia supernovae can be used to make Type Ia supernovae excellent standard candles for cosmology, so too, the relations between Epk, Eiso, and Egamma point toward a methodology for using GRBs as excellent standard candles for cosmology. In addition, GRBs occur over the broad redshift range from z=0.1 to at least z=4.5, and both they and their afterglows are easily detectable out to z > 8. Thus GRBs show great promise as cosmological ``yardsticks'' to measure the rate of expansion of the universe over time, and therefore the properties of dark energy (i.e., Omega_M, Omega_Lambda, w_0, and w_a). We describe a concept for a possible MIDEX-class mission dedicated to using GRBs to constrain the properties of dark energy that would obtain these quantities for > 800 bursts in the redshift range 0.1 >~ z <~ 10 during a 2-year mission. This burst sample would enable both Omega_M and w_0 to be determined to +/- 0.07 and +/- 0.06 (68% CL), respectively, and w_a to be significantly constrained.

It was argued in the past that bulges of galaxies cannot be formed through collisionless secular evolution because that would violate constraints on the phase-space density: the phase-space density in bulges is several times larger than in the inner parts of discs. We show that these arguments against secular evolution are not correct. Observations give estimates of the coarsely grained phase-space densities of galaxies, f'=rho_s/(sigma_R sigma_phi sigma_z), where rho_s is stellar density and sigma_R, sigma_phi, sigma_z are the radial, tangential, and vertical rms velocities of stars. Using high-resolution N-body simulations, we study the evolution of f' in discs of Galaxy-size models. During the secular evolution, the discs, which are embedded in live CDM haloes, form a bar and then a thick, dynamically hot, central mass concentration. In the course of evolution f' declines at all radii, not just in the central region. However, the decline is different in different parts of the disc. In the inner disc, f'(R) develops a valley with a minimum around the end of the central mass concentration. The final result is that the values of f' in the central regions are significantly larger than those in the inner disc. The minimum, which gets deeper with time, seems to be due to a large phase mixing produced by the outer bar. We find that the shape and the amplitude of f'(R) for different simulations agree qualitatively with the observed f'(R) in our Galaxy. Curiously enough, the fact that the coarsely grained phase-space density of the bulge is significantly larger than the one of the inner disc turns out to be an argument in favor of secular formation of bulges, not against it.

We report new calculations of the cooling rate of primordial gas by the HD molecule, taking into account its ro-vibrational structure. The HD cooling function is calculated including radiative and collisional transitions for J\le8 rotational levels, and for the vibrational levels v=0,1,2 and 3. The ro-vibrational level population is calculated from the balance equation assuming steady state. The cooling function is evaluated in the ranges of the kinetic temperatures, T_k, from 1e2 to 2e4 K and the number densities, n_H, from 1 to 1e8 cm**-3. We find that the inclusion of collisional ro-vibrational transitions increases significantly the HD cooling efficiency, in particular for high densities and temperatures. For n_H>~1e5 and T_k~1e4 K the cooling function becomes more than an order of magnitude higher than previously reported. We give also the HD cooling rate in the presence of the cosmic microwave radiation field for radiation temperatures of 30, 85 and 275 K (redshifts of 9, 29 and 99). The tabulated cooling functions are available at http://www.cifus.uson.mx/Personal\_Pages/anton/DATA/HD\_cooling/HD\_cool.html. We discuss the relevance to explore the effects of including our results into models and simulations of galaxy formation, especially in the regime when gas cools down from temperatures above ~3000 K.

The isotropic luminosity function (LF) and formation rate history (FRH) of long GRBs is by the first time constrained by using jointly both the observed GRB peak-flux and redshift distributions. Our results support an evolving LF and a FRH that keeps increasing after z=2. We discuss some interesting implications related to these results.

Apparently, Gamma-Ray Bursts (GRBs) are all but standard candles. Their emission is collimated into a cone and the received flux depends on the cone aperture angle. Fortunately we can derive the aperture angle through an achromatic steepening of the lightcurve of the afterglow, and thus we can measure the "true" energetics of the prompt emission. Ghirlanda et al. (2004) found that this collimation-corrected energy correlates tightly with thefrequency at which most of the radiation of the prompt is emitted. Through this correlation we can infer the burst energy accurately enough for a cosmological use. Using the best known 15 GRBs we find very encouraging results that emphasize the cosmological GRB role. Probing the universe with high accuracy up to high redshifts, GRBs establish a new insight on the cosmic expanding acceleration history and accomplish the role of "missing link" between the Cosmic Microwave Background and type Ia supernovae, motivating the most optimistic hopes for what can be obtained from the bursts detected by SWIFT.

We present a new method aimed to handle long Gamma-Ray Burst (GRBs) as cosmic rulers. The recent discovery of a tight correlation between the collimation corrected GRB energy and the peak of the gamma-ray spectrum has opened the possibility to use GRBs as a new category of standard candles. Unfortunately, because of the lack of low-z GRBs, up to now this correlation is obtained from high-z GRBs with the consequence that it depends on the cosmological parameters we pretend to constrain. Hopefully this circularity problem will be solved when, in a few years, the low-z GRB sample will be increased enough. In the meanwhile we present here a new Bayesian method that eases the aforesaid circularity problem, and allows to introduce new constrains on the cosmological (Om,OL) diagram as well as to explore the universe kinematics up to z~3. The method we propose offers the further advantage to make handy the problem of the (Om,OL) loitering line singularity which inevitably appears when standard candles with z>2 are used. The combination of GRB with SN Ia data makes the popular LambdaCDM cosmology more consistent with the Hubble diagram at a 68% confidence level. For a flat cosmology we find Om=0.28\pm0.03 for the combined GRB+SN Ia data set. Correspondingly, the transition redshift between cosmic deceleration and acceleration is z_T=0.73\pm0.09, slightly larger than the value found by considering SNe Ia alone. We briefly discuss our results also in terms of non--LambdaCDM dark energy models.

We present a structural analysis of nearby galaxies in spiral-spiral pairs in optical BVRI bands and compare with the structures of isolated spiral galaxies and galaxies in ongoing mergers. We use these comparisons to determine how galaxy structure changes during galaxy interactions and mergers. We analyze light concentration (C), asymmetry (A), and clumpiness (S) parameters, and use the projections of CAS parameter space to compare these samples. We find that the CAS parameters of paired galaxies are correlated with the projected separations of the pair. For the widest and closest pairs, the CAS parameters tend to be similar to those of isolated and ongoing major mergers (ULIRGs), respectively. Our results imply that galaxy CAS morphology is a robust quantity that only changes significantly during a strong interaction or major merger. The typical time-scale for this change in our paired sample, based on dynamical friction arguments, is short, t ~ 0.1 - 0.5 Gyr. We find average enhancement factors for the spiral pair asymmetries and clumpiness values of ~2.2 and 1.5. The S parameter, which is related to star formation activity, has a moderate level of enhancement suggesting that this activity in modern spirals depends more on internal processes than on external conditions. We furthermore test the statistical criterion for picking up interacting galaxies in an automated way by using the A-S projection plane. The diversity of our spiral pair sample in the CAS space suggests that structural/SF/morphological properties of interacting galaxies change abruptly only when the interaction becomes very strong and that the criteria for finding galaxies involved in major mergers from Conselice (2003) is effective.

We constrain the isotropic luminosity function (LF) and formation rate of long gamma-ray bursts (GRBs) by fitting models 'jointly' to both the observed differential peak flux and redshift distributions. We find evidence supporting an evolving LF, where the luminosity scales as (1+z)^delta with an optimal delta of 1.0\pm 0.2. For a single power-law LF, the best slope is ~ -1.57 with an upper luminosity of 10^53.3 erg/s, while the best slopes for a double power-law LF are approximately -1.6 and -2.6 with a break luminosity of 10^52.7 erg/s. Our finding implies a jet model intermediate between the universal structured epsilon(theta) proportional to theta^-2 model and the quasi-universal Gaussian structured model. For the uniform jet model our result is compatible with an angle distribution between 2 and 15 grades. Our best constrained GRB formation rate histories increase from z=0 to z=2 by a factor of ~30 and then continue increasing slightly. We connect these histories to that of the cosmic star formation history, and compare with observational inferences up to z~6. GRBs could be tracing the cosmic rates of both the normal and obscured star formation regimes. We estimate a current GRB event rate in the Milky Way of ~5 10^-5 yr^-1, and compare it with the birthrate of massive close WR+BH binaries with orbital periods of hours. The agreement is rather good suggesting that these systems could be the progenitors of the long GRBs.

We study the structure and substructure of halos obtained in N-body simulations for a Lambda Cold Dark Matter (LCDM) cosmology with non-Gaussian initial conditions (NGICs). The initial statistics are lognormal in the gravitational potential field with positive (LNp) and negative (LNn) skewness; the sign of the skewness is conserved by the density field, and the power spectrum is the same for all the simulations. Our aim is not to test a given non-Gaussian statistics, but to explore the generic effect of positive- and negative-skew statistics on halo properties. From our low-resolution simulations, we find that LNp (LNn) halos are systematically more (less) concentrated than their Gaussian counterparts. This result is confirmed by our Milky Way- and cluster-sized halos resimulated with high-resolution. In addition, they show inner density profiles that depend on the statistics: the innermost slopes of LNp (LNn) halos are steeper (shallower) than those obtained from the corresponding Gaussian halos. A subhalo population embedded in LNp halos is more susceptible to destruction than its counterpart inside Gaussian halos. On the other hand, subhalos in LNn halos tend to survive longer than subhalos in Gaussian halos. The spin parameter probability distribution of LNp (LNn) halos is skewed to smaller (larger) values with respect to the Gaussian case. Our results show how the statistics of the primordial density field can influence some halo properties, opening this the possibility to constrain, although indirectly, the primordial statistics at small scale.

High-resolution numerical simulations were performed to study the structure and substructure of Milky Way- and cluster-sized halos in a LCDM cosmology with self-interacting (SI) dark particles, where the particle cross section, +AFw-sig, is assumed constant or inversely proportional to the relative velocity. We conclude that the cuspy halo problem at galaxy scales of the LCDM cosmogony can be solved in the latter case. In this case, the inner density of galaxy-cluster halos is only slightly affected by the SI, resulting in agreement with observational constraints. At the same time, the subhalo population in galaxy and cluster SI-CDM halos remains roughly similar to that seen on collisionless CDM halos.

The study of formation and evolution of galaxies is reviewed, making emphasis on the physical factors which are important to understand the origin of the galaxy Hubble sequence. We concentrate on predictions of the hierarchical Cold Dark Matter (CDM) scenario and the confrontation with observations. The mass assembling of the CDM halos, the baryonic processes within them, and the evolution of disks and spheroids are described. The successes and shortcomings are discussed. Disk evolution seems to be a quiescent and extended process driven by the cosmological initial conditions, while spheroids are formed probably in violent events, where several astrophysical processes are competing.

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)

The kinetic energy dissipation rate in the turbulent ISM of disk galaxies is a key ingredient in galaxy evolution models since it determines the effectiveness of large-scale star formation (SF) feedback. Using magneto-hydro-dynamic simulations, we find that the ISM dissipates efficiently the turbulent kinetic energy injected by sources of stellar nature. Thus, the SF process may be self-regulated by an energy balance only at the level of the disk ISM. The use of the self-regulation SF mechanism in galaxy evolutionary models, where disks form inside growing Cold Dark Matter halos, allows to predict the SF history of disk galaxies, including the Milky Way and the solar neighborhood, as well as the contribution of the whole population of disk galaxies to the cosmic SF history. The results are encouraging.

We study the structure of Milky Way (MW)- and cluster-sized halos in a Lambda Cold Dark Matter (CDM) cosmology with self-interacting (SI) dark particles. The cross section per unit of particle mass has the form sigma = sig_0(1/v_100)^alpha, where sig_0 is a constant in units of cm^2/gr and v_100 is the relative velocity in units of 100 km/s. Different values for sigma with alpha= 0 or 1 were used. For small values of sigma = const. (sig_0<0.5), the core density of the halos at z=0 is typically higher at a given mass for lower values of sig_0 or, at a given sig_0, for lower masses. For values of sig_0 as high as 3.0, the halos may undergo the gravothermal catastrophe before z=0. When alpha = 1, the core density of cluster- and MW-sized halos is similar. Using sigma = 0.5-1.0x(1/v_100), our predictions agree with the central densities and the core scaling laws of halos both inferred from the observations of dwarf and LSB galaxies and clusters of galaxies. The cumulative Vmax-functions of subhalos in MW-sized halos with (sig_0,alpha) = (0.1,0.0), (0.5,0.0) and (0.5,1.0) agree roughly with observations (luminous satellites) for Vmax > 30 km/s, while at Vmax = 20 km/s the functions are a factor 5-8 higher, similar to the CDM predictions. The halos with SI have slightly more specific angular momentum at a given mass shell and are rounder than their CDM counterparts. We conclude that the introduction of SI particles with sigma \propto 1/v_100 may remedy the cuspy core problem of the CDM cosmogony, while the subhalo population number remains similar to that of the CDM halos.

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.

We find that the density distribution of high and low surface brightness (SB) galaxy disks formed within LCDM halos under the assumption of detailed angular momentum conservation is in rough agreement with observations. The luminous-to-dark matter ratio of the model galaxies increases with SB. The lowest SB models tend to be minimum disks, but the high SB models hardly attain the maximum disk solution. With the introduction of shallow cores in the halos, high SB models become maximum disks. The shallow cores also help to improve the inner SB profiles of bulge-less low SB models and the zero-point of the Tully-Fisher relation. The models predict well this relation and its scatter, as well as the small correlation among the residuals of this and the luminosity-radius relation, in spite of the dependence of the rotation curve shapes on SB.

We use a cosmological galactic evolutionary approach to model the Milky Way. A detailed treatment of the mass aggregation and dynamical history of the growing dark halo is included, together with a self consistent physical treatment for the star formation processes within the growing galactic disc. This allows us to calculate the temporal evolution of star and gas surface densities at all galactic radii, in particular, the star formation history (SFH) at the solar radius. A large range of cosmological mass aggregation histories (MAHs) is capable of producing a galaxy with the present day properties of the Milky Way. The resulting SFHs for the solar neighbourhood bracket the available observational data for this feature, the most probable MAH yielding the optimal comparison with these observations. We also find that the rotation curve for our Galaxy implies the presence of a constant density core in its dark matter halo.

(Abridged) Using high-resolution cosmological N-body simulations, we study how the density profiles of dark matter halos are affected by the filtering of the density power spectrum below a given scale length and by the introduction of a thermal velocity dispersion. In the warm dark matter (WDM) scenario, both the free-streaming scale, R_f, and the velocity dispersion, v_w, are determined by the mass m_w of the WDM particle. We find that v_w is too small to affect the density profiles of WDM halos. Down to the resolution attained in our simulations, there is not any significant difference in the density profiles and concentrations of halos obtained in simulations with and without the inclusion of v_w. The density profiles of halos with masses down to ~0.01 the filtering mass M_f can be described by the NFW shape; significant soft cores are not formed. Nevertheless, the concentrations of these halos are lower than those of the CDM counterparts and are approximately independent of mass. The lower concentrations of WDM halos with respect to their CDM counterparts can be accounted for their late formation epoch. Overall, our results point to a series of advantages of a WDM model over the CDM one. In addition to solving the substructure problem, a WDM model with R_f~0.16 Mpc (m_w~0.75 kev; flat cosmology with Omega_L=h=0.7) also predicts concentrations, a Tully-Fisher relation, and formation epochs for small halos which seems to be in better agreement with observations, relative to CDM predictions.

Using a cosmological galactic evolutionary approach to model the Milky Way, we calculate the star formation history (SFH) of the solar neighborhood. The good agreement we obtain with the observational inferences suggests that our physical model describes accurately the long term/large spatial trends of the local and global Milky Way SFH. In this model, star formation is triggered by disk gravitational instabilities and self-regulated by an energy balance in the ISM. The drivers of the SFH are the cosmological gas infall rate and the gas surface density determined by the primordial spin parameter. A LambdaCDM cosmology was used throughout.

We discuss the dissipation of turbulent kinetic energy Ek in the global ISM by means of 2-D, MHD, non-isothermal simulations in the presence of model radiative heating and cooling. We argue that dissipation in 2D is representative of that in three dimensions as long as it is dominated by shocks rather than by a turbulent cascade. Energy is injected at a few isolated sites in space, over relatively small scales, and over short time periods. This leads to the coexistence of forced and decaying regimes in the same flow. We find that the ISM-like flow dissipates its turbulent energy rapidly. In simulations with forcing, the input parameters are the radius l_f of the forcing region, the total kinetic energy e_k each source deposits into the flow, and the rate of formation of those regions, sfr_OB. The global dissipation time t_d depends mainly on l_f. In terms of measurable properties of the ISM, t_d >= Sigma_g u_rms^2/(e_k sfr_OB), where Sigma_g is the average gas surface density and u_rms is the rms velocity dispersion. For the solar neighborhood, t_d >= 1.5x10^7 yr. The global dissipation time is consistently smaller than the crossing time of the largest energy-containing scales. In decaying simulations, Ek decreases with time as t^-n, where n~0.8-0.9. This suggests a decay with distance d as Ek\propto d^{-2n/(2-n)} in the mixed forced+decaying case. If applicable to the vertical direction, our results support models of galaxy evolution in which stellar energy injection provides significant support for the gas disk thickness, but not models of galaxy formation in which this energy injection is supposed to reheat an intra-halo medium at distances of up to 10-20 times the optical galaxy size, as the dissipation occurs on distances comparable to the disk height.

High-resolution cosmological N-body simulations were performed in order to study the substructure of Milky Way-like galactic halos and the density profiles of halos in a warm dark matter scenario. The results favor this scenario with respect to the cold dark matter one.

We highlight some results from disk galaxy evolution models conceived within a cosmological context. When disk mergers and strong disk-halo feedback are omitted, several properties and correlations of disk galaxies seem to be related to initial conditions given by the CDM model, for example, the intensive galaxy properties, the disk Hubble sequence, and the Tully-Fisher relation in the infrared bands.

One of the predictions of the standard CDM is that dark haloes have centrally divergent density profiles. An extensive body of rotation curve observations of dwarf and low surface brightness galaxies shows the dark haloes of those systems to be characterized by soft constant density central cores. Several physical processes have been proposed to produce soft cores in dark haloes, each one with different scaling properties. With the aim of discriminating among them we have examined the rotation curves of dark matter dominated dwarf and low surface brightness galaxies and the inner mass profiles of two clusters of galaxies lacking a central cD galaxy and with evidence of soft cores in the centre. The core radii and central densities of these haloes scale in a well defined manner with the depth of their potential wells, as measured through the maximum circular velocity. As a result of our analysis we identify self-interacting CDM as a viable solution to the core problem, where a non-singular isothermal core is formed in the halo center surrounded by a Navarro, Frenk, & White profile in the outer parts. We show that this particular physical situation predicts core radii in agreement with observations. Furthermore, using the observed scalings, we derive an expression for the minimum cross section (\sigma) which has an explicit dependence with the halo dispersion velocity (v). If m_x is the mass of the dark matter particle: \sigma/m_x ~4 10^-25 (v/100 km s^-1)^-1 cm^2/Gev.

Using observational data for a sample of dark matter dominated galaxies and two cluster of galaxies, we have found that the central halo density does not depend on its mass, and the core radius is roughly proportional to the maximum rotation velocity. A good agreement with these scaling laws is obtained for CDM halos whose dense inner parts were expanded by gravothermal instabilities if the particles efficiently self-interact only in these parts. We find that the particle cross-section is inversely proportional to the velocity dispersion.

We study the evolution of disk galaxies using galaxy evolutionary models with initial and boundary conditions linked to the hierarchical formation scenario. Disks galaxies are modeled locally within growing cold dark matter halos and including a physical model for star formation. We focus our attention on predictions of the star formation history, size and surface brightness evolution and the evolution of the H- and B-band Tully-Fisher relations. Comparisons with available observational data are presented.

We performed a series of high-resolution simulations designed to study the substructure of Milky Way-size galactic halos (host halos) and the density profiles of halos in a warm dark matter (WDM) scenario with a non-vanishing cosmological constant. The virial masses of the host halos range from 3.5 x 10^12 to 1.7 x 10^12 solar masses and they have more than 10^5 particles each. A key feature of the WDM power spectrum is the free-streaming length R_f which fixes an additional parameter for the model of structure formation. We analyze the substructure of host halos using three R_f values: 0.2, 0.1, and 0.05 Mpc and compare results to the predictions of the cold dark matter (CDM) model. We find that guest halos (satellites) do form in the WDM scenario but are more easily destroyed by dynamical friction and tidal disruption than their counterparts in a CDM model. The small number of guest halos that we find within the virial radii of host halos at z = 0 in the WDM models is the result of a less efficient halo accretion and a higher satellite destruction rate. Under the assumption that each guest halo hosts a luminous galaxy, we find that the observed circular velocity function of satellites around the Milky Way and Andromeda is well described by the R_f = 0.1 Mpc WDM model. In the R_f = 0.1-0.2 Mpc models, the surviving subhalos at z=0 have an average concentration parameter c_1/5 which is approximately twice smaller than that of the corresponding CDM subhalos. This difference, very likely, produces the higher satellite destruction rate found in the WDM models. The density profile of host halos is well described by the NFW fit whereas guest halos show a wide variety of density profiles (abridged).

We study turbulent dissipation in the ISM and explore some implications for galaxy formation and evolution using 2D MHD numerical simulations of compressible fluids. The turbulent kinetic energy E_k is injected by stellar sources formed self-consistently in the simulation. In the ISM-like fluid, regimes of both forced and decaying turbulence coexist. In the active turbulent regions (forced regime), E_k is dissipated locally and efficiently. In the decaying regime (far from input sources), E_k(t) decays ~(1+t)^-0.8. The residual turbulent motions may propagate distances of the order of the observed disk height, suggesting that turbulence may be the responsible of vertical support and star formation self-regulation at the disk level, but not at the level of the whole cosmological halo, as would be required in some models of galaxy formation.

Within the framework of the cold dark matter (CDM) cosmogony, a central cusp in the density profiles of virialized dark haloes is predicted. This prediction disagrees with the soft inner halo mass distribution inferred from observations of dwarf and low surface brightness galaxies, and some clusters of galaxies. By analysing data for some of these objects, we find that the halo central density is nearly independent of the mass from galactic to galaxy cluster scales with an average value of around 0.02 M_sun/pc^3. We show that soft cores can be produced in the CDM haloes by introducing a lower cut-off in the power spectra of fluctuations and assuming high orbital thermal energies during halo formation. However, the scale invariance of the halo central density is not reproduced in these cases. The introduction of self-interaction in the CDM particles offers the most attractive alternative to the core problem. We propose gravothermal expansion as a possible mechanism to produce soft cores in the CDM haloes with self-interacting particles. A global thermodynamical equilibrium can explain the central density scale invariance. We find a minimum cross section capable of establishing isothermal cores in agreement with the observed shallow cores. If \sigma and m_x are the cross section and mass of the dark matter particle, and v is the halo velocity dispersion, then (\sigma/m_x) \~ 4 10^{-25} (100 km s^{-1}/v) cm^2/GeV.

We present results of extensive model calculations of disk galaxy evolution within an hierarchical inside-out formation scenario. We first compare properties of the dark halos identified in a cosmological N-body simulation with predictions of a seminumerical method based on an extended collapse model and find a good agreement. We also study the properties of the halos in dependence on their environment. We then describe detailed modelling of the formation and evolution of disks within the growing isolated cold dark matter halos and predictions for the main properties, correlations and evolutionary features of normal disk galaxies. The possible shortcomings of the scenario are discussed.

We propose and show that the dependence of the internal face-on dust absorption upon B-band luminosity reported by Wang & Heckman (1996) for a sample of late-type galaxies may explain the empirical color-magnitude and color Tully-Fisher relations of disk galaxies. Thus, in order to explain these relations, it is not necessary to evoke star formation and gas infall efficiencies dependent on the mass of the galaxy system. After applying the Wang & Heckman's luminosity-dependent extinction to hierarchical inside-out galaxy formation and evolution models where the star formation and gas infall efficiencies do not significantly depend on mass, we succesfully predict the observed Tully-Fisher relations in the H and B bands as well as the color-magnitude relation.

We used galaxy evolutionary models in a hierarchical inside-out formation scenario to study the origin of the main properties and correlations of disk galaxies. We found that most of these properties and correlations are the result of three (cosmological) initial factors and their dispersions: the virial mass, the halo mass aggregation history (MAH), and the angular momentum given through the spin parameter \lambda. The MAH determines mainly the halo structure and the color indexes while \lambda determines mainly the surface brightness and the bulge-to-disk ratio. We calculated star formation (SF) using a gravitational instability criterion and a self-regulation mechanism in the turbulent ISM. The efficiency of SF in this model is almost independent from the mass. We show that the luminosity-dependent dust absorption empirically determined by Wang & Heckman explains the observed color-magnitude and color Tully-Fisher (TF) relations without the necessity of introducing a mass-dependent SF efficiency. The disks in centrifugal equilibrium form within growing CDM halos with a gas accretion rate proportional to the MAH. The disks present exponential surface density and brightness profiles, negative radial color index gradients, and nearly flat rotation curves. We also calculated the secular formation of a bulge due to gravitational instabilities in the stellar disk. The intensive properties of our models agree with the observational data and the trends of the Hubble sequence are reproduced. The predicted infrared TF and luminosity-radius relations also agree with observations. The main shortcomings of our inside-out hierarchical models are the excessive radial color gradients and the dark halo dominion in the rotation curve decompositions.

We predict the structure and dynamics of disc galaxies using galaxy evolution models within a hierarchical formation scenario The halo mass aggregation histories, for a Lambda CDM model, were generated and used to calculate the virialization of dark matter (DM) haloes. A diversity of halo density profiles were obtained, the most typical one being close to the NFW profile. We modeled the formation of discs in centrifugal equilibrium within the evolving DM haloes using gas accretion rates proportional to the halo mass aggregation rates, and assuming detailed angular momentum conservation. We calculated the gravitational interactions between halo and disc, and the hydrodynamics, star formation, and evolution of the galaxy discs. We found that the slope and zero-point of the infrared Tully-Fisher relations (TFR) may be explained as a direct consequence of the cosmological initial conditions. This relation is almost independent of the assumed disc mass fraction. The rms scatter of the TFR originates mainly from the scatter in the DM halo structure and, to a minor extension, from the dispersion of the primordial spin parameter. The scatter obtained does not disagree with the observational estimates. Our models allow us to understand why the residuals of the TFR do not correlate significantly with disc size or surface brightness (SB), and why low and high SB galaxies have the same TFR. The correlations between gas fraction and SB, and between scale length and V_max agree with those observed. The discs present nearly exponential SB distributions. The shape of the rotation curves changes with the SB and is nearly flat for most cases. The rotation curve decompositions show a dominance of DM down to very small radii. The introduction of shallow cores in the DM halo attenuates this difficulty.

Using galaxy evolutionary models in a hierarchical formation scenario, we predict the structure, dynamics and evolution of disk galaxies in a LCDM universe. Our models include star formation and hydrodynamics of the ISM. We find that the Tully-Fisher relation (TFR) in the I and H bands is an imprint of the mass-velocity relation of the cosmological dark halos. The scatter of the TFR originates mainly from the scatter in the dark halo structure and, to a minor extension, from the dispersion of the primordial spin parameter lambda. Our models allow us to explain why low and high surface brightness galaxies have the same TFR. The disk gas fractions predicted agree with the observations. The disks formed within the growing halos have nearly exponential surface brightness and flat rotation curves. Towards high redshifts, the zero-point of the TFR in the H band increases while in the B-band it slightly decreases.

(Abridged) We study the outer density profiles of dark matter haloes predicted by a generalized secondary infall model and observed in a N-body cosmological simulation of a \Lambda CDM model. We find substantial systematic variations in shapes and concentrations of the halo profiles as well as a strong correlation of the profiles with the environment. In the N-body simulation, the average outer slope of the density profiles, \beta (\rho\propto r^{-\beta}), of isolated haloes is \approx 2.9; 68% of these haloes have values of \beta between 2.5 and 3.8. Haloes in dense environments of clusters are more concentrated and exhibit a broad distribution of \beta with values larger than for isolated haloes . Contrary to what one may expect, the haloes contained within groups and galaxy systems are less concentrated and have flatter outer density profiles than the isolated haloes. The concentration decreases with M_h, but its scatter for a given mass is substantial. The mass and circular velocity of the haloes are strongly correlated: M_h \propto V_m^{\alpha} with \alpha ~ 3.3 (isolated) and ~3.5 (haloes in clusters). For M_h=10^12M_sun the rms deviations from these relations are \Delta logM_h=0.12 and 0.18, respectively. Approximately 30% of the haloes are contained within larger haloes or have massive companions (larger than ~0.3 the mass of the current halo) within 3 virial radii. The remaining 70% of the haloes are isolated objects. The distribution of \beta as well as the concentration-mass and M_h-V_m relations for the isolated haloes agree very well with the predictions of our seminumerical approach which is based on a generalization of the secondary infall model and on the extended Press-Schechter formalism.

Making use of a seminumerical method we develop a scenario of disk galaxy formation and evolution in the framework of inflationary cold dark matter (CDM) cosmologies. Within the virializing dark matter halos, disks in centrifugal equilibrium are built-up and their galactic evolution is followed through an approach which considers the gravitational interactions among the galaxy components, the turbulence and energy balance of the ISM, the star formation (SF) process due to disk gravitational instabilities, the stellar evolution and the secular formation of a bulge. We find that the main properties and correlations of disk galaxies are determined by the mass, the hierarchical mass aggregation history and the primordial angular momentum. The models follow the same trends across the Hubble sequence than the observed galaxies. The predicted TF relation is in good agreement with the observations except for the standart CDM. While the slope of this relation remains almost constant up to intermediate redshifts, its zero-point decreases in the H-band and slightly increases in the B-band. A maximum in the SF rate for most of the models is attained at $z\sim 1.5-2.5$.

We study the evolution of disk galaxies within the frame of the cold dark matter (CDM) cosmologies. The hydrodynamics of a centrifugally supported gaseous disk and the growth of a stellar disk are calculated in detail taking into account the energy balance of the ISM and the gravitational instabilities that concern gas and stars. The halo density profile is derived from the primordial cosmological conditions and its gravitational contraction produced by the disk is included. Several features of the spiral galaxies at different redshifts are predicted, and the main factors which influence on these features are found. A strong evidence is provided that the Tully-Fisher relation is an imprint of the primordial cosmological conditions.

We study the formation of hot spheroidal systems within the frame of a scenario where galaxy formation and evolution is related to the gentle mass aggregation history and primordial angular momentum of protogalaxies, both defined by the cosmological initial conditions. We explore two cases: (1) the hot spheroidal system forms from the dynamical instabilities of the stellar disks, and (2) the spheroidal systems are formed during the dissipative collapse of the gas before falling to the disk in centrifugal equilibrium. In the former case a good agreement with observations for late type galaxies is found. In the second case, contrary to recent claims, we find that the tidal stability criterion is not easily reached. The gas that dissipatively collapses within the dark matter halos should be very clumpy, and the clumps very dense, in order to avoid the tidal destruction of the star formation unities.

We use a semianalytical approach and a CDM cosmological model to study the gravitational collapse and virialization, the structure, as well as the global and statistical properties of isolated dark matter galactic halos which emerge from primordial Gaussian fluctuations. From the statistical properties of the density fluctuation field the possible mass aggregation histories (MAHs) are generated, and these histories are used as the initial conditions of the gravitational collapse. To calculate the structure of the virialized systems we have generalized the secondary infall model. Although a range of halo structures is obtained, the average density profiles agree with the profile derived as a fitting formula to results of N-body cosmological simulations. The comparison of the density profiles with the observational data is disscused, and some possible solutions to the disagreement found in the inner regions are discussed. The results of our approach, after considering the gravitational dragging of the baryon matter that forms a central disk in centrifugal equilibrium, show that the Tully-Fisher relation and its scatter can be explained through the initial cosmological conditions. The sigma8=1 SCDM model produces galaxies with high velocities when compared to observations, but when the SCDM power spectrum is normalized to sigma8=0.57 an excellent agreement with the observable TF relation is found, suggesting that this relation is the natural extension to galactic scales of the observed galaxy distribution power spectrum. The theoretical TF scatter is close to the measured one. The slope of the TF relation is practically invariant with respect to the spin parameter lambda.

We present a scenario for the formation and evolution of disk galaxies within the framework of an inflationary CDM universe, and we compare the results with observations ranking from the present-day up to z~1. Galactic disks are built-up inside-out by gas infall with accretion rates driven by the cosmological mass aggregation history (MAH). We generate the MAHs for a Gaussian density fluctuation field, and we calculate the gravitational collapse and virialization of the dark halos. Assuming detailed angular momentum conservation, disks in centrifugal equilibrium are built-up within them. The disk galactic evolution is followed through a physically self-consistent approach. The main disk galaxy properties and their correlations are determined by the mass, the MAH, and the spin parameter. The models give exponential disk surface brightness (SB) profiles with realistic parameters (including LSB galaxies), nearly flat rotation curves, and negative gradients in the B-V radial profile. The main trends across the Hubble sequence are reproduced. A maximum in the star formation rate for most of the models is attained at z~ 1.5-2.5 (this rate is approximately 2.5-4.0 times larger than the present one). The scale radii and the bulge-to-total ratios decrease with z, while the central SB increase. The B-band TF relation remains almost the same at different redshifts. Our scenario reveals that the cosmological initial conditions are able to determine the main properties of disk galaxies and predict evolutionary features for the present-day dominant galaxy population that are in agreement with very recent deep field observational studies.

The galactic evolutionary theory is now sufficiently mature to provide information about the galaxy formation phase. From evolutionary models we establish a link between the present features of late type galaxies and the protogalactic density fluctuations. This link is used to estimate the fluctuation power spectrum amplitude at galactic scales. The agreement with the power spectrum derived from galaxy distribution and extrapolated to galactic scales, is satisfactory. Meanwhile, the comparison with respect to the standard cold dark matter power spectrum normalized to COBE satellite meassurements reveals only a marginal agreement. Among the most interesting results, the Tully-Fisher relation and the constant central brightness pointed out by Freeman (1970) appear to be the natural consequences of the initial cosmological conditions.

The formation and evolution of disk galaxies in the cosmological context is studied. We consider the observable properties of disk galaxies and treat the disk formation and galactic evolutionary processes in a self-consistent fashion. We find the matter accretion regime to be the dominant ingredient in establishing the Hubble sequence. The accretion regime is a phenomenon directly related to the statistical properties of the primordial density fluctuations from which disk galaxies emerged. In this paper we assumed a Gaussian density field and a standard CDM model.

The spherical gravitational collapse and virialization of arbitrary density fluctuations in an expanding universe is studied. In the context of the standard cosmological model and the peak $ansatz$, disk galaxies are supposed to be the product of "extended" gravitational collapses of regions surrounding local-maxima of the density fluctuation field. The substructure over the main profile (merging) is considered as a second-order phenomenon. The range and distribution of possible mass growth histories or accretion regimes for a given present-day mass are estimated and used as the initial condition for galaxy formation. Considering the effect of dark halo contraction due to disk formation, the final rotation curves are calculated. If the structural properties of dark halos are important in establishing the Hubble sequence's properties of visible galaxies, then this sequence is defined not only by the mass, but among other possible parameters, by the accretion regime.