Domanda:
Perché l'equazione di Schrödinger funziona così bene per l'atomo di idrogeno nonostante il confine relativistico al nucleo?
Paul Young
2019-06-01 20:23:24 UTC
view on stackexchange narkive permalink

Mi è stato insegnato che le condizioni al contorno sono importanti tanto quanto l'equazione differenziale stessa quando si risolvono problemi fisici reali.

Quando l'equazione di Schrödinger viene applicata all'atomo di idrogeno idealizzato, è separabile e le condizioni al contorno vengono applicate alla componente radiale.Sono preoccupato per il confine $ r = 0 $ vicino al nucleo.Vicino al protone, l'energia cinetica dell'elettrone sarà relativistica e guardare la stessa equazione di Schrödinger per come dovrebbe comportarsi questo confine sembra pericoloso perché il suo termine di energia cinetica è solo un'approssimazione non relativistica.

C'è qualche intuizione fisica, o matematica, che posso osservare che dovrebbe farmi sentire a mio agio con le condizioni al contorno in questa regione?

Quattro risposte:
RogerJBarlow
2019-06-01 21:32:38 UTC
view on stackexchange narkive permalink

Nella risoluzione dell'equazione radiale di Schroedinger non viene applicata alcuna condizione al contorno a $ r = 0 $ . A $ r = \ infty $ sì, $ R (r) $ deve tendere a zero, quindi rifiutare la soluzione esponenziale positiva; qualsiasi cambiamento in questo avrebbe conseguenze enormi. Ma non ci sono vincoli su $ R (r) $ o addirittura su $ R '(r) $ come $ r \ a 0 $ .

Quindi non c'è un cambiamento nella condizione al contorno. C'è un cambiamento nelle energie cinetiche e potenziali dovuto a effetti relativistici e al fatto che il protone non è una carica puntiforme. Questi hanno un effetto, ma molto piccolo, poiché il volume in questione è circa $ 10 ^ {- 15} $ del volume dell'atomo. (In realtà gli esperimenti dei fisici atomici possono rilevare questi effetti, almeno per i grandi $ Z $ atomi, grazie ad alcuni esperimenti ottici molto intelligenti e precisi.) Ma questo è un piccolo effetto , non il punto di svolta che una nuova condizione al contorno potrebbe dare.

