As smaller analogs of Active Galactic Nuclei (AGN), X-ray Binaries (XRBs) are also capable of launching jets that accelerate particles to high energies. In this work, we reexamine XRB jets as potential sources of high-energy cosmic rays (CRs) and explore whether they could provide a significant second Galactic component to the CR spectrum. In the most intriguing scenario, XRB-CRs could dominate the observed spectrum above the so-called "knee" feature at $\sim 3 \times 10^{15}$ eV, offering an explanation for several key issues in this transition zone from Galactic to extragalactic CRs. We discuss how such a scenario could be probed in the near future via multi-messenger observations of XRB jets, as well as diffuse Galactic neutrino flux measurements.

We explore the predictive power of cosmological, hydrodynamical simulations for stellar phase space substructure and velocity correlations with the Auriga simulations and Aurigaia mock-Gaia catalogues. We show that at the solar circle the Auriga simulations commonly host phase space structures in the stellar component that have constant orbital energies and arise from accreted subhaloes. These structures can persist for a few Gyrs, even after coherent streams in position space have been erased. We also explore velocity two-point correlation functions and find this diagnostic is not deterministic for particular clustering patterns in phase space. Finally, we explore these structure diagnostics with the Aurigaia catalogues and show that the current catalogues have the ability to recover some structures in phase space but careful consideration is required to separate physical structures from numerical structures arising from catalogue generation methods.

We present the status of the Dark Energy Spectroscopic Instrument (DESI) and its plans and opportunities for the coming decade. DESI construction and its initial five years of operations are an approved experiment of the US Department of Energy and is summarized here as context for the Astro2020 panel. Beyond 2025, DESI will require new funding to continue operations. We expect that DESI will remain one of the world's best facilities for wide-field spectroscopy throughout the decade. More about the DESI instrument and survey can be found at https://www.desi.lbl.gov.

