Final data for paper "A nonperturbative test of nucleation calculations for strong phase transitions"

In this Zenodo deposit we include the data used to make three key figures in our paper, showing our final nucleation rate results. The columns of the raw data files are labelled. In each case, the quantity being tabulated is the logarithm of the nucleation rate, typically labelled logGamma in our data files. We also include plotting scripts to generate the figures from the data.

• The directory a_limit contains the extrapolation to zero lattice spacing, with the raw data in the file a_limit_data. The temperature T is fixed to our benchmark point of 93.121 GeV. The linear size Nx and the lattice spacing a are varied to give a constant total volume. The file a_limit_lin_op_data contains our comparison measurement using the linear order parameter. The script plot_a_limit.py performs a nonlinear least squares fit to a cubic, writes the continuum extrapolated nuclation rate result to stdout, and plots the data and the fit.

• The directory vol_limit contains the extrapolation to infinite volume, with the raw data in the file vol_limit_data. Here the temperature T=93.121 GeV and lattice spacing a=1.5 are fixed, while the linear size Nx is varied. The script plot_vol_limit.py performs a nonlinear least squares fit to an exponential function, writes the infinite volume limit nucleation rate result to stdout, and again plots the data and the fit.

• The directory rate_reweighted contains the data needed to plot the nucleation rate as a function of temperature. This plot combines lattice data and various perturbative estimates of the nucleation rate. The lattice data are in three files:

  • rate_lattice_BM2_vol_limit contains the infinite volume limit result at the simulated lattice temperature T=93.121 GeV and with lattice spacing a=1.5.
  • rate_lattice_BM2 reweights the lattice data to different temperatures and then takes the infinite volume limit.
  • rate_lattice_BM2_Nx40_a1.5 contains data reweighted from the simulation at fixed volume Nx=40. In each lattice data file, the temperature, continuum potential parameters and nucleation rate are given.

  The perturbative results are in three files, rate_perturbative_BM2_rg_low, rate_perturbative_BM2 and rate_perturbative_BM2_rg_high, corresponding to the three RG scales mentioned in the text. These files again contain the temperature and continuum lattice parameters. The column eps gives the dimensionless loop expansion parameter around the metastable phase within the EFT. The column logGamma_0 gives the tree level result, logGamma_A and logGamma_B the LPA results with the two options for dealing with the imaginary parts described in the text, and logGamma_1 is the one loop result. The script plot_rate_reweighted.py plots all of these data together. For the tree level and one loop cases, the error bands plotted correspond to the minimum and maximum values of the nucleation rate across the three RG scales. For the LPA results, the bands include the extreme values for all six options including the two approaches to handling the imaginary parts.

In each case, the plots are saved to PDF and PNG plot files.

The figures used in this deposit and in the paper used the following package versions (obtained with the pipreqs package):

matplotlib==3.5.1
numpy==1.21.5
scipy==1.8.0
seaborn==0.11.2

Note that in Version v1 there was an error in the normalisation of our perturbative tree-level and LPA results which has now been fixed.