We present the first SkyMapper stellar population analysis of the Large Magellanic Cloud (hereafter LMC),including the identification of 3578 candidate Carbon Stars through their extremely red $g-r$ colours. Coupled with Gaia astrometry, we analyse the distribution and kinematics of this Carbon Star population, finding the LMC to be centred at $(R.A., Dec.) = (80.90^{\circ}\pm{0.29}, -68.74^{\circ}\pm{0.12})$, with a bulk proper motion of $(\mu_{\alpha},\mu_{\delta}) = (1.878\pm0.007,0.293\pm0.018) \mathrm{mas\ yr^{-1}}$ and a disk inclination of $i = 25.6^{\circ}\pm1.1$ at position angle $\theta = 135.6^{\circ}\pm 3.3^{\circ}$. We complement this study with the identification and analysis of additional stellar populations, finding that the dynamical centre for Red Giant Branch (RGB) stars is similar to that seen for the Carbon Stars, whereas for young stars the dynamical centre is significantly offset from the older populations. This potentially indicates that the young stars were formed as a consequence of a strong tidal interaction, probably with the Small Magellanic Cloud (SMC). In terms of internal dynamics, the tangential velocity profile increases linearly within $\sim3\ \mathrm{kpc}$, after which it maintains an approximately constant value of $V_{rot} = 83.6\pm 1.7 \mathrm{km\ s^{-1}}$ until $\sim7 \mathrm{kpc}$. With an asymmetric drift correction, we estimate the mass within $7 \mathrm{kpc}$ to be $M_{\rm LMC}(<7\mathrm{kpc}) = (2.5\pm0.1)\times10^{10}\ \mathrm{M}_{\odot}$ and within the tidal radius ($\sim 30\ \mathrm{kpc}$) to be $M_{\rm LMC}(<30\mathrm{kpc}) = (1.06 \pm 0.32)\times10^{11}\ \mathrm{M}_{\odot}$, consistent with other recent measurements.

There is compelling evidence for a highly energetic Seyfert explosion (10^{56-57} erg) that occurred in the Galactic Centre a few million years ago. The clearest indications are the x-ray/gamma-ray "10 kpc bubbles" identified by the Rosat and Fermi satellites. In an earlier paper, we suggested another manifestation of this nuclear activity, i.e. elevated H-alpha emission along a section of the Magellanic Stream due to a burst (or flare) of ionizing radiation from Sgr A*. We now provide further evidence for a powerful flare event: UV absorption line ratios (in particular CIV/CII, SiIV/SiII) observed by the Hubble Space Telescope reveal that some Stream clouds towards both galactic poles are highly ionized by a source capable of producing ionization energies up to at least 50 eV. We show how these are clouds caught in a beam of bipolar, radiative "ionization cones" from a Seyfert nucleus associated with Sgr A*. In our model, the biconic axis is tilted by about 15 deg from the South Galactic Pole with an opening angle of roughly 60 deg. For the Stream at such large Galactic distances (D > 75 kpc), nuclear activity is a plausible explanation for all of the observed signatures: elevated H-alpha emission and H ionization fraction (X_e > 0.5), enhanced CIV/CII and SiIV/SiII ratios, and high CIV and SiIV column densities. Wind-driven "shock cones" are ruled out because the Fermi bubbles lose their momentum and energy to the Galactic corona long before reaching the Stream. The nuclear flare event must have had a radiative UV luminosity close to the Eddington limit (f_E ~ 0.1-1). Our time-dependent Seyfert flare models adequately explain the observations and indicate the Seyfert flare event took place T_o = 3.5 +/- 1 Myr ago. The timing estimates are consistent with the mechanical timescales needed to explain the x-ray/gamma-ray bubbles in leptonic jet/wind models (2-8 Myr).

