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 generatethe 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 textfile (figN_*_data.txt) and, where applicable, a NumPy archive(figN_*_solutions.npz), so that the results can be reproduced withoutre-running the scripts. All scripts require Python 3 with numpy, scipy, and matplotlib.