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.

We study the nonlinear regime of large scale structure formation considering a dynamical dark energy (DE) component determined by a Steep Equation of State parametrization (SEoS) $w(z)=w_0+w_i\frac{(z/z_T)^q}{1+(z/z_T)^q}$. In order to perform the model exploration at low computational cost, we modified the public code L-PICOLA. We incorporate the DE model by means of the first and second-order matter perturbations in the Lagrangian frame and the expansion parameter. We analyze deviations of SEoS models with respect to $\Lambda$CDM in the non-linear matter power spectrum ($P_k$), the halo mass function (HMF), and the two-point correlation function (2PCF). On quantifying the nature of steep (SEoS-I) and smooth transitions in DE field (CPL-lim), no signature of steep transition is observed, rather found the overall impact of DE behaviors in $P_k$ at level of $\sim 2-3\%$ and $\sim 3-4\%$ differences w.r.t $\Lambda$CDM at $z=0$ respectively. HMF shows the possibility to distinguish between the models at the high mass ends. The best-fitted model assuming only background and linear perturbations dubbed as SEoS-II largely deviates from $\Lambda$CDM and current observations on studying the nonlinear growth. This large deviation in SEoS-II also quantified the combined effect of the dynamical DE and the larger amount of matter contained, $\Omega_{m0}$ and $H_{0}$ accordingly. 2PCF results are relatively robust with $\sim 1-2 \%$ deviation for SEoS-I and CPL-lim and a significant deviation for SEoS-II throughout $r$ from $\Lambda$CDM. Finally, we conclude that the search for viable DE models (like the SEoS) must include non-linear growth constraints.

Estimating the bar pattern speed (\Om{}) is one of the main challenges faced in understanding the role of stellar bars in galaxy dynamical evolution. This work aims to characterise different uncertainty sources affecting the Tremaine Weinberg (TW)-method to study the correlation between bar and galaxies physical parameters. We use a sample of 15 MaNGA SDSS-IV galaxies and 3 CALIFA galaxies from \cite{Aguerri2015}. We studied the errors related with (i) galaxy centre determination, (ii) disc position angle (PA) emphasising the difficulties triggered by outer non-axisymmetric structures besides the bar, (iii) the slits length and (iv) the spatial resolution. In average, the PA uncertainties range $\sim 15 \%$, the slit length $\sim 9 \%$ and the centring error $\sim 5 \%$. Reducing the spatial resolution increases the sensitivity to the PA error. Through Monte Carlo simulations, we estimate the probability distribution of the \R{} bar speed parameter. The present sample is composed of 7 slow, 4 fast and 7 ultrafast bars, with no trend with morphological types. Although uncertainties and low sample numbers may mask potential correlations between physical properties, we present a discussion of them: We observe an anti-correlation of \Om{} with the bar length and the stellar mass, suggesting that massive galaxies tend to host longer and slower bars. We also observe a correlation of the molecular gas fraction with \R{}, and a weak anti-correlation with \Om{}, suggesting that bars rotate slower in gaseous discs. Confirmation of such trends awaits future studies.

We explore two possible scenarios to explain the observed gamma-ray emission associated with the atypical globular cluster Omega-Centauri: emission from millisecond pulsars (MSP) and dark matter (DM) annihilation. In the first case the total number of MSPs needed to produce the gamma-ray flux is compatible with the known (but not confirmed) MSP candidates observed in X-rays. A DM interpretation is motivated by the possibility of Omega-Centauri being the remnant core of an ancient dwarf galaxy hosting a surviving DM component. At least two annihilation channels, light quarks and muons, can plausibly produce the observed gama-ray spectrum. We outline constraints on the parameter space of DM mass versus the product of the pair-annihilation cross section and integrated squared DM density (the so-called J-factor). We translate upper limits on the dark matter content of Omega-Centauri into lower limits on the annihilation cross section. This shows s-wave annihilation into muons to be inconsistent with CMB observations, while a small window for annihilation into light quarks is allowed. Further analysis of Omega-Centauri's internal kinematics, and/or additional information on the resident MSP population will yield much stronger constraints and shed light about the origin of this otherwise mysterious gamma-ray source.

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.

Recently, the Gaia data release 2 (DR2) showed us the richness in the kinematics of the Milky Way disk. Of particular interest is the presence of ridges covering the stellar velocity distribution, $V_{\phi}-R$; as shown by others, it is likely that these ridges are the signature of phase mixing, transient spirals, or the bar. Here, with a Galactic model containing both: bar and spirals, we found the same pattern of ridges extending from the inner to the outer disk. Interestingly, ridges in the $V_{\phi}-R$ plane correlate extremely well with wiggles in the computed rotation curve (RC). Hence, although the DR2 reveals (for the first time) such substructures in a wide spatial coverage, we notice that we have always seen such pattern of ridges, but projected into the form of wiggles in the RC. The separation and amplitude of the wiggles strongly depend on the extension and layout of ridges in the $V_{\phi}-R$ plane. This means that within the RC are encoded the kinematic state of the disk as well as information about the bar and spiral arms. The amplitude of the wiggles suggests that similar features currently observable in external galaxies RCs have similar origins, triggered by spirals and bars.

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.

