Globular Cluster Systems in Brightest Cluster Galaxies. III: Beyond Bimodality

We present new deep photometry of the rich globular cluster (GC) systems around the Brightest Cluster Galaxies UGC 9799 (Abell 2052) and UGC 10143 (Abell 2147), obtained with the HST ACS and WFC3 cameras. For comparison, we also present new reductions of similar HST/ACS data for the Coma supergiants NGC 4874 and 4889. All four of these galaxies have huge cluster populations (to the radial limits of our data, comprising from 12000 to 23000 clusters per galaxy). The metallicity distribution functions (MDFs) of the GCs can still be matched by a bimodal-Gaussian form where the metal-rich and metal-poor modes are separated by ~0.8 dex, but the internal dispersions of each mode are so large that the total MDF becomes very broad and nearly continuous from [Fe/H] = -2.4 to Solar. There are, however, significant differences between galaxies in the relative numbers of \emph{metal-rich} clusters, suggesting that they underwent significantly different histories of mergers with massive, gas-rich halos. Lastly, the proportion of metal-poor GCs rises especially rapidly outside projected radii R > 4 R_eff, suggesting the importance of accreted dwarf satellites in the outer halo. Comprehensive models for the formation of GCs as part of the hierarchical formation of their parent galaxies will be needed to trace the systematic change in structure of the MDF with galaxy mass, from the distinctly bimodal form in smaller galaxies up to the broad continuum that we see in the very largest systems.Comment: In press for Astrophysical Journa

Maximum likelihood estimation of photon number distribution from homodyne statistics

We present a method for reconstructing the photon number distribution from the homodyne statistics based on maximization of the likelihood function derived from the exact statistical description of a homodyne experiment. This method incorporates in a natural way the physical constraints on the reconstructed quantities, and the compensation for the nonunit detection efficiency.Comment: 3 pages REVTeX. Final version, to appear in Phys. Rev. A as a Brief Repor

The light curve of a transient X-ray source

The Ariel-V satellite monitored the X-ray light curve of A1524-62 almost continuously from 40 days prior to maximum light until its disappearance below the effective experimental sensitivity. The source exhibited maximum light on approximately 4 December 1974, at a level of 0.9 the apparent magnitude of the Crab Nebula in the energy band 3-6 keV. Although similar to previously reported transient sources with a decay time constant of approximately 2 months, the source exhibited an extended, variable pre-flare on-state of about 1 month at a level of greater than approximately 0.1 maximum light. The four bright (greater than 0.2 of the Crab Nebula) transient sources observed during the first half-year of Ariel-V operation are indicative of a galactic disk distribution, and a luminosity at maximum in excess of 10 to the 37th power ergs/sec

Poisson Cluster Process Models for Detecting Ultra-Diffuse Galaxies

We propose a novel set of Poisson Cluster Process (PCP) models to detect Ultra-Diffuse Galaxies (UDGs), a class of extremely faint, enigmatic galaxies of substantial interest in modern astrophysics. We model the unobserved UDG locations as parent points in a PCP, and infer their positions based on the observed spatial point patterns of their old star cluster systems. Many UDGs have somewhere from a few to hundreds of these old star clusters, which we treat as offspring points in our models. We also present a new framework to construct a marked PCP model using the marks of star clusters. The marked PCP model may enhance the detection of UDGs and offers broad applicability to problems in other disciplines. To assess the overall model performance, we design an innovative assessment tool for spatial prediction problems where only point-referenced ground truth is available, overcoming the limitation of standard ROC analyses where spatial Boolean reference maps are required. We construct a bespoke blocked Gibbs adaptive spatial birth-death-move MCMC algorithm to infer the locations of UDGs using real data from a \textit{Hubble Space Telescope} imaging survey. Based on our performance assessment tool, our novel models significantly outperform existing approaches using the Log-Gaussian Cox Process. We also obtained preliminary evidence that the marked PCP model improves UDG detection performance compared to the model without marks. Furthermore, we find evidence of a potential new ``dark galaxy'' that was not detected by previous methods.Comment: 47 pages, 8 figures, 3 tables; submitted to AoAS, comments are welcom

Comparison between methods for the determination of the primary cosmic ray mass composition from the longitudinal profile of atmospheric cascades

