ANEKS 2 - Dynamiczny licznik relacyjny i hipoteza kontaminacji tlenowej Cz-Si

Abstrakt [PL]:

W niniejszym aneksie do Zarzewia 1 usuwamy ostatnią stałą fenomenologiczną modelu Phinix-1 (100 nm) i zastępujemy ją dynamicznym licznikiem relacyjnym wyprowadzonym z geometrii mereologicznej. Wzmocnienie bliskiego pola przyjmuje postać

\[G = 1 + D^2/(L^2 + \delta^2)\]

gdzie

\[D_{linear} = \lambda_{res} / (2\pi^2 n) _ (1/3)\]

$\delta = a \cdot n \cdot 2$ (lub $(a+c)\cdot n \cdot 2$$ dla tetraedru SiO2), a czynnik $3^{-2}=1/9$ wynika z relacji mereologicznej synchronizacji na trzech osiach.

$\lambda_{res}$ jest długością fali nośnej informacji (rezonansu optycznego), a $n$ współczynnikiem załamania dla tej nośnej.

Model pokazuje, że dane Nature 2015 dla Si-Si (Cz-Si) nie odpowiadają czystemu krzemowi (rezonans wielofononowy ∼16-22 μm, $n=3.44$, predykcja ∼140x), lecz krzemowi zanieczyszczonemu tlenem międzywęzłowym (rezonans Si-O ∼9 μm, $n_{eff}=1.47$, predykcja ∼4000x, zgodna z pomiarem). SiO2 pozostaje zgodne z modelem tetraedrycznym $\delta=(a+c)\cdot n \cdot 2$.

Odkrycie implikuje falsyfikowalną predykcję: Float-Zone Si (beztlenowe) powinno wykazywać 3-5x mniejsze wzmocnienie NFRHT niż Czochralski Si i przesunięcie plateau odcięcia z 9 μm na 16-22 μm.

Abstract [EN]:

In this Annex 1 to Zarzewie 1, we remove the last phenomenological constant of the Phinix-1 model (100 nm) and replace it with a dynamic relational counter derived from mereological geometry. The near-field enhancement is

\[G = 1 + D^2/(L^2 + \delta^2)\]

where

\[D_{linear} = \lambda_{res} / (2\pi^2 n) _ (1/3)\]

$\delta = a \cdot n \cdot 2$ (or $(a+c)\cdot n \cdot 2$ for SiO2 tetrahedron), and the factor $3^{-2}=1/9$ arises from the mereological relation of synchronization on three axes.

$\lambda_{res}$ is the information carrier wavelength (optical resonance) and $n$ is the refractive index at the carrier.

We show that Nature 2015 Si-Si data (Cz-Si) does not correspond to pure Si (multiphonon resonance ∼16-22 μm, $n=3.44$, prediction ∼140x), but to oxygen-contaminated Si with interstitial Si-O resonance ∼9 μm ( $n_{eff}=1.47$, prediction ∼4000x, matching measurement). SiO2 remains consistent with tetrahedral $\delta=(a+c)\cdot n \cdot 2$.

The finding implies a falsifiable prediction: oxygen-free Float-Zone Si should exhibit 3-5x lower NFRHT enhancement than Czochralski Si with cutoff plateau shifted from 9 μm to 16-22 μm.

Draft for priority - full crystallographic analysis forthcoming.

PL

W temperaturze 300 K częstotliwości rezonansowe (mody optyczne sieci krystalicznej) dla Si oraz SiO₂ w zakresie podczerwieni różnią się zasadniczo ze względu na symetrię ich struktur krystalicznych. W spektroskopii podaje się je najczęściej w liczbach falowych ($\text{cm}^{-1}$).

1. Krzem (Si)

Czysty krzem krystaliczny (c-Si) ma strukturę diamentu o symetrii środkowosymetrycznej. Powoduje to, że jego podstawowy mod fononowy w centrum strefy Brillouina ($\sim 520\text{ cm}^{-1}$) jest nieaktywny w podczerwieni (jest aktywny tylko w spektroskopii Ramana). [1, 2]

W podczerwieni dla Si obserwuje się jednak rezonanse wynikające z absorpcji wielofononowej oraz zanieczyszczeń: [3]

$\sim 610 - 620\text{ cm}^{-1}$ ($\sim 16\text{ μm}$) – najsilniejszy pik rezonansowy związany z procesami wielofononowymi sieci (kombinacja fononów TO + TA).
$\sim 740\text{ cm}^{-1}$, $820\text{ cm}^{-1}$, $890\text{ cm}^{-1}$ – słabsze piki wielofononowe.
$\sim 1107\text{ cm}^{-1}$ ($\sim 9\text{ μm}$) – bardzo silny, charakterystyczny rezonans lokalny w krzemie hodowanym metodą Czochralskiego (CZ-Si), wywołany drganiami zanieczyszczeń tlenu międzywęzłowego (Si-O). [4, 5, 6]

2. Dwutlenek krzemu (SiO₂)