We used ionized gas and stellar kinematics for 667 spatially resolved galaxies publicly available from the Calar Alto Legacy Integral Field Area survey (CALIFA) 3rd Data Release with the aim of studying kinematic scaling relations as the Tully $\&$ Fisher (TF) relation using rotation velocity, $V_{rot}$, the Faber $\&$ Jackson (FJ) relation using velocity dispersion, $\sigma$, and also a combination of $V_{rot}$ and $\sigma$ through the $S_{K}$ parameter defined as $S_{K}^2 = KV_{rot}^2 + \sigma^2$ with constant $K$. Late-type and early-type galaxies reproduce the TF and FJ relations. Some early-type galaxies also follow the TF relation and some late-type galaxies the FJ relation, but always with larger scatter. On the contrary, when we use the $S_{K}$ parameter, all galaxies, regardless of the morphological type, lie on the same scaling relation, showing a tight correlation with the total stellar mass, $M_\star$. Indeed, we find that the scatter in this relation is smaller or equal to that of the TF and FJ relations. We explore different values of the $K$ parameter without significant differences (slope and scatter) in our final results with respect the case $K=0.5$ besides than a small change in the zero point. We calibrate the kinematic $S_{K}^2$ dynamical mass proxy in order to make it consistent with sophisticated published dynamical models within $0.15\ dex$. We show that the $S_{K}$ proxy is able to reproduce the relation between the dynamical mass and the stellar mass in the inner regions of galaxies. Our result may be useful in order to produce fast estimations of the central dynamical mass in galaxies and to study correlations in large galaxy surveys.

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 study the radial acceleration relation (RAR) for early-type galaxies (ETGs) in the SDSS MaNGA MPL5 dataset. The complete ETG sample show a slightly offset RAR from the relation reported by McGaugh et al. (2016) at the low-acceleration end; we find that the deviation is due to the fact that the slow rotators show a systematically higher acceleration relation than the McGaugh's RAR, while the fast rotators show a consistent acceleration relation to McGaugh's RAR. There is a 1\sigma significant difference between the acceleration relations of the fast and slow rotators, suggesting that the acceleration relation correlates with the galactic spins, and that the slow rotators may have a different mass distribution compared with fast rotators and late-type galaxies. We suspect that the acceleration relation deviation of slow rotators may be attributed to more galaxy merger events, which would disrupt the original spins and correlated distributions of baryons and dark matter orbits in galaxies.

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/

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 report the peculiar chemical abundance patterns of eleven atypical Milky Way (MW) field red giant stars observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE). These atypical giants exhibit strong Al and N enhancements accompanied by C and Mg depletions, strikingly similar to those observed in the so-called second-generation (SG) stars of globular clusters (GCs). Remarkably, we find low-Mg abundances ([Mg/Fe]$<$0.0) together with strong Al and N overabundances in the majority (5/7) of the metal-rich ([Fe/H]$\gtrsim - 1.0$) sample stars, which is at odds with actual observations of SG stars in Galactic CGs of similar metallicities. This chemical pattern is unique and unprecedented among MW stars, posing urgent questions about its origin. These atypical stars could be former SG stars of dissolved GCs formed with intrinsically lower abundances of Mg and enriched Al (subsequently self-polluted by massive AGB stars) or the result of exotic binary systems. We speculate that the stars Mg-deficiency as well as the orbital properties suggest that they could have an extragalactic origin. This discovery should guide future dedicated spectroscopic searches of atypical stellar chemical patterns in our Galaxy; a fundamental step forward to understand the Galactic formation and evolution.

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.

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.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

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.

We report on detection, from observations obtained with the APOGEE spectroscopic survey, of a metal-poor ([Fe/H] $= -1.3$ dex) field giant star with an extreme Mg-Al abundance ratio ([Mg/Fe] $= -0.31$ dex; [Al/Fe] $= 1.49$ dex). Such low Mg/Al ratios are seen only among the second-generation population of globular clusters, and are not present among Galactic disk field stars. The light element abundances of this star, 2M16011638-1201525, suggest that it could have been born in a globular cluster. We explore several origin scenarios, in particular studying the orbit of the star to check the probability of it being kinematically related to known globular clusters. We performed simple orbital integrations assuming the estimated distance of 2M16011638-1201525 and the available six-dimensional phase-space coordinates of 63 globular clusters, looking for close encounters in the past with a minimum distance approach within the tidal radius of each cluster. We found a very low probability that 2M16011638-1201525 was ejected from most globular clusters; however, we note that the best progenitor candidate to host this star is globular cluster $\omega$ Centauri (NGC 5139). Our dynamical investigation demonstrates that 2M16011638-1201525 reaches a distance $|Z_{max}| < 3 $ kpc from the Galactic plane and a minimum and maximum approach to the Galactic center of $R_{min}<0.62$ kpc and $R_{max}<7.26$ kpc in an eccentric ($e\sim0.53$) and retrograde orbit. Since the extreme chemical anomaly of 2M16011638-1201525 has also been observed in halo field stars, this object could also be considered a halo contaminant, likely been ejected into the Milky Way disk from the halo. We conclude that, 2M16011638-20152 is also kinematically consistent with the disk but chemically consistent with halo field stars.

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.

