Source author record

Raymond G. Carlberg

Raymond G. Carlberg appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

32works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

32 published item(s)

preprint2022arXiv

A stellar stream remnant of a globular cluster below the metallicity floor

Stellar ejecta gradually enrich the gas out of which subsequent stars form, making the least chemically enriched stellar systems direct fossils of structures formed in the early universe. Although a few hundred stars with metal content below one thousandth of the solar iron content are known in the Galaxy, none of them inhabit globular clusters, some of the oldest known stellar structures. These show metal content of at least ~0.2 percent of the solar metallicity ([Fe/H] > -2.7). This metallicity floor appears universal and it has been proposed that proto-galaxies that merge into the galaxies we observe today were simply not massive enough to form clusters that survived to the present day. Here, we report the discovery of a stellar stream, C-19, whose metallicity is less than 0.05 per cent the solar metallicity ([Fe/H]=-3.38 +/- 0.06 (stat.) +/- 0.20 (syst.)). The low metallicity dispersion and the chemical abundances of the C-19 stars show that this stream is the tidal remnant of the most metal-poor globular cluster ever discovered, and significantly below the purported metallicity floor: clusters with significantly lower metallicities than observed today existed in the past and contributed their stars to the Milky Way halo.

preprint2022arXiv

Infant-phase reddening by surface Fe-peak elements in a normal Type Ia Supernova

Type Ia Supernovae are thermonuclear explosions of white dwarf stars. They play a central role in the chemical evolution of the Universe and are an important measure of cosmological distances. However, outstanding questions remain about their origins. Despite extensive efforts to obtain natal information from their earliest signals, observations have thus far failed to identify how the majority of them explode. Here, we present infant-phase detections of SN 2018aoz from a brightness of -10.5 absolute AB magnitudes -- the lowest luminosity early Type Ia signals ever detected -- revealing a hitherto unseen plateau in the $B$-band that results in a rapid redward color evolution between 1.0 and 12.4 hours after the estimated epoch of first light. The missing $B$-band flux is best-explained by line-blanket absorption from Fe-peak elements in the outer 1% of the ejected mass. The observed $B-V$ color evolution of the SN also matches the prediction from an over-density of Fe-peak elements in the same outer 1% of the ejected mass, whereas bluer colors are expected from a purely monotonic distribution of Fe-peak elements. The presence of excess nucleosynthetic material in the extreme outer layers of the ejecta points to enhanced surface nuclear burning or extended sub-sonic mixing processes in some normal Type Ia Supernova explosions.

preprint2022arXiv

The Dark Matter Halo of M54

M54 is a prototype of a globular cluster embedded in a dark matter halo. Gaia EDR3 photometry and proper motions separate the old, metal-poor stars from the more metal rich and younger dwarf galaxy stars. The metal poor stars dominate the inner 50 pc, with a velocity dispersion profile that declines to a minimum around 30 pc then rises back to nearly the central velocity dispersion, as expected for a globular cluster at the center of a dark matter halo. The Jeans analysis of the three separate stellar populations give consistent masses that rise approximately linearly with radius to 1 kpc, implying a small core or cuspy halo. These data are compatible with an infalling CDM dark matter halo reduced to 3x10^8 M_sun at the 50 kpc apocenter 2.3 Gyr ago, with a central globular cluster surrounded by the remnant of a dwarf galaxy. Tides gradually remove material beyond 1 kpc but have little effect on the stars and dark matter within 300 pc of the center. M54 appears to be a transitional system between globular clusters with and without local dark halos, whose evolution within the galaxy depends sensitively on the time of accretion and orbital pericenter.

preprint2022arXiv

The Global Dynamical Atlas of the Milky Way mergers: Constraints from Gaia EDR3 based orbits of globular clusters, stellar streams and satellite galaxies