The determination of the primary cosmic ray mass composition from the longitudinal development of atmospheric cascades is still a debated issue. In this work we discuss several data analysis methods and show that if the entire information contained in the longitudinal profile is exploited, reliable results may be obtained. Among the proposed methods FCC ('Fit of the Cascade Curve'), MTA ('Multiparametric Topological Analysis') and NNA ('Neural Net Analysis') with conjugate gradient optimization algorithm give the best accuracy.Comment: 22 pages, 11 figures, accepted by Astroparticle Physics, minor misprints and an extra figure remove

A combined analysis of short-baseline neutrino experiments in the (3+1) and (3+2) sterile neutrino oscillation hypotheses

We investigate adding two sterile neutrinos to resolve the apparent tension existing between short-baseline neutrino oscillation results and CPT-conserving, four-neutrino oscillation models. For both (3+1) and (3+2) models, the level of statistical compatibility between the combined dataset from the null short-baseline experiments Bugey, CHOOZ, CCFR84, CDHS, KARMEN, and NOMAD, on the one hand; and the LSND dataset, on the other, is computed. A combined analysis of all seven short-baseline experiments, including LSND, is also performed, to obtain the favored regions in neutrino mass and mixing parameter space for both models. Finally, four statistical tests to compare the (3+1) and the (3+2) hypotheses are discussed. All tests show that (3+2) models fit the existing short-baseline data significantly better than (3+1) models.Comment: 16 pages, 15 figures. Added NOMAD data to the analysis, one statistical test, and two figures. References and text added. Version submitted to PR

Quantum homodyne tomography with a priori constraints

I present a novel algorithm for reconstructing the Wigner function from homodyne statistics. The proposed method, based on maximum-likelihood estimation, is capable of compensating for detection losses in a numerically stable way.Comment: 4 pages, REVTeX, 2 figure

Recent AEM Case Study Examples of a Full Waveform Time-Domain System for Near-Surface and Groundwater Applications

Early time or high frequency airborne electromagnetic data (AEM) are desirable for shallow sounding or mapping of resistive areas but this poses difficulties due to a variety of issues, such as system bandwidth, system calibration and parasitic loop capacitance. In an effort to address this issue, a continued system design strategy, aimed at improving its early-channel VTEM data, has achieved fully calibrated, quantitative measurements closer to the transmitter current turn-off, while maintaining reasonably optimal deep penetration characteristics. The new design implementation, known as “Full Waveform” VTEM was previously described by Legault et al. (2012). This paper presents some case-study examples of a Full Waveform helicopter time-domain EM system for near-surface application

Geo-neutrinos: A systematic approach to uncertainties and correlations

Geo-neutrinos emitted by heat-producing elements (U, Th and K) represent a unique probe of the Earth interior. The characterization of their fluxes is subject, however, to rather large and highly correlated uncertainties. The geochemical covariance of the U, Th and K abundances in various Earth reservoirs induces positive correlations among the associated geo-neutrino fluxes, and between these and the radiogenic heat. Mass-balance constraints in the Bulk Silicate Earth (BSE) tend instead to anti-correlate the radiogenic element abundances in complementary reservoirs. Experimental geo-neutrino observables may be further (anti)correlated by instrumental effects. In this context, we propose a systematic approach to covariance matrices, based on the fact that all the relevant geo-neutrino observables and constraints can be expressed as linear functions of the U, Th and K abundances in the Earth's reservoirs (with relatively well-known coefficients). We briefly discuss here the construction of a tentative "geo-neutrino source model" (GNSM) for the U, Th, and K abundances in the main Earth reservoirs, based on selected geophysical and geochemical data and models (when available), on plausible hypotheses (when possible), and admittedly on arbitrary assumptions (when unavoidable). We use then the GNSM to make predictions about several experiments ("forward approach"), and to show how future data can constrain - a posteriori - the error matrix of the model itself ("backward approach"). The method may provide a useful statistical framework for evaluating the impact and the global consistency of prospective geo-neutrino measurements and Earth models.Comment: 17 pages, including 4 figures. To appear on "Earth, Moon, and Planets," Special Issue on "Neutrino Geophysics," Proceedings of Neutrino Science 2005 (Honolulu, Hawaii, Dec. 2005