## METODE ELEMEN HINGGA UNTUK PENYELESAIAN PERSAMAAN SCHRÖDINGER ATOM HIDROGENIK
Paken Pandiangan (pakenp@mail.ut.ac.id) Universitas Terbuka
Supriyadi (adyph_0901@yahoo.com) A. Arkundato (a_arkundato@yahoo.com) Universitas Jember

## ABSTRACT
The research computed the energy levels and radial wave functions of the  Hydrogen Atom. The method used for computation was FEM (finite element method). Using the variational method approach, FEM was applied to the action integral of  Schrödinger equation. This lead to the eigenvalue equation in the form of  global matrix equation. The results of computation were depended on boundary of the action integral of Schrödinger equation and number of elements. For boundary 0 - 100a 0 and 100 elements,  they were the realistic and best choice of computation to the closed  analytic results. The computation of first five energy levels resulted E 1 = -0.99917211 R ∞ , E 2 = -0.24984445 R ∞ , E 3 = -0.11105532 R ∞ , E 4 = -0.06247405 R ∞ and  E 5 = -0.03998598 R ∞ where 1 R ∞ = 13.6 eV. They had relative error under 0.1% to the analytic results.
Key words: FEM, Hydrogen atom, Schrödinger equation .
Telaah mikroskopis atom hidrogen dengan menggunakan persamaan Schrödinger sejauh ini telah dikerjakan dengan baik secara analitik (Nikishkov, 2004). Untuk atom multi elektron, maka solusi analitik yang diinginkan sebagai pemecahan persamaan Schrödingernya menjadi makin rumit diperoleh dan memerlukan aproksimasi lebih lanjut. Salah satu aproksimasi yang dapat diadopsi adalah mencari solusi numerik menggunakan metode numerik yang ada. Metode numerik pemecahan persamaan diferensial (misalnya persamaan Schrödinger) yang umum digunakan adalah metode finite difference (FD).  Metode numerik yang lain adalah metode elemen hingga ( finite element methods, FEM) yang saat ini mulai banyak dieksplorasi dan berkembang cepat untuk diterapkan pada problem-problem fisika. Salah satu kelebihan FEM adalah dapat diterapkan untuk problem fisis dimana bentuk potensialnya (dalam persamaan Schrödinger) sangat kompleks dan kurang simetris (Brugioni, 1995). Pada riset yang kami lakukan, sebagai riset pendahuluan, kami mencoba menggunakan metode FEM ini untuk atom hidrogen yaitu menghitung tingkat-tingkat energinya secara first principle , yaitu memecahkan persamaan Schrödinger untuk atom hidrogen. Hasil numerik yang diperoleh kemudian dibandingkan dengan hasil analitik yang sudah ada.

