Within Rotation Curves (RC) is encoded the kinematical state of the stellar disc as well as information about the dynamical mechanisms driving the secular evolution of galaxies. To explain the characteristic features of RCs that arise by the influence of spiral patterns and bar, we study the kinematics of the stellar disc in a set of spiral galaxy models specifically tailored for this purpose. We find that, for our models, the induced non-circular motions are more prominent for spirals with larger pitch angle, the ones typical in late type galaxies. Moreover, inside corotation, stars rotate slower along the spiral arms than along the inter-arm, which translates into a local minima or maxima in the RC, respectively. We also see, from off-plane RC, that the rotation is faster for stars that at observed closer to the plane, and diminishes as one looks farther off plane; this trend is more noticeable in our Sa galaxy model than our Sc galaxy model. Additionally, in a previous work we found that the diagonal ridges in the $V_{\phi}-R$ plane, revealed through the Gaia DR2, have a resonant origin due to the spiral arms and bar and that these ridges project themselves as wiggles in the RC; here, we further notice that the development of these ridges, and the development of high orbital eccentricities in the stellar disc are the same. Hence, we conclude that, the following explanations of bumps and wiggles in RCs are equivalent: they are manifestations of diagonal ridges in the $V_{\phi}-R$ plane, or of the rearrangement of the orbital eccentricities in the stellar disc.

We explore the long-term evolution of a bias-free orbital representation of the cometary nuclei (with diameters above 2 km) of the Kuiper belt, using the so-called L7 synthetic model from CFEPS, which consists of three dynamical sub-populations: the Classical, the Resonant, and the Scattering. The dynamical evolution of belt particles is studied under the gravitational influence of the Sun and the four giant planets, as well as of the 34 largest known trans-Neptunian objects (TNOs with $H_V < 4$). Here we indistinctly call Dwarf Planets (DPs) to the full sample of 34 large TNOs. Over a 1 Gyr timescale, we analyze the secular influence of the DPs over Kuiper belt disk particles and their contribution to the injection rate of new visible Jupiter Family Comets (JFCs). We found that DPs globally increase the number of JFCs by 12.6% when compared with the comets produced by the giant planets alone. When considering each population separately, we find that the increment produced by DPs is 17%, 12%, and 3% for the Classical, Resonant, and Scattering populations, respectively. Given the rate of escapes from the Kuiper belt, we find upper limits to the number of objects in each population required to maintain the JFCs in steady-state; the results are $55.9\times10^6$, $78.5\times10^6$, and $274.3\times10^6$ for the Scattering, Resonant, and Classical populations, respectively. Finally, we find that the Plutinos are the most important source of comets which were originally in a resonant configuration, where the presence of Pluto alone enhances by 10% the number of JFCs.

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.

As part of our multi-observatory, multi-filter campaign, we present \rmi color observations of 82 Near-Earth Objects (NEOs) obtained with the RATIR instrument on the 1.5m robotic telescope at the San Pedro Martir's National Observatory in Mexico. Our project is particularly focused on rapid response observations of small ($\lesssim 850$ m) NEOs. The rapid response and the use of spectrophotometry allows us to constrain the taxonomic classification of NEOs with high efficiency. Here we present the methodology of our observations and our result, suggesting that the ratio of C-type to S-type asteroids in a size range of $\sim$30-850m is 1.1, which is in accordance with our previous results. We also find that 10$\%$ of all NEOs in our sample are neither C- nor S-type asteroids

