Quantitative holographic correspondence between the Peierls instability and the Yukawa-Sachdev-Ye-Kitaev model

This repository contains the Python scripts and raw data used to generate
the results discussed in the manuscript:

    "Quantitative holographic correspondence between the Peierls
     instability and the Yukawa-Sachdev-Ye-Kitaev model"

The four scripts correspond to the four figures of the manuscript:

    - bosonization.py
        Computes the Luttinger parameter K of the Peierls-Hubbard model
        as a function of U/t and lambda/t, and generates the
        K <-> M/N mapping. Produces Figure 1.

    - yukawa_syk.py
        Solves the Yukawa-SYK saddle-point equations and computes the
        boson scaling dimension Delta_phi(M/N) and the critical
        coupling g_c(M/N). Produces Figure 2.

    - susceptibility.py
        Computes the bond-order susceptibility of the Peierls chain
        and the order-parameter susceptibility of the Yukawa-SYK model,
        showing they share the same scaling exponent 2K - 1.
        Produces Figure 3.

    - conformal_higgs.py
        Solves the one-loop fixed-point equations for the minimal
        conformally extended Higgs sector, with electroweak gauge
        contributions. Produces Figure 4.

The raw data generated by each script is saved to a companion text
file (figN_*_data.txt) and, where applicable, a NumPy archive
(figN_*_solutions.npz), so that the results can be reproduced without
re-running the scripts.

All scripts require Python 3 with numpy, scipy, and matplotlib.