Use $ams2026

Run non-optimizing PES scans where the C coordinate is scaled (LatticeCRange)
from 6.0 to 9.0 in steps of 0.1 angstrom, and then another run where it is
continued from 9.0 to 16.0 in steps of 0.5 angstrom.

Initial graphite system:

```
System
   Atoms
      C 0 0 1.5
      C 0 0 4.5
      C 0 1.4202816622064793 1.5
      C 1.23 0.7101408311032397 4.5
   End
   Lattice
      2.46 0 0
      -1.23 2.130422493309719 0
      0 0 6
   End
End
```

Use QuantumESPRESSO with 12x12x4 k-points, gaussian smearing 0.001 Ry, energy
cutoff 100 Ry, density cutoff 1000 Ry.

Set `input_dft` explicitly.

Functionals:

- PBE with pslibrary-PAW
- BP86bPBE with pslibrary-PAW and XDM dispersion correction
- r2scan with Dojo
- r2scan with Dojo and D4 Engine addon dispersion correction.

In the report plot a figure with relative energy in meV/atom vs c/2 in angstrom
(c/2 = half c length, this is the interlayer graphite distance), where the
relative energy is energy relative to the energy of the longest c parameter. The
figure should have the plot for all tested functionals, and plot lines and small
points.

Plot the curves in the reverse order of minimum energy so that the order of the
colored lines in the legend matches the order in which they are shown in the
plot.