We present an innovative approach to the methodology of dynamical modelling, allowing practical reconstruction of the underlying dark matter mass without assuming both the density and anisotropy functions. With this, the mass-anisotropy degeneracy is reduced to simple model inference, incorporating the uncertainties inherent with observational data, statistically circumventing the mass-anisotropy degeneracy in spherical collisionless systems. We also tackle the inadequacy that the Jeans method of moments has on small datasets, with the aid of Generative Adversarial Networks: we leverage the power of artificial intelligence to reconstruct non-parametrically the projected line-of-sight velocity distribution. We show with realistic numerical simulations of dwarf spheroidal galaxies that we can distinguish between competing dark matter distributions and recover the anisotropy and mass profile of the system.

The Monoceros Ring (MRi) structure is an apparent stellar overdensity that has been postulated to entirely encircle the Galactic plane and has been variously described as being due to line-of-sight effects of the Galactic warp and flare or of extragalactic origin (via accretion). Despite being intensely scrutinised in the literature for more than a decade, no studies to-date have been able to definitively uncover its origins. Here we use $N$-body simulations and a genetic algorithm to explore the parameter space for the initial position, orbital parameters and, for the first time, the final location of a satellite progenitor. We fit our models to the latest Pan-STARRS data to determine whether an accretion scenario is capable of producing an in- Plane ring-like structure matching the known parameters of the MRi. Our simulations produce streams that closely match the location, proper motion and kinematics of the MRi structure. However, we are not able to reproduce the mass estimates from earlier studies based on Pan-STARRS data. Furthermore, in contrast with earlier studies our best-fit models are those for progenitors on retrograde orbits. If the MRi was produced by satellite accretion, we find that its progenitor has an initial mass upper limit of ~$10^{10}$M$_\odot$ and the remnant is likely located behind the Galactic bulge, making it diffcult to locate observationally. While our models produce realistic MRi-like structures we cannot definitively conclude that the MRi was produced by the accretion of a satellite galaxy.

We present results from deep and very spatially extended CTIO/DECam $g$ and $r$ photometry (reaching out to $\sim$ 2 mag below the oldest MSTO and covering $\sim$ 20 deg$^2$) around the Sextans dSph. We use this data-set to study the structural properties of Sextans overall stellar population and its different stellar evolutionary phases, as well as to search for signs of tidal disturbance from the MW, which would indicate departure from dynamical equilibrium. We perform the most accurate structural analysis to-date of Sextans' stellar components by applying Bayesian MCMC methods to the individual stars' positions. Surface density maps are built by decontaminating the sample through a matched filter analysis of the CMD, and then analysed for departures from axisymmetry. Sextans is found to be considerably less spatially extended than early studies suggested. No significant distortions or tidal disturbances are found down to a surface brightness of $\sim$ 31.8 mag/arcsec$^{-2}$ in V-band. We identify an overdensity in the central regions that may correspond to previously reported kinematic substructure(s). In agreement with previous findings, old $\&$ metal-poor stars such as BHB stars cover a much larger area than stars in other evolutionary phases, and bright BSs are less spatially extended than faint ones. However, the different spatial distribution of bright and faint BSs appears consistent with the general age/metallicity gradients found in Sextans' stellar component. This is compatible with BSs having formed by evolution of binaries and not necessarily due to the presence of a central disrupted globular cluster, as suggested in the literature. We provide structural parameters for the various populations analyzed and make publicly available the photometric catalogue of point-sources as well as a catalogue of literature spectroscopic measurements with updated membership probabilities.

The first and second moments of stellar velocities encode important information about the formation history of the Galactic halo. However, due to the lack of tangential motion and inaccurate distances of the halo stars, the velocity moments in the Galactic halo have largely remained "known unknowns". Fortunately, our off-centric position within the Galaxy allows us to estimate these moments in the galacto-centric frame using the observed radial velocities of the stars alone. We use these velocities coupled with the Hierarchical Bayesian scheme, which allows easy marginalisation over the missing data (the proper-motion, and uncertainty-free distance and line-of-sight velocity), to measure the velocity dispersions, orbital anisotropy ($\beta$) and streaming motion ($v_{\rm rot}$) of the halo main-sequence turn-off (MSTO) and K-giant (KG) stars in the inner stellar halo (r $\lesssim 15$ kpc). We study the metallicity bias in kinematics of the halo stars and observe that the comparatively metal-rich ([Fe/H]$>-1.4$) and the metal-poor ([Fe/H]$\leq - 1.4$) MSTO samples show a clear systematic difference in $v_{\rm rot} \sim 20-40$ km s$^{-1}$, depending on how restrictive the spatial cuts to cull the disk contamination are. The bias is also detected in KG samples but with less certainty. Both MSTO and KG populations suggest that the inner stellar halo of the Galaxy is radially biased i.e. $\sigma_r>\sigma_\theta$ or $\sigma_\phi$ and $\beta \simeq 0.5$. The apparent metallicity contrariety in the rotation velocity among the halo sub-populations supports the co-existence of multiple populations in the galactic halo that may have formed through distinct formation scenarios, i.e. in-situ versus accretion.