Quindi, immagina di adottare un approccio numerico per risolvere l'equazione differenziale radiale con $ l = 0 $.Io "sparerò" da $ r = \ infty $ verso l'interno a vari livelli di energia cercando i diversi autovalori per gli orbitali s.Se non esiste una condizione al contorno a $ r = 0 $, come faccio a sapere quali livelli di energia sono soluzioni?
Cosa intendo per metodo di "tiro": https://en.wikipedia.org/wiki/Shooting_method
Suggerimento: banalmente rilevabile per l'oro.La meccanica quantistica non relativistica produce un oro grigio, non un oro giallo.
@PaulYoung prova a prendere il limite di $ \ varepsilon \ to0 $ con il potenziale $ \ frac1 {\ sqrt {r ^ 2 + \ varepsilon ^ 2}} $ invece del solito $ 1 / r $.Questo ti darà il risultato senza imporre una condizione al contorno speciale.Vedi anche [questa mia domanda] (https://physics.stackexchange.com/q/292593/21441) e la mia auto-risposta lì.
@ruslan - perché non provare una risposta invece di un semplice commento?... mi sembra che tu mi capisca?
@PaulYoung Non sono sicuro di quello che stai veramente chiedendo: nel secondo paragrafo del tuo OP stai parlando dell'operatore di energia cinetica che è un'approssimazione non relativistica, che è irrilevante per le condizioni al contorno (per approssimare un po 'più da vicino userestiun'approssimazione di ordine superiore di $$ T = \ sqrt {m ^ 2c ^ 4 + p ^ 2c ^ 2} -mc ^ 2 = \ frac {p ^ 2} {2m} - \ frac {p ^ 4} {8m ^3c ^ 2} + \ frac {p ^ 6} {16m ^ 5c ^ 4} + ..., $$ non cambia una condizione al contorno) e nell'ultimo paragrafo stai chiedendo della condizione al contorno a $ r = 0$, che sembra essere affrontato dal mio link nel commento precedente.
Quello che penso che @PaulYoung stia ottenendo (e con cui sono d'accordo) è che l'affermazione "non c'è vincolo su $ R (r) $ o $ R '(r) $ come $ r \ a 0 $".La funzione d'onda deve essere finita come $ r \ a 0 $.Ciò funge effettivamente da condizione al contorno per l'ODE radiale, poiché la maggior parte delle soluzioni all'ODE radiale non sono finite all'origine.
@MichaelSeifert - in effetti ... so che non può divergere più velocemente di 1 / r (perché deve essere normalizzabile) ma questo lascia ancora aperto 1 / r, finito non zero e zero .. e poi forse succede qualcosa di strano primaandare troppo lontano - sebbene Ruslan anna_v e altri sostengano che la funzione d'onda si comporterà bene come testimoniato, ad esempio, dalle soluzioni dell'equazione di Dirac
@PaulYoung bene, se chiami illimitatezza comportati bene ...
tparker
2019-06-02 01:53:52 UTC
view on stackexchange narkive permalink

Bella domanda. La tua affermazione che

vicino al protone l'energia cinetica dell'elettrone sarà relativistica

non è così semplice come potrebbe sembrare. L'energia cinetica dell'elettrone $ \ langle \ hat {T} \ rangle = \ langle \ hat {p} ^ 2 \ rangle / (2m) $ è una quantità non locale che può essere espresso in modo equivalente come uno dei due integrali

$$ \ langle \ hat {T} \ rangle = \ frac {1} {2m} \ int d ^ 3x \ \ psi ^ * (x) \ left (- \ hbar ^ 2 \ nabla ^ 2 \ right) \ psi (x) = \ frac {1} {2m} \ int d ^ 3x \ | \ hbar \, {\ bf \ nabla} \ psi (x) | ^ 2 . $$

Quindi l'energia cinetica dell'elettrone "in" una posizione particolare non è ben definita; potrebbe essere il valore di uno dei due integrandi sopra in quel punto (o, in effetti, di qualsiasi altro integrando che si integra allo stesso valore su tutto lo spazio).

L'ultima espressione è la più naturale da usare, però, perché almeno è positiva-semidefinita. Abbiamo ancora il problema che $ \ hbar ^ 2 | \ nabla \ psi (0) | ^ 2 / (2m) $ è una "densità di energia cinetica" (qualunque cioè) piuttosto che un'energia cinetica reale, quindi non possiamo parlare di quanto sia relativistico l'elettrone "al" nucleo. (Potremmo integrare la dimensione empirica del nucleo, ma non credo che sia proprio questo il motivo per cui la tua domanda sta arrivando: non stai chiedendo quando l'elettrone è letteralmente all'interno del nucleo, ma quando è abbastanza vicino al potenziale centro che si muove intuitivamente molto rapidamente.)

Ma niente di tutto questo ha davvero importanza: il punto è che poiché l'integrando è definito positivo, il contributo all'energia cinetica su qualsiasi regione particolare è sempre inferiore (o uguale) all'energia cinetica totale su ogni regione.Quindi, per verificare in modo significativo se gli effetti relativistici devono essere presi in considerazione, è necessario calcolare l'energia cinetica totale su tutto lo spazio.Questo risulta essere $ \ hbar ^ 2 / (2m a ^ 2) = me ^ 4 / (2 \ hbar ^ 2) = (\ alpha ^ 2/2) mc ^2 $ , dove $ \ alpha $ è la costante di struttura fine.Gli effetti relativistici sono trascurabili se l'energia cinetica è molto inferiore all'energia a riposo dell'elettrone, che corrisponde alla condizione che $ \ alpha ^ 2/2 = 1/37538 \ ll 1 $ , il che, rassicurante, è vero.

Da ciò capisco che l'energia cinetica media o "attesa" dell'elettrone non è relativistica.Tuttavia, l'equazione di Schrödinger mi sembra "locale" e penso che sia ragionevole pensare che abbia un valore in ogni punto.Questo è ciò che rende la funzione d'onda "wigglier" vicino al nucleo, giusto?
Che sia più vicino al nucleo dipende dall'orbitale: per alcune scelte di $ n $ e $ l $, la funzione d'onda è piatta al nucleo.Ad ogni modo, come proponi di convertire la densità di energia cinetica data dal gradiente della funzione d'onda in un'energia effettivamente cinetica, che ti dirà se l'elettrone è relativistico?
Propongo di applicare localmente l'operatore di quantità di moto e dividere il risultato per la massa a riposo e confrontarlo con la velocità della luce
Considererei anche di applicare $ del ^ 2 $ e fare qualcosa di simile con $ mv ^ 2 $ - (scusa se non so come scrivere del in Tex)
Se lo fai, le unità non funzionano.Questo è quello che continuo a cercare di dirti;l'operatore quantità di moto restituisce solo una quantità con le unità di quantità di moto se la si integra nello spazio.(Inoltre, del è "\ nabla" in TeX.)
+1 su tutto quello che posso, e ci sto ancora lavorando ...
se non riesco ancora a ottenerlo aggiungerò una taglia
posso applicare l'operatore di quantità di moto alla funzione d'onda e quindi integrare il risultato su una piccola regione di spazio vicino al nucleo e quindi dividere quel risultato per la probabilità che l'elettrone si trovi in quella piccola regione?
@PaulYoung Il significato fisico di questo rapporto non è chiaro per me.Sarebbe molto grande in qualsiasi posizione in cui è molto improbabile che si trovi l'elettrone, e non vedo perché corrisponderebbe all'energia cinetica "in" quella posizione.Semplicemente non penso che esista un modo locale naturale per definire l'energia cinetica "in" una posizione all'interno della configurazione del QM non relativistico.(Nell'impostazione QFT, puoi parlare del tensore energia-momento.)
anna v
2019-06-01 20:46:31 UTC
view on stackexchange narkive permalink

Le condizioni al contorno che raccolgono le funzioni d'onda dell'idrogeno sono i "vincoli" posti alle soluzioni della funzione d'onda.Ricorda che l'osservabile è la distribuzione di probabilità da $ Ψ ^ * Ψ $ , non una posizione particolare.Si prega di leggere il collegamento.Dopo tutto, le soluzioni sono all'interno del postulato della meccanica quantistica.

Non esiste alcuna condizione al contorno relativistica, perché non esistono orbite, ma solo distribuzioni di probabilità.

Quindi le soluzioni non hanno una singolarità per r = 0, e in generale c'è una piccola probabilità di trovare l'elettrone all'origine, se i numeri quantici consentono un'interazione, come con cattura elettronica nei nuclei.Per l'atomo di idrogeno non c'è abbastanza energia per far apparire un neutrone.

Sono un po 'preoccupato per i collegamenti GSU perché sembrano essere basati su un presupposto di non relatività, ma questo è esattamente ciò di cui sono preoccupato
Ma le stesse soluzioni per quanto riguarda i livelli di energia ecc. Derivano dalla risoluzione dell'atomo di idrogeno con l'equazione di Dirac, che è relativistica.https://quantummechanics.ucsd.edu/ph130a/130_notes/node501.html.Ciò è dovuto alla natura probabilistica della meccanica quantistica, dove i vincoli ei postulati sono quelli che determinano le soluzioni.Non vi è alcun effetto di "avvicinarsi a r = 0" nelle distribuzioni di probabilità in entrambi i casi.
Non capisco l'equazione di Dirac, ma sembra dalla tua risposta e dal Prof Barlow che poiché la condizione al contorno $ r = 0 $ non deriva da nulla che abbia a che fare con l'equazione, non c'è molto dapreoccuparti oltre alla piccola correzione energetica che dovrebbe essere su una scala di "volume"
il collegamento che fornisco sopra dice "Questo risultato fornisce la stessa risposta del nostro calcolo non relativistico per ordinare a ^ 4 ma è anche corretto per un ordine superiore. È una soluzione esatta al problema della meccanica quantistica posto ma non include gli effetti diteoria dei campi, come lo spostamento di Lamb e il momento magnetico anomalo dell'elettrone ".
collegamento dall'aspetto impegnativo, ma vedo cosa posso fare con esso
my2cts
2019-06-01 22:52:25 UTC
view on stackexchange narkive permalink

La condizione al contorno per r = 0 è che la funzione d'onda dovrebbe essere finita.L'equazione di Schrödinger per gli atomi UC di idrogeno e probabilmente tutti gli atomi ha soluzioni con $ \ cal l $ negativo, che vengono rifiutate perché divergono in r = 0.Vedi ad esempio il libro di testo di Schiff sulla meccanica quantistica.

Per quanto riguarda gli effetti relativistici, potresti voler confrontare le espressioni di energia dell'idrogeno per Dirac, better, Klein-Gordon - no spin e Schrödinger.Dai un'occhiata a un altro ottimo testo, Itzykson e Zuber, per questi.

Controllerò i testi di riferimento.Schiff, in particolare, sembra che ci siano alcuni commenti che non ho mai visto.


Questa domanda e risposta è stata tradotta automaticamente dalla lingua inglese. Il contenuto originale è disponibile su stackexchange, che ringraziamo per la licenza cc by-sa 4.0 con cui è distribuito.
Loading...