The most accurate 6D phase-space information from the Radial Velocity Experiment (RAVE) was used to integrate the orbits of 105 stars around the galactic globular cluster 47 Tucanae, to look for close encounters between them in the past, with a minimum distance approach less than the cluster tidal radius. The stars are currently over the distance range 3.0 kpc $<$ d $<$ 5.5 kpc. Using the uncertainties in the current position and velocity vector for both, star and cluster, 105 pairs of star-cluster orbits were generated in a Monte Carlo numerical scheme, integrated over 2 Gyr and considering an axisymmetric and non-axisymmetric Milky-Way-like Galactic potential, respectively. In this scheme, we identified 20 potential cluster members that had close encounters with the globular cluster 47 Tucanae, all of which have a relative velocity distribution (V$_{rel}$) less than 200 km s$^{-1}$ at the minimum distance approach. Among these potential members, 9 had close encounters with the cluster with velocities less than the escape velocity of 47 Tucanae, therefore a scenario of tidal stripping seems likely. These stars have been classified with a 93\% confidence level, leading to the identification of extratidal cluster stars. For the other 11 stars, V$_{rel}$ exceeds the escape velocity of the cluster, therefore they were likely ejected or are unassociated interlopers.

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.

Using six-dimesional phase-space information from the Fourth Data release of the Radial Velocity Experiment (RAVE) over the range of Galactic longitude 240$^{\circ}< l <$ 360$^{\circ}$ and $V_{LSR} < -239$ kms$^{-1}$, we have computed orbits for 329 RAVE stars that were originally selected as chemically and kinematically related to $\omega$ Centauri. The orbits were integrated in a Milky-Way-like axisymmetric Galactic potential, ignoring the effects of the dynamical evolution of $\omega$ Centauri due to the tidal effects of the Galaxy disk on the cluster along time. We also ignored secular changes in the Milky Way potential over time. In a Monte Carlo scheme, and under the assumption that the stars may have been ejected with velocities greater than the escape velocity ($V_{rel}>V_{esc,0}$) from the cluster, we identified 15 stars as having close encounters with $\omega$ Centauri: (\textit{i}) 8 stars with relative velocities $V_{rel}< 200 $ kms$^{-1}$ may have been ejected $\sim$ 200 Myr ago from $\omega$ Centauri; (\textit{ii}) other group of 7 stars were identified with high relative velocity $V_{rel}> 200 $ kms$^{-1}$ during close encounters, and seems unlikely that they have been ejected from $\omega$ Centauri. We also confirm the link between J131340.4-484714 as potential member of $\omega$ Centauri, and probably ejected $\sim$ 2.0 Myr ago, with a relative velocity $V_{rel}\sim80$ kms$^{-1}$.

We present a method to identify Ultra Faint Dwarf Galaxy (UFDG) candidates in the halo of the Milky Way using the future Gaia catalogue and we explore its detection limits and completeness. The method is based on the Wavelet Transform and searches for over-densities in the combined space of sky coordinates and proper motions, using kinematics in the search for the first time. We test the method with a Gaia mock catalogue that has the Gaia Universe Model Snapshot (GUMS) as a background, and use a library of around 30 000 UFDGs simulated as Plummer spheres with a single stellar population. For the UFDGs we use a wide range of structural and orbital parameters that go beyond the range spanned by real systems, where some UFDGs may remain undetected. We characterize the detection limits as function of the number of observable stars by Gaia in the UFDGs with respect to that of the background and their apparent sizes in the sky and proper motion planes. We find that the addition of proper motions in the search improves considerably the detections compared to a photometric survey at the same magnitude limit. Our experiments suggest that Gaia will be able to detect UFDGs that are similar to some of the known UFDGs even if the limit of Gaia is around 2 magnitudes brighter than that of SDSS, with the advantage of having a full-sky catalogue. We also see that Gaia could even find some UFDGs that have lower surface brightness than the SDSS limit.