The PERLAS potential has been successfully used in many studies related with the dynamics of the spiral arms \textit{on} the equatorial plane of normal (non-barred) spiral galaxies. In the present work we extend these studies by investigating the three-dimensional dynamics of the spiral arms in the same type of potential. We consider a typical open, logarithmic, spiral pattern of pitch angle 25$^{\circ}$ and we examine the stellar orbits that can support it as the ratio of the masses of the spiral over the disk component ($M_{s}/M_{d}$) varies. We indicate the families of `three-dimensional' periodic orbits that act as the backbone of the spiral structure and we discuss their stability in the models we present. We study further the quasi-periodic and non-periodic orbits in general that follow spiral-supporting orbits as the $M_{s}/M_{d}$ ratio increases. We find that a bisymmetric spiral with 25$^{\circ}$ pitch angle is better supported by orbits in models with $0.03\lessapprox M_{s}/M_{d} \lessapprox 0.07$. In these cases a strong spiral pattern is supported between the radial 2:1 and 4:1 resonances, while local enhancements of the imposed spirals are encountered in some models between 4:1 and corotation. A characteristic bar-like structure is observed in all models at radii smaller than the radius of the 2:1 resonance.

In this work we explore a new dynamical path for the delivery of low-inclination comets. In a configuration formed by an interior giant planet and an exterior massive debris disk, where the mass is accounted for by the 50 largest objects in the disk, the strongest mean motion resonances of the giant, located along the belt, are replenished with new material (test particles) due to the influence of the 50 massive objects. Once in resonance, slow chaotic diffusion stirs the orbital elements of the cometary nuclei enough to encounter the giant and to be scattered by it. When the disk is massive enough, both resonant and non-resonant particles are stirred quickly to encounter the giant and form an scattered disk component, greatly increasing the rate for the delivery of cometary material to the inner part of the system. This mechanism is applicable both to the solar system and extrasolar systems in general. Preliminary results, using a disk as massive as the classical Kuiper belt, indicate that the mechanism here proposed can account for about a tenth of the required injection rate to maintain the population of ecliptic comets in steady state. In a more massive belt of 0.25 M$_\oplus$, an estimated rate of around 0.6 new comets per year is found. Such a high rate would pose a serious risk for the habitability of rocky interior planets, yet would resemble the late heavy bombardment that was present in the early solar system.

In this paper we report a new estimate of the absolute proper motion (PM) of the globular cluster NGC 5139 ($\omega$ Cen) as part of the HST large program GO-14118+14662. We analyzed a field 17 arcmin South-West of the center of $\omega$ Cen and computed PMs with an epoch span of $\sim$15.1 years. We employed 45 background galaxies to link our relative PMs to an absolute reference-frame system. The absolute PM of the cluster in our field is: $(\mu_\alpha \cos\delta , \mu_\delta) = (-3.341 \pm 0.028 , -6.557 \pm 0.043)$ mas yr$^{-1}$. Upon correction for the effects of viewing perspective and the known cluster rotation, this implies that for the cluster center of mass $(\mu_\alpha \cos\delta , \mu_\delta) = (-3.238 \pm 0.028, -6.716 \pm 0.043)$ mas yr$^{-1}$. This measurement is direct and independent, has the highest random and systematic accuracy to date, and will provide an external verification for the upcoming Gaia Data Release 2. It also differs from most reported PMs for $\omega$ Cen in the literature by more than 5$\sigma$, but consistency checks compared to other recent catalogs yield excellent agreement. We computed the corresponding Galactocentric velocity, calculated the implied orbit of $\omega$ Cen in two different Galactic potentials, and compared these orbits to the orbits implied by one of the PM measurements available in the literature. We find a larger (by about 500 pc) perigalactic distance for $\omega$ Cen with our new PM measurement, suggesting a larger survival expectancy for the cluster in the Galaxy.

NGC 6791 is one of the most studied open clusters, it is massive ($\sim5000\,M_{\odot}$), located at the solar circle, old ($~8\,$Gyr) and yet the most metal-rich cluster (${\rm [Fe/H]}\simeq0.4$) known in the Milky Way. By performing an orbital analysis within a Galactic model including spiral arms and a bar, we found that it is plausible that NGC 6791 formed in the inner thin disc or in the bulge, and later displaced by radial migration to its current orbit. We apply different tools to simulate NGC 6791, including direct $N$-body summation in time-varying potentials, to test its survivability when going through different Galactic environments. In order to survive the 8 Gyr journey moving on a migrating orbit, NGC 6791 must have been more massive, $M_0 \geq 5\times10^4 M_{\odot}$, when formed. We find independent confirmation of this initial mass in the stellar mass function, which is observed to be flat; this can only be explained if the average tidal field strength experienced by the cluster is stronger than what it is at its current orbit. Therefore, the birth place and journeys of NGC 6791 are imprinted in its chemical composition, in its mass loss, and in its flat stellar mass function, supporting its origin in the inner thin disc or in the bulge.

We provide new insight on the origin of the cold high-V$_{\rm los}$ peaks ($\sim$200 kms$^{-1}$) in the Milky Way bulge discovered in the APOGEE commissioning data \citep{Nidever2012}. Here we show that such kinematic behaviour present in the field regions towards the Galactic bulge is not likely associated with orbits that build the boxy/peanut (B/P) bulge. To this purpose, a new set of test particle simulations of a kinematically cold stellar disk evolved in a 3D steady-state barred Milky Way galactic potential, has been analysed in detail. Especially bar particles trapped into the bar are identified through the orbital Jacobi energy $E_{J}$, which allows us to identify the building blocks of the B/P feature and investigate their kinematic properties. Finally, we present preliminary results showing that the high-V$_{\rm los}$ features observed towards the Milky Way bulge are a natural consequence of a large-scale \textit{midplane} particle structure, which is unlikely associated with the Galactic bar.

Albeit radial migration must be a ubiquitous process in disc galaxies, its significance in the evolution of stellar discs is not always reflected through global trends. However, there are other key observables, such as the metallicity distribution function (MDF), that may shed some light in this matter. We argue that the shape of the MDF not only tells us whether the stellar disc experienced radial migration, but it also contains important clues on the structure that triggered it. Specifically, the MDF contains information about the dynamics and morphology of the spiral pattern. To constrain the spiral parameters, we have included a detailed chemical tagging in our simulations; this allows us to produce a restriction of the structural parameters of the spiral arms in the Milky Way as well as a method to constrain chemical evolution models towards the center of the Galactic disc, where no chemical model provides information. We also found that it is unlikely that the Sun was formed near its current galactocentric position, therefore it might be inaccurate to consider the Sun as representative of the chemical abundances in the solar neighborhood. We also show that a stellar disc of the Milky Way, after evolving dynamically and chemically for 5 Gyr, preserves 80% of its original global metallicity gradient despite having suffered important heating and radial migration; this means that the presence of a metallicity gradient in a given galaxy, does not guarantee that radial mixing has not played a role in its evolution.

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 have explored the evolution of a cold debris disk under the gravitational influence of dwarf planet sized objects (DPs), both in the presence and absence of an interior giant planet. Through detailed long-term numerical simulations, we demonstrate that, when the giant planet is not present, DPs can stir the eccentricities and inclinations of disk particles, in linear proportion to the total mass of the DPs; on the other hand, when the giant planet is included in the simulations, the stirring is approximately proportional to the mass squared. This creates two regimes: below a disk mass threshold (defined by the total mass of DPs), the giant planet acts as a stabilizing agent of the orbits of cometary nucleii, diminishing the effect of the scatterers; above the threshold, the giant contributes to the dispersion of the particles.

It is a common assumption that high-altitude open clusters live longer compared with clusters moving close to the Galactic plane. This is because at high altitudes, open clusters are far from the disruptive effects of in-plane substructures, such as spiral arms, molecular clouds and the bar. However, an important aspect to consider in this scenario is that orbits of high-altitude open clusters will eventually cross the Galactic plane, where the vertical tidal field of the disk is strong. In this work we simulate the interaction of open clusters with the tidal field of a detailed Milky Way Galactic model at different average altitudes and galactocentric radii. We find that the life expectancy of clusters decreases as the maximum orbital altitude increases and reaches a minimum at altitudes of approximately 600 pc. Clusters near the Galactic plane live longer because they do not experience strong vertical tidal shocks from the Galactic disk; then, for orbital altitudes higher than 600 pc, clusters start again to live longer due to the decrease in the number of encounters with the disk. With our study, we find that the compressive nature of the tides in the arms region and the bar have an important role on the survival of small clusters by protecting them from disruption: clusters inside the arms can live up to twice as long as those outside the arms at similar galactocentric distance.

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.

Recent observations show that the Milky Way's metallicity distribution function (MDF) changes its shape as a function of radius. This new evidence of radial migration within the stellar disc sets additional constraints on Galactic models. By performing controlled test particle simulations in a very detailed, observationally motivated model of the Milky Way, we demonstrate that, in the inner region of the disc, the MDF is shaped by the joint action of the bar and spiral arms, while at outer radii the MDF is mainly shaped by the spiral arms. We show that the spiral arms are able to imprint their signature in the radial migration, shaping the MDF in the outskirts of the Galactic disc with a minimal participation of the bar. Conversely, this work has the potential to characterise some structural and dynamical parameters of the spiral arms based on radial migration and the shape of the MDF. Finally, the resemblance obtained with this approximation to the MDF curves of the Galaxy as seen by APOGEE, show that a fundamental factor influencing their shape is the Galactic potential.

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.

Using the Wide Field Camera for the United Kingdom Infrared Telescope, we measure the near-infrared colors of near-Earth objects (NEOs) in order to put constraints on their taxonomic classifications. The rapid-response character of our observations allows us to observe NEOs when they are close to the Earth and bright. Here we present near-infrared color measurements of 86 NEOs, most of which were observed within a few days of their discovery, allowing us to characterize NEOs with diameters of only a few meters. Using machine-learning methods, we compare our measurements to existing asteroid spectral data and provide probabilistic taxonomic classifications for our targets. Our observations allow us to distinguish between S-complex, C/X-complex, D-type, and V-type asteroids. Our results suggest that the fraction of S-complex asteroids in the whole NEO population is lower than the fraction of ordinary chondrites in the meteorite fall statistics. Future data obtained with UKIRT will be used to investigate the significance of this discrepancy.

We present a dynamical study of the effect of the bar and spiral arms on the simulated orbits of open clusters in the Galaxy. Specifically, this work is devoted to the puzzling presence of high-altitude open clusters in the Galaxy. For this purpose we employ a very detailed observationally motivated potential model for the Milky Way and a careful set of initial conditions representing the newly born open clusters in the thin disk. We find that the spiral arms are able to raise an important percentage of open clusters (about one-sixth of the total employed in our simulations, depending on the structural parameters of the arms) above the Galactic plane to heights beyond 200 pc, producing a bulge-shaped structure toward the center of the Galaxy. Contrary to what was expected, the spiral arms produce a much greater vertical effect on the clusters than the bar, both in quantity and height; this is due to the sharper concentration of the mass on the spiral arms, when compared to the bar. When a bar and spiral arms are included, spiral arms are still capable of raising an important percentage of the simulated open clusters through chaotic diffusion (as tested from classification analysis of the resultant high-z orbits), but the bar seems to restrain them, diminishing the elevation above the plane by a factor of about two.

With the use of long-term numerical simulations, we study the evolution and orbital behavior of cometary nuclei in cold Kuiper belt-like debris disks under the gravitational influence of dwarf planets (DPs); we carry out these simulations with and without the presence of a Neptune-like giant planet. This exploratory study shows that in the absence of a giant planet, 10 DPs are enough to induce strong radial and vertical heating on the orbits of belt particles. On the other hand, the presence of a giant planet close to the debris disk, acts as a stability agent reducing the radial and vertical heating. With enough DPs, even in the presence of a Neptune-like giant planet some radial heating remains; this heating grows steadily, re-filling resonances otherwise empty of cometary nuclei. Specifically for the solar system, this secular process seems to be able to provide material that, through resonant chaotic diffusion, increase the rate of new comets spiraling into the inner planetary system, but only if more than the $\sim10$ known DP sized objects exist in the trans-Neptunian region.

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

Making use of a set of detailed potential models for normal spiral galaxies, we analyze the disk stellar orbital dynamics as the structural and dynamical parameters of the spiral arms (mass, pattern speed and pitch angle) are gradually modified. With this comprehensive study of ordered and chaotic behavior, we constructed an assemblage of orbitally supported galactic models and plausible parameters for orbitally self-consistent spiral arms models. We find that, to maintain orbital support for the spiral arms, the spiral arm mass, M$_{sp}$, must decrease with the increase of the pitch angle, $i$; if $i$ is smaller than $\sim10\deg$, M$_{sp}$ can be as large as $\sim7\%$, $\sim6\%$, $\sim5\%$ of the disk mass, for Sa, Sb, and Sc galaxies, respectively. If $i$ increases up to $\sim25\deg$, the maximum M$_{sp}$ is $\sim1\%$ of the disk mass independently in this case of morphological type. For values larger than these limits, spiral arms would likely act as transient features. Regarding the limits posed by extreme chaotic behavior, we find a strong restriction on the maximum plausible values of spiral arms parameters on disk galaxies beyond which, chaotic behavior becomes pervasive. We find that for $i$ smaller than $\sim20\deg$, $\sim25\deg$, $\sim30\deg$, for Sa, Sb, and Sc galaxies, respectively, M$_{sp}$ can go up to $\sim10\%$, of the mass of the disk. If the corresponding $i$ is around $\sim40\deg$, $\sim45\deg$, $\sim50\deg$, M$_{sp}$ is $\sim1\%$, $\sim2\%$, $\sim3\%$ of the mass of the disk. Beyond these values, chaos dominates phase space, destroying the main periodic and the neighboring quasi-periodic orbits.

With the use of a background Milky-Way-like potential model, we performed stellar orbital and magnetohydrodynamic (MHD) simulations. As a first experiment, we studied the gaseous response to a bisymmetric spiral arm potential: the widely employed cosine potential model and a self-gravitating tridimensional density distribution based model called PERLAS. Important differences are noticeable in these simulations, while the simplified cosine potential produces two spiral arms for all cases, the more realistic density based model produces a response of four spiral arms on the gaseous disk, except for weak arms -i.e. close to the linear regime- where a two-armed structure is formed. In order to compare the stellar and gas response to the spiral arms, we have also included a detailed periodic orbit study and explored different structural parameters within observational uncertainties. The four armed response has been explained as the result of ultra harmonic resonances, or as shocks with the massive bisymmetric spiral structure, among other. From the results of this work, and comparing the stellar and gaseous responses, we tracked down an alternative explanation to the formation of branches, based only on the orbital response to a self-gravitating spiral arms model. The presence of features such as branches, might be an indication of transiency of the arms.

With the use of a detailed Milky Way nonaxisymmetric potential, observationally and dynamically constrained, the effects of the bar and the spiral arms in the Galaxy are studied in the disc and in the stellar halo. Especially the trapping of stars in the disc and Galactic halo by resonances on the Galactic plane, induced by the Galactic bar, has been analysed in detail. To this purpose, a new method is presented to delineate the trapping regions using empirical diagrams of some orbital properties obtained in the Galactic potential. In these diagrams we plot in the inertial Galactic frame a characteristic orbital energy versus a characteristic orbital angular momentum, or versus the orbital Jacobi constant in the reference frame of the bar, when this is the only nonaxisymmetric component in the Galactic potential. With these diagrams some trapping regions are obtained in the disc and halo using a sample of disc stars and halo stars in the solar neighbourhood. We compute several families of periodic orbits on the Galactic plane, some associated with this resonant trapping. In particular, we find that the trapping effect of these resonances on the Galactic plane can extend several kpc from this plane, trapping stars in the Galactic halo. The purpose of our analysis is to investigate if the trapping regions contain some known moving groups in our Galaxy. We have applied our method to the Kapteyn group, a moving group in the halo, and we have found that this group appears not to be associated with a particular resonance on the Galactic plane.

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.

The first mechanism invoked to explain the existence of the thick disk in the Milky Way Galaxy, were the spiral arms. Up-to-date work summon several other possibilities that together seem to better explain this component of our Galaxy. All these processes must affect differently in distinct types of galaxies, but the contribution of each one has not been straightforward to quantify. In this work, we present a first comprehensive study of the effect of the spiral arms in the formation of thick disks, as going from early to late type disk galaxies, in an attempt to characterize and quantify this specific mechanism in galactic potentials. To this purpose, we perform numerical simulations of test particles in a three-dimensional spiral galaxy potential of normal spiral galaxies (from early to late types). By varying the parameters of the spiral arms we found that the vertical heating of the stellar disk becomes very important in some cases, and strongly depends on the galaxy morphology, pitch angle, arms mass and its pattern speed. The later the galaxy type, the larger is the effect on the disk heating. This study shows that the physical mechanism causing the vertical heating is different from simple resonant excitation. The spiral pattern induce chaotic behavior not linked necessarily to resonances but to direct scattering of disk stars, which leads to an increase of the velocity dispersion. We applied this study to the specific example of the Milky Way Galaxy, for which we have also added an experiment that includes the Galactic bar. From this study we deduce that the effect of spiral arms of a Milky-Way-like potential, on the dynamical vertical heating of the disk is negligible, unlike later galactic potentials for disks.

The orbital elements of comet Halley are known to a very high precision, suggesting that the calculation of its future dynamical evolution is straightforward. In this paper we seek to characterize the chaotic nature of the present day orbit of comet Halley and to quantify the timescale over which its motion can be predicted confidently. In addition, we attempt to determine the timescale over which its present day orbit will remain stable. Numerical simulations of the dynamics of test particles in orbits similar to that of comet Halley are carried out with the Mercury 6.2 code. On the basis of these we construct survival time maps to assess the absolute stability of Halley's orbit, frequency analysis maps, to study the variability of the orbit and we calculate the Lyapunov exponent for the orbit for variations in initial conditions at the level of the present day uncertainties in our knowledge of its orbital parameters. On the basis of our calculations of the Lyapunov exponent for comet Halley, the chaotic nature of its motion is demonstrated. The e-folding timescale for the divergence of initially very similar orbits is approximately 70 years. The sensitivity of the dynamics on initial conditions is also evident in the self-similarity character of the survival time and frequency analysis maps in the vicinity of Halley's orbit, which indicates that, on average, it is unstable on a timescale of hundreds of thousands of years. The chaotic nature of Halley's present day orbit implies that a precise determination of its motion, at the level of the present day observational uncertainty, is difficult to predict on a timescale of approximately 100 years. Furthermore, we also find that the ejection of Halley from the solar system or its collision with another body could occur on a timescale as short as 10,000 years.

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.

We calculate orbits, tidal radii, and bulge-bar and disk shocking destruction rates for 63 globular clusters in our Galaxy. Orbits are integrated in both an axisymmetric and a non-axisymmetric Galactic potential that includes a bar and a 3D model for the spiral arms. With the use of a Monte Carlo scheme, we consider in our simulations observational uncertainties in the kinematical data of the clusters. In the analysis of destruction rates due to the bulge-bar, we consider the rigorous treatment of using the real Galactic cluster orbit, instead of the usual linear trajectory employed in previous studies. We compare results in both treatments. We find that the theoretical tidal radius computed in the nonaxisymmetric Galactic potential compares better with the observed tidal radius than that obtained in the axisymmetric potential. In both Galactic potentials, bulge-shocking destruction rates computed with a linear trajectory of a cluster at its perigalacticons give a good approximation to the result obtained with the real trajectory of the cluster. Bulge-shocking destruction rates for clusters with perigalacticons in the inner Galactic region are smaller in the non-axisymmetric potential, as compared with those in the axisymmetric potential. For the majority of clusters with high orbital eccentricities (e > 0.5), their total bulge+disk destruction rates are smaller in the non-axisymmetric potential.

A general formulation to compute habitable zones (HZ) around binary stars is presented. A HZ in this context must satisfy two separate conditions: a radiative one and one of dynamical stability. For the case of single stars, the usual concept of circumstellar habitable zone is based on the radiative condition only, as the dynamical stability condition is taken for granted. For the radiative condition, we extend the simple formulation of the circumstellar habitable zone for single stars, to the case of eccentric stellar binary systems, where two sources of luminosity at different orbital phases contribute to the irradiance of their planetary circumstellar and circumbinary regions. Our approach considers binaries with eccentric orbits and guarantees that orbits in the computed habitable zone remain within it at all orbital phases. For the dynamical stability condition, we use the approach of invariant loops developed by Pichardo et al. 2005 to find regions of stable, non-intersecting orbits, which is a robust method to find stable regions in binary stars, as it is based in the existence of integrals of motion. We apply the combined criteria to calculate HZ for 64 binary stars in the solar neighborhood with known orbital parameters, including some with discovered planets. Formulae and interpolating tables are provided, so the reader can compute the boundaries of the HZ for an arbitrary binary system, using the stellar flux limits they prefer. Together with the formulae provided for stable zones, these allow the computation of both regions of stability and habitability around any binary stellar system. We found 56% of the cases we consider can satisfy both restrictions, this is a very important constriction to binary systems. Nevertheless, we conclude that these systems where a dynamical and radiative safe zone exists, must be considered strong candidates in the search for habitable planets

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.

Every Galactic environment is characterized by a stellar density and a velocity dispersion. With this information from literature, we simulated flyby encounters for several Galactic regions, numerically calculating stellar trajectories as well as orbits for particles in disks; our aim was to understand the effect of typical stellar flybys on planetary (debris) disks in the Milky Way Galaxy. For the Solar neighborhood, we examined nearby stars with known distance, proper motions, and radial velocities. We found occurrence of a disturbing impact to the Solar planetary disk within the next 8 Myr to be highly unlikely; perturbations to the Oort cloud seem unlikely as well. Current knowledge of the full phase space of stars in the Solar neighborhood, however, is rather poor, and thus we cannot rule out the existence of a star that is more likely to approach than those for which we have complete kinematic information. We studied the effect of stellar encounters on planetary orbits within the habitable zones of stars in more crowded stellar environments, such as stellar clusters. We found that in open clusters habitable zones are not readily disrupted; this is true if they evaporate in less than 108 years. For older clusters the results may not be the same. We specifically studied the case of Messier 67, one of the oldest open clusters known, and show the effect of this environment on debris disks. We also considered the conditions in globular clusters, the Galactic nucleus, and the Galactic bulge-bar. We calculated the probability of whether Oort clouds exist in these Galactic environments.

We built a family of non-axisymmetric potential models for normal non-barred or weakly-barred spiral galaxies as defined in the simplest classification of galaxies: the Hubble sequence. For this purpose a three-dimensional self-gravitating model for spiral arms PERLAS is superimposed to the galactic axisymmetric potentials. We analyze the stellar dynamics varying only the pitch angle of the spiral arms, from 4$\deg$ to 40$\deg$, for an Sa galaxy, from 8$\deg$ to 45$\deg$, for an Sb galaxy, and from 10$\deg$ to 60$\deg$, for an Sc galaxy. Self-consistency is indirectly tested through periodic orbital analysis, and through density response studies for each morphological type. Based on ordered behavior, periodic orbits studies show that for pitch angles up to approximately $15\deg$, $18\deg$, and $20\deg$ for Sa, Sb and Sc galaxies, respectively, the density response supports the spiral arms potential, a requisite for the existence of a long-lasting large-scale spiral structure. Beyond those limits, the density response tends to "avoid" the potential imposed by mantaining lower pitch angles in the density response; in that case the spiral arms may be explained as transient features rather than long-lasting large-scale structures. In a second limit, from a phase space orbital study based on chaotic behavior, we found that for pitch angles larger than $\sim30\deg$, $\sim40\deg$ and $\sim50\deg$ for Sa, Sb, and Sc galaxies, respectively, chaotic orbits dominate all phase space prograde region that surrounds the periodic orbits sculpting the spiral arms and even destroying them. This result seems to be in good agreement with observations of pitch angles in typical isolated normal spiral galaxies.

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.

It is generally assumed that gas in a galactic disk follows closely non self-intersecting periodic stellar orbits. In order to test this common assumption, we have performed MHD simulations of a galactic-like disk under the influence of a spiral galactic potential. We also have calculated the actual orbit of a gas parcel and compared it to stable periodic stellar orbits in the same galactic potential and position. We found that the gaseous orbits approach periodic stellar orbits far from the major orbital resonances only. Gas orbits initialized at a given galactocentric distance but at different azimuths can be different, and scattering is conspicuous at certain galactocentric radii. Also, in contrast to the stellar behaviour, near the 4:1 (or higher order) resonance the gas follows nearly circular orbits, with much shorter radial excursions than the stars. Also, since the gas does not settle into a steady state, the gaseous orbits do not necessarily close on themselves.

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.

Using the results of Pichardo \et (2005,2008), we determine regions of dynamical stability where planets (or discs in general) could survive in stable orbits around binary stellar systems. We produce this study for 161 binary stars in the Solar neighborhood with known orbital parameters. Additionally, we constructed numerically the discs (invariant loops) around five binary systems with known orbital parameters and with confirmed planets: HIP 10138, HIP 4954, HIP 67275, HIP 116727 and Kepler 16, as a test to the approximation of Pichardo et al. (2005,2008). In each single case, the reported position of the planets lay within our calculated stability regions. This study intends to provide a guide in the search for planets around binary systems with well know orbital parameters, since our method defines precise limits for the stable regions, where discs may have established and planets formed.

Using the most recent proper-motion determination of the old, Solar-metallicity, Galactic open cluster M 67, in orbital computations in a non-axisymmetric model of the Milky Way, including a bar and 3D spiral arms, we explore the possibility that the Sun once belonged to this cluster. We have performed Monte Carlo numerical simulations to generate the present-day orbital conditions of the Sun and M 67, and all the parameters in the Galactic model. We compute 3.5 \times 10^5 pairs of orbits Sun-M 67 looking for close encounters in the past with a minimum distance approach within the tidal radius of M 67. In these encounters we find that the relative velocity between the Sun and M 67 is larger than 20 km/s. If the Sun had been ejected from M 67 with this high velocity by means of a three-body encounter, this interaction would destroy an initial circumstellar disk around the Sun, or disperse its already formed planets. We also find a very low probability, much less than 10^-7, that the Sun was ejected from M 67 by an encounter of this cluster with a giant molecular cloud. This study also excludes the possibility that the Sun and M 67 were born in the same molecular cloud. Our dynamical results convincingly demonstrate that M67 could not have been the birth cluster of our Solar System.

We built models for low bulge mass spiral galaxies (late type as defined by the Hubble classification) using a 3-D self-gravitating model for spiral arms, and analyzed the orbital dynamics as a function of pitch angle, going from 10$\deg$ to 60$\deg$. Testing undirectly orbital self-consistency, we search for the main periodic orbits and studied the density response. For pitch angles up to approximately $\sim 20\deg$, the response supports closely the potential permitting readily the presence of long lasting spiral structures. The density response tends to "avoid" larger pitch angles in the potential, by keeping smaller pitch angles in the corresponding response. Spiral arms with pitch angles larger than $\sim 20\deg$, would not be long-lasting structures but rather transient. On the other hand, from an extensive orbital study in phase space, we also find that for late type galaxies with pitch angles larger than $\sim 50\deg$, chaos becomes pervasive destroying the ordered phase space surrounding the main stable periodic and quasi-periodic orbits and even destroying them. This result is in good agreement with observations of late type galaxies, where the maximum observed pitch angle is $\sim 50\deg$.

In Papers I and II of this series, the existence of two distinct halo populations of stars have been found in the solar neighborhood. Precise relative ages and orbital parameters are determined for 67 halo and 16 thick-disk stars having metallicities in the range -1.4 < [Fe/H] < -0.4 to better understand the context of the two halo populations in the formation and evolution of the Galaxy. Ages are derived by comparing the positions of stars in the logT_{eff}-log(g) diagram with isochrones from the Y^2 models interpolated to the exact [Fe/H] and [alpha/Fe] values of each star. Possible systematic errors in T_{eff} and log(g) are considered and corrected. With space velocities from Paper I as initial conditions, orbital integrations have been carried out using a detailed, observationally constrained Milky Way model including a bar and spiral arms. The `high-alpha' halo stars have ages 2-3 Gyr larger than the `low-alpha' ones. The orbital parameters show very distinct differences between the `high-alpha' and `low-alpha' halo stars. The `low-alpha' ones have r_{max}'s to 30-40 kpc, z_{max}'s to approx. 18 kpc, and e_{max}'s clumped at values greater than 0.85, while the `high-alpha' ones, r_{max}'s to about 16 kpc, z_{max}'s to 6-8 kpc, and e_{max} more or less uniformly distributed over 0.4-1.0. A dual in situ-plus-accretion formation scenario best explains the existence and characteristics of these two halo populations, but one remaining defect is that this model is not consistent regarding the r_{max}'s obtained for the in situ `high-alpha' component; the predicted values are too small. It appears that omega Cen may have contributed in a significant way to the existence of the `low-alpha' component; recent models, including dynamical friction and tidal stripping, have produced orbital parameters as great as those of the `low-alpha' component.

