Researcher profile

Michael J. Irwin

Michael J. Irwin contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
18works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2022arXiv

Forward and Back: Kinematics of the Palomar 5 Tidal Tails

The tidal tails of Palomar 5 (Pal 5) have been the focus of many spectroscopic studies in an attempt to identify individual stars lying along the stream and characterise their kinematics. The well-studied trailing tail has been explored out to a distance of 15^\text{o} from the cluster centre, while less than four degrees have been examined along the leading tail. In this paper, we present results of a spectroscopic study of two fields along the leading tail that we have observed with the AAOmega spectrograph on the Anglo-Australian telescope. One of these fields lies roughly 7^\text{o} along the leading tail, beyond what has been previously been explored spectroscopically. Combining our measurements of kinematics and line strengths with Pan-STARRS1 photometric data and Gaia EDR3 astrometry, we adopt a probabilistic approach to identify 16 stars with high probability of belonging to the Pal 5 stream. Eight of these stars lie in the outermost field and their sky positions confirm the presence of ``fanning'' in the leading arm. We also revisit previously-published radial velocity studies and incorporate Gaia EDR3 astrometry to remove interloping field stars. With a final sample of 109 {\it bona fide} Pal 5 cluster and tidal stream stars, we characterise the 3D kinematics along the the full extent of the system. We provide this catalogue for future modeling work.

preprint2016arXiv

Major Substructure in the M31 Outer Halo: Distances and Metallicities along the Giant Stellar Stream

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.

preprint2014arXiv

Comparing the observable properties of dwarf galaxies on and off the Andromeda plane

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.

preprint2014arXiv

Precise time-series photometry for the Kepler-2.0 mission

The recently approved NASA K2 mission has the potential to multiply by an order of magnitude the number of short-period transiting planets found by Kepler around bright and low-mass stars, and to revolutionise our understanding of stellar variability in open clusters. However, the data processing is made more challenging by the reduced pointing accuracy of the satellite, which has only two functioning reaction wheels. We present a new method to extract precise light curves from K2 data, combining list-driven, soft-edged aperture photometry with a star-by-star correction of systematic effects associated with the drift in the roll-angle of the satellite about its boresight. The systematics are modelled simultaneously with the stars' intrinsic variability using a semi-parametric Gaussian process model. We test this method on a week of data collected during an engineering test in January 2014, perform checks to verify that our method does not alter intrinsic variability signals, and compute the precision as a function of magnitude on long-cadence (30-min) and planetary transit (2.5-hour) timescales. In both cases, we reach photometric precisions close to the precision reached during the nominal Kepler mission for stars fainter than 12th magnitude, and between 40 and 80 parts per million for brighter stars. These results confirm the bright prospects for planet detection and characterisation, asteroseismology and stellar variability studies with K2. Finally, we perform a basic transit search on the light curves, detecting 2 bona fide transit-like events, 7 detached eclipsing binaries and 13 classical variables.

preprint2014arXiv

The PAndAS Field of Streams: stellar structures in the Milky Way halo toward Andromeda and Triangulum

We reveal the highly structured nature of the Milky Way stellar halo within the footprint of the PAndAS photometric survey from blue main sequence and main sequence turn-off stars. We map no fewer than five stellar structures within a heliocentric range of ~5 to 30 kpc. Some of these are known (the Monoceros Ring, the Pisces/Triangulum globular cluster stream), but we also uncover three well-defined stellar structures that could be, at least partly, responsible for the so-called Triangulum/Andromeda and Triangulum/Andromeda 2 features. In particular, we trace a new faint stellar stream located at a heliocentric distance of ~17 kpc. With a surface brightness of Σ_V ~ 32-32.5 mag/arcsec^2, it follows an orbit that is almost parallel to the Galactic plane north of M31 and has so far eluded surveys of the Milky Way halo as these tend to steer away from regions dominated by the Galactic disk. Investigating our follow-up spectroscopic observations of PAndAS, we serendipitously uncover a radial velocity signature from stars that have colors and magnitudes compatible with the stream. From the velocity of eight likely member stars, we show that this stellar structure is dynamically cold, with an unresolved velocity dispersion that is lower than 7.1 km/s at the 90-percent confidence level. Along with the width of the stream (300-650 pc), its dynamics points to a dwarf-galaxy-accretion origin. The numerous stellar structures we can map in the Milky Way stellar halo between 5 and 30 kpc and their varying morphology is a testament to the complex nature of the stellar halo at these intermediate distances.

preprint2013arXiv

A kinematic study of the Andromeda dwarf spheroidal system

