Showing 1–38 of 38 results for author: Mateu, C

Structure, kinematics and time evolution of the Galactic Warp from Classical Cepheids

Authors: Mauro Cabrera-Gadea, Cecilia Mateu, Pau Ramos, Mercè Romero-Gómez, Teresa Antoja, Luis Aguilar, ., UdelaR, Uruguay., National Astronomical Observatory of Japan, ICCUB, FQA, IEEC, IA-UNAM Ensenada

Abstract: The warp is a well-known undulation of the Milky Way disc. Its structure has been widely studied, but only since Gaia DR2 has it been possible to reveal its kinematic signature beyond the solar neighbourhood. In this work we present an analysis of the warp traced by Classical Cepheids by means of a Fourier decomposition of their height ($Z$) and, for the first time, of their vertical velocity (… ▽ More The warp is a well-known undulation of the Milky Way disc. Its structure has been widely studied, but only since Gaia DR2 has it been possible to reveal its kinematic signature beyond the solar neighbourhood. In this work we present an analysis of the warp traced by Classical Cepheids by means of a Fourier decomposition of their height ($Z$) and, for the first time, of their vertical velocity ($V_z$). We find a clear but complex signal that in both variables reveals an asymmetrical warp. In $Z$ we find the warp to be almost symmetric in amplitude at the disc's outskirts, with the two extremes never being diametrically opposed at any radius and the line of nodes presenting a twist in the direction of stellar rotation for $R>11$ kpc. For $V_z$, in addition to the usual $m=1$ mode, an $m=2$ mode is needed to represent the kinematic signal of the warp, reflecting its azimuthal asymmetry. The line of maximum vertical velocity is similarly twisted as the line of nodes and trails behind by $\approx 25^\circ$. We develop a new formalism to derive the pattern speed and change in amplitude with time $\dot{A}$ of each Fourier mode at each radius, via a joint analysis of the Fourier decomposition in $Z$ and $V_z$. By applying it to the Cepheids we find, for the $m=1$ mode, a constant pattern speed in the direction of stellar rotation of $9.2\pm3.1$ km/s/kpc, a negligible $\dot{A}$ up to $R\approx 14$ kpc and a slight increase at larger radii, in agreement with previous works. △ Less

Submitted 24 January, 2024; originally announced January 2024.

Comments: Accepted for publication in MNRAS. 15 pages, 10 figures, 4 tables (+ Appendices). Full version of tables 1 and 2 included

Fusing Feature Engineering and Deep Learning: A Case Study for Malware Classification

Authors: Daniel Gibert, Carles Mateu, Jordi Planes, Quan Le

Abstract: Machine learning has become an appealing signature-less approach to detect and classify malware because of its ability to generalize to never-before-seen samples and to handle large volumes of data. While traditional feature-based approaches rely on the manual design of hand-crafted features based on experts knowledge of the domain, deep learning approaches replace the manual feature engineering p… ▽ More Machine learning has become an appealing signature-less approach to detect and classify malware because of its ability to generalize to never-before-seen samples and to handle large volumes of data. While traditional feature-based approaches rely on the manual design of hand-crafted features based on experts knowledge of the domain, deep learning approaches replace the manual feature engineering process by an underlying system, typically consisting of a neural network with multiple layers, that perform both feature learning and classification altogether. However, the combination of both approaches could substantially enhance detection systems. In this paper we present an hybrid approach to address the task of malware classification by fusing multiple types of features defined by experts and features learned through deep learning from raw data. In particular, our approach relies on deep learning to extract N-gram like features from the assembly language instructions and the bytes of malware, and texture patterns and shapelet-based features from malwareś grayscale image representation and structural entropy, respectively. These deep features are later passed as input to a gradient boosting model that combines the deep features and the hand-crafted features using an early-fusion mechanism. The suitability of our approach has been evaluated on the Microsoft Malware Classification Challenge benchmark and results show that the proposed solution achieves state-of-the-art performance and outperforms gradient boosting and deep learning methods in the literature. △ Less

Submitted 12 June, 2022; originally announced June 2022.

galstreams: A Library of Milky Way Stellar Stream Footprints and Tracks

Authors: Cecilia Mateu

Abstract: Nearly a hundred stellar streams have been found to date around the Milky Way and the number keeps growing at an ever faster pace. Here we present the galstreams library, a compendium of angular position, distance, proper motion and radial velocity track data for nearly a hundred (95) Galactic stellar streams. The information published in the literature has been collated and homogenised in a consi… ▽ More Nearly a hundred stellar streams have been found to date around the Milky Way and the number keeps growing at an ever faster pace. Here we present the galstreams library, a compendium of angular position, distance, proper motion and radial velocity track data for nearly a hundred (95) Galactic stellar streams. The information published in the literature has been collated and homogenised in a consistent format and used to provide a set of features uniformly computed throughout the library: e.g. stream length, end points, mean pole, stream's coordinate frame, polygon footprint, and pole and angular momentum tracks. We also use the information compiled to analyse the distribution of several observables across the library and to assess where the main deficiencies are found in the characterisation of individual stellar streams, as a resource for future follow-up efforts. The library is intended to facilitate keeping track of new discoveries and to encourage the use of automated methods to characterise and study the ensemble of known stellar streams by serving as a starting point. The library is publicly available as a Python package and served at the galstreams GitHub repository. △ Less

Submitted 26 January, 2023; v1 submitted 21 April, 2022; originally announced April 2022.

Comments: Accepted for publication at MNRAS. 29 pages, 19 figures, 2 tables, updated to match final version (new figures added! conclusions unchanged). The library of stellar stream tracks is publicly available at https://github.com/cmateu/galstreams

The Sagittarius stream in Gaia eDR3 and the origin of the bifurcations

Authors: P. Ramos, T. Antoja, Z. Yuan, A. Arentsen, P. -A. Oria, B. Famaey, R. Ibata, C. Mateu, J. A. Carballo-Bello

Abstract: The Sagittarius dwarf spheroidal (Sgr) is a dissolving galaxy being tidally disrupted by the Milky Way (MW). Its stellar stream still poses serious modelling challenges, which hinders our ability to use it effectively as a probe of the MW gravitational potential. Our goal is to construct the best possible sample of stars with which we can advance our understanding of the Sgr-MW interaction, focusi… ▽ More The Sagittarius dwarf spheroidal (Sgr) is a dissolving galaxy being tidally disrupted by the Milky Way (MW). Its stellar stream still poses serious modelling challenges, which hinders our ability to use it effectively as a probe of the MW gravitational potential. Our goal is to construct the best possible sample of stars with which we can advance our understanding of the Sgr-MW interaction, focusing on the characterisation of the bifurcations. We improve on previous methods based on the use of the wavelet transform to systematically search for the kinematic signature of the Sgr stream in the Gaia data. We refine our selection via the use of a clustering algorithm on the statistical properties of the colour-magnitude diagrams. Our final sample contains > 700k candidate stars, 3x larger than previous samples. We have been able to detect the bifurcation of the stream in both the northern and southern hemispheres, requiring four branches to fully describe this system. We present the detailed proper motion distribution of the trailing arm as a function of Lambda showing the presence of a sharp edge (on the side of the small proper motions) beyond which there are no Sgr stars. We also characterise the correlation between kinematics and distance. The chemical analysis of our sample shows a significant difference between the faint and bright branches. We provide analytical descriptions for the proper motion trends as well as for the sky distribution of the four branches of the stream. We interpret the bifurcations as the misaligned overlap of the material stripped at the antepenultimate pericentre (faint branches) with the stars ejected at the penultimate pericentre (bright branch). The source of this misalignment is still unknown but we argue that models with some internal rotation in the progenitor, at least during the time of stripping of the stars in the faint branches, are worth exploring. △ Less

Submitted 3 December, 2021; originally announced December 2021.

Comments: Submitted to A&A on the 03/12/2021. Online content available upon publication

Journal ref: A&A 666, A64 (2022)

Enhancing the Insertion of NOP Instructions to Obfuscate Malware via Deep Reinforcement Learning

Authors: Daniel Gibert, Matt Fredrikson, Carles Mateu, Jordi Planes, Quan Le

Abstract: Current state-of-the-art research for tackling the problem of malware detection and classification is centered on the design, implementation and deployment of systems powered by machine learning because of its ability to generalize to never-before-seen malware families and polymorphic mutations. However, it has been shown that machine learning models, in particular deep neural networks, lack robus… ▽ More Current state-of-the-art research for tackling the problem of malware detection and classification is centered on the design, implementation and deployment of systems powered by machine learning because of its ability to generalize to never-before-seen malware families and polymorphic mutations. However, it has been shown that machine learning models, in particular deep neural networks, lack robustness against crafted inputs (adversarial examples). In this work, we have investigated the vulnerability of a state-of-the-art shallow convolutional neural network malware classifier against the dead code insertion technique. We propose a general framework powered by a Double Q-network to induce misclassification over malware families. The framework trains an agent through a convolutional neural network to select the optimal positions in a code sequence to insert dead code instructions so that the machine learning classifier mislabels the resulting executable. The experiments show that the proposed method significantly drops the classification accuracy of the classifier to 56.53% while having an evasion rate of 100% for the samples belonging to the Kelihos_ver3, Simda, and Kelihos_ver1 families. In addition, the average number of instructions needed to mislabel malware in comparison to a random agent decreased by 33%. △ Less

Submitted 18 November, 2021; originally announced November 2021.

Report number: 0167-4048

Journal ref: Journal Computers & Security, Volume 113, 2022, 102543

Phenotypic-dependent variability and the emergence of tolerance in bacterial populations

Authors: José Camacho Mateu, Matteo Sireci, Miguel A. Muñoz

Abstract: Ecological and evolutionary dynamics have been historically regarded as unfolding at broadly separated timescales. However, these two types of processes are nowadays well documented to much more tightly than traditionally assumed, especially in communities of microorganisms. With this motivation in mind, here we scrutinize recent experimental results showing evidence of rapid evolution of toleranc… ▽ More Ecological and evolutionary dynamics have been historically regarded as unfolding at broadly separated timescales. However, these two types of processes are nowadays well documented to much more tightly than traditionally assumed, especially in communities of microorganisms. With this motivation in mind, here we scrutinize recent experimental results showing evidence of rapid evolution of tolerance by lag in bacterial populations that are periodically exposed to antibiotic stress in laboratory conditions. In particular, the distribution of single-cell lag times evolves its average value to approximately fit the antibiotic-exposure time. Moreover, the distribution develops right skewed heavy tails, revealing the presence of individuals with anomalously large lag times. Here, we develop a parsimonious individual-based model mimicking the actual demographic processes of the experimental setup. Individuals are characterized by a single phenotypic trait: their intrinsic lag time, which is transmitted with variation to the progeny. The model (in a version in which the amplitude of phenotypic variations grows with the parentś lag time) is able to reproduce quite well the key empirical observations. Furthermore, we develop a general mathematical framework allowing us to describe with good accuracy the properties of the stochastic model by means of a macroscopic equation, which generalizes the Crow Kimura equation in population genetics. From a broader perspective, this work represents a benchmark for the mathematical framework designed to tackle much more general eco-evolutionary problems, thus paving the road to further research avenues. △ Less

Submitted 6 October, 2021; v1 submitted 17 September, 2021; originally announced September 2021.

The RR Lyrae Delay-Time Distribution: A Novel Perspective on Models of Old Stellar Populations

Authors: Sumit K. Sarbadhicary, Mairead Heiger, Carles Badenes, Cecilia Mateu, Jeffrey Newman, Robin Ciardullo, Na'ama Hallakoun, Dan Maoz, Laura Chomiuk

Abstract: The delay-time distribution (DTD) is the occurrence rate of a class of objects as a function of time after a hypothetical burst of star formation. DTDs are mainly used as a statistical test of stellar evolution scenarios for supernova progenitors, but they can be applied to many other classes of astronomical objects. We calculate the first DTD for RR Lyrae variables using 29,810 RR Lyrae from the… ▽ More The delay-time distribution (DTD) is the occurrence rate of a class of objects as a function of time after a hypothetical burst of star formation. DTDs are mainly used as a statistical test of stellar evolution scenarios for supernova progenitors, but they can be applied to many other classes of astronomical objects. We calculate the first DTD for RR Lyrae variables using 29,810 RR Lyrae from the OGLE-IV survey and a map of the stellar-age distribution (SAD) in the Large Magellanic Cloud (LMC). We find that $\sim 46\%$ of the OGLE-IV RR Lyrae are associated with delay-times older than 8 Gyr (main-sequence progenitor masses less than 1 M$_{\odot}$), and consistent with existing constraints on their ages, but surprisingly about $51\%$ of RR Lyrae appear have delay times $1.2-8$ Gyr (main-sequence masses between $1 - 2$ M$_{\odot}$ at LMC metallicity). This intermediate-age signal also persists outside the Bar-region where crowding is less of a concern, and we verified that without this signal, the spatial distribution of the OGLE-IV RR Lyrae is inconsistent with the SAD map of the LMC. Since an intermediate-age RR Lyrae channel is in tension with the lack of RR Lyrae in intermediate-age clusters (noting issues with small-number statistics), and the age-metallicity constraints of LMC stars, our DTD result possibly indicates that systematic uncertainties may still exist in SAD measurements of old-stellar populations, perhaps stemming from the construction methodology or the stellar evolution models used. We described tests to further investigate this issue. △ Less

Submitted 27 January, 2021; originally announced January 2021.

Comments: 21 pages, 11 figures. Accepted to ApJ

Revealing the Structure and Internal Rotation of the Sagittarius Dwarf Spheroidal Galaxy with Gaia and Machine Learning

Authors: Andrés del Pino, Mark A. Fardal, Roeland P. van der Marel, Ewa L. Łokas, Cecilia Mateu, Sangmo Tony Sohn

Abstract: We present a detailed study of the internal structure and kinematics of the core of the Sagittarius dwarf spheroidal galaxy (Sgr). Using machine-learning techniques, we have combined the information provided by 3300 RR Lyrae stars, more than 2000 spectroscopically observed stars, and the Gaia second data release to derive the full phase space, i.e. 3D positions and kinematics, of more than… ▽ More We present a detailed study of the internal structure and kinematics of the core of the Sagittarius dwarf spheroidal galaxy (Sgr). Using machine-learning techniques, we have combined the information provided by 3300 RR Lyrae stars, more than 2000 spectroscopically observed stars, and the Gaia second data release to derive the full phase space, i.e. 3D positions and kinematics, of more than $1.2\times10^5$ member stars in the core of the galaxy. Our results show that Sgr has a bar structure $\sim 2.5$ kpc long, and that tidal tails emerge from its tips to form what it is known as the Sgr stream. The main body of the galaxy, strongly sheared by tidal forces, is a triaxial (almost prolate) ellipsoid with its longest principal axis of inertia inclined $43^\circ\pm6^\circ$ with respect to the plane of the sky and axis ratios of 1:0.67:0.60. Its external regions are expanding mainly along its longest principal axis, yet the galaxy conserves an inner core of about $500\times330\times300$ pc that shows no net expansion and is rotating at $v_{\rm rot} = 4.13 \pm 0.16$ ${\rm{ km \ s^{-1}}}$. The internal angular momentum of Sgr forms an angle $θ= 18^\circ\pm6^\circ$ with respect to its orbital angular momentum, meaning that the galaxy is in an inclined prograde orbit around the Milky Way. We compared our results with predictions from $N$-body models with spherical, pressure-supported progenitors and a model whose progenitor is a flattened rotating disk. Only the rotating model, based on preexisting simulations aimed at reproducing the line-of-sight velocity gradients observed in Sgr, was able to reproduce the observed properties in the core of the galaxy. △ Less

Submitted 15 February, 2021; v1 submitted 4 November, 2020; originally announced November 2020.

Comments: 37 pages, 26 figures, 2 tables. Accepted for publication in the Astrophysical Journal

The outer disc in shambles: blind detection of Monoceros and ACS with Gaia's astrometric sample

Authors: P. Ramos, T. Antoja, C. Mateu, F. Anders, C. F. P. Laporte, J. A. Carballo-Bello, B. Famaey, R. Ibata

Abstract: The astrometric sample of Gaia allows us to study the outermost Galactic disc, the halo and their interface. It is precisely at the very edge of the disc where the effects of external perturbations are expected to be the most noticeable. Our goal is to detect the kinematic substructure present in the halo and at the edge of the Milky Way (MW) disc, and provide observational constraints on their ph… ▽ More The astrometric sample of Gaia allows us to study the outermost Galactic disc, the halo and their interface. It is precisely at the very edge of the disc where the effects of external perturbations are expected to be the most noticeable. Our goal is to detect the kinematic substructure present in the halo and at the edge of the Milky Way (MW) disc, and provide observational constraints on their phase-space distribution. We download, one HEALpix at a time, the proper motion histogram of distant stars, to which we apply a Wavelet Transformation to reveal the significant overdensities. We then analyse the large coherent structures that appear in the sky. We reveal a sharp yet complex anticentre dominated by Monoceros (MNC) and the Anticentre Stream (ACS) in the north, which we find with an intensity comparable to the Magellanic clouds and the Sagittarius stream, and by MNC south and TriAnd at negative latitudes. Our method allows us to perform a morphological analysis of MNC and ACS, both spanning more than 100$^\circ$ in longitude, and to provide a high purity sample of giants with which we track MNC down to latitudes as low as $\sim$5$^\circ$. Their colour-magnitude diagram is consistent with extended structures at a distance of $\sim$10-11 kpc originated in the disc, with a very low ratio of RR Lyrae over M giants, and kinematics compatible with the rotation curve at those distances or only slightly slower. We present a precise characterisation of MNC and ACS, two previously known structures that our method reveals naturally, allowing us to detect them without limiting ourselves to a particular stellar type and, for the first time, using only kinematics. Our results allow future studies to model their chemo-dynamics and evolution, thus constraining some of the most influential processes that shaped the MW. △ Less

Submitted 2 November, 2020; originally announced November 2020.

Comments: Submitted to Astronomy and Astrophysics. Online material available upon acceptance

Journal ref: A&A 646, A99 (2021)

Discovery of Extended Tidal Tails around the Globular Cluster Palomar 13

Authors: Nora Shipp, Adrian Price-Whelan, Kiyan Tavangar, Cecilia Mateu, Alex Drlica-Wagner

Abstract: We use photometry from the DECam Legacy Survey to detect candidate tidal tails extending ~5 deg on either side of the Palomar 13 globular cluster. The tails are aligned with the proper motion of Palomar 13 and are consistent with its old, metal-poor stellar population. We identify three RR Lyrae stars that are plausibly associated with the tails, in addition to four previously known in the cluster… ▽ More We use photometry from the DECam Legacy Survey to detect candidate tidal tails extending ~5 deg on either side of the Palomar 13 globular cluster. The tails are aligned with the proper motion of Palomar 13 and are consistent with its old, metal-poor stellar population. We identify three RR Lyrae stars that are plausibly associated with the tails, in addition to four previously known in the cluster. From these RR Lyrae stars, we find that the mean distance to the cluster and tails is $23.6 \pm 0.2$ kpc and estimate the total (initial) luminosity of the cluster to be $L_V=5.1^{+9.7}_{-3.4}\times 10^3 L_\odot$, consistent with previous claims that its initial luminosity was higher than its current luminosity. Combined with previously-determined proper motion and radial velocity measurements of the cluster, we find that Palomar 13 is on a highly eccentric orbit ($e\sim 0.8$) with a pericenter of ~9 kpc and an apocenter of ~69 kpc, and a recent pericentric passage of the cluster ~75 Myr ago. We note a prominent linear structure in the interstellar dust map that runs parallel to the candidate tidal features, but conclude that reddening due to dust is unlikely to account for the structure that we observe. If confirmed, the Palomar 13 stellar stream would be one of very few streams with a known progenitor system, making it uniquely powerful for studying the disruption of globular clusters, the formation of the stellar halo, and the distribution of matter within our Galaxy. △ Less

Submitted 22 June, 2020; originally announced June 2020.

Comments: 8 pages, 4 figures, 1 table. Submitted to AAS journals

Empirical completeness assessment of the Gaia DR2, Pan-STARRS 1 and ASAS-SN-II RR Lyrae catalogues

Authors: Cecilia Mateu, Berry Holl, Joris de Ridder, Lorenzo Rimoldini

Abstract: RR Lyrae stars are an important and widely used tracer of the most ancient populations of our Galaxy, mainly due to their standard candle nature. The availability of large scale surveys of variable stars is allowing us to trace the structure of our entire Galaxy, even in previously inaccessible areas like the Galactic disc. In this work we aim to provide an empirical assessment of the completeness… ▽ More RR Lyrae stars are an important and widely used tracer of the most ancient populations of our Galaxy, mainly due to their standard candle nature. The availability of large scale surveys of variable stars is allowing us to trace the structure of our entire Galaxy, even in previously inaccessible areas like the Galactic disc. In this work we aim to provide an empirical assessment of the completeness of the three largest RR Lyrae catalogues available: Gaia DR2, PanSTARRS-1 and ASAS-SN-II. Using a joint probabilistic analysis of the three surveys we compute 2D and 3D completeness maps in each survey's full magnitude range. At the bright end (G<13) we find ASAS-SN-II and Gaia are near 100% complete in RRab at high latitude (|b|>20deg); ASAS-SN-II has the best completeness at low latitude for RRab and at all latitudes for RRc. At the faint end (G>13), Gaia DR2 is the most complete catalogue for both RR Lyrae types, at any latitude, with median completeness rates of 95% (RRab) and >85% (RRc) outside the ecliptic plane (|β|>25deg). We confirm a high and uniform completeness of PanSTARRS-1 RR Lyrae at 91% (RRab) and 82% (RRc) down to G~18, and provide the first estimate of its completeness at low galactic latitude (|b|<20deg) at an estimated median 65% (RRab) and 50-60% (RRc). Our results are publicly available as 2D and 3D completeness maps, and as functions to evaluate each survey's completeness versus distance or per line-of sight. △ Less

Submitted 16 June, 2020; originally announced June 2020.

Comments: 15 pages, 14 figures, 4+10 tables. Accepted for publication at MNRAS. Completeness maps and Python utils to plot and/or recompute them, available at https://github.com/cmateu/rrl_completeness

Full 5D characterisation of the Sagittarius stream with Gaia DR2 RR Lyrae

Authors: P. Ramos, C. Mateu, T. Antoja, A. Helmi, A. Castro-Ginard, E. Balbinot, J. M. Carrasco

Abstract: The Sagittarius stream is one of the best tools that we currently have to estimate the mass and shape of our Galaxy. However, assigning membership and obtaining the phase-space distribution of the stars that form the tails is quite challenging. Our goal is to produce a catalogue of RR Lyrae stars of Sagittarius and obtain an empiric measurement of the trends along the stream in sky position, dista… ▽ More The Sagittarius stream is one of the best tools that we currently have to estimate the mass and shape of our Galaxy. However, assigning membership and obtaining the phase-space distribution of the stars that form the tails is quite challenging. Our goal is to produce a catalogue of RR Lyrae stars of Sagittarius and obtain an empiric measurement of the trends along the stream in sky position, distance and tangential velocities. We generate two initial samples from the Gaia DR2 RR Lyrae catalogue: one, selecting only the stars within \pm20deg of the orbital plane of Sagittarius (Strip) and the other, the result of applying the Pole Count Map (nGC3) algorithm. We then use the model-independent, deterministic method developed in this work to remove most of the contamination by detecting and isolating the stream in distance and proper motions. The output is two empiric catalogues: the Strip sample (higher-completeness, lower-purity) which contains 11 677 stars, and the nGC3 sample (higher-purity, lower-completeness) with 6 608 stars. We characterise the changes along the stream in all the available dimensions, the 5 astrometric ones plus the metallicity, covering more than 2pi rad in the sky and obtain new estimates for the apocentres and the mean [Fe/H] of the RR Lyrae population. Also, we show the first map of the two components of the tangential velocity, thanks to the combination of distances and proper motions. Finally, we detect the bifurcation in the leading arm and report no significant difference between the two branches, either in metallicity, kinematics or distance. We provide the largest sample of RR Lyrae candidates of Sagittarius, which can be used as an input for a spectroscopic follow-up or as a reference for the new generation of models of the stream through the interpolators in distance and velocity that we have constructed. △ Less

Submitted 23 July, 2020; v1 submitted 25 February, 2020; originally announced February 2020.

Comments: 19 pages. Published in Astronomy and Astrophysics. Related on-line material available at https://services.fqa.ub.edu/sagittarius

Journal ref: A&A 638, A104 (2020)

An all-sky proper motion map of the Sagittarius stream using Gaia DR2

Authors: T. Antoja, P. Ramos, C. Mateu, A. Helmi, F. Anders, C. Jordi, J. A. Carballo-Bello

Abstract: We aim to measure the proper motion along the Sagittarius stream that is the missing piece to determine its full 6D phase space coordinates. We conduct a blind search of over-densities in proper motion from Gaia DR2 in a broad region around the Sagittarius stream by applying wavelet transform techniques. We find that for most of the sky patches, the highest intensity peaks delineate the path of th… ▽ More We aim to measure the proper motion along the Sagittarius stream that is the missing piece to determine its full 6D phase space coordinates. We conduct a blind search of over-densities in proper motion from Gaia DR2 in a broad region around the Sagittarius stream by applying wavelet transform techniques. We find that for most of the sky patches, the highest intensity peaks delineate the path of the Sagittarius stream. The 1500 peaks identified depict a continuous sequence spanning almost $2π$ in the sky, only obscured when the stream crosses the Galactic disk. Altogether, around $100\,000$ stars potentially belong to the stream as indicated by a coarse inspection of the colour-magnitude diagrams. From these stars, we determine the proper motion along the Sagittarius stream, making it the proper motion sequence with the largest span and continuity ever measured for a stream. A first comparison with existing N-body models of the stream reveals some discrepancies, especially near the pericentre of the trailing arm and an overestimation of the total proper motion for the leading arm. Our study can be the starting point for determining the variation of the population of stars along the stream, the distance to the stream with red clump stars, and the solar motion. It will also allow a much better measurement of the Milky Way potential. △ Less

Submitted 31 January, 2020; v1 submitted 27 January, 2020; originally announced January 2020.

Comments: Accepted for publication in A&A Letters. See additional material at https://services.fqa.ub.edu/sagittarius, v2 with minor editing

Journal ref: A&A 635, L3 (2020)

Kinematics of the Palomar 5 stellar stream from RR Lyrae stars

Authors: Adrian M. Price-Whelan, Cecilia Mateu, Giuliano Iorio, Sarah Pearson, Ana Bonaca, Vasily Belokurov

Abstract: Thin stellar streams, formed from the tidal disruption of globular clusters, are important gravitational tools, sensitive to both global and small-scale properties of dark matter. The Palomar 5 stellar stream (Pal 5) is an exemplar stream within the Milky Way: Its $\sim 20^\circ$ tidal tails connect back to the progenitor cluster, and the stream has been used to study the shape, total mass, and su… ▽ More Thin stellar streams, formed from the tidal disruption of globular clusters, are important gravitational tools, sensitive to both global and small-scale properties of dark matter. The Palomar 5 stellar stream (Pal 5) is an exemplar stream within the Milky Way: Its $\sim 20^\circ$ tidal tails connect back to the progenitor cluster, and the stream has been used to study the shape, total mass, and substructure fraction of the dark matter distribution of the Galaxy. However, most details of the phase-space distribution of the stream are not fully explained, and dynamical models that use the stream for other inferences are therefore incomplete. Here we aim to measure distance and kinematic properties along the Pal 5 stream in order to motivate improved models of the system. We use a large catalog of RR Lyrae-type stars (RRLs) with astrometric data from the Gaia mission to probabilistically identify RRLs in the Pal 5 stream. RRLs are useful because they are intrinsically-luminous standard candles and their distances can be inferred with small relative precision ($\sim3\%$). By building a probabilistic model of the Pal 5 cluster and stream in proper motion and distance, we find 27 RRLs consistent with being members of the cluster (10) and stream (17). Using these RRLs, we detect gradients in distance and proper motion along the stream, and provide an updated measurement of the distance to the Pal 5 cluster using the RRLs, $d = 20.6 \pm 0.2~\textrm{kpc}$. We provide a catalog of Pal 5 RRLs with inferred membership probabilities for future modeling work. △ Less

Submitted 5 February, 2020; v1 submitted 1 October, 2019; originally announced October 2019.

Comments: 13 pages, 4 figures. Published in AJ

QUEST Eclipsing Binaries and ELLISA: Eclipsing binary search in the QUEST low latitude catalogue and the ELLISA light curve simulator

Authors: Bolivia Cuevas-Otahola, Cecilia Mateu, Fabiola Hernández-Pérez, Juan José Downes, A. Katherina Vivas, César Briceño

Abstract: The realistic simulation of variable star populations is fundamental to determine the selection function and contamination in existing and upcoming multi-epoch surveys. We present \ellisa, a simulator that produces an ensemble of mock light curves for a population of eclipsing binaries obtained from physical and orbital parameters consistent with different galactic populations, and which considers… ▽ More The realistic simulation of variable star populations is fundamental to determine the selection function and contamination in existing and upcoming multi-epoch surveys. We present \ellisa, a simulator that produces an ensemble of mock light curves for a population of eclipsing binaries obtained from physical and orbital parameters consistent with different galactic populations, and which considers user-supplied time sampling and photometric errors to represent any given survey. We carried out a search for eclipsing binaries in the QUEST low-galactic latitude catalogue of variable stars, spanning an area of 476 sq. deg at $-25^\circ \la b \la 30^\circ$ and $190^ \circ \leqslant l \leqslant 230^ \circ$ towards the galactic anti-centre, and use \ellisa~to characterise the completeness of the resulting catalogue in terms of amplitudes and periods of variation as well as eclipsing binary type. The resulting catalogue consists of $1,125$ eclipsing binaries, out of which $179$, $60$ and $886$ are EA, EB and EW types respectively. We estimate, on average, $30\%$ completeness in the period range $0.25 \la P/d \la 1$ for EB+EW binaries and $15\%$ completeness for EA binaries with periods $2 \la P/d \la 10$, being the time sampling the primary factor determining the completeness of each type of eclipsing binary. This is one of few eclipsing binary catalogues reported with an estimate of the selection function. Mock eclipsing binary light curve libraries produced with \ellisa~can be used to estimate the selection function and optimise eclipsing binary searches in upcoming multi-epoch surveys such as Gaia, the Panoramic Survey Telescope and Rapid Response System, the Zwicky Transient Factory or the Large Synoptic Survey Telescope. △ Less

Submitted 11 January, 2019; originally announced January 2019.

Comments: 17 pages, 13 figures

The total mass of the Large Magellanic Cloud from its perturbation on the Orphan stream

Authors: D. Erkal, V. Belokurov, C. F. P. Laporte, S. E. Koposov, T. S. Li, C. J. Grillmair, N. Kallivayalil, A. M. Price-Whelan, N. W. Evans, K. Hawkins, D. Hendel, C. Mateu, J. F. Navarro, A. del Pino, C. T. Slater, S. T. Sohn

Abstract: In a companion paper by Koposov et al., RR Lyrae from \textit{Gaia} Data Release 2 are used to demonstrate that stars in the Orphan stream have velocity vectors significantly misaligned with the stream track, suggesting that it has received a large gravitational perturbation from a satellite of the Milky Way. We argue that such a mismatch cannot arise due to any realistic static Milky Way potentia… ▽ More In a companion paper by Koposov et al., RR Lyrae from \textit{Gaia} Data Release 2 are used to demonstrate that stars in the Orphan stream have velocity vectors significantly misaligned with the stream track, suggesting that it has received a large gravitational perturbation from a satellite of the Milky Way. We argue that such a mismatch cannot arise due to any realistic static Milky Way potential and then explore the perturbative effects of the Large Magellanic Cloud (LMC). We find that the LMC can produce precisely the observed motion-track mismatch and we therefore use the Orphan stream to measure the mass of the Cloud. We simultaneously fit the Milky Way and LMC potentials and infer that a total LMC mass of $1.38^{+0.27}_{-0.24} \times10^{11}\,\rm{M_\odot}$ is required to bend the Orphan Stream, showing for the first time that the LMC has a large and measurable effect on structures orbiting the Milky Way. This has far-reaching consequences for any technique which assumes that tracers are orbiting a static Milky Way. Furthermore, we measure the Milky Way mass within 50 kpc to be $3.80^{+0.14}_{-0.11}\times10^{11} M_\odot$. Finally, we use these results to predict that, due to the reflex motion of the Milky Way in response to the LMC, the outskirts of the Milky Way's stellar halo should exhibit a bulk, upwards motion. △ Less

Submitted 14 May, 2019; v1 submitted 19 December, 2018; originally announced December 2018.

Comments: 17 pages, 11 figures. Updated to version accepted to MNRAS after minor revision

Piercing the Milky Way: an all-sky view of the Orphan Stream

Authors: S. E. Koposov, V. Belokurov, T. S. Li, C. Mateu, D. Erkal, C. J. Grillmair, D. Hendel, A. M. Price-Whelan, C. F. P. Laporte, K. Hawkins, S. T. Sohn, A. del Pino, N. W. Evans, C. T. Slater, N. Kallivayalil, J. F. Navarro

Abstract: We use astrometry, broad-band photometry and variability information from the Data Release 2 of ESA's Gaia mission (GDR2) to identify members of the Orphan Stream (OS) across the whole sky. The stream is traced above and below the celestial equator and in both Galactic hemispheres, thus increasing its visible length to ~ 210 degrees equivalent to ~150 kpc in physical extent. Taking advantage of th… ▽ More We use astrometry, broad-band photometry and variability information from the Data Release 2 of ESA's Gaia mission (GDR2) to identify members of the Orphan Stream (OS) across the whole sky. The stream is traced above and below the celestial equator and in both Galactic hemispheres, thus increasing its visible length to ~ 210 degrees equivalent to ~150 kpc in physical extent. Taking advantage of the large number of RR Lyrae stars in the OS, we extract accurate distances and proper motions across the entire stretch of the tidal debris studied. As delineated by the GDR2 RR Lyrae, the stream exhibits two prominent twists in its shape on the sky which are accompanied by changes in the tangential motion. We complement the RR Lyrae maps with those created using GDR2 Red Giants and the DECam Legacy Survey Main Sequence Turn-Off stars. The behavior of the OS track on the sky is consistent across all three tracers employed. We detect a strong non-zero motion in the across-stream direction for a substantial portion of the stream. Such a misalignment between the debris track and the streaming velocity cannot be reproduced in a static gravitational potential and signals an interaction with a massive perturber. △ Less

Submitted 11 February, 2019; v1 submitted 19 December, 2018; originally announced December 2018.

Comments: 18 pages; accepted to MNRAS

Gaia kinematics reveal a complex lopsided and twisted Galactic disc warp

Authors: M. Romero-Gómez, C. Mateu, L. Aguilar, F. Figueras, A. Castro-Ginard

Abstract: 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 dep… ▽ More 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] △ Less

