Проверка 4D-модели квантовой памяти: сходится ли математика по уравнениям Аллена-Кана?

Question

[Александр Петров]

Привет, инженеры, программисты и физики! Я слесарь КИПиА. На смене пришла идея: спроектировал контур 4D-памяти в капле сверхчистой воды через ТГц-лазер (2.4 ТГц), уставку угла 104,5° и дифракцию кварков в возбужденном вакууме. Набросал простейший код симуляции на Python. Подскажите, корректно ли работает алгоритм компенсации теплового шума или я где-то ошибся в калибровочном контуре? Кто может, протестируйте софт у себя на компьютерах.

import numpy as np
import plotly.graph_objects as go

# ==========================================
# КОНТУР АВТОМАТИЗАЦИИ (КИПиА) КХД-ЯЧЕЙКИ
# ==========================================
N_GRID = 30                 # Матрица пространства-времени 30х30х30
VALENCE_SETPOINT = 104.5    # Уставка: валентный угол H-O-H (градусы)
LASER_FREQ_THZ = 2.4       # Частота исполнительного механизма (лазера)
THERMAL_DEVIATION = 1.8     # Возмущающее воздействие (тепловой шум среды)

# Инициализируем 3D-координаты пространства капли
x = np.linspace(-1, 1, N_GRID)
y = np.linspace(-1, 1, N_GRID)
z = np.linspace(0, 1, N_GRID)
X, Y, Z = np.meshgrid(x, y, z)

# Расчет геометрии: берем только пограничную EZ-зону у подложки (Z близко к 0)
ez_layer_mask = (X**2 + Y**2 <= 0.8) & (Z >= 0) & (Z <= 0.4)

X_cell = X[ez_layer_mask]
Y_cell = Y[ez_layer_mask]
Z_cell = Z[ez_layer_mask]

# ==========================================
# ФИЗИЧЕСКИЙ РАСЧЕТ ДИФРАКЦИИ И ВОЗБУЖДЕНИЯ ВАКУУМА
# ==========================================
# Симулируем интерференционную решетку возбужденного лазером вакуума
vacuum_excitation = (np.sin(5 * X_cell) * np.cos(5 * Y_cell)) ** 2

# Кварковая дифракция заполняет канавки трафарета вакуума
# Выпрямляем знак случайной погрешности через np.abs
np.random.seed(1045)
vibration_noise = np.random.normal(0, THERMAL_DEVIATION, size=Z_cell.shape)

# Расчет отклонения от уставки. Чем сильнее возбужден вакуум, тем меньше ошибка!
angle_error = np.abs(vibration_noise * np.exp(-3 * vacuum_excitation))
current_angles = VALENCE_SETPOINT - angle_error

# Фильтруем данные: оставляем точки, где структура удержала уставку
stable_bits_mask = angle_error < 5.0
X_fin = X_cell[stable_bits_mask]
Y_fin = Y_cell[stable_bits_mask]
Z_fin = Z_cell[stable_bits_mask]
Color_fin = current_angles[stable_bits_mask]

# ==========================================
# ВИЗУАЛИЗАЦИЯ ДЛЯ ПРЕЗЕНТАЦИИ (СТИЛЬ СТАРКА)
# ==========================================
fig = go.Figure(data=go.Scatter3d(
    x=X_fin, y=Y_fin, z=Z_fin,
    mode='markers',
    marker=dict(
        size=3,
        color=Color_fin,
        colorscale='Hot',  # Светлые точки - идеальное удержание уставки 104.5°
        opacity=0.8,
        colorbar=dict(title="Текущий угол связи (Уставка = 104.5°)")
    )
))

fig.update_layout(
    title=f"Прототип 4D-памяти: Дифракция кварков в возбужденном КХД-вакууме капли<br>Частота лазера: {LASER_FREQ_THZ} ТГц",
    scene=dict(
        xaxis_title='X (Нанометры)', 
        yaxis_title='Y (Нанометры)', 
        zaxis_title='Z (Высота EZ-слоя)',
        backgroundcolor="black"
    ),
    template="plotly_dark"
)

# Сохраняем готовый 3D-интерфейс
fig.write_html("quantum_drop_prototype.html")
print("ПРОТОТИП СГЕНЕРИРОВАН! Открывай файл 'quantum_drop_prototype.html' в браузере.")