We present a homogeneous kinematic analysis of red giant branch stars within 18 of the 28 Andromeda dwarf spheroidal (dSph) galaxies, obtained using the Keck I LRIS and Keck II DEIMOS spectrographs. Based on their g-i colors (taken with the CFHT MegaCam imager), physical positions on the sky, and radial velocities, we assign probabilities of dSph membership to each observed star. Using this information, the velocity dispersions, central masses and central densities of the dark matter halos are calculated for these objects, and compared with the properties of the Milky Way dSph population. We also measure the average metallicity ([Fe/H]) from the co-added spectra of member stars for each M31 dSph and find that they are consistent with the trend of decreasing [Fe/H] with luminosity observed in the Milky Way population. We find that three of our studied M31 dSphs appear as significant outliers in terms of their central velocity dispersion, And XIX, XXI and XXV, all of which have large half-light radii (>700 pc) and low velocity dispersions (sigma_v<5 km/s). In addition, And XXV has a mass-to-light ratio within its half-light radius of just [M/L]_{half}=10.3^{+7.0}_{-6.7}, making it consistent with a simple stellar system with no appreciable dark matter component within its 1 sigma uncertainties. We suggest that the structure of the dark matter halos of these outliers have been significantly altered by tides.

preprint2013arXiv

A Vast Thin Plane of Co-rotating Dwarf Galaxies Orbiting the Andromeda Galaxy

Dwarf satellite galaxies are thought to be the remnants of the population of primordial structures that coalesced to form giant galaxies like the Milky Way. An early analysis noted that dwarf galaxies may not be isotropically distributed around our Galaxy, as several are correlated with streams of HI emission, and possibly form co-planar groups. These suspicions are supported by recent analyses, and it has been claimed that the apparently planar distribution of satellites is not predicted within standard cosmology, and cannot simply represent a memory of past coherent accretion. However, other studies dispute this conclusion. Here we report the existence (99.998% significance) of a planar sub-group of satellites in the Andromeda galaxy, comprising approximately 50% of the population. The structure is vast: at least 400 kpc in diameter, but also extremely thin, with a perpendicular scatter <14.1 kpc (99% confidence). Radial velocity measurements reveal that the satellites in this structure have the same sense of rotation about their host. This finding shows conclusively that substantial numbers of dwarf satellite galaxies share the same dynamical orbital properties and direction of angular momentum, a new insight for our understanding of the origin of these most dark matter dominated of galaxies. Intriguingly, the plane we identify is approximately aligned with the pole of the Milky Way&#39;s disk and is co-planar with the Milky Way to Andromeda position vector. The existence of such extensive coherent kinematic structures within the halos of massive galaxies is a fact that must be explained within the framework of galaxy formation and cosmology.

preprint2013arXiv

The large-scale structure of the halo of the Andromeda Galaxy Part I: global stellar density, morphology and metallicity properties

We present an analysis of the large-scale structure of the halo of the Andromeda galaxy, based on the Pan-Andromeda Archeological Survey (PAndAS), currently the most complete map of resolved stellar populations in any galactic halo. Despite copious substructure, the global halo populations follow closely power law profiles that become steeper with increasing metallicity. We divide the sample into stream-like populations and a smooth halo component. Fitting a three-dimensional halo model reveals that the most metal-poor populations ([Fe/H]<-1.7) are distributed approximately spherically (slightly prolate with ellipticity c/a=1.09+/-0.03), with only a relatively small fraction (42%) residing in discernible stream-like structures. The sphericity of the ancient smooth component strongly hints that the dark matter halo is also approximately spherical. More metal-rich populations contain higher fractions of stars in streams (86% for [Fe/H]>-0.6). The space density of the smooth metal-poor component has a global power-law slope of -3.08+/-0.07, and a non-parametric fit shows that the slope remains nearly constant from 30kpc to 300kpc. The total stellar mass in the halo at distances beyond 2 degrees is 1.1x10^10 Solar masses, while that of the smooth component is 3x10^9 Solar masses. Extrapolating into the inner galaxy, the total stellar mass of the smooth halo is plausibly 8x10^9 Solar masses. We detect a substantial metallicity gradient, which declines from [Fe/H]=-0.7 at R=30kpc to [Fe/H]=-1.5 at R=150kpc for the full sample, with the smooth halo being 0.2dex more metal poor than the full sample at each radius. While qualitatively in-line with expectations from cosmological simulations, these observations are of great importance as they provide a prototype template that such simulations must now be able to reproduce in quantitative detail.

preprint2013arXiv

