Using Hubble Space Telescope ACS/WFC data we present the photometry and spatial distribution of resolved stellar populations of four fields within the extended ultraviolet disk (XUV disk) of M83. These observations show a clumpy distribution of main-sequence stars and a mostly smooth distribution of red giant branch stars. We constrain the upper-end of the initial mass function (IMF) in the outer disk using the detected population of main-sequence stars and an assumed constant star formation rate (SFR) over the last 300 Myr. By comparing the observed main-sequence luminosity function to simulations, we determine the best-fitting IMF to have a power law slope $\alpha=-2.35 \pm 0.3$ and an upper-mass limit $\rm M_{u}=25_{-3}^{+17} \, M_\odot$. This IMF is consistent with the observed H$\alpha$ emission, which we use to provide additional constraints on the IMF. We explore the influence of deviations from the constant SFR assumption, finding that our IMF conclusions are robust against all but strong recent variations in SFR, but these are excluded by causality arguments. These results, along with our similar studies of other nearby galaxies, indicate that some XUV disks are deficient in high-mass stars compared to a Kroupa IMF. There are over one hundred galaxies within 5 Mpc, many already observed with HST, thus allowing a more comprehensive investigation of the IMF, and how it varies, using the techniques developed here.

Large galaxies grow through the accumulation of dwarf galaxies. In principle it is possible to trace this growth history using the properties of a galaxy's stellar halo. Previous investigations of the galaxy M31 (Andromeda) have shown that outside a radius of 25 kpc the population of halo globular clusters is rotating in alignment with the stellar disk, as are more centrally located clusters. The M31 halo also contains coherent stellar substructures, along with a smoothly distributed stellar component. Many of the globular clusters outside 25 kpc are associated with the most prominent substructures, while others are part of the smooth halo. Here we report a new analysis of the kinematics of these globular clusters. We find that the two distinct populations are rotating with perpendicular orientations. The rotation axis for the population associated with the smooth halo is aligned with the rotation axis for the plane of dwarf galaxies that encircles M31. We interpret these separate cluster populations as arising from two major accretion epochs, likely separated by billions of years. Stellar substructures from the first epoch are gone, but those from the more recent second epoch still remain.

We use Hyper Suprime-Cam on the Subaru Telescope to investigate the structural and photometric properties of early-type dwarf galaxies and young stellar systems at the center of the M81 Group. We have mapped resolved stars to $\sim2$ magnitudes below the tip of the red giant branch over almost 6.5 square degrees, corresponding to a projected area of $160\times160 \rm{kpc}$ at the distance of M81. The resulting stellar catalogue enables a homogeneous analysis of the member galaxies with unprecedented sensitivity to low surface brightness emission. The radial profiles of the dwarf galaxies are well-described by Sersic and King profiles, and show no obvious signatures of tidal disruption. The measured radii for most of these systems are larger than the existing literature values and we find the total luminosity of IKN ($\rm{M_{V,0}}=-14.29$) to be almost 3 magnitudes brighter than previously-thought. We identify new dwarf satellite candidates, d1006+69 and d1009+68, which we estimate to lie at a distance of $4.3\pm0.2$ Mpc and $3.5\pm0.5$ Mpc. With $\rm{M_{V,0}}=-8.91\pm0.40$ and $\rm{[M/H]}=-1.83\pm0.28$, d1006+69 is one of the faintest and most metal-poor dwarf satellites currently-known in the M81 Group. The luminosity functions of young stellar systems in the outlying tidal HI debris imply continuous star formation in the recent past and the existence of populations as young as 30 Myr old. We find no evidence for old RGB stars coincident with the young MS/cHeB stars which define these objects, supporting the idea that they are genuinely new stellar systems resulting from triggered star formation in gaseous tidal debris.

We present the star formation histories (SFHs) of 20 faint M31 satellites ($-12 \lesssim M_V \lesssim -6$) that were measured by modeling sub-horizontal branch (HB) depth color-magnitude diagrams constructed from Hubble Space Telescope (HST) imaging. Reinforcing previous results, we find that virtually all galaxies quenched between 3 and 9 Gyr ago, independent of luminosity, with a notable concentration $3-6$ Gyr ago. This is in contrast to the Milky Way (MW) satellites, which are generally either faint with ancient quenching times or luminous with recent ($<3$ Gyr) quenching times. We suggest that systematic differences in the quenching times of M31 and MW satellites may be a reflection of the varying accretion histories of M31 and the MW. This result implies that the formation histories of low-mass satellites may not be broadly representative of low-mass galaxies in general. Among the M31 satellite population we identify two distinct groups based on their SFHs: one with exponentially declining SFHs ($\tau \sim 2$ Gyr) and one with rising SFHs with abrupt quenching. We speculate how these two groups could be related to scenarios for a recent major merger involving M31. The Cycle 27 HST Treasury survey of M31 satellites will provide well-constrained ancient SFHs to go along with the quenching times we measure here. The discovery and characterization of M31 satellites with $M_V \gtrsim -6$ would help quantify the relative contributions of reionization and environment to quenching of the lowest-mass satellites.

We present new horizontal branch (HB) distance measurements to 17 of the faintest known M31 satellites ($-6 \lesssim M_{V} \lesssim -13$) based on deep Hubble Space Telescope (HST) imaging. The color-magnitude diagrams extend $\sim$1-2 magnitudes below the HB, which provides for well-defined HBs, even for faint galaxies in which the tip of the red giant branch (TRGB) is sparsely populated. We determine distances across the sample to an average precision of 4% ($\sim 30$~kpc at $800$~kpc). We find that the majority of these galaxies are in good agreement, though slightly farther (0.1-0.2 mag) when compared to recent ground-based TRGB distances. Two galaxies (And~IX and And~XVII) have discrepant HST and ground-based distances by $\sim 0.3$ mag ($\sim 150$~kpc), which may be due to contamination from Milky Way foreground stars and/or M31 halo stars in sparsely populated TRGB regions. We use the new distances to update the luminosities and structural parameters for these 17 M31 satellites. The new distances do not substantially change the spatial configuration of the M31 satellite system. We comment on future prospects for precise and accurate HB distances for faint galaxies in the Local Group and beyond.

We present the first stellar spectroscopy in the low-luminosity (M_V ~-9.3 mag), dwarf galaxy Leo P. Its significantly low oxygen abundance (3% solar) and relative proximity (~1.6 Mpc) make it a unique galaxy to investigate the properties of massive stars with near-primordial compositions akin to those in the early Universe. From our VLT-MUSE spectroscopy we find the first direct evidence for an O-type star in the prominent HII region, providing an important test case to investigate the potential environmental dependence of the upper end of the initial mass function in the dwarf galaxy regime. We classify 14 further sources as massive stars (and 17 more as candidate massive stars), most likely B-type objects. From comparisons with published evolutionary models we argue that the absolute visual magnitudes of massive stars in very metal-poor systems such as Leo P and I Zw 18 may be fainter by ~0.5 mag compared to Galactic stars. We also present spectroscopy of two carbon stars identified previously as candidate asymptotic-giant-branch stars. Two of three further candidate asymptotic-giant-branch stars display CaII absorption, confirming them as cool, evolved stars; we also recover CaII absorption in the stacked data of the next brightest 16 stars in the upper red giant branch. These discoveries will provide targets for future observations to investigate the physical properties of these objects and to calibrate evolutionary models of luminous stars at such low metallicity. The MUSE data also reveal two 100pc-scale ring structures in Halpha emission, with the HII region located on the northern edge of the southern ring. Lastly, we report serendipitous observations of 20 galaxies, with redshifts ranging from z=0.39, to a close pair of star-forming galaxies at z=2.5.