We investigate how strong gravitational lensing can test contemporary models of massive elliptical (ME) galaxy formation, by combining a traditional decomposition of their visible stellar distribution with a lensing analysis of their mass distribution. As a proof of concept, we study a sample of three ME lenses, observing that all are composed of two distinct baryonic structures, a `red' central bulge surrounded by an extended envelope of stellar material. Whilst these two components look photometrically similar, their distinct lensing effects permit a clean decomposition of their mass structure. This allows us to infer two key pieces of information about each lens galaxy: (i) the stellar mass distribution (without invoking stellar populations models) and (ii) the inner dark matter halo mass. We argue that these two measurements are crucial to testing models of ME formation, as the stellar mass profile provides a diagnostic of baryonic accretion and feedback whilst the dark matter mass places each galaxy in the context of LCDM large scale structure formation. We also detect large rotational offsets between the two stellar components and a lopsidedness in their outer mass distributions, which hold further information on the evolution of each ME. Finally, we discuss how this approach can be extended to galaxies of all Hubble types and what implication our results have for studies of strong gravitational lensing.

The Great Andromeda Galaxy (M31) is the nexus of the near-far galaxy evolution connection and a principal data point for near-field cosmology. Due to its proximity (780 kpc), M31 can be resolved into individual stars like the Milky Way (MW). Unlike the MW, we have the advantage of a global view of M31, enabling M31 to be observed with techniques that also apply to more distant galaxies. Moreover, recent evidence suggests that M31 may have survived a major merger within the last several Gyr, shaping the morphology of its stellar halo and triggering a starburst, while leaving the stellar disk largely intact. The MW and M31 thus provide complementary opportunities for in-depth studies of the disks, halos, and satellites of L* galaxies. Our understanding of the M31 system will be transformed in the 2020s if they include wide field facilities for both photometry (HST-like sensitivity and resolution) and spectroscopy (10-m class telescope, >1 sq. deg. field, highly multiplexed, R~ 3000 to 6000). We focus here on the power of these facilities to constrain the past, present, and future merger history of M31, via chemo-dynamical analyses and star formation histories of phase-mixed stars accreted at early times, as well as stars in surviving tidal debris features, M31's extended disk, and intact satellite galaxies that will eventually be tidally incorporated into the halo. This will yield an unprecedented view of the hierarchical formation of the M31 system and the subhalos that built it into the L* galaxy we observe today.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

We present and analyse mock stellar catalogues that match the selection criteria and observables (including uncertainties) of the Gaia satellite data release 2 (DR2). The source are six cosmological high-resolution magneto-hydrodynamic $\Lambda$CDM zoom simulations of the formation of Milky Way analogues from the AURIGA project. Mock data are provided for stars with $V < 16$ mag, and $V < 20$ mag at $|b|>20$ degrees. The mock catalogues are made using two different methods: the public SNAPDRAGONS code, and a method based on that of Lowing et al. that preserves the phase-space distribution of the model stars. These publicly available catalogues contain 5-parameter astrometry, radial velocities, multi-band photometry, stellar parameters, dust extinction values, and uncertainties in all these quantities. In addition, we provide the gravitational potential and information on the origin of each star. By way of demonstration, we apply the mock catalogues to analyses of the young stellar disc and the stellar halo. We show that: i) the young outer stellar disc exhibits a flared distribution that is detectable in the height and vertical velocity distribution of A- and B-dwarf stars up to radii of ~15 kpc; and ii) the spin of the stellar halo out to 100 kpc can be accurately measured with Gaia DR2 RR Lyrae stars. These catalogues are well suited for comparisons with observations and should help to: i) develop and test analysis methods for the Gaia DR2 data; ii) gauge the limitations and biases of the data and iii) interpret the data in the light of theoretical predictions from realistic $ab$ $initio$ simulations of galaxy formation in the $\Lambda$CDM cosmological model.

We present results from a pilot observation of nearby (~20 Mpc) galaxies with mass similar to that of the Milky Way (MW) to address the missing satellite problem. This is the first paper from an on-going project to address the problem with a statistical sample of galaxies outside of the Local Group (LG) without employing an assumption that the LG is a typical halo in the Universe. Thanks to the close distances of our targets, dwarf galaxies around them can be identified as extended, diffuse galaxies. By applying a surface brightness cut together with a careful visual screening to remove artifacts and background contamination, we construct a sample of dwarf galaxies. The luminosity function (LF) of one of the targets is broadly consistent with that of the MW, but the other has a more abundant dwarf population. Numerical simulations by Okamoto (2013) seem to overpredict the number of dwarfs on average, while more recent predictions from Copernicus Complexio are in a better agreement. In both observations and simulations, there is a large diversity in the LFs, demonstrating the importance of addressing the missing satellite problem with a statistically representative sample. We also characterize the projected spatial distributions of the satellites and do not observe strong evidence for alignments around the central galaxies. Based on this successful pilot observation, we are carrying out further observations to increase the sample of nearby galaxies, which we plan to report in our future paper.

In hierarchical models of galaxy formation, stellar tidal streams are expected around most, if not all, galaxies. Although these features may provide useful diagnostics of the $\Lambda$CDM model, their observational properties remain poorly constrained because they are challenging to detect and interpret and have been studied in detail for only a sparse sampling of galaxy population. More quantitative, systematic approaches are required. We advocate statistical analysis of the counts and properties of such features in archival wide-field imaging surveys for a direct comparison against results from numerical simulations. Thus, in this paper we aim to study systematically the frequency of occurrence and other observational properties of tidal features around nearby galaxies. The sample we construct will act as a foundational dataset for statistical comparison with cosmological models of galaxy formation. Our approach is based on a visual classification of diffuse features around a volume-limited sample of nearby galaxies, using a post-processing of Sloan Digital Sky Survey imaging optimized for the detection of stellar structure with low surface brightness. At a limiting surface brightness of $28\ \mathrm{mag~arcsec^{-2}}$, 14% of the galaxies in our sample exhibit evidence of diffuse features likely to have arisen from minor merging events. Our technique recovers all previously known streams in our sample and yields a number of new candidates. Consistent with previous studies, coherent arc-like features and shells are the most common type of tidal structures found in this study. We conclude that although some detections are ambiguous and could be corroborated or refuted with deeper imaging, our technique provides a reliable foundation for the statistical analysis of diffuse circumgalactic features in wide-area imaging surveys, and for the identification of targets for follow-up studies.

We present the extended GALEX Arecibo SDSS Survey (xGASS), a gas fraction-limited census of the atomic (HI) gas content of 1179 galaxies selected only by stellar mass ($M_\star =10^{9}-10^{11.5} M_\odot$) and redshift ($0.01<z<0.05$). This includes new Arecibo observations of 208 galaxies, for which we release catalogs and HI spectra. In addition to extending the GASS HI scaling relations by one decade in stellar mass, we quantify total (atomic+molecular) cold gas fractions and molecular-to-atomic gas mass ratios, $R_{mol}$, for the subset of 477 galaxies observed with the IRAM 30 m telescope. We find that atomic gas fractions keep increasing with decreasing stellar mass, with no sign of a plateau down to $\log M_\star/M_\odot = 9$. Total gas reservoirs remain HI-dominated across our full stellar mass range, hence total gas fraction scaling relations closely resemble atomic ones, but with a scatter that strongly correlates with $R_{mol}$, especially at fixed specific star formation rate. On average, $R_{mol}$ weakly increases with stellar mass and stellar surface density $\mu_\star$, but individual values vary by almost two orders of magnitude at fixed $M_\star$ or $\mu_\star$. We show that, for galaxies on the star-forming sequence, variations of $R_{mol}$ are mostly driven by changes of the HI reservoirs, with a clear dependence on $\mu_\star$. Establishing if galaxy mass or structure plays the most important role in regulating the cold gas content of galaxies requires an accurate separation of bulge and disk components for the study of gas scaling relations.

Particle tagging is an efficient, but approximate, technique for using cosmological N-body simulations to model the phase-space evolution of the stellar populations predicted, for example, by a semi-analytic model of galaxy formation. We test the technique developed by Cooper et al. (which we call STINGS here) by comparing particle tags with stars in a smooth particle hydrodynamic (SPH) simulation. We focus on the spherically averaged density profile of stars accreted from satellite galaxies in a Milky Way (MW)-like system. The stellar profile in the SPH simulation can be recovered accurately by tagging dark matter (DM) particles in the same simulation according to a prescription based on the rank order of particle binding energy. Applying the same prescription to an N-body version of this simulation produces a density profile differing from that of the SPH simulation by <10 per cent on average between 1 and 200 kpc. This confirms that particle tagging can provide a faithful and robust approximation to a self-consistent hydrodynamical simulation in this regime (in contradiction to previous claims in the literature). We find only one systematic effect, likely due to the collisionless approximation, namely that massive satellites in the SPH simulation are disrupted somewhat earlier than their collisionless counterparts. In most cases this makes remarkably little difference to the spherically averaged distribution of their stellar debris. We conclude that, for galaxy formation models that do not predict strong baryonic effects on the present-day DM distribution of MW-like galaxies or their satellites, differences in stellar halo predictions associated with the treatment of star formation and feedback are much more important than those associated with the dynamical limitations of collisionless particle tagging.

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

Current data broadly support trends of galaxy surface brightness profile amplitude and shape with total stellar mass predicted by state-of-the-art Lambda-CDM cosmological simulations, although recent results show signs of interesting discrepancies, particularly for galaxies less massive than the Milky Way. Here I discuss how perhaps the largest contribution to such discrepancies can be inferred almost directly from how well a given model agrees with the observed present-day galaxy stellar mass function.

We use a pair of high resolution N-body simulations implementing two dark matter models, namely the standard cold dark matter (CDM) cosmogony and a warm dark matter (WDM) alternative where the dark matter particle is a 1.5keV thermal relic. We combine these simulations with the GALFORM semi-analytical galaxy formation model in order to explore differences between the resulting galaxy populations. We use GALFORM model variants for CDM and WDM that result in the same z=0 galaxy stellar mass function by construction. We find that most of the studied galaxy properties have the same values in these two models, indicating that both dark matter scenarios match current observational data equally well. Even in under-dense regions, where discrepancies in structure formation between CDM and WDM are expected to be most pronounced, the galaxy properties are only slightly different. The only significant difference in the local universe we find is in the galaxy populations of "Local Volumes", regions of radius 1 to 8Mpc around simulated Milky Way analogues. In such regions our WDM model provides a better match to observed local galaxy number counts and is five times more likely than the CDM model to predict sub-regions within them that are as empty as the observed Local Void. Thus, a highly complete census of the Local Volume and future surveys of void regions could provide constraints on the nature of dark matter.

Observations of diffuse starlight in the outskirts of galaxies are thought to be a fundamental source of constraints on the cosmological context of galaxy assembly in the $\Lambda$CDM model. Such observations are not trivial because of the extreme faintness of such regions. In this work, we investigate the photometric properties of six massive early type galaxies (ETGs) in the VEGAS sample (NGC 1399, NGC 3923, NGC 4365, NGC 4472, NGC 5044, and NGC 5846) out to extremely low surface brightness levels, with the goal of characterizing the global structure of their light profiles for comparison to state-of-the-art galaxy formation models. We carry out deep and detailed photometric mapping of our ETG sample taking advantage of deep imaging with VST/OmegaCAM in the g and i bands. By fitting the light profiles, and comparing the results to simulations of elliptical galaxy assembly, we identify signatures of a transition between "relaxed" and "unrelaxed" accreted components and can constrain the balance between in situ and accreted stars. The very good agreement of our results with predictions from theoretical simulations demonstrates that the full VEGAS sample of $\sim 100$ ETGs will allow us to use the distribution of diffuse light as a robust statistical probe of the hierarchical assembly of massive galaxies.

Using the exquisite depth of the Hubble Ultra Deep Field (HUDF12 programme) dataset, we explore the ongoing assembly of the outermost regions of the most massive galaxies ($\rm M_{\rm stellar}\geq$ 5$\times$10$^{10}$ M$_{\odot}$) at $z \leq$ 1. The outskirts of massive objects, particularly Early-Types Galaxies (ETGs), are expected to suffer a dramatic transformation across cosmic time due to continuous accretion of small galaxies. HUDF imaging allows us to study this process at intermediate redshifts in 6 massive galaxies, exploring the individual surface brightness profiles out to $\sim$25 effective radii. We find that 5-20\% of the total stellar mass for the galaxies in our sample is contained within 10 $< R <$ 50 kpc. These values are in close agreement with numerical simulations, and higher than those reported for local late-type galaxies ($\lesssim$5\%). The fraction of stellar mass stored in the outer envelopes/haloes of Massive Early-Type Galaxies increases with decreasing redshift, being 28.7\% at $< z > =$ 0.1, 15.1\% at $< z > =$ 0.65 and 3.5\% at $< z > =$ 2. The fraction of mass in diffuse features linked with ongoing minor merger events is $>$ 1-2\%, very similar to predictions based on observed close pair counts. Therefore, the results for our small albeit meaningful sample suggest that the size and mass growth of the most massive galaxies have been solely driven by minor and major merging from $z =$ 1 to today.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

We use the APOSTLE cosmological hydrodynamic simulations to examine the effects of tidal stripping on cold dark matter (CDM) sub haloes that host three of the most luminous Milky Way (MW) dwarf satellite galaxies: Fornax, Sculptor, and Leo I. We identify simulated satellites that match the observed spatial and kinematic distributions of stars in these galaxies, and track their evolution after infall. We find $\sim$ 30$\%$ of subhaloes hosting satellites with present-day stellar mass $10^6$-$10^8$ $M_{\odot}$ experience $>20\%$ stellar mass loss after infall. Fornax analogues have earlier infall times compared to Sculptor and Leo I analogues. Star formation in Fornax analogues continues for $\sim3$-$6$ Gyr after infall, whereas Sculptor and Leo I analogues stop forming stars $< 2$-$3$ Gyr after infall. Fornax analogues typically show more significant stellar mass loss and exhibit stellar tidal tails, whereas Sculptor and Leo I analogues, which are more deeply embedded in their host DM haloes at infall, do not show substantial mass loss due to tides. When additionally comparing the orbital motion of the host subaloes to the measured proper motion of Fornax we find the matching more difficult; host subhaloes tend to have pericentres smaller than that measured for Fornax itself. From the kinematic and orbital data, we estimate that Fornax has lost $10-20\%$ of its infall stellar mass. Our best estimate for the surface brightness of a stellar tidal stream associated with Fornax is $\Sigma \sim$ 32.6 mag $ {\rm arcsec^{-2}}$, which may be detectable with deep imaging surveys such as DES and LSST.

We use the Wide-field Infrared Survey Explorer (WISE) and the Sloan Digital Sky Survey (SDSS) to confirm a connection between dust-obscured active galactic nuclei (AGNs) and galaxy merging. Using a new, volume-limited (z <= 0.08) catalog of visually-selected major mergers and galaxy-galaxy interactions from the SDSS, with stellar masses above 2x10^10 Msun, we find that major mergers (interactions) are 5-17 (3-5) times more likely to have red [3.4]-[4.6] colors associated with dust-obscured or `dusty' AGNs, compared to non-merging galaxies with similar masses. Using published fiber spectral diagnostics, we map the [3.4]-[4.6] versus [4.6]-[12] colors of different emission-line galaxies and find one-quarter of Seyferts have colors indicative of a dusty AGN. We find that AGNs are five times more likely to be obscured when hosted by a merging galaxy, half of AGNs hosted by a merger are dusty, and we find no enhanced frequency of optical AGNs in merging over non-merging galaxies. We conclude that undetected AGNs missed at shorter wavelengths are at the heart of the ongoing AGN-merger connection debate. The vast majority of mergers hosting dusty AGNs are star-forming and located at the centers of Mhalo<10^13 Msun groups. Assuming plausibly short duration dusty-AGN phases, we speculate that a large fraction of gas-rich mergers experience a brief obscured AGN phase, in agreement with the strong connection between central star formation and black hole growth seen in merger simulations.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

