Analiza Procesu Obliczeniowego Modelu Phinix-1 (Test 1)
Niniejszy dokument zawiera rygorystyczny opis procesu obliczeniowego dla energii emisyjnej ciała doskonale czarnego w ujęciu mereologicznym (bez punktów Cantora) z wprowadzeniem sztywnych geometrycznych granic sumowania.
1. Pełen Opis Procesu Obliczeniowego
W modelu Phinix-1 całkowita gęstość energii we wnęce o rozmiarze liniowym \(L\) nie jest całką rozciągającą się od częstotliwości \(0\) do nieskończoności \(\infty\). Sumowanie stanów jest zablokowane geometrycznie.
Krok 1: Wyznaczenie granic operacyjnych
Dolna granica częstotliwości \(\nu_{\text{min}}\) wynika ze skwantowania najdłuższego dopuszczalnego modu fali (harmonicznej podstawowej), który fizycznie może stać się częścią wnęki:
Górna granica częstotliwości \(\nu_{\text{max}}\) jest zdeterminowana przez dno strukturalne czasoprzestrzeni – atom mereologiczny \(\delta\):
Krok 2: Całkowanie numeryczne (Skończona suma)
W temperaturze pokojowej (\(T = 300\text{ K}\)) gęstość energii dla wysokich częstotliwości wygasza się samoczynnie w sposób wykładniczy na długo przed osiągnięciem granicy \(\nu_{\text{max}}\). Dlatego rzeczywistą, nieliniową modyfikację rozkładu opisuje ucięcie dolnej granicy pasma \(\nu_{\text{min}}\).
Obliczamy gęstość energii mereologicznej \(E_{\text{mereo}}\) za pomocą numerycznej metody trapezów:
Gdzie geometryczny współczynnik skali \(\alpha\) odpowiada stałej Plancka \(h\).
Krok 3: Wyznaczenie punktu odniesienia Cantora (\(E_{\text{classic}}\))
Klasyczny model ciągły zakłada, że wnęka może pomieścić nieskończenie długie fale (\(\nu_{\text{min}} = 0\)). Całkowita gęstość energii zintegrowana w pełnym kontinuum wynosi:
2. Wyniki Cząstkowe i Tabela Porównawcza
Poniższe obliczenia numeryczne wykonano dla stałych zweryfikowanych laboratoryjnie: * \(\alpha = 6.62607015 \times 10^{-34}\text{ J}\cdot\text{s}\) * \(c = 299792458\text{ m/s}\) * \(k = 1.380649 \times 10^{-23}\text{ J/K}\) * \(T = 300\text{ K}\) (Temperatura pokojowa)
Dla tych warunków klasyczna gęstość energii wynosi: \(E_{\text{classic}} = 6.1282 \times 10^{-6}\text{ J/m}^3\).
Tabela: Wpływ skrótu geometrycznego na energię układu (T = 300 K)
| Rozmiar \(L\) (m) | Rozmiar \(L\) (nm / µm) | Odcięcie \(\nu_{\text{min}}\) (Hz) | Granica \(\lambda_{\text{max}}\) | \(E_{\text{mereo}}\) (\(\text{J/m}^3\)) | Stosunek \(\frac{E_{\text{mereo}}}{E_{\text{classic}}}\) | Status fizyczny układu |
|---|---|---|---|---|---|---|
| \(1.00 \times 10^{-4}\) | \(100\text{ \µm}\) | \(1.49 \times 10^{12}\) | \(200000\text{ nm}\) | \(6.1242 \times 10^{-6}\) | 0.9993 | Zachowanie klasyczne (Skala Makro) |
| \(2.66 \times 10^{-5}\) | \(26.6\text{ \µm}\) | \(5.63 \times 10^{12}\) | \(53200\text{ nm}\) | \(5.9657 \times 10^{-6}\) | 0.9735 | Początek blokady długofalowej |
| \(6.58 \times 10^{-6}\) | \(6.58\text{ \µm}\) | \(2.27 \times 10^{13}\) | \(13160\text{ nm}\) | \(2.8880 \times 10^{-6}\) | 0.4712 | Silna anomalia emisyjna |
| \(1.63 \times 10^{-6}\) | \(1.63\text{ \µm}\) | \(9.19 \times 10^{13}\) | \(3260\text{ nm}\) | \(1.5103 \times 10^{-9}\) | 0.0002 | Krytyczne wygaszenie promieniowania |
| \(4.04 \times 10^{-7}\) | \(404\text{ nm}\) | \(3.71 \times 10^{14}\) | \(808\text{ nm}\) | \(3.3226 \times 10^{-27}\) | \(5.42 \times 10^{-22}\) | Pełna blokada (Stan \(\frac{1}{2}\) Łukasiewicza) |
| \(1.00 \times 10^{-7}\) | \(100\text{ nm}\) | \(1.49 \times 10^{15}\) | \(200\text{ nm}\) | \(0.0000\) | 0.0000 | Zapadnięcie continuum (Skala Atomowa) |
3. Kod Skryptu Wizualizacyjnego (Matplotlib na Kubuntu)
Poniższy skrypt napisany w czystym Pythonie 3 wygeneruje i zapisze na Twoim dysku plik wykresu porównawczego w wysokiej rozdzielczości (300 DPI). Kod jest w pełni otwarty (FOSS) i zoptymalizowany pod standardowe środowisko Linux.
#!/usr/bin/env python3
import os
import numpy as np
import matplotlib.pyplot as plt
# 1. Definicja stałych fizycznych
alpha = 6.62607015e-34 # Geometryczny współczynnik skali (h)
c = 299792458 # Prędkość światła (m/s)
k = 1.380649e-23 # Stała Boltzmanna (J/K)
T = 300 # Temperatura pokojowa (K)
# 2. Obliczenie referencyjnej energii klasycznej Cantora
E_classic = (8 * np.pi**5 * (k * T)**4) / (15 * c**3 * alpha**3)
# 3. Generowanie dziedziny rozmiarów wnęki L (od 100 nm do 100 um)
L_plot = np.logspace(-7, -4, 200)
ratios = []
# 4. Pętla obliczeniowa dla modelu mereologicznego
for L in L_plot:
nu_min = c / (2 * L)
# Bezpiecznik numeryczny dla skrajnie małych wartości (zapobieganie overflow w exp)
if (alpha * nu_min) / (k * T) > 100:
ratios.append(0.0)
else:
# Generowanie siatki częstotliwości od nu_min do bezpiecznego odcięcia termicznego
freqs = np.linspace(nu_min, 1e15, 10000)
integrand = (8 * np.pi * alpha * freqs**3) / (c**3 * (np.exp((alpha * freqs) / (k * T)) - 1))
E_mereo = np.trapezoid(integrand, freqs)
ratios.append(float(E_mereo / E_classic))
# 5. Tworzenie wykresu za pomocą matplotlib
plt.style.use('seaborn-v0_8-whitegrid' if 'seaborn-v0_8-whitegrid' in plt.style.available else 'default')
fig, ax = plt.subplots(figsize=(10, 6))
ax.semilogx(L_plot * 1e6, ratios, label='Model Mereologiczny Phinix-1', color='firebrick', lw=2.5)
ax.axhline(1.0, color='darkslategrey', linestyle='--', alpha=0.7, label='Model Ciągły Cantora (Limit Plancka)')
# Stylizacja osi i siatki (Kubuntu Desktop Ready)
ax.set_title('Złamanie Limitu Plancka (Blokada Geometryczna) w T = 300K', fontsize=13, pad=15, fontweight='bold')
ax.set_xlabel('Rozmiar liniowy wnęki L [µm]', fontsize=11)
ax.set_ylabel('Stosunek całkowitej energii emisyjnej (E_mereo / E_classic)', fontsize=11)
ax.set_xlim(0.1, 100)
ax.set_ylim(-0.05, 1.05)
ax.grid(True, which="both", ls=":", alpha=0.6)
ax.legend(loc='lower right', frameon=True, facecolor='white', framealpha=0.9, fontsize=10)
# 6. Zapis do katalogu mkdocs (dostosuj ścieżkę jeśli docs/img istnieje)
output_filename = 'blokada_geometryczna_phinix1.png'
plt.tight_layout()
plt.savefig(output_filename, dpi=300)
print(f"[OK] Wykres zapisany pomyślnie jako: {output_filename}")