We present the very first wide-field, 11 by 11 arcmin, optical spectral mapping of M16, one of the most famous star-forming regions in the Galaxy. The data were acquired with the new imaging Fourier transform spectrograph SITELLE mounted on the Canada-France-Hawaii Telescope (CFHT). We obtained three spectral cubes with R=10'000 (SN1 filter), 1500 (SN2 filter) and 600 (SN3 filter), centered on the Pillars of Creation and the HH216 flow, covering the main optical nebular emission lines: [O II]3726,29 (SN1), Hb, [O III]4959,5007 (SN2), [N II]6548,84, Ha, and [S II]6717,31 (SN3). We validate the performance, calibration, and data reduction of SITELLE, and analyze the structures in the large field-of-view in terms of their kinematics and nebular emission. We compared the SITELLE data to MUSE integral field observations and other spectroscopic and narrow-band imaging data to validate the performance of SITELLE. We computed gas-phase metallicities via the strong-line method, performed a pixel-by-pixel fit to the main emission lines to derive kinematics of the ionized gas, computed the mass-loss rate of the Eastern pillar (the Spire), and combined the SITELLE data with near-infrared narrow-band imaging to characterize the HH216 flow. The comparison with previously published fluxes demonstrates very good agreement. We disentangle the dependence of the gas-phase metallicities (derived via abundance-tracing line ratios) on the degree of ionization and obtain metallicities that are in excellent agreement with the literature. We confirm the bipolar structure of HH216, find evidence for episodic accretion from the source of the flow, and identify its likely driving source. We compute the mass-loss rate of the Spire pillar on the East side of the H II region and find excellent agreement with the correlation between the mass-loss rate and the ionizing photon flux from the nearby cluster NGC6611.

There are few warp kinematic models of the Galaxy able to characterise structure and kinematics. These models are necessary to study the lopsidedness of the warp and the twisting of the line-of-nodes of the stellar warp, already seen in gas and dust. We use the \Gaia~Data Release 2 astrometric data up to $G=20$mag to characterise the structure of the Galactic warp, the vertical motions and the dependency on the age. We use two populations up to galactocentric distances of $16$kpc, a young (OB-type) and an old (Red Giant Branch, RGB). We use the nGC3 PCM and LonKin methods based on the Gaia observables, together with 2D projections of the positions and proper motions in the Galactic plane. We confirm the age dependency of the Galactic warp, both in positions and kinematics, being the height of the Galactic warp of about $0.2$kpc for the OB sample and of $1.$kpc for the RGB at a galactocentric distance of $14$kpc. Both methods find that the onset radius is $12\sim 13$kpc for the OB sample and $10\sim 11$kpc for the RGB. From the RGB sample, we find from galactocentric distances larger than $10$kpc the line-of-nodes twists away from the Sun-anticentre line towards galactic azimuths $\sim 180-200^{\circ}$ increasing with radius, though possibly influenced by extinction. The RGB sample reveals a slightly lopsided stellar warp with $\sim 250$pc between the up and down sides. The line of maximum of proper motions in latitude is systematically offset from the line-of-nodes estimated from the spatial data, which our models predict as a kinematic signature of lopsidedness. We also show a prominent wave-like pattern of a bending mode different in the OB and RGB, and substructures that might not be related to the Galactic warp nor to a bending mode. GDR2 triggers the need for complex kinematic models, flexible enough to combine both wave-like patterns and an S-shaped lopsided warp.[abridged]

The purpose of this study is the identification of young ($1< age < 100$ Myr), nearby ($d \leqslant100$ pc) moving groups (YNMGs) through their kinematic signature. YNMGs could be the result of the recent dispersal of young embedded clusters, such that they still represent kinematically cold groups, carrying the residual motion of their parental cloud. Using the fact that a large number ($\sim$ 14000) of the RAVE sources with evidence of chromospheric activity, also present signatures of stellar youth, we selected a sample of solar type sources with the highest probability of chromospheric activity to look for common kinematics. We made use of radial velocity information from RAVE and astrometric parameters from GAIA DR2 to construct a 6-dimension position-velocity vector catalog for our full sample. We developed a method based on the grouping of stars with similar orientation of their velocity vectors, which we call the Cone Method Sampling. Using this method, we detected 646 sources with high significance in the velocity space, with respect to the average orientation of artificial distributions made from a purely Gaussian velocity ellipsoid with null vertex deviation. We compared this sample of highly significant sources with a catalog of YNMGs reported in previous studies, which yield 75 confirmed members. From the remaining sample, about 50% of the sources have ages younger than 100 Myr, which indicate they are highly probable candidates to be new members of identified or even other YNMGs in the solar neighborhood.