Submitted 7 June, 2019; v1 submitted 18 December, 2018; originally announced December 2018.

Comments: 14 pages (+7 pages of appendix), matches the accepted version in A&A after referee comments (June 5th 2019)

Journal ref: A&A 627, A150 (2019)

A Gaia DR2 search for dwarf galaxies towards Fermi-LAT sources: implications for annihilating dark matter

Authors: Ioana Ciucă, Daisuke Kawata, Shin'ichiro Ando, Francesca Calore, Justin I. Read, Cecilia Mateu

Abstract: We make a first attempt to find dwarf galaxies in eight \Fermi-LAT extended, unassociated, source fields using \Gaia\ DR2. We probe previously unexplored heliocentric distances of $d<20$~kpc with an extreme-deconvolution (XD) technique. We find no signature of a dwarf galaxy in any of these fields despite \Gaia's excellent astrometric accuracy. We estimate our detection limits by applying the XD m… ▽ More We make a first attempt to find dwarf galaxies in eight \Fermi-LAT extended, unassociated, source fields using \Gaia\ DR2. We probe previously unexplored heliocentric distances of $d<20$~kpc with an extreme-deconvolution (XD) technique. We find no signature of a dwarf galaxy in any of these fields despite \Gaia's excellent astrometric accuracy. We estimate our detection limits by applying the XD method to mock data, obtaining a conservative limit on the stellar mass of $M_* < 10^4$~M$_{\sun}$ for $d < 20$\, kpc. Such a low stellar mass implies either a low-mass subhalo, or a massive stripped-down subhalo. We use an analytic model for stripped subhalos to argue that, given the sizes and fluxes of the \Fermi-LAT sources, we can reject the hypothesis that they owe to dark matter annihilation. △ Less

