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

<p>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&nbsp;<em>logGamma</em>&nbsp;in our data files. We also include plotting scripts to generate the figures from the data.</p> <ul> <li> <p>The directory&nbsp;<code>a_limit</code>&nbsp;contains the extrapolation to zero lattice spacing, with the raw data in the file&nbsp;<code>a_limit_data</code>. The temperature&nbsp;<em>T</em>&nbsp;is fixed to our benchmark point of 93.121 GeV. The linear size&nbsp;<em>Nx</em>&nbsp;and the lattice spacing&nbsp;<em>a</em>&nbsp;are varied to give a constant total volume. The file&nbsp;<code>a_limit_lin_op_data</code>&nbsp;contains our comparison measurement using the linear order parameter. The script&nbsp;<code>plot_a_limit.py</code>&nbsp;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.</p> </li> <li> <p>The directory&nbsp;<code>vol_limit</code>&nbsp;contains the extrapolation to infinite volume, with the raw data in the file&nbsp;<code>vol_limit_data</code>. Here the temperature&nbsp;<em>T</em>=93.121 GeV and lattice spacing&nbsp;<em>a</em>=1.5 are fixed, while the linear size&nbsp;<em>Nx</em>&nbsp;is varied. The script&nbsp;<code>plot_vol_limit.py</code>&nbsp;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.</p> </li> <li> <p>The directory&nbsp;<code>rate_reweighted</code>&nbsp;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:</p> <ul> <li><code>rate_lattice_BM2_vol_limit</code>&nbsp;contains the infinite volume limit result at the simulated lattice temperature&nbsp;<em>T</em>=93.121 GeV and with lattice spacing&nbsp;<em>a</em>=1.5.</li> <li><code>rate_lattice_BM2</code>&nbsp;reweights the lattice data to different temperatures and then takes the infinite volume limit.</li> <li><code>rate_lattice_BM2_Nx40_a1.5</code>&nbsp;contains data reweighted from the simulation at fixed volume&nbsp;<em>Nx</em>=40. In each lattice data file, the temperature, continuum potential parameters and nucleation rate are given.</li> </ul> <p>The perturbative results are in three files,&nbsp;<code>rate_perturbative_BM2_rg_low</code>,&nbsp;<code>rate_perturbative_BM2</code>&nbsp;and&nbsp;<code>rate_perturbative_BM2_rg_high</code>, corresponding to the three RG scales mentioned in the text. These files again contain the temperature and continuum lattice parameters. The column&nbsp;<em>eps</em>&nbsp;gives the dimensionless loop expansion parameter around the metastable phase within the EFT. The column&nbsp;<em>logGamma_0</em>&nbsp;gives the tree level result,&nbsp;<em>logGamma_A</em>&nbsp;and&nbsp;<em>logGamma_B</em>&nbsp;the LPA results with the two options for dealing with the imaginary parts described in the text, and&nbsp;<em>logGamma_1</em>&nbsp;is the one loop result. The script&nbsp;<code>plot_rate_reweighted.py</code>&nbsp;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.</p> </li> </ul> <p>In each case, the plots are saved to PDF and PNG plot files.</p> <p>The figures used in this deposit and in the paper used the following package versions (obtained with the&nbsp;<code>pipreqs</code>&nbsp;package):</p> <pre><code>matplotlib==3.5.1 numpy==1.21.5 scipy==1.8.0 seaborn==0.11.2</code></pre> <p>Note that in <a href="https://doi.org/10.5281/zenodo.10891524">Version v1</a> there was an error in the normalisation of our perturbative tree-level and LPA results which has now been fixed.</p>