We utilise the final catalogue from the Pan-Andromeda Archaeological Survey to investigate the links between the globular cluster system and field halo in M31 at projected radii $R_p=25-150$ kpc. In this region the cluster radial density profile exhibits a power-law decline with index $\Gamma=-2.37\pm0.17$, matching that for the stellar halo component with [Fe/H] $<-1.1$. Spatial density maps reveal a striking correspondence between the most luminous substructures in the metal-poor field halo and the positions of many globular clusters. By comparing the density of metal-poor halo stars local to each cluster with the azimuthal distribution at commensurate radius, we reject the possibility of no correlation between clusters and field overdensities with high confidence. We use our stellar density measurements and previous kinematic data to demonstrate that $\approx35-60\%$ of clusters exhibit properties consistent with having been accreted into the outskirts of M31 at late times with their parent dwarfs. Conversely, at least $\sim40\%$ of remote clusters show no evidence for a link with halo substructure. The radial density profile for this subgroup is featureless and closely mirrors that observed for the apparently smooth component of the metal-poor stellar halo. We speculate that these clusters are associated with the smooth halo; if so, their properties appear consistent with a scenario where the smooth halo was built up at early times via the destruction of primitive satellites. In this picture the entire M31 globular cluster system outside $R_p=25$ kpc comprises objects accumulated from external galaxies over a Hubble time of growth.

Faint tidal features around galaxies record their merger and interaction histories over cosmic time. Due to their low surface brightnesses and complex morphologies, existing automated methods struggle to detect such features and most work to date has heavily relied on visual inspection. This presents a major obstacle to quantitative study of tidal debris features in large statistical samples, and hence the ability to be able to use these features to advance understanding of the galaxy population as a whole. This paper uses convolutional neural networks (CNNs) with dropout and augmentation to identify galaxies in the CFHTLS-Wide Survey that have faint tidal features. Evaluating the performance of the CNNs against previously-published expert visual classifications, we find that our method achieves high (76%) completeness and low (20%) contamination, and also performs considerably better than other automated methods recently applied in the literature. We argue that CNNs offer a promising approach to effective automatic identification of low surface brightness tidal debris features in and around galaxies. When applied to forthcoming deep wide-field imaging surveys (e.g. LSST, Euclid), CNNs have the potential to provide a several order-of-magnitude increase in the sample size of morphologically-perturbed galaxies and thereby facilitate a much-anticipated revolution in terms of quantitative low surface brightness science.