Recent studies have presented evidence that the Milky Way global potential may be nonspherical. In this case, the assembling process of the Galaxy may have left long lasting stellar halo kinematic fossils due to the shape of the dark matter halo, potentially originated by orbital resonances. We further investigate such possibility, considering now potential models further away from $\Lambda$CDM halos, like scalar field dark matter halos, MOND, and including several other factors that may mimic the emergence and permanence of kinematic groups, such as, a spherical and triaxial halo with an embedded disk potential. We find that regardless of the density profile (DM nature), kinematic groups only appear in the presence of a triaxial halo potential. For the case of a MOND like gravity theory no kinematic structure is present. We conclude that the detection of these kinematic stellar groups could confirm the predicted triaxiality of dark halos in cosmological galaxy formation scenarios.

We present a survey for RR Lyrae stars in an area of 50 deg$^2$ around the globular cluster $\omega$ Centauri, aimed to detect debris material from the alleged progenitor galaxy of the cluster. We detected 48 RR Lyrae stars of which only 11 have been previously reported. Ten among the eleven previously known stars were found inside the tidal radius of the cluster. The rest were located outside the tidal radius up to distances of $\sim 6$ degrees from the center of the cluster. Several of those stars are located at distances similar to that of $\omega$ Centauri. We investigated the probability that those stars may have been stripped off the cluster by studying their properties (mean periods), calculating the expected halo/thick disk population of RR Lyrae stars in this part of the sky, analyzing the radial velocity of a sub-sample of the RR Lyrae stars, and finally, studying the probable orbits of this sub-sample around the Galaxy. None of these investigations support the scenario that there is significant tidal debris around $\omega$ Centauri, confirming previous studies in the region. It is puzzling that tidal debris have been found elsewhere but not near the cluster itself.

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

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

We map the kinematics of stars in simulated galaxy disks with spiral arms using the velocity ellipsoid vertex deviation (l$_v$). We use test particle simulations, and for the first time, fully self-consistent high resolution N-body models. We compare our maps with the Tight Winding Approximation model analytical predictions. We see that for all barred models spiral arms rotate closely to a rigid body manner and the vertex deviation values correlate with the density peaks position bounded by overdense and underdense regions. In such cases, vertex deviation sign changes from negative to positive when crossing the spiral arms in the direction of disk rotation, in regions where the spiral arms are in between corotation (CR) and the Outer Lindblad Resonance (OLR). By contrast, when the arm sections are inside the CR and outside the OLR, l$_v$ changes from negative to positive.We propose that measurements of the vertex deviations pattern can be used to trace the position of the main resonances of the spiral arms. We propose that this technique might exploit future data from Gaia and APOGEE surveys. For unbarred N-body simulations with spiral arms corotating with disk material at all radii, our analysis suggests that no clear correlation exists between l$_v$ and density structures.

We present the result of an extensive search for non-axisymmetric structures in the dwarf satellite galaxy of M81: NGC 2976, using multiwavelength archival observations. The galaxy is known to present kinematic evidence for a bysimmetric distortion, anyhow the stellar bar presence is controversial, due to the possible interpretation of NGC 2976 presenting an elliptical disk triggered by a prolate dark matter halo. We applied diagnostics in order to detect stellar bars or spiral arms. The m=2 fourier phase has a jump around 60 arcsecs consistent with a central bar and bisymmetric arms. The CO, 3.6 $\mu$ surface brightness and the dust lanes are consistent with a gas rich central bar and possibly with gaseous spiral arms. The kinematic jumps related with the dust lanes suggest that the bar perturbation in the disk kinematics is non-negligible and the reported non-circular motions, the central gas excess and the nuclear X-ray source (AGN/Starburst) might be produced by the central bar. SPH simulations of disks inside triaxial dark halos suggest that the two symmetric spots at 130 arcsecs and the narrow arms may be produced by gas at turning points in an elliptical disk, alternatively the potential ellipticity can be originated by tidally induced strong stellar bar/arms. We conclude that, the gas rich bar and spiral arms triggered by tidal stirring, and primordial halo triaxiality, can explain most of the non-circular motions, mass redistribution and nuclear activity. Distinguishing between them requires detailed modeling of environmental effects. An analysis similar to this may reveal such kind of structures in other nearby dwarf satellite galaxies, and if confirmed, the same evolutionary scenario will be applicable to them. This implies biases constraining their dark matter distribution and also making comparison against theoretical predictions for isolated galaxies.

We present a set of hydrodynamical/Nbody controlled simulations of isolated gas rich galaxies that self-consistently include SN feedback and a detailed chemical evolution model, both tested in cosmological simulations. The initial conditions are motivated by the observed star forming galaxies at z ~ 2-3. We find that the presence of a multiphase interstellar media in our models promotes the growth of disc instability favouring the formation of clumps which in general, are not easily disrupted on timescales compared to the migration time. We show that stellar clumps migrate towards the central region and contribute to form a classical-like bulge with a Sersic index, n > 2. Our physically-motivated Supernova feedback has a mild influence on clump survival and evolution, partially limiting the mass growth of clumps as the energy released per Supernova event is increased, with the consequent flattening of the bulge profile. This regulation does not prevent the building of a classical-like bulge even for the most energetic feedback tested. Our Supernova feedback model is able to establish a self-regulated star formation, producing mass-loaded outflows and stellar age spreads comparable to observations. We find that the bulge formation by clumps may coexit with other channels of bulge assembly such as bar and mergers. Our results suggest that galactic bulges could be interpreted as composite systems with structural components and stellar populations storing archaeological information of the dynamical history of their galaxy.

