Recent news about our research group:


2016/11/06: "La vida intima de las particulas" ("The inner life of particles"), my first popular science book about high energy physics, is now available (in Spanish).

2016/11/01: New publication: "Precision determination of the small-x gluon from charm production at LHCb"

2016/10/25: The 4th instalment of the "Handbook of LHC Cross-Sections", "Deciphering the Higgs sector", is now available on the arXiv. The updated cross-sections have been computed using the new PDF4LHC working group recommendations for Run II.

My research group is currently composed of the following members:

Postdoctoral researchers

  • Valerio Bertone (Oxford and VU/Nikhef, from 2015)
  • Nathan P. Hartland (Oxford and VU/Nikhef, from 2014)
  • Emanuele R. Nocera (Oxford, from 2015)

PhD students

  • Francesco Giuli (Oxford, from 2014, co-supervised with Amanda Cooper-Sarkar)
  • Luca Rottoli (Oxford, from 2014)
  • Emma Slade (Oxford, from 2016)

Previous members of my research group, with their next destination

  • Marco Bonvini (in Oxford, postdoc between 2014 and 2016), now postdoc at Roma I

Here I provide some background material and references about my research interests.

With the discovery of the Higgs boson at the LHC, awarded with the Nobel prize in 2013, particle physics has entered a new era. The Higgs boson, responsible for giving their mass to all known particles, is the last missing element of the extremely successful Standard Model (SM), but at the same time its discovery raises a number of new fascinating questions that need to be addressed. The challenge for particle physics in the coming years is to understand in detail the properties of this new particle, the first ever fundamental scalar: any deviation with respect to the SM predictions would be a smoking gun for New Physics beyond the SM (BSM). In addition, the LHC will continue the exploration of the high energy frontier with the search for exotic heavy particles, new forces or additional space-time dimensions. The LHC program has also profound implications for exciting open problems in astronomy and cosmology, for example, many scenarios predict the production of dark matter particles at the LHC.

Parton Distribution Functions

It is therefore of utmost importance to provide accurate theoretical predictions for the most relevant processes at the LHC, including Higgs production, as well as for a variety of New Physics scenarios. Crucial ingredients of these predictions are the Parton Distribution Functions (PDFs) of the proton, which encode the dynamics determining how the energy of the protons that collide at the LHC is split among its fundamental constituents, the quarks and gluons. Being intrinsically non-perturbative, PDFs cannot be computed from first principles, and thus need to be extracted from experimental data, in the context of the so-called global QCD analysis framework.

Our limited knowledge of parton distributions degrades the scientific output of many important LHC analyses. To begin with, PDFs are one of the dominant theory systematics for Higgs production cross-sections, limiting the best precision that can be achieved the extraction of Higgs couplings. PDF uncertainties also affect high-mass particle production at the LHC, ubiquitous in BSM scenarios: for example, the cross-section for the production of a pair of the supersymmetric parters of the gluon, the gluinos, can be affected by up to 100% uncertainties. Moreover, PDF uncertainties are also the largest theory uncertainty in precision SM measurements such as the determination of the W mass at hadron colliders.

I am among the developers of a novel approach to the determination of PDFs based on Artificial Neural Networks (ANN) and Maching Learning techniques, the so-called NNPDF approach. In this framework, PDFs learn the underlying physical law from experimental data without the need of imposing any prior knowledge. The ANN training is performed using Genetic Algorithms, non-deterministic minimization strategies suitable for the solution of complex optimization problems, for instance when a very large number of quasi-equivalent minima are present. Non-linear propagation of uncertainties from the fitted data to the PDFs is achieved by means of Monte Carlo replicas, in contrast to the usual Hessian method whose validity is restricted to the gaussian (linear) regime of PDF uncertainties.

One recent breakthrough in PDF analysis has been the constraints of LHC Run I data in PDF fits. I have studied the inclusion of LHC data into PDF fits, from direct photons to top quark pair production cross-sections and charmed meson differential distributions. LHC data is now a central ingredient of PDF fits, providing vital information on poorly-known PDFs such as the large and small-x gluon or the large-x antiquarks. I have also explored novel ways to exploit the LHC data to provide PDF information, for instance with the proposal of using the ratios of cross-sections between different center of mass energies, which benefits from high PDF sensitivity and reduced experimental and theory systematics due to cancellations in the ratio.

In the figure below I show a snapshot of our current understanding of the proton structure in terms of unpolarized (NNPDF3.0, left plot) and polarized (NNPDFpol1.1, right plot) parton distributuons, from the Particle Data Group review.

Unpolarized and polarized PDFs

Higgs pair production and jet substructure.

Another of my research lines is the exploitation of perturbative QCD calculations for the internal structure of jets in hadronic collisions. Jet substructure methods allow to disentangle QCD jets from jets arising from the decays of heavy resonances, and at the LHC they are crucial to optimize the exploration of the Higgs sector and in searches for BSM physics. In the paradigmatic example, jet substructure methods can turn Higgs associate production in the bbfinal state into a discovery channel despite the overwhelming QCD background. At the LHC, many searches for heavy new physics at the TeV scale rely heavily on jet substructure methods.

One of the cornerstones of the LHC program for the next years will be Higgs pair production, a process for which jet substructure techniques are crucial, in particular when one of the Higgs (or the two) decay into bbpairs. The measurement of the production of Higgs pairs is necessary in order to be able to reconstruct the entire electroweak symmetry breaking potential. While in single Higgs production we can only access the local region around the minimum, in double Higgs production we can reconstruct the full Higgs sector potential. It is difficult to underestimate the importance of this process: not only it would provide the first evidence that Higgs bosons have self-interactions, but also it allows to probe the more ad-hoc sector of the Standard Model, the Higgs potential.  Therefore, measuring Higgs pair production could be the key to unveil the origin of the mechanism of electroweak symmetry breaking and open a new window to physics beyond the SM.

