-
The Effect of Age and Stellar Model Choice on Globular Cluster Color-to-Metallicity Conversions
Authors:
Kate Hartman,
William E. Harris
Abstract:
The photometric colors of globular clusters (GCs) act as effective proxies for metallicity, since all normally used optical/IR color indices exhibit a nonlinear but monotonic relation between their integrated color and their metallicity. One color index, (g - z) or (F475W - F850LP), has been spectroscopically calibrated in several studies, providing leverage to define color-to-metallicity conversi…
▽ More
The photometric colors of globular clusters (GCs) act as effective proxies for metallicity, since all normally used optical/IR color indices exhibit a nonlinear but monotonic relation between their integrated color and their metallicity. One color index, (g - z) or (F475W - F850LP), has been spectroscopically calibrated in several studies, providing leverage to define color-to-metallicity conversions for other indices. In this paper, building on the work of arXiv:2307.01863, we study the GC color-metallicity relation in more detail by testing the dependence of the relations on different suites of stellar models and different assumed GC ages. Though noticeable differences between models exist, we find that the net effect on the derived GCS metallicity distributions is small.
△ Less
Submitted 13 June, 2024;
originally announced June 2024.
-
The turbulent life of NGC 4696 as told by its globular cluster system
Authors:
S. Federle,
M. Gómez,
S. Mieske,
W. E. Harris,
M. Hilker,
I. A. Yegorova,
G. L. H. Harris
Abstract:
In this work we perform the photometric analysis of the globular cluster system (GCS) of the giant elliptical NGC4696, which is the brightest member of Centaurus, a rich and dynamically young galaxy cluster. We obtained deep Magellan 6.5 m/MegaCam (g', r', i') photometry, with which we identified a sample of 3818 stellar clusters around NGC4696 that were analyzed in the context of possible interac…
▽ More
In this work we perform the photometric analysis of the globular cluster system (GCS) of the giant elliptical NGC4696, which is the brightest member of Centaurus, a rich and dynamically young galaxy cluster. We obtained deep Magellan 6.5 m/MegaCam (g', r', i') photometry, with which we identified a sample of 3818 stellar clusters around NGC4696 that were analyzed in the context of possible interactions and its assembly history. After carefully modeling and subtracting the galaxy light, we used selection criteria based on the shape, colors, and magnitudes to identify GC candidates. We find a number of features that indicate a disturbed GCS that points toward a complex evolution with other neighboring members of Centaurus. Formally, two subpopulations could be found at (g'-i')_0 = 0.763 $\pm$ 0.004 and (g'-i')_0=1.012 $\pm$ 0.004. Moreover, the color distribution does not show the presence of a significant blue tilt, but it presents a trend with the radius, where at small galactocentric distances a unimodal distribution is preferable to a bimodal one, suggesting the presence of an intermediate GC population. Besides the color distribution, the metallicity distribution also shows a bimodal trend, with peaks at [Fe/H]=-1.363 $\pm$ 0.010 and [Fe/H]=-0.488 $\pm$ 0.012. The radial density profiles show different slopes for the blue and red populations and the azimuthal distributions are well fitted by an asymmetrical sinusoidal function, with peaks projecting toward two nearby galaxies, NGC4696B and NGC4709, indicating past interactions among these three galaxies. Finally, we derived a GC specific frequency of S_N=6.8 $\pm$ 0.9, in good agreement with the values obtained for other giant ellipticals and with previously estimated S_N of NGC4696. All these results point toward a complex GCS, strongly influenced by the interaction history of NGC4696 with the other galaxies of the Centaurus cluster.
△ Less
Submitted 12 June, 2024;
originally announced June 2024.
-
Utilizing Voronoi Tessellations to Determine Radial Density Profiles
Authors:
Veronika Dornan,
William E. Harris
Abstract:
We have developed a novel method of determining 2D radial density profiles for astronomical systems of discrete objects using Voronoi tessellations. This Voronoi-based method was tested against the standard annulus-based method on 5 simulated systems of objects, following known Hubble density profiles of varying parameters and sizes. It was found that the Voronoi-based method returned radial densi…
▽ More
We have developed a novel method of determining 2D radial density profiles for astronomical systems of discrete objects using Voronoi tessellations. This Voronoi-based method was tested against the standard annulus-based method on 5 simulated systems of objects, following known Hubble density profiles of varying parameters and sizes. It was found that the Voronoi-based method returned radial density fits with lower uncertainties on the fitting parameters across all 5 systems compared to the annulus-based method. The Voronoi-based method also consistently returned more accurate estimates of the total number of objects in each system than the annulus-based method, and this accuracy increased with increasing system size. Finally, the Voronoi-based method was applied to two observed globular cluster systems around brightest cluster galaxies ESO 444-G046 and 2MASX J13272961-3123237 and the results were compared to previous results for these galaxies obtained with the annulus-based method. Again, it was found that the Voronoi-based method returned fits with lower uncertainties on the fitting parameters, and the total number of globular clusters returned are within errors of the annulus-based method estimates, however also with lower uncertainties.
△ Less
Submitted 4 June, 2024;
originally announced June 2024.
-
Photometric Completeness Modelled With Neural Networks
Authors:
William E. Harris,
Joshua S. Speagle
Abstract:
In almost any study involving optical/NIR photometry, understanding the completeness of detection and recovery is an essential part of the work. The recovery fraction is, in general, a function of several variables including magnitude, color, background sky noise, and crowding. We explore how completeness can be modelled, {with the use of artificial-star tests,} in a way that includes all of these…
▽ More
In almost any study involving optical/NIR photometry, understanding the completeness of detection and recovery is an essential part of the work. The recovery fraction is, in general, a function of several variables including magnitude, color, background sky noise, and crowding. We explore how completeness can be modelled, {with the use of artificial-star tests,} in a way that includes all of these parameters \emph{simultaneously} within a neural network (NN) framework. The method is able to manage common issues including asymmetric completeness functions and the bilinear dependence of the detection limit on color index. We test the method with two sample HST (Hubble Space Telescope) datasets: the first involves photometry of the star cluster population around the giant Perseus galaxy NGC 1275, and the second involves the halo-star population in the nearby elliptical galaxy NGC 3377. The NN-based method achieves a classification accuracy of $>$\,94\%, and produces results entirely consistent with more traditional techniques for determining completeness. Additional advantages of the method are that none of the issues arising from binning of the data are present, and that a recovery probability can be assigned to every individual star in the real photometry. Our data, models, and code (called COINTOSS) can be found online on Zenodo at the following link: https://doi.org/10.5281/zenodo.8306488.
△ Less
Submitted 29 May, 2024;
originally announced May 2024.
-
Deep HST/UVIS imaging of the candidate dark galaxy CDG-1
Authors:
Pieter van Dokkum,
Dayi David Li,
Roberto Abraham,
Shany Danieli,
Gwendolyn M. Eadie,
William E. Harris,
Aaron J. Romanowsky
Abstract:
CDG-1 is a tight grouping of four likely globular clusters in the Perseus cluster, and a candidate dark galaxy with little or no diffuse light. Here we provide new constraints on the luminosity of any underlying stellar emission, using HST/UVIS F200LP imaging. No diffuse emission is detected, with a 2$σ$ upper limit of F200LP>28.1 mag/arcsec$^2$ on the 5'' scale of CDG-1. This surface brightness l…
▽ More
CDG-1 is a tight grouping of four likely globular clusters in the Perseus cluster, and a candidate dark galaxy with little or no diffuse light. Here we provide new constraints on the luminosity of any underlying stellar emission, using HST/UVIS F200LP imaging. No diffuse emission is detected, with a 2$σ$ upper limit of F200LP>28.1 mag/arcsec$^2$ on the 5'' scale of CDG-1. This surface brightness limit corresponds to a 2$σ$ lower limit of >0.5 for the fraction of the total luminosity that is in the form of globular clusters. The most likely alternative, although improbable, is that CDG-1 is a chance grouping of four globular clusters in the halo of the Perseus galaxy IC312.
△ Less
Submitted 21 May, 2024;
originally announced May 2024.
-
JWST Photometry of Globular Clusters in Abell 2744. II: luminosity and color distributions
Authors:
William E. Harris,
Marta Reina-Campos
Abstract:
Deep JWST imaging of the giant galaxy cluster Abell 2744, at redshift $z=0.308$, is used to explore the features of its rich population of globular clusters (GCs), building on our initial survey of the system (Harris & Reina-Campos 2023). We use the photometry of more than $10,000$ GCs over a three-magnitude range to construct the GC luminosity function (GCLF) and color distribution (CDF). These r…
▽ More
Deep JWST imaging of the giant galaxy cluster Abell 2744, at redshift $z=0.308$, is used to explore the features of its rich population of globular clusters (GCs), building on our initial survey of the system (Harris & Reina-Campos 2023). We use the photometry of more than $10,000$ GCs over a three-magnitude range to construct the GC luminosity function (GCLF) and color distribution (CDF). These results now specifically account for photometric incompleteness as a function of location relative to the five giant galaxies that dominate the gravitational potential of A2744. The total GC population in A2744 is estimated at $N_{\rm GC} \geq 1.1 \times 10^5$, consistent with its high total mass. We also directly compare the observed distributions with theoretical predictions for GC populations drawn from the recent EMP-Pathfinder simulations (Reina-Campos et al. 2022), viewed at the same 3.5 Gyr lookback time as the cluster. The simulations match the observations well, with the only notable disagreement being that the simulations predict larger numbers of GCs at high metallicity formed after $z\simeq2$ than are seen in the data.
△ Less
Submitted 16 April, 2024;
originally announced April 2024.
-
Imaging dark matter at the smallest scales with $z\approx1$ lensed stars
Authors:
J. M. Diego,
Sung Kei Li,
Alfred Amruth,
Ashish K. Meena,
Tom J. Broadhurst,
Patrick L. Kelly,
Alexei V. Filippenko,
Liliya L. R. Williams,
Adi Zitrin,
William E. Harris,
Marta Reina-Campos,
Carlo Giocoli,
Liang Dai,
Mitchell F. Struble,
Tommaso Treu,
Yoshinobu Fudamoto,
Daniel Gilman,
Anton M. Koekemoer,
Jeremy Lim,
J. M. Palencia,
Fengwu Sun,
Rogier A. Windhorst
Abstract:
Observations of caustic-crossing galaxies at redshift $0.7<z<1$ show a wealth of transient events. Most of them are believed to be microlensing events of highly magnified stars. Earlier work predicted such events should be common near the critical curves (CCs) of galaxy clusters, but some are found relatively far away from these CCs. We consider the possibility that substructure on milliarcsecond…
▽ More
Observations of caustic-crossing galaxies at redshift $0.7<z<1$ show a wealth of transient events. Most of them are believed to be microlensing events of highly magnified stars. Earlier work predicted such events should be common near the critical curves (CCs) of galaxy clusters, but some are found relatively far away from these CCs. We consider the possibility that substructure on milliarcsecond scales (few parsecs in the lens plane) is boosting the microlensing signal. We study the combined magnification from the macrolens, millilenses, and microlenses (3M-lensing). After considering realistic populations of millilenses and microlenses, we conclude that the enhanced microlensing rate around millilenses is not sufficient to explain the high fraction of observed events in the far region. Instead we find a that the shape of the luminosity function (LF) of the lensed stars combined with the amount of substructure in the lens plane determines the number of mcirolensing events found near and far from the CC. By measuring $β$ (the exponent of the LF), and the number density of microlensing events at each location, one can create a pseudoimage of the underlying distribution of mass on small scales. We identify two regimes: (i) positive imaging regime where $β>2$ and the number density of events is greater around substructureand the number density of events is greater around substructures, and (ii) negative imaging regime where $β<2$. We study the particular case of seven microlensing events found by HST in the Dragon arc (at z=0.725). We find that a population of supergiant stars with a steep LF with $β=2.55$ fits the distribution of these events. We identify a small region of high density of microlensing events, and interpret it as evidence of a possible invisible substructure, for which we derive a mass of $\sim 1.3 \times 10^8\,\Msun$ (within its Einstein radius).
△ Less
Submitted 22 April, 2024; v1 submitted 11 April, 2024;
originally announced April 2024.
-
Polarization of circumstellar debris disk light echoes
Authors:
Austin J. King,
Benjamin C. Bromley,
Preston W. Harris,
Scott J. Kenyon
Abstract:
Light echoes of debris disks around active stars can reveal disk structure and composition even when disks are not spatially resolved. Unfortunately, distinguishing reflected light from quiescent starlight and unexpected post-peak flare structure is challenging, especially for edge-on geometries where the time delay between observed flare photons and light scattered from the near side of the disk…
▽ More
Light echoes of debris disks around active stars can reveal disk structure and composition even when disks are not spatially resolved. Unfortunately, distinguishing reflected light from quiescent starlight and unexpected post-peak flare structure is challenging, especially for edge-on geometries where the time delay between observed flare photons and light scattered from the near side of the disk is short. Here, we take advantage of the fact that scattered light from a dusty disk is polarized, depending on the location of the scattering site and the orientation of the disk relative to a distant observer. Filtering reflected light into its polarized components allows echoes to stand out in predictable ways. We test this idea with a simple model for a disk around an active M dwarf. Our results demonstrate that the use of polarimetric data of flaring stars can significantly enhance echo signals relative to starlight and yield more robust and accurate fits to disk parameters compared to analyses based on the total intensity alone.
△ Less
Submitted 9 February, 2024;
originally announced February 2024.
-
RESCUER: Cosmological K-corrections for star clusters
Authors:
Marta Reina-Campos,
William E. Harris
Abstract:
The advent of JWST (the James Webb Space Telescope) now allows entire star cluster populations to be imaged in galaxies at cosmologically significant redshifts, bringing with it the need to apply K-corrections to their magnitudes and colour indices. Since the stellar populations within star clusters can be well approximated by a single age and metallicity, their spectral energy distributions are v…
▽ More
The advent of JWST (the James Webb Space Telescope) now allows entire star cluster populations to be imaged in galaxies at cosmologically significant redshifts, bringing with it the need to apply K-corrections to their magnitudes and colour indices. Since the stellar populations within star clusters can be well approximated by a single age and metallicity, their spectral energy distributions are very different from those of galaxies or supernovae, and their K-corrections behave differently. We derive the photometric K-corrections versus redshift for model star clusters that cover a wide range of ages and metallicities, illustrating the results particularly for the broadband filters on the HST/ACS and the JWST/NIRCam cameras that are most commonly being used for imaging of populations of star clusters in distant galaxies. In an Appendix, we introduce a simple webtool called RESCUER that can generate K-values for any user-defined combination of cluster properties.
△ Less
Submitted 3 October, 2023;
originally announced October 2023.
-
JWST Photometry of Globular Cluster Populations in Abell 2744 at $z=0.3$
Authors:
William E. Harris,
Marta Reina-Campos
Abstract:
JWST imaging of the rich galaxy cluster Abell 2744 at $z=0.308$ has been used by the UNCOVER team (Bezanson et al. 2022) to construct mosaic images in the NIRCAM filters. The exceptionally deep images in the ($F115W$, $F150W$, $F200W$) bands reveal a large population of unresolved pointlike sources across the field, the vast majority of which are globular clusters (GCs). To the limits of our photo…
▽ More
JWST imaging of the rich galaxy cluster Abell 2744 at $z=0.308$ has been used by the UNCOVER team (Bezanson et al. 2022) to construct mosaic images in the NIRCAM filters. The exceptionally deep images in the ($F115W$, $F150W$, $F200W$) bands reveal a large population of unresolved pointlike sources across the field, the vast majority of which are globular clusters (GCs). To the limits of our photometry, more than 10,000 such objects were measured, most of which are in the halos of the five largest A2744 galaxies but which also include GCs around some satellite galaxies and throughout the IntraCluster Medium. Their luminosity function follows a lognormal shape, with the data reaching to within one magnitude of the classic GCLF turnover point. The colour index ($F115W-F200W$) in particular covers a range of $0.5$ mag, clearly resolving the expected internal spread of GC metallicities. The estimated GC masses are systematically higher than in present-day galaxies, consistent with a large, normal GC population seen at a $3.5~$Gyr earlier stage of dynamical evolution. Lastly, the spatial distribution of the bluer (more metal-poor) GCs resembles the gravitational lensing map of the cluster, consistent with recent theoretical suggestions.
△ Less
Submitted 22 September, 2023; v1 submitted 26 July, 2023;
originally announced July 2023.
-
Comparing Globular Cluster System Properties with Host Galaxy Environment
Authors:
Kate Hartman,
William E. Harris,
John P. Blakeslee,
Chung-Pei Ma,
Jenny E. Greene
Abstract:
We present Hubble Space Telescope photometry in optical (F475X) and near-infrared (F110W) bands of the globular cluster (GC) systems of the inner halos of a sample of 15 massive elliptical galaxies. The targets are selected from the volume-limited MASSIVE survey, and chosen to sample a range of environments from sparsely populated groups to BCGs in dense clusters. We also present a quantitative mo…
▽ More
We present Hubble Space Telescope photometry in optical (F475X) and near-infrared (F110W) bands of the globular cluster (GC) systems of the inner halos of a sample of 15 massive elliptical galaxies. The targets are selected from the volume-limited MASSIVE survey, and chosen to sample a range of environments from sparsely populated groups to BCGs in dense clusters. We also present a quantitative model of the relation between (F475X - F110W) colour and cluster metallicity [M/H], using simulated GCs. Because much of the GC population in such galaxies is built up through accretion, the metallicity distribution of the GC systems might be expected to vary with galaxy environment. The photometry is used to create a completeness-corrected metallicity distribution for each galaxy in the sample, and to fit a double Gaussian curve to each histogram in order to model the two standard red and blue subpopulations. Finally, the properties of the GC metallicity distribution are correlated against galaxy environment. We find that almost no GCS properties and host galaxy environmental properties are correlated, with the exception of a weak but consistent correlation between blue fraction and nth-nearest neighbour surface density. The results suggest that the systemic properties of the GCS, at least in the inner to mid-halo regions, are influenced more strongly by the local environment at early times, rather than by the environmental properties we see today.
△ Less
Submitted 4 July, 2023;
originally announced July 2023.
-
The HERBAL model: A hierarchical errors-in-variables Bayesian lognormal hurdle model for galactic globular cluster populations
Authors:
Samantha C. Berek,
Gwendolyn M. Eadie,
Joshua S. Speagle,
William E. Harris
Abstract:
Galaxy stellar mass is known to be monotonically related to the size of the galaxy's globular cluster (GC) population for Milky Way sized and larger galaxies. However, the relation becomes ambiguous for dwarf galaxies, where there is some evidence for a downturn in GC population size at low galaxy masses. Smaller dwarfs are increasingly likely to have no GCs, and these zeros cannot be easily incor…
▽ More
Galaxy stellar mass is known to be monotonically related to the size of the galaxy's globular cluster (GC) population for Milky Way sized and larger galaxies. However, the relation becomes ambiguous for dwarf galaxies, where there is some evidence for a downturn in GC population size at low galaxy masses. Smaller dwarfs are increasingly likely to have no GCs, and these zeros cannot be easily incorporated into linear models. We introduce the Hierarchical ERrors-in-variables Bayesian lognormAL hurdle (HERBAL) model to represent the relationship between dwarf galaxies and their GC populations, and apply it to the sample of Local Group galaxies where the luminosity range coverage is maximal. This bimodal model accurately represents the two populations of dwarf galaxies: those that have GCs and those that do not. Our model thoroughly accounts for all uncertainties, including measurement uncertainty, uncertainty in luminosity to stellar mass conversions, and intrinsic scatter. The hierarchical nature of our Bayesian model also allows us to estimate galaxy masses and individual mass-to-light ratios from luminosity data within the model. We find that 50% of galaxies are expected to host globular cluster populations at a stellar mass of $\log_{10}(M_*)=6.996$, and that the expected mass of GC populations remains linear down to the smallest galaxies. Our hierarchical model recovers an accurate estimate of the Milky Way stellar mass. Under our assumed error model, we find a non-zero intrinsic scatter of $0.59_{-0.21}^{+0.30}$ (95% credible interval) that should be accounted for in future models.
△ Less
Submitted 22 September, 2023; v1 submitted 26 June, 2023;
originally announced June 2023.
-
Investigating the $M_{GCS}-M_h$ Relation in the Most Massive Galaxies
Authors:
Veronika Dornan,
William E. Harris
Abstract:
The relation between the total mass contained in the globular clusters of a galaxy and the mass of its dark matter halo has been found observationally to be nearly linear over five decades of mass. However, the high-mass end of this relation is not well determined from previous data and shows large scatter. We analyze the globular cluster systems (GCSs) of a homogeneous sample of 11 brightest clus…
▽ More
The relation between the total mass contained in the globular clusters of a galaxy and the mass of its dark matter halo has been found observationally to be nearly linear over five decades of mass. However, the high-mass end of this relation is not well determined from previous data and shows large scatter. We analyze the globular cluster systems (GCSs) of a homogeneous sample of 11 brightest cluster galaxies (BCGs) through DOLPHOT photometry of their deep Hubble Space Telescope (HST) images in the F814W filter. We standardize the definition of $M_{GCS}$, the total GCS mass, by using the GC total population within a limiting radius of $0.1 R_{virial}$, while the dark-matter halo mass $M_h$ is determined from the weak-lensing calibration of $M_h$ versus $M_{bary}$. When these 11 BCGs are added to the previously studied homogeneous catalogue of Virgo member galaxies, a total value for $η= M_{GCS}/M_h$ is found to be $(3.0\pm1.8_{internal})\times10^{-5}$, slightly higher than previous estimates but with much reduced uncertainty. Perhaps more importantly, the results suggest that the relation continues to have a near-linear shape at the highest galaxy masses, strongly reinforcing the conclusion that accreted GCs make a major contribution to the GC populations at high galaxy mass.
△ Less
Submitted 21 April, 2023;
originally announced April 2023.
-
Ice, Ice, Maybe? Investigating 46P/Wirtanen's Inner Coma For Icy Grains
Authors:
Theodore Kareta,
John W. Noonan,
Walter M. Harris,
Alessondra Springmann
Abstract:
The release of volatiles from comets is usually from direct sublimation of ices on the nucleus, but for very or hyper-active comets other sources have to be considered to account for the total production rates. In this work, we present new near-infrared imaging and spectroscopic observations of 46P/Wirtanen taken during its close approach to the Earth on 2018 December 19 with the MMIRS instrument…
▽ More
The release of volatiles from comets is usually from direct sublimation of ices on the nucleus, but for very or hyper-active comets other sources have to be considered to account for the total production rates. In this work, we present new near-infrared imaging and spectroscopic observations of 46P/Wirtanen taken during its close approach to the Earth on 2018 December 19 with the MMIRS instrument at the MMT Observatory to search for signatures of icy or ice-rich grains in its inner coma that might explain its previously reported excess water production. The morphology of the images does not suggest any change in grain properties within the field of view, and the NIR spectra do not show the characteristic absorption features of water ice. Using a new MCMC-based implementation of the spectral modeling approach of Protopapa et al. (2018), we estimate the areal water ice fraction of the coma to be less than 0.6%. When combined with slit-corrected Afrho values for the J, H, and K bands and previously measured dust velocities for this comet, we estimate an icy grain production rate of less than 4.6 kg/s. This places a strict constraint on the water production rate from pure icy grains in the coma, and in turn we find that for the 2018-2019 apparition approximately 64% of 46P's surface was sublimating water near perihelion. WE then discuss 46P's modern properties within the context of other (formerly) hyper-active comets to understand how these complex objects evolve.
△ Less
Submitted 12 April, 2023;
originally announced April 2023.
-
Early Evolution and 3D Structure of Embedded Star Clusters
Authors:
Claude Cournoyer-Cloutier,
Alison Sills,
William E. Harris,
Sabrina M. Appel,
Sean C. Lewis,
Brooke Polak,
Aaron Tran,
Martijn J. C. Wilhelm,
Mordecai-Mark Mac Low,
Stephen L. W. McMillan,
Simon Portegies Zwart
Abstract:
We perform simulations of star cluster formation to investigate the morphological evolution of embedded star clusters in the earliest stages of their evolution. We conduct our simulations with Torch, which uses the AMUSE framework to couple state-of-the-art stellar dynamics to star formation, radiation, stellar winds, and hydrodynamics in FLASH. We simulate a suite of $10^4$ M$_{\odot}$ clouds at…
▽ More
We perform simulations of star cluster formation to investigate the morphological evolution of embedded star clusters in the earliest stages of their evolution. We conduct our simulations with Torch, which uses the AMUSE framework to couple state-of-the-art stellar dynamics to star formation, radiation, stellar winds, and hydrodynamics in FLASH. We simulate a suite of $10^4$ M$_{\odot}$ clouds at 0.0683 pc resolution for $\sim$ 2 Myr after the onset of star formation, with virial parameters $α$ = 0.8, 2.0, 4.0 and different random samplings of the stellar initial mass function and prescriptions for primordial binaries. Our simulations result in a population of embedded clusters with realistic morphologies (sizes, densities, and ellipticities) that reproduce the known trend of clouds with higher initial $α$ having lower star formation efficiencies. Our key results are as follows: (1) Cluster mass growth is not monotonic, and clusters can lose up to half of their mass while they are embedded. (2) Cluster morphology is not correlated with cluster mass and changes over $\sim$ 0.01 Myr timescales. (3) The morphology of an embedded cluster is not indicative of its long-term evolution but only of its recent history: radius and ellipticity increase sharply when a cluster accretes stars. (4) The dynamical evolution of very young embedded clusters with masses $\lesssim$ 1000 M$_{\odot}$ is dominated by the overall gravitational potential of the star-forming region rather than by internal dynamical processes such as two- or few-body relaxation.
△ Less
Submitted 16 February, 2023;
originally announced February 2023.
-
A Photometric Survey of Globular Cluster Systems in Brightest Cluster Galaxies
Authors:
William E. Harris
Abstract:
Hubble Space Telescope imaging for 26 giant early-type galaxies, all drawn from the MAST archive, is used to carry out photometry of their surrounding globular cluster (GC) systems. Most of these targets are Brightest Cluster Galaxies (BCGs) and their distances range from 24 to 210 Mpc. The catalogs of photometry, completed with DOLPHOT, are publicly available. The GC color indices are converted t…
▽ More
Hubble Space Telescope imaging for 26 giant early-type galaxies, all drawn from the MAST archive, is used to carry out photometry of their surrounding globular cluster (GC) systems. Most of these targets are Brightest Cluster Galaxies (BCGs) and their distances range from 24 to 210 Mpc. The catalogs of photometry, completed with DOLPHOT, are publicly available. The GC color indices are converted to [Fe/H] through a combination of 12-Gyr SSP (Single Stellar Population) models and direct spectroscopic calibration of the fiducial color index (F475W-F850LP). All the resulting metallicity distribution functions (MDFs) can be accurately matched by bimodal Gaussian functions. The GC systems in all the galaxies also exhibit shallow metallicity gradients with projected galactocentric distance that average $Z \sim R_{gc}^{-0.3}$. Several parameters of the MDFs including the means, dispersions, and blue/red fractions are summarized. Perhaps the most interesting new result is the trend of blue/red GC fraction with galaxy mass, which connects with predictions from recent simulations of GC formation within hierarchical assembly of large galaxies. The observed trend reveals two major transition stages: for low-mass galaxies, the metal-rich (red) GC fraction $f(red)$ increases steadily with galaxy mass, until halo mass $M_h \simeq 3 \times 10^{12} M_{\odot}$. Above this point, more than half the metal-poor (blue) GCs come from accreted satellites and $f(red)$ starts declining. But above a still higher transition point near $M_h \simeq 10^{14} M_{\odot}$, the data hint that $f(red)$ may start to increase again because the metal-rich GCs also become dominated by accreted systems.
△ Less
Submitted 12 December, 2022;
originally announced December 2022.
-
29P/Schwassmann-Wachmann: A Rosetta Stone for Amorphous Water Ice and CO <-> CO2 Conversion in Centaurs and Comets?
Authors:
C. M. Lisse,
J. K. Steckloff,
D. Prialnik,
M. Womack,
O. Harrington-Pinto,
G. Sarid,
Y. R. Fernandez,
C. A. Schambeau,
T. Kareta,
N. H. Samarasinha,
W. Harris,
K. Volk,
L. M. Woodney,
D. P. Cruikshank,
S. A. Sandford
Abstract:
Centaur 29P/Schwassmann-Wachmann 1 (SW1) is a highly active object orbiting in the transitional Gateway region (Sarid et al. 2019) between the Centaur and Jupiter Family Comet regions. SW1 is unique among the Centaurs in that it experiences quasi-regular major outbursts and produces CO emission continuously; however, the source of the CO is unclear. We argue that due to its very large size (approx…
▽ More
Centaur 29P/Schwassmann-Wachmann 1 (SW1) is a highly active object orbiting in the transitional Gateway region (Sarid et al. 2019) between the Centaur and Jupiter Family Comet regions. SW1 is unique among the Centaurs in that it experiences quasi-regular major outbursts and produces CO emission continuously; however, the source of the CO is unclear. We argue that due to its very large size (approx. 32 km radius), SW1 is likely still responding, via amorphous water ice (AWI) conversion to crystalline water ice (CWI), to the rapid change in its external thermal environment produced by its dynamical migration from the Kuiper belt to the Gateway Region at the inner edge of the Centaur region at 6 au. It is this conversion process that is the source of the abundant CO and dust released from the object during its quiescent and outburst phases. If correct, these arguments have a number of important predictions testable via remote sensing and in situ spacecraft characterization, including: the quick release on Myr timescales of CO from AWI conversion for any few km-scale scattered disk KBO transiting into the inner system; that to date SW1 has only converted between 50 to 65% of its nuclear AWI to CWI; that volume changes upon AWI conversion could have caused subsidence and cave-ins, but not significant mass wasting or crater loss on SW1; that SW1s coma should contain abundant amounts of CWI CO2-rich icy dust particles; and that when SW1 transits into the inner system within the next 10,000 years, it will be a very different kind of JFC comet.
△ Less
Submitted 27 October, 2022; v1 submitted 19 September, 2022;
originally announced September 2022.
-
Light from the Darkness: Detecting Ultra-Diffuse Galaxies in the Perseus Cluster through Over-densities of Globular Clusters with a Log-Gaussian Cox Process
Authors:
Dayi Li,
Gwendolyn M. Eadie,
Roberto G. Abraham,
Patrick E. Brown,
William E. Harris,
Steven R. Janssens,
Aaron J. Romanowsky,
Pieter van Dokkum,
Shany Danieli
Abstract:
We introduce a new method for detecting ultra-diffuse galaxies by searching for over-densities in intergalactic globular cluster populations. Our approach is based on an application of the log-Gaussian Cox process, which is a commonly used model in the spatial statistics literature but rarely used in astronomy. This method is applied to the globular cluster data obtained from the PIPER survey, a \…
▽ More
We introduce a new method for detecting ultra-diffuse galaxies by searching for over-densities in intergalactic globular cluster populations. Our approach is based on an application of the log-Gaussian Cox process, which is a commonly used model in the spatial statistics literature but rarely used in astronomy. This method is applied to the globular cluster data obtained from the PIPER survey, a \textit{Hubble Space Telescope} imaging program targeting the Perseus cluster. We successfully detect all confirmed ultra-diffuse galaxies with known globular cluster populations in the survey. We also identify a potential galaxy that has no detected diffuse stellar content. Preliminary analysis shows that it is unlikely to be merely an accidental clump of globular clusters or other objects. If confirmed, this system would be the first of its kind. Simulations are used to assess how the physical parameters of the globular cluster systems within ultra-diffuse galaxies affect their detectability using our method. We quantify the correlation of the detection probability with the total number of globular clusters in the galaxy and the anti-correlation with increasing half-number radius of the globular cluster system. The Sérsic index of the globular cluster distribution has little impact on detectability.
△ Less
Submitted 16 June, 2022; v1 submitted 11 April, 2022;
originally announced April 2022.
-
Near-Infrared Spectroscopy Of The Nucleus Of Low-Activity Comet P/2016 BA$_{14}$ During Its 2016 Close Approach
Authors:
Theodore Kareta,
Vishnu Reddy,
Juan A. Sanchez,
Walter M. Harris
Abstract:
The Near-Earth Comet P/2016 BA$_{14}$ (PanSTARRS) is a slow-rotatating nearly-dormant object, a likely dynamical twin of 252P/LINEAR, and was recently shown to have a mid-infrared spectrum very dissimilar to other comets. BA$_{14}$ also recently selected one of the back-up targets for the ESA's \textit{Comet Interceptor}, so a clearer understanding of BA$_{14}$'s modern properties would not just i…
▽ More
The Near-Earth Comet P/2016 BA$_{14}$ (PanSTARRS) is a slow-rotatating nearly-dormant object, a likely dynamical twin of 252P/LINEAR, and was recently shown to have a mid-infrared spectrum very dissimilar to other comets. BA$_{14}$ also recently selected one of the back-up targets for the ESA's \textit{Comet Interceptor}, so a clearer understanding of BA$_{14}$'s modern properties would not just improve our understanding of how comets go dormant, but could also aid planning for a potential spacecraft visit. We present observations of BA$_{14}$ taken on two dates during its 2016 Earth close approach with the NASA Infrared Telescope Facility, both of which are consistent with direct observations of its nucleus. The reflectance spectrum of BA$_{14}$ is similar to 67P/Churyumov-Gerasimenko, albeit highly phase-reddened. Thermal emission contaminates the reflectance spectrum at longer wavelengths, which we correct with a new Markov Chain Monte Carlo thermal modeling code. The models suggest $BA_{14}$'s visible geometric albedo is $p_V=0.01-0.03$, consistent with radar observations, its beaming parameter is typical for NEOs observed in its geometry, and its reflectrance spectrum is red and linear throughout H and K band. It appears very much like a "normal" comet nucleus, despite its mid-infrared oddities. A slow loss of fine grains as the object's activity diminished might help to reconcile some of the lines of evidence, and we discuss other possibilities. A spacecraft flyby past BA$_{14}$ could get closer to the nucleus than with a more active target, and we highlight some science questions that could be addressed with a visit to a (nearly-)dormant comet.
△ Less
Submitted 3 April, 2022;
originally announced April 2022.
-
Operating Spacecraft Around Comets: Evaluation of the Near-Nucleus Environment
Authors:
C. M. Lisse,
M. R. Combi,
T. L. Farnham,
N. Dello Russo,
S. Sandford,
A. F. Cheng,
U. Fink,
W. M. Harris,
J. McMahon,
D. J. Scheeres,
H. A. Weaver,
J. Leary
Abstract:
We present a study of the current state of knowledge concerning spacecraft operations and potential hazards while operating near a comet nucleus. Starting from simple calculations comparing the cometary coma environment to benign conditions on Earth, we progress to sophisticated engineering models of spacecraft behavior, and then confront these models with recent spacecraft proximity operations ex…
▽ More
We present a study of the current state of knowledge concerning spacecraft operations and potential hazards while operating near a comet nucleus. Starting from simple calculations comparing the cometary coma environment to benign conditions on Earth, we progress to sophisticated engineering models of spacecraft behavior, and then confront these models with recent spacecraft proximity operations experience. Finally, we make recommendations from lessons learned for future spacecraft missions that enter into orbit around a comet for long-term operations. All of these considerations indicate that, with a proper spacecraft design and operations planning, the near-nucleus environment can be a relatively safe region in which to operate, even for an active short period comet near perihelion with gas production rates as high as 1e29 molecules/s. With gas densities similar to those found in good laboratory vacuums, dust densities similar to Class 100 cleanrooms, dust particle velocities of 10s of m/s, and microgravity forces that permit slow and deliberate operations, the conditions around a comet are generally more benign than a typical day on Mars. Even in strong dust jets near the nucleus surface, dust densities tend to be only a few grains/cm3, about the same as in a typical interior room on Earth. Stochastic forces on a modern spacecraft with tens of square meters of projected surface area can be accounted for using modern Attitude Control Systems to within tens of meters navigation error; surface contamination issues are only important for spacecraft spending months to years within a few kilometers of the nucleus surface; and the issues the Rosetta spacecraft faced, confusion of celestial star trackers by sunlit dust particles flying past the spacecraft, will be addressed using the next generation of star trackers implementing improved transient rejection algorithms.
△ Less
Submitted 26 January, 2022;
originally announced January 2022.
-
Radial Distribution of the Dust Comae of Comets 45P/Honda-Mrkos-Pajduusáková and 46P/Wirtanen
Authors:
C. Lejoly,
W. Harris,
N. Samarasinha,
B. E. A. Mueller,
E. Howell,
J. Bodnarik,
A. Springmann,
T. Kareta,
B. Sharkey,
J. Noonan,
L. R. Bedin,
J. -G. Bosch,
A. Brosio,
E. Bryssinck,
J. -B. de Vanssay,
F. -J. Hambsch,
O. Ivanova,
V. Krushinsky,
Z. -Y. Lin,
F. Manzini,
A. Maury,
N. Moriya,
P. Ochner,
V. Oldani
Abstract:
There was an unprecedented opportunity to study the inner dust coma environment, where the dust and gas are not entirely decoupled, of comets 45P/Honda-Mrkos-Pajduusáková (45P/HMP) from Dec. 26, 2016 - Mar. 15, 2017, and 46P/Wirtanen from Nov. 10, 2018 - Feb. 13, 2019, both in visible wavelengths. The radial profile slopes of these comets were measured in the R and HB-BC filters most representativ…
▽ More
There was an unprecedented opportunity to study the inner dust coma environment, where the dust and gas are not entirely decoupled, of comets 45P/Honda-Mrkos-Pajduusáková (45P/HMP) from Dec. 26, 2016 - Mar. 15, 2017, and 46P/Wirtanen from Nov. 10, 2018 - Feb. 13, 2019, both in visible wavelengths. The radial profile slopes of these comets were measured in the R and HB-BC filters most representative of dust, and deviations from a radially expanding coma were identified as significant. The azimuthally averaged radial profile slope of comet 45P/HMP gradually changes from -1.81 $\pm$ 0.20 at 5.24 days pre-perihelion to -0.35 $\pm$ 0.16 at 74.41 days post perihelion. Contrastingly, the radial profile slope of 46P/Wirtanen stays fairly constant over the observed time period at -1.05 $\pm$ 0.05. Additionally, we find that the radial profile of 46P/Wirtanen is azimuthally dependent on the skyplane-projected solar position angle, while that of 45P/HMP is not. These results suggest that comet 45P/HMP and 46P/Wirtanen have vastly different coma dust environments and that their dust properties are distinct. As evident from these two comets, well-resolved inner comae are vital for detailed characterization of dust environments.
△ Less
Submitted 20 December, 2021;
originally announced December 2021.
-
The outermost stellar halo of NGC 5128 (Centaurus A): Radial structure
Authors:
M. Rejkuba,
W. E. Harris,
L. Greggio,
D. Crnojević,
G. L. H. Harris
Abstract:
The extended stellar halos of galaxies contain important clues for investigating their assembly history and evolution. We investigate the resolved stellar content and the extended halo of NGC 5128 as a function of galactocentric distance. We used HST images to resolve individual red giant branch (RGB) stars in 28 independent pointings. Star counts from deep VI color-magnitude diagrams reaching at…
▽ More
The extended stellar halos of galaxies contain important clues for investigating their assembly history and evolution. We investigate the resolved stellar content and the extended halo of NGC 5128 as a function of galactocentric distance. We used HST images to resolve individual red giant branch (RGB) stars in 28 independent pointings. Star counts from deep VI color-magnitude diagrams reaching at least 1.5 mag below the tip of the RGB are used to derive the surface density distribution of the halo. The contamination by Milky Way stars is assessed with a new control field, with models, and by combining optical and near-IR photometry. We present a new calibration of the WFC3 F606W+F814W photometry to the ground-based VI photometric system. The photometry shows that the stellar halo of NGC 5128 is dominated by old RGB stars that are present in all fields. The V-band surface brightness changes from 23 to 32 mag/arcsec$^2$ between 8.3 kpc from the galaxy center to our outermost halo fields located 140 kpc away from the center along the major axis and 92 kpc along the minor axis. Within ~30 kpc, we also find evidence for a 2-3 Gyr old population traced by bright asymptotic giant branch stars. This population contributes only up to 10% in total stellar mass if it is 2 Gyr old, but a larger fraction of 30-40% is required if its age is 3 Gyr. The stellar surface density profile is well fit by a r$^{1/4}$ curve or a power-law $\sim r^{-3.1}$ over the full radial range, with no obvious break in the slope, but with large field-to-field scatter. The ellipticity measured from integrated-light photometry in the inner parts, $e=(b/a)=0.77$, flattens to $e=0.54 \pm 0.02$ beyond 30 kpc. Considering the flattening of the outer halo, the projection of the elliptical isophote on the semimajor axis for our most distant field reaches nearly 30 effective radii. [abridged]
△ Less
Submitted 4 November, 2021;
originally announced November 2021.
-
Clearing the hurdle: The mass of globular cluster systems as a function of host galaxy mass
Authors:
Gwendolyn M. Eadie,
William E. Harris,
Aaron Springford
Abstract:
Current observational evidence suggests that all large galaxies contain globular clusters (GCs), while the smallest galaxies do not. Over what galaxy mass range does the transition from GCs to no GCs occur? We investigate this question using galaxies in the Local Group, nearby dwarf galaxies, and galaxies in the Virgo Cluster Survey. We consider four types of statistical models: (1) logistic regre…
▽ More
Current observational evidence suggests that all large galaxies contain globular clusters (GCs), while the smallest galaxies do not. Over what galaxy mass range does the transition from GCs to no GCs occur? We investigate this question using galaxies in the Local Group, nearby dwarf galaxies, and galaxies in the Virgo Cluster Survey. We consider four types of statistical models: (1) logistic regression to model the probability that a galaxy of stellar mass $M_{\star}$ has any number of GCs; (2) Poisson regression to model the number of GCs versus $M_{\star}$, (3) linear regression to model the relation between GC system mass ($\log{M_{gcs}}$) and host galaxy mass ($\log{M_{\star}}$), and (4) a Bayesian lognormal hurdle model of the GC system mass as a function of galaxy stellar mass for the entire data sample. From the logistic regression, we find that the 50% probability point for a galaxy to contain GCs is $M_{\star}=10^{6.8}M_{\odot}$. From post-fit diagnostics, we find that Poisson regression is an inappropriate description of the data. Ultimately, we find that the Bayesian lognormal hurdle model, which is able to describe how the mass of the GC system varies with $M_{\star}$ even in the presence of many galaxies with no GCs, is the most appropriate model over the range of our data. In an Appendix, we also present photometry for the little-known GC in the Local Group dwarf Ursa Major II.
△ Less
Submitted 28 October, 2021;
originally announced October 2021.
-
Investigating the Relationship between (3200) Phaethon and (155140) 2005 UD through Telescopic and Laboratory Studies
Authors:
Theodore Kareta,
Vishnu Reddy,
Neil Pearson,
Juan A. Sanchez,
Walter M. Harris
Abstract:
The relationship between the Near-Earth Objects (3200) Phaethon and (155140) 2005 UD is unclear. While both are parents to Meteor Showers, (the Geminids and Daytime Sextantids, respectively), have similar visible-wavelength reflectance spectra and orbits, dynamical investigations have failed to find any likely method to link the two objects in the recent past. Here we present the first near-infrar…
▽ More
The relationship between the Near-Earth Objects (3200) Phaethon and (155140) 2005 UD is unclear. While both are parents to Meteor Showers, (the Geminids and Daytime Sextantids, respectively), have similar visible-wavelength reflectance spectra and orbits, dynamical investigations have failed to find any likely method to link the two objects in the recent past. Here we present the first near-infrared reflectance spectrum of 2005 UD, which shows it to be consistently linear and red-sloped unlike Phaethon's very blue and concave spectrum. Searching for a process that could alter some common starting material to both of these end states, we hypothesized that the two objects had been heated to different extents, motivated by their near-Sun orbits, the composition of Geminid meteoroids, and previous models of Phaethon's surface. We thus set about building a new laboratory apparatus to acquire reflectance spectra of meteoritic samples after heating to higher temperatures than available in the literature to test this hypothesis and were loaned a sample of the CI Chondrite Orgueil from the Vatican Meteorite Collection for testing. We find that while Phaethon's spectrum shares many similarities with different CI Chondrites, 2005 UD's does not. We thus conclude that the most likely relationship between the two objects is that their similar properties are only by coincidence as opposed to a parent-fragment scenario, though the ultimate test will be when JAXA's DESTINY+ mission visits one or both objects later this decade. We also discuss possible paths forward to understanding Phaethon's properties from dynamical and compositional grounds.
△ Less
Submitted 2 September, 2021;
originally announced September 2021.
-
Repeating gas ejection events from comet 45P/Honda-Mrkos-Pajdušáková
Authors:
Alessondra Springmann,
Walter M. Harris,
Erin L. Ryan,
Cassandra Lejoly,
Ellen S. Howell,
Beatrice E. A. Mueller,
Nalin H. Samarasinha,
Laura M. Woodney,
Jordan K. Steckloff
Abstract:
Studying materials released from Jupiter-family comets (JFCs) -- as seen in their inner comæ, the envelope of gas and dust that forms as the comet approaches the Sun -- improves the understanding of their origin and evolutionary history. As part of a coordinated, multi-wavelength observing campaign, we observed comet 45P/Honda-Mrkos-Pajdušáková during its close approach to Earth in February 2017.…
▽ More
Studying materials released from Jupiter-family comets (JFCs) -- as seen in their inner comæ, the envelope of gas and dust that forms as the comet approaches the Sun -- improves the understanding of their origin and evolutionary history. As part of a coordinated, multi-wavelength observing campaign, we observed comet 45P/Honda-Mrkos-Pajdušáková during its close approach to Earth in February 2017. Narrowband observations were taken using the Bok 90" telescope at Kitt Peak National Observatory on February 16 and 17 UT, revealing gas and dust structures. We observed different jet directions for different volatile species, implying source region heterogeneity, consistent with other ground-based and \textit{in situ} observations of other comet nuclei. A repeating feature visible in CN and C$_2$ images on February 16 was also observed on February 17 with an interval of $7.6\pm0.1$ hours, consistent with the rotation period of the comet derived from Arecibo Observatory radar observations. The repeating feature's projected gas velocity away from the nucleus is 0.8 km s$^{-1}$, with an expansion velocity of 0.5 km s$^{-1}$. A bright compact spot adjacent to the nucleus provides a lower limit of the amount of material released in one cycle of $\sim$9.2 kg, depending on composition -- a quantity small enough to be produced by repeated exposure of nucleus ices to sunlight. This repeating CN jet, forming within 400 km of the nucleus, may be typical of inner coma behavior in JFCs; however, similar features could be obscured by other processes and daughter product species when viewed from distances further than the scale length of CN molecules.
△ Less
Submitted 28 November, 2021; v1 submitted 21 July, 2021;
originally announced July 2021.
-
Complex Water Ice Mixtures on NII Nereid: Constraints from NIR Reflectance
Authors:
Benjamin N. L. Sharkey,
Vishnu Reddy,
Juan A. Sanchez,
Matthew R. M. Izawa,
Walter M. Harris
Abstract:
Nereid, Neptune's third largest satellite, lies in an irregular orbit and is the only outer satellite in the system (apart from Triton) that can be spectroscopically characterized with the current generation of Earth-based telescopes. We report our results on spectral characterization of Nereid using its reflectance spectrum from 0.8-2.4 $μm$, providing the first measurements over the range of 0.8…
▽ More
Nereid, Neptune's third largest satellite, lies in an irregular orbit and is the only outer satellite in the system (apart from Triton) that can be spectroscopically characterized with the current generation of Earth-based telescopes. We report our results on spectral characterization of Nereid using its reflectance spectrum from 0.8-2.4 $μm$, providing the first measurements over the range of 0.8-1.4 $μm$. We detect spectral absorption features of crystalline water ice in close agreement with previous measurements. We show that model fits of simple intimate mixtures including water ice do not provide simultaneous matches to absorption band depths at 1.5 and 2.0 $μm$ when accounting for the spectral continuum. Possible solutions include invoking a more complex continuum, including both crystalline and amorphous water ice, and allowing for sub-micron sized grains. We show that mixtures including magnetite and the CM2 chondrite Murchison provide a flexible framework for interpreting spectral variation of bodies with neutral-sloped spectra like that of Nereid. Magnetite in particular provides a good match to the spectral continuum without requiring the presence of Tholin-like organics. We note that carbonaceous chondrites and their components may be useful analogs for the non-ice components of outer solar system bodies, consistent with recent findings by Fraser et al. (2019). Comparison to spectra of large TNOs and satellites of Uranus show that Nereid's low albedo, deep water bands, and neutral color is distinct from many other icy objects, but such comparisons are limited by incomplete understanding of spectral variability among $\sim$100km-sized icy bodies.
△ Less
Submitted 22 June, 2021;
originally announced June 2021.
-
FUV Observations of the Inner Coma of 46P/Wirtanen
Authors:
John W. Noonan,
Walter M. Harris,
Steven Bromley,
Davide Farnocchia,
Jian-Yang Li,
Kathleen E. Mandt,
Joel Wm. Parker,
Kumar Venkataramani,
Dennis Bodewits
Abstract:
Far ultraviolet observations of comets yield information about the energetic processes that dissociate the sublimated gases from their primitive surfaces. Understanding which emission processes are dominant, their effects on the observed cometary spectrum, and how to properly invert the spectrum back to composition of the presumably pristine surface ices of a comet nuclei are all critical componen…
▽ More
Far ultraviolet observations of comets yield information about the energetic processes that dissociate the sublimated gases from their primitive surfaces. Understanding which emission processes are dominant, their effects on the observed cometary spectrum, and how to properly invert the spectrum back to composition of the presumably pristine surface ices of a comet nuclei are all critical components for proper interpretation and analysis of comets. The close approach of comet 46P/Wirtanen in 2018-2019 provided a unique opportunity to study the inner most parts of a cometary coma with the Hubble Space Telescope Cosmic Origins Spectrograph, rarely accessible with remote observations, at length scales (100's of km) and wavelengths (900-1430 Angstroms) previously probed only by the European Space Agency's Rosetta spacecraft. Our observations show a complex picture for the inner coma; atomic production rates for H and O that show water is the dominant source of both, an abundance of atomic sulfur that is difficult to explain with the lifetimes of common sulfur parent molecules, and a density distribution that is poorly fit with both Haser and vectorial models.
△ Less
Submitted 8 December, 2020;
originally announced December 2020.
-
Contemporaneous Multi-Wavelength and Precovery Observations of Active Centaur P/2019 LD2 (ATLAS
Authors:
Theodore Kareta,
Laura M. Woodney,
Charles Schambeau,
Yanga Fernandez,
Olga Harrington Pinto,
Kacper Wierzchos,
M. Womack,
S. J. Bus,
Jordan Steckloff,
Gal Sarid,
Kathryn Volk,
Walter M. Harris,
Vishnu Reddy
Abstract:
Gateway Centaur and Jupiter co-orbital P/2019 LD2 (ATLAS) (Sarid et al. 2019) provides the first opportunity to observe the migration of a Solar System small body from a Centaur orbit to a Jupiter Family Comet (JFC) four decades from now (Kareta et al., 2020; Hsieh et al. 2020). The Gateway transition region is beyond where water ice can power cometary activity, and coma production there is as poo…
▽ More
Gateway Centaur and Jupiter co-orbital P/2019 LD2 (ATLAS) (Sarid et al. 2019) provides the first opportunity to observe the migration of a Solar System small body from a Centaur orbit to a Jupiter Family Comet (JFC) four decades from now (Kareta et al., 2020; Hsieh et al. 2020). The Gateway transition region is beyond where water ice can power cometary activity, and coma production there is as poorly understood as in all Centaurs. We present contemporaneous multi-wavelength observations of LD2 from 2020 July 2-4: Gemini-North visible imaging, NASA IRTF near-infrared spectroscopy, and ARO SMT millimeter-wavelength spectroscopy. Precovery DECam images limit the nucleu's effective radius to <=1.2 km and no large outbursts were seen in archival Catalina Sky Survey observations. LD2's coma has g'-r'=0.70+/-0.07, r'-i'=0.26+/-0.07, a dust production rate of ~10-20 kg/s, and an outflow velocity between v~0.6-3.3 m/s. We did not detect CO towards LD2 on 2020 July 2-3, with a 3-sigma upper limit of Q(CO) < 4.4 * 10^27 mol/s (<200 kg/s). Near-infrared spectra show evidence for water ice at the 1-10% level depending on grain size. Spatial profiles and archival data are consistent with sustained activity. The evidence supports the hypothesis that LD2 is a typical small Centaur that will become a typical JFC, and thus it is critical to understanding the transition between these two populations. Finally, we discuss potential strategies for a community-wide, long baseline monitoring effort.
△ Less
Submitted 11 February, 2021; v1 submitted 19 November, 2020;
originally announced November 2020.
-
Implementing Primordial Binaries in Simulations of Star Cluster Formation with a Hybrid MHD and Direct N-Body Method
Authors:
Claude Cournoyer-Cloutier,
Aaron Tran,
Sean Lewis,
Joshua E. Wall,
William E. Harris,
Mordecai-Mark Mac Low,
Stephen L. W. McMillan,
Simon Portegies Zwart,
Alison Sills
Abstract:
The fraction of stars in binary systems within star clusters is important for their evolution, but what proportion of binaries form by dynamical processes after initial stellar accretion remains unknown. In previous work, we showed that dynamical interactions alone produced too few low-mass binaries compared to observations. We therefore implement an initial population of binaries in the coupled M…
▽ More
The fraction of stars in binary systems within star clusters is important for their evolution, but what proportion of binaries form by dynamical processes after initial stellar accretion remains unknown. In previous work, we showed that dynamical interactions alone produced too few low-mass binaries compared to observations. We therefore implement an initial population of binaries in the coupled MHD and direct N-body star cluster formation code Torch. We compare simulations with, and without, initial binary populations and follow the dynamical evolution of the binary population in both sets of simulations, finding that both dynamical formation and destruction of binaries take place. Even in the first few million years of star formation, we find that an initial population of binaries is needed at all masses to reproduce observed binary fractions for binaries with mass ratios above the $q \geq 0.1$ detection limit. Our simulations also indicate that dynamical interactions in the presence of gas during cluster formation modify the initial distributions towards binaries with smaller primary masses, larger mass ratios, smaller semi-major axes and larger eccentricities. Systems formed dynamically do not have the same properties as the initial systems, and systems formed dynamically in the presence of an initial population of binaries differ from those formed in simulations with single stars only. Dynamical interactions during the earliest stages of star cluster formation are important for determining the properties of binary star systems.
△ Less
Submitted 11 November, 2020;
originally announced November 2020.
-
Measuring Dark Matter in Galaxies: The Mass Fraction Within 5 Effective Radii
Authors:
William E. Harris,
Rhea-Silvia Remus,
Gretchen L. H. Harris,
Iurii V. Babyk
Abstract:
Large galaxies may contain an "atmosphere" of hot interstellar X-ray gas, and the temperature and radial density profile of this gas can be used to measure the total mass of the galaxy contained within a given radius r. We use this technique for 102 early-type galaxies (ETGs) with stellar masses M_* > 10^10 M_Sun, to evaluate the mass fraction of dark matter (DM) within the fiducial radius r = 5 r…
▽ More
Large galaxies may contain an "atmosphere" of hot interstellar X-ray gas, and the temperature and radial density profile of this gas can be used to measure the total mass of the galaxy contained within a given radius r. We use this technique for 102 early-type galaxies (ETGs) with stellar masses M_* > 10^10 M_Sun, to evaluate the mass fraction of dark matter (DM) within the fiducial radius r = 5 r_e, denoted f_5 = f_{DM}(5r_e). On average, these systems have a median f_5 = 0.8 - 0.9 with a typical galaxy-to-galaxy scatter +-0.15. Comparisons with mass estimates made through the alternative techniques of satellite dynamics (e.g. velocity distributions of globular clusters, planetary nebulae, satellite dwarfs) as well as strong lensing show encouraging consistency over the same range of stellar mass. We find that many of the disk galaxies (S0/SA0/SB0) have a significantly higher mean $f_5$ than do the pure ellipticals, by Delta f_5 = 0.1. We suggest that this higher level may be a consequence of sparse stellar haloes and quieter histories with fewer major episodes of feedback or mergers. Comparisons are made with the Magneticum Pathfinder suite of simulations for both normal and centrally dominant "Brightest Cluster" galaxies. Though the observed data exhibit somewhat larger scatter at a given galaxy mass than do the simulations, the mean level of DM mass fraction for all classes of galaxies is in good first-order agreement with the simulations. Lastly, we find that the group galaxies with stellar masses near M_* ~ 10^11 M_Sun have relatively more outliers at low $f_5$ than in other mass ranges, possibly the result of especially effective AGN feedback in that mass range leading to expansion of their dark matter halos.
△ Less
Submitted 27 October, 2020;
originally announced October 2020.
-
Asteroid Thermal Inertia Estimates from Remote Infrared Observations: The Effects of Surface Roughness and Rotation Rate
Authors:
Alan W. Harris,
Line Drube
Abstract:
The thermal inertia of an asteroid's surface can provide insight into regolith properties, such as the presence of a layer of fine dust, the density and thermal conductivity of a rocky surface, and, together with other observational data, mineralogy. Knowledge of the surface characteristics of asteroids is important for planetary defense initiatives and the extraction of resources ("asteroid minin…
▽ More
The thermal inertia of an asteroid's surface can provide insight into regolith properties, such as the presence of a layer of fine dust, the density and thermal conductivity of a rocky surface, and, together with other observational data, mineralogy. Knowledge of the surface characteristics of asteroids is important for planetary defense initiatives and the extraction of resources ("asteroid mining"). A simple means of estimating asteroid thermal inertia has been proposed by Harris & Drube, which is suitable for application to large sets of thermal-infrared observational data, such as those obtained by infrared space telescopes. We compare results from the Harris-Drube estimator with recently published values of asteroid thermal inertia from detailed thermophysical modeling, and provide an explanation in terms of reduced surface roughness for some discrepant results. Smooth surfaces covered in fine dust may provide an explanation for the unexpectedly low values of thermal inertia derived from thermophysical modeling for some slowly rotating main-belt asteroids (MBAs). In the case of near-Earth objects (NEOs) we show that results from the estimator are in good agreement with those from thermophysical modeling, with just a few exceptions. We discuss the special cases of the NEOs (101955) Bennu, (162173) Ryugu, and (29075) 1950 DA in the context of results from our estimator. Given the data requirements and complexity of thermophysical modeling, data-analysis tools based on relatively simple concepts can play an important role in allowing "quick-look" assessment of thermal-infrared data of asteroids, especially NEOs.
△ Less
Submitted 5 October, 2020;
originally announced October 2020.
-
Strength In Diversity: Small Bodies as the Most Important Objects in Planetary Sciences
Authors:
Laura M. Woodney,
Andrew S. Rivkin,
Walter Harris,
Barbara A. Cohen,
Gal Sarid,
Maria Womack,
Olivier Barnouin,
Kat Volk,
Rachel Klima,
Yanga R. Fernandez,
Jordan K. Steckloff,
Paul A. Abell
Abstract:
Small bodies, the unaccreted leftovers of planetary formation, are often mistaken for the leftovers of planetary science in the sense that they are everything else after the planets and their satellites (or sometimes just their regular satellites) are accounted for. This mistaken view elides the great diversity of compositions, histories, and present-day conditions and processes found in the small…
▽ More
Small bodies, the unaccreted leftovers of planetary formation, are often mistaken for the leftovers of planetary science in the sense that they are everything else after the planets and their satellites (or sometimes just their regular satellites) are accounted for. This mistaken view elides the great diversity of compositions, histories, and present-day conditions and processes found in the small bodies, and the interdisciplinary nature of their study. Understanding small bodies is critical to planetary science as a field, and we urge planetary scientists and our decision makers to continue to support science-based mission selections and to recognize that while small bodies have been grouped together for convenience, the diversity of these objects in terms of composition, mass, differentiation, evolution, activity, dynamical state, physical structure, thermal environment, thermal history, and formation vastly exceeds the observed variability in the major planets and their satellites. Treating them as a monolithic group with interchangeable members does a grave injustice to the range of fundamental questions they address. We advocate for a deep and ongoing program of missions, telescopic observations, R and A funding, and student support that respects this diversity.
△ Less
Submitted 14 August, 2020;
originally announced August 2020.
-
P/2019 LD2 (ATLAS): An Active Centaur in Imminent Transition to the Jupiter Family
Authors:
Jordan Steckloff,
Gal Sarid,
Kathryn Volk,
Theodore Kareta,
Maria Womack,
Walter Harris,
Laura Woodney,
Charles Schambeau
Abstract:
The recently discovered object P/2019 LD2 (ATLAS) was initially thought to be a Jupiter Trojan asteroid, until dynamical studies and the appearance of persistent cometary activity revealed that this object is actually an active Centaur. However, the dynamical history, thermal environment, and impact of such environments on the activity of 2019 LD2 are poorly understood. Here we conduct dynamical s…
▽ More
The recently discovered object P/2019 LD2 (ATLAS) was initially thought to be a Jupiter Trojan asteroid, until dynamical studies and the appearance of persistent cometary activity revealed that this object is actually an active Centaur. However, the dynamical history, thermal environment, and impact of such environments on the activity of 2019 LD2 are poorly understood. Here we conduct dynamical simulations to constrain its orbital history and resulting thermal environment over the past 3000 years. We find that 2019 LD2 is currently in the vicinity of a dynamical "Gateway" that facilitates the majority of transitions from the Centaur population into the Jupiter Family of Comets (JFC population; Sarid et al. 2019). Our calculations show that it is unlikely to have spent significant amounts of time in the inner solar system, suggesting that its nucleus is relatively pristine in terms of physical, chemical, and thermal processing through its history. This could explain its relatively high level of distant activity as a recently activated primordial body. Finally, we find that the median frequency of transition from the Gateway population into the JFC population varies from once every ~3 years to less than once every 70 years, if 2019 LD2's nucleus is ~1 km in radius or greater than 3 km in radius. Forward modeling of 2019 LD2 shows that it will transition into the JFC population in 2063, representing the first known opportunity to observe the evolution of an active Centaur nucleus as it experiences this population-defining transition.
△ Less
Submitted 9 October, 2020; v1 submitted 6 August, 2020;
originally announced August 2020.
-
The carbon monoxide-rich interstellar comet 2I/Borisov
Authors:
D. Bodewits,
J. W. Noonan,
P. D. Feldman,
M. T. Bannister,
D. Farnocchia,
W. M. Harris,
J. -Y. Li,
K. E. Mandt,
J. Wm. Parker,
Z. Xing
Abstract:
Interstellar comets offer direct samples of volatiles from distant protoplanetary disks. 2I/Borisov is the first notably active interstellar comet discovered in our solar system[1]. Comets are condensed samples of the gas, ice, and dust that were in a star's protoplanetary disk during the formation of its planets and inform our understanding on how chemical compositions and abundances vary with di…
▽ More
Interstellar comets offer direct samples of volatiles from distant protoplanetary disks. 2I/Borisov is the first notably active interstellar comet discovered in our solar system[1]. Comets are condensed samples of the gas, ice, and dust that were in a star's protoplanetary disk during the formation of its planets and inform our understanding on how chemical compositions and abundances vary with distance from the central star. Their orbital migration moves volatiles[2], organic material, and prebiotic chemicals in their host system[3]. In our solar system, hundreds of comets have been observed remotely, and a few have been studied up close by space missions[4]. However, knowledge of extrasolar comets has been limited to what could be gleaned from distant, unresolved observations of cometary regions around other stars, with only one detection of carbon monoxide[5]. Here we report that the coma of 2I/Borisov contains significantly more CO than H2O gas, with abundances of at least 173%, more than three times higher than previously measured for any comet in the inner (<2.5 au) solar system[4]. Our ultraviolet observations of 2I/Borisov provide the first glimpse into the ice content and chemical composition of the protoplanetary disk of another star that is substantially different from our own.
△ Less
Submitted 19 April, 2020;
originally announced April 2020.
-
Planetary Science with Astrophysical Assets: Defining the Core Capabilities of Platforms
Authors:
James Bauer,
Stefanie Milam,
Gordon Bjoraker,
Sean Carey,
Doris Daou,
Leigh Fletcher,
Walt Harris,
Paul Hartogh,
Christine Hartzell,
Amanda Hendrix,
Carrie Nugent,
Andy Rivkin,
Timothy Swindle,
Cristina Thomas,
Matthew Tiscareno,
Geronimo Villanueva,
Scott Wolk
Abstract:
We seek to compile a uniform set of basic capabilities and needs to maximize the yield of Solar System science with future Astrophysics assets while allowing those assets to achieve their Astrophysics priorities. Within considerations of cost and complexity, inclusion of capabilities that make a particular platform useable to planetary science provide a critical advantage over platforms lacking su…
▽ More
We seek to compile a uniform set of basic capabilities and needs to maximize the yield of Solar System science with future Astrophysics assets while allowing those assets to achieve their Astrophysics priorities. Within considerations of cost and complexity, inclusion of capabilities that make a particular platform useable to planetary science provide a critical advantage over platforms lacking such capabilities.
△ Less
Submitted 13 February, 2020;
originally announced February 2020.
-
The PIPER Survey: I. An Initial Look at the Intergalactic Globular Cluster Population in the Perseus Cluster
Authors:
William E. Harris,
Rachel A. Brown,
Patrick R. Durrell,
Aaron J. Romanowsky,
John Blakeslee,
Jean Brodie,
Steven Janssens,
Thorsten Lisker,
Sakurako Okamoto,
Carolin Wittmann
Abstract:
We describe the goals and first results of a Program for Imaging of the PERseus cluster of galaxies (PIPER). The first phase of the program builds on imaging of fields obtained with the Hubble Space Telescope (HST) ACS/WFC and WFC3/UVIS cameras. Our PIPER target fields with HST include major early-type galaxies including the active central giant NGC 1275; known Ultra-Diffuse Galaxies; and the Intr…
▽ More
We describe the goals and first results of a Program for Imaging of the PERseus cluster of galaxies (PIPER). The first phase of the program builds on imaging of fields obtained with the Hubble Space Telescope (HST) ACS/WFC and WFC3/UVIS cameras. Our PIPER target fields with HST include major early-type galaxies including the active central giant NGC 1275; known Ultra-Diffuse Galaxies; and the Intracluster Medium. The resulting two-color photometry in F475W and F814W reaches deep enough to resolve and measure the globular cluster (GC) populations in the Perseus member galaxies. Here we present initial results for eight pairs of outer fields that confirm the presence of Intergalactic GCs (IGCs) in fields as distant as 740 kpc from the Perseus center (40\% of the virial radius of the cluster). Roughly 90% of these IGCs are identifiably blue (metal-poor) but there is a clear trace of a red (metal-rich) component as well, even at these very remote distances.
△ Less
Submitted 17 January, 2020;
originally announced January 2020.
-
The Strikingly Metal-Rich Halo of the Sombrero Galaxy
Authors:
Roger E. Cohen,
Paul Goudfrooij,
Matteo Correnti,
Oleg Y. Gnedin,
William E. Harris,
Rupali Chandar,
Thomas H. Puzia,
Ruben Sanchez-Janssen
Abstract:
The nature of the Sombrero galaxy (M 104 = NGC 4594) has remained elusive despite many observational studies at a variety of wavelengths. Here we present Hubble Space Telescope imaging of two fields at $\sim$16 and 33 kpc along the minor axis to examine stellar metallicity gradients in the extended spheroid. We use this imaging, extending more than 2 mag below the tip of the red giant branch (TRGB…
▽ More
The nature of the Sombrero galaxy (M 104 = NGC 4594) has remained elusive despite many observational studies at a variety of wavelengths. Here we present Hubble Space Telescope imaging of two fields at $\sim$16 and 33 kpc along the minor axis to examine stellar metallicity gradients in the extended spheroid. We use this imaging, extending more than 2 mag below the tip of the red giant branch (TRGB), in combination with artificial star tests to forward model observed color-magnitude diagrams (CMDs), measuring metallicity distribution functions (MDFs) at different radii along the minor axis. An important and unexpected result is that the halo of the Sombrero is strikingly metal-rich: even the outer field, located at $\sim$17 effective radii of the bulge, has a median metallicity [Z/H]$\sim$-0.15 and the fraction of stars with [Z/H]<-1.0 is negligible. This is unprecedented among massive galaxy halos studied to date, even among giant ellipticals. We find significant radial metallicity gradients, characterized by an increase in the fraction of metal-poor stars with radius and a gradient in median metallicity of $\sim$-0.01 dex/kpc. The density profile is well fit by power laws with slopes that exhibit a dependence on metallicity, with flatter slopes for more metal-poor stars. We discuss our results in the context of recent stellar MDF studies of other nearby galaxies and potential formation scenarios for the Sombrero galaxy.
△ Less
Submitted 7 January, 2020; v1 submitted 6 January, 2020;
originally announced January 2020.
-
Globular Cluster Systems and X-Ray Atmospheres in Galaxies
Authors:
Gretchen L. H. Harris,
Iurii V. Babyk,
William E. Harris,
Brian R. McNamara
Abstract:
We compare the empirical relationships between the mass of a galaxy's globular system M_GCS, the gas mass in the hot X-ray atmosphere M_X within a fiducial radius of 5 r_e, the total gravitational mass M_grav within 5 r_e, and lastly the total halo mass M_h calibrated from weak lensing. We use a sample of 45 early-type galaxies (ETGs) for which both GCS and X-ray data are available; all the galaxi…
▽ More
We compare the empirical relationships between the mass of a galaxy's globular system M_GCS, the gas mass in the hot X-ray atmosphere M_X within a fiducial radius of 5 r_e, the total gravitational mass M_grav within 5 r_e, and lastly the total halo mass M_h calibrated from weak lensing. We use a sample of 45 early-type galaxies (ETGs) for which both GCS and X-ray data are available; all the galaxies in our sample are relatively high-mass ones with M_h > 10^12 M_sun. We find that M_X ~ M_h^1.0, similar to the previously known scaling relation M_GCS ~ M_h^1.0. Both components scale much more steeply than the more well known dependence of total stellar mass M_star ~ M_h^0.3 for luminous galaxies. These results strengthen previous suggestions that feedback had little effect on formation of the globular cluster system. The current data are also used to measure the relative mass fractions of baryonic matter and dark matter (DM) within 5 r_e. We find a strikingly uniform mean of <f_DM> = 0.83 with few outliers and an rms scatter of +-0.07. This result is in good agreement with two recent suites of hydrodynamic galaxy formation models.
△ Less
Submitted 19 November, 2019;
originally announced November 2019.
-
Physical Characterization of the December 2017 Outburst of the Centaur 174P/Echeclus
Authors:
Theodore Kareta,
Benjamin Sharkey,
John Noonan,
Kathryn Volk,
Vishnu Reddy,
Walter Harris,
Richard Miles
Abstract:
The Centaurs are the small solar system bodies intermediate between the active inner solar system Jupiter Family Comets and their inactive progenitors in the trans-Neptunian region. Among the fraction of Centaurs which show comet-like activity, 174P/Echeclus is best known for its massive 2005 outburst in which a large apparently active fragment was ejected above the escape velocity from the primar…
▽ More
The Centaurs are the small solar system bodies intermediate between the active inner solar system Jupiter Family Comets and their inactive progenitors in the trans-Neptunian region. Among the fraction of Centaurs which show comet-like activity, 174P/Echeclus is best known for its massive 2005 outburst in which a large apparently active fragment was ejected above the escape velocity from the primary nucleus. We present visible imaging and near-infrared spectroscopy of Echeclus during the first week after its December 2017 outburst taken at the Faulkes North & South Telescopes and the NASA IRTF, the largest outburst since 2005. The coma was seen to be highly asymmetric. A secondary peak was seen in the near-infrared 2D spectra, which is strongly hinted at in the visible images, moving hyperbolically with respect to the nucleus. The retrieved reflectance spectrum of Echelcus is consistent with the unobscured nucleus but becomes bluer when a wider extraction aperture is used. We find that Echeclus's coma is best explained as dominated by large blue dust grains, which agrees with previous work. We also conducted a high-resolution orbital integration of Echeclus's recent evolution and found no large orbital changes that could drive its modern evolution. We interpret the second peak in the visible and near-infrared datasets as a large cloud of larger-than-dust debris ejected at the time of outburst. If Echeclus is typical of the Centaurs, there may be several debris ejection or fragmentation events per year on other Centaurs that are going unnoticed.
△ Less
Submitted 21 October, 2019;
originally announced October 2019.
-
Carbon Chain Depletion of 2I/Borisov
Authors:
Theodore Kareta,
Jennifer Andrews,
John W. Noonan,
Walter M. Harris,
Nathan Smith,
Patrick O'Brien,
Benjamin N. L. Sharkey,
Vishnu Reddy,
Alessondra Springmann,
Cassandra Lejoly,
Kathryn Volk,
Albert Conrad,
Christian Veillet
Abstract:
The composition of comets in the Solar System come in multiple groups thought to encode information about their formation in different regions of the outer protosolar disk. The recent discovery of the second interstellar object, 2I/Borisov, allows for spectroscopic investigations into its gas content and a preliminary classification of it within the Solar System comet taxonomies to test the applic…
▽ More
The composition of comets in the Solar System come in multiple groups thought to encode information about their formation in different regions of the outer protosolar disk. The recent discovery of the second interstellar object, 2I/Borisov, allows for spectroscopic investigations into its gas content and a preliminary classification of it within the Solar System comet taxonomies to test the applicability of planetesimal formation models to other stellar systems. We present spectroscopic and imaging observations from 2019 September 20th to October 26th at the Bok, MMT, and LBT telescopes. We identify CN in the comet's spectrum and set precise upper limits on the abundance of C2 on all dates. We use a Haser model to convert our integrated fluxes to production rates and find Q(CN) = 5.0 +/- 2.0 * 10^24 mol/s on September 20th and Q(CN) = 1.1 - 1.9 * 10^24 mol/s on later dates, both consistent with contemporaneous observations. We set our lowest upper limit on a C2 production rate, Q(C2) < 1.6 * 10^23 mol/s, on October 10th. The measured ratio upper limit for that date, Q(C2)/Q(CN) < 0.095 indicates that 2I/Borisov is strongly in the (carbon chain) 'depleted' taxonomic group. The only comparable Solar System comets have detected ratios near this limit, making 2I/Borisov statistically likely to be more depleted than any known comet. Most 'depleted' comets are Jupiter Family Comets, perhaps indicating a similiarity in formation conditions between the most depleted of the JFCs and 2I/Borisov. More work is needed to understand the applicability of our knowledge of Solar System comet taxonomies onto interstellar objects, and we discuss future work that could help clarify the usefulness of the approach.
△ Less
Submitted 25 November, 2019; v1 submitted 8 October, 2019;
originally announced October 2019.
-
29P/Schwassmann-Wachmann 1, A Centaur in the Gateway to the Jupiter-Family Comets
Authors:
Gal Sarid,
Kathryn Volk,
Jordan Steckloff,
Walter Harris,
Maria Womack,
Laura Woodney
Abstract:
Jupiter-family comets (JFCs) are the evolutionary products of trans-Neptunian objects (TNOs) that evolve through the giant planet region as Centaurs and into the inner solar system. Through numerical orbital evolution calculations following a large number of TNO test particles that enter the Centaur population, we have identified a short-lived dynamical Gateway, a temporary low-eccentricity region…
▽ More
Jupiter-family comets (JFCs) are the evolutionary products of trans-Neptunian objects (TNOs) that evolve through the giant planet region as Centaurs and into the inner solar system. Through numerical orbital evolution calculations following a large number of TNO test particles that enter the Centaur population, we have identified a short-lived dynamical Gateway, a temporary low-eccentricity region exterior to Jupiter through which the majority of JFCs pass. We apply an observationally based size distribution function to the known Centaur population and obtain an estimated Gateway region population. We then apply an empirical fading law to the rate of incoming JFCs implied by the the Gateway region residence times. Our derived estimates are consistent with observed population numbers for the JFC and Gateway populations. Currently, the most notable occupant of the Gateway region is 29P/Schwassmann-Wachmann 1 (SW1), a highly active, regularly outbursting Centaur. SW1's present-day, very-low-eccentricity orbit was established after a 1975 Jupiter conjunction and will persist until a 2038 Jupiter conjunction doubles its eccentricity and pushes its semi-major axis out to its current aphelion. Subsequent evolution will likely drive SW1's orbit out of the Gateway region, perhaps becoming one of the largest JFCs in recorded history. The JFC Gateway region coincides with a heliocentric distance range where the activity of observed cometary bodies increases significantly. SW1's activity may be typical of the early evolutionary processing experienced by most JFCs. Thus, the Gateway region, and its most notable occupant SW1, are critical to both the dynamical and physical transition between Centaurs and JFCs.
△ Less
Submitted 12 September, 2019; v1 submitted 12 August, 2019;
originally announced August 2019.
-
Spitzer Albedos of Near-Earth Objects
Authors:
Annika Gustafsson,
David E. Trilling,
Michael Mommert,
Andrew McNeill,
Joseph L. Hora,
Howard A. Smith,
Stephan Hellmich,
Stefano Mottola,
Alan W. Harris
Abstract:
Thermal infrared observations are the most effective way to measure asteroid diameter and albedo for a large number of near-Earth objects. Major surveys like NEOWISE, NEOSurvey, ExploreNEOs, and NEOLegacy find a small fraction of high albedo objects that do not have clear analogs in the current meteorite population. About 8% of Spitzer-observed near-Earth objects have nominal albedo solutions grea…
▽ More
Thermal infrared observations are the most effective way to measure asteroid diameter and albedo for a large number of near-Earth objects. Major surveys like NEOWISE, NEOSurvey, ExploreNEOs, and NEOLegacy find a small fraction of high albedo objects that do not have clear analogs in the current meteorite population. About 8% of Spitzer-observed near-Earth objects have nominal albedo solutions greater than 0.5. This may be a result of lightcurve variability leading to an incorrect estimate of diameter or inaccurate absolute visual magnitudes. For a sample of 23 high albedo NEOs we do not find that their shapes are significantly different from the McNeill et al. (2019) near-Earth object shape distribution. We performed a Monte Carlo analysis on 1505 near-Earth objects observed by Spitzer, sampling the visible and thermal fluxes of all targets to determine the likelihood of obtaining a high albedo erroneously. Implementing the McNeill shape distribution for near-Earth objects, we provide an upper-limit on the geometric albedo of 0.5+/-0.1 for the near-Earth population.
△ Less
Submitted 17 June, 2019;
originally announced June 2019.
-
New Insights into the Physical Conditions and Internal Structure of a Candidate Proto-Globular Cluster
Authors:
Molly K. Finn,
Kelsey E. Johnson,
Crystal L. Brogan,
Christine D. Wilson,
Remy Indebetouw,
William E. Harris,
Julia Kamenetzky,
Ashley Bemis
Abstract:
We present $\sim$0.1" resolution ($\sim$10 pc) ALMA observations of a molecular cloud identified in the merging Antennae galaxies with the potential to form a globular cluster, nicknamed the ``Firecracker.' Since star formation has not yet begun at an appreciable level, this cloud provides an example of what the birth environment of a globular cluster may have looked like before stars form and dis…
▽ More
We present $\sim$0.1" resolution ($\sim$10 pc) ALMA observations of a molecular cloud identified in the merging Antennae galaxies with the potential to form a globular cluster, nicknamed the ``Firecracker.' Since star formation has not yet begun at an appreciable level, this cloud provides an example of what the birth environment of a globular cluster may have looked like before stars form and disrupt the natal conditions. Using emission from $^{12}$CO(2-1), $^{12}$CO(3-2), $^{13}$CO(2-1), HCN(4-3), and HCO$^+$(4-3), we are able to resolve the cloud's structure and find that it has a characteristic radius of 22 pc and a mass of 1--9$\times10^6 M_\odot$. We also constrain the abundance ratios of $^{12}$CO/$^{13}$CO and H$_2$/\twelveCO. Based on the calculated mass, we determine that the commonly used CO-to-H$_2$ conversion factor varies spatially, with average values in the range $X_{CO}=(0.12-1.1)\times10^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$. We demonstrate that if the cloud is bound (as is circumstantially suggested by its bright, compact morphology), an external pressure of $P/k > 10^8$ K cm$^{-3}$ is required. This would be consistent with theoretical expectations that globular cluster formation requires high pressure environments. The position-velocity diagram of the cloud and its surrounding material suggests that this pressure may be produced by the collision of filaments. The radial profile of the column density can be fit with both a Gaussian and Bonnor-Ebert profile. The relative line strengths of HCN and HCO$^+$ in this region suggest that these molecular lines can be used as tracers for the evolutionary stage of a cluster.
△ Less
Submitted 20 March, 2019;
originally announced March 2019.
-
The large-scale structure of the halo of the Andromeda galaxy II. Hierarchical structure in the Pan-Andromeda Archaeological Survey
Authors:
Alan W. McConnachie,
Rodrigo Ibata,
Nicolas Martin,
Annette M. N. Ferguson,
Michelle Collins,
Stephen Gwyn,
Mike Irwin,
Geraint F. Lewis,
A. Dougal Mackey,
Tim Davidge,
Veronica Arias,
Anthony Conn,
Patrick Cote,
Denija Crnojevic,
Avon Huxor,
Jorge Penarrubia,
Chelsea Spengler,
Nial Tanvir,
David Valls-Gabaud,
Arif Babul,
Pauline Barmby,
Nicholas F. Bate,
Edouard Bernard,
Scott Chapman,
Aaron Dotter
, et al. (7 additional authors not shown)
Abstract:
The Pan-Andromeda Archaeological Survey is a survey of $>400$ square degrees centered on the Andromeda (M31) and Triangulum (M33) galaxies that has provided the most extensive panorama of a $L_\star$ galaxy group to large projected galactocentric radii. Here, we collate and summarise the current status of our knowledge of the substructures in the stellar halo of M31, and discuss connections betwee…
▽ More
The Pan-Andromeda Archaeological Survey is a survey of $>400$ square degrees centered on the Andromeda (M31) and Triangulum (M33) galaxies that has provided the most extensive panorama of a $L_\star$ galaxy group to large projected galactocentric radii. Here, we collate and summarise the current status of our knowledge of the substructures in the stellar halo of M31, and discuss connections between these features. We estimate that the 13 most distinctive substructures were produced by at least 5 different accretion events, all in the last 3 or 4 Gyrs. We suggest that a few of the substructures furthest from M31 may be shells from a single accretion event. We calculate the luminosities of some prominent substructures for which previous estimates were not available, and we estimate the stellar mass budget of the outer halo of M31. We revisit the problem of quantifying the properties of a highly structured dataset; specifically, we use the OPTICS clustering algorithm to quantify the hierarchical structure of M31's stellar halo, and identify three new faint structures. M31's halo, in projection, appears to be dominated by two `mega-structures', that can be considered as the two most significant branches of a merger tree produced by breaking M31's stellar halo into smaller and smaller structures based on the stellar spatial clustering. We conclude that OPTICS is a powerful algorithm that could be used in any astronomical application involving the hierarchical clustering of points. The publication of this article coincides with the public release of all PAndAS data products.
△ Less
Submitted 25 October, 2018; v1 submitted 18 October, 2018;
originally announced October 2018.
-
Estimating the Milky Way's Mass via Hierarchical Bayes: A Blind Test on MUGS2 Simulated Galaxies
Authors:
Gwendolyn M. Eadie,
Benjamin W. Keller,
William E. Harris
Abstract:
In a series of three papers, Eadie et al. developed a hierarchical Bayesian method to estimate the Milky Way Galaxy's mass given a physical model for the potential, a measurement model, and kinematic data of test particles such as globular clusters (GCs) or halo stars in the Galaxy's halo. The Galaxy's virial mass was found to have a 95\% Bayesian credible region (c.r.) of…
▽ More
In a series of three papers, Eadie et al. developed a hierarchical Bayesian method to estimate the Milky Way Galaxy's mass given a physical model for the potential, a measurement model, and kinematic data of test particles such as globular clusters (GCs) or halo stars in the Galaxy's halo. The Galaxy's virial mass was found to have a 95\% Bayesian credible region (c.r.) of $(0.67, 1.09) \times 10^{12} M_{\odot}$. In the present study, we test the hierarchical Bayesian method against simulated galaxies created in the McMaster Unbiased Galaxy Simulations 2 (MUGS2), for which the true mass is known. We estimate the masses of MUGS2 galaxies using GC analogs from the simulations as tracers. The analysis, completed as a blind test, recovers the true $M_{200}$ of the MUGS2 galaxies within 95\% Bayesian c.r. in 8 out of 18 cases. Of the 10 galaxy masses that were not recovered within the 95\% c.r., a large subset have posterior distributions that occupy extreme ends of the parameter space allowed by the priors. A few incorrect mass estimates are explained by the exceptional evolution history of the galaxies. We also find evidence that the model cannot describe both the galaxies' inner and outer structure simultaneously in some cases. After removing the GC analogs associated with the galactic disks, the true masses were found more reliably (13 out of 18 were predicted within the c.r.). Finally, we discuss how representative the GC analogs are of the real GC population in the Milky Way.
△ Less
Submitted 4 October, 2018;
originally announced October 2018.
-
On the Origin of Multiple Populations During Massive Star Cluster Formation
Authors:
Corey S. Howard,
Ralph E. Pudritz,
Alison Sills,
William E. Harris
Abstract:
We investigate the possibility that multiple populations in globular clusters arise as a natural by-product of massive star-cluster formation. We use 3D radiative hydrodynamics simulations for the formation of young massive clusters to track their chemical self-enrichment during their first 5 Myr. These clusters form embedded within filamentary Giant Molecular Clouds by a combination of gas accret…
▽ More
We investigate the possibility that multiple populations in globular clusters arise as a natural by-product of massive star-cluster formation. We use 3D radiative hydrodynamics simulations for the formation of young massive clusters to track their chemical self-enrichment during their first 5 Myr. These clusters form embedded within filamentary Giant Molecular Clouds by a combination of gas accretion and rapid merging of protoclusters. Chemical enrichment is a dynamic process happening as the young cluster assembles, so that the original (1P) and enriched (2P) subpopulations of stars form almost simultaneously. Here we test two simple and opposite extremes for the injection of enriched material into the intracluster gas: we assume either continuous injection in a way that tracks the star formation rate; or sudden injection by a single instantaneous event. Using helium abundance as a proxy for the enrichment, we find that realistic multiple population features can be reproduced by injecting a total helium mass amounting to a few percent of the cluster's total mass. The differences in individual growth histories can lead to widely differing 1P/2P outcomes. These models suggest that dual or multiple populations can emerge rapidly in massive star clusters undergoing the typical mode of star cluster formation.
△ Less
Submitted 5 April, 2019; v1 submitted 21 August, 2018;
originally announced August 2018.
-
A universal route for the formation of massive star clusters in giant molecular clouds
Authors:
Corey S. Howard,
Ralph E. Pudritz,
William E. Harris
Abstract:
Young massive star clusters (YMCs, with M $\geq$10$^4$ M$_{\odot}$) are proposed modern-day analogues of the globular clusters (GCs) that were products of extreme star formation in the early universe. The exact conditions and mechanisms under which YMCs form remain unknown -- a fact further complicated by the extreme radiation fields produced by their numerous massive young stars. Here we show tha…
▽ More
Young massive star clusters (YMCs, with M $\geq$10$^4$ M$_{\odot}$) are proposed modern-day analogues of the globular clusters (GCs) that were products of extreme star formation in the early universe. The exact conditions and mechanisms under which YMCs form remain unknown -- a fact further complicated by the extreme radiation fields produced by their numerous massive young stars. Here we show that GC-sized clusters are naturally produced in radiation-hydrodynamic simulations of isolated 10$^7$ M$_{\odot}$ Giant Molecular Clouds (GMCs) with properties typical of the local universe, even under the influence of radiative feedback. In all cases, these massive clusters grow to GC-level masses within 5 Myr via a roughly equal combination of filamentary gas accretion and mergers with several less massive clusters. Lowering the heavy-element abundance of the GMC by a factor of 10 reduces the opacity of the gas to radiation and better represents the high-redshift formation conditions of GCs. This results in higher gas accretion leading to a mass increase of the largest cluster by a factor of ~4. When combined with simulations of less massive GMCs (10$^{4-6}$ M$_{\odot}$), a clear relation emerges between the maximum YMC mass and the mass of the host GMC. Our results demonstrate that YMCs, and potentially GCs, are a simple extension of local cluster formation to more massive clouds and do not require suggested exotic formation scenarios.
△ Less
Submitted 21 August, 2018;
originally announced August 2018.
-
Evolution of H$_2$O Production in Comet C/2012 S1 (ISON) as Inferred from Forbidden Oxygen and OH Emission
Authors:
Adam J. McKay,
Anita L. Cochran,
Michael A. DiSanti,
Neil Dello Russo,
Harold Weaver,
Ronald J. Vervack Jr.,
Walter M. Harris,
Hideyo Kawakita
Abstract:
We present H$_2$O production rates for comet C/2012 S1 (ISON) derived from observations of [OI] and OH emission during its inbound leg, covering a heliocentric distance range of 1.8-0.44 AU. Our production rates are in agreement with previous measurements using a variety of instruments and techniques and with data from the various observatories greatly differing in their projected fields of view.…
▽ More
We present H$_2$O production rates for comet C/2012 S1 (ISON) derived from observations of [OI] and OH emission during its inbound leg, covering a heliocentric distance range of 1.8-0.44 AU. Our production rates are in agreement with previous measurements using a variety of instruments and techniques and with data from the various observatories greatly differing in their projected fields of view. The consistent results across all data suggest the absence of an extended source of H$_2$O production, for example sublimation of icy grains in the coma, or a source with spatial extent confined to the dimensions of the smallest projected field of view (in this case $<$ 1,000 km). We find that ISON had an active area of around 10 km$^2$ for heliocentric distances R$_h$ > 1.2 AU, which then decreased to about half this value from R$_h$=1.2-0.9 AU. This was followed by a rapid increase in active area at about R$_h$=0.6 AU, corresponding to the first of three major outbursts ISON experienced inside of 1 AU. The combination of a detected outburst in the light curve and rapid increase in active area likely indicates a major nucleus fragmentation event. The 5-10 km$^2$ active area observed outside of R$_h$=0.6 AU is consistent with a 50-100% active fraction for the nucleus, larger than typically observed for cometary nuclei. Although the absolute value of the active area is somewhat dependent on the thermal model employed, the changes in active area observed are consistent among models. The conclusion of a 50-100+% active fraction is robust for realistic thermal models of the nucleus. However the possibility of a contribution of a spatially unresolved distribution of icy grains cannot be discounted.
△ Less
Submitted 16 February, 2018;
originally announced February 2018.
-
A Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE).I. Introduction to the Survey
Authors:
A. Boselli,
M. Fossati,
L. Ferrarese,
S. Boissier,
G. Consolandi,
A. Longobardi,
P. Amram,
M. Balogh,
P. Barmby,
M. Boquien,
F. Boulanger,
J. Braine,
V. Buat,
D. Burgarella,
F. Combes,
T. Contini,
L. Cortese,
P. Cote,
S. Cote,
J. C. Cuillandre,
L. Drissen,
B. Epinat,
M. Fumagalli,
S. Gallagher,
G. Gavazzi
, et al. (34 additional authors not shown)
Abstract:
The Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) is a blind narrow-band Halpha+[NII] imaging survey carried out with MegaCam at the Canada-France-Hawaii Telescope. The survey covers the whole Virgo cluster region from its core to one virial radius (104 deg^2). The sensitivity of the survey is of f(Halpha) ~ 4 x 10^-17 erg sec-1 cm^-2 (5 sigma detection limit) for point sources…
▽ More
The Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) is a blind narrow-band Halpha+[NII] imaging survey carried out with MegaCam at the Canada-France-Hawaii Telescope. The survey covers the whole Virgo cluster region from its core to one virial radius (104 deg^2). The sensitivity of the survey is of f(Halpha) ~ 4 x 10^-17 erg sec-1 cm^-2 (5 sigma detection limit) for point sources and Sigma (Halpha) ~ 2 x 10^-18 erg sec^-1 cm^-2 arcsec^-2 (1 sigma detection limit at 3 arcsec resolution) for extended sources, making VESTIGE the deepest and largest blind narrow-band survey of a nearby cluster. This paper presents the survey in all its technical aspects, including the survey design, the observing strategy, the achieved sensitivity in both the narrow-band Halpha+[NII] and in the broad-band r filter used for the stellar continuum subtraction, the data reduction, calibration, and products, as well as its status after the first observing semester. We briefly describe the Halpha properties of galaxies located in a 4x1 deg^2 strip in the core of the cluster north of M87, where several extended tails of ionised gas are detected. This paper also lists the main scientific motivations of VESTIGE, which include the study of the effects of the environment on galaxy evolution, the fate of the stripped gas in cluster objects, the star formation process in nearby galaxies of different type and stellar mass, the determination of the Halpha luminosity function and of the Halpha scaling relations down to ~ 10^6 Mo stellar mass objects, and the reconstruction of the dynamical structure of the Virgo cluster. This unique set of data will also be used to study the HII luminosity function in hundreds of galaxies, the diffuse Halpha+[NII] emission of the Milky Way at high Galactic latitude, and the properties of emission line galaxies at high redshift.
△ Less
Submitted 8 February, 2018;
originally announced February 2018.
-
The LUVOIR Ultraviolet Multi-Object Spectrograph (LUMOS): Instrument Definition and Design
Authors:
Kevin France,
Brian Fleming,
Garrett West,
Stephan R. McCandliss,
Matthew R. Bolcar,
Walter Harris,
Leonidas Moustakas,
John M. O'Meara,
Ilaria Pascucci,
Jane Rigby,
David Schiminovich,
Jason Tumlinson,
Jean-Claude Bouret,
Christopher J. Evans,
Miriam Garcia
Abstract:
The Large Ultraviolet / Optical / Infrared Surveyor (LUVOIR) is one of four large mission concepts currently undergoing community study for consideration by the 2020 Astronomy and Astrophysics Decadal Survey. The LUVOIR Ultraviolet Multi-Object Spectrograph, LUMOS, is being designed to support all of the UV science requirements of LUVOIR, from exoplanet host star characterization to tomography of…
▽ More
The Large Ultraviolet / Optical / Infrared Surveyor (LUVOIR) is one of four large mission concepts currently undergoing community study for consideration by the 2020 Astronomy and Astrophysics Decadal Survey. The LUVOIR Ultraviolet Multi-Object Spectrograph, LUMOS, is being designed to support all of the UV science requirements of LUVOIR, from exoplanet host star characterization to tomography of circumgalactic halos to water plumes on outer solar system satellites. LUMOS offers point source and multi-object spectroscopy across the UV bandpass, with multiple resolution modes to support different science goals. The instrument will provide low (R = 8,000-18,000) and medium (R = 30,000-65,000) resolution modes across the far-ultraviolet (FUV: 100-200 nm) and near-ultraviolet (NUV: 200-400 nm) windows, and a very low resolution mode (R = 500) for spectroscopic investigations of extremely faint objects in the FUV. Imaging spectroscopy will be accomplished over a 3 x 1.6 arcminute field-of-view by employing holographically-ruled diffraction gratings to control optical aberrations, microshutter arrays (MSA), advanced optical coatings for high-throughput in the FUV, and next generation large-format photon-counting detectors. The spectroscopic capabilities of LUMOS are augmented by an FUV imaging channel (100-200nm, 13 milliarcsecond angular resolution, 2 x 2 arcminute field-of-view) that will employ a complement of narrow and medium-band filters. We present an overview of LUMOS' observing modes and estimated performance curves for effective area, spectral resolution, and imaging performance. Example "LUMOS 100-hour Highlights" observing programs are presented to demonstrate the potential power of LUVOIR's ultraviolet spectroscopic capabilities.
△ Less
Submitted 18 September, 2017;
originally announced September 2017.