We have computed trajectories, distances and times of closest approaches to the Sun by stars in the Solar neighbourhood with known position, radial velocity and proper motions. For this purpose we have used a full potential model of the Galaxy that reproduces the local z-force, the Oort constants, the local escape velocity, and the rotation curve of the Galaxy. From our sample we constructed initial conditions, within observational uncertainties, with a Monte Carlo scheme for the twelve most suspicious candidates because of their small tangential motion. We find that the star Gliese 710 will have the closest approach to the Sun, with a distance of approximately 0.34 pc at 1.36 Myr in the future. We show that the effect of a flyby with the characteristics of Gliese 710 on a 100 AU test particle disk representing the Solar system is negligible. However, since there is a lack of 6D data for a large percentage of stars in the Solar neighbourhood, closer approaches may exist. We calculate parameters of passing stars that would cause noticeable effects on the Solar disk. Regarding the birth cloud of the Sun, we performed experiments to reproduce roughly the observed orbital parameters such as eccentricities and inclinations of the Kuiper Belt. It is known now that in Galactic environments, such as stellar formation regions, the stellar densities of new born stars, are high enough to produce close encounters within 200 AU. Moreover, in these Galactic environments, the velocity dispersion is relatively low, typically approximately 1-3 km s-1.We find that with a velocity dispersion of approximately 1 km s-1 and an approach distance of about 150 AU, typical of these regions, we obtain approximately the eccentricities and inclinations seen in the current Solar system. Simple analytical calculations of stellar encounters effects on the Oort cloud are presented.

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.

