Source author record

O. Valenzuela

O. Valenzuela appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2022arXiv

Assessing the Hierarchical Hamiltonian Splitting Integrator for Collisionless N-body Simulations

The large dynamic range in some astrophysical N-body problems led to the use of adaptive multi-time-steps; however, the search for optimal strategies is still challenging. We numerically quantify the performance of the hierarchical Hamiltonian Splitting (HHS) integrator for collisionless simulations using a direct summation code. We compare HHS with the constant global time-step leapfrog integrator, and with the adaptive one (AKDK). We find that HHS is approximately reversible, whereas AKDK not. Therefore, it is possible to find a combination of parameters where the energy drift is considerably milder for HHS, resulting in a better performance. We conclude that HHS is an attractive alternative to AKDK, and it is certainly advantageous for direct summation and P3M codes. Also, we find advantages with GADGET4 (Tree/FMM) HHS implementation that are worth exploring further.

preprint2022arXiv

SDSS-IV MaNGA: The MaNGA Dwarf Galaxy Sample

We present the MaNGA Dwarf galaxy, MaNDala, Value-Added-Catalog, VAC, from the final release of the Sloan Digital Sky Survey-IV program. MaNDala consists of 136 randomly selected bright dwarf galaxies with $M_{*}<10^{9.1}\odot$ and $M_{g}>-18.5$ making it the largest Integral Field Spectroscopy homogeneous sample of dwarf galaxies. We release a photometric analysis of the $g,r$ and $z$ broadband imaging based on the DESI Legacy Imaging Surveys as well as the spectroscopic analysis based on the Pipe3D SDSS-IV VAC. Our release includes the surface brightness (SB), geometric parameters and color profiles, Sérsic fits as well as stellar population properties (such as, stellar ages, metallicities, star formation histories), and emission lines fluxes within the FOV and the effective radii of the galaxies. We find that the majority of the MaNDala galaxies are star forming late-type galaxies with $\langle{}n_{\text{Sersic,r}}\rangle\sim1.6$ that are centrals (central/satellite dichotomy). MaNDala covers a large range of SB values (we find 11 candidates of ultra diffuse galaxies and 3 compact ones), filling the gap between classical dwarfs and low-mass galaxies in the Kormendy Diagram and in the size-mass/luminosity relation, whichseems to flatten at $10^8<M_{*}/\odot<10^{9}$ with $\langle{}R_{e,r}\rangle\sim2.7$ kpc. A large fraction of MaNDala galaxies formed from an early low-metallicity burst of SF but also of late SF events from more metal-enriched gas: half of the MaNDala galaxies assembled $50\%$ of their mass at $\langle{}z\rangle>2$, while the last $20\%$ was at $\langle{}z\rangle<0.3$. Finally, a bending of the sSFR-$M_{*}$ relation at $M_{*}\sim10^{9}\odot$ for the main sequence galaxies seems to be supported by MaNDala.

preprint2022arXiv

The Lockman-SpReSO project. Description, target selection, observations and preliminary results

Context. Extragalactic surveys are a key tool for better understanding the evolution of galaxies. Both deep and wide-field surveys serve to provide a clearer emerging picture of the physical processes that take place in and around galaxies, and to identify which of these processes are the most important in shaping the properties of galaxies. Aims. The Lockman Spectroscopic Redshift Survey using Osiris (Lockman-SpReSO) aims to provide one of the most complete optical spectroscopic follow-ups of the far-infrared (FIR) sources detected by the \textit{Herschel} Space Observatory in the Lockman Hole (LH) field. The optical spectroscopic study of the FIR-selected galaxies supplies valuable information about the relation between fundamental FIR and optical parameters, including extinction, star formation rate, and gas metallicity. In this article, we introduce and provide an in-depth description of the Lockman-SpReSO project and of its early results. Methods. We selected FIR sources from \textit{Herschel} observations of the central 24 arcmin $\times$ 24 arcmin of the LH field with an optical counterpart up to 24.5 $R_{\rm C}$(AB). The sample comprises 956 \textit{Herschel} FIR sources, plus 188 additional interesting objects in the field. These are point X-ray sources, cataclysmic variable star candidates, high-velocity halo star candidates, radio sources, very red quasi-stellar objects, and optical counterparts of sub-millimetre galaxies. The faint component of the catalogue ($R_{\rm C}(\mathrm{AB})\geq20$) was observed using the OSIRIS instrument on the 10.4 m Gran Telescopio Canarias in multi-object spectroscopy (MOS) mode. The bright component was observed using two multi-fibre spectrographs: the AF2-WYFFOS at the William Herschel Telescope and the HYDRA instrument at the WYIN telescope.

preprint2020arXiv

A Universal fundamental plane and the $M_{dyn}-M_{\star}$ relation for galaxies with CALIFA and MaNGA

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}+σ_{\star_e}^{2}$, which combines rotation velocity $V_{R_e}$, and velocity dispersion $σ_{\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}-σ_{\star_e}$ ones. We also explore a more general Universal Fundamental Plane described by the equation $log(Υ_{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, $Υ_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, $Υ_{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.

preprint2019arXiv

On the origin of the gamma-ray emission from Omega Centauri: Milisecond pulsars and dark matter annihilation

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.

preprint2019arXiv

SDSS-IV MaNGA: Bar pattern speed estimates with the Tremaine-Weinberg method and their error sources

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.

preprint2016arXiv

Analysis of the Very Inner Milky Way Dark Matter Distribution and Gamma-Ray Signals

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.

preprint2016arXiv

Close encounters involving RAVE stars beyond the 47 Tucanae tidal radius

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.

preprint2016arXiv

Discovery of a Metal-Poor Field Giant with a Globular Cluster Second-Generation Abundance Pattern

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 $ω$ 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.

preprint2015arXiv

RAVE stars tidally stripped/ejected from $ω$ Centauri globular cluster

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 $ω$ Centauri. The orbits were integrated in a Milky-Way-like axisymmetric Galactic potential, ignoring the effects of the dynamical evolution of $ω$ 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 $ω$ Centauri: (\textit{i}) 8 stars with relative velocities $V_{rel}< 200 $ kms$^{-1}$ may have been ejected $\sim$ 200 Myr ago from $ω$ 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 $ω$ Centauri. We also confirm the link between J131340.4-484714 as potential member of $ω$ Centauri, and probably ejected $\sim$ 2.0 Myr ago, with a relative velocity $V_{rel}\sim80$ kms$^{-1}$.

preprint2014arXiv

Searching for tidal tails around $ω$ Centauri using RR Lyrae Stars

We present a survey for RR Lyrae stars in an area of 50 deg$^2$ around the globular cluster $ω$ 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 $ω$ 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 $ω$ Centauri, confirming previous studies in the region. It is puzzling that tidal debris have been found elsewhere but not near the cluster itself.

preprint2014arXiv

Simulations of isolated dwarf galaxies formed in dark matter halos with different mass assembly histories

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.

preprint2013arXiv

Clumpy Disc and Bulge Formation

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.

preprint2013arXiv

Non-Axisymmetric Structure in the Satellite Dwarf Galaxy NGC2976: Implications for its Dark/Bright Mass Distribution and Evolution

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 $μ$ 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.

preprint2011arXiv

The Bulgeless Seyfert/LINER Galaxy NGC 3367: Disk, Bar, Lopsidedness and Environment

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.

preprint2011arXiv

The specific star formation rate and stellar mass fraction of low-mass central galaxies in cosmological simulations

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

preprint2011arXiv

Understanding the spiral structure of the Milky Way using the local kinematic groups

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.

preprint2008arXiv

Dynamics of barred galaxies: effects of disk height

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.