The masses of Local Group dwarf spheroidal galaxies: The death of the universal mass profile

We investigate the claim that all dwarf spheroidal galaxies (dSphs) reside within halos that share a common, universal mass profile as has been derived for dSphs of the Galaxy. By folding in kinematic information for 25 Andromeda dSphs, more than doubling the previous sample size, we find that a singular mass profile can not be found to fit all the observations well. Further, the best-fit dark matter density profile measured for solely the Milky Way dSphs is marginally discrepant (at just beyond the 1 sigma level) with that of the Andromeda dSphs, where a profile with lower maximum circular velocity, and hence mass, is preferred. The agreement is significantly better when three extreme Andromeda outliers, And XIX, XXI and XXV, all of which have large half-light radii (>600pc) and low velocity dispersions (sigma_v < 5km/s) are omitted from the sample. We argue that the unusual properties of these outliers are likely caused by tidal interactions with the host galaxy.

preprint2013arXiv

The PAndAS view of the Andromeda satellite system - I. A Bayesian search for dwarf galaxies using spatial and color-magnitude information

We present a generic algorithm to search for dwarf galaxies in photometric catalogs and apply it to the Pan-Andromeda Archaeological Survey (PAndAS). The algorithm is developed in a Bayesian framework and, contrary to most dwarf-galaxy-search codes, makes use of both the spatial and color-magnitude information of sources in a probabilistic approach. Accounting for the significant contamination from the Milky Way foreground and from the structured stellar halo of the Andromeda galaxy, we recover all known dwarf galaxies in the PAndAS footprint with high significance, even for the least luminous ones. Some Andromeda globular clusters are also recovered and, in one case, discovered. We publish a list of the 143 most significant detections yielded by the algorithm. The combined properties of the 39 most significant isolated detections show hints that at least some of these trace genuine dwarf galaxies, too faint to be individually detected. Follow-up observations by the community are mandatory to establish which are real members of the Andromeda satellite system. The search technique presented here will be used in an upcoming contribution to determine the PAndAS completeness limits for dwarf galaxies. Although here tuned to the search of dwarf galaxies in the PAndAS data, the algorithm can easily be adapted to the search for any localised overdensity whose properties can be modeled reliably in the parameter space of any catalog.

preprint2013arXiv

The Three-Dimensional Structure of the M31 Satellite System; Strong Evidence for an Inhomogeneous Distribution of Satellites

We undertake an investigation into the spatial structure of the M31 satellite system utilizing the distance distributions presented in a previous publication. These distances make use of the unique combination of depth and spatial coverage of the Pan-Andromeda Archaeological Survey (PAndAS) to provide a large, homogeneous sample consisting of 27 of M31&#39;s satellites, as well as M31 itself. We find that the satellite distribution, when viewed as a whole, is no more planar than one would expect from a random distribution of equal size. A disk consisting of 15 of the satellites is however found to be highly significant, and strikingly thin, with a root-mean-square thickness of just $12.34^{+0.75}_{-0.43}$ kpc. This disk is oriented approximately edge on with respect to the Milky Way and almost perpendicular to the Milky Way disk. It is also roughly orthogonal to the disk like structure regularly reported for the Milky Way satellite system and in close alignment with M31&#39;s Giant Stellar Stream. A similar analysis of the asymmetry of the M31 satellite distribution finds that it is also significantly larger than one would expect from a random distribution. In particular, it is remarkable that 20 of the 27 satellites most likely lie on the Milky Way side of the galaxy, with the asymmetry being most pronounced within the satellite subset forming the aforementioned disk. This lopsidedness is all the more intriguing in light of the apparent orthogonality observed between the satellite disk structures of the Milky Way and M31.

preprint2012arXiv

A Deep, Wide-Field Study of Holmberg II with Suprime-Cam: Evidence for Ram Pressure Stripping

We present a deep, wide-field optical study of the M81 group dwarf galaxy Holmberg II (HoII) based on Subaru/Suprime-Cam imaging. Individual stars are resolved down to I~25.2, i.e. about 1.5 mag below the tip of the red giant branch (RGB). We use resolved star counts in the outskirts of the galaxy to measure the radial surface brightness profile down to μ_V~32 mag arcsec^-2, from which we determine a projected exponential scalelength of 0.70&#39;+-0.01&#39; (i.e. 0.69+-0.01 kpc). The composite profile, ranging from the cored centre out to R=7&#39;, is best fit by an EFF profile which gives a half-light radius of 1.41&#39;+-0.04&#39; (i.e. 1.39+-0.04 kpc), and an absolute magnitude M_V=-16.3. The low surface-brightness stellar component of HoII is regular and symmetric and has an extent much smaller than the vast HI cloud in which it is embedded. We compare the spatial distribution of the young, intermediate age, and old stellar populations, and find that the old RGB stars are significantly more centrally concentrated than the young stellar populations, contrary to what is observed in most dwarf galaxies of the Local Universe. We discuss these properties in the context of the comet-like distribution of HI gas around HoII, and argue for the presence of a hot intragroup medium in the vicinity of HoII to explain the contrasting morphologies of the gas and stars.