The Pan-Andromeda Archaeological Survey is a survey of $>400$ square degrees centered on the Andromeda (M31) and Triangulum (M33) galaxies that has provided the most extensive panorama of a $L_\star$ galaxy group to large projected galactocentric radii. Here, we collate and summarise the current status of our knowledge of the substructures in the stellar halo of M31, and discuss connections between these features. We estimate that the 13 most distinctive substructures were produced by at least 5 different accretion events, all in the last 3 or 4 Gyrs. We suggest that a few of the substructures furthest from M31 may be shells from a single accretion event. We calculate the luminosities of some prominent substructures for which previous estimates were not available, and we estimate the stellar mass budget of the outer halo of M31. We revisit the problem of quantifying the properties of a highly structured dataset; specifically, we use the OPTICS clustering algorithm to quantify the hierarchical structure of M31's stellar halo, and identify three new faint structures. M31's halo, in projection, appears to be dominated by two `mega-structures', that can be considered as the two most significant branches of a merger tree produced by breaking M31's stellar halo into smaller and smaller structures based on the stellar spatial clustering. We conclude that OPTICS is a powerful algorithm that could be used in any astronomical application involving the hierarchical clustering of points. The publication of this article coincides with the public release of all PAndAS data products.

We present a 3x3 degrees, 105-pointing, high-resolution neutral hydrogen (HI) mosaic of the M81 galaxy triplet (including the galaxies M81, M82 and NGC 3077, as well as dwarf galaxy NGC 2976) obtained with the Very Large Array (VLA) C and D arrays. This uniformly covers the entire area and velocity range of the triplet with a resolution of ~20'' or ~420 pc. The data reveal many small-scale anomalous velocity features highlighting the complexity of the interacting M81 triplet. We compare our data with Green Bank Telescope (GBT) observations of the same area. This provides evidence for a substantial reservoir of low-column density gas in the northern part of the triplet, probably associated with M82. Such a reservoir is not found in the southern part. We report a number of kpc-sized low-mass HI clouds with HI masses of a few times 10^6 Msun. Their dynamical masses are much larger than their baryonic masses, which could indicate the presence of dark matter if the clouds are rotationally supported. However, due to their spatial and kinematical association with HI tidal features, it is more likely that the velocity widths indicate tidal effects or streaming motions. We do not find any clouds not associated with tidal features down to an HI mass limit of a few times 10^4 Msun. We compare the HI column densities with resolved stellar density maps and find a star formation threshold around 3-6 10^20 cm-2$. We find that extreme velocity dispersions can be explained by a superposition of multiple components along the line of sight near M81 as well as winds or outflows around M82. The velocity dispersions found are high enough that these processes could explain the linewidths of Damped-Lyman-alpha absorbers observed at high redshift.

We discuss some of the key open questions regarding the formation and evolution of globular clusters (GCs) during galaxy formation and assembly within a cosmological framework. The current state-of-the-art for both observations and simulations is described, and we briefly mention directions for future research. The oldest GCs have ages $\ge$ 12.5 Gyr and formed around the time of reionisation. Resolved colour-magnitude diagrams of Milky Way GCs and direct imaging of lensed proto-GCs at z $\sim$ 6 with JWST promise further insight. Globular clusters are known to host multiple populations of stars with variations in their chemical abundances. Recently, such multiple populations have been detected in $\sim$2 Gyr old compact, massive star clusters. This suggests a common, single pathway for the formation of GCs at high and low redshift. The shape of the initial mass function for GCs remains unknown, however for massive galaxies a power-law mass function is favoured. Significant progress has been made recently modelling GC formation in the context of galaxy formation, with success in reproducing many of the observed GC-galaxy scaling relations.

We present homogeneous, sub-horizontal branch photometry of twenty dwarf spheroidal satellite galaxies of M31 observed with the Hubble Space Telescope. Combining our new data for sixteen systems with archival data in the same filters for another four, we show that Andromeda dwarf spheroidal galaxies favor strikingly red horizontal branches or red clumps down to ~10^{4.2} Lsun (M_V ~ -5.8). The age-sensitivity of horizontal branch stars implies that a large fraction of the M31 dwarf galaxies have extended star formation histories (SFHs), and appear inconsistent with early star formation episodes that were rapidly shutdown. Systems fainter than ~10^{5.5} Lsun show the widest range in the ratios and morphologies of red and blue horizontal branches, indicative of both complex SFHs and a diversity in quenching timescales and/or mechanisms, which is qualitatively different from what is currently known for faint Milky Way (MW) satellites of comparable luminosities. Our findings bolster similar conclusions from recent deeper data for a handful of M31 dwarf galaxies. We discuss several sources for diversity of our data such as varying halo masses, patchy reionization, mergers/accretion, and the environmental influence of M31 and the Milky Way on the early evolution of their satellite populations. A detailed comparison between the histories of M31 and MW satellites would shed signifiant insight into the processes that drive the evolution of low-mass galaxies. Such a study will require imaging that reaches the oldest main sequence turnoffs for a significant number of M31 companions.

The Large Synoptic Survey Telescope (LSST) will enable revolutionary studies of galaxies, dark matter, and black holes over cosmic time. The LSST Galaxies Science Collaboration has identified a host of preparatory research tasks required to leverage fully the LSST dataset for extragalactic science beyond the study of dark energy. This Galaxies Science Roadmap provides a brief introduction to critical extragalactic science to be conducted ahead of LSST operations, and a detailed list of preparatory science tasks including the motivation, activities, and deliverables associated with each. The Galaxies Science Roadmap will serve as a guiding document for researchers interested in conducting extragalactic science in anticipation of the forthcoming LSST era.

We confirm the reality of the recently discovered Milky Way stellar cluster $\textit{Gaia}$ 1 using spectra acquired with the HERMES and AAOmega spectrographs of the Anglo-Australian Telescope. This cluster had been previously undiscovered due to its close angular proximity to Sirius, the brightest star in the sky at visual wavelengths. Our observations identified 41 cluster members, and yielded an overall metallicity of [Fe/H]$=-0.13\pm0.13$ and barycentric radial velocity of $v_r=58.30\pm0.22$ km/s. These kinematics provide a dynamical mass estimate of $12.9^{+4.6}_{-3.9}\times10^3$ M$_{\odot}$. Isochrone fits to $\textit{Gaia}$, 2MASS, and Pan-STARRS1 photometry indicate that $\textit{Gaia}$ 1 is an intermediate age ($\sim3$ Gyr) stellar cluster. Combining the spatial and kinematic data we calculate $\textit{Gaia}$ 1 has a circular orbit with a radius of about 12~kpc, but with a large out of plane motion: $z_\textrm{max}=1.1^{+0.4}_{-0.3}$ kpc. Clusters with such orbits are unlikely to survive long due to the number of plane passages they would experience.

We show that hard encounters in the central regions of globular clusters embedded in dark matter (DM) haloes necessarily lead to the formation of gravitationally-bound stellar envelopes that extend far beyond the nominal tidal radius of the system. Using statistical arguments and numerical techniques we derive the equilibrium distribution function of stars ejected from the centre of a non-divergent spherical potential. Independently of the velocity distribution with which stars are ejected, GC envelopes have density profiles that approach asymptotically $\rho\sim r^{-4}$ at large distances and become isothermal towards the centre. Adding a DM halo component leaves two clear-cut observational signatures: (i) a flattening, or slightly increase of the projected velocity dispersion profile at large distances, and (ii) an outer surface density profile that is systematically shallower than in models with no dark matter.

Gas infall and outflow are critical for determining the star formation rate and chemical evolution of galaxies but direct measurements of gas flows are diffcult to make. Young massive stars and HII regions in the halos of galaxies are potential tracers for accretion and/or outflows of gas. Gas phase abundances of three HII regions in the lower halos of the edge-on galaxies NGC 3628 and NGC 4522 are determined by analysing optical long-slit spectra. The observed regions have projected distances to the midplane of their host from 1.4 to 3 kpc. With the measured flux densities of the optical nebular emission lines, we derive the oxygen abundance 12 + log(O/H) for the three extraplanar HII regions. The analysis is based on one theoretical and two empirical strong-line calibration methods. The resulting oxygen abundances of the extraplanar HII regions are comparable to the disk HII regions in one case and a little lower in the other case. Since our results depend on the accuracy of the metallicity determinations, we critically discuss the difference of the calibration methods we applied and confirm previously noted offsets. From our measurements, we argue that these three extraplanar HII regions were formed in the disk or at least from disk material. We discuss the processes that could transport disk material into the lower halo of these systems and conclude that gravitational interaction with a companion galaxy is most likely for NGC 3628 while ram pressure is favoured in the case of NGC 4522.

We present K-band Multi-Object Spectrograph (KMOS) observations of 18 Red Supergiant (RSG) stars in the Sculptor Group galaxy NGC 55. Radial velocities are calculated and are shown to be in good agreement with previous estimates, confirming the supergiant nature of the targets and providing the first spectroscopically confirmed RSGs in NGC 55. Stellar parameters are estimated for 14 targets using the $J$-band analysis technique, making use of state-of-the-art stellar model atmospheres. The metallicities estimated confirm the low-metallicity nature of NGC 55, in good agreement with previous studies. This study provides an independent estimate of the metallicity gradient of NGC 55, in excellent agreement with recent results published using hot massive stars. In addition, we calculate luminosities of our targets and compare their distribution of effective temperatures and luminosities to other RSGs, in different environments, estimated using the same technique.

Ultrahigh energy cosmic ray air showers probe particle physics at energies beyond the reach of accelerators. Here we introduce a new method to test hadronic interaction models without relying on the absolute energy calibration, and apply it to events with primary energy 6-16 EeV (E_CM = 110-170 TeV), whose longitudinal development and lateral distribution were simultaneously measured by the Pierre Auger Observatory. The average hadronic shower is 1.33 +- 0.16 (1.61 +- 0.21) times larger than predicted using the leading LHC-tuned models EPOS-LHC (QGSJetII-04), with a corresponding excess of muons.

We present a comprehensive analysis of the structural properties and luminosities of the 23 dwarf spheroidal galaxies that fall within the footprint of the Pan-Andromeda Archaeological Survey (PAndAS). These dwarf galaxies represent the large majority of Andromeda's known satellite dwarf galaxies and cover a wide range in luminosity ($-11.6<M_V<-5.8$ or $10^{4.2}< L <10^{6.5} L_\odot$) and surface brightness ($25.1<\mu_0<29.3$ mag/arcsec$^2$). We confirm most previous measurements, but find And XIX to be significantly larger than before ($r_h=3065^{+1065}_{-935}$ pc, $M_V=-10.1^{+0.8}_{-0.4}$) and cannot derive parameters for And XXVII as it is likely not a bound stellar system. We also significantly revise downward the luminosities of And~XV and And~XVI, which are now $M_V\sim-7.5$ or $L\sim10^5 L_\odot$. Finally, we provide the first detailed analysis of Cas II/And XXX, a fairly faint system ($M_V=-8.0^{+0.4}_{-0.3}$) of typical size ($r_h=270\pm50$ pc), located in close proximity to the two bright elliptical dwarf galaxies NGC 147 & 185. Combined with the set of homogeneous distances published in an earlier contribution, our analysis dutifully tracks all relevant sources of uncertainty in the determination of the properties of the dwarf galaxies from the PAndAS photometric catalogue. We further publish the posterior probability distribution functions of all the parameters we fit for in the form of MCMC chains available online; these inputs should be used in any analysis that aims to remain truthful to the data and properly account for covariance between parameters.

We present a panoramic map of the entire Milky Way halo north of dec~-30 degrees (~30,000 deg^2), constructed by applying the matched-filter technique to the Pan-STARRS1 3Pi Survey dataset. Using single-epoch photometry reaching to g~22, we are sensitive to stellar substructures with heliocentric distances between 3.5 and ~35 kpc. We recover almost all previously-reported streams in this volume and demonstrate that several of these are significantly more extended than earlier datasets have indicated. In addition, we also report five new candidate stellar streams. One of these features appears significantly broader and more luminous than the others and is likely the remnant of a dwarf galaxy. The other four streams are consistent with a globular cluster origin, and three of these are rather short in projection (<10 degrees), suggesting that streams like Ophiuchus may not be that rare. Finally, a significant number of more marginal substructures are also revealed by our analysis; many of these features can also be discerned in matched-filter maps produced by other authors from SDSS data, and hence they are very likely to be genuine. However, the extant 3Pi data is currently too shallow to determine their properties or produce convincing CMDs. The global view of the Milky Way provided by Pan-STARRS1 provides further evidence for the important role of both globular cluster disruption and dwarf galaxy accretion in building the Milky Way's stellar halo.

The stellar halos of large galaxies represent a vital probe of the processes of galaxy evolution. They are the remnants of the initial bouts of star formation during the collapse of the proto-galactic cloud, coupled with imprint of ancient and on-going accretion events. Previously, we have reported the tentative detection of a possible, faint, extended stellar halo in the Local Group spiral, the Triangulum Galaxy (M33). However, the presence of substructure surrounding M33 made interpretation of this feature difficult. Here, we employ the final data set from the Pan-Andromeda Archaeological Survey (PAndAS), combined with an improved calibration and a newly derived contamination model for the region to revisit this claim. With an array of new fitting algorithms, fully accounting for contamination and the substantial substructure beyond the prominent stellar disk in M33, we reanalyse the surrounds to separate the signal of the stellar halo and the outer halo substructure. Using more robust search algorithms, we do not detect a large scale smooth stellar halo and place a limit on the maximum surface brightness of such a feature of ${\mu}_V$ = 35.5 mags per square arcsec, or a total halo luminosity of $L < 10^6L_{\odot}$.

We measure the energy emitted by extensive air showers in the form of radio emission in the frequency range from 30 to 80 MHz. Exploiting the accurate energy scale of the Pierre Auger Observatory, we obtain a radiation energy of 15.8 \pm 0.7 (stat) \pm 6.7 (sys) MeV for cosmic rays with an energy of 1 EeV arriving perpendicularly to a geomagnetic field of 0.24 G, scaling quadratically with the cosmic-ray energy. A comparison with predictions from state-of-the-art first-principle calculations shows agreement with our measurement. The radiation energy provides direct access to the calorimetric energy in the electromagnetic cascade of extensive air showers. Comparison with our result thus allows the direct calibration of any cosmic-ray radio detector against the well-established energy scale of the Pierre Auger Observatory.

The Auger Engineering Radio Array (AERA) is part of the Pierre Auger Observatory and is used to detect the radio emission of cosmic-ray air showers. These observations are compared to the data of the surface detector stations of the Observatory, which provide well-calibrated information on the cosmic-ray energies and arrival directions. The response of the radio stations in the 30 to 80 MHz regime has been thoroughly calibrated to enable the reconstruction of the incoming electric field. For the latter, the energy deposit per area is determined from the radio pulses at each observer position and is interpolated using a two-dimensional function that takes into account signal asymmetries due to interference between the geomagnetic and charge-excess emission components. The spatial integral over the signal distribution gives a direct measurement of the energy transferred from the primary cosmic ray into radio emission in the AERA frequency range. We measure 15.8 MeV of radiation energy for a 1 EeV air shower arriving perpendicularly to the geomagnetic field. This radiation energy -- corrected for geometrical effects -- is used as a cosmic-ray energy estimator. Performing an absolute energy calibration against the surface-detector information, we observe that this radio-energy estimator scales quadratically with the cosmic-ray energy as expected for coherent emission. We find an energy resolution of the radio reconstruction of 22% for the data set and 17% for a high-quality subset containing only events with at least five radio stations with signal.

The azimuthal asymmetry in the risetime of signals in Auger surface detector stations is a source of information on shower development. The azimuthal asymmetry is due to a combination of the longitudinal evolution of the shower and geometrical effects related to the angles of incidence of the particles into the detectors. The magnitude of the effect depends upon the zenith angle and state of development of the shower and thus provides a novel observable, $(\sec \theta)_\mathrm{max}$, sensitive to the mass composition of cosmic rays above $3 \times 10^{18}$ eV. By comparing measurements with predictions from shower simulations, we find for both of our adopted models of hadronic physics (QGSJETII-04 and EPOS-LHC) an indication that the mean cosmic-ray mass increases slowly with energy, as has been inferred from other studies. However, the mass estimates are dependent on the shower model and on the range of distance from the shower core selected. Thus the method has uncovered further deficiencies in our understanding of shower modelling that must be resolved before the mass composition can be inferred from $(\sec \theta)_\mathrm{max}$.

The Pierre Auger Observatory has begun a major Upgrade of its already impressive capabilities, with an emphasis on improved mass composition determination using the surface detectors of the Observatory. Known as AugerPrime, the upgrade will include new 4 m$^2$ plastic scintillator detectors on top of all 1660 water-Cherenkov detectors, updated and more flexible surface detector electronics, a large array of buried muon detectors, and an extended duty cycle for operations of the fluorescence detectors. This Preliminary Design Report was produced by the Collaboration in April 2015 as an internal document and information for funding agencies. It outlines the scientific and technical case for AugerPrime. We now release it to the public via the arXiv server. We invite you to review the large number of fundamental results already achieved by the Observatory and our plans for the future.

Tidal streams from existing and destroyed satellite galaxies populate the outer regions of the Andromeda galaxy (M31). This inhomogeneous debris can be studied without many of the obstacles that plague Milky Way research. We review the history of tidal stream research in M31, and in its main satellite galaxies. We highlight the numerous tidal streams observed around M31, some of which reside at projected distances of up to 120 kpc from the center of this galaxy. Most notable is the Giant Stellar Stream, a signature of the most recent significant accretion event in the M31 system. This event involved an early-type progenitor of ~10^9 solar masses that came within a few kpc of M31's center roughly a gigayear ago; almost all of the inner halo debris (within 50 kpc) in M31 can be tied either directly or indirectly to this event. We draw attention to the fact that most of M31's outer halo globular clusters lie preferentially on tidal streams and discuss the potential this offers to use these systems as probes of the accretion history. Tidal features observed around M33, M32, NGC 205 and NGC 147 are also reviewed. We conclude by discussing future prospects for this field.

We present a renewed look at M31's Giant Stellar Stream along with the nearby structures Stream C and Stream D, exploiting a new algorithm capable of fitting to the red giant branch (RGB) of a structure in both colour and magnitude space. Using this algorithm, we are able to generate probability distributions in distance, metallicity and RGB width for a series of subfields spanning these structures. Specifically, we confirm a distance gradient of approximately 20 kpc per degree along a 6 degree extension of the Giant Stellar Stream, with the farthest subfields from M31 lying ~ 120 kpc more distant than the inner-most subfields. Further, we find a metallicity that steadily increases from -0.7^{+0.1}_{-0.1} dex to -0.2^{+0.2}_{-0.1} dex along the inner half of the stream before steadily dropping to a value of -1.0^{+0.2}_{-0.2} dex at the farthest reaches of our coverage. The RGB width is found to increase rapidly from 0.4^{+0.1}_{-0.1} dex to 1.1^{+0.2}_{-0.1} dex in the inner portion of the stream before plateauing and decreasing marginally in the outer subfields of the stream. In addition, we estimate Stream C to lie at a distance between 794 and 862 kpc and Stream D between 758 kpc and 868 kpc. We estimate the median metallicity of Stream C to lie in the range -0.7 to -1.6 dex and a metallicity of -1.1^{+0.3}_{-0.2} dex for Stream D. RGB widths for the two structures are estimated to lie in the range 0.4 to 1.2 dex and 0.3 to 0.7 dex respectively. In total, measurements are obtained for 19 subfields along the Giant Stellar Stream, 4 along Stream C, 5 along Stream D and 3 general M31 spheroid fields for comparison. We thus provide a higher resolution coverage of the structures in these parameters than has previously been available in the literature.

This paper presents the results of different searches for correlations between very high-energy neutrino candidates detected by IceCube and the highest-energy cosmic rays measured by the Pierre Auger Observatory and the Telescope Array. We first consider samples of cascade neutrino events and of high-energy neutrino-induced muon tracks, which provided evidence for a neutrino flux of astrophysical origin, and study their cross-correlation with the ultrahigh-energy cosmic ray (UHECR) samples as a function of angular separation. We also study their possible directional correlations using a likelihood method stacking the neutrino arrival directions and adopting different assumptions on the size of the UHECR magnetic deflections. Finally, we perform another likelihood analysis stacking the UHECR directions and using a sample of through-going muon tracks optimized for neutrino point-source searches with sub-degree angular resolution. No indications of correlations at discovery level are obtained for any of the searches performed. The smallest of the p-values comes from the search for correlation between UHECRs with IceCube high-energy cascades, a result that should continue to be monitored.

We present the first spectroscopic analysis of the faint and compact stellar system Draco II (Dra II, M_V=-2.9+/-0.8, r_h=19^{+8}_{-6} pc), recently discovered in the Pan-STARRS1 3{\pi} survey. The observations, conducted with DEIMOS on the Keck II telescope, establish some of its basic characteristics: the velocity data reveal a narrow peak with 9 member stars at a systemic heliocentric velocity <v_r>=-347.6^{+1.7}_{-1.8} km/s, thereby confirming Dra II is a satellite of the Milky Way; we infer a velocity dispersion with \sigma_{vr}=2.9+/-2.1 km/s (<8.4 km/s at the 95% confidence level), which implies log_{10}(M_{1/2})=5.5^{+0.4}_{-0.6} and log_{10}((M/L)_{1/2})=2.7^{+0.5}_{-0.8}, in Solar units; furthermore, very weak Calcium triplet lines in the spectra of the high signal-to-noise member stars imply [Fe/H]<-2.1, whilst variations in the line strengths of two stars with similar colours and magnitudes suggest a metallicity spread in Dra II. These new data cannot clearly discriminate whether Draco II is a star cluster or amongst the faintest, most compact, and closest dwarf galaxies. However, the sum of the three --- individually inconclusive --- pieces of evidence presented here, seems to favour the dwarf galaxy interpretation.

We present a study of the recently discovered compact stellar system Triangulum II. From observations conducted with the DEIMOS spectrograph on Keck II, we obtained spectra for 13 member stars that follow the CMD features of this very faint stellar system and include two bright red giant branch stars. Tri II has a very negative radial velocity (<v_r>=-383.7^{+3.0}_{-3.3} km/s) that translates to <v_{r,gsr}> ~ -264 km/s and confirms it is a Milky Way satellite. We show that, despite the small data set, there is evidence that Tri II has complex internal kinematics. Its radial velocity dispersion increases from 4.4^{+2.8}_{-2.0} km/s in the central 2' to 14.1^{+5.8}_{-4.2} km/s outwards. The velocity dispersion of the full sample is inferred to be \sigma_{vr}=9.9^{+3.2}_{-2.2} km/s. From the two bright RGB member stars we measure an average metallicity <[Fe/H]>=-2.6+/-0.2, placing Tri II among the most metal-poor Milky Way dwarf galaxies. In addition, the spectra of the fainter member stars exhibit differences in their line widths that could be the indication of a metallicity dispersion in the system. All these properties paint a complex picture for Tri II, whose nature and current state are largely speculative. The inferred metallicity properties of the system however lead us to favor a scenario in which Tri II is a dwarf galaxy that is either disrupting or embedded in a stellar stream.

The Ophiuchus stellar stream presents a dynamical puzzle: its old stellar populations ($\sim 12$ Gyr) cannot be reconciled with (1) its orbit in a simple model for the Milky Way potential and (2) its short angular extent, both of which imply that the observed stream formed within the last $<1$ Gyr. Recent theoretical work has shown that streams on chaotic orbits may abruptly fan out near their apparent ends; stars in these fans are dispersed in both position and velocity and may be difficult to associate with the stream. Here we present the first evidence of such stream-fanning in the Ophiuchus stream, traced by four blue horizontal branch (BHB) stars beyond the apparent end of the stream. These stars stand out from the background by their high velocities ($v_{\rm los} > 230$ km s$^{-1}$) against $\sim 40$ other stars: their velocities are comparable to those of the stream, but would be exceptional if they were unrelated halo stars. Their positions and velocities are, however, inconsistent with simple extrapolation of the observed cold, high-density portion of the stream. These observations suggest that stream-fanning may be a real, observable effect and, therefore, that Ophiuchus may be on a chaotic orbit. They also show that the Ophiuchus stream is more extended and hence dynamically older than previously thought, easing the stellar population vs. dynamical age tension.

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.

Neutrinos in the cosmic ray flux with energies near 1 EeV and above are detectable with the Surface Detector array of the Pierre Auger Observatory. We report here on searches through Auger data from 1 January 2004 until 20 June 2013. No neutrino candidates were found, yielding a limit to the diffuse flux of ultra-high energy neutrinos that challenges the Waxman-Bahcall bound predictions. Neutrino identification is attempted using the broad time-structure of the signals expected in the SD stations, and is efficiently done for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for "Earth-skimming" neutrino interactions in the case of tau neutrinos. In this paper the searches for downward-going neutrinos in the zenith angle bins $60^\circ-75^\circ$ and $75^\circ-90^\circ$ as well as for upward-going neutrinos, are combined to give a single limit. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high energy neutrinos with an $E^{-2}$ spectrum in the energy range $1.0 \times 10^{17}$ eV - $2.5 \times 10^{19}$ eV is $E_\nu^2 dN_\nu/dE_\nu < 6.4 \times 10^{-9}~ {\rm GeV~ cm^{-2}~ s^{-1}~ sr^{-1}}$.

A measurement of the cosmic-ray spectrum for energies exceeding $4{\times}10^{18}$ eV is presented, which is based on the analysis of showers with zenith angles greater than $60^{\circ}$ detected with the Pierre Auger Observatory between 1 January 2004 and 31 December 2013. The measured spectrum confirms a flux suppression at the highest energies. Above $5.3{\times}10^{18}$ eV, the "ankle", the flux can be described by a power law $E^{-\gamma}$ with index $\gamma=2.70 \pm 0.02 \,\text{(stat)} \pm 0.1\,\text{(sys)}$ followed by a smooth suppression region. For the energy ($E_\text{s}$) at which the spectral flux has fallen to one-half of its extrapolated value in the absence of suppression, we find $E_\text{s}=(5.12\pm0.25\,\text{(stat)}^{+1.0}_{-1.2}\,\text{(sys)}){\times}10^{19}$ eV.

We present the first detailed analysis of the East Cloud, a highly disrupted diffuse stellar substructure in the outer halo of M31. The core of the substructure lies at a projected distance of $\sim100$ kpc from the centre of M31 in the outer halo, with possible extensions reaching right into the inner halo. Using Pan-Andromeda Archaeological Survey photometry of red giant branch stars, we measure the distance, metallicity and brightness of the cloud. Using Hubble Space Telescope data, we independently measure the distance and metallicity to the two globular clusters coincident with the East Cloud core, PA-57 and PA-58, and find their distances to be consistent with the cloud. Four further globular clusters coincident with the substructure extensions are identified as potentially associated. Combining the analyses, we determine a distance to the cloud of $814^{+20}_{-9}$ kpc, a metallicity of $[Fe/H] = -1.2\pm0.1$, and a brightness of $M_V = -10.7\pm0.4$ mag. Even allowing for the inclusion of the potential extensions, this accounts for less than $20$ per cent of the progenitor luminosity implied by the luminosity-metallicity relation. Using the updated techniques developed for this analysis, we also refine our estimates of the distance and brightness of the South-West Cloud, a separate substructure analyzed in the previous work in this series.

We have conducted three searches for correlations between ultra-high energy cosmic rays detected by the Telescope Array and the Pierre Auger Observatory, and high-energy neutrino candidate events from IceCube. Two cross-correlation analyses with UHECRs are done: one with 39 cascades from the IceCube `high-energy starting events' sample and the other with 16 high-energy `track events'. The angular separation between the arrival directions of neutrinos and UHECRs is scanned over. The same events are also used in a separate search using a maximum likelihood approach, after the neutrino arrival directions are stacked. To estimate the significance we assume UHECR magnetic deflections to be inversely proportional to their energy, with values $3^\circ$, $6^\circ$ and $9^\circ$ at 100 EeV to allow for the uncertainties on the magnetic field strength and UHECR charge. A similar analysis is performed on stacked UHECR arrival directions and the IceCube sample of through-going muon track events which were optimized for neutrino point-source searches.

