We report mass measurements of neutron-rich Ga isotopes $^{80-85}$Ga with TRIUMF's Ion Trap for Atomic and Nuclear science (TITAN). The measurements determine the masses of $^{80-83}$Ga in good agreement with previous measurements. The masses of $^{84}$Ga and $^{85}$Ga were measured for the first time. Uncertainties between $25-48$ keV were reached. The new mass values reduce the nuclear uncertainties associated with the production of A $\approx$ 84 isotopes by the \emph{r}-process for astrophysical conditions that might be consistent with a binary neutron star (BNS) merger producing a blue kilonova. Our nucleosynthesis simulations confirm that BNS merger may contribute to the first abundance peak under moderate neutron-rich conditions with electron fractions $Y_e=0.35-0.38$.

Several anomalous elemental abundance ratios have been observed in the metal-poor star HD94028. We assume that its high [As/Ge] ratio is a product of a weak intermediate (i) neutron-capture process. Given that observational errors are usually smaller than predicted nuclear physics uncertainties, we have first set up a benchmark one-zone i-process nucleosynthesis simulation results of which provide the best fit to the observed abundances. We have then performed Monte Carlo simulations in which 113 relevant (n,$\gamma$) reaction rates of unstable species were randomly varied within Hauser-Feshbach model uncertainty ranges for each reaction to estimate the impact on the predicted stellar abundances. One of the interesting results of these simulations is a double-peaked distribution of the As abundance, which is caused by the variation of the $^{75}$Ga (n,$\gamma$) cross section. This variation strongly anti-correlates with the predicted As abundance, confirming the necessity for improved theoretical or experimental bounds on this cross section. The $^{66}$Ni (n,$\gamma$) reaction is found to behave as a major bottleneck for the i-process nucleosynthesis. Our analysis finds the Pearson product-moment correlation coefficient $r_\mathrm{P} > 0.2$ for all of the i-process elements with $32 \leq Z \leq 42$, with significant changes in their predicted abundances showing up when the rate of this reaction is reduced to its theoretically constrained lower bound. Our results are applicable to any other stellar nucleosynthesis site with the similar i-process conditions, such as Sakurai's object (V4334 Sagittarii) or rapidly-accreting white dwarfs.

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.

Scalar singlet dark matter is one of the simplest and most predictive realisations of the WIMP (weakly-interacting massive particle) idea. Although the model is constrained from all directions by the latest experimental data, it still has viable regions of parameter space. Another compelling aspect of scalar singlets is their ability to stabilise the electroweak vacuum. Indeed, models of scalar dark matter are not low-energy effective theories, but can be valid all the way to the Planck scale. Using the GAMBIT framework, we present the first global fit to include both the low-energy experimental constraints and the theoretical constraints from UV physics, considering models with a scalar singlet charged under either a $\mathbb{Z}_2$ or a $\mathbb{Z}_3$ symmetry. We show that if the model is to satisfy all experimental constraints, completely stabilise the electroweak vacuum up to high scales, and also remain perturbative to those scales, one is driven to a relatively small region of parameter space. This region has a Higgs-portal coupling slightly less than 1, a dark matter mass of 1 to 2 TeV and a spin-independent nuclear scattering cross-section around 10$^{-45}$ cm$^2$.

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.

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.

The extension of the standard model (SM) by a Higgs portal scalar field is one the simplest dark matter theories. We present here the first results for a global fit to this model using the global and beyond the SM inference tool (GAMBIT). This software enables the combination of dark matter constraints in a statistically consistent manner. In total 15 parameters are varied and the parameter space explored using four different scanning algorithms. The viable parameter space is reduced from previous studies of this model due to the inclusion of the latest direct detection constraints.

We characterize a cosmic rest frame in which the monopole variation of the spherically averaged nonlinear Hubble expansion is most uniform, under arbitrary local Lorentz boosts of the central observer. Using the COMPOSITE sample of 4534 galaxies, we identify a degenerate set of candidate minimum nonlinear variation frames, which includes the rest frame of the Local Group (LG) of galaxies, but excludes the standard Cosmic Microwave Background (CMB) frame. Candidate rest frames defined by a boost from the LG frame close to the plane of the galaxy have a statistical likelihood similar to the LG frame. This may result from a lack of constraining data in the Zone of Avoidance. We extend our analysis to the Cosmicflows-2 (CF2) sample of 8162 galaxies. While the signature of a systematic boost offset between the CMB and LG frame averages is still detected, the spherically averaged nonlinear expansion variation in all rest frames is significantly larger in the CF2 sample than would be reasonably expected. We trace this to the CF2 distances being reported without a correction for inhomogeneous distribution Malmquist bias. Systematic differences in the inclusion of the large SFI++ subsample into the COMPOSITE and CF2 catalogues are analysed. Our results highlight the importance of a careful treatment of Malmquist biases for future peculiar velocities studies, including tests of the hypothesis of Wiltshire et al [Phys. Rev. D 88 (2013) 083529; arXiv:1201.5371] that a significant fraction of the CMB temperature dipole may be nonkinematic in origin.