Higgs pair production

My studies have concentrated on the 4b final state, characterised by the higher signal yields, but where the huge QCD background needs to be reduced by several orders of magnitude. For the analysis of the resonant case, where di-Higgs production is mediated by a heavy resonance, as in scenarios with extra dimensions, We developed the scale-invariant resonance tagging method. The basic idea of this technique is the consistent combination of all possible final-state topologies, boosted, resolved and intermediate, into a common optimised analysis. This study allowed to probe a wide range of the BSM parameter space, in particular from low to very high radion and graviton masses.

The same techniques have also been used in a feasibility study which demonstrates how, using the 4b final state, it might be possible to measure the Higgs self-coupling already at Run II, earlier than what previous estimates suggested. This result was achieved thanks to the use of advanced Multivariate Analysis (MVA) techniques. In particular we used Artificial Neural Networks (ANN) to enhance the discrimination of the signal over background events. The MVA is trained on the kinematics of signal and background events (including jet substructure variables), allowing to derive a non-trivial set of highly-correlated cuts that lead to a high signal significance at the High-Luminosity LHC.

Astrophysical implications of collider measurements.

I believe that it is always fascinating to cross boundaries between research fields, and I am always keen to explore the implications of high-precision QCD calculations beyond the collider framework. My first excursion into astro-particle physics was a direct determination of the conventional neutrino flux from Super-Kamiokande data using Artificial Neural Networks, applying in a completely new context a methodology that was being developed for proton PDF fits. This calculation has been widely used by the neutrino community.

More recently, we have completed an updated calculation of the prompt neutrino flux at IceCube. The IceCube experiment at the South Pole claimed in 2013 evidence of the first ever detection of astrophysical neutrinos, opening a whole new window to the Universe. In these measurements, it is important to provide careful estimates of all relevant backgrounds, which in this case are dominated by the prompt flux of neutrinos from charm decays produced in cosmic rays collisions in the atmosphere. The key ingredient of this calculation was the exploitation of the constraints from the LHCb measurements of forward charm production to reduce the PDF uncertainties of the poorly-know gluon at very small-x. Our results, together with recent IceCube bounds, indicate that a first direct measurement of this prompt flux is now within reach, which would allow a final vindication of the cosmological nature of the highest energy neutrinos ever measured.

The figure below shows the comparison between our updated predictions for the prompt neutrino flux as a function of the neutrino energy, with theoretical uncertainties (blue band) compared with a previous benchmark calculation (ERS) and with a direct exclusion bound from the analysis of the IceCube data (dotted red line).

The prompt neutrino flux at IceCube



Postdoctoral research positions

The theory group of Nikhef is advertising the opening of two postdoctoral research positions, one of which within our group. The deadline for these applications is the 1st of December 2017. Full details about these positions can be found hereApplications for these positions should be submitted via

Moreover, we welcome applications for individual postdoctoral fellowships including:

PhD student positions

There are currently no openings for PhD positions in the group.


UvA/VU Master in Physics and Astronomy

Research projects for 2nd year students of the joint UvA/VU Master programme in Physics and Astronomy are available. These projects are suitable both for those students following the Theory track as well as the GRAPPA track. Illustrative examples of some of the available projects can be found here:

  • A three-dimensional imaging of the proton with Neural Networks (description)
  • Constraining the proton structure with Run II LHC data (description)
  • Probing electroweak symmetry breaking with Higgs pair production at the LHC and beyond (description)

Interested students should contact me at This email address is being protected from spambots. You need JavaScript enabled to view it.  by April or May at the latest, to arrange a meeting to discuss the details of the project.

 Junior research positions

Candidates interested in applying to five-year junior research positions with Nikhef/VU as host institution are encouraged to contact any of the staff members with enough time in advance. These include


One of the main activities of PDF4BSM is to establish a strong visitor program in Oxford to develop synergies between researchers of different LHC communities with a common interest in parton distribution functions. A list of Oxford visitors supported by PDF4BSM can be found below:



  • Carlota Andres Casas (Santiago de Compostela), April 2015 - June 2015.

Carlota did her undergraduate studies at the Universidade de Santiago de Compostela, followed by a Master degree under the supervision of Carlos Salgado. Carlota's research interests focus on initial stage effects and jet quenching in high energy collisions at the LHC and they include nuclear Parton Distributions and medium induced gluon radiation.


  • Simone Marzani (Massachusets Institute of Technology - MIT)

Simone did his undergraduate studies at the University of Milan, followed by a doctoral studies at the University of Edinburgh, where he obtained his PhD in 2008. He then went on for postdoc positions at the University of Manchester (2008-2011) and Durham University (2011-2014) and he currently holds an LHC-TI fellowship at MIT. His research focus his in LHC phenomenology with particular focus on applying all-order techniques in QCD to Higgs production, jet substructure and small-x physics. 


  • Richard Ruiz (Pittsburgh), September 2014

Richard E. Ruiz is a K. Leroy Irvis Graduate Fellow and PhD Candidate in the Department of Physics & Astronomy at the University of Pittsburgh. His primary research addresses new tests of TeV-scale neutrino mass-generating mechanisms (Seesaw models), extended gauge sectors (LR-Symmetric, GUT models), and extended scalar sectors (SU(2)L doublets and triplets), with an emphasis on hadron collider signatures involving final-state scalar, vector bosons, and top quarks.