The mass of the dark matter halo of the Milky Way can be estimated by fitting analytical models to the phase-space distribution of dynamical tracers. We test this approach using realistic mock stellar halos constructed from the Aquarius N-body simulations of dark matter halos in the $\Lambda$CDM cosmology. We extend the standard treatment to include a Navarro-Frenk-White (NFW) potential and use a maximum likelihood method to recover the parameters describing the simulated halos from the positions and velocities of their mock halo stars. We find that the estimate of halo mass is highly correlated with the estimate of halo concentration. The best-fit halo masses within the virial radius, $R_{200}$, are biased, ranging from a 40\% underestimate to a 5\% overestimate in the best case (when the tangential velocities of the tracers are included). There are several sources of bias. Deviations from dynamical equilibrium can potentially cause significant bias; deviations from spherical symmetry are relatively less important. Fits to stars at different galactocentric radii can give different mass estimates. By contrast, the model gives good constraints on the mass within the half-mass radius of tracers even when restricted to tracers within 60kpc. The recovered velocity anisotropies of tracers, $\beta$, are biased systematically, but this does not affect other parameters if tangential velocity data are used as constraints.

We study the formation of stellar haloes in three Milky Way-mass galaxies using cosmological smoothed particle hydrodynamics simulations, focusing on the subset of halo stars that form in situ, as opposed to those accreted from satellites. In situ stars in our simulations dominate the stellar halo out to 20 kpc and account for 30-40 per cent of its total mass. We separate in situ halo stars into three straightforward, physically distinct categories according to their origin: stars scattered from the disc of the main galaxy ('heated disc'), stars formed from gas smoothly accreted on to the halo ('smooth' gas) and stars formed in streams of gas stripped from infalling satellites ('stripped' gas). We find that most belong to the stripped gas category. Those originating in smooth gas outside the disc tend to form at the same time and place as the stripped-gas population, suggesting that their formation is associated with the same gas-rich accretion events. The scattered disc star contribution is negligible overall but significant in the Solar neighbourhood, where ~90 per cent of stars on eccentric orbits once belonged to the disc. However, the distinction between halo and thick disc in this region is highly ambiguous. The chemical and kinematic properties of the different components are very similar at the present day, but the global properties of the in situ halo differ substantially between the three galaxies in our study. In our simulations, the hierarchical buildup of structure is the driving force behind not only the accreted stellar halo, but also those halo stars formed in situ.

