Dlaczego mechanika kwantowa potrzebuje symulacji
Mechanika klasyczna pozwala wyobrazić sobie piłkę toczącą się ze wzgórza. Mechanika kwantowa opisuje elektron jako zespoloną funkcję falową ψ(x,t), której kwadrat modułu |ψ|² daje gęstość prawdopodobieństwa. Nie ma dla tego intuicyjnego obrazu myślowego. Funkcja falowa nie żyje w zwykłej przestrzeni — żyje w przestrzeni Hilberta. Gdy „mierzysz” elektron, funkcja falowa zapada się do określonej wartości w sposób, którego żaden proces klasyczny nie jest w stanie odtworzyć.
Symulacja pomaga, renderując to, co inaczej pozostałoby niewidoczne. Interaktywny pakiet falowy Schrödingera pozwala bezpośrednio obserwować superpozycję i interferencję. Tunelowanie kwantowe staje się wizualnie oczywiste, gdy widzisz, jak amplituda prawdopodobieństwa przecieka przez barierę, której cząstka klasycznie nie mogłaby przekroczyć. Algorytm Grovera, który na papierze wygląda jak czysta abstrakcyjna algebra, staje się zrozumiały, gdy obserwujesz, jak amplitudy wzmacniają się iteracyjnie w kierunku oznaczonego stanu.
Ten przewodnik przechodzi przez pięć kolejnych warstw teorii kwantowej, z których każda jest zakotwiczona w symulacji, którą możesz otworzyć i przetestować już teraz.
Część 1 — Funkcja falowa i równanie Schrödingera
Zależne od czasu równanie Schrödingera
Centralnym równaniem mechaniki kwantowej jest zależne od czasu równanie Schrödingera (TDSE). Rządzi ono ewolucją stanu kwantowego ψ(x,t) w czasie w obecności potencjału V(x). W przeciwieństwie do F=ma Newtona — równania różniczkowego dla trajektorii — TDSE jest równaniem różniczkowym dla amplitudy prawdopodobieństwa. Funkcja falowa ψ jest wartością zespoloną; kwadrat jej modułu |ψ(x,t)|² daje prawdopodobieństwo znalezienia cząstki w położeniu x w chwili t.
Zależne od czasu równanie Schrödingera
iℏ · ∂ψ/∂t = Ĥψ = −(ℏ²/2m)·∂²ψ/∂x² + V(x)·ψ
Gęstość prawdopodobieństwa: ρ(x,t) = |ψ(x,t)|² = ψ*·ψ
Normalizacja: ∫|ψ|² dx = 1 (całkowite prawdopodobieństwo = 1)
Stan stacjonarny: ψₙ(x,t) = φₙ(x)·exp(−iEₙt/ℏ)
Dla cząstki w pudełku: Eₙ = n²π²ℏ²/(2mL²), n = 1, 2, 3, …
Gaussowski pakiet falowy: ψ(x,0) = A·exp(−x²/4σ² + ik₀x)
Symulacja Schrödingera rozwiązuje numerycznie jednowymiarowe TDSE metodą Cranka-Nicolson, która jest bezwarunkowo stabilna i zachowuje normę (utrzymuje ∫|ψ|² = 1 na każdym kroku). Możesz umieścić gaussowski pakiet falowy — zlokalizowane wybrzuszenie prawdopodobieństwa o średnim pędzie ℏk₀ — i obserwować, jak się rozchodzi, rozpływa wskutek dyspersji i oddziałuje z barierami.
Pakiet falowy Schrödingera
1D TDSE rozwiązane metodą Cranka-Nicolson. Wizualizuje |ψ|², Re(ψ), Im(ψ). Umieszczaj bariery, studnie i potencjały harmoniczne. Obserwuj dyspersję, odbicie i superpozycję stanów własnych energii.
Eksperyment z dwiema szczelinami
Interferencja pojedynczej cząstki poprzez propagację funkcji falowej. Przełącz detekcję „którą drogą”, by zniszczyć wzór interferencyjny. Demonstruje dualizm korpuskularno-falowy i problem pomiaru.
Część 2 — Tunelowanie kwantowe
Obszary klasycznie zabronione
Mechanika klasyczna jest bezwzględna: cząstka w położeniu x o energii całkowitej E nie może wejść w obszar, gdzie energia potencjalna V(x) > E. Energia kinetyczna K = E − V byłaby ujemna, co jest niemożliwe dla cząstki klasycznej. Mechanika kwantowa nie zna takiego zakazu. Funkcja falowa może być niezerowa wewnątrz klasycznie zabronionego obszaru — zanika wykładniczo, ale nie znika całkowicie. Jeśli bariera jest wystarczająco cienka, funkcja falowa pojawia się po drugiej stronie ze skończoną amplitudą. To jest tunelowanie kwantowe.
Tunelowanie nie jest rzadkim, egzotycznym efektem. To mechanizm leżący u podstaw rozpadu alfa jąder promieniotwórczych, działania diod tunelowych, skaningowej mikroskopii tunelowej (STM) oraz enzymatycznego transferu wodoru w biochemii. Samo Słońce zależy od tunelowania kwantowego: bez niego fuzja proton-proton byłaby niemożliwa w temperaturach jądra Słońca (≈15 milionów K to o wiele za mało, by klasycznie pokonać barierę kulombowską).
Tunelowanie kwantowe — współczynnik transmisji
Bariera: V(x) = V₀ dla 0 ≤ x ≤ a, V = 0 poza nią
Stała zaniku: κ = √(2m(V₀−E)) / ℏ (E < V₀)
Transmisja: T ≈ exp(−2κa) (cienka bariera, κa ≫ 1)
Przybliżenie WKB: T ≈ exp(−2∫κ(x) dx) (dowolny kształt)
Czynnik Gamowa (rozpad alfa): Γ ∝ exp(−2∫√(2m(V(r)−E))/ℏ dr)
Prąd STM: I ∝ exp(−2κd) (d = odległość ostrze-próbka)
Tunelowanie kwantowe
Pakiet falowy uderza w barierę potencjału. Animowane |ψ|² pokazuje składowe odbitą i przepuszczoną. Współczynnik transmisji T w funkcji szerokości i wysokości bariery. Nakładka przybliżenia WKB.
Schrödinger: presety bariery i studni
Skorzystaj z presetów stanu związanego i bariery schodkowej, by zbadać kwantowanie energii w studni potencjału — oraz przejście od bariery odbijającej do warunku rezonansowej transmisji.
Skaningowy mikroskop tunelowy osiąga rozdzielczość przestrzenną poniżej angstremu, ponieważ prąd tunelowy spada mniej więcej o rząd wielkości na każdy 1 Å wzrostu odległości ostrza od próbki. Ta niezwykła czułość — wykładnicza zależność od szerokości szczeliny — zamienia topografię każdego atomu bezpośrednio w fluktuacje prądu, które można odczytać wzmacniaczem synchronicznym. Symulacja tunelowania kwantowego uchwytuje tę fizykę w jednym wymiarze z regulowaną szerokością i wysokością bariery.
Część 3 — Struktura atomowa: orbitale wodoru
Liczby kwantowe wodoru
Atom wodoru to jedyny atom, dla którego równanie Schrödingera ma dokładne rozwiązanie analityczne. Rozwiązanie trójwymiarowego TDSE we współrzędnych sferycznych daje zbiór funkcji falowych ψ_nlm(r,θ,φ) oznaczonych trzema liczbami kwantowymi: n (główna, ustala energię), l (orbitalny moment pędu, 0 ≤ l ≤ n−1) oraz m (magnetyczna liczba kwantowa, −l ≤ m ≤ l). Każda unikalna trójka (n,l,m) to orbital — przestrzenny rozkład prawdopodobieństwa dla elektronu.
Funkcje falowe wodoru
ψₙₗₘ(r,θ,φ) = Rₙₗ(r) · Yₗᵐ(θ,φ)
Energia: Eₙ = −13,6 eV / n² (n = 1, 2, 3, …)
Rₙₗ(r): część radialna — wielomiany Laguerre'a × exp(−r/na₀)
Yₗᵐ(θ,φ): harmoniki sferyczne (część kątowa)
Promień Bohra: a₀ = ℏ²/(mₑe²) ≈ 0,529 Å
Orbital 1s (n=1,l=0,m=0): ψ = (1/√π)·a₀^{−3/2}·exp(−r/a₀)
Orbitale wodoru
Wizualizacja 3D orbitali ψₙₗₘ do n=4. Mapy kolorów |ψ|² (gęstość prawdopodobieństwa) i arg(ψ) (faza). Suwaki dla n, l, m. Przekrój poprzeczny ujawnia powierzchnie węzłowe.
Tunelowanie kwantowe
Tunelowanie leży u podstaw rozpadu alfa: jądro helu formuje się wewnątrz jądra uranu i tuneluje przez barierę kulombowską. Ustaw parametry bariery zgodnie z czynnikiem Gamowa i obserwuj reakcję prawdopodobieństwa transmisji.
Kształty orbitali atomowych determinują wiązania chemiczne. Orbital 2p ma dwa płaty rozdzielone płaszczyzną węzłową — ta geometria wymusza tworzenie wiązań sigma wzdłuż osi międzyjądrowej. Złożone geometrie płatowe orbitali 3d tworzą rozszczepienie pola krystalicznego odpowiedzialne za barwy związków metali przejściowych. Zrozumienie kształtów orbitali jako rozkładów prawdopodobieństwa (a nie „powłok”) to koncepcyjny krok od modelu Bohra do współczesnej chemii kwantowej.
Część 4 — Obliczenia kwantowe: kubity i bramki
Od bitów klasycznych do kubitów
Bit klasyczny to albo 0, albo 1. Kubit to układ kwantowy, który może istnieć w superpozycji |0⟩ i |1⟩: |ψ⟩ = α|0⟩ + β|1⟩, gdzie |α|² + |β|² = 1. Stan kubitu to wektor jednostkowy w dwuwymiarowej zespolonej przestrzeni Hilberta — wizualizowany jako punkt na sferze Blocha. Bramki kwantowe to operacje unitarne obracające ten wektor i są kwantowym odpowiednikiem klasycznych bramek logicznych.
Stan kubitu i bramki jednokubitowe
Stan: |ψ⟩ = α|0⟩ + β|1⟩, gdzie |α|² + |β|² = 1
Wektor Blocha: r = (2Re(α*β), 2Im(α*β), |α|²−|β|²)
Hadamard: H = (1/√2)[[1,1],[1,−1]] → |0⟩ → (|0⟩+|1⟩)/√2
Pauli-X (NOT): X = [[0,1],[1,0]] → |0⟩ ↔ |1⟩
Faza: S = [[1,0],[0,i]] → dodaje fazę 90° do |1⟩
Bramka T: T = [[1,0],[0,exp(iπ/4)]] → faza 45° do |1⟩
CNOT (2-kubitowa): odwraca kubit docelowy, jeśli kontrolny = |1⟩
Kubit i sfera Blocha
Interaktywna sfera Blocha 3D. Stosuj bramki jednokubitowe (H, X, Y, Z, S, T) i obserwuj obrót wektora stanu. Dziennik bramek pokazuje sekwencję. Wizualizuje superpozycję, fazę i kolaps przy pomiarze.
Splątanie kwantowe
Stany Bella (Φ+, Φ−, Ψ+, Ψ−) z animowaną wiązką splątania. Pomiar probabilistyczny, test nierówności CHSH (S = 2√2 ≈ 2,83). Demonstruje nielokalne korelacje, których systemy klasyczne nie są w stanie odtworzyć.
Moc obliczeń kwantowych pochodzi ze splątania: dwóch splątanych kubitów nie można opisać niezależnie — ich wspólny stan ma korelacje, których żaden klasyczny rozkład prawdopodobieństwa nie potrafi odtworzyć. To właśnie dowodzi twierdzenie Bella: nierówność CHSH S ≤ 2 musi obowiązywać dla dowolnej teorii lokalnych zmiennych ukrytych, ale mechanika kwantowa przewiduje — a eksperyment potwierdza — S = 2√2 ≈ 2,83. Symulacja splątania kwantowego uruchamia test CHSH na żywo, który możesz zweryfikować, zmieniając kąty pomiaru.
Część 5 — Kwantowy algorytm przeszukiwania Grovera
Kwadratowe przyspieszenie dla przeszukiwania nieustrukturyzowanego
Mając nieuporządkowaną bazę danych N elementów z jednym oznaczonym elementem, komputer klasyczny wymaga średnio O(N) zapytań, by go znaleźć. Algorytm Grovera rozwiązuje ten sam problem w O(√N) zapytań kwantowych — kwadratowe przyspieszenie osiągnięte dzięki wykorzystaniu kwantowego wzmocnienia amplitudy. Algorytm wymaga jedynie około π√N/4 iteracji, by osiągnąć niemal pewne prawdopodobieństwo sukcesu.
Każda iteracja algorytmu Grovera składa się z dwóch operacji: (1) wyroczni, która odwraca znak (fazę) amplitudy oznaczonego stanu, i (2) operatora dyfuzji, który odbija wszystkie amplitudy względem ich średniej. Efekt netto to zwiększenie amplitudy oznaczonego stanu o w przybliżeniu 2/√N na krok. Po ~π√N/4 krokach oznaczony stan ma niemal jednostkową amplitudę, a pomiar znajduje go z wysokim prawdopodobieństwem.
Algorytm Grovera
Stan początkowy: |s⟩ = (1/√N)·Σ|x⟩ (jednolita superpozycja)
Wyrocznia Uω: |x⟩ → −|x⟩ jeśli x=ω, |x⟩ → |x⟩ w przeciwnym razie
Dyfuzja D = 2|s⟩⟨s| − I (inwersja względem średniej)
Iteracja: |ψₜ₊₁⟩ = D·Uω·|ψₜ⟩
Po T krokach: P(znalezienia ω) = sin²((2T+1)·arcsin(1/√N))
Optymalne T ≈ π√N/4 → P ≈ 1 − O(1/N)
Złożoność: O(√N) w porównaniu z klasycznym O(N) — przyspieszenie kwadratowe
Algorytm przeszukiwania Grovera
Animacja krok po kroku dla N=16 stanów. Wyrocznia oznacza cel; każda iteracja Grovera pokazana jest jako histogram amplitud z animacją lerp. Porównaj O(√N) zapytań kwantowych z O(N) klasycznymi.
Kubit i bramki — wiedza wstępna
Zanim uruchomisz algorytm Grovera, zbuduj tutaj intuicję dla bramek kwantowych. Bramka Hadamarda, która generuje początkową jednolitą superpozycję, to bramka H zastosowana do |0⟩ — zobacz jej działanie bezpośrednio na sferze Blocha.
Algorytm Grovera demonstruje kluczową zasadę przewagi kwantowej: nie wymaga, by baza danych miała jakąkolwiek strukturę — działa na dowolnej wyroczni, która potrafi zidentyfikować oznaczony element. Technika wzmocnienia amplitudy, którą wprowadza, jest także podprogramem w potężniejszych algorytmach kwantowych. W symulacji możesz zmieniać rozmiar bazy danych (N), obserwować sinusoidalną ewolucję prawdopodobieństwa i zobaczyć, że zatrzymanie się dokładnie na optymalnej liczbie iteracji ma znaczenie — wykonaj zbyt wiele kroków, a prawdopodobieństwo zacznie ponownie spadać.
Powiązania kwantowe w całej kolekcji
Mechanika kwantowa leży u podstaw symulacji w wielu innych kategoriach platformy:
- Promieniowanie ciała doskonale czarnego — kwantowa hipoteza Plancka (E=hν) rozwiązała katastrofę w nadfiolecie fizyki klasycznej i położyła fundament pod mechanikę kwantową.
- Masa-energia (E=mc²) — rozszczepienie i fuzja jądrowa to procesy kwantowe; krzywa energii wiązania jest obliczana z półempirycznego wzoru masowego Bethego-Weizsäckera.
- Fuzja jądrowa — łańcuch proton-proton w Słońcu wymaga tunelowania kwantowego przez barierę kulombowską.
- Ruchy Browna — w nanoskali fluktuacje termiczne stają się porównywalne z kwantową energią stanu podstawowego; długość fali de Broglie'a λ = h/mv zaczyna mieć znaczenie.
- Struktury krystaliczne — typy wiązań metalicznych, jonowych i kowalencyjnych w symulacji struktur krystalicznych mają wszystkie pochodzenie mechaniczno-kwantowe.
Algorytmy i metody
Dalej: Devlog #27 — Aktualizacja platformy Q2 2027 omawia sprint treści fali 8, nowe symulacje i to, co dalej w planach rozwoju platformy.