Numerical simulation of quark-gluon plasma formed in the collision of two gold atoms.

Calculations used to produce this figure are from this work.

Research

Quarks and gluons are fundamental building blocks of matter. Their interactions are governed by Quantum Chromodynamics, one of the core parts of the Standard Model of particle physics.

Microseconds after the Big Bang, a plasma of quarks and gluons filled the Universe. As the Universe expanded, its temperature decreased until quarks and gluons could combine into protons, neutrons, nuclei and other bound states.

Aerial view of the Relativistic Heavy Ion Collider at Brookhaven National Laboratory

Quark-gluon plasma can now be produced in laboratory, by colliding large nuclei at velocities close to the speed of light. The two colliders that can achieve this feat are the Relativistic Heavy Ion Collider at Brookhaven National Laboratory (New York) and the Large Hadron Collider at CERN (Switzerland/France).



Relativistic collisions of two lead nuclei as seen by the CMS detector circa 2011

Quark-gluon plasma is never seen directly: its existence and properties must be inferred from the shower of subatomic particles that hits the detector (image on the right). I study the properties of quark-gluon plasma using multi-stage numerical simulations whose core is relativistic viscous fluid dynamics.

I also use Bayesian inference to quantify the properties of quark-gluon plasma with reliable uncertainties. The specific shear viscosity of quark-gluon plasma is among the smallest ever measured in a liquid (an honour shared with cold atomic gases); it is orders of magnitude smaller than that of liquid helium, for example.

Temperature dependence of the specific shear viscosity

The figure on the left shows a recent constraint on the specific shear viscosity η/s of quark-gluon plasma as a function of the plasma's temperature (the orange band is a Bayesian model average, the bottom panel is the posterior-prior Kullback–Leibler divergence). This result is the product of a multi-institution collaboration, obtained from millions of core-hours of calculations on US supercomputers.

While the viscosity of strongly-coupled quark-gluon plasma is famously difficult to calculate, the long-term objective is to find agreement between measurements of the viscosity and ab initio calculations.

I also study Apparent temperature of quark-gluon plasma as seen through the photon energy spectrum, as a function of the true maximum temperature of the plasma and of the energy range of the measured photons the production of photons by the quark-gluon plasma. Because the quark-gluon plasma produced in nuclear collisions is small, photons can escape with negligible rescattering, and they can be used to probe deep into the hottest and densest parts of the plasma. Recently, I have used analytical solutions of relativistic fluid dynamics (the "Gubser solution") to better understand the relation between the energy spectrum of photons, the temperature of the plasma and the Doppler shift on photons resulting from the plasma's expansion, finding \(T_{\textrm{eff}} \approx T_0/(1+\frac{5}{2} \frac{T_0}{E})\). The figure on the left shows this relation between the actual maximum temperature of the plasma and the apparent ("effective") temperature extracted from the photon spectrum in different energy ranges.

See also our latest work on the relation between the speed of sound of quark-gluon plasma and the final particles' energy and multiplicity.



A high-level overview of ultrarelativistic heavy-ion collisions

(adapted from a topical review I wrote in 2024)

            We aim to describe nucleus collisions from the initial impact of the nuclei to the final shower of particles observed by detectors. The problem is broadly approached piecewise. Nuclear structure properties of heavy nuclei, measured in scattering experiments, can be combined with effective field theory calculations to understand the state of the nuclei shortly after their impact. Many-body quantum field theory is used to study the interaction of energetic quarks and gluons with their lower-energy counterparts. In certain regimes, the collective motion of strongly-interacting quarks and gluons is described with relativistic fluid dynamics, using an equation of state calculated numerically from lattice quantum field theory. By combining these different theoretical approaches and more, one can obtain a multi-stage model that is a functional description of the successive stages of ultrarelativistic nuclear collisions. Any quantity that is unknown or only partly constrained by theoretical guidance can be considered a model parameter that enters the description of the collision. For example, the shear and bulk viscosities of strongly-coupled nuclear matter are generally parametrized as a function of temperature, and the same can be done for the equation of state.