We show that diffusion due to chaotic mixing in the Neighbourhood of the Sun may not be as relevant as previously suggested in erasing phase space signatures of past Galactic accretion events. For this purpose, we analyse Solar Neighbourhood-like volumes extracted from cosmological simulations that naturally account for chaotic orbital behaviour induced by the strongly triaxial and cuspy shape of the resulting dark matter haloes, among other factors. In the approximation of an analytical static triaxial model, our results show that a large fraction of stellar halo particles in such local volumes have chaos onset times (i.e., the timescale at which stars commonly associated with chaotic orbits will exhibit their chaotic behaviour) significantly larger than a Hubble time. Furthermore, particles that do present a chaotic behaviour within a Hubble time do not exhibit significant diffusion in phase space.

We simulate the phase-space distribution of stellar mass in 9 massive Lambda-CDM galaxy clusters by applying the semi-analytic particle tagging method of Cooper et al. to the Phoenix suite of high-resolution N-body simulations (M200 = 7.5 to 33 x 10^14 Msol). The resulting surface brightness (SB) profiles of brightest cluster galaxies (BCGs) match well to observations. On average, stars formed in galaxies accreted by the BCG account for ~90 per cent of its total mass (the remainder is formed in situ). In circular BCG-centred apertures, the superposition of multiple debris clouds (each ~10 per cent of the total BCG mass) from different progenitors can result in an extensive outer diffuse component, qualitatively similar to a 'cD envelope'. These clouds typically originate from tidal stripping at z < 1 and comprise both streams and the extended envelopes of other massive galaxies in the cluster. Stars at very low SB contribute a significant fraction of the total cluster stellar mass budget: in the central 1 Mpc^2 of a z ~ 0.15 cluster imaged at SDSS-like resolution, our fiducial model predicts 80-95 per cent of stellar mass below a SB of mu_V = 26.5 mag arcsec^2 is associated with accreted stars in the envelope of the BCG. The ratio of BCG stellar mass (including this diffuse component) to total cluster stellar mass is ~30 per cent.

We apply a semi-analytic galaxy formation model to two high resolution cosmological N-body simulations to investigate analogues of the Milky Way system. We select these according to observed properties of the Milky Way rather than by halo mass as in most previous work. For disk-dominated central galaxies with stellar mass (5--7) x 10d10Msun, the median host halo mass is 1.4 x 10d12Msun, with 1 sigma dispersion in the range [0.86, 3.1] x 10d12Msun, consistent with dynamical measurements of the Milky Way halo mass. For any given halo mass, the probability of hosting a Milky Way system is low, with a maximum of ~20% in haloes of mass ~10d12Msun. The model reproduces the V-band luminosity function and radial profile of the bright (MV < -9) Milky Way satellites. Galaxy formation in low mass haloes is found to be highly stochastic, resulting in an extremely large scatter in the relation between MV (or stellar mass) for satellites and the depth of the subhalo potential well in which they live, as measured by the maximum of the rotation curve, Vmax. We conclude that the "too big to fail" problem is an artifact of selecting satellites in N-body simulations according to subhalo properties: in 10% of cases we find that three or fewer of the brightest (or most massive) satellites have Vmax > 30 km/s. Our model predicts that around half of the dark matter subhaloes with Vmax > 20 km/s host satellites fainter than MV = -9 and so may be missing from existing surveys.

We establish a controlled comparison between the properties of galactic stellar halos obtained with hydrodynamical simulations and with `particle tagging'. Tagging is a fast way to obtain stellar population dynamics: instead of tracking gas and star formation, it `paints' stars directly onto a suitably defined subset of dark matter particles in a collisionless, dark-matter-only simulation.Our study shows that there are conditions under which particle tagging generates good fits to the hydrodynamical stellar density profiles of a central Milky-Way-like galaxy and its most prominent substructure. Phase-space diffusion processes are crucial to reshaping the distribution of stars in infalling spheroidal systems and hence the final stellar halo. We conclude that the success of any particular tagging scheme hinges on this diffusion being taken into account, at a minimum by making use of `live' tagging schemes, in which particles are regularly tagged throughout the evolution of a galaxy.

