Source author record

Gregory D. Martinez

Gregory D. Martinez appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

16 published item(s)

preprint2022arXiv

Simple and statistically sound recommendations for analysing physical theories

Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by multiple experiments, and random or grid sampling of model parameters. Whilst these methods are easy to apply, they exhibit pathologies even in low-dimensional parameter spaces, and quickly become problematic to use and interpret in higher dimensions. In this article we give clear guidance for going beyond these procedures, suggesting where possible simple methods for performing statistically sound inference, and recommendations of readily-available software tools and standards that can assist in doing so. Our aim is to provide any physicists lacking comprehensive statistical training with recommendations for reaching correct scientific conclusions, with only a modest increase in analysis burden. Our examples can be reproduced with the code publicly available at https://doi.org/10.5281/zenodo.4322283.

preprint2020arXiv

A population of dust-enshrouded objects orbiting the Galactic black hole

The central 0.1 parsecs of the Milky Way host a supermassive black hole identified with the position of the radio and infrared source Sagittarius A*, a cluster of young, massive stars (the S stars) and various gaseous features. Recently, two unusual objects have been found to be closely orbiting Sagittarius A*: the so-called G sources, G1 and G2. These objects are unresolved (having a size of the order of 100 astronomical units, except at periapse, where the tidal interaction with the black hole stretches them along the orbit) and they show both thermal dust emission and line emission from ionized gas. G1 and G2 have generated attention because they appear to be tidally interacting with the supermassive Galactic black hole, possibly enhancing its accretion activity. No broad consensus has yet been reached concerning their nature: the G objects show the characteristics of gas and dust clouds but display the dynamical properties of stellar-mass objects. Here we report observations of four additional G objects, all lying within 0.04 parsecs of the black hole and forming a class that is probably unique to this environment. The widely varying orbits derived for the six G objects demonstrate that they were commonly but separately formed.

preprint2020arXiv

Global analyses of Higgs portal singlet dark matter models using GAMBIT

We present global analyses of effective Higgs portal dark matter models in the frequentist and Bayesian statistical frameworks. Complementing earlier studies of the scalar Higgs portal, we use GAMBIT to determine the preferred mass and coupling ranges for models with vector, Majorana and Dirac fermion dark matter. We also assess the relative plausibility of all four models using Bayesian model comparison. Our analysis includes up-to-date likelihood functions for the dark matter relic density, invisible Higgs decays, and direct and indirect searches for weakly-interacting dark matter including the latest XENON1T data. We also account for important uncertainties arising from the local density and velocity distribution of dark matter, nuclear matrix elements relevant to direct detection, and Standard Model masses and couplings. In all Higgs portal models, we find parameter regions that can explain all of dark matter and give a good fit to all data. The case of vector dark matter requires the most tuning and is therefore slightly disfavoured from a Bayesian point of view. In the case of fermionic dark matter, we find a strong preference for including a CP-violating phase that allows suppression of constraints from direct detection experiments, with odds in favour of CP violation of the order of 100:1. Finally, we present DDCalc 2.0.0, a tool for calculating direct detection observables and likelihoods for arbitrary non-relativistic effective operators.

preprint2019arXiv

Combined collider constraints on neutralinos and charginos

Searches for supersymmetric electroweakinos have entered a crucial phase, as the integrated luminosity of the Large Hadron Collider is now high enough to compensate for their weak production cross-sections. Working in a framework where the neutralinos and charginos are the only light sparticles in the Minimal Supersymmetric Standard Model, we use gambit to perform a detailed likelihood analysis of the electroweakino sector. We focus on the impacts of recent ATLAS and CMS searches with 36 fb$^{-1}$ of 13 TeV proton-proton collision data. We also include constraints from LEP and invisible decays of the $Z$ and Higgs bosons. Under the background-only hypothesis, we show that current LHC searches do not robustly exclude any range of neutralino or chargino masses. However, a pattern of excesses in several LHC analyses points towards a possible signal, with neutralino masses of $(m_{\tildeχ_1^0}, m_{\tildeχ_2^0}, m_{\tildeχ_3^0}, m_{\tildeχ_4^0})$ = (8-155, 103-260, 130-473, 219-502) GeV and chargino masses of $(m_{\tildeχ_1^{\pm}}, m_{\tildeχ_2^{\pm}})$ = (104-259, 224-507) GeV at the 95% confidence level. The lightest neutralino is mostly bino, with a possible modest Higgsino or wino component. We find that this excess has a combined local significance of $3.3σ$, subject to a number of cautions. If one includes LHC searches for charginos and neutralinos conducted with 8 TeV proton-proton collision data, the local significance is lowered to 2.9$σ$. We briefly consider the implications for dark matter, finding that the correct relic density can be obtained through the Higgs-funnel and $Z$-funnel mechanisms, even assuming that all other sparticles are decoupled. All samples, gambit input files and best-fit models from this study are available on Zenodo.