Joint contributions of the Pierre Auger Collaboration and the Telescope Array Collaboration to the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands.

We present the discovery of three new Milky Way satellites from our search for compact stellar overdensities in the photometric catalog of the Panoramic Survey Telescope and Rapid Response System 1 (Pan-STARRS 1, or PS1) 3pi survey. The first satellite, Laevens 3, is located at a heliocentric distance of d=67+/-3 kpc. With a total magnitude of Mv=-4.4+/-0.3 and a half-light radius rh=7+/-2 pc, its properties resemble those of outer halo globular clusters. The second system, Draco II/Laevens 4 (Dra II), is a closer and fainter satellite (d~20 kpc, Mv =-2.9+/-0.8), whose uncertain size (rh = 19 +8/-6 pc) renders its classification difficult without kinematic information; it could either be a faint and extended globular cluster or a faint and compact dwarf galaxy. The third satellite, Sagittarius II/Laevens 5 (Sgr II), has an ambiguous nature as it is either the most compact dwarf galaxy or the most extended globular cluster in its luminosity range (rh = 37 +9/-8 pc and Mv=-5.2+/-0.4). At a heliocentric distance of 67+/-5 kpc, this satellite lies intriguingly close to the expected location of the trailing arm of the Sagittarius stellar stream behind the Sagittarius dwarf spheroidal galaxy (Sgr dSph). If confirmed through spectroscopic follow up, this connection would locate this part of the trailing arm of the Sagittarius stellar stream that has so far gone undetected. It would further suggest that Sgr II was brought into the Milky Way halo as a satellite of the Sgr dSph.