Dwarf galaxies, among the most dark matter dominated structures of our universe, are excellent test-beds for dark matter theories. Unfortunately, mass modelling of these systems suffers from the well documented mass-velocity anisotropy degeneracy. For the case of spherically symmetric systems, we describe a method for non-parametric modelling of the radial and tangential velocity moments. The method is a numerical velocity anisotropy "inversion", with parametric mass models, where the radial velocity dispersion profile, $\sigma_{\mathrm{rr}}^2$ is modeled as a B-spline, and the optimization is a three step process that consists of: (i) an Evolutionary modelling to determine the mass model form and the best B-spline basis to represent $\sigma_{\mathrm{rr}}^2$; (ii) an optimization of the smoothing parameters; (iii) a Markov chain Monte Carlo analysis to determine the physical parameters. The mass-anisotropy degeneracy is reduced into mass model inference, irrespective of kinematics. We test our method using synthetic data. Our algorithm constructs the best kinematic profile and discriminates between competing dark matter models. We apply our method to the Fornax dwarf spheroidal galaxy. Using a King brightness profile and testing various dark matter mass models, our model inference favours a simple mass-follows-light system. We find that the anisotropy profile of Fornax is tangential ($\beta(r) < 0$) and we estimate a total mass of $M_{\text{tot}} = 1.613 ^{+0.050}_{-0.075} \times 10^8 \, \text{M}_{\odot}$, and a mass-to-light ratio of $\Upsilon_V = 8.93 ^{+0.32}_{-0.47} \, (\text{M}_{\odot}/\text{L}_{\odot})$. The algorithm we present is a robust and computationally inexpensive method for non-parametric modelling of spherical clusters independent of the mass-anisotropy degeneracy.

The recently discovered vast thin plane of dwarf satellites orbiting the Andromeda Galaxy (M31) adds to the mystery of the small scale distribution of the Local Group's galaxy population. Such well defined planar structures are apparently rare occurrences in cold dark matter cosmological simulations, and we lack a coherent explanation of their formation and existence. In this paper, we explore the long-term survivability of thin planes of dwarfs in galactic halos, focusing, in particular, on systems mimicking the observed Andromeda distribution. The key results show that, in general, planes of dwarf galaxies are fragile, sensitive to the shape of the dark matter halo and other perturbing effects. In fact, long lived planes of satellites only exist in polar orbits in spherical dark matter halos, presenting a challenge to the observed Andromeda plane which is significantly tilted with respect to the optical disk. Our conclusion is that, in standard cosmological models, planes of satellites are generally short lived, and hence we must be located at a relatively special time in the evolution of the Andromeda Plane, lucky enough to see its coherent pattern.

Astronomy is undergoing through a methodological revolution triggered by an unprecedented wealth of complex and accurate data. DAMEWARE (DAta Mining & Exploration Web Application and REsource) is a general purpose, Web-based, Virtual Observatory compliant, distributed data mining framework specialized in massive data sets exploration with machine learning methods. We present the DAMEWARE (DAta Mining & Exploration Web Application REsource) which allows the scientific community to perform data mining and exploratory experiments on massive data sets, by using a simple web browser. DAMEWARE offers several tools which can be seen as working environments where to choose data analysis functionalities such as clustering, classification, regression, feature extraction etc., together with models and algorithms.