preprint2018arXiv

A global fit of the MSSM with GAMBIT

We study the seven-dimensional Minimal Supersymmetric Standard Model (MSSM7) with the new GAMBIT software framework, with all parameters defined at the weak scale. Our analysis significantly extends previous weak-scale, phenomenological MSSM fits, by adding more and newer experimental analyses, improving the accuracy and detail of theoretical predictions, including dominant uncertainties from the Standard Model, the Galactic dark matter halo and the quark content of the nucleon, and employing novel and highly-efficient statistical sampling methods to scan the parameter space. We find regions of the MSSM7 that exhibit co-annihilation of neutralinos with charginos, stops and sbottoms, as well as models that undergo resonant annihilation via both light and heavy Higgs funnels. We find high-likelihood models with light charginos, stops and sbottoms that have the potential to be within the future reach of the LHC. Large parts of our preferred parameter regions will also be accessible to the next generation of direct and indirect dark matter searches, making prospects for discovery in the near future rather good.

preprint2018arXiv

Global fits of GUT-scale SUSY models with GAMBIT

We present the most comprehensive global fits to date of three supersymmetric models motivated by grand unification: the Constrained Minimal Supersymmetric Standard Model (CMSSM), and its Non-Universal Higgs Mass generalisations NUHM1 and NUHM2. We include likelihoods from a number of direct and indirect dark matter searches, a large collection of electroweak precision and flavour observables, direct searches for supersymmetry at LEP and Runs I and II of the LHC, and constraints from Higgs observables. Our analysis improves on existing results not only in terms of the number of included observables, but also in the level of detail with which we treat them, our sampling techniques for scanning the parameter space, and our treatment of nuisance parameters. We show that stau co-annihilation is now ruled out in the CMSSM at more than 95\% confidence. Stop co-annihilation turns out to be one of the most promising mechanisms for achieving an appropriate relic density of dark matter in all three models, whilst avoiding all other constraints. We find high-likelihood regions of parameter space featuring light stops and charginos, making them potentially detectable in the near future at the LHC. We also show that tonne-scale direct detection will play a largely complementary role, probing large parts of the remaining viable parameter space, including essentially all models with multi-TeV neutralinos.

preprint2018arXiv

Status of the scalar singlet dark matter model

One of the simplest viable models for dark matter is an additional neutral scalar, stabilised by a $\mathbb{Z}_2$ symmetry. Using the GAMBIT package and combining results from four independent samplers, we present Bayesian and frequentist global fits of this model. We vary the singlet mass and coupling along with 13 nuisance parameters, including nuclear uncertainties relevant for direct detection, the local dark matter density, and selected quark masses and couplings. We include the dark matter relic density measured by Planck, direct searches with LUX, PandaX, SuperCDMS and XENON100, limits on invisible Higgs decays from the Large Hadron Collider, searches for high-energy neutrinos from dark matter annihilation in the Sun with IceCube, and searches for gamma rays from annihilation in dwarf galaxies with the Fermi-LAT. Viable solutions remain at couplings of order unity, for singlet masses between the Higgs mass and about 300 GeV, and at masses above $\sim$1 TeV. Only in the latter case can the scalar singlet constitute all of dark matter. Frequentist analysis shows that the low-mass resonance region, where the singlet is about half the mass of the Higgs, can also account for all of dark matter, and remains viable. However, Bayesian considerations show this region to be rather fine-tuned.

preprint2017arXiv

Comparison of statistical sampling methods with ScannerBit, the GAMBIT scanning module