Submitted 20 July, 2018; v1 submitted 7 May, 2018; originally announced May 2018.

Comments: 9 pages, 6 figures, accepted for publication in MNRAS

The CIDA Variability Survey of Orion OB1 II: demographics of the young, low-mass stellar populations

Authors: Cesar Briceno, Nuria Calvet, Jesus Hernandez, Anna K. Vivas, Cecilia Mateu, Juan Jose Downes, Jaqueline Loerincs, Alice Perez-Blanco, Perry Berlind, Catherine Espaillat, Lori Allen, Lee Hartmann, Mario Mateo, John Bailey III

Abstract: We present results of our large scale, optical, multi-epoch photometric survey of ~180 square degrees across the Orion OB1 association, complemented with extensive follow up spectroscopy. We map and characterize in an uniform way the off-cloud, low-mass, pre-main sequence populations. We report 2064, mostly K and M-type, confirmed T Tauri members. Most (59%) are located in the OB1a subassociation,… ▽ More We present results of our large scale, optical, multi-epoch photometric survey of ~180 square degrees across the Orion OB1 association, complemented with extensive follow up spectroscopy. We map and characterize in an uniform way the off-cloud, low-mass, pre-main sequence populations. We report 2064, mostly K and M-type, confirmed T Tauri members. Most (59%) are located in the OB1a subassociation, 27% in OB1b, and 14% within the confines of the A and B molecular clouds. There is significant structure in the spatial distribution of the young stars. We characterize two new clusterings of T Tauri stars, HD 35762 and HR 1833 in the OB1a subassociation, and two stellar overdensities in OB1b. There is indication of two populations of young stars in the OB1b region, located at two different distances, which may be due to the OB1a subassociation overlapping on front of OB1b. The various groups and regions can be ordered in an age sequence that agrees with the long standing picture of star formation starting in Orion OB1a some 10-15 Myr ago. We define a new type of T Tauri star, the C/W class, objects we propose may be nearing the end of their accretion phase. We detect the observational signature of Li depletion in young K and M stars with a timescale of 8.5 Myr. The decline of the accretion fraction from ~2 - 10 Myr, implies an accretion e-folding timescale of 2.1 Myr. Finally, the median amplitude of the V-band variability shows the decline of stellar activity, from accreting Classical T Tauri stars to the least active field dwarfs. △ Less

