Mechanika kwantowa i funkcje falowe — Schrödinger, tunelowanie, dwie szczeliny i obwody kwantowe

Mechanika kwantowa jest najprecyzyjniej sprawdzoną teorią fizyczną w historii — i jedną z najtrudniejszych do zinternalizowania, ponieważ jej przewidywania są fundamentalnie probabilistyczne, a jej obiekty (funkcje falowe) nie mają klasycznego odpowiednika. Ten przewodnik łączy każde kluczowe pojęcie kwantowe z działającą symulacją, budując intuicję krok po kroku od eksperymentu z dwiema szczelinami po kwantowy algorytm przeszukiwania Grovera.

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.

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)

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₀)

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⟩

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 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:

Algorytmy i metody

Solver TDSE metodą Cranka-Nicolson Gaussowski pakiet falowy Przybliżenie WKB Metoda macierzy transferu Współczynniki Clebscha-Gordana Harmoniki sferyczne Yₗᵐ Stowarzyszone wielomiany Laguerre'a Parametryzacja sfery Blocha Macierze bramek unitarnych (SU(2)) Przygotowanie stanów Bella Test nierówności CHSH Fazowy odrzut wyroczni Operator dyfuzji Grovera Wzmocnienie amplitudy Reguła Borna dla pomiaru kwantowego

Dalej: Devlog #27 — Aktualizacja platformy Q2 2027 omawia sprint treści fali 8, nowe symulacje i to, co dalej w planach rozwoju platformy.