We introduce ScannerBit, the statistics and sampling module of the public, open-source global fitting framework GAMBIT. ScannerBit provides a standardised interface to different sampling algorithms, enabling the use and comparison of multiple computational methods for inferring profile likelihoods, Bayesian posteriors, and other statistical quantities. The current version offers random, grid, raster, nested sampling, differential evolution, Markov Chain Monte Carlo (MCMC) and ensemble Monte Carlo samplers. We also announce the release of a new standalone differential evolution sampler, Diver, and describe its design, usage and interface to ScannerBit. We subject Diver and three other samplers (the nested sampler MultiNest, the MCMC GreAT, and the native ScannerBit implementation of the ensemble Monte Carlo algorithm T-Walk) to a battery of statistical tests. For this we use a realistic physical likelihood function, based on the scalar singlet model of dark matter. We examine the performance of each sampler as a function of its adjustable settings, and the dimensionality of the sampling problem. We evaluate performance on four metrics: optimality of the best fit found, completeness in exploring the best-fit region, number of likelihood evaluations, and total runtime. For Bayesian posterior estimation at high resolution, T-Walk provides the most accurate and timely mapping of the full parameter space. For profile likelihood analysis in less than about ten dimensions, we find that Diver and MultiNest score similarly in terms of best fit and speed, outperforming GreAT and T-Walk; in ten or more dimensions, Diver substantially outperforms the other three samplers on all metrics.

preprint2017arXiv

GAMBIT: The Global and Modular Beyond-the-Standard-Model Inference Tool

We describe the open-source global fitting package GAMBIT: the Global And Modular Beyond-the-Standard-Model Inference Tool. GAMBIT combines extensive calculations of observables and likelihoods in particle and astroparticle physics with a hierarchical model database, advanced tools for automatically building analyses of essentially any model, a flexible and powerful system for interfacing to external codes, a suite of different statistical methods and parameter scanning algorithms, and a host of other utilities designed to make scans faster, safer and more easily-extendible than in the past. Here we give a detailed description of the framework, its design and motivation, and the current models and other specific components presently implemented in GAMBIT. Accompanying papers deal with individual modules and present first GAMBIT results. GAMBIT can be downloaded from gambit.hepforge.org.

preprint2015arXiv

A Robust Determination of Milky Way Satellite Properties using Hierarchical Mass Modeling

We introduce a new methodology to robustly determine the mass profile, as well as the overall distribution, of Local Group satellite galaxies. Specifically we employ a statistical multilevel modelling technique, Bayesian hierarchical modelling, to simultaneously constrain the properties of individual Local Group Milky Way satellite galaxies and the characteristics of the Milky Way satellite population. We show that this methodology reduces the uncertainty in individual dwarf galaxy mass measurements up to a factor of a few for the faintest galaxies. We find that the distribution of Milky Way satellites inferred by this analysis, with the exception of the apparent lack of high-mass haloes, is consistent with the Lambda cold dark matter (Lambda-CDM) paradigm. In particular we find that both the measured relationship between the maximum circular velocity and the radius at this velocity, as well as the inferred relationship between the mass within 300 pc and luminosity, match the values predicted by Lambda-CDM simulations for halos with maximum circular velocities below 20 km/sec. Perhaps more striking is that this analysis seems to suggest a more cusped "average" halo shape that is shared by these galaxies. While this study reconciles many of the observed properties of the Milky Way satellite distribution with that of Lambda-CDM simulations, we find that there is still a deficit of satellites with maximum circular velocities of 20-40 km/sec.

preprint2013arXiv

Stellar Populations in the Central 0.5 pc of the Galaxy I: A New Method for Constructing Luminosity Functions and Surface-Density Profiles