Submitted 2 May, 2018; originally announced May 2018.

Comments: 41 pages, 30 figures, submitted to the Astronomical Journal on May 02, 2018. Abridged abstract

The Galactic Thick Disc density profile traced with RR Lyrae stars

Authors: Cecilia Mateu, A. Katherina Vivas

Abstract: We used a combination of public RR Lyrae star catalogs and a Bayesian methodology to derive robust structural parameters of the inner Halo (<25 kpc) and Thick Disc of the Milky Way. RR Lyrae stars are an unequivocal tracer of old metal-poor populations, for which accurate distances and extinctions can be individually estimated and so, are a reliable independent means of tracing the population of t… ▽ More We used a combination of public RR Lyrae star catalogs and a Bayesian methodology to derive robust structural parameters of the inner Halo (<25 kpc) and Thick Disc of the Milky Way. RR Lyrae stars are an unequivocal tracer of old metal-poor populations, for which accurate distances and extinctions can be individually estimated and so, are a reliable independent means of tracing the population of the old high-[α/Fe] disc usually associated to the Thick Disc. In particular, the chosen RR Lyrae sample spans regions at low galactic latitude toward the anti-center direction, allowing to probe the outermost parts of the disc. Our results favour a Thick Disc with short scale height and short scale length, $h_z=0.65_{-0.05}^{+0.09}$ kpc, $h_R=2.1_{-0.25}^{+0.82}$ kpc, for a model in which the inner Halo has a constant flattening of $q=0.90_{-0.03}^{+0.05}$ and a power law index of $n=-2.78_{-0.05}^{+0.05}$. Similar short scales for the Thick Disc are also found when considering an inner Halo with flattening dependent on radius. We also explored a model in which the Thick Disc has a flare and, although this is only mildly constrained by our data, a flare onset in the inner $\sim11$ kpc is highly disfavoured. △ Less