We derived the absolute proper motion (PM) of the old, solar-metallicity Galactic open cluster M67 using observations collected with CFHT (1997) and with LBT (2007). About 50 galaxies with relatively sharp nuclei allow us to determine the absolute PM of the cluster. We find (mu_alpha cos(delta),mu_delta)_J2000.0 = (-9.6+/-1.1,-3.7+/-0.8) mas/yr. By adopting a line-of-sight velocity of 33.8+/-0.2 km/s, and assuming a distance of 815+/-50 pc, we explore the influence of the Galactic potential, with and without the bar and/or spiral arms, on the galactic orbit of the cluster.

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.

Only a few molecular clouds in the Solar Neighborhood exhibit the formation of only low-mass stars. Traditionally, these clouds have been assumed to be supported against more vigorous collapse by magnetic fields. The existence of strong magnetic fields in molecular clouds, however, poses serious problems for the formation of stars and of the clouds themselves. In this {\em Letter}, we review the three-dimensional structure and kinematics of Taurus --the archetype of a region forming only low-mass stars-- as well as its orientation within the Milky way. We conclude that the particularly low star-formation efficiency in Taurus may naturally be explained by tidal forces from the Galaxy, with no need for magnetic regulation or stellar feedback.

(Abridged) The gravitational term for clouds and cores entering in the virial theorem is usually assumed to be equal to the gravitational energy, since the contribution to the gravitational force from the mass distribution outside the volume of integration is assumed to be negligible. Such approximation may not be valid in the presence of an important external net potential. In the present work we analyze the effect of an external gravitational field on the gravitational budget of a density structure. Our cases under analysis are (a) a giant molecular cloud (GMC) with different aspect ratios embedded within a galactic net potential, and (b) a molecular cloud core embedded within the gravitational potential of its parent molecular cloud. We find that for roundish GMCs, the tidal tearing due to the shear in the plane of the galaxy is compensated by the tidal compression in the z direction. The influence of the external effective potential on the total gravitational budget of these clouds is relatively small, although not necessarily negligible. However, for more filamentary GMCs, the external effective potential can be dominant and can even overwhelm self-gravity, regardless of whether its main effect on the cloud is to disrupt it or compress it. This may explain the presence of some GMCs with few or no signs of massive star formation, such as the Taurus or the Maddalena's clouds. In the case of dense cores embedded in their parent molecular cloud, we found that the gravitational content due to the external field may be more important than the gravitational energy of the cores themselves. This effect works in the same direction as the gravitational energy, i.e., favoring the collapse of cores. We speculate on the implications of these results for star formation models.