High resolution N-body simulations using different codes and initial condition techniques reveal two different behaviours for the rotation frequency of transient spiral arms like structures. Whereas unbarred disks present spiral arms nearly corotatingwith disk particles, strong barred models (bulged or bulge-less) quickly develop a bar-spiral structure dominant in density, with a pattern speed almost constant in radius. As the bar strength decreases the arm departs from bar rigid rotation and behaves similar to the unbarred case. In strong barred models we detect in the frequency space other subdominant and slower modes at large radii, in agreement with previous studies, however we also detect them in the configuration space. We propose that the distinctive behaviour of the dominant spiral modes can be exploited in order to constraint the nature of Galactic spiral arms by the astrometric survey GAIA and by 2-D spectroscopic surveys like CALIFA and MANGA in external galaxies.

The accuracy of the measurements of some astrophysical dynamical systems allows to constrain the existence of incredibly small gravitational perturbations. In particular, the internal Solar System dynamics (planets, Earth-Moon) opens up the possibility, for the first time, to prove the abundance, mass and size, of dark sub-structures at the Earth vicinity. We find that adopting the standard dark matter density, its local distribution can be composed by sub-solar mass halos with no currently measurable dynamical consequences, regardless of the mini-halo fraction. On the other hand, it is possible to exclude the presence of dark streams with linear mass densities higher than $\lambda_{\rm st}> 10^{-10} \Msun/\AU$ (about the Earth mass spread along the diameter of the SS up to the Kuiper belt). In addition, we review the dynamics of wide binaries inside the dwarf spheroidal galaxies in the MW. The dynamics of such kind of binaries seem to be compatible with the presence of a huge fraction of dark sub-structure, thus their existence is not a sharp discriminant of the dark matter hypothesis as been claimed before. However, there are regimes where the constraints from different astrophysical systems may reveal the sub-structure mass function cut-off scale.

A massive, self-interacting scalar field has been considered as a possible candidate for the dark matter in the universe. We present an observational constraint to the model arising from strong lensing observations in galaxies. The result points to a discrepancy in the properties of scalar field dark matter halos for dwarf and lens galaxies, mainly because halo parameters are directly related to physical quantities in the model. This is an important indication that it becomes necessary to have a better understanding of halo evolution in scalar field dark matter models, where the presence of baryons can play an important role.

Assuming the dark matter halo of the Milky Way as a non-spherical potential (i.e. triaxial, prolate, oblate), we show how the assembling process of the Milky Way halo, may have left long lasting stellar halo kinematic fossils only due to the shape of the dark matter halo. In contrast with tidal streams, associated with recent satellite accretion events, these stellar kinematic groups will typically show inhomogeneous chemical and stellar population properties. However, they may be dominated by a single accretion event for certain mass assembling histories. If the detection of these peculiar kinematic stellar groups is confirmed, they would be the smoking gun for the predicted triaxiality of dark halos in cosmological galaxy formation scenarios.

We study the spiral arm influence on the solar neighbourhood stellar kinematics. As the nature of the Milky Way (MW) spiral arms is not completely determined, we study two models: the Tight-Winding Approximation (TWA) model, which represents a local approximation, and a model with self-consistent material arms named PERLAS. This is a mass distribution with more abrupt gravitational forces. We perform test particle simulations after tuning the two models to the observational range for the MW spiral arm properties. We explore the effects of the arm properties and find that a significant region of the allowed parameter space favours the appearance of kinematic groups. The velocity distribution is mostly sensitive to the relative spiral arm phase and pattern speed. In all cases the arms induce strong kinematic imprints for pattern speeds around 17 km/s/kpc (close to the 4:1 inner resonance) but no substructure is induced close to corotation. The groups change significantly if one moves only ~0.6 kpc in galactocentric radius, but ~2 kpc in azimuth. The appearance time of each group is different, ranging from 0 to more than 1 Gyr. Recent spiral arms can produce strong kinematic structures. The stellar response to the two potential models is significantly different near the Sun, both in density and kinematics. The PERLAS model triggers more substructure for a larger range of pattern speed values. The kinematic groups can be used to reduce the current uncertainty about the MW spiral structure and to test whether this follows the TWA. However, groups such as the observed ones in the solar vicinity can be reproduced by different parameter combinations. Data from velocity distributions at larger distances are needed for a definitive constraint.

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