Submitted 22 May, 2018; v1 submitted 21 February, 2018; originally announced February 2018.

Comments: 14 pages, 8 figures. Updated to match MNRAS accepted version. No results or conclusions changed

Fourteen candidate RR Lyrae star streams in the inner Galaxy

Authors: Cecilia Mateu, Justin I. Read, Daisuke Kawata

Abstract: We apply the GC3 stream-finding method to RR Lyrae stars (RRLS) in the Catalina survey. We find two RRLS stream candidates at $>4σ$ confidence and another 12 at $>3.5σ$ confidence over the Galactocentric distance range $4 < D/{\rm kpc} < 26$. Of these, only two are associated with known globular clusters (NGC 1261 and Arp2). The remainder are candidate `orphan' streams, consistent with the idea th… ▽ More We apply the GC3 stream-finding method to RR Lyrae stars (RRLS) in the Catalina survey. We find two RRLS stream candidates at $>4σ$ confidence and another 12 at $>3.5σ$ confidence over the Galactocentric distance range $4 < D/{\rm kpc} < 26$. Of these, only two are associated with known globular clusters (NGC 1261 and Arp2). The remainder are candidate `orphan' streams, consistent with the idea that globular cluster streams are most visible close to dissolution. Our detections are likely a lower bound on the total number of dissolving globulars in the inner galaxy, since many globulars have few RRLS while only the brightest streams are visible over the Galactic RRLS background, particularly given the current lack of kinematical information. We make all of our candidate streams publicly available and provide a new GALSTREAMS Python library for the footprints of all known streams and overdensities in the Milky Way. △ Less