In a previous work (Pichardo et al. 2005), we studied stable configurations for circumstellar discs in eccentric binary systems. We searched for "invariant loops": closed curves (analogous to stable periodic orbits in time-independent potentials) that change shape with the binary orbital phase, as test particles in them move under the influence of the binary potential. This approach allows us to identify stable configurations when pressure forces are unimportant, and dissipation acts only to prevent gas clouds from colliding with one another. We now extend this work to study the main geometrical properties of circumbinary discs. We have studied more than 100 cases with a range in eccentricity 0 .le. e .le. 0.9, and mass ratio 0.1 .le. q .le. 0.9. Although gas dynamics may impose further restrictions, our study sets lower stable bounds for the size of the central hole in a simple and computationally cheap way, with a relation that depends on the eccentricity and mass ratio of the central binary. We extend our previous studies and focus on an important component of these systems: circumbinary discs. The radii for stable orbits that can host gas in circumbinary discs are sharply constrained as a function of the binary's eccentricity. The circumbinary disc configurations are almost circular, with eccentricity e_d < 0.15, but if the mass ratio is unequal the disk is offset from the center of mass of the system. We compare our results with other models, and with observations of specific systems like GG Tauri A, UY Aurigae, HD 98800 B, and Fomalhaut, restricting the plausible parameters for the binary.