## Solusi Analitik Persamaan Schrödinger Atom Hidrogen
Persamaan Schrödinger bebas waktu untuk atom Hidrogen dimana potensialnya dinyatakan dalam satuan CGS adalah
) ( ) ( ) ( 2 2 2 2 r E r r e r m e            (1)
Dalam sistem koordinat bola, persamaan diatas dapat dinyatakan dengan
) , , ( sin 1 sin sin 1 1 2 2 2 2 2 2 2 2 φ r E ψ(r,θ,φ) r e φ θ r θ θ θ θ r r r r r e m                                           
(2)
Penyelesaian lengkap secara analitik secara panjang lebar telah banyak dibahas dalam buku teks (Mohan, 2002), yang menghasilkan solusi untuk atom hidrogen adalah:
,... 3 , 2 , 1 6 , 13 2 2 2 0 2      n eV n n a e E n (3)

## Integral Aksi
Dinamika sistem fisis, baik klasik maupun kuantum sebenarnya mengikuti rumusan persamaan diferensial yang dapat diturunkan dari prinsip aksi stasioner yang bersifat umum. Dalam hal ini, Euler dan Lagrange telah memberikan suatu framework untuk merealisasikan prinsip ini dalam memecahkan problem fisika. Di dalam mekanika klasik, maka azas aksi stasioner menyatakan untuk sistem konservatif integral berikut
dt t q q L b a t t ) , , (     (4)
bersifat stasioner (Brugioni, 1995). L adalah Langrangian dan variabel Λ disebut dengan integral aksi atau singkatnya disebut aksi. Berdasarkan syarat batas (BC = B oundary Conditions ) integral aksi Λ diambil pada nilai ektremum menurut syarat
0    (5)
Evaluasi integral aksi sangat membantu kita dalam mengkaji dinamika sistem kuantum, dalam hal ini persamaan Schrödinger, dengan menggunakan FEM.

## Metode Elemen Hingga (FEM)
Konsep dasar yang melandasi FEM adalah prinsip diskritisasi. Secara umum, diskritisasi dapat diartikan sebagai upaya untuk membagi sistem dari problem yang akan diselesaikan (obyek) menjadi bagian-bagian yang lebih kecil. Bagian-bagian yang lebih kecil tersebut selanjutnya disebut sebagai elemen hingga. Diskritisasi ini muncul karena adanya kesulitan untuk mempelajari sistem secara keseluruhan.
Pendekatan yang sering digunakan untuk merumuskan suatu problem fisis ke dalam analogi diskritisasi FEM yaitu metode Galerkin dan metode variasi (Gasiorowicz, 1995). Perumusan FEM dengan menggunakan metode Galerkin didasarkan pada minimasi residu (sisa) yang tertinggal setelah suatu solusi pendekatan disubsitusikan ke dalam persamaan diferensial dari suatu problem fisis. Sedangkan metode variasi didasarkan pada nilai-nilai stasioner dari suatu fungsi atau besaran skalar yang berkaitan dengan suatu problem fisis untuk mendapatkan solusinya (Thankappan, 1985).
Perumusan FEM dalam penelitian ini dipilih pendekatan metode variasi dengan menggunakan integral aksi yang diperoleh dari prinsip aksi stasioner. Dengan pendekatan ini, perumusan FEM berangkat dari definisi integral fungsi. Misalkan kita mempunyai fungsi f ( x ) yang akan diintegralkan dalam selang dari x a ,x b . Dengan metode trapezium kita dapat menghitung integral tertentu ini dalam bentuk jumlahan dari elemen-elemen luas secara keseluruhan, seperti pada gambar (1a). Elemen-elemen luas diaproksimasi dengan sebuah konstanta atau lebih umum sebuah polinomial (Gasiorowicz, 1995).
Untuk problem fisis satu dimensi dalam daerah [ x min , x mak ] maka perhitungan integral aksi dapat dilakukan dalam bentuk diskritisasi region. Daerah integrasi dibagi menjadi elemen-elemen yang lebih kecil yang disebut elemen-elemen hingga, seperti pada gambar (1b). Integral aksi merupakan jumlahan integral-integral aksi itu sendiri pada masing-masing elemen sehingga dapat dituliskan menjadi
    nelem iel iel ) ( ) ( (6) (a) (b)
Gambar 1. (a) Integral f ( x ) merupakan jumlahan elemen-elemen luas di bawah kurva (b) Daerah diskretisasi
Penerapan azas variasi nodal yaitu dengan memvariasikan integral aksi terhadap variabel nodal misalkan ψ α * , maka didapatkan persamaan simultan untuk koefisien ψ β yang belum diketahui             0 * M (7)
Beberapa koefisien ψ β diketahui dari syarat batas (BC) dan dimasukkan ke dalam persamaan simultan ini. Koefisien-koefisien yang lainnya dapat diperoleh dengan metode matriks standar setelah nilai batas diterapkan untuk koefisien yang bersangkutan (Gasiorowicz, 1995).
METODOLOGI PENELITIAN
Penelitian yang dilakukan ini adalah penelitian dalam bidang fisika teori-komputasi. Metode yang digunakan untuk mendapatkan eigen nilai atom hidrogen dalam penelitian ini adalah
x min x max
pendekatan numerik ( numerical approach ) dengan menggunakan metode elemen hingga ( finite element method, FEM ). Dengan menggunakan kesalahan relatif, hasil numerik berupa eigen energi atom hidrogen dibandingkan terhadap hasil analitik.
kesalahan relatif
0 0 100 N T T E E x E  
(8)
E N mewakili tingkat-tingkat energi atom hidrogen secara numerik sedangkan E T mewakili tingkat-tingkat energi atom hidrogen secara analitik. Eigen energi atom hidrogen dikatakan akurat jika kesalahan relatifnya kurang dari 1%. Sedangkan analisis terhadap fungsi gelombang radial dilakukan secara kualitatif melalui visualisasi grafik. Hasil perhitungan numerik fungsi gelombang radial atom hidrogen mempunyai kesesuaian yang baik jika simpangan terhadap hasil analitiknya cukup kecil.

## HASIL DAN PEMBAHASAN
Perumusan Metode Elemen Hingga Untuk Problem Atom Hidrogen Bagian Radial Dengan mengambil bagian radial dari persamaan (2) dengan batasan untuk keadaan ns, maka dapat dituliskan
) ( 1 2 2 2 2 2 r E ψ(r) r e dr d r dr d r m e                        (9)
Untuk tujuan komputasi masalah atomik, persamaan di atas dapat dituliskan kembali dalam satuan yang lebih sesuai yaitu jari-jari Bohr a o sebagai satuan panjang dan Rydberg R ∞ (R ∞ = 13.6 eV) sebagai satuan energi.
) ( ) ( 2 ) ( 1 2 2 r r r r dr d r dr d r             (10)
Prinsip aksi stasioner dapat diterapkan untuk menurunkan bentuk diskret dari persamaan Schrödinger melalui pendekatan metode variasi[2,5,7].
          0 2 ) ( ) 2 )( ( * ) ( ' ) ( *' r r r r r r dr      (11)
Pada daerah asimtotik yaitu di r = ∞ , fungsi gelombang dan turunannya menghilang secara eksponensial, sehingga agar dapat diselesaikan secara numerik, maka batas atas integral aksi harus diubah berhingga, r = r c . Nilai r c dipilih sedemikian rupa sehingga fungsi gelombangnya mendekati nol di daerah asimtotik jauh.
         c r r r r r r r dr
0 2 ) ( ) 2 )( ( * ) ( ' ) ( *'      (12)
Dengan menggunakan prinsip daerah diskritisasi, maka integral aksi di atas dapat diperkirakan dengan memecah daerah integrasi menjadi nelem buah elemen kecil satu dimensi. Untuk daerah fisis [ 0, r c ] seperti gambar berikut ini.
Gambar 2. Elemen di dalam daerah 0, r c
Ukuran setiap elemen yaitu
nelem r h c /  (13)
Gambar 3. Koordinat r dalam elemen i el (Desai, 1988)
Persamaan (12) dievaluasi dalam diskret dan untuk elemen ke- i el dapat dinyatakan

        h r r r r r r dr l e i 0 2 ) ( ) ( ) 2 )( ( * ) ( ' ) ( *'      (14)
Fungsi gelombang ψ ( r ) pada persamaan ini dapat diinterpolasi dengan menggunakan interpolasi polinomial Lagrange. Interpolasi polinomial Lagrange yang digunakan memerlukan dua buah node ( n=2 ) untuk setiap elemen, seperti ditunjukkan pada dibawah ini.
Gambar 4. Fungsi Lagrange untuk interpolasi linear-dua node tiap elemen (Brugioni, 1995) node node node R = (i el –1)h+ηh r h Elemen ke i el η h N 2 (η) N 1 (η) 0 1 η 1 N(η) i el
Fungsi gelombang ψ ( r ) menjadi
                            el i r el i r el i r r el i el i r el i r r el i r l e i r el i N el i r el i N el i r 1 2 1 2 1 2 2 1 ) ( 2 2 ) ( 1 1 ) (      (15)
el i r 1 merupakan titik awal dan el i r 2 merupakan titik akhir. Pada Gambar 4 digunakan koordinat lokal η dalam sebuah elemen yang mempunyai batas [0,1] sehingga titik-titik nodal- nya
h h el i el i r h el i el i r      ) 1 ( 2 ) 1 ( 1 (16)
Subsitusi titik-titik nodal (16) ke dalam fungsi gelombang (15) didapatkan
) ( 2 ) 1 ( 1 ) (      el i el i r    (17)
Dari persamaan ini dapat diketahui bahwa fungsi interpolasi polinomial lokal pada titik awal N 1 ( η ) = 1- η dan pada titik akhir N 2 ( η ) =  η . Turunan dari fungsi gelombang ψ ( r ) adalah
2 1 1 ) (          el i el i h dr r d   
(18)
Faktor h sebagai denominator muncul akibat transformasi koordinat dari r menjadi η . Jika integral aksi pada sebuah elemen (persamaan 14) dituliskan dalam notasi matriks bujur sangkar 2x2, maka
                        1 1 1 1 2 1 ) ( * 2 , ) ( * 1 2 1 0 ) 1 ( 3 ) ( h el i el i el i d h el i                                                         ) ( 2 ) ( 1 2 1 1 2 1 2 1 1 2 1 ) 1 ( 2 el i el i el i h                 (19)
Integral aksi di atas dapat dituliskan dalam bentuk
                                                 ) ( 2 ) ( 1 ) ( 22 ) ( 21 ) ( 12 ) ( 11 ) ( 22 ) ( 21 ) ( 12 ) ( 11 ) ( * 2 , ) ( * 1 ) ( el i el i el i U el i U el i U el i U el i H el i H el i H el i H el i el i el i      (20)
di mana
            n h n n h H n h n n h H H n h n n h H n n n n 3 2 6 1 2 2 3 1 ) ( 22 3 1 6 1 2 2 3 1 ) ( 21 ) ( 12 3 2 2 1 2 2 3 1 ) ( 11                 
,
      2 3 1 6 1 30 1 3 ) ( 22 2 6 1 6 1 20 1 3 ) ( 21 ) ( 12 2 3 1 2 1 5 1 3 ) ( 11 n n h U n n h U U n n h U n n n n           (21) ( n = i el )
Matriks global dapat dibentuk dengan menjumlahkan keseluruhan kontribusi integral aksi dari nelem buah elemen dengan syarat bahwa antar elemen kontinu. Dalam representasi matriks, kontinuitas antar elemen mengijinkan penjumlahan elemen matriks ) ( 11 el i H pada ) 1 ( 22  el i H , hal ini berlaku juga untuk matriks U . Untuk tiga elemen dapat digambarkan
    
            ) 3 ( 22 ) 3 ( 21 ) 3 ( 12 ) 3 ( 11 ) 2 ( 22 ) 2 ( 21 ) 2 ( 12 ) 2 ( 11 ) 1 ( 22 ) 1 ( 21 ) 1 ( 12 ) 1 ( 11 H H H H H H H H H H H H
H ij                  ) 3 ( 22 ) 3 ( 21 ) 3 ( 12 ) 3 ( 11 ) 2 ( 22 ) 2 ( 21 ) 2 ( 12 ) 2 ( 11 ) 1 ( 22 ) 1 ( 21 ) 1 ( 12 ) 1 ( 11 U U U U U U U U U U U U
U ij
Gambar 5. Matriks global H dan U untuk tiga elemen
Dengan menggunakan indeks global sebagai pengganti indeks lokal maka node 1 dan 2 dalam elemen i el dinyatakan kembali sebagai node α = i el dan α +1. Demikian juga nilai nodal fungsi gelombang menjadi ψ α dan ψ α+1 . Total integral aksi untuk seluruh elemen
  , 1 ,..., 2 , 1 , ;      
nelem U H            (22)
Matriks global H dan U memiliki dimensi (nelem+1) x (nelem+1) .
Berdasarkan prinsip aksi stasioner dengan variasi integral aksi terhadap *   maka
             ) ( * * * 0 U H       (23)
Karena variasi *   adalah sembarang maka bentuk umum persamaan nilai eigen
0          U H (24)
Jadi masalah nilai eigen akan dipecahkan sebagai masalah matriks melalui persamaan ini. Variabel yang dicari adalah nilai eigen ε yang merupakan energi atom hidrogen dalam satuan Rydberg R ∞ .
Analisa Perhitungan Tingkat-tingkat Energi Atom Hidrogen Dengan menggunakan kode komputer yang telah dibuat didapatkan nilai eigen yang menyatakan tingkat-tingkat energi atom hidrogen. Pemilihan kombinasi dari parameter-parameret input yang tepat  menghasilkan tingkat-tingkat energi atom hidrogen dengan akurasi yang tinggi. Penulis menggunakan beberapa pilihan kombinasi parameter-parameter input yaitu batas atas integrasi ( r c ) 20-1000 a 0 dan jumlah elemen ( nelem ) 20-1000 elemen.
Gambar 6. Tampilan input - output kode komputer
Potensial Coulomb dalam integral aksi akan mendekati nol untuk r yang besar sehingga mempunyai nilai eigen positif. Untuk atom hidrogen dimana elektron terikat pada inti ( bound state ) maka hanya nilai eigen negatif yang diijinkan. Nilai eigen negatif tersebut merupakan tingkat-tingkat energi atom hidrogen. Tingkat-tingkat energi atom hidrogen secara analitik dapat dirumuskan
   R n n E 2 1
(25)
dimana n adalah bilangan kuantum utama. Hasil perhitungan tingkat-tingkat energi atom hidrogen secara numerik dan secara analitik dalam satuan Rydberg R ∞ serta kesalahan relatifnya dapat diketahui secara langsung dengan kode komputer yang telah dibuat.
Analisa secara terpisah parameter input batas atas integrasi r c dilakukan dengan memberikan input batas atas integrasi 100 a 0 , 200 a 0 , 500 a 0 dan 1000 a 0 dengan jumlah elemen 1000. Dapat dilihat bahwa semakin besar batas atas integrasi maka semakin banyak tingkat-tingkat energi atom hidrogen yang dapat dihitung. Beberapa tingkat energi tertinggi mempunyai kesalahan relatif yang cukup besar. Hal ini disebabkan karena pengaruh syarat batas, dimana pada r = r c fungsi gelombang radial ψ ( r c ) diatur bernilai nol.
Pengaruh jumlah elemen yang digunakan selama perhitungan ( nelem ) terhadap akurasi data yang akan dihasilkan dapat diketahui dengan menganalisis hasil perhitungan tingkat-tingkat energi atom hidrogen secara numerik untuk input r c = 100 a 0 dan nelem = 100, 600 dan 1000. Jika jumlah elemen yang digunakan dalam perhitungan semakin besar maka akurasi hasil perhitungan akan semakin baik. Hal ini dibuktikan dengan kesalahan relatif untuk masing-masing tingkat energi semakin kecil. Sebagai tambahan, jumlah elemen sangat berpengaruh terhadap lamanya running time dari kode komputer yang dibuat. Hal ini disebabkan karena jumlah elemen yang digunakan selama perhitungan menentukan besarnya orde matriks global yang harus ditangani. Semakin besar jumlah elemen maka orde matriks global juga semakin besar, akibatnya running time yang dibutuhkan semakin lama.
Batas atas integrasi dan jumlah elemen secara bersama-sama, dinyatakan dengan ukuran elemen (persamaan 13) juga mempengaruhi akurasi hasil perhitungan. Jika h dalam perhitungan sama, maka akan mempunyai keadaan dasar ( ground state ) yang sama pula. Misalnya untuk h = 1 yaitu pasangan batas atas integrasi r c dan jumlah elemen nelem 20 a 0 : 20, 100 a 0 :100, dan 1000 a 0 : 1000 mempunyai keadaan dasar yang tepat sama sebesar -0.94170913 R ∞ dengan kesalahan relatif 5.82908663%. Hal ini disebabkan karena 20 elemen pertama dalam perhitungan baik untuk 20 elemen, 100 elemen maupun 1000 elemen, mempunyai distribusi elemen yang sama pada seluruh daerah kritis yang berhubungan dengan keadaan dasar. Demikian juga untuk h=1/6 yaitu pasangan batas atas integrasi r c dan jumlah elemen nelem 20 a 0 : 120, dan 100 a 0 : 600  juga mempunyai keadaan dasar yang sama sebesar -0.99772506 R ∞ dengan kesalahan relatif  0.22749426%. Dari dua  contoh nilai h di atas, dapat dilihat bahwa h = 1/6 mempunyai kesalahan relatif  yang lebih kecil dibandingkan h = 1 . Hal ini berarti perhitungan keadaan dasar untuk h = 1/6 mempunyai akurasi yang lebih baik. Berdasarkan nilai kesalahan relatif masing-masing tingkat energi untuk dua nilai h tersebut dan nilai h yang lain yaitu 1/2, 1/5 dan 1/10 dapat disimpulkan bahwa semakin kecil ukuran elemen maka akurasi perhitungan secara numerik juga semakin baik.
Telaah atom hidrogen biasanya hanya melibatkan beberapa tingkat energi pertama yaitu yaitu keadaan dasar dan beberapa keadaan eksitasi ( excitation states ) pertama ( n < 10). Pasangan input batas atas itegrasi r c dan jumlah elemen nelem 100 a 0 : 600, 100 a 0 : 1000, dan  200 a 0 : 1000 menghasilkan eigen energi untuk enam tingkat pertama dengan kesalahan relatif kurang dari 1%. Bahkan untuk pasangan input 200 a 0 : 1000 sampai sembilan tingkat pertama mempunyai kesalahan relatif kurang dari 1%. Jika ditinjau hanya pada lima tingkat energi pertama saja yaitu keadaan dasar dan empat keadaan eksitasi pertama, maka kombinasi parameter input batas atas integrasi r c 100 a 0 dan jumlah elemen sebanyak 1000 merupakan pilihan realistis dan memberikan  cukup informasi mengenai tingkat-tingkat energi atom hidrogen dengan akurasi yang baik dimana  kesalahan relatifnya kurang dari 0,1%.
E 1 =  -0.99917211 R ∞ dengan  kesalahan relatif  0,08278922%, E 2 =  -0.24984445 R ∞ dengan  kesalahan relatif  0,06222133%, E 3 =  -0.11105532 R ∞ dengan  kesalahan relatif  0,05021519%, E 4 =  -0.06247405 R ∞ dengan  kesalahan relatif  0,04151635%, E 5 =  -0.03998598 R ∞ dengan  kesalahan relatif  0,03504057%.
Secara umum, parameter input yang digunakan selama perhitungan dapat diatur sedemikian hingga untuk mendapatkan tingkat-tingkat energi atom hidrogen dengan akurasi yang tinggi. Tentunya, dengan tetap memperhatikan spesifikasi dari komputer yang digunakan.
Analisa Fungsi Gelombang Radial dan Probabilitas Radial Dari perhitungan numerik yang telah dilakukan, juga didapatkan fungsi eigen yang menyatakan fungsi gelombang radial atom hidrogen untuk keadaan ns . Secara analitik, fungsi gelombang radial atom hidrogen untuk tiga keadaan ns yang pertama dalam satuan atomik a.u dapat dirumuskan (Brasden & Joachain, 1995).
) 3 / exp( ) 27 / 2 2 3 / 2 1 ( 2 / 3 ) 3 / 1 ( 2 ) ( 30 ) 2 / exp( ) 2 / 1 ( 2 / 3 ) 2 / 1 ( 2 ) ( 20 ) exp( 2 ) ( 10 r r r r r r r r r             (26)
Visualisasi fungsi gelombang radial atom hidrogen secara numerik dan secara analitik untuk tiga keadaan pertama ( 1s , 2s dan 3s ) untuk pilihan kombinasi parameter input Rc = 100 a 0 dan nelem 100 ditunjukkan pada Gambar 7 sampai Gambar 9.
Kesesuaian antara fungsi gelombang radial atom hidrogen secara numerik dan secara analitik dapat diketahui dengan membandingkan secara langsung melalui visualisasi grafik.
Seperti halnya dalam analisis tingkat-tingkat energi atom hidrogen, batas atas integrasi r c dan jumlah elemen nelem yang digunakan selama perhitungan sangat mempengaruhi akurasi fungsi gelombang radial yang dihasilkan. Fungsi gelombang radial atom hidrogen ) ( r  tidak sama dengan nol pada r = 0 . Begitu juga dengan turunannya 0 / )) ( (  dr r d  pada r = 0. Hal ini disebabkan karena energi potensial untuk atom hidrogen bernilai tak hingga pada titik r = 0. Untuk nilai-nilai n yang besar (keadaan-keadaan eksitasi), fungsi gelombang radial dipengaruhi oleh pemberian syarat batas dimana pada r  =  r c fungsi gelombang radial diatur bernilai nol,
0 ) (   c r r  . Secara umum, untuk n yang besar  tingkat kesesuaian antara fungsi gelombang radial atom hidrogen secara numerik dan secara analitik lebih kecil. Misalnya untuk batas integrasi 20 a 0 , terlihat bahwa ) ( 30 r  mempunyai simpangan yang lebih besar jila dibandingkan dengan ) ( 20 r  dan ) ( 20 r  juga mempunyai simpangan yang lebih besar bila dibandingkan dengan ) ( 10 r  . Berdasarkan visualisasi fungsi gelombang radial atom hidrogen untuk beberapa pilihan kombinasi parameter input, semakin besar batas atas integrasi yang diberikan maka semakin baik kesesuaian antara fungsi gelombang radial atom hidrogen secara numerik dan secara analitik. Hal ini terutama terlihat untuk keadaan- keadaan eksitasi. Jumlah elemen yang digunakan juga mempengaruhi akurasi fungsi gelombang radial atom hidrogen yang dihasilkan dimana semakin banyak jumlah elemen yang digunakan maka semakin baik kesesuaian antara fungsi gelombang radial atom hidrogen secara numerik dan secara analitik. Hal ini berlaku baik untuk keadaan dasar maupun keadaan-keadaan eksitasi.
Kita telah menghitung fungsi gelombang radial atom hidrogen ) ( r  dimana kuadrat dari fungsi gelombang radial 2 ) ( r  menyatakan kerapatan elektron sebagai fungsi dari koordinat radial r . Probabilitas per satuan panjang bahwa elektron dapat ditemukan pada jarak r dari inti disebut fungsi distribusi radial (Brasden & Joachain, 1995).
2 ) ( 2 ) ( r ns r r ns
D  
(27)
Gambar 7.   Fungsi Gelombang Radial Atom Hidrogen untuk n = 1 Batas Integrasi 100 dan Jumlah Elemen 100
Gambar 8.   Fungsi Gelombang Radial Atom Hidrogen untuk n = 2 Batas Integrasi 100 dan Jumlah Elemen 100
Gambar 9.   Fungsi Gelombang Radial Atom Hidrogen untuk n = 3 Batas Integrasi 100 dan Jumlah Elemen 100
Dari visualisasi fungsi distribusi radial dapat diketahui bahwa probabilitas terbesar untuk menemukan elektron pada jarak tertentu r sebanding dengan n 2 . Nilai maksimum fungsi distribusi radial semakin jauh dari inti jika n semakin besar.
Analisis kesesuaian antara fungsi distribusi radial secara numerik dan analitik tidak dibahas dengan alasan fungsi distribusi radial merupakan kuadrat dari fungsi gelombang radial. Artinya jika fungsi gelombang radial secara numerik mempunyai kesesuaian yang baik dengan hasil secara analitik maka fungsi distribusi radial secara numerik juga mempunyai kesesuaian yang baik dengan hasil analitik.

## KESIMPULAN
Telaah problem atom hidrogen secara numerik dengan menggunakan FEM menghasilkan tingkat- tingkat energi dan fungsi gelombang radial atom hidrogen. Hasil perhitungan yang didapatkan sangat bergantung pada parameter input yang digunakan. Batas integrasi dari integral aksi persamaan Schrödinger dan jumlah elemen, baik secara terpisah maupun secara bersama sama,  sangat mempengaruhi hasil perhitungan numerik. Pilihan kombinasi batas atas integrasi 100 a 0 dengan jumlah elemen sebanyak 1000 memberikan tingkat-tingkat energi atom hidrogen dengan akurasi yang baik. Untuk lima tingkat energi pertama yaitu E 1 = -0.99917211 R ∞ , E 2 = -0.24984445 R ∞ , E 3 = - 0.11105532 R ∞ , E 4 = -0.06247405 R ∞ dan E 5 = -0.03998598 R ∞ dimana R ∞ = 13.6 eV. Masing- masing tingkat energi tersebut mempunyai kesalahan relatif  kurang dari 0.1 % terhadap hasil analitik.

## REFERENSI
Brasden, B. H. & Joachain, C. J. (1995). Physics of atoms and molecules . New York: John Wiley & Sons, Inc. Brugioni, J. P. (1995) . Finite element method in quantum mechanics. http://www.science.gmu.edu/~jbrugion/RTF/qm1.htm . Diambil 19 Oktober 2004 . Gasiorowicz, S. (1995). Quantum physics . USA: John Wiley & Sons, Inc. Mohan, L. R. (2002). Finite element and boundary element applications in quantum mechanics. New York: Oxford University Press.
Nikishkov, G. P. (2004). Introduction to the finite elemen method . Lecture Note. University of Aizau, Japan.
Thankappan, V. K. (1985). Quantum mechanics . New Delhi: Wiley Eastern Limited.