Submitted 10 November, 2017; originally announced November 2017.

Comments: 18 pages, 4 figures. Accepted for publication at MNRAS. GALSTREAMS Milky Way Streams Footprint Library are available at https://github.com/cmateu/galstreams . All RRL data and code used in the paper are available at https://cmateu.github.io/Cecilia_Mateu_WebPage/CatalinaGC3_Streams.html

The Astropy Problem

Authors: Demitri Muna, Michael Alexander, Alice Allen, Richard Ashley, Daniel Asmus, Ruyman Azzollini, Michele Bannister, Rachael Beaton, Andrew Benson, G. Bruce Berriman, Maciej Bilicki, Peter Boyce, Joanna Bridge, Jan Cami, Eryn Cangi, Xian Chen, Nicholas Christiny, Christopher Clark, Michelle Collins, Johan Comparat, Neil Cook, Darren Croton, Isak Delberth Davids, Éric Depagne, John Donor , et al. (129 additional authors not shown)

Abstract: The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical… ▽ More The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software. △ Less

Submitted 10 October, 2016; originally announced October 2016.

Predictions for the detection of Tidal Streams with Gaia using Great Circle Methods

Authors: Cecilia Mateu, Andrew Cooper, Andreea S. Font, Luis Aguilar, Carlos Frenk, Shaun Cole, Wenting Wang, Ian G. McCarthy

Abstract: 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 sim… ▽ More 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. △ Less

Submitted 31 May, 2017; v1 submitted 27 September, 2016; originally announced September 2016.

Comments: 25 pages, 15 figures, MNRAS in press. Updated to match accepted version, includes new discussion of distance distribution of recovered progenitors. Stream candidates and full mock catalogues available at https://cmateu.github.io/Cecilia_Mateu_WebPage/Gaia_Halo_Mocks.html

The rotation period distributions of 4--10 Myr T Tauri stars in Orion OB1: New constraints on pre-main-sequence angular momentum evolution

Authors: Md. Tanveer Karim, Keivan G. Stassun, Cesar Briceno, A. Katherina Vivas, Stefanie Raetz, Cecilia Mateu, Juan Jose Downes, Nuria Calvet, Jesus Hernandez, Ralph Neuhauser, Markus Mugrauer, Hidenori Takahashi, Kengo Tachihara, Rolf Chini, Gustavo A. Cruz-Dias, Alicia Aarnio, David J. James, Moritz Hackstein

Abstract: Most existing studies of the angular momentum evolution of young stellar populations have focused on the youngest (1-3 Myr) T Tauri stars. In contrast, the angular momentum distributions of older T Tauri stars (4-10 Myr) have been less studied, even though they hold key insight to understanding stellar angular momentum evolution at a time when protoplanetary disks have largely dissipated and when… ▽ More Most existing studies of the angular momentum evolution of young stellar populations have focused on the youngest (1-3 Myr) T Tauri stars. In contrast, the angular momentum distributions of older T Tauri stars (4-10 Myr) have been less studied, even though they hold key insight to understanding stellar angular momentum evolution at a time when protoplanetary disks have largely dissipated and when models therefore predict changes in the rotational evolution that can in principle be tested. We present a study of photometric variability among 1,974 confirmed T Tauri members of various sub-regions of the Orion OB1 association, and with ages spanning 4-10 Myr, using optical time-series from three different surveys. For 564 of the stars (~32% of the weak-lined T Tauri stars and ~13% of the classical T Tauri stars in our sample) we detect statistically significant periodic variations which we attribute to the stellar rotation periods, making this one of the largest samples of T Tauri star rotation periods yet published. We observe a clear change in the overall rotation period distributions over the age range 4-10 Myr, with the progressively older sub-populations exhibiting systematically faster rotation. This result is consistent with angular momentum evolution model predictions of an important qualitative change in the stellar rotation periods starting at ~5 Myr, an age range for which very few observational constraints were previously available. △ Less

Submitted 15 August, 2016; v1 submitted 13 May, 2016; originally announced May 2016.

Comments: Accepted for publication in the Astronomical Journal. 16 pages, full appendices will be available in the electronic journal

Mapping optically variable quasars towards the galactic plane

Authors: J. G. Fernandez-Trincado, T. Verdugo, C. Reylé, A. C Robin, J. A. de Diego, V. Motta, L. Vega, J. J. Downes, C. Mateu, A. K. Vivas, C. Briceño, C. Abad, K. Vieira, J. Hernández, A. Nuñez, E. Gatuzz

Abstract: We present preliminary results of the CIDA Equatorial Variability Survey (CEVS), looking for quasar (hereafter QSO) candidates near the Galactic plane. The CEVS contains photometric data from extended and adjacent regions of the Milky Way disk ($\sim$ 500 sq. deg.). In this work 2.5 square degrees with moderately high temporal sampling in the CEVS were analyzed. The selection of QSO candidates was… ▽ More We present preliminary results of the CIDA Equatorial Variability Survey (CEVS), looking for quasar (hereafter QSO) candidates near the Galactic plane. The CEVS contains photometric data from extended and adjacent regions of the Milky Way disk ($\sim$ 500 sq. deg.). In this work 2.5 square degrees with moderately high temporal sampling in the CEVS were analyzed. The selection of QSO candidates was based on the study of intrinsic optical photometric variability of 14,719 light curves. We studied samples defined by cuts in the variability index (Vindex $>$ 66.5), periodicity index (Q $>$ 2), and the distribution of these sources in the plane (AT , $γ$), using a slight modification of the first-order of the structure function for the temporal sampling of the survey. Finally, 288 sources were selected as QSO candidates. The results shown in this work are a first attempt to develop a robust method to detect QSO towards the Galactic plane in the era of massive surveys such as VISTA and Gaia. △ Less

Submitted 3 October, 2015; originally announced October 2015.

Comments: 7 pages, 2 figures, to appear in the proceedings of SF2A 2015

Detection of satellite remnants in the Galactic Halo with Gaia III. Detection limits for Ultra Faint Dwarf Galaxies

Authors: Teresa Antoja, Cecilia Mateu, Luis Aguilar, Francesca Figueras, Erika Antiche, Fabiola Hernandez-Perez, Anthony Brown, Octavio Valenzuela, Antonio Aparicio, Sebastian Hidalgo, Hector Velazquez

Abstract: 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 wi… ▽ More 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. △ Less

Submitted 15 July, 2015; originally announced July 2015.

Comments: Accepted for publication in MNRAS

The number fraction of discs around brown dwarfs in Orion OB1a and the 25 Orionis group