Contributions of the Pierre Auger Collaboration to the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands

We present deep Hubble Space Telescope Advanced Camera for Surveys observations of the stellar populations in two fields lying at 20 and 23 kpc from the centre of M31 along the south-west semi-major axis. These data enable the construction of colour-magnitude diagrams reaching the oldest main-sequence turn-offs (~13 Gyr) which, when combined with another field at 25 kpc from our previous work, we use to derive the first precision constraints on the spatially-resolved star formation history of the M31 disc. The star formation rates exhibit temporal as well as field-to-field variations, but are generally always within a factor of two of their time average. There is no evidence of inside-out growth over the radial range probed. We find a median age of ~7.5 Gyr, indicating that roughly half of the stellar mass in the M31 outer disc was formed before z ~ 1. We also find that the age-metallicity relations (AMRs) are smoothly increasing from [Fe/H]~-0.4 to solar metallicity between 10 and 3 Gyr ago, contrary to the flat AMR of the Milky Way disc at a similar number of scale lengths. Our findings provide insight on the roles of stellar feedback and radial migration in the formation and evolution of large disc galaxies.

The Ophiuchus stream is a recently discovered stellar tidal stream in the Milky Way. We present high-quality spectroscopic data for 14 stream member stars obtained using the Keck and MMT telescopes. We confirm the stream as a fast moving ($v_{los}\sim290$ km s$^{-1}$), kinematically cold group ($\sigma_{v_{los}}\lesssim1$ km s$^{-1}$) of $\alpha$-enhanced and metal-poor stars (${\rm [\alpha/Fe]\sim0.4}$ dex, ${\rm [Fe/H]\sim-2.0}$ dex). Using a probabilistic technique, we model the stream simultaneously in line-of-sight velocity, color-magnitude, coordinate, and proper motion space, and so determine its distribution in 6D phase-space. We find that that the stream extends in distance from 7.5 to 9 kpc from the Sun; it is 50 times longer than wide, merely appearing highly foreshortened in projection. The analysis of the stellar population contained in the stream suggests that it is $\sim12$ Gyr old, and that its initial stellar mass was $\sim2\times10^4$ $M_{\odot}$ (or at least $\gtrsim7\times10^3$ $M_{\odot}$). Assuming a fiducial Milky Way potential, we fit an orbit to the stream which matches the observed phase-space distribution, except for some tension in the proper motions: the stream has an orbital period of $\sim350$ Myr, and is on a fairly eccentric orbit ($e\sim0.66$) with a pericenter of $\sim3.5$ kpc and an apocenter of $\sim17$ kpc. The phase-space structure and stellar population of the stream show that its progenitor must have been a globular cluster that was disrupted only $\sim240$ Myr ago. We do not detect any significant overdensity of stars along the stream that would indicate the presence of a progenitor, and conclude that the stream is all that is left of the progenitor.