NGC3367 is a nearby isolated active galaxy that shows a radio jet, a strong bar and evidence of lopsidedness. We present a quantitative analysis of the stellar and gaseous structure of the galaxy disk and a search for evidence of recent interaction based on new UBVRI Halpha and JHK images and on archival Halpha Fabry-Perot and HI VLA data. From a coupled 1D/2D GALFIT bulge/bar/disk decomposition an (B/D ~ 0.07-0.1) exponential pseudobulge is inferred in all the observed bands. A NIR estimate of the bar strength <Q_T{max}(R)> = 0.44 places NGC 3367 bar among the strongest ones. The asymmetry properties were studied using (1) optical and NIR CAS indexes (2) the stellar (NIR) and gaseous (Halpha, HI) A_1 Fourier mode amplitudes and (3) the HI integrated profile and HI mean intensity distribution. While the average stellar component shows asymmetry values close to the average found in the Local Universe for isolated galaxies, the young stellar component and gas values are largely decoupled showing significantly larger A_1 mode amplitudes suggesting that the gas has been recently perturbed. Our search for (1) faint stellar structures in the outer regions (up to u_R ~ 26 mag arcsec^{-2}), (2) (Halpha) star-forming satellite galaxies and (3) regions with different colors (stellar populations) along the disk all failed. Such an absence is interpreted using recent numerical simulations to constrain a tidal event with an LMC like galaxy to some dynamical times in the past or to a current very low mass, gas rich accretion. We conclude that a cold accretion mode (gas and small/dark galaxies) may be responsible of the nuclear activity and peculiar (young stars and gas) morphology regardless of the highly isolated environment. Black hole growth in bulgeless galaxies may be triggered by cosmic smooth mass accretion.

We propose an extension of the GC3 streamer finding method of Johnston et al. 1996 that can be applied to the future Gaia database. The original method looks for streamers along great circles in the sky, our extension adds the kinematical restriction that velocity vectors should also be constrained to lie along these great circles, as seen by a Galactocentric observer. We show how to use these combined criteria starting from heliocentric observables. We test it by using the mock Gaia catalogue of Brown et al. 2005, which includes a realistic Galactic background and observational errors, but with the addition of detailed star formation histories for the simulated satellites. We investigate its success rate as a function of initial satellite luminosity, star formation history and orbit. We find that the inclusion of the kinematical restriction vastly enhances the contrast between a streamer and the background, even in the presence of observational errors, provided we use only data with good astrometric quality (fractional errors of 30 per cent or better). The global nature of the method diminishes the erasing effect of phase mixing and permits the recovery of merger events of reasonable dynamical age. Satellites with a star formation history different to that of the Galactic background are also better isolated. We find that satellites in the range of 10^8-10^9 Lsun can be recovered even for events as old as ~10 Gyr. Even satellites with 4-5x10^7 Lsun can be recovered for certain combinations of dynamical ages and orbits.

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.

We study the phase space available to the local stellar distribution using a Galactic potential consistent with several recent observational constraints. We find that the induced phase space structure has several observable consequences. The spiral arm contribution to the kinematic structure in the solar neighborhood may be as important as the one produced by the Galactic bar. We suggest that some of the stellar kinematic groups in the solar neighborhood, like the Hercules structure and the kinematic branches, can be created by the dynamical resonances of self-gravitating spiral arms and not exclusively by the Galactic bar. A structure coincident with the Arcturus kinematic group is developed when a hot stellar disk population is considered, which introduces a new perspective on the interpretation of its extragalactic origin. A bar-related resonant mechanism can modify this kinematic structure.We show that particles in the dark matter disk-like structure predicted by recent LCDM galaxy formation experiments, with similar kinematics to the thick disk, are affected by the same resonances, developing phase space structures or dark kinematic groups that are independent of the Galaxy assembly history and substructure abundance. We discuss the possibility of using the stellar phase space groups as constraints to non-axisymmetric models of the Milky Way structure.

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.

We study dynamics of bars in models of disk galaxies embeded in realistic dark matter halos. We find that disk thickness plays an important, if not dominant, role in the evolution and structure of the bars. We also make extensive numerical tests of different N-body codes used to study bar dynamics. Models with thick disks typically used in this type of modeling (height-to-length ratio hz/Rd=0.2) produce slowly rotating, and very long, bars. In contrast, more realistic thin disks with the same parameters as in our Galaxy (hz/Rd= 0.1) produce bars with normal length Rbar approx R_d, which rotate quickly with the ratio of the corotation radius to the bar radius 1.2-1.4 compatible with observations. Bars in these models do not show a tendency to slow down, and may lose as little as 2-3 percent of their angular momentum due to dynamical friction with the dark matter over cosmological time. We attribute the differences between the models to a combined effect of high phase-space density and smaller Jeans mass in the thin disk models, which result in the formation of a dense central bulge. Special attention is paid to numerical effects such as the accuracy of orbital integration, force and mass resolution. Using three N-body codes -- Gadget, ART, and Pkdgrav -- we find that numerical effects are very important and, if not carefully treated, may produce incorrect and misleading results. Once the simulations are performed with sufficiently small time-steps and with adequate force and mass resolution, all the codes produce nearly the same results: we do not find any systematic deviations between the results obtained with TREE codes (Gadget and Pkdgrav) and with the Adaptive-Mesh-Refinement (ART) code.