Authors: Juan José Downes, Carlos Román-Zúñiga, Javier Ballesteros-Paredes, Cecilia Mateu, César Briceño, Jesús Hernández, Monika G. Petr-Gotzens, Nuria Calvet, Lee Hartmann, Karina Mauco

Abstract: We present a study of 15 new brown dwarfs belonging to the $\sim7$ Myr old 25 Orionis group and Orion OB1a sub-association with spectral types between M6 and M9 and estimated masses between $\sim0.07$M$_\odot$ and $\sim0.01$ M$_\odot$. By comparing them through a Bayesian method with low mass stars ($0.8\lesssim$ M/M$_\odot\lesssim0.1$) from previous works in the 25 Orionis group, we found statist… ▽ More We present a study of 15 new brown dwarfs belonging to the $\sim7$ Myr old 25 Orionis group and Orion OB1a sub-association with spectral types between M6 and M9 and estimated masses between $\sim0.07$M$_\odot$ and $\sim0.01$ M$_\odot$. By comparing them through a Bayesian method with low mass stars ($0.8\lesssim$ M/M$_\odot\lesssim0.1$) from previous works in the 25 Orionis group, we found statistically significant differences in the number fraction of classical T Tauri stars, weak T Tauri stars, class II, evolved discs and purely photospheric emitters at both sides of the sub-stellar mass limit. Particularly we found a fraction of $3.9^{+2.4}_{-1.6}~\%$ low mass stars classified as CTTS and class II or evolved discs, against a fraction of $33.3^{+10.8}_{-9.8}~\%$ in the sub-stellar mass domain. Our results support the suggested scenario in which the dissipation of discs is less efficient for decreasing mass of the central object. △ Less

Submitted 20 April, 2015; originally announced April 2015.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 15 pages, 9 figures and 5 tables

On the recovery of galaxy properties from SED fitting solutions

Authors: Gladis Magris C., Juan Mateu P., Cecilia Mateu, Gustavo Bruzual A., Ivan Cabrera-Ziri, Alfredo Mejía-Narváez

Abstract: We explore the ability of four different inverse population synthesis codes to recover the physical properties of galaxies from their spectra by SED fitting. Three codes, DynBaS, TGASPEX, and GASPEX, have been implemented by the authors and are described in detail in the paper. STARLIGHT, the fourth code, is publicly available. DynBaS selects dynamically a different spectral basis to expand the sp… ▽ More We explore the ability of four different inverse population synthesis codes to recover the physical properties of galaxies from their spectra by SED fitting. Three codes, DynBaS, TGASPEX, and GASPEX, have been implemented by the authors and are described in detail in the paper. STARLIGHT, the fourth code, is publicly available. DynBaS selects dynamically a different spectral basis to expand the spectrum of each target galaxy; TGASPEX uses an unconstrained age basis, whereas GASPEX and STARLIGHT use for all fits a fixed spectral basis selected a priori by the code developers. Variable and unconstrained basis reflect the peculiarities of the fitted spectrum and allow for simple and robust solutions to the problem of extracting galaxy parameters from spectral fits. We assemble a Synthetic Spectral Atlas of Galaxies (SSAG), comprising 100,000 galaxy spectra corresponding to an equal number of star formation histories based on the recipe of Chen et al. (2012). We select a subset of 120 galaxies from SSAG with a colour distribution similar to that of local galaxies in the seventh data release (DR7) of the Sloan Digital Sky Survey (SDSS) and produce 30 random noise realisations for each of these spectra. For each spectrum we recover the mass, mean age, metallicity, internal dust extinction, and velocity dispersion characterizing the dominant stellar population in the problem galaxy. All methods produce almost perfect fits to the target spectrum, but the recovered physical parameters can differ significantly. Our tests provide a quantitative measure of the accuracy and precision with which these parameters are recovered by each method. From a statistical point of view all methods yield similar precisions, whereas DynBaS produces solutions with minimal systematic biases in the distributions of residuals for all of these parameters. △ Less

Submitted 7 January, 2015; v1 submitted 25 November, 2014; originally announced November 2014.

Comments: Minor typos fixed. Accepted for publication at PASP. Synthetic Spectral Atlas of Galaxies (SSAG) available at http://www.astro.ljmu.ac.uk/~asticabr/SSAG.html

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

Authors: J. G. Fernández-Trincado, A. K. Vivas, C. E. Mateu, R. Zinn, A. C. Robin, O. Valenzuela, E. Moreno, B. Pichardo

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

Submitted 11 November, 2014; originally announced November 2014.

Comments: 11 pages, 11 figures, Accepted

Report number: AA/2014/24899

Journal ref: A&A 574, A15 (2015)

The low mass star and sub-stellar populations of the 25 Orionis group

Authors: Juan José Downes, César Briceño, Cecilia Mateu, Jesús Hernández, Anna Katherina Vivas, Nuria Calvet, Lee Hartmann, Monika G. Petr-Gotzens, Lori Allen

Abstract: We present the results of a survey of the low mass star and brown dwarf population of the 25 Orionis group. Using optical photometry from the CIDA Deep Survey of Orion, near IR photometry from the Visible and Infrared Survey Telescope for Astronomy and low resolution spectroscopy obtained with Hectospec at the MMT, we selected 1246 photometric candidates to low mass stars and brown dwarfs with est… ▽ More We present the results of a survey of the low mass star and brown dwarf population of the 25 Orionis group. Using optical photometry from the CIDA Deep Survey of Orion, near IR photometry from the Visible and Infrared Survey Telescope for Astronomy and low resolution spectroscopy obtained with Hectospec at the MMT, we selected 1246 photometric candidates to low mass stars and brown dwarfs with estimated masses within $0.02 \lesssim M/M_\odot \lesssim 0.8$ and spectroscopically confirmed a sample of 77 low mass stars as new members of the cluster with a mean age of $\sim$7 Myr. We have obtained a system initial mass function of the group that can be well described by either a Kroupa power-law function with indices $α_3=-1.73\pm0.31$ and $α_2=0.68\pm0.41$ in the mass ranges $0.03\leq M/M_\odot\leq 0.08$ and $0.08\leq M/M_\odot\leq0.5$ respectively, or a Scalo log-normal function with coefficients $m_c=0.21^{+0.02}_{-0.02}$ and $σ=0.36\pm0.03$ in the mass range $0.03\leq M/M_\odot\leq0.8$. From the analysis of the spatial distribution of this numerous candidate sample, we have confirmed the East-West elongation of the 25 Orionis group observed in previous works, and rule out a possible southern extension of the group. We find that the spatial distributions of low mass stars and brown dwarfs in 25 Orionis are statistically indistinguishable. Finally, we found that the fraction of brown dwarfs showing IR excesses is higher than for low mass stars, supporting the scenario in which the evolution of circumstellar discs around the least massive objects could be more prolonged. △ Less

Submitted 31 July, 2014; originally announced August 2014.

Comments: Accepted for publication in MNRAS

Characterizing the Galactic warp with Gaia: I. The tilted ring model with a twist

Authors: H. Abedi, C. Mateu, L. A. Aguilar, F. Figueras, M. Romero-Gomez

Abstract: 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… ▽ More 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. △ Less

Submitted 23 May, 2014; originally announced May 2014.

Comments: 18 pages, 16 figures, accepted for publication in MNRAS

The QUEST RR Lyrae Survey: III. The Low Galactic Latitude Catalogue

Authors: Cecilia Mateu, A. Katherina Vivas, Juan José Downes, César Briceño, Robert Zinn, Gustavo Cruz-Diaz