The Gaia astrometric mission may offer an unprecedented opportunity to discover new tidal streams in the Galactic halo. To test this, we apply nGC3, a great-circle-cell count method that combines position and proper motion data to identify streams, to ten mock Gaia catalogues of K giants and RR Lyrae stars constructed from cosmological simulations of Milky Way analogues. We analyse two sets of simulations, one using a combination of $N$-body and semi-analytical methods which has extremely high resolution, the other using hydro-dynamical methods, which captures the dynamics of baryons, including the formation of an in situ halo. These ten realisations of plausible Galactic merger histories allow us to assess the potential for the recovery of tidal streams in different Milky Way formation scenarios. We include the Gaia~selection function and observational errors in these mock catalogues. We find that the nGC3 method has a well-defined detection boundary in the space of stream width and projected overdensity, that can be predicted based on direct observables alone. We predict that about 4-13 dwarf galaxy streams can be detected in a typical Milky Way-mass halo with Gaia+nGC3, with an estimated efficiency of $>$80\% inside the detection boundary. The progenitors of these streams are in the mass range of the classical dwarf galaxies and may have been accreted as early as redshift $\sim3$. Finally, we analyse how different possible extensions of the Gaia mission will improve the detection of tidal streams.

We present a novel positive potential-density pair expansion for modelling galaxies, based on the Miyamoto-Nagai (MN) disc. By using three sets of such discs, each one of them aligned along each symmetry axis, we are able to reconstruct a broad range of potentials that correspond to density profiles from exponential discs to 3D power law models with varying triaxiality (henceforth simply "twisted" models). We increase the efficiency of our expansion by allowing the scale length parameter of each disc to be negative. We show that, for suitable priors on the scale length and height parameters, these "MNn discs" have just one negative density minimum. This allows us to ensure global positivity by demanding that the total density at the global minimum is positive. We find that at better than 10\% accuracy in our density reconstruction, we can represent a radial and vertical exponential disc over $0.1-10$ scale lengths/heights with 4 MNn discs, an NFW profile over $0.1-10$ scale lengths with 4 MNn discs, and a twisted triaxial NFW profile with 3 MNn discs per symmetry axis. Our expansion is efficient, fully analytic, and well-suited to reproducing the density distribution and gravitational potential of galaxies from discs to ellipsoids.

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

In this first paper we simulate the population of disc Red Clump stars to be observed by Gaia. We generate a set of test particles and we evolve it in a 3D barred Milky Way like galactic potential. We assign physical properties of the Red Clump trace population and a realistic 3D interstellar extinction model. We add Gaia observational constraints and an error model according to the pre-commissioning scientific performance assessments. We present and analyse two mock catalogues, offered to the community, that are an excellent test bed for testing tools being developed for the future scientific exploitation of Gaia data. The first catalogue contains stars up to Gaia G 20, while the second is the subset containing Gaia radial velocity data with a maximum error of sigmaVr=10 kms. Here we present first attempts to characterise the density structure of the Galactic bar in the Gaia space of observables. The Gaia large errors in parallax and the high interstellar extinction in the inner parts of the Galactic disc prevent us to model the bar overdensity. This result suggests the need to combine Gaia and IR data to undertake such studies. We find that IR photometric distances for this Gaia sample allow us to recover the Galactic bar orientation angle with an accuracy of approximately 5 degrees.

We explore the possibility of detecting and characterizing the warp of the stellar disc of our Galaxy using synthetic Gaia data. The availability of proper motions and, for the brightest stars radial velocities, adds a new dimension to this study. A family of Great Circle Cell Counts (GC3) methods is used. They are ideally suited to find the tilt and twist of a collection of rings, which allow us to detect and measure the warp parameters. To test them, we use random realizations of test particles which evolve in a realistic Galactic potential warped adiabatically to various final configurations. In some cases a twist is introduced additionally. The Gaia selection function, its errors model and a realistic 3D extinction map are applied to mimic three tracer populations: OB, A and Red Clump stars. We show how the use of kinematics improves the accuracy in the recovery of the warp parameters. The OB stars are demonstrated to be the best tracers determining the tilt angle with accuracy better than $\sim0.5$ up to Galactocentric distance of $\sim16$ kpc. Using data with good astrometric quality, the same accuracy is obtained for A type stars up to $\sim13$ kpc and for Red Clump up to the expected stellar cut-off. Using OB stars the twist angle is recovered to within $< 3^\circ$ for all distances.

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

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

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.

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