We present the first results of a wide-field mapping survey of the M81 group conducted with Hyper Suprime-Cam on the Subaru Telescope. Our deep photometry reaches $\sim2$ magnitudes below the tip of the red giant branch (RGB) and reveals the spatial distribution of both old and young stars over an area of $\sim 100\times115$ kpc at the distance of M81. The young stars ($\sim30-160$ Myr old) closely follow the neutral hydrogen distribution and can be found in a stellar stream between M81 and NGC\,3077 and in numerous outlying stellar associations, including the known concentrations of Arp's Loop, Holmberg\,IX, an arc in the halo of M82, BK3N, and the Garland. Many of these groupings do not have counterparts in the RGB maps, suggesting they may be genuinely young systems. Our survey also reveals for the first time the very extended ($\geq 2\times \rm{R_{25}}$) halos of RGB stars around M81, M82 and NGC\,3077, as well as faint tidal streams that link these systems. The halos of M82 and NGC\,3077 exhibit highly disturbed morphologies, presumably a consequence of the recent gravitational encounter and their ongoing disruption. While the halos of M81, NGC\,3077 and the inner halo of M82 have the similar $(g-i)_{0}$ colors, the outer halo of M82 is significantly bluer indicating it is more metal-poor. Remarkably, our deep panoramic view of the M81 group demonstrates that the complexity long-known to be present in HI is equally matched in the low surface brightness stellar component.