Absolute proper motions for six new globular clusters have recently been determined. This motivated us to obtain the Galactic orbits of these six clusters both in an axisymmetric Galactic potential and in a barred potential, such as the one of our Galaxy. Orbits are also obtained for a Galactic potential that includes spiral arms. The orbital characteristics are compared and discussed for these three cases. Tidal radii and destruction rates are also computed and discussed.

An analytical approximation to periodic orbits in the circular restricted three-body problem is provided. The formulation given in this work is based in calculations known from classical mechanics, but with the addition of the necessary terms to give a fairly good approximation that we compare with simulations, resulting in a simple set of analytical expressions that solve periodic orbits on discs of binary systems without the need of solving the motion equations by numerical integrations.

The effect of a barred potential (such as the one of the Milky Way) on the galactic orbits of forty-eight globular clusters for which absolute proper motions are known is studied. The orbital characteristics are compared with those obtained for the case of an axisymmetric galactic potential. Tidal radii are computed and discussed for both the better known axisymmetric case and that including a bar. The destruction rates due to bulge and disk shocking are calculated and compared in both galactic potentials.

We explore test particle orbits in the orbital plane of eccentric stellar binary systems, searching for ``invariant loops'': closed curves that change shape periodically as a function of binary orbital phase as the test particles in them move under the stars' gravity. Stable invariant loops play the same role in this periodically-varying potential as stable periodic orbits do in stationary potentials; in particular, when dissipation is weak, gas will most likely follow the non-intersecting loops, while nearby particle orbits librate around them. We use this method to set bounds on the sizes of disks around the stars, and on the gap between those and the inner edge of a possible circumbinary disk. Gas dynamics may impose further restrictions, but our study sets upper bounds for the size of circumstellar disks, and a lower bound for the inner radius of a circumbinary disk. We find that circumstellar disks are sharply reduced as the binary's eccentricity grows. Disks change in size and shape only marginally with the binary phase, with no strong preference to increase or decrease at any particular phase. The circumstellar disks in particular can be quite asymmetric. We compare our results with other numerical and theoretical results and with observations of the $\alpha$ Centauri and L1551 systems, finding very good agreement. The calculated changes in the shapes and crowding of the circumstellar orbits can be used to predict how the disk luminosity and mass inflow should vary with binary phase.