We present here results of numerical calculations demonstrating that the very long periodic variability detected in the Cataclysmic Variable FS Aur can be the result of eccentricity modulation of a close binary (CB) orbit induced by the presence of a third body on a circumbinary orbit. A third component with a substellar mass on a circular, relatively close orbit, modulates the mass transfer rate of the binary on much longer time scales than periods within the triple system. We also report here preliminary results of X-ray observations of FS Aur, providing further evidence that it contains magnetic and freely precessing white dwarf. These two findings allow us to incorporate the new and previously stressed hypothesis on the nature of FS Aur into one consistent model.

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.

We address the problem of identifying remnants of satellite galaxies in the halo of our galaxy with Gaia data. The remnants have to be extracted from a very large data set (of order 10^9 stars) in the presence of observational errors and against a background population of Galactic stars. We address this issue through a numerical simulation of the Gaia catalogue with a realistic number of entries. The Galaxy model includes separate light distributions and kinematics for the bulge, disc and stellar halo. Our catalogue contains 3.5x10^8 stars. We perform tree code 10^6-body simulations of satellite dwarf galaxies in orbit around a rigid mass model of the Galaxy. The shrinking satellite cores and tidal tails are added to the Monte Carlo simulation of the Gaia catalogue. When combining the Milky Way and dwarf galaxy models we include the complication that the luminosity function of the satellite is probed at various depths as a function of position along the tidal tails. The combined models are converted to a synthetic Gaia catalogue using a detailed model for the expected astrometric and radial velocity errors, depending on magnitude, colour and sky position of the stars. We explore the feasibility of detecting tidal streams in the halo using the energy versus angular momentum plane. A straightforward search in this plane will be very challenging. The background population and the observational errors make it difficult to detect tidal streams as discrete structures in the E-Lz plane. The propagation of observational errors leads to apparent caustic structures in the integrals of motion space that may be mistaken for physical entities. Any practical search strategy will have to use a combination of pre-selection of high-quality data and complementary searches using the photometric data from Gaia. (abridged)

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.

We find a deficiency of binary stars moving on strongly retrograde Galactic orbits. No binary deficiencies are seen for U or W velocities, however. From theoretical analyses, we rule out preferential disruption of pre-existing binary stars due to encounters with massive perturbers. We also rule out globular clusters as the source of the effect since prograde motions are more likely to create such an effect. We search for star streams and find one possible candidate, but it is not on a retrograde orbit and probably represents the remains of a cluster that has passed too near the Galactic center. Based on a very small number of stars, we find that about the right fraction of stars on retrograde Galactic orbits share some chemical similarities to the cluster omega Cen, suggesting that its parent galaxy could be the explanation.

We address the problem of identifying remnants of satellite galaxies in the halo of our galaxy with Gaia data. We make use of N-body simulations of dwarf galaxies being disrupted in the halo of our galaxy combined with a Monte Carlo model of the Milky Way galaxy. The models are converted to a simulated Gaia catalogue containing a realistic number ($\sim10^8$--$10^9$) of stars. The simulated catalogue can be used to study how to handle the large data set that Gaia will provide and to study issues such as how to best retrieve information on substructure in the Galactic halo. The techniques described are applicable to any set of N-body simulations of (parts of) the Galaxy.

A new method to identify coherent structures in velocity space --- moving groups --- in astrometric catalogues is presented: the Spaghetti method. It relies on positions, parallaxes, and proper motions and is ideally suited to search for moving groups in the Hipparcos Catalogue. No radial velocity information is required. The method has been tested extensively on synthetic data, and applied to the Hipparcos measurements for the Hyades and IC2602 open clusters. The resulting lists of members agree very well with those of Perryman et al. for the Hyades and of Whiteoak and Braes for IC2602.

We present a new procedure for the identification of moving groups. It is a combination of two independent member selection methods. One is a modern implementation of a classical convergent point method for proper motion data. The other objectively identifies moving groups in velocity space using proper motions and parallaxes. We briefly describe both methods, and illustrate their power by applying them to Hipparcos measurements of a field containing the Upper Scorpius subgroup of the Scorpio-Centaurus association (Sco OB2). We show how our membership selection procedure not only improves the list of previously known B and A-type members, but also identifies many new members, including a significant number of F stars. We apply our procedure to other nearby OB associations elsewhere in these proceedings (Hoogerwerf et al.; de Zeeuw et al.).