A large amount of data is available to constrain the model's parameters, a result of the comprehensive heavy-ion-collision program at the Relativistic Heavy Ion Collider and the Large Hadron Collider. The detectors measure the number and momentum distributions of different species of particles produced after each collision, as well as various correlations between these particles. The final shower of particles fluctuates from one collision to the next, due to the variable degree of overlap of nuclei --- which cannot be controlled --- and the quantum nature of nuclear collisions. Measurements are defined by averaging over a large number of collisions. The variable degree of overlap of the nuclei --- the centrality --- provides a natural parameter that changes the system's size, geometry and density. By further colliding different nuclei species and varying the collision energy of the nuclei, one can probe the system in a large number of different configurations.

The number of measurements available is large, and the type of measurements varies from simple averages (e.g., particle numbers) to intricate correlation observables. The measurements are stochastic, with statistical uncertainties ranging from a fraction of a percent to tens of percent. There are also systematic uncertainties that can be correlated across different data sets, although information about these correlations is often not readily available. On the theory side, multistage models provide a general description of the collisions to compare with these measurements. There is uncertainty in the physics description of the collisions, represented in the literature by different models used to describe the same stage of the collision. For example, a given stage of the collision can be parametrized with a flexible ansatz, or one can use an ab initio model with fewer parameters but additional physics assumptions.

JFP, "Applications of emulation and Bayesian methods in heavy-ion physics" (2024) J.Phys.G.

Group

J.-F. Paquet

Gabriel Soares Rocha
(postdoc)

Teerthal Patel
(postdoc)

Mayank Singh
(postdoc)

Andi Mankolli
(PhD student)

McKenna Sleeth
(PhD student)

Harvir Dhindsa
(PhD candidate)

Joseph Sanders
(PhD candidate)

Jeseleth Van Rose
(PhD candidate)

Collaborations

I am a member of the Vanderbilt Initiative for Gravity, Waves, and Fluids (VandyGRAF), which brings together 20+ faculty, postdocs and graduate students from physics, astrophysics and math for weekly seminars and discussions. We have been hosting a stellar lineup of visitors across disciplines, which you can see here.

I am also a member of the JETSCAPE Collaboration. Since 2018, I have been co-convener of the Simulation and Distribution Computing Working Group, one of the four working groups of the collaboration, which focuses on the soft sector of heavy-ion collisions, and in particular on Bayesian studies using low-energy hadron measurements from the Relativistic Heavy Ion Collider and the Large Hadron Collider.

Undergraduate research opportunities

Reach out if you are interested in an undergraduate research project: jean-francois.paquet@vanderbilt.edu


Projects of interest include
  • Relativistic fluid dynamics
  • Numerical simulations of relativistic nuclear collisions and data analysis
  • Perturbative field theory
  • Bayesian inference, emulation, machine learning


Sample of previous undergraduate research projects:

Teaching

Fall 2022, 2023 & 2024

PHYS 8010: Particle and Continuum Mechanics


I use Goldstein's classic, and Hand and Finch's Analytical Mechanics, as references. [Note: Hand and Finch does appear to have some typos].

There are excellent lecture notes available online, including David Tong's.

Spring 2024

PHYS 8040: Statistical Mechanics


As an undergrad, I did a term project on the life and work of Josiah Willard Gibbs: what a fascinating story. He studied at Yale and got his PhD in engineering (apparently the fifth PhD ever granted by a U.S. institution).                   He worked as a tutor for a few years, then visited Europe and attended lectures in France and Germany by some of the leading physicists of the time.                   Sometime after his return to New Haven, he secured a faculty position at Yale, and then proceeded to revolutionize statistical mechanics...
                 
                  It may be of some historical interest to recall the peculiar circumstances connected with the publication of the famous memoirs of Professor Gibbs, on which his great reputation rests. They were published in the Transactions of the Connecticut Academy of Sciences, which at that time was not connected with Yale University and received no support from it. When Professor Gibbs presented his papers, the academy had no regular publication funds [...] They were expensive to set up, owing to the complex mathematical formulae. Our funds were small. On nearly every occasion, we had to go out and raise a subscription to pay the cost, partly among college men and partly among the business and professional men of New Haven. Long discussions took place as to our ability to print the articles. Two able mathematical professors were on the committee Loomis and Newton. Both protested that they did not understand Gibbs's papers at all. One insisted that no man ever lived who could except Maxwell and he was dead. Yet we all believed that what Gibbs wrote must be of intrinsic value in his branch of science. Therefore we raised the money and printed each paper as it came in.