We present new high angular resolution near-infrared spectroscopic observations of the nuclear star cluster surrounding the Milky Way's central supermassive black hole. Using the integral-field spectrograph OSIRIS on Keck II behind the laser-guide-star adaptive optics system, this spectroscopic survey enables us to separate early-type (young, 4-6 Myr) and late-type (old, >1 Gyr) stars with a completeness of 50% down to K' = 15.5 mag, which corresponds to ~10 \msun for the early-type stars. This work increases the radial extent of reported OSIRIS/Keck measurements by more than a factor of 3 from 4" to 14" (0.16 pc to 0.56 pc), along the projected disk of young stars. For our analysis, we implement a new method of completeness correction using a combination of star-planting simulations and Bayesian inference. We assign probabilities for the spectral type of every source detected in deep imaging down to K' = 15.5 mag using information from spectra, simulations, number counts, and the distribution of stars. The inferred radial surface-density profiles, $Σ(R) \propto R^{-Γ}$, for the young stars and late-type giants are consistent with earlier results ($Γ_{early} = 0.93 \pm 0.09$, $Γ_{late} = 0.16 \pm 0.07$). The late-type surface-density profile is approximately flat out to the edge of the survey. While the late-type stellar luminosity function is consistent with the Galactic bulge, the completeness-corrected luminosity function of the early-type stars has significantly more young stars at faint magnitudes compared to previous surveys with similar depth. This luminosity function indicates that the corresponding mass function of the young stars is likely less top-heavy than that inferred from previous surveys.

preprint2012arXiv

Evidence for Substructure in Ursa Minor Dwarf Spheroidal Galaxy using a Bayesian Object Detection Method

We present a method for identifying localized secondary populations in stellar velocity data using Bayesian statistical techniques. We apply this method to the dwarf spheroidal galaxy Ursa Minor and find two secondary objects in this satellite of the Milky Way. One object is kinematically cold with a velocity dispersion of $4.25 \pm 0.75\ \kms$ and centered at $(9.1\arcmin \pm 1.5, 7.2\arcmin \pm 1.2)$ in relative RA and DEC with respect to the center of Ursa Minor. The second object has a large velocity offset of $-12.8^{+1.75}_{-1.5}\ \kms$ compared to Ursa Minor and centered at $(-14.0\arcmin^{+2.4}_{-5.8}, -2.5\arcmin^{+0.4}_{-1.0})$. The kinematically cold object has been found before using a smaller data set but the prediction that this cold object has a velocity dispersion larger than $2.0\ \kms$ at 95% C.L. differs from previous work. We use two and three component models along with the information criteria and Bayesian evidence model selection methods to argue that Ursa Minor has one or two localized secondary populations. The significant probability for a large velocity dispersion in each secondary object raises the intriguing possibility that each has its own dark matter halo, that is, it is a satellite of a satellite of the Milky Way.

preprint2011arXiv

A Complete Spectroscopic Survey of the Milky Way satellite Segue 1: Dark matter content, stellar membership and binary properties from a Bayesian analysis

We introduce a comprehensive analysis of multi-epoch stellar line-of-sight velocities to determine the intrinsic velocity dispersion of the ultrafaint satellites of the Milky Way. Our method includes a simultaneous Bayesian analysis of both membership probabilities and the contribution of binary orbital motion to the observed velocity dispersion within a 14-parameter likelihood. We apply our method to the Segue 1 dwarf galaxy and conclude that Segue 1 is a dark-matter-dominated galaxy at high probability with an intrinsic velocity dispersion of 3.7^{+1.4}_{-1.1} km/sec. The dark matter halo required to produce this dispersion must have an average density of 2.5^{+4.1}_{-1.9} solar mass/pc^3 within a sphere that encloses half the galaxy's stellar luminosity. This is the highest measured density of dark matter in the Local Group. Our results show that a significant fraction of the stars in Segue 1 may be binaries with the most probable mean period close to 10 years, but also consistent with the 180 year mean period seen in the solar vicinity at about 1 sigma. Despite this binary population, the possibility that Segue 1 is a bound star cluster with the observed velocity dispersion arising from the orbital motion of binary stars is disfavored by the multi-epoch stellar velocity data at greater than 99% C.L. Finally, our treatment yields a projected (two-dimensional) half-light radius for the stellar profile of Segue 1 of 28^{+5}_{-4} pc, in excellent agreement with photometric measurements.

preprint2011arXiv

A Complete Spectroscopic Survey of the Milky Way Satellite Segue 1: The Darkest Galaxy

