Оптимизировать код для работы над квантовым хаосомPython

Программы на Python
Ответить
Anonymous
 Оптимизировать код для работы над квантовым хаосом

Сообщение Anonymous »

Я работаю над воспроизведением рисунков 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) - η=0. Как это обычно бывает в области квантового хаоса, Изменение η характеризует эволюцию P(s). Изменение η относительно J/∆0 представлено на рис. Зависимость J_{c} от n приведена на рис.2"

Упомянутый гамильтониан имеет вид:
$H = \sum_i Gamma_i\sigma_i + \sum_{i

Подробнее здесь: https://stackoverflow.com/questions/798 ... s-disorder
Ответить

Быстрый ответ

Изменение регистра текста: 
Смайлики
:) :( :oops: :roll: :wink: :muza: :clever: :sorry: :angel: :read: *x)
Ещё смайлики…
   
К этому ответу прикреплено по крайней мере одно вложение.

Если вы не хотите добавлять вложения, оставьте поля пустыми.

Максимально разрешённый размер вложения: 15 МБ.

Вернуться в «Python»