We describe the public release of the Cluster Monte Carlo (CMC) code, a parallel, star-by-star N-body code for modeling dense star clusters. CMC treats collisional stellar dynamics using Henon's method, where the cumulative effect of many two-body encounters is statistically reproduced as a single effective encounter between nearest-neighbor particles on a relaxation timescale. The star-by-star approach allows for the inclusion of additional physics, including strong gravitational three- and four-body encounters, two-body tidal and gravitational-wave captures, mass loss in arbitrary galactic tidal fields, and stellar evolution for both single and binary stars. The public release of CMC is pinned directly to the COSMIC population synthesis code, allowing dynamical star cluster simulations and population synthesis studies to be performed using identical assumptions about the stellar physics and initial conditions. As a demonstration, we present two examples of star cluster modeling: first, we perform the largest (N = 10(8)) star-by-star N-body simulation of a Plummer sphere evolving to core collapse, reproducing the expected self-similar density profile over more than 15 orders of magnitude; second, we generate realistic models for typical globular clusters, and we show that their dynamical evolution can produce significant numbers of black hole mergers with masses greater than those produced from isolated binary evolution (such as GW190521, a recently reported merger with component masses in the pulsational pair-instability mass gap).

Ultralight bosons are a proposed solution to outstanding problems in cosmology and particle physics: they provide a dark-matter candidate while potentially explaining the strong charge-parity problem. If they exist, ultralight bosons can interact with black holes through the superradiant instability. In this work we explore the consequences of this instability on the evolution of hierarchical black holes within dense stellar clusters. By reducing the spin of individual black holes, superradiance reduces the recoil velocity of merging binary black holes, which, in turn, increases the retention fraction of hierarchical merger remnants. We show that the existence of ultralight bosons with mass 2 x 10(-14) less than or similar to mu/eV less than or similar to 2 x 10(-13) would lead to an increased rate of hierarchical black hole mergers in nuclear star clusters. An ultralight boson in this energy range would result in up to approximate to 60% more present-day nuclear star clusters supporting hierarchical growth. The presence of an ultralight boson can also double the rate of intermediate-mass black hole mergers to approximate to 0.08 Gpc(-3) yr(-1) in the local universe. These results imply that a select range of ultralight boson masses can have far-reaching consequences for the population of black holes in dense stellar environments. Future studies into black hole cluster populations and the spin distribution of hierarchically formed black holes will test this scenario.

The globular cluster 47 Tucanae (47 Tuc) is one of the most massive star clusters in the Milky Way and is exceptionally rich in exotic stellar populations. For several decades it has been a favorite target of observers, and yet it is computationally very challenging to model because of its large number of stars (N greater than or similar to 10(6)) and high density. Here we present detailed and self-consistent 47 Tuc models computed with the Cluster Monte Carlo code (CMC). The models include all relevant dynamical interactions coupled to stellar and binary evolution, and reproduce various observations, including the surface brightness and velocity dispersion profiles, pulsar accelerations, and numbers of compact objects. We show that the present properties of 47 Tuc are best reproduced by adopting an initial stellar mass function that is both bottom-heavy and top-light relative to standard assumptions (as in, e.g., Kroupa 2001), and an initial Elson profile (Elson et al. 1987) that is overfilling the cluster's tidal radius. We include new prescriptions in CMC for the formation of binaries through giant star collisions and tidal captures, and we show that these mechanisms play a crucial role in the formation of neutron star binaries and millisecond pulsars in 47 Tuc; our best-fit model contains similar to 50 millisecond pulsars, 70% of which are formed through giant collisions and tidal captures. Our models also suggest that 47 Tuc presently contains up to similar to 200 stellar-mass black holes, similar to 5 binary black holes, similar to 15 low-mass X-ray binaries, and similar to 300 cataclysmic variables.