We present the results of a comprehensive Keck/DEIMOS spectroscopic survey of the ultra-faint Milky Way satellite galaxy Segue 1. We have obtained velocity measurements for 98.2% of the stars within 67 pc (10 arcmin, or 2.3 half-light radii) of the center of Segue 1 that have colors and magnitudes consistent with membership, down to a magnitude limit of r=21.7. Based on photometric, kinematic, and metallicity information, we identify 71 stars as probable Segue 1 members, including some as far out as 87 pc. After correcting for the influence of binary stars using repeated velocity measurements, we determine a velocity dispersion of 3.7^{+1.4}_{-1.1} km/s, with a corresponding mass within the half-light radius of 5.8^{+8.2}_{-3.1} x 10^5 Msun. The stellar kinematics of Segue 1 require very high mass-to-light ratios unless the system is far from dynamical equilibrium, even if the period distribution of unresolved binary stars is skewed toward implausibly short periods. With a total luminosity less than that of a single bright red giant and a V-band mass-to-light ratio of 3400 Msun/Lsun, Segue 1 is the darkest galaxy currently known. We critically re-examine recent claims that Segue 1 is a tidally disrupting star cluster and that kinematic samples are contaminated by the Sagittarius stream. The extremely low metallicities ([Fe/H] < -3) of two Segue 1 stars and the large metallicity spread among the members demonstrate conclusively that Segue 1 is a dwarf galaxy, and we find no evidence in favor of tidal effects. We also show that contamination by the Sagittarius stream has been overestimated. Segue 1 has the highest measured dark matter density of any known galaxy and will therefore be a prime testing ground for dark matter physics and galaxy formation on small scales.

preprint2010arXiv

Accurate masses for dispersion-supported galaxies

We derive an accurate mass estimator for dispersion-supported stellar systems and demonstrate its validity by analyzing resolved line-of-sight velocity data for globular clusters, dwarf galaxies, and elliptical galaxies. Specifically, by manipulating the spherical Jeans equation we show that the dynamical mass enclosed within the 3D deprojected half-light radius r_1/2 can be determined with only mild assumptions about the spatial variation of the stellar velocity dispersion anisotropy. We find M_1/2 = 3 σ_los^2 r_1/2 / G ~ 4 σ_los^2 R_eff / G, where σ_los^2 is the luminosity-weighted square of the line-of-sight velocity dispersion and R_eff is the 2D projected half-light radius. While deceptively familiar in form, this formula is not the virial theorem, which cannot be used to determine accurate masses unless the radial profile of the total mass is known a priori. We utilize this finding to show that all of the Milky Way dwarf spheroidal galaxies (MW dSphs) are consistent with having formed within a halo of mass approximately 3 x 10^9 M_sun in Lambda CDM cosmology. The faintest MW dSphs seem to have formed in dark matter halos that are at least as massive as those of the brightest MW dSphs, despite the almost five orders of magnitude spread in luminosity. We expand our analysis to the full range of observed dispersion-supported stellar systems and examine their I-band mass-to-light ratios (M/L). The M/L vs. M_1/2 relation for dispersion-supported galaxies follows a U-shape, with a broad minimum near M/L ~ 3 that spans dwarf elliptical galaxies to normal ellipticals, a steep rise to M/L ~ 3,200 for ultra-faint dSphs, and a more shallow rise to M/L ~ 800 for galaxy cluster spheroids.

preprint2010arXiv

The Case Against Warm or Self-Interacting Dark Matter as Explanations for Cores in Low Surface Brightness Galaxies

Warm dark matter (WDM) and self-interacting dark matter (SIDM) are often motivated by the inferred cores in the dark matter halos of low surface brightness (LSB) galaxies. We test thermal WDM, non-thermal WDM, and SIDM using high-resolution rotation curves of nine LSB galaxies. We fit these dark matter models to the data and determine the halo core radii and central densities. While the minimum core size in WDM models is predicted to decrease with halo mass, we find that the inferred core radii increase with halo mass and also cannot be explained with a single value of the primordial phase space density. Moreover, if the core size is set by WDM particle properties, then even the smallest cores we infer would require primordial phase space density values that are orders of magnitude smaller than lower limits obtained from the Lyman alpha forest power spectra. We also find that the dark matter halo core densities vary by a factor of about 30 from system to system while showing no systematic trend with the maximum rotation velocity of the galaxy. This strongly argues against the core size being directly set by large self-interactions (scattering or annihilation) of dark matter. We therefore conclude that the inferred cores do not provide motivation to prefer WDM or SIDM over other dark matter models.