Researcher profile

Roberto Trotta

Roberto Trotta contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

5 published item(s)

preprint2026arXiv

Simulation-based population inference of LISA's Galactic binaries: Bypassing the global fit

The Laser Interferometer Space Antenna (LISA) is expected to detect thousands of individually resolved gravitational wave sources, overlapping in time and frequency, on top of unresolved astrophysical and/or primordial backgrounds. Disentangling resolved sources from backgrounds and extracting their parameters in a computationally intensive "global fit" is normally regarded as a necessary step toward reconstructing the properties of the underlying astrophysical populations. Here, we show that it is in principle feasible to infer the population properties of the most numerous of LISA sources -- Galactic double white dwarfs -- directly from the frequency (or, equivalently, time) strain series by adopting a simulation-based approach, without extracting and estimating the parameters of each single source. By training a normalizing flow on a custom-designed compression of simulated LISA frequency series from the Galactic double white dwarf population, we demonstrate how to infer the posterior distribution of population parameters (e.g., mass function, frequency, and spatial distributions). This allows for extracting information on the population parameters from both resolved and unresolved sources simultaneously and in a computationally efficient manner. This approach can be extended to other source classes (e.g., massive and stellar-mass black holes, extreme mass ratio inspirals) and to scenarios involving non-Gaussian or non-stationary noise (e.g., data gaps), provided that fast and accurate simulations are available.

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 robust estimate of the Milky Way mass from rotation curve data

We present a new estimate of the mass of the Milky Way, inferred via a Bayesian approach by making use of tracers of the circular velocity in the disk plane and stars in the stellar halo, as from the publicly available {\tt galkin} compilation. We use the rotation curve method to determine the dark matter distribution and total mass under different assumptions for the dark matter profile, while the total stellar mass is constrained by surface stellar density and microlensing measurements. We also include uncertainties on the baryonic morphology via Bayesian model averaging, thus converting a potential source of systematic error into a more manageable statistical uncertainty. We evaluate the robustness of our result against various possible systematics, including rotation curve data selection, uncertainty on the Sun's velocity $V_0$, dependence on the dark matter profile assumptions, and choice of priors. We find the Milky Way's dark matter virial mass to be $\log_{10}M_{200}^{\rm DM}/ {\rm M_\odot} = 11.92^{+0.06}_{-0.05}{\rm(stat)}\pm{0.28}\pm0.27{\rm(syst)}$ ($M_{200}^{\rm DM}=8.3^{+1.2}_{-0.9}{\rm(stat)}\times10^{11}\,{\rm M_\odot}$). We also apply our framework to Gaia DR2 rotation curve data and find good statistical agreement with the above results.

preprint2020arXiv

STACCATO: A Novel Solution to Supernova Photometric Classification with Biased Training Sets

We present a new solution to the problem of classifying Type Ia supernovae from their light curves alone given a spectroscopically confirmed but biased training set, circumventing the need to obtain an observationally expensive unbiased training set. We use Gaussian processes (GPs) to model the supernovae's (SN) light curves, and demonstrate that the choice of covariance function has only a small influence on the GPs ability to accurately classify SNe. We extend and improve the approach of Richards et al (2012} -- a diffusion map combined with a random forest classifier -- to deal specifically with the case of biassed training sets. We propose a novel method, called STACCATO (SynThetically Augmented Light Curve ClassificATiOn') that synthetically augments a biased training set by generating additional training data from the fitted GPs. Key to the success of the method is the partitioning of the observations into subgroups based on their propensity score of being included in the training set. Using simulated light curve data, we show that STACCATO increases performance, as measured by the area under the Receiver Operating Characteristic curve (AUC), from 0.93 to 0.96, close to the AUC of 0.977 obtained using the 'gold standard' of an unbiased training set and significantly improving on the previous best result of 0.88. STACCATO also increases the true positive rate for SNIa classification by up to a factor of 50 for high-redshift/low brightness SNe.

preprint2020arXiv

Structure Formation Models Weaken Limits on WIMP Dark Matter from Dwarf Spheroidal Galaxies

Dwarf spheroidal galaxies that form in halo substructures provide stringent constraints on dark matter annihilation. Many ultrafaint dwarfs discovered with modern surveys contribute significantly to these constraints. At present, because of the lack of abundant stellar kinematic data for the ultrafaints, non-informative prior assumptions are usually made for the parameters of the density profiles. Based on semi-analytic models of dark matter subhalos and their connection to satellite galaxies, we present more informative and realistic satellite priors. We show that our satellite priors lead to constraints on the annihilation rate that are between a factor of 2 and a factor of 7 weaker than under non-informative priors. As a result, the thermal relic cross section can at best only be excluded (with 95% probability) for dark matter masses of $\lesssim 40$ GeV from dwarf spheroidal data, assuming annihilation into $b\bar{b}$.