Many recent observational and theoretical studies suggest that globular clusters (GCs) host compact object populations large enough to play dominant roles in their overall dynamical evolution. Yet direct detection, particularly of black holes and neutron stars, remains rare and limited to special cases, such as when these objects reside in close binaries with bright companions. Here we examine the potential of microlensing detections to further constrain these dark populations. Based on state-of-the-art GC models from the CMC Cluster Catalog, we estimate the microlensing event rates for black holes, neutron stars, white dwarfs (WDs), and, for comparison, also for M dwarfs in Milky Way GCs, as well as the effects of different initial conditions on these rates. Among compact objects, we find that WDs dominate the microlensing rates, simply because they largely dominate by numbers. We show that microlensing detections are in general more likely in GCs with higher initial densities, especially in clusters that undergo core collapse. We also estimate microlensing rates in the specific cases of M22 and 47 Tuc using our best-fitting models for these GCs. Because their positions on the sky lie near the rich stellar backgrounds of the Galactic bulge and the Small Magellanic Cloud, respectively, these clusters are among the Galactic GCs best suited for dedicated microlensing surveys. The upcoming 10 yr survey with the Rubin Observatory may be ideal for detecting lensing events in GCs.

The total integrated emission from galaxies, known as the Extragalactic Background Light (EBL), is an important observable for understanding the history of star formation over the history of the universe. Spatial fluctuations in the infrared EBL as measured by the Cosmic Infrared Background ExpeRiment (CIBER), Spitzer and AKARI exceed the predicted signal from galaxy clustering alone. The CIBER-2 project seeks to extend CIBER observations of the EBL throughout the near infrared into the optical, through measurements above Earth's atmosphere during a suborbital sounding rocket flight. The experiment has a LN2-cooled 28.5 cm Cassegrain telescope along with three optical paths and dichroic beamsplitters, which are used to obtain three wide-field images in six broad spectral bands between 0.5-2.0 mu m. The three focal planes also contain linear variable filters (LVFs) which simultaneously take spectra with resolution R=20 across the same range. CIBER-2 is scheduled to fly multiple times on a Black Brant IX sounding rocket from White Sands Missile Range in the New Mexico desert. For the first flight, scheduled for early 2021, we have completed a variety of pre-flight optical tests, which we use to make focus adjustments, spectral response measurements, and absolute photometric calibrations. In this paper, we describe the methods behind these tests and present their results for pre-flight performance evaluation. In particular, we present measurements of the PSF for each broad spectral band, along with absolute calibration factors for each band and the LVF. Through monochromator scans, we also measure the spectral responsivity of each LVF as a function of position.

Scattered sunlight from the interplanetary dust (IPD) cloud in our solar system presents a serious foreground challenge for spectrophotometric measurements of the extragalactic background light (EBL). In this work, we report on inferred measurements of the absolute intensity of the zodiacal light (ZL) using the novel technique of Fraunhofer line spectroscopy on the deepest 8542 angstrom line of the near-infrared Ca ii absorption triplet. The measurements are performed with the narrow band spectrometer (NBS) on board the Cosmic Infrared Background Experiment sounding rocket instrument. We use the NBS data to test the accuracy of two ZL models widely cited in the literature, the Kelsall and Wright models, which have been used in foreground removal analyses that produce high and low EBL results respectively. We find a mean reduced chi (2) = 3.5 for the Kelsall model and chi (2) = 2.0 for the Wright model. The best description of our data is provided by a simple modification to the Kelsall model, which includes a free ZL offset parameter. This adjusted model describes the data with a reduced chi (2) = 1.5 and yields an inferred offset amplitude of 46 +/- 19 nW m(-2) sr(-1) extrapolated to 12500 angstrom. These measurements elude to the potential existence of a dust cloud component in the inner solar system whose intensity does not strongly modulate with the Earth's motion around the Sun.

We report the first measurement of the zodiacal light (ZL) polarization spectrum in the near-infrared between 0.8 and 1.8 mu m. Using the low-resolution spectrometer on board the Cosmic Infrared Background Experiment, calibrated for absolute spectrophotometry and spectropolarimetry, we acquire long-slit polarization spectral images of the total diffuse sky brightness toward five fields. To extract the ZL spectrum, we subtract the contribution of other diffuse radiation, such as the diffuse galactic light, the integrated starlight, and the extragalactic background light. The measured ZL polarization spectrum shows little wavelength dependence in the near-infrared, and the degree of polarization clearly varies as a function of the ecliptic coordinates and solar elongation. Among the observed fields, the North Ecliptic Pole shows the maximum degree of polarization of similar to 20%, which is consistent with an earlier observation from the Diffuse Infrared Background Experiment on board on the Cosmic Background Explorer. The measured degree of polarization and its solar elongation dependence are reproduced by an empirical scattering model in the visible band and also by a Mie scattering model for large absorptive particles, while a Rayleigh scattering model is ruled out. All of our results suggest that the interplanetary dust is dominated by large particles.