We present a new technique for creating mock catalogues of the individual stars that make up the accreted component of stellar haloes in cosmological simulations and show how the catalogues can be used to test and interpret observational data. The catalogues are constructed from a combination of methods. A semi-analytic galaxy formation model is used to calculate the star formation history in haloes in an N-body simulation and dark matter particles are tagged with this stellar mass. The tags are converted into individual stars using a stellar population synthesis model to obtain the number density and evolutionary stage of the stars, together with a phase-space sampling method that distributes the stars while ensuring that the phase-space structure of the original N-body simulation is maintained. A set of catalogues based on the $\Lambda$CDM Aquarius simulations of Milky Way mass haloes have been created and made publicly available on a website. Two example applications are discussed that demonstrate the power and flexibility of the mock catalogues. We show how the rich stellar substructure that survives in the stellar halo precludes a simple measurement of its density profile and demonstrate explicitly how pencil-beam surveys can return almost any value for the slope of the profile. We also show that localized variations in the abundance of particular types of stars, a signature of differences in the composition of stellar populations, allow streams to be easily identified.

Recent work has suggested that the amplitude of the size mass relation of massive early type galaxies evolves with redshift. Here we use a semi-analytical galaxy formation model to study the size evolution of massive early type galaxies. We find this model is able to reproduce the amplitude of present day amplitude and slope of the relation between size and stellar mass for these galaxies, as well as its evolution. The amplitude of this relation reflects the typical compactness of dark halos at the time when most of the stars are formed. This link between size and star formation epoch is propagated in galaxy mergers. Mergers of high or moderate mass ratio (less than 1:3) become increasingly important with increasing present day stellar mass for galaxies more massive than $10^{11.4}M_{\odot}$. At lower masses, low mass ratio mergers play a more important role. In situ star formation contribute more to the size growth than it does to stellar mass growth. We also find that, for ETGs identified at $z=2$, minor mergers dominate subsequent growth both for stellar mass and in size, consistent with earlier theoretical results.