Dwutlenek krzemu (zarówno w formie kwarcu, jak i bezpostaciowej krzemionki/szkła kwarcowego) posiada silne, aktywne w podczerwieni mody optyczne (tzw. pasmo Reststrahlen). Wyróżnia się w nich rezonanse poprzeczne (TO) oraz podłużne (LO): [7, 8, 9]

$\sim 1075 - 1100\text{ cm}^{-1}$ ($\sim 9.2\text{ μm}$) – główny, najsilniejszy rezonans asymetrycznego rozciągania wiązań Si-O-Si (mod TO). Towarzyszy mu mod LO przy $\sim 1250\text{ cm}^{-1}$.
$\sim 800 - 810\text{ cm}^{-1}$ ($\sim 12.4\text{ μm}$) – rezonans symetrycznego rozciągania/zginania wiązań Si-O (mod TO).
$\sim 450 - 460\text{ cm}^{-1}$ ($\sim 22\text{ μm}$) – rezonans drgań kołyszących (rocking) oraz zginających Si-O-Si (mod TO). [7, 8, 9, 10, 11]

Wniosek: czysty c-Si nie posiada silnego rezonansu IR na 9μm. Pik 1107 cm⁻¹ jest sygnaturą tlenową Cz-Si. To podstawa hipotezy kontaminacji.

EN

At 300 K, the resonant frequencies (optical lattice modes) for Si and SiO2 in the infrared differ fundamentally due to crystal symmetry. They are given in wavenumbers ($\text{cm}^{-1}$).

1. Silicon (Si)

Pure crystalline Si (c-Si) has a diamond structure with centrosymmetric symmetry. Its fundamental phonon at Γ (~520 cm⁻¹) is IR-inactive (Raman-only) [1,2]. IR resonances are observed due to multiphonon absorption and impurities [3]:

~610-620 cm⁻¹ (~16 μm) – strongest multiphonon peak (TO+TA combination)
~740, 820, 890 cm⁻¹ – weaker multiphonon peaks
~1107 cm⁻¹ (~9 μm) – very strong local resonance in Czochralski-grown Si (Cz-Si), caused by interstitial oxygen (Si-O) [4,5,6]

2. Silicon Dioxide (SiO2)

SiO2 (quartz and amorphous silica) has strong IR-active optical modes (Reststrahlen band) [7,8,9]: * ~1075-1100 cm⁻¹ (~9.2 μm) – main TO resonance of asymmetric Si-O-Si stretching, LO at ~1250 cm⁻¹ * ~800-810 cm⁻¹ (~12.4 μm) – symmetric stretching/bending TO * ~450-460 cm⁻¹ (~22 μm) – rocking/bending TO [7,8,9,10,11]

Conclusion: Pure c-Si has no strong IR resonance at 9μm. The 1107 cm⁻¹ peak is the oxygen signature of Cz-Si. This is the basis for the O-contamination hypothesis for Nature 2015 Si-Si data.

References

[1-11] - to be completed (Cz-Si 1107 cm⁻¹, Reststrahlen SiO2)

[1] https://arxiv.org
[2] https://pubs.acs.org
[3] https://www.rp-photonics.com
[4] https://www.researchgate.net
[5] https://www.researchgate.net
[6] https://www.crystran.com
[7] https://www.researchgate.net
[8] https://www.researchgate.net
[9] https://www.researchgate.net
[10] https://www.researchgate.net
[11] https://infoscience.epfl.ch

Python

import numpy as np
import matplotlib.pyplot as plt

# 1. Definicja stałych fizycznych bazowych
alpha = 6.62607015e-34  # Stała Plancka h
c = 299792458          # Prędkość światła (m/s)
k = 1.380649e-23       # Stała Boltzmanna (J/K)
T = 300                # Temperatura (K)

E_classic = (8 * np.pi**5 * (k * T)**4) / (15 * c**3 * alpha**3)

def planck_integrand(freqs, alpha, c, k, T):
    exponent = (alpha * freqs) / (k * T)
    integrand = np.zeros_like(freqs)
    mask_standard = exponent <= 300
    if np.any(mask_standard):
        f = freqs[mask_standard]
        integrand[mask_standard] = (8 * np.pi * alpha * f**3) / (c**3 * (np.exp(exponent[mask_standard]) - 1))
    mask_wien = (exponent > 300) & (exponent < 700)
    if np.any(mask_wien):
        f = freqs[mask_wien]
        integrand[mask_wien] = (8 * np.pi * alpha * f**3 / c**3) * np.exp(-exponent[mask_wien])
    return integrand

# Dziedzina szerokości szczeliny L
L_plot = np.logspace(-10, -4, 500)

# Surowe stałe krystalograficzne i widmowe dla wyliczenia dynamicznego licznika i mianownika
# lambda_res to długość fali głównego rezonansu w próżni (m)
mats_data = {
    'Si':   {'a': 0.5431e-9,  'c': 0.5431e-9,   'n': 3.44, 'lambda_res': 9.0e-6},
    'Si_':  {'a': 0.5431e-9,  'c': 0.5431e-9,   'n': 1.47, 'lambda_res': 9.0e-6}, # Cz-Si: siatka Si + rezonans Si-O
    'Si16': {'a': 0.5431e-9,  'c': 0.5431e-9,   'n': 3.44, 'lambda_res': 16.0e-6},
    'SiO2': {'a': 0.4914e-9,  'c': 0.5405e-9,   'n': 1.47, 'lambda_res': 9.3e-6},
    'SiC':  {'a': 0.3076e-9,  'c': 1.0053e-9,   'n': 2.55, 'lambda_res': 10.6e-6},
    'GaN':  {'a': 0.3189e-9,  'c': 0.5186e-9,   'n': 2.35, 'lambda_res': 21.4e-6}
}

