Я работаю над воспроизведением рисунков 1 и 2 из статьи.
Эти рисунки:

и

Для Рис.2 только основной график, а не вставки.
В документе говорится следующее:
"Большая часть наших данных отображается для середины энергетического спектра, где начинается переход. Модель имеет два класса симметрии, характеризующиеся нечетным или четным числом кубитов, и данные приведены для одного класса симметрии. Чтобы уменьшить статистические флуктуации, мы используем $5 ≤ ND ≤ 4 × 10^4$ случайных реализаций Γi и J_{ij} , как это обычно делается в теории случайных матриц. Собственные значения и собственные векторы вычисляются путем точной диагонализации матрицы Гамильтона для каждой реализации. общее количество интервалов составляет 10^4 < NS ≤ 1,6 × 10^5 (NS ∝ ND NH ). Пример перехода в спектральной статистике показан на рис.1. Для анализа эволюции P(s) при связи J удобно использовать параметр η = [∫₀^{s₀} P(s)ds - I_W] / [I_P - I_W] где I_P = ∫₀^{s₀} P_P(s)ds, I_W = ∫₀^{s₀} P_W(s)ds и s₀=0,4729... является точкой пересечения P_{P}(s) и P_{W}(s). Таким образом, P_{P}(s) соответствует η=1, а P_{W}(s) Как обычно в области квантового хаоса, изменение η характеризует эволюцию P(s). Изменение η относительно J/∆0 представлено на рис. 2 для δ = 0, показывая, что действительно η падает от 1 до 0 с увеличением силы связи. условие η(J_{c})=0,3. Зависимость J_{c} от n представлена на рис.2"
Упомянутый гамильтониан имеет вид:
$H = \sum_i Gamma_i\sigma_i + \sum_{i
Подробнее здесь: https://stackoverflow.com/questions/798 ... estimation
Мобильная версия