preprint2012arXiv

PAndAS in the mist: The stellar and gaseous mass within the halos of M31 and M33

Large scale surveys of the prominent members of the Local Group have provided compelling evidence for the hierarchical formation of massive galaxies, revealing a wealth of substructure that is thought to be the debris from ancient and on-going accretion events. In this paper, we compare two extant surveys of the M31-M33 subgroup of galaxies; the Pan-Andromeda Archaeological Survey (PAndAS) of the stellar structure, and a combination of observations of the HI gaseous content, detected at 21cm. Our key finding is a marked lack of spatial correlation between these two components on all scales, with only a few potential overlaps between stars and gas.The paucity of spatial correlation significantly restricts the analysis of kinematic correlations, although there does appear to the HI kinematically associated with the Giant Stellar Stream where it passes the disk of M31. These results demonstrate that that different processes must significantly influence the dynamical evolution of the stellar and HI components of substructures, such as ram pressure driving gas away from a purely gravitational path. Detailed modelling of the offset between the stellar and gaseous substructure will provide a determination of the properties of the gaseous halo of M31 and M33.

preprint2011arXiv

The Star Formation History and Dust Content in the Far Outer Disc of M31

We present a detailed analysis of two fields located 26 kpc (~5 scalelengths) from the centre of M31. One field samples the major axis populations--the Outer Disc field--while the other is offset by ~18&#39; and samples the Warp in the stellar disc. The CMDs based on HST/ACS imaging reach old main-sequence turn-offs (~12.5 Gyr). We apply the CMD-fitting technique to the Warp field to reconstruct the star formation history (SFH). We find that after undergoing roughly constant SF until about 4.5 Gyr ago, there was a rapid decline in activity and then a ~1.5 Gyr lull, followed by a strong burst lasting 1.5 Gyr and responsible for 25% of the total stellar mass in this field. This burst appears to be accompanied by a decline in metallicity which could be a signature of the inflow of metal-poor gas. The onset of the burst (~3 Gyr ago) corresponds to the last close passage of M31 and M33 as predicted by detailed N-body modelling, and may have been triggered by this event. We reprocess the deep M33 outer disc field data of Barker et al. (2011) in order to compare consistently-derived SFHs. This reveals a similar duration burst that is exactly coeval with that seen in the M31 Warp field, lending further support to the interaction hypothesis. The complex SFHs and the smoothly-varying age-metallicity relations suggest that the stellar populations observed in the far outer discs of both galaxies have largely formed in situ rather than migrated from smaller galactocentric radii. The strong differential reddening affecting the CMD of the Outer Disc field prevents derivation of the SFH. Instead, we quantify this reddening and find that the fine-scale distribution of dust precisely follows that of the HI gas. This indicates that the outer HI disc of M31 contains a substantial amount of dust and therefore suggests significant metal enrichment in these parts, consistent with inferences from our CMD analysis.

preprint2010arXiv

Modelling of the Complex CASSOWARY/SLUGS Gravitational Lenses

We present the first high-resolution images of CSWA 31, a gravitational lens system observed as part of the SLUGS (Sloan Lenses Unravelled by Gemini Studies) program. These systems exhibit complex image structure with the potential to strongly constrain the mass distribution of the massive lens galaxies, as well as the complex morphology of the sources. In this paper, we describe the strategy used to reconstruct the unlensed source profile and the lens galaxy mass profiles. We introduce a prior distribution over multi-wavelength sources that is realistic as a representation of our knowledge about the surface brightness profiles of galaxies and groups of galaxies. To carry out the inference computationally, we use Diffusive Nested Sampling, an efficient variant of Nested Sampling that uses Markov Chain Monte Carlo (MCMC) to sample the complex posterior distributions and compute the normalising constant. We demonstrate the efficacy of this approach with the reconstruction of the group-group gravitational lens system CSWA 31, finding the source to be composed of five merging spiral galaxies magnified by a factor of 13.

preprint2010arXiv

The photometric properties of a vast stellar substructure in the outskirts of M33