A.E. Verrill, "How the works of Professor Willard Gibbs were published" (1925) Science

Spring 2022 (Duke) & 2023

Introduction to Fluid Dynamics (for nuclear collisions)



Too few people have had the chance to watch this series of videos from 1961: https://web.mit.edu/hml/ncfmf.html
See for example, Sir G.I. Taylor himself explain [non-relativistic] fluid dynamics at low Reynolds numbers.
These videos are a reminder of the decades of experimentation it took to establish the foundations of non-relativistic fluid dynamics.

Selected recent presentations

Nuclear collisions as seen through photons
Nuclear physics seminar (Nov 2023, Institute for Nuclear Theory, University of Washington)

Flow and transport properties: Assessment and outlook
Invited talk at the 2022 Hot and Cold QCD Town Hall meeting (Sept 2022, MIT)

Matching conditions across time evolution stages of ultrarelativistic heavy-ion collisions
Plenary talk prepared for Initial Stages 2021

Multi-system Bayesian constraints on the transport coefficients of QCD
Parallel talk prepared for Quark Matter 2019 (Wuhan), on behalf of the JETSCAPE Collaboration

Publications

A complete list of my publications is available on INSPIRE [ INSPIRE profile ], with links to the published versions and the (free) arXiv versions. Google Scholar is also eerily good at finding my work.



In the summer of 2018, I wrote a manual for the heavy-ion hydrodynamic simulation code MUSIC as support material for lectures I gave at a workshop at the University of São Paulo. Since then, a new official version of MUSIC has been published, and I updated the manual accordingly for this new version: here. The older version of MUSIC is still usable, of course; its manual remains available here. While it is MUSIC-oriented, the manual should provide a good overview of the physics behind hydrodynamic simulations of heavy-ion collisions in general.

Theses

Master's: Photon and neutral pion production in d+Au collisions at RHIC
PhD: Characterizing the non-equilibrium quark-gluon plasma with photons and hadrons

Scientific software

MUSIC

We use relativistic viscous fluid dynamics to describe the spacetime expansion of quark-gluon plasma. MUSIC is a full 3+1D code developed in McGill that has been widely used by us and others to study heavy-ion collisions.

The official website is http://www.physics.mcgill.ca/music/.

The latest version of the code is available here: https://github.com/MUSIC-fluid/MUSIC/

A legacy version of the code is hosted on Sourceforce: https://sourceforge.net/projects/music-hydro/.

In 2018, I wrote a manual for the code, which I have updated since to keep up with changes in the code: https://webhome.phy.duke.edu/~jp401/music_manual/


KøMPøST

KøMPøST: a kinetic theory propagator for the early stage of heavy-ion collisions

Other codes: My Github

https://github.com/j-f-paquet

A few highlights:

About

Jean-François Paquet
[ jean-francois.paquet@vanderbilt.edu ]

Assistant Professor of Physics and Mathematics
Vanderbilt University
Nashville, TN


Previous position:


Postdoctoral Researcher and Research Scientist
Department of Physics
Duke University
Durham, NC

Postdoctoral Researcher
Department of Physics & Astronomy
Stony Brook University
Stony Brook, NY


Education:


Master's and PhD in Physics
Under the supervision of Prof. Charles Gale
McGill University, Montréal

Bachelor's degree in Physics
Université de Sherbrooke, Sherbrooke