We present a comprehensive analysis of the globular cluster (GC) system of the Local Group dwarf irregular galaxy NGC 6822. Our study is based on homogeneous optical and near-IR photometry, as well as long-slit spectroscopic observations which are used to determine new radial velocities for 6 GCs, two of which had no previous spectroscopic information. We construct optical-near IR colour-colour diagrams and through comparison to simple stellar population models infer that the GCs have old ages consistent with being 9 Gyr or older, while their metallicities are in the range between -1.6 < [Fe/H] < -0.4. We conduct a kinematic analysis of the GC population and find tentative evidence for weak net rotation of the GC system, in the same sense as that exhibited by the underlying spheroid. The most likely amplitude of rotation is ~10 km/s, approximately half the magnitude of the observed velocity dispersion. Finally, we use the GCs to estimate the dynamical mass of NGC 6822 within 11 kpc and we formally find it to be in the range between (3-4)10^9 Msun. This implies an overall mass-to-light ratio in the range of ~ 30-40 and indicates that NGC 6822 is highly dark matter dominated. The mass and the corresponding mass-to-light ratio estimates are affected by various additional systematic effects due to limitations of the data and the model that are not necessary reflected in the formal uncertainties.

The Gaia-ESO Survey (GES) is a large public spectroscopic survey at the European Southern Observatory Very Large Telescope. A key aim is to provide precise radial velocities (RVs) and projected equatorial velocities (v sin i) for representative samples of Galactic stars, that will complement information obtained by the Gaia astrometry satellite. We present an analysis to empirically quantify the size and distribution of uncertainties in RV and v sin i using spectra from repeated exposures of the same stars. We show that the uncertainties vary as simple scaling functions of signal-to-noise ratio (S/N) and v sin i, that the uncertainties become larger with increasing photospheric temperature, but that the dependence on stellar gravity, metallicity and age is weak. The underlying uncertainty distributions have extended tails that are better represented by Student's t-distributions than by normal distributions. Parametrised results are provided, that enable estimates of the RV precision for almost all GES measurements, and estimates of the v sin i precision for stars in young clusters, as a function of S/N, v sin i and stellar temperature. The precision of individual high S/N GES RV measurements is 0.22-0.26 km/s, dependent on instrumental configuration.

This White Paper presents the scientific motivations for a multi-object spectrograph (MOS) on the European Extremely Large Telescope (E-ELT). The MOS case draws on all fields of contemporary astronomy, from extra-solar planets, to the study of the halo of the Milky Way and its satellites, and from resolved stellar populations in nearby galaxies out to observations of the earliest 'first-light' structures in the partially-reionised Universe. The material presented here results from thorough discussions within the community over the past four years, building on the past competitive studies to agree a common strategy toward realising a MOS capability on the E-ELT. The cases have been distilled to a set of common requirements which will be used to define the MOSAIC instrument, entailing two observational modes ('high multiplex' and 'high definition'). When combined with the unprecedented sensitivity of the E-ELT, MOSAIC will be the world's leading MOS facility. In analysing the requirements we also identify a high-multiplex MOS for the longer-term plans for the E-ELT, with an even greater multiplex (>1000 targets) to enable studies of large-scale structures in the high-redshift Universe. Following the green light for the construction of the E-ELT the MOS community, structured through the MOSAIC consortium, is eager to realise a MOS on the E-ELT as soon as possible. We argue that several of the most compelling cases for ELT science, in highly competitive areas of modern astronomy, demand such a capability. For example, MOS observations in the early stages of E-ELT operations will be essential for follow-up of sources identified by the James Webb Space Telescope (JWST). In particular, multi-object adaptive optics and accurate sky subtraction with fibres have both recently been demonstrated on sky, making fast-track development of MOSAIC feasible.

We present the discovery of a faint Milky Way satellite, Laevens 2/Triangulum II, found in the Panoramic Survey Telescope And Rapid Response System (Pan-STARRS 1) 3 pi imaging data and confirmed with follow-up wide-field photometry from the Large Binocular Cameras. The stellar system, with an absolute magnitude of M_V=-1.8 +/-0.5, a heliocentric distance of 30 +2/-2 kpc, and a half-mass radius of 34 +9/-8 pc, shows remarkable similarity to faint, nearby, small satellites such as Willman 1, Segue 1, Segue 2, and Bo\"otes II. The discovery of Laevens 2/Triangulum II further populates the region of parameter space for which the boundary between dwarf galaxies and globular clusters becomes tenuous. Follow-up spectroscopy will ultimately determine the nature of this new satellite, whose spatial location hints at a possible connection with the complex Triangulum-Andromeda stellar structures.

A small fraction of the halo field is made up of stars that share the light element (Z<=13) anomalies characteristic of second generation globular cluster (GC) stars. The ejected stars shed light on the formation of the Galactic halo by tracing the dynamical history of the clusters, which are believed to have once been more massive. Some of these ejected stars are expected to show strong Al enhancement at the expense of shortage of Mg, but until now no such star has been found. We search for outliers in the Mg and Al abundances of the few hundreds of halo field stars observed in the first eighteen months of the Gaia-ESO public spectroscopic survey. One halo star at the base of the red giant branch, here referred to as 22593757-4648029 is found to have [Mg/Fe]=-0.36+-0.04 and [Al/Fe]=0.99+-0.08, which is compatible with the most extreme ratios detected in GCs so far. We compare the orbit of 22593757-4648029 to GCs of similar metallicity and find it unlikely that this star has been tidally stripped with low ejection velocity from any of the clusters. However, both chemical and kinematic arguments render it plausible that the star has been ejected at high velocity from the anomalous GC omega Centauri within the last few billion years. We cannot rule out other progenitor GCs, because some may have disrupted fully, and the abundance and orbital data are inadequate for many of those that are still intact.

We present near-IR spectroscopy of red supergiant (RSG) stars in NGC 6822, obtained with the new VLT-KMOS instrument. From comparisons with model spectra in the J-band we determine the metallicity of 11 RSGs, finding a mean value of [Z] = -0.52 $\pm$ 0.21 which agrees well with previous abundance studies of young stars and HII regions. We also find an indication for a low-significance abundance gradient within the central 1 kpc. We compare our results with those derived from older stellar populations and investigate the difference using a simple chemical evolution model. By comparing the physical properties determined for RSGs in NGC 6822 with those derived using the same technique in the Galaxy and the Magellanic Clouds, we show that there appears to be no significant temperature variation of RSGs with respect to metallicity, in contrast with recent evolutionary models.

