43. DFT. Wprawka I. Energia stanu podstawowego atomu wodoru z funkcjonału Weizsäckera

    W książce [1] znajduje się wyprowadzenie energii stanu podstawowego atomu wodoru i jonu H₂+ w formie zależności tej energii od gęstości elektronowej. Jest to zatem element teorii funkcjonałów gęstości (DFT). Po wyprowadzeniu równanie uzyskało postać:

    W innych pracach pierwsza całka ma jednak nieco inną postać, którą nazywa się funkcjonałem energii kinetycznej Weizsäckera T_W:

    Różnica polegająca na tym, że we wzorze pierwszym mianownik jest pod kwadratem powoduje kłopoty z całkowaniem. Wystarczy zauważyć, że dla energii stanu podstawowego atomu wodoru znamy ścisły wzór na unormowaną gęstość elektronową w jednostkach atomowych:

   Łatwo się przekonać, że omawiany stosunek z pierwszego równania, będący kwadratem ilorazu gradientu gęstości przez gęstość, musi dać, dla dokładnej postaci równania gęstości, wartość stałą, a to powoduje, że całkowanie, tak zapisanego funkcjonału Weizsäckera, staje się rozbieżne. Należy przyjąć, że kwadrat dotyczy tylko licznika. Poprawiony wzór ma zatem postać:

    Stosunek kwadratu gradientu gęstości elektronowej do gęstości ma postać:

   Pamiętając, że całkowanie odbywa się we współrzędnych sferycznych (jakobian przekształcenia), całkowanie po kątach daje współczynnik 4pi. Przyjmując, że potencjał V_N(r) jest dany wzorem -1/r, ostateczny wzór na energię całkowitą stanu podstawowego atomu wodoru przyjmuje następującą postać:

    

    Wynik jest prawidłowy, gdyż dokładna wartość energii, obliczona dla porównania z równania Schrödingera wynosi dokładnie tyle samo (czyli -1/2) w jednostkach hartree (atomowych).

Literatura: 1. N. H. March, Applications of density functional theory in chemistry, in Density Functional Theory, Ed. E.K.U. Gross and R. M. Dreizler, Plenum Press, New York, 1995 p. 299.

42. Wartości minimalne i maksymalne zmiennej chaotycznej dla kwadratu funcji logarytmicznej w zależności od stałej c

    Ciekawiło mnie, jak zmieniają się wartości minimalne i maksymalne zmiennej chaotycznej x dla szeregu danego wzorem 1.

w zakresie chaotycznym dla stałej c od 0.1 do 1.84. W tym celu napisałem program realizujący to zagadnienie:

'xmin i x max dla obszaru chaotycznego f. log^2

CLS : SCREEN 12

x = .7

testmax = 0

testmin = 100000

OPEN "minmax.txt" FOR APPEND AS #1

FOR c = .1 TO 1.84 STEP .01

 llos = 10000  'liczba iteracji

     'dla c=1.2 min=1.15E-9, max=507.92

 FOR i = 1 TO llos

  x = c * LOG(x) ^ 2  'funkcja logarytmiczna

  IF x >= testmax THEN testmax = x

  IF x <= testmin THEN testmin = x

 NEXT i

 PRINT #1, c, testmin, testmax

NEXT c

    Wyniki przedstawiam w postaci dwóch wykresów ze względu na różnicę wielu rzędów wielkości dla wartości maksymalnej i minimalnej. Pierwszy posiada skalę liniową i pokazuje skokową zmianę wartości maksymalnej przy przechodzeniu od jednej wartości zmiennej c do innej. 

    Drugi wykres ze skalą logarytmiczną eksponuje podobne zmiany dla wartości minimalnej. Kolor niebieski oznacza wartości zmiennej c. Liczby na osi odciętych pokazują numer zmiennej c a nie jej wartość (ot Excel).

    Zmiany te również są skokowe, ale niekoniecznie w tych  samych miejscach, co skoki wartości maksymalnej.