The gas response to a proposed spiral stellar pattern for our Galaxy is presented here as calculated via 2D hydrodynamic calculations utilizing the ZEUS code in the disk plane. The locus is that found by Drimmel (2000) from emission profiles in the K band and at 240 $\mu m$. The self-consistency of the stellar spiral pattern was studied in previous work (see Martos et al. 2004). It is a sensitive function of the pattern rotation speed, $\Omega_p$, among other parameters which include the mass in the spiral and its pitch angle. Here we further discuss the complex gaseous response found there for plausible values of $\Omega_p$ in our Galaxy, and argue that its value must be close to $20 km s^{-1} kpc^{-1}$ from the strong self-consistency criterion and other recent, independent studies which depend on such parameter. However, other values of $\Omega_p$ that have been used in the literature are explored to study the gas response to the stellar (K band) 2-armed pattern. For our best fit values, the gaseous response to the 2-armed pattern displayed in the K band is a four-armed pattern with complex features in the interarm regions. This response resembles the optical arms observed in the Milky Way and other galaxies with the smooth underlying two-armed pattern of the old stellar disk populations in our interpretation.

We report calculations of the stellar and gaseous response to a Milky Way mass distribution model including a spiral pattern with a locus as traced by K-band observations, over imposed on the axisymmetric components in the plane of the disk. The stellar study extends calculations from previous work concerning the self-consistency of the pattern. The stellar response to the imposed spiral mass is studied via computations of the central family of periodic and nearby orbits as a function of the pattern rotation speed, $\Omega_p$, among other parameters. A fine grid of values of $\Omega_p$ was explored ranging from 12 to 25 $km s^{-1} kpc^{-1}$. Dynamical self-consistency is highly sensitive to $\Omega_p$, with the best fit appearing at 20 $km s^{-1} kpc^{-1}$. We give an account of recent independent pieces of theoretical and observational work that are dependent on the value of $\Omega_p$, all of which are consistent with the value found here; the recent star formation history of the Milky Way, local inferences of cosmic ray flux variations and Galactic abundance patterns. The gaseous response, which is also a function of $\Omega_p$, was calculated via 2D hydrodynamic simulations with the ZEUS code. For $\Omega_p = 20 km s^{-1} kpc^{-1}$, the response to a two-armed pattern is a structured pattern of 4 arms, with bifurcations along the arms and interarm features. The pattern resembles qualitatively the optical arms observed in our Galaxy and other galaxies. The complex gaseous pattern appears to be linked to resonances in stellar orbits. Among these, the 4:1 resonance plays an important role, as it also determines the extent over which the imposed spiral pattern is dynamically self-consistent.