We have surveyed $\sim40$sq.degrees surrounding M33 with CFHT MegaCam in the g and i filters, as part of the Pan-Andromeda Archaeological Survey. Our observations are deep enough to resolve the top 4mags of the red giant branch population in this galaxy. We have previously shown that the disk of M33 is surrounded by a large, irregular, low-surface brightness substructure. Here, we quantify the stellar populations and structure of this feature using the PAndAS data. We show that the stellar populations of this feature are consistent with an old population with $<[Fe/H]>\sim-1.6$dex and an interquartile range in metallicity of $\sim0.5$dex. We construct a surface brightness map of M33 that traces this feature to $μ_V\simeq33$mags\,arcsec$^{-2}$. At these low surface brightness levels, the structure extends to projected radii of $\sim40$kpc from the center of M33 in both the north-west and south-east quadrants of the galaxy. Overall, the structure has an &#34;S-shaped&#34; appearance that broadly aligns with the orientation of the HI disk warp. We calculate a lower limit to the integrated luminosity of the structure of $-12.7\pm0.5$mags, comparable to a bright dwarf galaxy such as Fornax or AndII and slightly less than $1\$ of the total luminosity of M33. Further, we show that there is tentative evidence for a distortion in the distribution of young stars near the edge of the HI disk that occurs at similar azimuth to the warp in HI. The data also hint at a low-level, extended stellar component at larger radius that may be a M33 halo component. We revisit studies of M33 and its stellar populations in light of these new results, and we discuss possible formation scenarios for the vast stellar structure. Our favored model is that of the tidal disruption of M33 in its orbit around M31.

preprint2009arXiv

IPHAS A-type Stars with Mid-IR Excesses in Spitzer Surveys

We have identified 17 A-type stars in the Galactic Plane that have mid-IR excesses at 8 micron. From the observed colors in the (r&#39;-H_alpha)-(r&#39;-i&#39;) plane, we first identified 23050 early A-type main sequence (MS) star candidates in the Isaac Newton Photometric H-Alpha Survey (IPHAS) point source database that are located in Spitzer GLIMPSE Galactic Plane fields. Imposing the requirement that they be detected in all seven 2MASS and IRAC bands led to a sample of 2692 candidate A-type stars with fully sampled 0.6 to 8 micron SEDs. Optical classification spectra of 18 of the IPHAS candidate A-type MS stars showed that all but one could be well fitted using main sequence A-type templates, with the other being an A-type supergiant. Out of the 2692 A-type candidates 17 (0.6%) were found to have 8-micron excesses above the expected photospheric values. Taking into account non-A-Type contamination estimates, the 8-micron excess fraction is adjusted to ~0.7%. The distances to these sources range from 0.7-2.5 kpc. Only 10 out of the 17 excess stars had been covered by Spitzer MIPSGAL survey fields, of which 5 had detectable excesses at 24 micron. For sources with excesses detected in at least two mid-IR wavelength bands, blackbody fits to the excess SEDs yielded temperatures ranging from 270 to 650 K, and bolometric luminosity ratios L_IR/L* from 2.2x10^{-3}-1.9x10^{-2}, with a mean value of 7.9x10^{-3} (these bolometric luminosities are lower limits as cold dust is not detectable by this survey). Both the presence of mid-IR excesses and the derived bolometric luminosity ratios are consistent with many of these systems being in the planet-building transition phase between the early protoplanetary disk phase and the later debris disk phase.

preprint2009arXiv

The remnants of galaxy formation from a panoramic survey of the region around M31

In hierarchical cosmological models, galaxies grow in mass through the continual accretion of smaller ones. The tidal disruption of these systems is expected to result in loosely bound stars surrounding the galaxy, at distances that reach $10 - 100$ times the radius of the central disk. The number, luminosity and morphology of the relics of this process provide significant clues to galaxy formation history, but obtaining a comprehensive survey of these components is difficult because of their intrinsic faintness and vast extent. Here we report a panoramic survey of the Andromeda galaxy (M31). We detect stars and coherent structures that are almost certainly remnants of dwarf galaxies destroyed by the tidal field of M31. An improved census of their surviving counterparts implies that three-quarters of M31&#39;s satellites brighter than $M_V < -6$ await discovery. The brightest companion, Triangulum (M33), is surrounded by a stellar structure that provides persuasive evidence for a recent encounter with M31. This panorama of galaxy structure directly confirms the basic tenets of the hierarchical galaxy formation model and reveals the shared history of M31 and M33 in the unceasing build-up of galaxies.