ratios_coupled = {mat: [] for mat in mats_data.keys()}

# Główna pętla numeryczna silnika Phinix-1
for L in L_plot:
    nu_limit_thermal = (700 * k * T) / alpha
    freqs_c = np.linspace(1e11, nu_limit_thermal, 15000)
    nu_cutoff_coupled = c / L
    base_integrand = planck_integrand(freqs_c, alpha, c, k, T)
    mask_near_field = freqs_c <= nu_cutoff_coupled

    for mat in mats_data.keys():
        d = mats_data[mat]

        # MIANOWNIK: Bezpiecznik atomowy wyliczony z krystalografii
        if mat == 'Si':
            delta = d['a'] * d['n'] * 2
        if mat == 'Si9':
            delta = d['a'] * d['n'] * 2
        elif mat == 'SiO2' or mat == 'SiC':
            delta = (d['a'] + d['c']) * d['n'] * 2
        elif mat == 'GaN':
            delta = d['c'] * d['n'] * 2 # Wybrany mod osiowy o najwyższym suficie

        # LICZNIK: Dynamiczny próg D zredukowany mereologicznym czynnikiem relacji 3^-2 (1/9)
        D_linear = d['lambda_res'] / (2 * (np.pi**2) * d['n']) * (1.0 / 3.0)
        D_squared = D_linear ** 2

        # Profil wzmocnienia wolny od stałych fenomenologicznych
        gain = 1 + D_squared / (L**2 + delta**2)

        gain_profile = np.ones_like(freqs_c)
        if np.any(mask_near_field):
            gain_profile[mask_near_field] = gain

        coupled_integrand = base_integrand * gain_profile
        ratios_coupled[mat].append(float(np.trapezoid(coupled_integrand, freqs_c) / E_classic))

# Dane empiryczne do wykresu walidacyjnego
L_mit = np.array([0.030, 0.050, 0.100])
ratio_mit = np.array([80.3, 32.1, 5.4])
L_natureSiO = np.array([0.002, 0.005, 0.010, 0.020, 0.040, 0.100])
ratio_natureSiO = np.array([6100.0/6.1, 3300.0/6.1, 1250.0/6.1, 410.0/6.1, 115.0/6.1, 32.1/6.1])
L_natureSi = np.array([0.002, 0.005, 0.010, 0.020, 0.040, 0.100])
ratio_natureSi = np.array([4200.0/6.1, 2100.0/6.1, 910.0/6.1, 310.0/6.1, 95.0/6.1, 28.5/6.1])

# Renderowanie wykresu
plt.figure(figsize=(9, 6.5))
colors = {'Si_': 'lightblue', 'Si16': 'royalblue', 'Si': 'blue', 'SiO2': 'cyan', 'SiC': 'orange', 'GaN': 'magenta'}
styles = {'Si_': 'dotted', 'Si16': '-.', 'Si': '-', 'SiO2': '-', 'SiC': '--', 'GaN': '--'}
labels = {'Si_': 'Six µm', 'Si16': 'Si 16 µm', 'Si': 'Si Theory', 'SiO2': 'SiO₂ Theory', 'SiC': 'SiC Prediction', 'GaN': 'GaN Prediction'}

for mat in mats_data.keys():
    plt.loglog(L_plot * 1e6, ratios_coupled[mat], color=colors[mat], ls=styles[mat], lw=2.0, label=labels[mat])

plt.scatter(L_mit, ratio_mit, color='black', marker='^', s=70, zorder=6, label='MIT Data')
plt.scatter(L_natureSiO, ratio_natureSiO, color='darkcyan', marker='o', s=45, zorder=5, label='Nature 2015 (SiO₂)')
plt.scatter(L_natureSi, ratio_natureSi, color='darkblue', marker='s', s=35, zorder=5, label='Nature 2015 (Si)')
plt.axhline(1.0, color='gray', ls='--', alpha=0.5, label='Planck Limit')

plt.title('ZRELATYWIZOWANA MACIERZ DYNAMICZNA PHINIX-1 (CZYNNIK 3^-2)', fontsize=11, fontweight='bold')
plt.xlabel('Gap size d [µm]')
plt.ylabel('Planck Limit Exceedance Factor')
plt.legend(loc='upper right', frameon=True, facecolor='white')
plt.grid(True, which="both", ls=":", alpha=0.5)
plt.xlim(0.0005, 5)
plt.ylim(0.5, 5000)
plt.tight_layout()

# Zapis pliku
output_filename = 'phinix1_dynamiczny_licznik_relacyjny.png'
plt.show()
# plt.savefig(output_filename, dpi=300)
# print(f"[OK] Wykres zapisany jako: {output_filename}")