We built three models for the gravitational field of the Galactic bar. These models are an inhomogeneous ellipsoid, an inhomogeneous prolate spheroid, and a superposition of four inhomogeneous ellipsoids. Among the three models, the superposition provides our best approximation to the observed boxy mass distribution of the Galactic bar. Adding the bar component to an axisymmetric Galactic model, we have calculated stellar midplane orbits and orbits of some globular clusters with known kinematical data. For all models we find a secular dispersion effect upon the orbital energy and angular momentum, as measured in the Galactic inertial frame. This effect might be relevant to explain the orbital prograde-retrograde distribution of globular clusters. For the stellar kinematics, we study the connection between the sense of orbital motion in the midplane and the onset of chaos in the presence of the bar. In the inner region of the bar, chaos is induced by an axisymmetric central component (bulge) and it arises in orbits that change its orbital sense from prograde to retrograde and vice versa as seen from an inertial reference frame. Outside the bar region, chaos appears only in prograde orbits. Our results concerning such connection are consistent and extend those obtained for midplane orbits in the presence of only a spiral pattern in the axisymmetric Galactic model.

With the aim of studying the nonlinear stellar and gaseous response to the gravitational potential of a galaxy such as the Milky Way, we have modeled 3D galactic spiral arms as a superposition of inhomogeneous oblate spheroids and added their contribution to an axisymmetric model of the Galactic mass distribution. Three spiral loci are proposed here, based in different sets of observations. A comparison of our model with a tight-winding approximation shows that the self-gravitation of the whole spiral pattern is important in the middle and outer galactic regions. As a first step to full 3D calculations the model is suitable for, we have explored the stellar orbital structure in the midplane of the Galaxy. We present the standard analysis in the pattern rotating frame, and complement this analysis with orbital information from the Galactic inertial frame. Prograde and retrograde orbits are defined unambiguously in the inertial frame, then labeled as such in the Poincar\'e diagrams of the non-inertial frame. In this manner we found a sharp separatrix between the two classes of orbits. Chaos is restricted to the prograde orbits, and its onset occurs for the higher spiral perturbation considered plausible in our Galaxy.

In this paper we test the results of a recent analytical study by Lazarian and Pogosyan, on the statistics of emissivity in velocity channel maps, in the case of realistic density and velocity fields obtained from numerical simulations of MHD turbulence in the interstellar medium (ISM). To compensate for the lack of well-developed inertial ranges in the simulations due to the limited resolution, we apply a procedure for modifying the spectral slopes of the fields while preserving the spatial structures. We find that the density and velocity are moderately correlated in space and prove that the analytical results by Lazarian and Pogosyan hold in the case when these fields obey the fluid conservation equations. Our results imply that the spectra of velocity and density can be safely recovered from the position-position-velocity (PPV) data cubes available through observations, and confirm that the relative contributions of the velocity and density fluctuations to those of the emissivity depend on the velocity resolution used and on the steepness of the density spectral index. Furthermore, this paper supports previous reports that an interpretation of the features in the PPV data cubes as simple density enhancements (i.e., ``clouds'') can be often erroneous, as we observe that changes in the velocity statistics substantially modify the statistics of emissivity within the velocity data cubes.

Among the various phenomena observed in interacting galaxies is the ejection due to tidal forces of material into the intergalactic medium and its subsequent rearranging which can lead to the formation of self-gravitating tidal dwarf galaxies (TDGs). We investigate this process with a detailed multiwavelength study of the interacting system Arp 245 (NGC 2992/3) and a numerical model of the collision. The system is observed at an early stage of the interaction, about 100 Myr after perigalacticon, though at a time when tidal tails have already developed. Our HI map shows prominent counterparts to the optical tails. Whereas the stellar and gaseous components of the plume that originates from NGC 2992 match, the stellar and HI tails emanating from NGC 2993 have a different morphology. In particular, the HI forms a ring, a feature that has been successfully reproduced by our numerical simulations. The HI emission in the system as a whole peaks at the tip of the NGC 2992 tail where a gas reservoir of about 10^9 Msun, about 60% of the HI towards NGC 2992, coincides with a star-forming optical condensation, A245N. The latter tidal object exhibits properties ranging between those of dwarf irregular galaxies (structural parameters, gas content, star formation rate) and those of spiral disks (metallicity, star formation efficiency, stellar population). Prompted by the questions raised for this particular object, we discuss some issues related to the definition and identification of TDGs (abridged)

We study the effects of projection of three-dimensional (3D) data onto the plane of the sky by means of numerical simulations of turbulence in the interstellar medium including the magnetic field, parameterized cooling and diffuse and stellar heating, self-gravity and rotation. We compare the physical-space density and velocity distributions with their representation in position-position-velocity (PPV) space (``channel maps''), noting that the latter can be interpreted in two ways: either as maps of the column density's spatial distribution (at a given line-of-sight (LOS) velocity), or as maps of the spatial distribution of a given value of the LOS velocity (weighted by density). This ambivalence appears related to the fact that the spatial and PPV representations of the data give significantly different views. First, the morphology in the channel maps more closely resembles that of the spatial distribution of the LOS velocity component than that of the density field, as measured by pixel-to-pixel correlations between images. Second, the channel maps contain more small-scale structure than 3D slices of the density and velocity fields, a fact evident both in subjective appearance and in the power spectra of the images. This effect may be due to a pseudo-random sampling (along the LOS) of the gas contributing to the structure in a channel map: the positions sampled along the LOS (chosen by their LOS velocity) may vary significantly from one position in the channel map to the next.

Whereas it is well known that gas may be driven into the core of merging disks, fueling a central AGN or a nuclear starburst, recent studies have shown that a significant fraction of the stellar/gaseous components is expelled into the intergalactic medium along tidal tails. The tidal debris might be dispersed in the intergalactic/intracluster medium where it adds to the diffuse background light such as that observed in the Coma cluster, might fall back towards the merger, or regroup to form a new generation of galaxies, the so-called tidal dwarf galaxies (TDGs). We illustrate in this paper these various phenomena with a detailed multi-wavelength study of the interacting system NGC 2992/3 (Arp 245).