We present, using a novel technique, a study of the angular distribution of satellite galaxies around a sample of isolated, blue host galaxies selected from the sixth data release of the Sloan Digital Sky Survey. As a complement to previous studies we subdivide the sample of galaxies into bins of differing inclination and use the systematic differences that would exist between the different bins as the basis for our approach. We parameterize the cumulative distribution function of satellite galaxies and apply a maximum likelihood, Monte-Carlo technique to determine allowable distributions, which we show as an exclusion plot. We find that the allowed distributions of the satellites of spiral hosts are very nearly isotropic. We outline our formalism and our analysis and discuss how this technique may be refined for future studies and future surveys.

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 used fully cosmological, high resolution N-body + SPH simulations to follow the formation of disk galaxies with rotational velocities between 135 and 270 km/sec in a Lambda CDM universe. The simulations include gas cooling, star formation, the effects of a uniform UV background and a physically motivated description of feedback from supernovae. The host dark matter halos have a spin and last major merger redshift typical of galaxy sized halos as measured in recent large scale N--Body simulations. The simulated galaxies form rotationally supported disks with realistic exponential scale lengths and fall on both the I-band and baryonic Tully Fisher relations. An extended stellar disk forms inside the Milky Way sized halo immediately after the last major merger. The combination of UV background and SN feedback drastically reduces the number of visible satellites orbiting inside a Milky Way sized halo, bringing it in fair agreement with observations. Our simulations predict that the average age of a primary galaxy's stellar population decreases with mass, because feedback delays star formation in less massive galaxies. Galaxies have stellar masses and current star formation rates as a function of total mass that are in good agreement with observational data. We discuss how both high mass and force resolution and a realistic description of star formation and feedback are important ingredients to match the observed properties of galaxies.

Two well studied dwarf galaxies -- NGC 3109 and NGC 6822 -- present some of the strongest observational support for a flat core at the center of galactic dark matter (DM) halos. We use detailed cosmologically motivated numerical models to investigate the systematics and the accuracy of recovering parameters of the galaxies. Some of our models match the observed structure of the two galaxies remarkably well. Our analysis shows that the rotation curves of these two galaxies are instead quite compatible with their DM halos having steep cuspy density profiles. The rotation curves in our models are measured using standard observational techniques. The models reproduce the rotation curves of both galaxies, the disk surface brightness profiles as well as the profile of isophotal ellipticity and position angle. The models are centrally dominated by baryons; however, the dark matter component is globally dominant. The simulated disk mass is marginally consistent with a stellar mass-to-light ratio in agreement with the observed colors. We show that non-circular motions combined with gas pressure support and projection effects results in a large underestimation of the circular velocity in the central $\sim 1$ kpc region, creating the illusion of a constant density core. Although the systematic effects mentioned above are stronger in barred systems, they are also present in axisymetric disks. Our results strongly suggest that there is no contradiction between the observed rotation curves in dwarf galaxies and the cuspy central dark matter density profiles predicted by Cold Dark Matter models.

The central part of a dark matter halo reacts to the presence and evolution of a bar. Not only does the halo absorb angular momentum from the disk, it can also be compressed and have its shape modified. We study these issues in a series of cosmologically motivated, highly resolved N-body simulations of barred galaxies run under different initial conditions. In all models we find that the inner halo's central density increases. We model this density increase using the standard adiabatic approximation and the modified formula by Gnedin et al. and find that halo mass profiles are better reproduced by this latter. In models with a strong bar, the dark matter in the central region forms a bar-like structure (``dark matter bar''), which rotates together with the normal bar formed by the stellar component (``stellar bar''). The minor-to-major axial ratio of a halo bar changes with radius with a typical value 0.7 in the central disk region. DM bar amplitude is mostly a function of the stellar bar strength. Models in which the bar amplitude increases or stays roughly constant with time, initially large (40%-60%) misalignment between the halo and disk bars quickly decreases with time as the bar grows. The halo bar is nearly aligned with the stellar bar (~10 degrees lag for the halo) after ~2 Gyr. The torque, which the halo bar exerts on the stellar bar, can serve as a mechanism to regulate the angular momentum transfer from the disk to the halo.

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 address the issue of accuracy in recovering density profiles from observations of rotation curves of galaxies. We ``observe'' and analyze our models in much the same way as observers do the real galaxies. We find that the tilted ring model analysis produces an underestimate of the central rotational velocity. In some cases the galaxy halo density profile seems to have a flat core, while in reality it does not. We identify three effects, which explain the systematic biases: (1) inclination (2), small bulge, and (3) bar. The presence of even a small non-rotating bulge component reduces the rotation velocity. In the case of a disk with a bar, the underestimate of the circular velocity is larger due to a combination of non-circular motions and random velocities. Signatures of bars can be difficult to detect in the surface brightness profiles of the model galaxies. The variation of inclination angle and isophote position angle with radius are more reliable indicators of bar presence than the surface brightness profiles. The systematic biases in the central ~ 1 kpc of galaxies are not large. Each effect separately gives typically a few kms error, but the effects add up. In some cases the error in circular velocity was a factor of two, but typically we get about 20 percent. The result is the false inference that the density profile of the halo flattens in the central parts. Our observations of real galaxies show that for a large fraction of galaxies the velocity of gas rotation (as measured by emission lines) is very close to the rotation of stellar component (as measured by absorption lines). This implies that the systematic effects discussed in this paper are also applicable both for the stars and emission-line gas.