NGC147, NGC185 and CassiopeiaII (CassII) have similar positions in the sky, distances and measured line of sight velocities. This proximity in phase space suggests that these three satellites of M31 form a subgroup within the Local Group. Nevertheless, the differences in their star formation history and interstellar medium, and the recent discovery of a stellar stream in NGC~147, combined with the lack of tidal features in the other two satellites, are all indications of complex and diverse interactions between M31 and these three satellites. We use a genetic algorithm to explore the different orbits that these satellites can have and select six sets of orbits that could best explain the observational features of the NGC147, NGC185 and CassII satellites. The parameters of these orbits are then used as a starting point for N-body simulations. We present models for which NGC147, NGC185 and CassII are a bound group for a total time of at least one Gyr but still undergo different interactions with M31 and as a result NGC147 has a clear stellar stream whereas the other two satellites have no significant tidal features. This result shows that it is possible to find solutions that reproduce the contrasting properties of the satellites and for which NGC147-NGC185-CassII have been gravitationally bound.

We present deep optical photometry from the DECam imager on the 4m Blanco telescope of over 12 deg$^2$ around the Carina dwarf spheroidal, with complete coverage out to 1 degree and partial coverage extending out to 2.6 degrees. Using a Poisson-based matched filter analysis to identify stars from each of the three main stellar populations, old, intermediate, and young, we confirm the previously identified radial age gradient, distance, tidal radius, stellar radial profiles, relative stellar population sizes, ellipticity, and position angle. We find an angular offset between the three main elliptical populations of Carina, and find only tentative evidence for tidal debris, suggesting that past tidal interactions could not have significantly influenced the Carina dwarf. We detect stars in the vicinity of, but distinct to, the Carina dwarf, and measure their distance to be 46$\pm$2 kpc. We determine this population to be part of the halo of the Large Magellanic Cloud at an angular radius of over 20 degrees. Due to overlap in colour-magnitude space with Magellanic stars, previously detected tidal features in the old population of Carina are likely weaker than previously thought.

The history of the Magellanic Clouds is investigated using N-body hydrodynamic simulations where the initial conditions are set by a genetic algorithm. This technique allows us to identify possible orbits for the Magellanic Clouds around the Milky Way, by directly comparing the simulations with observational constraints. We explore the parameter space of the interaction between the Magellanic Clouds and the Milky Way, considering as free parameters the proper motions of the Magellanic Clouds, the virial mass and the concentration parameter ($c$) of the Galactic dark matter halo. The best orbital scenarios presented here are considered with two different sets of parameters for the Milky Way disc and bulge components. The total circular velocity at the Sun's position ($\rm{R_{\odot}}=8.5\,\rm {kpc}$) is directly calculated from the rotation curve of the corresponding Galactic mass model. Our analysis suggests that the Magellanic Clouds have orbited inside the virial radius of the Milky Way for at least $3\,\rm{Gyr}$, even for low mass haloes. However, this is possible only with high values for the concentration parameter ($c\ge20$). In both orbital models presented here, the mutual interaction between the Magellanic Clouds is able to reproduce the observed features of the Magellanic System.

Astronomy is undergoing through a methodological revolution triggered by an unprecedented wealth of complex and accurate data. The new panchromatic, synoptic sky surveys require advanced tools for discovering patterns and trends hidden behind data which are both complex and of high dimensionality. We present DAMEWARE (DAta Mining & Exploration Web Application REsource): a general purpose, web-based, distributed data mining environment developed for the exploration of large datasets, and finely tuned for astronomical applications. By means of graphical user interfaces, it allows the user to perform classification, regression or clustering tasks with machine learning methods. Salient features of DAMEWARE include its capability to work on large datasets with minimal human intervention, and to deal with a wide variety of real problems such as the classification of globular clusters in the galaxy NGC1399, the evaluation of photometric redshifts and, finally, the identification of candidate Active Galactic Nuclei in multiband photometric surveys. In all these applications, DAMEWARE allowed to achieve better results than those attained with more traditional methods. With the aim of providing potential users with all needed information, in this paper we briefly describe the technological background of DAMEWARE, give a short introduction to some relevant aspects of data mining, followed by a summary of some science cases and, finally, we provide a detailed description of a template use case.