Abstract: We present results for the QUEST RR Lyrae Survey at low galactic latitude, conducted entirely with observations obtained with the QUEST mosaic camera and the 1.0/1.5m Jürgen Stock Schmidt telescope at the National Observatory of Venezuela. The survey spans an area of 476 sq. deg on the sky, with multi-epoch observations in the V, R and I photometric bands for 6.5x10^6 stars in the galactic latitud… ▽ More We present results for the QUEST RR Lyrae Survey at low galactic latitude, conducted entirely with observations obtained with the QUEST mosaic camera and the 1.0/1.5m Jürgen Stock Schmidt telescope at the National Observatory of Venezuela. The survey spans an area of 476 sq. deg on the sky, with multi-epoch observations in the V, R and I photometric bands for 6.5x10^6 stars in the galactic latitude range 30<= b(deg) <=+25, in a direction close to the Galactic Anticenter 190<= l(deg)<= 230. The variability survey has a typical number of 30 observations per object in V and I and ~25 in R, with up to ~120-150 epochs in V and I and up to ~100 in R in the best sampled regions. The completeness magnitudes of the survey are V=R=18.5 mag, and I=18.0 mag. We identified 211 RR Lyrae stars, 160 bona fide stars of type ab and 51 candidates of type c, ours being the first deep RR Lyrae survey conducted at low galactic latitude.The completeness of the RR Lyrae survey was estimated in >95 per cent and ~85 per cent for RRab and RRc stars respectively. Photometric metallicities were computed based on the light curves and individual extinctions calculated from minimum light colours for each RRab star. Distances were obtained with typical errors ~7 per cent. The RR Lyrae survey simultaneously spans a large range of heliocentric distances 0.5<= R_hel(kpc)<=40 and heights above the plane -15<=z(kpc)<=+20, with well known completeness across the survey area, making it an ideal set for studying the struc ture of the Galactic thick disk. △ Less

Submitted 22 August, 2012; originally announced August 2012.

Comments: 24 pages, 17 figures, 8 tables. Accepted for publication in MNRAS

Detection of Satellite Remnants in the Galactic Halo with Gaia - II. A modified Great Circle Cell Method

Authors: Cecilia Mateu, Gustavo Bruzual, Luis Aguilar, Anthony Brown, Octavio Valenzuela, Leticia Carigi, Fabiola Hernandez, Hector Velazquez

Abstract: 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 com… ▽ More 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. △ Less

Submitted 1 March, 2011; originally announced March 2011.

Comments: 13 pages, 12 figures. Accepted for publication in MNRAS

No Excess of RR Lyrae Stars in the Canis Major Overdensity

Authors: Cecilia Mateu, A. Katherina Vivas, Robert Zinn, Lissa Miller, Carlos Abad

Abstract: Our multi-epoch survey of ~20 sq. deg. of the Canis Major overdensity has detected only 10 RR Lyrae stars (RRLS). We show that this number is consistent with the number expected from the Galactic halo and thick disk populations alone, leaving no excess that can be attributed to the dwarf spheroidal (dSph) galaxy that some authors have proposed as the origin of the CMa overdensity. If this galaxy… ▽ More Our multi-epoch survey of ~20 sq. deg. of the Canis Major overdensity has detected only 10 RR Lyrae stars (RRLS). We show that this number is consistent with the number expected from the Galactic halo and thick disk populations alone, leaving no excess that can be attributed to the dwarf spheroidal (dSph) galaxy that some authors have proposed as the origin of the CMa overdensity. If this galaxy resembles the dSph satellites of the Milky Way and of M31 and has the putative Mv~-14.5, our survey should have detected several tens of RRLS. Even if Mv<-12, the expected excess is >10, which is not observed. Either the old stellar population of this galaxy has unique properties or, as others have argued before, the CMa overdensity is produced by the thin and thick disk and spiral arm populations of the Milky Way and not by a collision with a dSph satellite galaxy. △ Less

Submitted 2 March, 2009; originally announced March 2009.

Comments: 35 pages, 9 figures, 5 tables. Accepted for publication at the Astronomical Journal

Spectroscopy of Bright QUEST RR Lyrae Stars: Velocity Substructures toward Virgo

Authors: A. Katherina Vivas, Yara Jaffe, Robert Zinn, Rebeccah Winnick, Sonia Duffau, Cecilia Mateu

Abstract: Using a sample of 43 bright (V<16.1, distance <13 kpc) RR Lyrae stars (RRLS) from the QUEST survey with spectroscopic radial velocities and metallicities, we find that several separate halo substructures contribute to the Virgo overdensity (VOD). While there is little evidence for halo substructure in the spatial distribution of these stars, their distribution in radial velocity reveals two movi… ▽ More Using a sample of 43 bright (V<16.1, distance <13 kpc) RR Lyrae stars (RRLS) from the QUEST survey with spectroscopic radial velocities and metallicities, we find that several separate halo substructures contribute to the Virgo overdensity (VOD). While there is little evidence for halo substructure in the spatial distribution of these stars, their distribution in radial velocity reveals two moving groups. These results are reinforced when the sample is combined with a sample of blue horizontal branch stars that were identified in the SDSS, and the combined sample provides evidence for one additional moving group. These groups correspond to peaks in the radial velocity distribution of a sample of F type main-sequence stars that was recently observed in the same directon by SEGUE, although in one case the RRLS and F star groups may not lie at the same distance. One of the new substructures has a very narrow range in metallicity, which is more consistent with it being the debris from a destroyed globular cluster than from a dwarf galaxy. A small concentration of stars have radial velocities that are similar to the Virgo Stellar Stream (VSS) that was identified previously in a fainter sample of RRLS. Our results suggest that this feature extends to distances as short as ~12 kpc from its previous detection at ~19 kpc. None of the new groups and only one star in the sample have velocities that are consistent with membership in the leading tidal stream from the Sagittarius Dwarf Spheroidal Galaxy, which some authors have suggested is the origin of the VOD. △ Less

Submitted 10 July, 2008; originally announced July 2008.

Comments: Accepted for publication in the AJ

Journal ref: Astron.J.136:1645-1655,2008

Searching for RR Lyrae stars in the Canis Major Overdensity

Authors: Cecilia Mateu, A. Katherina Vivas, Robert Zinn, Lissa Miller, .

Abstract: The Canis Major overdensity (CMa) was initially proposed to be the remnant of a tidally disrupting dSph galaxy. Since its nature is still subject of debate, the goal of the present work was to conduct a large-scale RR Lyrae survey in CMa, in order to see if there is an overdensity of these stars. The survey spans a total area of ~34 sq. deg. with observations in V and R filters, made with the 1.… ▽ More The Canis Major overdensity (CMa) was initially proposed to be the remnant of a tidally disrupting dSph galaxy. Since its nature is still subject of debate, the goal of the present work was to conduct a large-scale RR Lyrae survey in CMa, in order to see if there is an overdensity of these stars. The survey spans a total area of ~34 sq. deg. with observations in V and R filters, made with the 1.0m Jurgen Stock Schmidt telescope at the National Astronomical Observatory of Venezuela. Current results in a subregion, including spectroscopic observations, show that the small number of RR Lyrae stars found can be accounted for by the halo and thick disk components of our Galaxy. △ Less

Submitted 6 February, 2007; v1 submitted 5 February, 2007; originally announced February 2007.

Comments: 2 pages, 1 figure, to appear in the proceedings of IAU Symposium No241 "Stellar Populations as Building Blocks of Galaxies"

RR Lyrae search and stellar populations study in Canis Major: Preliminary Results

Authors: Cecilia E. Mateu, A. Katherina Vivas, Robert Zinn, Lissa Miller

Abstract: We present preliminary results of a RR Lyrae star search and stellar populations study performed in the Canis Major overdensity, spanning an area of 8.35sq deg. The observations were made in R and V bands, with the QUEST camera installed in the 1m Jurgen Stock Telescope, at the Venezuela National Observatory. The resulting Hess diagram shows a possible, but weak, red giant branch and no obvious… ▽ More We present preliminary results of a RR Lyrae star search and stellar populations study performed in the Canis Major overdensity, spanning an area of 8.35sq deg. The observations were made in R and V bands, with the QUEST camera installed in the 1m Jurgen Stock Telescope, at the Venezuela National Observatory. The resulting Hess diagram shows a possible, but weak, red giant branch and no obvious horizontal branch, red clump or main sequence turnoff. After a multi--epoch photometric search, 6 RR Lyrae stars were confirmed with further observations obtained at the 1.0 and 1.3m telescopes of the SMARTS consortium at CTIO. Of these confirmed RR Lyrae stars, 5 have heliocentric distances between 5 and 7 kpc. Confirmation of their physsical association with the Canis Major system awaits for a study of their radial velocities. △ Less

Submitted 14 April, 2005; v1 submitted 14 April, 2005; originally announced April 2005.

Comments: 4 pages, 10 figures. To appear in the proceedings of the International Workshop Astronomia Dinamica en Latino-America (ADeLA 2004), Merida, Venezuela, Nov 22-26, 2004, RMAA Conf Ser., C. Abad, A. Bongiovanni & Y. Guillen (eds)