We study properties of dark matter halos at high redshifts z=2-10 for a vast range of masses with the emphasis on dwarf halos with masses 10^7-10^9 Msun/h. We find that the density profiles of relaxed dwarf halos are well fitted by the NFW profile and do not have cores. We compute the halo mass function and the halo spin parameter distribution and find that the former is very well reproduced by the Sheth & Tormen model while the latter is well fitted by a lognormal distribution with lambda_0 = 0.042 and sigma_lambda = 0.63. We estimate the distribution of concentrations for halos in mass range that covers six orders of magnitude from 10^7 Msun/h to 10^13} Msun/h, and find that the data are well reproduced by the model of Bullock et al. The extrapolation of our results to z = 0 predicts that present-day isolated dwarf halos should have a very large median concentration of ~ 35. We measure the subhalo circular velocity functions for halos with masses that range from 4.6 x 10^9 Msun/h to 10^13 Msun/h and find that they are similar when normalized to the circular velocity of the parent halo. Dwarf halos studied in this paper are many orders of magnitude smaller than well-studied cluster- and Milky Way-sized halos. Yet, in all respects the dwarfs are just down-scaled versions of the large halos. They are cuspy and, as expected, more concentrated. They have the same spin parameter distribution and follow the same mass function that was measured for large halos.

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.

Using high resolution N-body simulations of stellar disks embedded in cosmologically motivated dark matter halos, we study the evolution of bars and the transfer of angular momentum between halos and bars. We find that dynamical friction results in some transfer of angular momentum to the halo, but the effect is much smaller than previously found in low resolution simulations and is incompatible with early analytical estimates. In simulations with millions of particles reaching a resolution of 20-40 pc, the pattern speed may not change over billions of years. Our high resolution models produce bars which are fast rotators, where the ratio of the corotation radius to the bar major semi-axis lies in the range R =1.2-1.7, marginally compatible with observational results. In contrast to many previous simulations, we find that bars are relatively short. As in many observed cases, the bar major semi-axis is close to the exponential length of the disk. The transfer of angular momentum between inner and outer parts of the disk plays a very important role in the secular evolution of the disk and the bar. The bar formation increases the exponential length of the disk by a factor of 1.2 -1.5. The transfer substantially increases the stellar mass in the centre of the galaxy and decreases the dark matter-to-baryons ratio. As the result, the central 2 kpc region is always strongly dominated by the baryonic component. At intermediate 3-10kpc scales the disk is sub-dominant. In summary, realistic models produce bar structure in striking agreement with observational results.

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.

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

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

Using published data, we have compiled the circular velocity (Vc) distribution function (VDF) of galaxy satellites in the Local Group. We find that within the volumes of radius of 570 kpc (400/h kpc for h=0.7) centered on the Milky Way and Andromeda, the average VDF is roughly approximated as n(>Vc)~ 45(Vc/10 km/s)^{-1} h^3 Mpc^{-3} for Vc in the range ~10-70 km/s. The observed VDF is compared with results of high-resolution cosmological simulations. We find that the VDF in models is very different from the observed one: n(>Vc)~1200(Vc/10 km/s)^{-2.75}h^3 Mpc^{-3}. Cosmological models thus predict that a halo of the size of our Galaxy should have about 50 dark matter satellites with circular velocity >20 km/s and mass >3x10^8/h Msun within a 570 kpc radius. This number is significantly higher than the approximate dozen satellites actually observed around our Galaxy. The observed and predicted VDFs cross at ~50 km/s, indicating that the predicted abundance of satellites with Vc> 50 km/s is in reasonably good agreement with observations. We conclude, therefore, that unless a large fraction of the Local Group satellites has been missed in observations, there is a dramatic discrepancy between observations and hierarchical models, regardless of the model parameters. We discuss several possible explanations for this discrepancy including identification of some satellites with the High Velocity Clouds observed in the Local Group, and the existence of dark satellites that failed to accrete gas and form stars due either to the expulsion of gas in the supernovae-driven winds or to gas heating by the intergalactic ionizing background. (Abridged)