This paper describes the analysis of UVES and GIRAFFE spectra acquired by the Gaia-ESO Public Spectroscopic Survey in the fields of young clusters whose population includes pre-main sequence (PMS) stars. Both methods that have been extensively used in the past and new ones developed in the contest of the Gaia-ESO survey enterprise are available and used. The internal precision of these quantities is estimated by inter-comparing the results obtained by such different methods, while the accuracy is estimated by comparison with independent external data, like effective temperature and surface gravity derived from angular diameter measurements, on a sample of benchmarks stars. Specific strategies are implemented to deal with fast rotation, accretion signatures, chromospheric activity, and veiling. The analysis carried out on spectra acquired in young clusters' fields during the first 18 months of observations, up to June 2013, is presented in preparation of the first release of advanced data products. Stellar parameters obtained with the higher resolution and larger wavelength coverage from UVES are reproduced with comparable accuracy and precision using the smaller wavelength range and lower resolution of the GIRAFFE setup adopted for young stars, which allows us to provide with confidence stellar parameters for the much larger GIRAFFE sample. Precisions are estimated to be $\approx$ 120 K r.m.s. in Teff, $\approx$0.3 dex r.m.s. in logg, and $\approx$0.15 dex r.m.s. in [Fe/H], for both the UVES and GIRAFFE setups.

Detailed chemical abundances are presented for seven M31 outer halo globular clusters (with projected distances from M31 greater than 30 kpc), as derived from high resolution integrated light spectra taken with the Hobby Eberly Telescope. Five of these clusters were recently discovered in the Pan-Andromeda Archaeological Survey (PAndAS)---this paper presents the first determinations of integrated Fe, Na, Mg, Ca, Ti, Ni, Ba, and Eu abundances for these clusters. Four of the target clusters (PA06, PA53, PA54, and PA56) are metal-poor ([Fe/H] < -1.5), alpha-enhanced (though they are possibly less alpha-enhanced than Milky Way stars at the 1 sigma level), and show signs of star-to-star Na and Mg variations. The other three globular clusters (H10, H23, and PA17) are more metal rich, with metallicities ranging from [Fe/H] = -1.4 to -0.9. While H23 is chemically similar to Milky Way field stars, Milky Way globular clusters, and other M31 clusters, H10 and PA17 have moderately low [Ca/Fe], compared to Milky Way field stars and clusters. Additionally, PA17's high [Mg/Ca] and [Ba/Eu] ratios are distinct from Milky Way stars, and are in better agreement with the stars and clusters in the Large Magellanic Cloud (LMC). None of the clusters studied here can be conclusively linked to any of the identified streams from PAndAS; however, based on their locations, kinematics, metallicities, and detailed abundances, the most metal-rich PAndAS clusters H23 and PA17 may be associated with the progenitor of the Giant Stellar Stream, H10 may be associated with the SW Cloud, and PA53 and PA56 may be associated with the Eastern Cloud.

The thin, extended planes of satellite galaxies detected around both the Milky Way and Andromeda are not a natural prediction of the LCDM paradigm. Galaxies in these distinct planes may have formed and evolved in a different way (e.g., tidally) to their off-plane neighbours. If this were the case, one would expect the on- and off-plane dwarf galaxies in Andromeda to have experienced different evolutionary histories, which should be reflected by the chemistries, dynamics, and star formation histories of the two populations. In this work, we present new, robust kinematic observations for 2 on-plane M31 dSphs (And XVI and XVII) and compile and compare all available observational metrics for the on- and off-plane dwarfs to search for a signal that would corroborate such a hypothesis. We find that, barring their spatial alignment, the on- and off-plane Andromeda dwarf galaxies are indistinguishable from one another, arguing against vastly different formative and evolutionary histories for these two populations.

While wide-field surveys of M31 have revealed much substructure at large radii, understanding the nature and origin of this material is not straightforward from morphology alone. Using deep HST/ACS data, we have derived further constraints in the form of quantitative star formation histories (SFHs) for 14 inner halo fields which sample diverse substructures. In agreement with our previous analysis of colour-magnitude diagram morphologies, we find the resultant behaviours can be broadly separated into two categories. The SFHs of 'disc-like' fields indicate that most of their mass has formed since z~1, with one quarter of the mass formed in the last 5 Gyr. We find 'stream-like' fields to be on average 1.5 Gyr older, with <10 percent of their stellar mass formed within the last 5 Gyr. These fields are also characterised by an age--metallicity relation showing rapid chemical enrichment to solar metallicity by z=1, suggestive of an early-type progenitor. We confirm a significant burst of star formation 2 Gyr ago, discovered in our previous work, in all the fields studied here. The presence of these young stars in our most remote fields suggests that they have not formed in situ but have been kicked-out from through disc heating in the recent past.

We exploit data from the Pan-Andromeda Archaeological Survey (PAndAS) to study the extended structures of M31's dwarf elliptical companions, NGC147 and NGC185. Our wide-field, homogeneous photometry allows to construct deep colour-magnitude diagrams (CMDs) which reach down to $\sim3$ mag below the red giant branch (RGB) tip. We trace the stellar components of the galaxies to surface brightness of $\mu_g \sim 32$ mag arcsec$^{-2}$ and show they have much larger extents ($\sim5$ kpc radii) than previously recognised. While NGC185 retains a regular shape in its peripheral regions, NGC147 exhibits pronounced isophotal twisting due to the emergence of symmetric tidal tails. We fit single Sersic models to composite surface brightness profiles constructed from diffuse light and star counts and find that NGC147 has an effective radius almost 3 times that of NGC185. In both cases, the effective radii that we calculate are larger by a factor of $\sim2$ compared to most literature values. We also calculate revised total magnitudes of $M_g=-15.36\pm0.04$ for NGC185 and $M_g=-16.36\pm0.04$ for NGC147. Using photometric metallicities computed for RGB stars, we find NGC185 to exhibit a metallicity gradient of [Fe/H]$\sim-0.15$ dex/kpc over the radial range 0.125 to 0.5 deg. On the other hand, NGC147 exhibits almost no metallicity gradient, $\sim-0.02$ dex/kpc from 0.2 to 0.6 deg. The differences in the structure and stellar populations in the outskirts of these systems suggest that tidal influences have played an important role in governing the evolution of NGC147.

Spherical harmonic moments are well-suited for capturing anisotropy at any scale in the flux of cosmic rays. An unambiguous measurement of the full set of spherical harmonic coefficients requires full-sky coverage. This can be achieved by combining data from observatories located in both the northern and southern hemispheres. To this end, a joint analysis using data recorded at the Telescope Array and the Pierre Auger Observatory above $10^{19}$ eV is presented in this work. The resulting multipolar expansion of the flux of cosmic rays allows us to perform a series of anisotropy searches, and in particular to report on the angular power spectrum of cosmic rays above $10^{19}$ eV. No significant deviation from isotropic expectations is found throughout the analyses performed. Upper limits on the amplitudes of the dipole and quadrupole moments are derived as a function of the direction in the sky, varying between 7% and 13% for the dipole and between 7% and 10% for a symmetric quadrupole.