We use the Sloan Digital Sky Survey to investigate the properties of massive elliptical galaxies in the local Universe (z\leq0.08) that have unusually blue optical colors. Through careful inspection, we distinguish elliptical from non-elliptical morphologies among a large sample of similarly blue galaxies with high central light concentrations (c_r\geq2.6). These blue ellipticals comprise 3.7 per cent of all c_r\geq2.6 galaxies with stellar masses between 10^10 and 10^11 h^{-2} {\rm M}_{\sun}. Using published fiber spectra diagnostics, we identify a unique subset of 172 non-star-forming ellipticals with distinctly blue urz colors and young (< 3 Gyr) light-weighted stellar ages. These recently quenched ellipticals (RQEs) have a number density of 2.7-4.7\times 10^{-5}\,h^3\,{\rm Mpc}^{-3} and sufficient numbers above 2.5\times10^{10} h^{-2} {\rm M}_{\sun} to account for more than half of the expected quiescent growth at late cosmic time assuming this phase lasts 0.5 Gyr. RQEs have properties that are consistent with a recent merger origin (i.e., they are strong `first-generation' elliptical candidates), yet few involved a starburst strong enough to produce an E+A signature. The preferred environment of RQEs (90 per cent reside at the centers of < 3\times 10^{12}\,h^{-1}{\rm M}_{\sun} groups) agrees well with the `small group scale' predicted for maximally efficient spiral merging onto their halo center and rules out satellite-specific quenching processes. The high incidence of Seyfert and LINER activity in RQEs and their plausible descendents may heat the atmospheres of small host halos sufficiently to maintain quenching.

We use the very high resolution, fully cosmological simulations from the Aquarius project, coupled to a semi-analytical model of galaxy formation, to study the phase-space distribution of halo stars in "solar neighbourhood"-like volumes. We find that this distribution is very rich in substructure in the form of stellar streams for all five stellar haloes we have analysed. These streams can be easily identified in velocity space, as well as in spaces of pseudo-conserved quantities such as E vs. Lz. In our best-resolved local volumes, the number of identified streams ranges from ~ 300 to 600, in very good agreement with previous analytical predictions, even in the presence of chaotic mixing. The fraction of particles linked to (massive) stellar streams in these volumes can be as large as 84%. The number of identified streams is found to decrease as a power-law with galactocentric radius. We show that the strongest limitation to the quantification of substructure in our poorest-resolved local volumes is particle resolution rather than strong diffusion due to chaotic mixing.

We present the final data release from the GALEX Arecibo SDSS Survey (GASS), a large Arecibo program that measured the HI properties for an unbiased sample of ~800 galaxies with stellar masses greater than 10^10 Msun and redshifts 0.025<z<0.05. This release includes new Arecibo observations for 250 galaxies. We use the full GASS sample to investigate environmental effects on the cold gas content of massive galaxies at fixed stellar mass. The environment is characterized in terms of dark matter halo mass, obtained by cross-matching our sample with the SDSS group catalog of Yang et al. Our analysis provides, for the first time, clear statistical evidence that massive galaxies located in halos with masses of 10^13-10^14 Msun have at least 0.4 dex less HI than objects in lower density environments. The process responsible for the suppression of gas in group galaxies most likely drives the observed quenching of the star formation in these systems. Our findings strongly support the importance of the group environment for galaxy evolution, and have profound implications for semi-analytic models of galaxy formation, which currently do not allow for stripping of the cold interstellar medium in galaxy groups.

We have combined the semi-analytic galaxy formation model of Guo et al. (2011) with the particle-tagging technique of Cooper et al. (2010) to predict galaxy surface brightness profiles in a representative sample of ~1900 massive dark matter haloes (10^12--10^14 M_sol) from the Millennium II Lambda-CDM N-body simulation. Here we present our method and basic results focusing on the outer regions of galaxies, consisting of stars accreted in mergers. These simulations cover scales from the stellar haloes of Milky Way-like galaxies to the 'cD envelopes' of groups and clusters, and resolve low surface brightness substructure such as tidal streams. We find that the surface density of accreted stellar mass around the central galaxies of dark matter haloes is well described by a Sersic profile, the radial scale and amplitude of which vary systematically with halo mass (M_200). The total stellar mass surface density profile breaks at the radius where accreted stars start to dominate over stars formed in the galaxy itself. This break disappears with increasing M_200 because accreted stars contribute more of the total mass of galaxies, and is less distinct when the same galaxies are averaged in bins of stellar mass, because of scatter in the relation between M_star and M_200. To test our model we have derived average stellar mass surface density profiles for massive galaxies at z~0.08 by stacking SDSS images. Our model agrees well with these stacked profiles and with other data from the literature, and makes predictions that can be more rigorously tested by future surveys that extend the analysis of the outer structure of galaxies to fainter isophotes. We conclude that it is likely that the outer structure of the spheroidal components of galaxies is largely determined by collisionless merging during their hierarchical assembly

ASTRA (ASTrometric and phase-Referencing Astronomy) is an upgrade to the existing Keck Interferometer which aims at providing new self-phase referencing (high spectral resolution observation of YSOs), dual-field phase referencing (sensitive AGN observations), and astrometric (known exoplanetary systems characterization and galactic center general relativity in strong field regime) capabilities. With the first high spectral resolution mode now offered to the community, this contribution focuses on the progress of the dual field and astrometric modes.

Ground-based long baseline interferometers have long been limited in sensitivity by the short integration periods imposed by atmospheric turbulence. The first observation fainter than this limit was performed on January 22, 2011 when the Keck Interferometer observed a K=11.5 target, about one magnitude fainter than its K=10.3 limit. This observation was made possible by the Dual Field Phase Referencing instrument of the ASTRA project: simultaneously measuring the real-time effects of the atmosphere on a nearby bright guide star, and correcting for it on the faint target, integration time longer than the turbulence time scale are made possible. As a prelude to this demonstration, we first present the implementation of Dual Field Phase Referencing on the interferometer. We then detail its on-sky performance focusing on the accuracy of the turbulence correction, and on the resulting fringe contrast stability. We conclude with a presentation of early results obtained with Laser Guide Star AO and the interferometer.

We present the second data release from the GALEX Arecibo SDSS Survey (GASS), an ongoing large Arecibo program to measure the HI properties for an unbiased sample of ~1000 galaxies with stellar masses greater than 10^10 Msun and redshifts 0.025<z<0.05. GASS targets are selected from the Sloan Digital Sky Survey (SDSS) spectroscopic and Galaxy Evolution Explorer (GALEX) imaging surveys, and are observed until detected or until a gas mass fraction limit of a few per cent is reached. This second data installment includes new Arecibo observations of 240 galaxies, and marks the 50% of the complete survey. We present catalogs of the HI, optical and ultraviolet parameters for these galaxies, and their HI-line profiles. Having more than doubled the size of the sample since the first data release, we also revisit the main scaling relations of the HI mass fraction with galaxy stellar mass, stellar mass surface density, concentration index, and NUV-r color, as well as the gas fraction plane introduced in our earlier work.

We investigate the spatial distribution of galactic satellites in high resolution simulations of structure formation in the LCDM model: the Aquarius dark matter simulations of individual halos and the Millennium II simulation of a large cosmological volume. To relate the simulations to observations of the Milky Way we use two alternative models to populate dark halos with "visible" galaxies: a semi-analytic model of galaxy formation and an abundance matching technique. We find that the radial density profile of massive satellites roughly follows that of the dark matter halo (unlike the distribution of dark matter subhalos). Furthermore, our two galaxy formation models give results consistent with the observed profile of the 11 classical satellites of the Milky Way. Our simulations predict that larger, fainter samples of satellites should still retain this profile at least up to samples of 100 satellites. The angular distribution of the classical satellites of the Milky Way is known to be highly anisotropic. Depending on the exact measure of flattening, 5--10 per cent of satellite systems in our simulations are as flat as the Milky Way's and this fraction does not change when we correct for possible obscuration of satellites by the Galactic disk. A moderate flattening of satellite systems is a general property of LCDM, best understood as the consequence of preferential accretion along filaments of the cosmic web. Accretion of a single rich group of satellites can enhance the flattening due to such anisotropic accretion. We verify that a typical Milky Way-mass CDM halo does not acquire its 11 most massive satellites from a small number of rich groups. Single--group accretion becomes more likely for less massive satellites. Our model predictions should be testable with forthcoming studies of satellite systems in other galaxies and surveys of fainter satellites in the Milky Way.

The flux-ratio anomalies observed in multiply-lensed quasar images are most plausibly explained as the result of perturbing structures superposed on the underlying smooth matter distribution of the primary lens. The cold dark matter cosmological model predicts that a large number of substructures should survive inside larger halos but, surprisingly, this population alone has been shown to be insufficient to explain the observed distribution of the flux ratios of quasar's multiple images. Other halos (and their own subhalos) projected along the line of sight to the primary lens have been considered as additional source of perturbation. In this work, we use ray tracing through the Millennium II simulation to investigate the importance of projection effects due to halos and subhalos of mass m>1E8 Msun/h and extend our analysis to lower masses, m>1E6 Msun/h, using Monte-Carlo halo distributions. We find that the magnitude of the violation depends strongly on the density profile and concentration of the intervening halos, but clustering plays only a minor role. For a typical lensing geometry (lens at redshift 0.6 and source at redshift 2), background haloes (behind the main lens) are more likely to cause a violation than foreground halos. We conclude that line-of-sight structures can be as important as intrinsic substructures in causing flux-ratio anomalies. The combined effect of perturbing structures within the lens and along the line of sight in the LCDM universe results in a cusp-violation probability of 20-30%. This alleviates the discrepancy between models and current data, but a larger observational sample is required for a stronger test of the theory.

We present new deep observations of 'shell' structures in the halo of the nearby elliptical galaxy NGC 7600, alongside a movie of galaxy formation in a cold dark matter universe (available at http://www.virgo.dur.ac.uk/shell-galaxies). The movie, based on an ab initio cosmological simulation, shows how continuous accretion of clumps of dark matter and stars creates a swath of diffuse circumgalactic structures. The disruption of a massive clump on a near-radial orbit creates a complex system of transient concentric shells which bare a striking resemblance to those of NGC 7600. With the aid of the simulation we interpret NGC 7600 in the context of the CDM model.

We characterize substructure in the simulated stellar halos of Cooper et al. (2010) which were formed by the disruption of satellite galaxies within the cosmological N-body simulations of galactic halos of the Aquarius Project. These stellar halos exhibit a wealth of tidal features: broad overdensities and very narrow faint streams akin to those observed around the Milky Way. The substructures are distributed anisotropically on the sky, a characteristic that should become apparent in the next generation of photometric surveys. The normalized RMS of the density of stars on the sky appears to be systematically larger for our halos compared to the value estimated for the Milky Way from main sequence turn-off stars in the Sloan Digital Sky Survey. We show that this is likely to be due in part to contamination by faint QSOs and redder main sequence stars, and might suggest that ~10% of the Milky Way halo stars have formed in-situ.

We present dynamical scaling relations for a homogeneous and representative sample of ~500 massive galaxies, selected only by stellar mass (>10^10 Msun) and redshift (0.025<z<0.05) as part of the ongoing GALEX Arecibo SDSS Survey. We compare baryonic Tully-Fisher (BTF) and Faber-Jackson (BFJ) relations for this sample, and investigate how galaxies scatter around the best fits obtained for pruned subsets of disk-dominated and bulge-dominated systems. The BFJ relation is significantly less scattered than the BTF when the relations are applied to their maximum samples, and is not affected by the inclination problems that plague the BTF. Disk-dominated, gas-rich galaxies systematically deviate from the BFJ relation defined by the spheroids. We demonstrate that by applying a simple correction to the stellar velocity dispersions that depends only on the concentration index of the galaxy, we are able to bring disks and spheroids onto the same dynamical relation -- in other words, we obtain a generalized BFJ relation that holds for all the galaxies in our sample, regardless of morphology, inclination or gas content, and has a scatter smaller than 0.1 dex. We find that disks and spheroids are offset in the stellar dispersion-size relation, and that the offset is removed when corrected dispersions are used instead. The generalized BFJ relation represents a fundamental correlation between the global dark matter and baryonic content of galaxies, which is obeyed by all (massive) systems regardless of morphology. [abridged]

We test whether the masses of old globular clusters (GCs) in dwarf galaxies are consistent with the same initial mass spectrum as young massive clusters (YMCs) in nearby star-forming galaxies. The most massive GCs of dwarf galaxies are compared to their expected masses when drawing from the Schechter-type ICMF of YMCs. It is found that the most massive GCs of galaxies in the stellar mass range M_*,gal=10^7-10^9 Msun are consistent with the same initial mass spectrum as YMCs in about 90% of the cases, suggesting that their formation mechanisms were the same. For the remaining 10%, the most massive clusters are nuclear GCs, which have been able to grow to higher masses through further merging after their initial formation ended. Because the effects of cluster disruption are weaker for more massive clusters, we estimate that up to one third of the metal-poor GCs in the Milky Way may have a nuclear origin, while the remaining two thirds formed through the same process as YMCs in the local universe. A log-normal ICMF is inconsistent with observed GCs at a 99.6% confidence level.

We describe a correlation function statistic that quantifies the amount of spatial and kinematic substructure in the stellar halo. We test this statistic using model stellar halo realizations constructed from the Aquarius suite of six high-resolution N-body simulations in combination with the Galform semi-analytic galaxy formation model. These simulations show considerable scatter in the properties of stellar haloes. We find that our statistic can distinguish between these plausible alternatives for the global structure of the Milky Way stellar halo. We compare with observational data and show that pencil beam surveys of ~100 tracer stars (such as the Spaghetti Survey) are not sufficient to constrain the degree of structure in the Milky Way halo with this statistic. Larger area surveys with >1000 tracer stars (such as BHB stars in the Sloan Digital Sky Survey) provide much tighter constraints on comparisons between models and data. In our simulations, we find examples of haloes with spatial and kinematic substructure consistent with the available Milky Way data.

We use cosmological simulations from the Aquarius Project to study the orbital history of the Large Magellanic Cloud (LMC) and its potential association with other satellites of the Milky Way (MW). We search for dynamical analogs to the LMC and find a subhalo that matches the LMC position and velocity at either of its two most recent pericentric passages. This suggests that the LMC is not necessarily on its first approach to the MW, provided that the virial mass of the Milky Way is M_200 ~1.8e12 M_sun. The simulation results yield specific predictions for the position and velocity of systems associated with the LMC prior to infall. If on first approach, most should lie close to the LMC because the Galactic tidal field has not yet had enough time to disperse them. If on second approach, the list of potential associates increases substantially. Interestingly, our analysis rules out an LMC association for Draco and Ursa Minor, two of the dwarf spheroidals suggested by Lynden-Bell & Lynden-Bell to form part of the "Magellanic Ghostly Stream". Our results also indicate that the direction of the orbital angular momentum is a powerful test of LMC association. This test, however, requires precise proper motions, which are unavailable for most MW satellites. Of the 4 satellites with published proper motions, only the Small Magellanic Cloud is clearly associated with the LMC. Taken at face value, the proper motions of Carina, Fornax and Sculptor rule them out as potential associates, but this conclusion should be revisited when better data become available. The dearth of satellites clearly associated with the Clouds might be solved by wide-field imaging surveys that target its surroundings, a region that may prove a fertile hunting ground for faint, previously unnoticed MW satellites.

We present a model for the satellites of the Milky Way in which galaxy formation is followed using semi-analytic techniques applied to the six high-resolution N-body simulations of galactic halos of the Aquarius project. The model, calculated using the Galform code, incorporates improved treatments of the relevant physics in the LambdaCDM cosmogony, particularly a self-consistent calculation of reionization by UV photons emitted by the forming galaxy population, including the progenitors of the central galaxy. Along the merger tree of each halo, the model calculates gas cooling (by Compton scattering off cosmic microwave background photons, molecular hydrogen and atomic processes), gas heating (from hydrogen photoionization and supernova energy), star formation and evolution. The evolution of the intergalactic medium is followed simultaneously with that of the galaxies. Star formation in the more massive progenitor subhalos is suppressed primarily by supernova feedback, while for smaller subhalos it is suppressed primarily by photoionization due to external and internal sources. The model is constrained to match a wide range of properties of the present day galaxy population as a whole, but at high redshift it requires an escape fraction of UV photons near unity in order completely to reionize the universe by redshift z ~ 8. In the most successful model the local sources photoionize the pre-galactic region completely by z ~ 10. In addition to the luminosity function of Milky Way satellites, the model matches their observed luminosity-metallicity relation, their radial distribution and the inferred values of the mass within 300 pc, which in the models increase slowly but significantly with luminosity. There is a large variation in satellite properties from halo to halo, with the luminosity function, for example, varying by a factor of ~ 2 among the six simulations.

We report on the successful science verification phase of a new observing mode at the Keck interferometer, which provides a line-spread function width and sampling of 150km/s at K'-band, at a current limiting magnitude of K'~7mag with spatial resolution of lam/2B ~2.7mas and a measured differential phase stability of unprecedented precision (3mrad at K=5mag, which represents 3uas on sky or a centroiding precision of 10^-3). The scientific potential of this mode is demonstrated by the presented observations of the circumstellar disk of the evolved Be-star 48Lib. In addition to indirect methods such as multi-wavelength spectroscopy and polaritmetry, the here described spectro-interferometric astrometry provides a new tool to directly constrain the radial density structure in the disk. We resolve for the first time several Pfund emission lines, in addition to BrGam, in a single interferometric spectrum, and with adequate spatial and spectral resolution and precision to analyze the radial disk structure in 48Lib. The data suggest that the continuum and Pf-emission originates in significantly more compact regions, inside of the BrGam emission zone. Thus, spectro-interferometric astrometry opens the opportunity to directly connect the different observed line profiles of BrGam and Pfund in the total and correlated flux to different disk radii. The gravitational potential of a rotationally flattened Be star is expected to induce a one-armed density perturbation in the circumstellar disk. Such a slowly rotating disk oscillation has been used to explain the well known periodic V/R spectral profile variability in these stars, as well as the observed V/R cycle phase shifts between different disk emission lines. The differential line properties and linear constraints set by our data lend support to the existence of a radius-dependent disk density perturbation.

Lensing flux-ratio anomalies have been frequently observed and taken as evidence for the presence of abundant dark matter substructures in lensing galaxies, as predicted by the cold dark matter (CDM) model of cosmogony. In previous work, we examined the cusp-caustic relations of the multiple images of background quasars lensed by galaxy-scale dark matter haloes, using a suite of high-resolution N-body simulations (the Aquarius simulations). In this work, we extend our previous calculations to incorporate both the baryonic and diffuse dark components in lensing haloes. We include in each lensing simulation: (1) a satellite galaxy population derived from a semi-analytic model applied to the Aquarius haloes, (2) an empirical Milky-Way globular cluster population and (3) satellite streams (diffuse dark component) identified in the simulations. Accounting for these extra components, we confirm our earlier conclusion that the abundance of intrinsic substructures (dark or bright, bound or diffuse) is not sufficient to explain the observed frequency of cusp-caustic violations in the CLASS survey. We conclude that the observed effect could be the result of the small number statistics of CLASS, or intergalactic haloes along the line of sight acting as additional sources of lensing flux anomalies. Another possibility is that this discrepancy signals a failure of the CDM model.

We present six simulations of galactic stellar haloes formed by the tidal disruption of accreted dwarf galaxies in a fully cosmological setting. Our model is based on the Aquarius project, a suite of high resolution N-body simulations of individual dark matter haloes. We tag subsets of particles in these simulations with stellar populations predicted by the Galform semi-analytic model. Our method self-consistently tracks the dynamical evolution and disruption of satellites from high redshift. The luminosity function and structural properties of surviving satellites, which agree well with observations, suggest that this technique is appropriate. We find that accreted stellar haloes are assembled between 1<z<7 from less than 5 significant progenitors. These progenitors are old, metal-rich satellites with stellar masses similar to the brightest Milky Way dwarf spheroidals. In contrast to previous stellar halo simulations, we find that several of these major contributors survive as self-bound systems to the present day. Both the number of these significant progenitors and their infall times are inherently stochastic. This results in great diversity among our stellar haloes, which amplifies small differences between the formation histories of their dark halo hosts. The masses and density/surface-brightness profiles of the stellar haloes are consistent with expectations from the Milky Way and M31. Each halo has a complex structure, consisting of well-mixed components, tidal streams, shells and other subcomponents. This structure is not adequately described by smooth models. We find one example of an accreted thick disk. Contrasts in age and metallicity between halo stars and those in surviving satellites are in broad agreement with recent observations. [Abridged]

The sensitivity and astrometry upgrade ASTRA of the Keck Interferometer is introduced. After a brief overview of the underlying interferometric principles, the technology and concepts of the upgrade are presented. The interferometric dual-field technology of ASTRA will provide the KI with the means to observe two objects simultaneously, and measure the distance between them with a precision eventually better than 100 uas. This astrometric functionality of ASTRA will add a unique observing tool to fields of astrophysical research as diverse as exo-planetary kinematics, binary astrometry, and the investigation of stars accelerated by the massive black hole in the center of the Milky Way as discussed in this contribution.

We explore the form of extragalactic reddening toward quasars using a sample of 9566 quasars with redshifts 0<z<2.2, and accurate optical colors from the Sloan Digital Sky Survey (SDSS). We confirm that dust reddening is the primary explanation for the red ``tail'' of the color distribution of SDSS quasars. Our fitting to 5-band photometry normalized by the modal quasar color as a function of redshift shows that this ``tail'' is well described by SMC-like reddening but not by LMC-like, Galactic, or Gaskell et al. (2004) reddening. Extension to longer wavelengths using a subset of 1886 SDSS-2MASS matches confirms these results at high significance. We carry out Monte-Carlo simulations that match the observed distribution of quasar spectral energy distributions using a Lorentzian dust reddening distribution; 2% of quasars selected by the main SDSS targeting algorithm (i.e., which are not extincted out of the sample) have E_{B-V} > 0.1; less than 1% have E_{B-V} > 0.2, where the extinction is relative to quasars with modal colors. Reddening is uncorrelated with the presence of intervening narrow-line absorption systems, but reddened quasars are much more likely to show narrow absorption at the redshift of the quasar than are unreddened quasars. Thus the reddening towards quasars is dominated by SMC-like dust at the quasar redshift.

There exists a class of cosmic strings that turn matter into antimatter (Alice strings). In a GUT where the unbroken gauge group contains charge conjugation ($C$), such strings form when a phase transition renders $C$ a discrete symmetry. They become boundaries of domain walls at a later, $C$-breaking transition. These `Alice walls' are cosmologically harmless, but can play an important role in baryogenesis. We present a three-generation toy model with scalar baryons, where a quasi-static Alice wall (or a gas of such walls) temporarily gives rise to net baryogenesis of uniform sign everywhere in space. This becomes a permanent baryon excess if the wall shrinks away early enough. We comment on the possible relevance of a similar mechanism to baryogenesis in a realistic $\soten$ unification model, where Alice walls would form at the scale of left-right symmetry breaking.