We report on the discovery of an infrared cometary nebula around PSR J1549-4848 in our Spitzer survey of a few middle-aged radio pulsars. Following the discovery, multi-wavelength imaging and spectroscopic observations of the nebula were carried out. We detected the nebula in Spitzer Infrared Array Camera 8.0, Multiband Imaging Photometer for Spitzer 24 and 70 mu m imaging, and in Spitzer IRS 7.5-14.4 mu m spectroscopic observations, and also in the Wide-field Infrared Survey Explorer all-sky survey at 12 and 22 mu m. These data were analyzed in detail, and we find that the nebula can be described with a standard bow shock shape, and that its spectrum contains polycyclic aromatic hydrocarbon and H-2 emission features. However, it is not certain which object drives the nebula. We analyze the field stars and conclude that none of them can be the associated object because stars with a strong wind or mass ejection that usually produce bow shocks are much brighter than the field stars. The pulsar is approximately 15 '' away from the region in which the associated object is expected to be located. In order to resolve the discrepancy, we suggest that a highly collimated wind could be emitted from the pulsar and produce the bow shock. X-ray imaging to detect the interaction of the wind with the ambient medium-and high-spatial resolution radio imaging to determine the proper motion of the pulsar should be carried out, which will help verify the association of the pulsar with the bow shock nebula.

A key tracer of the elusive progenitor systems of Type Ia supernovae (SNe Ia) is the detection of narrow blueshifted time-varying Na I D absorption lines, interpreted as evidence of circumstellar material surrounding the progenitor system. The origin of this material is controversial, but the simplest explanation is that it results from previous mass-loss in a system containing a white dwarf and a non-degenerate companion star. We present new single-epoch intermediate-resolution spectra of 17 low-redshift SNe Ia taken with XShooter on the European Southern Observatory Very Large Telescope. Combining this sample with events from the literature, we confirm an excess (similar to 20 per cent) of SNe Ia displaying blueshifted narrow Na I D absorption features compared to redshifted Na I D features. The host galaxies of SNe Ia displaying blueshifted absorption profiles are skewed towards later-type galaxies, compared to SNe Ia that show no Na I D absorption and SNe Ia displaying blueshifted narrow Na I D absorption features have broader light curves. The strength of the Na I D absorption is stronger in SNe Ia displaying blueshifted Na I D absorption features than those without blueshifted features, and the strength of the blueshifted Na I D is correlated with the B - V colour of the SN at maximum light. This strongly suggests the absorbing material is local to the SN. In the context of the progenitor systems of SNe Ia, we discuss the significance of these findings and other recent observational evidence on the nature of SN Ia progenitors. We present a summary that suggests that there are at least two distinct populations of normal, cosmologically useful SNe Ia.

We present the second installment of GOSSS, a massive spectroscopic survey of Galactic O stars, based on new homogeneous, high signal-to-noise ratio, R similar to 2500 digital observations from both hemispheres selected from the Galactic O-Star Catalog (GOSC). In this paper we include bright stars and other objects drawn mostly from the first version of GOSC, all of them south of delta = -20 degrees, for a total number of 258 O stars. We also revise the northern sample of Paper I to provide the full list of spectroscopically classified Galactic O stars complete to B = 8, bringing the total number of published GOSSS stars to 448. Extensive sequences of exceptional objects are given, including the early Of/WN, O Iafpe, Ofc, ON/OC, Onfp, Of?p, and Oe types, as well as double/triple-lined spectroscopic binaries. The new spectral subtype O9.2 is also discussed. The magnitude and spatial distributions of the observed sample are analyzed. We also present new results from OWN, a multi-epoch high-resolution spectroscopic survey coordinated with GOSSS that is assembling the largest sample of Galactic spectroscopic massive binaries ever attained. The OWN data combined with additional information on spectroscopic and visual binaries from the literature indicate that only a very small fraction (if any) of the stars with masses above 15-20M(circle dot) are born as single systems. In the future we will publish the rest of the GOSSS survey, which is expected to include over 1000 Galactic O stars.