The Milky Way halo was predominantly formed by the merging of numerous progenitor galaxies. However, our knowledge of this process is still incomplete, especially in regard to the total number of mergers, their global dynamical properties and their contribution to the stellar population of the Galactic halo. Here, we uncover the Milky Way mergers by detecting groupings of globular clusters, stellar streams and satellite galaxies in action ($\mathbf{J}$) space. While actions fully characterize the orbits, we additionally use the redundant information on their energy ($\textit{E}$) to enhance the contrast between groupings. For this endeavour, we use $\textit{Gaia}$ EDR3 based measurements of $170$ globular clusters, $41$ streams and $46$ satellites to derive their $\mathbf{J}$ and $\textit{E}$. To detect groups, we use the $\texttt{ENLINK}$ software, coupled with a statistical procedure that accounts for the observed phase-space uncertainties of these objects. We detect a total of $N=6$ groups, including the previously known mergers $\textit{Sagittarius}$, $\textit{Cetus}$, $\textit{Gaia-Sausage/Enceladus}$, $\textit{LMS-1/Wukong}$, $\textit{Arjuna/Sequoia/I'itoi}$ and one new merger that we call $\textit{Pontus}$. All of these mergers, together, comprise $62$ objects ($\approx 25\%$ of our sample). We discuss their members, orbital properties and metallicity distributions. We find that the three most metal-poor streams of our Galaxy -- "C-19" ([Fe/H]$=-3.4$ dex), "Sylgr" ([Fe/H]$=-2.9$ dex) and "Phoenix" ([Fe/H]$=-2.7$ dex) -- are associated with $\textit{LMS-1/Wukong}$; showing it to be the most metal-poor merger. The global dynamical atlas of Milky Way mergers that we present here provides a present-day reference for galaxy formation models.

preprint2021arXiv

Simulating Globular Clusters in Dark Matter Sub-Halos

A cosmological zoom-in simulation which develops into a Milky Way-like halo is started at redshift 7. The initial dark matter distribution is seeded with dense star clusters, median mass $5\times 10^5 M_\sun$, placed in the largest sub-halos present, which have a median peak circular velocity of 25 \kms. Three simulations are initialized using the same dark matter distribution, with the star clusters started on approximately circular orbits having initial median radii 6.8 kpc, 0.14 kpc, and, at the exact center of the sub-halos. The simulations are evolved to the current epoch at which time the median galactic orbital radii of the three sets of clusters are 30, 5 and 16 kpc, with the clusters losing about 2, 50 and 15\% of their mass, respectively. Clusters started at small orbital radii have so much tidal forcing that they are often not in equilibrium. Clusters started at larger sub-halo radii have a velocity dispersion that declines smoothly to $\simeq$20\% of the central value at $\simeq$20 half mass radii. The clusters started at the sub-halo centers can show a rise in velocity dispersion beyond 3-5 half mass radii. That is, the clusters formed without local dark matter always have stellar mass dominated kinematics at all radii, whereas about 25\% of the clusters started at sub-halo centers have remnant local dark matter.

preprint2020arXiv

A Deep CFHT Optical Search for a Counterpart to the Possible Neutron Star -- Black Hole Merger GW190814

We present a wide-field optical imaging search for electromagnetic counterparts to the likely neutron star - black hole (NS-BH) merger GW190814/S190814bv. This compact binary merger was detected through gravitational waves by the LIGO/Virgo interferometers, with masses suggestive of a NS-BH merger. We imaged the LIGO/Virgo localization region using the MegaCam instrument on the Canada-France-Hawaii Telescope. We describe our hybrid observing strategy of both tiling and galaxy-targeted observations, as well as our image differencing and transient detection pipeline. Our observing campaign produced some of the deepest multi-band images of the region between 1.7 and 8.7 days post-merger, reaching a 5sigma depth of g > 22.8 (AB mag) at 1.7 days and i > 23.1 and i > 23.9 at 3.7 and 8.7 days, respectively. These observations cover a mean total integrated probability of 67.0% of the localization region. We find no compelling candidate transient counterparts to this merger in our images, which suggests that either the lighter object was tidally disrupted inside of the BH's innermost stable circular orbit, the transient lies outside of the observed sky footprint, or the lighter object is a low-mass BH. We use 5sigma source detection upper limits from our images in the NS-BH interpretation of this merger to constrain the mass of the kilonova ejecta to be Mej < 0.015Msun for a 'blue' (kappa = 0.5 cm^2 g^-1) kilonova, and Mej < 0.04Msun for a 'red' (kappa = 5-10 cm^2 g^-1) kilonova. Our observations emphasize the key role of large-aperture telescopes and wide-field imagers such as CFHT MegaCam in enabling deep searches for electromagnetic counterparts to gravitational wave events.

preprint2020arXiv

Dynamical Simulations of the First Globular Clusters

Tidally limited star clusters are started on nearly circular orbits in the dark matter sub-halos present at high redshift and evolved with an n-body code augmented with gravitational interactions in the clusters. The tidally removed stars and the remaining high redshift clusters from a redshift 8 start are more concentrated than the dark matter, as expected. However, the subset of stars from the clusters that began in the lower mass sub-halos have a distribution somewhat more extended than the dark matter halo, with a mean galactic radius of about 60 kpc inside 150 kpc. The clusters from low mass sub-halos, those with a peak circular velocity of $ 12-18$ \kms, also produce most of the population's thin stellar streams. The dependence of the stellar population distribution on sub-halo mass is not seen in simulations that start clusters at lower redshift. The half mass radii of the clusters are set by the tidal fields of the initial cluster orbits, which leads to the average half mass radius decreasing with increased redshift of formation. Starting clusters at greater than redshift 4 leads to cluster half mass radii approximately as seen in the Milky Way, with lower redshifts leading to clusters that are too large. If extremely metal poor globular clusters were preferentially formed in high redshift, low mass, dwarf galaxies, then the clusters should be spread nearly uniformly over the dark halo, along with their tidally removed stars.

preprint2020arXiv

The Likelihood of Undiscovered Globular Clusters in the Outskirts of the Milky Way

The currently known Galactic globular cluster population extends out to a maximum galactocentric distance of $\sim$ 145 kpc, with the peculiarity that the outermost clusters predominantly have an inward velocity. Orbit averaging finds that this configuration occurs by chance about $6\%$ of the time, suggesting that several globular clusters with positive radial velocities remain undiscovered. We evaluate the expected number of undiscovered clusters at large distances under the assumption that the cluster population has a smooth radial distribution and is in equilibrium within the Milky Way's virial radius. By comparing the present day kinematic properties of outer clusters to random orbital configurations of the Galactic globular cluster system through orbit averaging, we estimate a likelihood of $73\%$ of there being at least one undiscovered globular cluster within the Milky Way. This estimate assumes the current population is complete out to 50 kpc, and increases to $91\%$ if the population is complete out to 150 kpc. The likelihood of there being two undiscovered clusters is between $60\%$ and $70\%$, with the likelihood of there being three undiscovered clusters being on the order of $50\%$. The most likely scenario is that the undiscovered clusters are moving outwards, which results in the outer cluster population being consistent with an equilibrium state. Searches for distant and possibly quite low concentration and very low metallicity globular clusters will be enabled with upcoming deep imaging surveys.

preprint2020arXiv

The Pristine Inner Galaxy Survey (PIGS) II: Uncovering the most metal-poor populations in the inner Milky Way

Metal-poor stars are important tools for tracing the early history of the Milky Way, and for learning about the first generations of stars. Simulations suggest that the oldest metal-poor stars are to be found in the inner Galaxy. Typical bulge surveys, however, lack low metallicity ([Fe/H] < -1.0) stars because the inner Galaxy is predominantly metal-rich. The aim of the Pristine Inner Galaxy Survey (PIGS) is to study the metal-poor and very metal-poor (VMP, [Fe/H] < -2.0) stars in this region. In PIGS, metal-poor targets for spectroscopic follow-up are selected from metallicity-sensitive CaHK photometry from the CFHT. This work presents the ~250 deg^2 photometric survey as well as intermediate-resolution spectroscopic follow-up observations for ~8000 stars using AAOmega on the AAT. The spectra are analysed using two independent tools: ULySS with an empirical spectral library, and FERRE with a library of synthetic spectra. The comparison between the two methods enables a robust determination of the stellar parameters and their uncertainties. We present a sample of 1300 VMP stars -- the largest sample of VMP stars in the inner Galaxy to date. Additionally, our spectroscopic dataset includes ~1700 horizontal branch stars, which are useful metal-poor standard candles. We furthermore show that PIGS photometry selects VMP stars with unprecedented efficiency: 86%/80% (lower/higher extinction) of the best candidates satisfy [Fe/H] < -2.0, as do 80%/63% of a larger, less strictly selected sample. We discuss future applications of this unique dataset that will further our understanding of the chemical and dynamical evolution of the innermost regions of our Galaxy.

preprint2019arXiv

Milky Way Halo Vibrations and Incommensurate Stream Velocities

Collisionless dark matter galactic halos are expected to exhibit damped oscillations as a result of ongoing late time accretion. An n-body model of the cosmological assembly of a Milky Way-like halo is used to quantify the time dependence of its gravitational field. The simulation contains stellar streams whose incommensurate perpendicular velocities are found to have an approximately exponential distribution with a scale of 10-20\kms, depending on how the stars are selected, comparable to those reported for the Orphan stream. The fluctuations in the quadrupole moment of the dark matter halo are sufficient to largely explain the tangential velocities. If velocity measurements of a larger sample of Milky Way streams finds (or does not find) the expected distribution of transverse velocities it will lead to limits on the cross-section of self-interacting dark matter, in which kinetic viscosity can damp the oscillations more rapidly than the mixing processes of collisionless dark matter alone.

preprint2019arXiv

Phase-space correlation in stellar streams of the Milky Way halo: The clash of Kshir and GD-1

We report the discovery of a 70-degree long stellar stream in the Milky Way halo, which criss-crosses the well known "GD-1" stream. We show that this new stellar structure ("Kshir") and GD-1 lie at similar distance, and are remarkably correlated in kinematics. We propose several explanations for the nature of this new structure and its possible association with GD-1. However, a scenario in which these two streams were accreted onto the Milky Way within the same dark matter sub-halo seems to provide a natural explanation for their phase-space entanglement, and other complexities of this coupled-system.

preprint2019arXiv

The Canada-France Imaging Survey: Reconstructing the Milky Way Star Formation History from its White Dwarf Population

As the remnants of stars with initial masses $\lesssim$ 8 M$_{\odot}$, white dwarfs contain valuable information on the formation histories of stellar populations. In this paper, we use deep, high-quality, u-band photometry from the Canada France Imaging Survey (CFIS), griz photometry from Pan-STARRS 1 (PS1), as well as proper motions from Gaia DR2, to select 25,156 white dwarf candidates over $\sim$4500 deg$^2$ using a reduced proper motion diagram. We develop a new white dwarf population synthesis code that returns mock observations of the Galactic field white dwarf population for a given star formation history, while simultaneously taking into account the geometry of the Milky Way, survey parameters, and selection effects. We use this model to derive the star formation histories of the thin disk, thick disk, and stellar halo. Our results show that the Milky Way disk began forming stars (11.3 $\pm$ 0.5) Gyr ago, with a peak rate of (8.8 $\pm$ 1.4) M$_{\odot}$yr$^{-1}$ at (9.8 $\pm$ 0.4) Gyr, before a slow decline to a constant rate until the present day --- consistent with recent results suggesting a merging event with a satellite galaxy. Studying the residuals between the data and best-fit model shows evidence for a slight increase in star formation over the past 3 Gyr. We fit the local fraction of helium-atmosphere white dwarfs to be (21 $\pm$ 3) %. Incorporating this methodology with data from future wide-field surveys such as LSST, Euclid, CASTOR, and WFIRST should provide an unprecedented view into the formation of the Milky Way at its earliest epoch through its white dwarfs.

preprint2019arXiv

The Density Structure of Simulated Stellar Streams

Simulations are run with and without a normal cold dark matter sub-halo population below 4x10^8 M_sun to examine the role of the lower mass sub-halos in the creation of density variations, "gaps", within thin tidal star streams. Dense star clusters composed of stellar mass particles are started within a redshift 3 sub-galactic dark matter halo distribution which assembles to create a Milky Way like halo. Hierarchical assembly causes the early-time tidal star streams from the clusters to later be dispersed as thick streams around the thin streams that develop in the main late-time halo as it settles towards equilibrium. The two simulations are set up to be essentially dynamically identical and result in about 140 star clusters within a radius of 50 kpc, of which 40 some have thin tidal streams longer than 40\degr\ with a separation of more than 20\degr\ from the cluster as viewed from the halo center. Averaged over about 20 separate snapshots, the fraction of streams with significant density variations inside 50 kpc is $0.58\pm0.02$ in the simulation with low mass sub-halos and $0.52\pm0.02$ in the simulation without low mass halos. That is, only $12\pm5$\% of the stream density variations can be attributed to the sub-halos below 4x10^8 M_sun. The dominant source of density variations in streams is the varying tidal fields along the complex trajectories of the clusters and their streams as the galaxy assembles.

preprint2019arXiv

The Pristine survey X: a large population of low-metallicity stars permeates the Galactic disk

The orbits of the least chemically enriched stars open a window on the formation of our Galaxy when it was still in its infancy. The common picture is that these low-metallicity stars are distributed as an isotropic, pressure-supported component since these stars were either accreted from the early building blocks of the assembling Milky Way, or were later brought by the accretion of faint dwarf galaxies. Combining the metallicities and radial velocities from the Pristine and LAMOST surveys and Gaia DR2 parallaxes and proper motions for an unprecedented large and unbiased sample of very metal-poor stars at $[Fe/H]\leq-2.5$ we show that this picture is incomplete. This sample shows strong statistical evidence (at the $5.0σ$ level) of asymmetry in their kinematics, favouring prograde motion. Moreover, we find that $31\%$ of the stars that currently reside in the disk do not venture outside of the disk plane throughout their orbit. The discovery of this population implies that a significant fraction of stars with iron abundances $[Fe/H]\leq-2.5$ formed within or concurrently with the Milky Way disk and that the history of the disk was quiet enough to allow them to retain their disk-like orbital properties.

preprint2016arXiv

Dispersal of tidal debris in a Milky-Way-sized dark matter halo

We simulate the tidal disruption of a collisionless N-body globular star cluster in a total of 300 different orbits selected to have galactocentric radii between 10 and 30 kpc in four dark matter halos: (a) a spherical halo with no subhalos, (b) a spherical halo with subhalos, (c) a realistic halo with no subhalos, and (d) a realistic halo with subhalos. This allows us to isolate and study how the halo's (lack of) dynamical symmetry and substructures affect the dispersal of tidal debris. The realistic halos are constructed from the snapshot of the Via Lactea II simulation at redshift zero. We find that the overall halo's lack of dynamical symmetry disperses tidal debris to make the streams fluffier, consistent with previous studies of tidal debris of dwarf galaxies in larger orbits than ours in this study. On the other hand, subhalos in realistic potentials can locally enhance the densities along streams, making streams denser than their counterparts in smooth potentials. We show that many long and thin streams can survive in a realistic and lumpy halo for a Hubble time. This suggests that upcoming stellar surveys will likely uncover more thin streams which may contain density gaps that have been shown to be promising probes for dark matter substructures.

preprint2016arXiv

Modeling GD-1 Gaps in a Milky-Way Potential

The GD-1 star stream is currently the best available for identifying density fluctuations, "gaps", along its length as a test of the LCDM prediction of large numbers of dark matter sub-halos orbiting in the halo. Density variations of some form are present, since the variance of the density along the stream is three times that expected from the empirically estimated variation in the filtered mean star counts. The density variations are characterized with filters that approximate the shape of sub-halo gravitationally induced stream gaps. The filters locate gaps and measures their amplitude, leading to a measurement of the distribution of gap widths. To gain understanding of the factors influencing the gap width distribution, a suite of collisionless n-body simulations for a GD-1 like orbit in a Milky Way-like potential provides a dynamically realistic statistical prediction of the gap distribution. The simulations show that every location in the stream has been disturbed to some degree by a sub-halo. The small gaps found via the filtering are largely noise. Larger gaps, those longer than 1 kpc, or 10 degrees for GD-1, are the source of the excess variance. The suite of stream simulations shows that sub-halos at the predicted inner halo abundance or possibly somewhat higher can produce the required large sale density variations.

preprint2016arXiv

The DESI Experiment Part I: Science,Targeting, and Survey Design

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-$α$ 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.

preprint2016arXiv

The DESI Experiment Part II: Instrument Design

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= λ/Δλ$ 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.

preprint2016arXiv

Tidal Debris as a Dark Matter Probe

Tidal debris streams from galaxy satellites can provide insight into the dark matter distribution in halos. This is because we have more information about stars in a debris structure than about a purely random population of stars: we know that in the past they were all bound to the same dwarf galaxy; and we know that they form a dynamically cold population moving on similar orbits. They also probe a different region of the matter distribution in a galaxy than many other methods of mass determination, as their orbits take them far beyond the typical extent of those for the bulk of stars. Although conclusive results from this information have yet to be obtained, significant progress has been made in developing the methodologies for determining both the global mass distribution of the Milky Way's dark matter halo and the amount of dark matter substructure within it. Methods for measuring the halo shape are divided into "predictive methods," which predict the tidal debris properties from the progenitor satellite's mass and orbit, given an assumed parent galaxy mass distribution; and "fundamental methods," which exploit properties fundamental to the nature of tidal debris as global potential constraints. Methods for quantifying the prevalence of dark matter subhalos within halos through the analysis of the gaps left in tidal streams after these substructures pass through them are reviewed.

preprint2016arXiv

What a Tangled Web We Weave: Hermus as the Northern Extension of the Phoenix Stream

We investigate whether the recently discovered Phoenix stream may be part of a much longer structure that includes the previously discovered Hermus stream. Using a simple model of the Galaxy with a disk, bulge, and a spherical dark matter halo, we show that a nearly circular orbit, highly inclined with respect to the disk, can be found that fits the positions, orientations, and distances of both streams. While the two streams are somewhat misaligned in the sense that they do not occupy the same plane, nodal precession due to the Milky Way disk potential naturally brings the orbit into line with each stream in the course of half an orbit. We consequently consider a common origin for the two streams as plausible. Based on our best fitting orbit, we make predictions for the positions, distances, radial velocities, and proper motions along each stream. If our hypothesis is borne out by measurements, then at ~183 degrees (~235 degrees with respect to the Galactic center) and ~76 kpc in length, Phoenix-Hermus would become the longest cold stream yet found. This would make it a particularly valuable new probe of the shape and mass of the Galactic halo out to ~20 kpc.

preprint2015arXiv

An Orphan No Longer? Detection of the Southern Orphan Stream and a Candidate Progenitor

Using a shallow, two-color survey carried out with the Dark Energy Camera, we detect the southern, possibly trailing arm of the Orphan Stream. The stream is reliably detected to a declination of $-38^\circ$, bringing the total known length of the Orphan stream to $108^\circ$. We find a slight offset or "S" shape in the stream at $δ\simeq -14^\circ$ that would be consistent with the transition from leading to trailing arms. This coincides with a moderate concentration of $137 \pm 25$ stars (to $g = 21.6$) that we consider a possible remnant of the Orphan progenitor. The position of this feature is in agreement with previous predictions.

preprint2015arXiv

Simulating tidal streams in a high resolution dark matter halo

We simulate tidal streams in the presence and absence of substructures inside the zero redshift snapshot of the Via Lactea II (VL-2) simulation. A halo finder is used to remove and isolate the subhalos found inside the high resolution dark matter halo of VL-2, and the potentials for both the main halo and all the subhalos are constructed individually using the self-consistent field (SCF) method. This allows us to make direct comparison of tidal streams between a smooth halo and a lumpy halo without assuming idealized profiles or triaxial fits. We simulate the kinematics of a star cluster starting with the same orbital position but two different velocities. Although these two orbits are only moderately eccentric and have similar apo- and pericentric distances, we find that the two streams have very different morphologies. We conclude that our model of the potential of VL-2 can provide insights about tidal streams that have not been explored by previous studies using idealized or axisymmetric models.

preprint2015arXiv

Star Streams in Triaxial Isochrone Potentials with Sub-Halos

The velocity, position, and action variable evolution of a tidal stream drawn out of a star cluster in a triaxial isochrone potential containing a sub-halo population reproduces many of the orbital effects of more general cosmological halos but allows easy calculation of orbital actions. We employ a spherical shell code which we show accurately reproduces the results of a tree gravity code for a collisionless star cluster. Streams from clusters on high eccentricity orbits, $e\gtrsim 0.6$, can spread out so much that the amount of material at high enough surface density to stand out on the sky may be only a few percent of the stream's total mass. Low eccentricity streams remain more spatially coherent, but sub-halos both broaden the stream and displace the centerline with details depending on the orbits allowed within the potential. Overall, the majority of stream particles have changes in their total actions of only 1-2\%, leaving the mean stream relatively undisturbed. A halo with 1\% of the mass in sub-halos typically spreads the velocity distribution about a factor of two wider than would be expected for a smooth halo. Strong density variations, "gaps", along with mean velocity offsets, are clearly detected in low eccentricity streams for even a 0.2\% sub-halo mass fraction. Around one hundred velocity measurements per kiloparsec of stream will enable tests for the presence of a local sub-halo density as small as 0.2-0.5\% of the local mass density, with about 1\% predicted for 30 kiloparsec orbital radii streams.

preprint2014arXiv

Dynamical Properties of Collisionless Star Streams

A sufficiently extended satellite in the tidal field of a host galaxy loses mass to create nearly symmetric leading and trailing tidal streams. We study the case in which tidal heating drives mass loss from a low mass satellite. The stream effectively has two dynamical components, a common angular momentum core superposed with episodic pulses with a broader angular momentum distribution. The pulses appear as spurs on the stream, oscillating above and below the stream centerline, stretching and blurring in configuration space as they move away from the cluster. Low orbital eccentricity streams are smoother and have less differential motion than high eccentricity streams. The tail of a high eccentricity stream can develop a fan of particles which wraps around at apocenter in a shell feature. We show that scaling the essentially stationary action-angle variables with the cube root of the satellite mass allows a low mass satellite stream to accurately predict the features in the stream from a satellite a thousand times more massive. As a practical astrophysical application, we demonstrate that narrow gaps in a moderate eccentricity stream, such as GD-1, blur out to 50\% contrast over approximately 6 radial periods. A high eccentricity stream, such as Pal~5, will blur small gaps in only radial 2 orbits as can be understood from the much larger dispersion of angular momentum in the stream.

preprint2014arXiv

Using gaps in N-body tidal streams to probe missing satellites

We use N-body simulations to model the tidal disruption of a star cluster in a Milky-Way-sized dark matter halo, which results in a narrow stream comparable to (but slightly wider than) Pal-5 or GD-1. The mean Galactic dark matter halo is modeled by a spherical Navarro-Frenk-White (NFW) potential with subhalos predicted by the LCDM cosmological model. The distribution and mass function of the subhalos follow the results from the Aquarius simulation. We use a matched filter approach to look for "gaps" in tidal streams at 12 length scales from 0.1 kpc to 5 kpc, which appear as characteristic dips in the linear densities along the streams. We find that, in addition to the subhalos' perturbations, the epicyclic overdensities (EOs) due to the coherent epicyclic motions of particles in a stream also produce gap-like signals near the progenitor. We measure the gap spectra - the gap formation rates as functions of gap length - due to both subhalo perturbations and EOs, which have not been accounted for together by previous studies. Finally, we project the simulated streams onto the sky to investigate issues when interpreting gap spectra in observations. In particular, we find that gap spectra from low signal-to-noise observations can be biased by the orbital phase of the stream. This indicates that the study of stream gaps will benefit greatly from high-quality data from future missions.

preprint2013arXiv

Gaps in the GD-1 Star Stream

GD-1 is a long, thin, Milky Way star stream that has readily visible density variations along its length. We quantify the locations, sizes and statistical significance of the density structure, i.e. gaps, using a set of scaled filters. The shapes of the filters are based on the gaps that develop in simulations of dark matter sub-halos crossing a star stream. The high Galactic latitude 8.4 kpc long segment of GD-1 that we examine has 8+/-3 gaps of 99% significance or greater, with the error estimated on the basis of tests of the gap-filtering technique. The cumulative distribution of gaps more than 3 times the width of the stream is in good agreement with predictions for dark matter sub-halo encounters with cold star streams. The number of gaps narrower than 3 times the width of the GD-1 stream falls well below the cold stream prediction which is taken into account for the gap creation rate integrated over all sizes. Simple warm stream simulations scaled to GD-1 show that the falloff in gaps is expected for sub-halos below a mass of 10^6 M_sun. The GD-1 gaps requires 100 sub-halos >10^6 M_sun within 30 kpc, the apocenter of GD-1 orbit. These results are consistent with LCDM sub-halo predictions but further improvements in stream signal-to-noise and gap modeling will be welcome.

preprint2013arXiv

The Next Generation Virgo Cluster Survey - Infrared (NGVS-IR): I. A new Near-UV/Optical/Near-IR Globular Cluster selection tool

The NGVS-IR project (Next Generation Virgo Survey - Infrared) is a contiguous near-infrared imaging survey of the Virgo cluster of galaxies. It complements the optical wide-field survey of Virgo (NGVS). The current state of NGVS-IR consists of Ks-band imaging of 4 deg^2 centered on M87, and J and Ks-band imaging of 16 deg^2 covering the region between M49 and M87. In this paper, we present the observations of the central 4 deg^2 centered on Virgo's core region. The data were acquired with WIRCam on the Canada-France-Hawaii Telescope and the total integration time was 41 hours distributed in 34 contiguous tiles. A survey-specific strategy was designed to account for extended galaxies while still measuring accurate sky brightness within the survey area. The average 5σlimiting magnitude is Ks=24.4 AB mag and the 50% completeness limit is Ks=23.75 AB mag for point source detections, when using only images with better than 0.7" seeing (median seeing 0.54"). Star clusters are marginally resolved in these image stacks, and Virgo galaxies with μ_Ks=24.4 AB mag arcsec^-2 are detected. Combining the Ks data with optical and ultraviolet data, we build the uiK color-color diagram which allows a very clean color-based selection of globular clusters in Virgo. This diagnostic plot will provide reliable globular cluster candidates for spectroscopic follow-up campaigns needed to continue the exploration of Virgo's photometric and kinematic sub-structures, and will help the design of future searches for globular clusters in extragalactic systems. Equipped with this powerful new tool, future NGVS-IR investigations based on the uiK diagram will address the mapping and analysis of extended structures and compact stellar systems in and around Virgo galaxies.

preprint2012arXiv

Counting Dark Sub-halos with Star Stream Gaps

The Cold Dark Matter paradigm predicts vast numbers of dark matter sub-halos to be orbiting in galactic halos. The sub-halos are detectable through the gaps they create gaps in stellar streams. The gap-rate is an integral over the density of sub-halos, their mass function, velocity distribution and the dynamical age of the stream. The rate of visible gap creation is a function of the width of the stream. The available data for four streams: the NW stream of M31, the Pal~5 stream, the Orphan Stream and the Eastern Banded Structure, are compared to the LCDM predicted relation. We find a remarkably good agreement, although there remains much to be done to improve the quality of the result. The narrower streams require that there is a total population of order 10^5 sub-halos above 10^5 M_sun to create the gaps.

preprint2012arXiv

Super-luminous supernovae at redshifts of 2.05 and 3.90

A rare class of `super-luminous' supernovae that are about ten or more times more luminous at their peaks than other types of luminous supernovae has recently been found at low to intermediate redshifts. A small subset of these events have luminosities that evolve slowly and result in radiated energies of around 10^51 ergs or more. Therefore, they are likely examples of `pair-instability' or `pulsational pair-instability' supernovae with estimated progenitor masses of 100 - 250 times that of the Sun. These events are exceedingly rare at low redshift, but are expected to be more common at high redshift because the mass distribution of the earliest stars was probably skewed to high values. Here we report the detection of two super-luminous supernovae, at redshifts of 2.05 and 3.90, that have slowly evolving light curves. We estimate the rate of events at redshifts of 2 and 4 to be approximately ten times higher than the rate at low redshift. The extreme luminosities of super-luminous supernovae extend the redshift limit for supernova detection using present technology, previously 2.36, and provide a way of investigating the deaths of the first generation of stars to form after the Big Bang.

preprint2011arXiv

Dark Matter Sub-Halo Counts via Star Stream Crossings

Dark matter sub-halos create gaps in the stellar streams orbiting in the halos of galaxies. We evaluate the sub-halo stream crossing integral with the guidance of simulations to find that the linear rate of gap creation, R_gap, in a typical Cold Dark Matter (CDM) galactic halo at 100 kpc is R_gap = 0.0066 M_8^-0.35 kpc^-1 Gyr^-1$, where M_8 = M/10^8 M_sun, is the minimum mass halo that creates a visible gap. The relation can be recast entirely in terms of observables, as R_gap= 0.059 w^-0.85 kpc^-1 Gyr^-1, for w in kpc, normalized at 100 kpc. Using published data, the density of gaps is estimated for M31's NW stream and the Milky Way Pal 5 stream, Orphan stream, and Eastern Banded Structure. The estimated rates of gap creation all have errors of 50% or more due to uncertain dynamical ages and the relatively noisy stream density measurements. The gap rate-width data are in good agreement with the CDM predicted relation. The high density of gaps in the narrow streams require a total halo population of 10^5 sub-halos above a minimum mass of 10^5 M_sun.

preprint2011arXiv

Red Nuggets at High Redshift: Structural Evolution of Quiescent Galaxies Over 10 Gyr of Cosmic History

We present an analysis of the size growth seen in early-type galaxies over 10 Gyr of cosmic time. Our analysis is based on a homogeneous synthesis of published data from 17 spectroscopic surveys observed at similar spatial resolution, augmented by new measurements for galaxies in the Gemini Deep Deep Survey. In total, our sample contains structural data for 465 galaxies (mainly early-type) in the redshift range 0.2<z<2.7. The size evolution of passively-evolving galaxies over this redshift range is gradual and continuous, with no evidence for an end or change to the process around z~1, as has been hinted at by some surveys which analyze subsets of the data in isolation. The size growth appears to be independent of stellar mass, with the mass-normalized half-light radius scaling with redshift as R_e (1+z)^(-1.62 +/- 0.34). Surprisingly, this power law seems to be in good agreement with the recently reported continuous size evolution of UV-bright galaxies in the redshift range z~0.5-3.5. It is also in accordance with the predictions from recent theoretical models.

preprint2010arXiv

The CFHTLS Deep Catalog of Interacting Galaxies I. Merger Rate Evolution to z=1.2

We present the rest-frame optical galaxy merger fraction between 0.2<z<1.2, as a function of stellar mass and optical luminosity, as observed by the Canada-France-Hawaii Telescope Legacy Deep Survey (CFHTLS-Deep). We developed a new classification scheme to identify major galaxy-galaxy mergers based on the presence of tidal tails and bridges. These morphological features are signposts of recent and ongoing merger activity. Through the visual classification of all galaxies, down to i_vega<22.2 (~27,000 galaxies) over 2 square degrees, we have compiled the CFHTLS Deep Catalog of Interacting Galaxies, with ~1600 merging galaxies. We find the merger fraction to be 4.3% +/-0.3% at z~0.3 and 19.0% +/-2.5% at z~1, implying evolution of the merger fraction going as (1+z)^m, with m=2.25 +/-0.24. This result is inconsistent with a mild or non-evolving (m<1.5) scenario at a >4sigma level of confidence. A mild trend, where massive galaxies with M>10^10.7 M_sun are undergoing fewer mergers than less massive systems M~10^10 M_sun), consistent with the expectations of galaxy assembly downsizing is observed. Our results also show that interacting galaxies have on average SFRs double that found in non-interacting field galaxies. We conclude that (1) the optical galaxy merger fraction does evolve with redshift, (2) the merger fraction depends mildly on stellar mass, with lower mass galaxies having higher merger fractions at z<1, and (3) star formation is triggered at all phases of a merger, with larger enhancements at later stages, consistent with N-body simulations.