DINAMIKA DAN INTERAKSI SOLITON DNA MODEL PEYRARD-BISHOP-DAUXOIS
SWITENIA WANA PUTRI
DEPARTEMEN FISIKA FAKULTAS MATEMATIKA DAN ILMU PENGETAHUAN ALAM
INSTITUT PERTANIAN BOGOR BOGOR 2012
Switenia Wana Putri. Dinamika dan Interaksi Soliton DNA Model PeyrardBishop-Dauxois. Dibimbing oleh : Dr. Husin Alatas.
Abstrak Soliton DNA model PBD menggambarkan dinamika peregangan (denaturasi) pada ikatan hidrogen antar nukleotida dalam rantai DNA yang berbeda yang direpresentasikan oleh potensial Morse. Dengan menggunakan metode beda hingga (finite-difference) dibantu dengan interpolasi Lagrange, maka diperoleh solusi numerik dari dinamika dan interaksi soliton DNA model PBD dengan memberikan gangguan pada solusi stabilnya. Hasil yang diperoleh menunjukkan bahwa solusi numerik dari dinamika dan interaksi soliton DNA model PBD yang diberikan efek gangguan akan menunjukkan hasil yang berbeda dengan solusi eksak tanpa gangguan. Ada tiga gangguan yang diberikan pada solusi stabil ini. Gangguan pertama dan kedua berupa gangguan pada amplitudo menunjukkan bahwa soliton DNA yang pada awalnya stabil mengalami peristiwa undulasi, dimana terjadi penyempitan soliton dan kenaikan amplitudo. Pada gangguan ketiga diberikan gangguan untuk tiga beda fase yang berbeda, dimana akan terbentuk dua soliton yang mengalami superposisi dengan interaksi antar soliton yang berbeda-beda pada setiap beda fasenya. Untuk beda fase ΞΈ=0, kedua soliton akan saling tarik menarik, π lalu kedua soliton akan saling tolak menolak dengan perubahan asimetrik pada ΞΈ= , dan 2 saling tolak menolak dengan perubahan simetrik pada ΞΈ=π. Kata kunci: soliton DNA, PBD, denaturasi, metode beda hingga, interpolasi Lagrange. .
Judul
: Dinamika dan Interaksi Soliton DNA Model Peyrard-Bishop-Dauxois
Nama
: Switenia Wana Putri
NRP
: G74070026
Menyetujui, Pembimbing
Dr. Husin Alatas NIP : 19710604 199802 1 001
Mengetahui, Ketua Departemen Fisika FMIPA IPB
Dr. Akhiruddin Maddu NIP : 19660907 199802 1 000
Disetujui tanggal :
DINAMIKA DAN INTERAKSI SOLITON DNA MODEL PEYRARD-BISHOP-DAUXOIS
Switenia Wana Putri G74070026
Skripsi Sebagai salah satu syarat untuk memperoleh gelar Sarjana Sains pada Fakultas Matematika dan Ilmu Pengetahuan Alam Institut Pertanian Bogor
DEPARTEMEN FISIKA FAKULTAS MATEMATIKA DAN ILMU PENGETAHUAN ALAM
INSTITUT PERTANIAN BOGOR BOGOR 2012
i
KATA PENGANTAR Puji syukur penulis panjatkan ke hadirat Allah swt. karena berkat rahmat, kasih sayang, dan karunia yang tak terhingga dari-Nya penulis dapat menyelesaikan penulisan laporan tugas akhir yang berjudul βDinamika dan Interaksi Soliton DNA Model PeyrardBishop-Dauxoisβ. Laporan ini disusun sebagai salah satu syarat untuk memperoleh gelar Sarjana Sains pada Fakultas Matematika dan Ilmu Pengetahuan Alam Institut Pertanian Bogor. Dalam skripsi ini tertuang penjelasan tentang hasil penelitian tugas akhir yang telah dilakukan oleh penulis dalam beberapa bulan ke belakang. Penulis juga ingin menyampaikan ucapan terima kasih kepada pihak-pihak yang telah membantu dalam penyelesaian penulisan usulan penelitian ini, yaitu kepada : 1. Bapak Dr. Husin Alatas selaku dosen pembimbing yang telah memberikan bimbingan, masukan, dan semangat kepada penulis. 2. Ibu Mersi Kurniati, M.Si dan Bapak Dr. Ir. Irmansyah selaku dosen penguji serta staf dosen Departemen Fisika FMIPA IPB lainnya yang telah memberikan bimbingan, masukan, dan ilmu-ilmu yang bermanfaat selama penulis menempuh pendidikan S1 ini. 3. Mama dan papa selaku orang tua, adik Erin dan Melly, serta keluarga besar di Palembang yang selalu memberi dorongan baik secara materi maupun spiritual. 4. Izzatu Yazidah dan Dede Hermanudin sebagai partners satu topik penelitian ini, terima kasih untuk bantuan dan dukungan yang luar biasa dari kalian. 5. Keluarga di rumah kedua βArsida 4β, Erna Piantari, Mery Purnamasarie, Hesti Paramita Sari, dan Rithoh Yahya. Thanks for being so lovely and take care of me here. 6. Teman-teman seperjuangan di Fisika 44. 7. Keluarga besar Departemen Fisika FMIPA IPB. 8. Teman-teman Ikamusi IPB, teman-teman dari TK sampai SMA, dan teman-teman IPB lainnya yang tidak dapat disebutkan satu per satu, terima kasih untuk doa dan semangatnya. Penulis menyadari bahwa skripsi ini masih jauh dari sempurna. Oleh karena itu, saran dan kritik yang bersifat membangun sangat diperlukan bagi penulis. Semoga tulisan ini bermanfaat bagi semuanya.
Bogor,
Juli 2012
Penulis
ii
DAFTAR RIWAYAT HIDUP Penulis dilahirkan di Palembang pada tanggal 16 Desember 1989 dari pasangan Isman Prabujaya dan Ilya Rosdiana. Penulis menyelesaikan pendidikan Taman Kanak-Kanak hingga Sekolah Menengah Atas di Sumatera Selatan, yaitu TK Tri Daqa Kartika II Sriwijaya Karang Endah, SDN 1 Karang Endah, SMP Negeri 2 Gelumbang, dan SMA Negeri 1 Prabumulih. Pada tahun 2007, penulis diterima menjadi mahasiswa IPB di Departemen Fisika Fakultas Matematika dan Ilmu Pengetahuan Alam melalui jalur Undangan Seleksi Masuk IPB (USMI). Selain aktif mengikuti kegiatan akademik, penulis juga aktif di beberapa organisasi kemahasiswaan di IPB seperti di Ikatan Keluarga Mahasiswa Bumi Sriwijaya, Himpunan Mahasiswa Fisika, Komunitas Seni Budaya Masyarakat Roempoet, dan dalam kepanitiaan pada beberapa acara yang diselenggarakan oleh mahasiswa IPB. Penulis juga berkesempatan menjadi asisten praktikum mata kuliah Fisika Dasar TPB IPB dan mata kuliah Eksperimen Fisika II di Departemen Fisika IPB.
iii
DAFTAR ISI Halaman DAFTAR GAMBAR .............................................................................................
iv
DAFTAR LAMPIRAN .........................................................................................
vi
BAB 1 PENDAHULUAN .....................................................................................
1
1.1 Latar Belakang ....................................................................................
1
1.2 Tujuan Penelitian ................................................................................
1
1.3 Perumusan Masalah ............................................................................
2
1.4 Hipotesis .............................................................................................
2
BAB 2 TINJAUAN PUSTAKA ............................................................................
2
2.1 Deoxyribose Nucleic Acid (DNA) ......................................................
2
2.2 Solusi Soliton DNA Model PBD Teraproksimasi ..............................
3
BAB 3 METODE PENELITIAN ..........................................................................
7
3.1 Tempat dan Waktu Penelitian .............................................................
7
3.2 Peralatan ..............................................................................................
7
3.3 Metode Penelitian ...............................................................................
7
3.3.1 Studi Pustaka ..............................................................................
8
3.3.2 Pembuatan Program Simulasi ....................................................
8
3.3.2.1 Metode Beda Hingga (Finite-Difference) ......................
8
3.3.2.2 Interpolasi Lagrange .......................................................
9
3.3.3 Skala Ulang Parameter ...............................................................
11
BAB 4 HASIL DAN PEMBAHASAN .................................................................
12
4.1 Simulasi Perambatan Soliton akibat Gangguan pada Amplitudo .......
13
4.2 Simulasi Interaksi Dua Buah Soliton ..................................................
15
BAB 5 KESIMPULAN .........................................................................................
18
SARAN β SARAN ................................................................................................
19
DAFTAR PUSTAKA ............................................................................................
19
LAMPIRAN ..........................................................................................................
21
iv
DAFTAR GAMBAR Halaman Gambar 1. (a) Struktur untai komplementer DNA yang membentuk DNA beruntai ganda ...........................................................................
2
(b) Struktur molekul DNA. Atom karbon berwarna hitam, oksigen merah, nitrogen biru, fosfor hijau, dan hidrogen putih .............
2
Gambar 2. Proses replikasi DNA .....................................................................
3
Gambar 3
Solusi soliton pada DNA model PBD.............................................
4
Gambar 4. Representasi grafis model pegas sederhana untuk rantai DN..........
4
Gambar 5. Kurva suatu fungsi f(x) yang dibagi sama besar berjarak h ...........
8
Gambar 6. Karekteristik solusi persamaan NLS soliton DNA model PBD hingga orde-3 stabil ........................................................................
12
(a) profil soliton DNA dalam tiga dimensi ...................................
12
(b) plot hubungan yn (pm) terhadap nl (pm), dimana grafik berwarna merah menunjukkan grafik pada saat Tawal, dan grafik biru menunjukkan grafik pada saat Takhir ......................
12
Gambar 7. Karekteristik solusi persamaan NLS soliton DNA model PBD hingga orde-3 Gangguan I ..............................................................
13
(a) profil soliton DNA dalam tiga dimensi ...................................
13
(b) plot hubungan yn (pm) terhadap nl (pm), dimana grafik berwarna merah menunjukkan grafik pada saat Tawal, grafik biru menunjukkan grafik pada saat Takhir, dan grafik hijau menunjukkan saat terjadinya undulasi ....................................
13
Gambar 8. Karekteristik solusi persamaan NLS soliton DNA model PBD hingga orde-3 Gangguan II .............................................................
14
(a) profil soliton DNA dalam tiga dimensi ...................................
14
(b) plot hubungan yn (pm) terhadap nl (pm), dimana grafik berwarna merah menunjukkan grafik pada saat Tawal, dan grafik biru menunjukkan grafik pada saat Takhir ......................
14
Gambar 9. Karekteristik solusi persamaan NLS soliton DNA model PBD hingga orde-3 Gangguan III ...........................................................
15
(a) profil soliton DNA dalam tiga dimensi untuk ΞΈ = 0 ...............
15
v
(b) plot hubungan yn (pm) terhadap nl (pm), dimana grafik berwarna merah menunjukkan grafik pada saat Tawal, dan grafik biru menunjukkan grafik pada saat Takhir ...................... Gambar 10.
15
Karekteristik solusi persamaan NLS soliton DNA model PBD hingga orde-3 Gangguan III .................................................. (a) profil soliton DNA dalam tiga dimensi untuk ΞΈ =
π 2
...............
16 16
(b) plot hubungan yn (pm) terhadap nl (pm), dimana grafik berwarna merah menunjukkan grafik pada saat Tawal, dan grafik biru menunjukkan grafik pada saat Takhir .................. Gambar 11.
16
Karekteristik solusi persamaan NLS soliton DNA model PBD hingga orde-3 Gangguan III .........................................................
17
(a) profil soliton DNA dalam tiga dimensi untuk ΞΈ = Ο ...............
17
(b) plot hubungan yn (pm) terhadap nl (pm), dimana grafik berwarna merah menunjukkan grafik pada saat Tawal, dan grafik biru menunjukkan grafik pada saat Takhir .....................
17
vi
DAFTAR LAMPIRAN Halaman Lampiran A. Diagram Alir Penelitian .................................................................
22
Lampiran B. SkalaUlang Parameter Persamaan NLS Tahap 1 ...........................
23
Lampiran C. Skala Ulang Parameter Persamaan NLS Tahap 2 ..........................
23
Lampiran D. Source Code Program dengan software Matlab.............................
25
BAB 1 PENDAHULUAN 1.1 Latar Belakang Dinamika transkripsi Deoxyribo Nucleic Acid (DNA) adalah salah satu gejala fisis yang paling menarik dalam biofisika modern karena terkait dengan dasar kehidupan.1 Pada tahun-tahun terakhir gagasan bahwa eksitasi nonlinier memainkan peran dalam dinamika DNA telah menjadi semakin populer, sehingga selanjutnya semakin banyak penelitianpenelitian terkait dengan dinamika soliton DNA ini sendiri, seperti penelitian Englander et al2, Yomosa3, Takeno dan Homma4, Peyrard-Bishop1, PeyrardBishop-Dauxois5, Zhang6 dan lain-lain. Ini merupakan suatu hal yang sangat menarik, menelisik bagaimana βcetak biruβ ini dipandang dari sudut pandang seorang fisikawan, tidak hanya di bidang Biologi atau Biokimia. Banyak penelitian Fisika berkaitan dengan gerak internal amplitudo besar DNA sampai pada kesimpulan bahwa molekul dapat ditinjau sebagai sistem dinamika nonlinier yang terkait dengan gejala soliton.1-6 Struktur heliks ganda DNA mengalami dinamika yang sangat kompleks. Pengetahuan tentang dinamika tersebut dapat memberikan wawasan yang dalam mengenai berbagai fenomena biologis terkait seperti transkripsi, translasi dan mutasi. Diyakini dengan pemahaman akan mekanisme organisme elementer (DNA), kelak diharapkan fenomena makroskopis organisme hidup bisa dijelaskan dan diprediksi dengan akurat dan mudah. Permasalahan utama dalam biofisika sendiri adalah bagaimana menghubungkan sifat fungsional DNA dengan struktur dan karakteristik fisik dinamik-nya. Dari sudut pandang fisikawan, molekul DNA tak lain adalah sistem yang terdiri dari banyak atom berinteraksi, tersusun dengan cara tertentu dalam ruang.7 Sebuah hirarki dari model yang paling penting bagi dinamika nonlinier DNA disajikan oleh Yakushevich.8
Salah satu model dinamika DNA telah diajukan oleh Peyrard dan Bishop1 yang kemudian dikembangkan oleh Dauxois.9,10 Selanjutnya dalam laporan tugas akhir ini, model tersebut disingkat menjadi model PBD. Model PBD terutama ditujukan untuk menggambarkan proses peregangan pasangan basa DNA (denaturasi). Berbeda dengan dua paper Peyrard-Bishop11 sebelumnya, dalam model PBD juga diperhitungkan efek helikoidal yang menggambarkan interaksi antara dua rantai DNA akibat struktur heliks. Soliton DNA model PBD menggambarkan dinamika peregangan pada ikatan hidrogen antar nukleotida dalam rantai yang berbeda yang direpresentasikan oleh potensial Morse, yang terlokalisasi dan bersifat stabil. Dalam model pendekatan potensial Morse hingga orde 4 dikendalikan oleh persamaan Nonlinear Schrodinger 12 (NLS). Profil sebaran rapat energinya menyerupai gundukan yang terpusat dalam rentang ruang berhingga. Setiap soliton dicirikan oleh sifat tidak berubahnya topologi yang menunjukkan sifat kestabilannya. Pada penelitian ini, dilakukan analisis numerik yaitu membuat simulasi model menggunakan pemrograman MATLAB dari solusi analitik persamaan NLS yang telah diperoleh pada penelitian Hermanudin12 sebelumnya. Untuk menentukan solusi numerik dari persamaan NLS ini, metode yang digunakan adalah metode beda hingga (finite- difference) dibantu dengan interpolasi Lagrange.
1.2 Tujuan Penelitian Penelitian ini bertujuan untuk menentukan solusi numerik dinamika dan interaksi soliton DNA model PBD yang diberi efek gangguan dengan menggunakan metode beda hingga (finitedifference) dan interpolasi Lagrange, serta menjelaskan pengertian fisis yang terjadi dari solusi numerik yang diperoleh. Dalam perhitungannya digunakan software MATLAB.
1.3 Perumusan Masalah Pada penelitian sebelumnya, telah diperoleh solusi analitik untuk dinamika DNA model PBD hingga orde ke-4.12 Lalu ada penelitian ini akan ditinjau pemecahan solusi numerik persamaan soliton DNA model PBD dengan memberikan efek gangguan menggunakan metode beda hingga untuk dinamika DNA model PBD hingga orde ke-3.
1.4 Hipotesis Solusi numerik yang akan diperoleh dari dinamika dan interaksi DNA model PBD yang diberikan efek gangguan akan menunjukkan hasil yang berbeda dengan solusi eksak tanpa gangguan.
BAB 2 TINJAUAN PUSTAKA 2.1 Deoxyribose Nucleic Acid (DNA) DNA atau Deoxyribo Nucleic Acid merupakan substansi penurunan sifat yang erat kaitannya dengan hampir semua aktivitas biologi. DNA adalah asam nukleotida, biasanya dalam bentuk heliks ganda yang mengandung instruksi genetik yang menentukan perkembangan biologis dari seluruh bentuk kehidupan sel. DNA berbentuk polimer panjang nukleotida, mengkode barisan residu asam amino
(a) Gambar 1.
13
dalam protein dengan menggunakan kode genetik, sebuah kode nukleotida triplet.13 DNA pada pengertian sempit adalah suatu cetakan untuk membuat protein.14 DNA membawa informasi genetik sel dan terdiri dari ribuan gen. Setiap gen berfungsi sebagai cetakan untuk membangun sebuah molekul protein. Protein melakukan tugastugas penting untuk fungsi sel atau berfungsi sebagai blok bangunan. Aliran informasi dari gen ini menentukan komposisi protein dan fungsi dari sel itu sendiri.15 DNA tersusun sebagai untai komplementer dengan ikatan hidrogen di antara mereka. Masing-masing untai DNA adalah rantai kimia βbatu bata penyusunβ, yakni nukleotida, yang terdiri dari empat tipe: Adenine (A), Cytosine (C), Guanine (G) dan Thymine (T).13 DNA mengandung informasi genetika yang diwariskan oleh keturunan dari suatu organisme; informasi ini ditentukan oleh barisan pasangan basa. Sebuah untai DNA mengandung gen, sebagai βcetak biruβ organisme. Rantai DNA memiliki lebar 22-24 Γ
, sementara panjang satu unit nukleotida 3,3 Γ
.13 Struktur komplementer dan molekul DNA ditunjukkan oleh Gambar 1. Replikasi DNA adalah proses penggandaan rantai ganda DNA. Proses replikasi DNA ini ditunjukkan oleh Gambar 2.
(b)
(a) Struktur untai komplementer DNA yang membentuk DNA beruntai ganda. (b) Struktur molekul DNA. Atom karbon berwarna hitam, oksigen merah, nitrogen biru, fosfor hijau, dan hidrogen putih.
3
Replikasi DNA
DNA Induk
Tahap 1 Pembukaan rantai DNA
Tahap 2 Melengkapi pasangan basa
Dua molekul DNA semikonservatif baru
Tahap 3 Pembentukan gula fosfat yang menghubungkan nukleotidanukleotida pada dua rantai baru
Gambar 2. Proses Replikasi DNA.16 Satu pengamatan penting yang menuntun Watson dan Crick ke penemuan terkenal dari struktur heliks ganda DNA adalah basa cenderung berpasangan melalui ikatan hidrogen dan pasangan yang terbentuk oleh purin dan pirimidin memiliki ukuran yang hampir sama sehingga pasangan-pasangan demikian membentuk struktur dua rantai nukleotida berikatan hidrogen yang sangat teratur. Barisan pasangan basa dalam DNA mengkode informasi untuk mensintesa protein, adalah polimer yang terdiri dari 20 asam amino berbeda. Ikatan hidrogen pada pasangan basa dalam DNA ini pada suatu saat dapat mengalami pemutusan yang disebut dengan proses denaturasi. Selama replikasi DNA, yang terjadi selama pembelahan sel, misalnya molekul disalin dengan cara membukanya seperti ritsleting.Transkripsi DNA adalah pembacaan gen tunggal untuk mensintesa protein.5 Pada sel, replikasi DNA terjadi sebelum pembelahan sel. Sel prokariot terus-menerus melakukan replikasi DNA, sedangkan sel eukariot waktu terjadinya replikasi DNA sangat teratur, yaitu pada fase S daur sel, sebelum mitosis atau meiosis I. Penggandaan tersebut memanfaatkan enzim DNA polimerase yang membantu pembentukan ikatan antara nukleotidanukleotida penyusun polimer DNA. Proses replikasi DNA dapat pula dilakukan in vitro dalam proses yang disebut reaksi berantai polimerase (PCR). Replikasi DNA dikatakan semikonservatif
jika masing-masing strand DNA berlaku sebagai template untuk sintesa rantai baru, dua molekul DNA baru akan dihasilkan, masing-masing dengan satu rantai baru dan satu rantai lama. Proses ini disebut replikasi semikonservatif. 17
2.2 Solusi Soliton DNA Model PBD Teraproksimasi Untuk mendeskripsikan gerak internal DNA dimulai dengan memilih model aproksimasi yang bersesuaian. Biasanya pilihan model bergantung pada masalah yang ditinjau dan akurasi deskripsi yang diperlukan.18 Dalam hal ini untuk gejala denaturasi DNA digunakan model PBD. Sedang untuk mempelajari dinamika model nonlinier DNA model PBD tersebut, diperlukan suatu algoritma yang mencakup beberapa unsur, yaitu memilih gerakan yang dominan, mengkonstruksikan persamaan diferensial parsial terkait dan mencari solusi untuk persamaan, serta menginterpretasikan solusi yang diperoleh.8 Gejala denaturasi pada dinamika DNA ditunjukkan oleh adanya solusi soliton. Representasi grafis dari solusi persamaan gerak yang memiliki solusi soliton ditunjukkan oleh Gambar 3 yang ditafsirkan sebagai keadaan peregangan untai ganda DNA. Dalam dinamika replikasi DNA, DNA diasumsikan mengalami tiga gerakan, yaitu gerak longitudinal (searah sumbu x), gerak bending (searah sumbu y), dan gerak memutar (torsional).
4
un
un-1
s1
m
k
yn (pm)
un+1 Morse Potensial
k
s2
Bentuk B-DNA dalam model Watson-Crick merupakan heliks ganda, terdiri atas dua alur s1 dan s2 (Gambar 4) yang dihubungkan oleh interaksi secara harmonik dengan tetangga terdekat di sepanjang rantai. Untai dipasangkan satu sama lain oleh suatu ikatan hidrogen, yang bersesuaian untuk perpindahan transversal nukleotida.5 Berdasarkan Gambar 4, diasumsikan bahwa massa nukleotida pada DNA sebagai m dan kopling k konstan di sepanjang untai yang sama.5 Struktur helikoid dari untai DNA menunjukkan bahwa nukleotida dari untai yang berbeda menjadi cukup dekat, sehingga nukleotida-nukleotida tersebut dapat berinteraksi melalui filamen air, artinya bahwa nukleotida n salah satu rantai berinteraksi dengan (n+h) dan (n-h) nukleotida lainya. Nukleotida selanjutnya mengalami perpindahan gerak transversal un dan vn dari posisi kesetimbangannya di 4 sepanjang arah ikatan hidrogen, Hamiltonian untuk rantai DNA menjadi 9,10 π
π»=
2
π’2 + π£ 2 +
π 2
2
π· π βπ
π’ π βπ£π
β1
2
Gambar 4. Representasi grafis model pegas sederhana untuk rantai DNA.5
Dengan membuat transformasi ke koordinat pusat massa yang mewakili gerak tranversal ke dalam dan keluar, maka akan lebih mudah menggambarkan gerakan dua alur, yaitu π₯π =
π’ π +π£π
, π¦π =
2
π’ π βπ£π 2
................ (2)
Bentuk Persamaan (1) setelah di transformasi ke koordinat pusat massa Persamaan (2) : π
π» π₯π =
2 πΎ
π₯π2 +
π 2
2 π₯π β π₯πβ1
2 βπ₯π + π₯π+β
2
2 βπ₯π + π₯πββ π
π» π¦π = πΎ 2
2
π¦π2 +
π 2
π¦π + π¦π+β
π· π βπ
2π¦π
2
2
............. (3)
β1
2
+ π¦π + π¦πββ
2
+
+
2 π¦π β π¦πβ1 2
2
+ 2
+
..................... (4)
dari hubungan Persamaan (3) dan (4) akan diperoleh persamaan dinamika DNA yang menggambarkan gelombang linear dan non linear ππ₯π = π π₯π+1 + π₯π β1 β 2π₯π + πΎ π₯π +β + π₯πββ β 2π₯π ....... (5)
π’π β
π’πβ1 2 + π£π β π£πβ1 2 + πΎ π’π β π£π+β 2 + π’π β π£πββ
vn+1
vn
vn-1
nl (pm)
Gambar 3. Solusi soliton pada DNA model PBD
m
2
+
................... (1)
dimana k adalah konstanta harmonik helicoid untuk untai yang sama; K adalah konstanta harmonik helicoid untuk untai yang berbeda; H adalah Hamiltonian potensial morse yang mendekati potensial ikatan hidrogen. D adalah kedalaman potensial morse; dan a adalah jarak antar nukleotida pada rantai yang berbeda.
dan ππ¦ = π π¦π +1 + π¦π β1 β2π¦π β πΎ π¦π +β +π¦π ββ + 2π¦π 2 2ππ· π βπ 2π¦π β 1 π βπ
2π¦π
(6)
Seperti dijelaskan dengan asumsi DNA bergerak dengan amplitudo sangat kecil dalam beberapa artikel dapat menerapkan transformasi π¦ = πΞ¦n atau Ξ¦π =
π¦π πΊ
; (π βͺ 1)...... (7)
5
Substitusi Persamaan (6) ke Persamaan (7) (yang artinya diasumsikan bahwa osilasi nukleotida cukup besar untuk anharmonik, tetapi masih cukup kecil sehingga partikel berosilasi di sekitar bagian bawah potensial Morse) menjadi Ξ¦π =
π π πΎ π
Ξ¦π+1 + Ξ¦πβ1 β2Ξ¦π β Ξ¦π+β + Ξ¦πββ + 2Ξ¦π β
ππ2 (Ξ¦π dimana: ππ2 =
2
2
,πΌ=β
3π
3
+ ππΌΞ¦π + π π½Ξ¦π ) ... (8) 4π 2 π· π
2
,π½=
7π 2 3
(9)
agar Persamaan (8) dapat dipecahkan, harus diterapkan pendekatan semi-diskrit yang artinya bahwa mencari solusi dalam bentuk gelombang soliton DNA model PBD. Ξ¦π π‘ = πΉ1 πππ, ππ‘ π ππ π + π πΉ0 πππ, ππ‘ + πΉ2 πππ, ππ‘ π π2ππ π π3ππ +π. π + π π 3 ................... (10)
Dari Persamaan (10) - (14) diperoleh keadaan untuk semua persamaan dalam Persamaan (8). Kasus kontinu dari Persamaan (8) adalah menyamakan koefisien untuk berbagai keadaan. Sebagai contoh, dengan menyamakan koefisien untuk π ππ π kemudian diperoleh hubungan dispersi π2 =
2π
2πΎ
1 β cos ππ + cos πβπ + π 1+ππ2 .................................... (15) π
Selanjutnya, dengan menurunkan Persamaan (15) secara implisit terhadap q ππ ( ) maka akan diperoleh kecepatan grup ππ
ππ =
π
π sin(ππ) β πΎ sin(πβπ) ... (16)
ππ
Dengan cara yang sama yaitu menyamakan koefisien untuk π π0 = 1 yang artinya cos ππ = cos πβπ = 1, diperoleh πΉ0 =
dengan π β‘ ππ = πππ β ππ‘, .............. (11) dimana l adalah jarak antara dua nukleotida tetangga pada untai yang sama; Ο adalah frekuensi optik dari getaran pendekatan linier; q adalah bilangan gelombang soliton DNA; dan c.c adalah istilah conjugate-complex dari fungsi F0.5 Untuk memecahkan Persamaan (10), fungsi Fi akan ubah ke bentuk kontinu. dengan mengambil batas ππ β π§ dan menerapkan transformasi π = ππ§, π = ππ‘ .............. (12) sehingga menghasilkan kontinu berikut
pendekatan
β2πΌ
πΉ1
4πΎ 2 πππ
1+
dengan π = β2πΌ 1 +
4πΎ
2
.......... (17) β1
π π π2
maka diperoleh hubungan πΉ0 = π πΉ1 Asumsikan π πΉ2 =
π2ππ
2
............... (18)
untuk koefisien πΉ2 π π2 πΌ
2π π
cos 2ππ β1
2πΎ β π
cos 2πβπ +1 +4π 2
πΉ12
............................................................ (19) dengan πΏ =
π π2 πΌ 2π π
2πΎ π
cos 2ππ β1 β
cos 2πβπ +1 +4π 2
πΉ π π Β± β π β πΉ π, π Β± πΉπ π, π ππβ + 1 πΉππ (π, π)π 2 π 2 β2 ... (13)
hubungannya menjadi
dimana FZ dan FZZ berarti turunan 5 terhadap Z. Hal ini memungkinkan untuk mendapatkan ekspresi baru untuk fungsi Ξ¦π (π‘) :
Fungsi F0 dan F2 kemudian dapat dinyatakan ke dalam bentuk fungsi F1 dengan menggunakan transformasi koordinat yang baru lagi
Ξ¦π π‘ = πΉ1 π, π π ππ π + π πΉ0 π, π + πΉ2 π, π π π2ππ + π. π = πΉ1 π ππ π + π πΉ0 + πΉ2 π π2ππ + πΉ1β π βππ π + ππΉ2β π βπ2ππ .......... (14)
π = π β ππ π atau π = ππ π + π, π = ππ ............................................................ (21)
2
πΉ2 = πΏπΉ1 2 ................. (20)
Untuk koefisien π ππ π kasus kontinu
6
π 2 πΉ1ππ β 2ππππΉ1π β π2 πΉ1 =
π
{2πΉ1 cos ππ β 1 + 2ππππΉ1π sin ππ + πΎ π 2 π 2 πΉ1ππ cos ππ }π ππ π β {2πΉ1 cos πβπ + 1 + π 2πππβ πΉ1π sin πβπ + π 2 π 2 β2 πΉ1ππ cos(πβπ)}π πππ β ππ2 πΉ1 + π 2 2πΌ π + πΏ + 3π½ πΉ1 1 πΉ1 π ππ π ................ (22) π
Berdasarkan Persamaan (22) persamaan NLS dapat ditulis ππΉ1π +
π2
1
maka
1 π
π cos ππ β πΎ β2 cos(πβπ) β
2π π
π2
2π π
ππ2
.......... (25)
2π
ππ’ 2 ππ
ππΉ1π + ππΉ1ππ + π πΉ1 2 πΉ1 = 0 ............ (26)
βππ’ + π
ππ 2
Misalkan
ππ
+π
π 2 π’ ππ’
......... (27)
ππ 2 ππ
ππ
π
ππ
=
2π/π
2π
π’=
π
......... (32b)
dan
π ππ π
πππ π ππ , maka ruas kanan
ππ’ π’
2π βπ’ 2 π
ππ
β
2π π
......... (33)
π ππ π
Integrasikan ruas kanan Persamaan (33) terhadap variable Ο, diperoleh
+ ππ’ = 0 ..... (28)
ππ’
2
Persamaan (32b) dalam variable Ο dapat ditulis menjadi
3
+ ππ’3
2π
ππ’ =
ππ 2π π
ππ
ππ’
π
ππ’2 β π’4 ... (32a)
2π βπ’ 2 π
π’
Selanjutnya Persamaan (28) dikali dengan ππ’ , diperoleh βππ’
1 π
ππ’
dengan Ο merupakan frekuensi gelombang soliton DNA yang berperan sebagai variabel bebas dan π’(π) merupakan fungsi real. Kemudian Persamaan (27) disubstitusikan ke Persamaan (26) π2π’
=
atau
Untuk menyelesaikan Persamaan (26) harus menggunakan ansatz (tebakan), sehingga dapat diperoleh solusi eksak soliton DNA. Persamaan ansatz dari persamaan NLS tersebut diberikan πΉ1 π, π = π’(π)π
2
Dari hasil ini dapat ditulis kembali Persamaan (31) menjadi
Maka persamaan NLS pada Persamaan (23) dapat ditulis menjadi5
πππ
π
+ π’4 = π ... (31)
ππ
dimana c merupakan sebuah konstanta. Selanjutnya kembali akan diberikan batasan pada solusi yang memiliki kondisi ππ’ β 0 dan u β 0 pada S β Β± β dan ππ mengimplikasikan c = 0.
π cos ππ β πΎ β2 cos(πβπ) β .............................................
(24) π 2 π π=β 2πΌ π + πΏ + 3π½
2
ππ’ 2
βππ’2 + π
............................................................ (23) dengan koefisien dispersi 1
π
+ π’4 = 0 ... (30)
ππ
yang mengindikasikan bahwa
ππ2πΉ1ππβππ22π2πΌπ+πΏ+3π½πΉ11πΉ1=0
π=
ππ’ 2
βππ’2 + π
2 ππ
π ππ π
=
1 2π π
ππ
1 π πππ
β
πππ π π ππ π
Ruas kanan Persamaan (34) variabel u dapat dinyatakan
= 0 ...... (29)
..... (34) dalam
Persamaan (29) dapat dituliskan kembali dalam bentuk 1 2π π
ππ
1 π πππ
β
πππ π π ππ π
=
1 2π π
ππ
2π π
π’ 2π 2π
1β 1β π’
........................ (35)
Sedangkan integrasi ruas kanan untuk Persamaan (32b) diperoleh ππ 2π/π
π
=
2π π
............... (36)
ππ
π’ 2π 2π
1β 1β
π
π
=π
π’
π
(37)
Dari Persamaan (37) diperoleh persamaan untuk u 2π π
2
π’ π =
π π
ππ₯π π
........ (38)
π π
1+ππ₯π 2π
Bentuk solusi Persamaan (38) dapat diubah menjadi bentuk trigonometri berikut π’ π =
2π π
π
sππβ π
...... (39)
π
Substitusi Persamaan (39) ke Persamaan (27) maka diperoleh 2π
πΉ1 π, π =
π
π
sππβ π
π πππ (40)
π
Solusi eksak untuk Ξ¦π (π‘) dalam ekspresi F1 berdasarkan Persamaan (14), (18) dan (20) adalah Ξ¦π π‘ = πΉ1 π ππ π + π π πΉ1 2 + πΏπΉ12 π π2ππ + π. π .................. (41) Substitusi Persamaan (38) ke Persamaan (39) 2π
Ξ¦π π‘ =
π
sππβ π
2πππ
π
sππβ π
π
π
π ππ +ππ
+
2
π
sππβ π
π 2πππΏ
π
+
π 2
π
π 2π(ππ +ππ ) +
π
π. π ....................................... (42a) Ξ¦π π‘ =
2π π
sππβ π
ππ
+
2πππΏ π
2πππΏ π
π
2cos ππ +
π
sππβ π
sππβ π
2π π
sππβ π
π
π
2
π
2
π
2 cos 2(ππ + ππ )
........ (42b)
π π
,
maka substitusi Persamaan (42b) ke Persamaan (7) diperoleh
Dengan demikian pesamaan yang harus dipecahkan adalah 2π
dengan F =
12
π¦π = ππΉ 2cos ππ + ππ + π 2 πΉ 2 π + 2πΏ cos 2(ππ + ππ ) (43) Penurunan untuk solusi analitik ini diperoleh dan dapat dilihat secara lengkap pada penelitian Hermanudin12. Dalam perhitungan pada penelitian ini juga, digunakan nilai variabel-variabel karakteristik DNA, yaitu sebagai berikut k = 3K = 24 N/m, l = 3,4 x 10β10 m, m = 5,1 x 10β25 kg, a = 2 x 1010 mβ1, D = 0,1 eV. Massa yang digunakan dalam hal ini merupakan massa rata-rata untuk empat buah nukleotida (h = 4).9,10
BAB 3 METODE PENELITIAN 3.1 Tempat dan Waktu Penelitian Penelitian ini dilakukan di Laboratorium Fisika Teori dan Komputasi, Departemen Fisika, Fakultas Matematika dan Ilmu Pengetahuan Alam, Institut Pertanian Bogor dari bulan Februari 2011 sampai dengan Januari 2012.
3.2 Peralatan Peralatan yang digunakan adalah perangkat komputer dengan software yang digunakan MATLAB R2010b dan Microsoft Office 2007. Sebagai pendukung penulis menggunakan sumber literatur, yaitu jurnal-jurnal ilmiah, bukubuku, dan sumber-sumber lain yang terkait.
3.3 Metode Penelitian Metode penelitian ini adalah mencari solusi numerik dari dinamika gangguan soliton DNA model Peyrard-BishopDauxois (PBD) orde ke-3 dengan menggunakan metode finite-difference
8
dibantu dengan interpolasi Lagrange. Langkah-langkah yang dilakukan secara garis besar adalah studi pustaka dan pembuatan program simulasi.
3.3.1 Studi Pustaka Studi pustaka dilakukan dengan membaca dan memahami jurnal-jurnal dan buku-buku yang berkaitan dengan penelitian ini. Studi pustaka juga membantu dalam menganalisis hasil simulasi yang dilakukan.
3.3.2 Pembuatan Program Simulasi Membuat sintak simulasi dengan bantuan software MATLAB menggunakan interpolasi Lagrange dan metode finite-difference. Termasuk dengan mengubah persamaan dinamika DNA model PBD ke dalam persamaan finite-difference dan interpolasi Lagrange. 3.3.2.1 Metode Beda Hingga (Finite Difference) Sebuah aplikasi penting dari persamaan diferensial adalah dalam analisis numerik, terutama dalam persamaan diferensial numerik, tujuannya untuk menentukan solusi numerik dari persamaan diferensial biasa dan parsial. Idenya adalah dengan menggantikan turunan yang muncul dalam persamaan diferensial dengan persamaan finitedifference yang teraproksimasi. Metode yang dihasilkan kemudian disebut metode beda hingga.19 Atau dalam kata lain, metode beda hingga (finite-difference) diterapkan untuk persamaan diferensial dengan melibatkan pergantian semua turunan dengan formula perbedaan (difference).20 Metode ini digunakan untuk membantu dalam menentukan solusi numerik dari Persamaan NLS (26) sebelumnya [halaman 11]. Persamaan tersebut kemudian dapat diubah ke dalam bentuk persamaan differensial biasa, menjadi π
ππΉ ππ‘
+π
π2πΉ ππ₯ 2
+ π πΉ 2 πΉ = 0 ..... (44)
Maka persamaan tersebut dapat diselesaikan dengan melakukan pendekatan numerik beda hingga untuk masing-masing turunan parsial. Untuk melakukannya pertama kali pilih angka integer sembarang yaitu N dimana N > 0 dan membagi interval [a, b] dengan (N + 1) sebagaimana diilustrasikan pada Gambar 5, sedemikan rupa sehingga β=
πβπ π+1
.................. (45)
Dengan demikian maka titik-titik x yang merupakan sub-interval antara a dan b dapat dinyatakan sebagai21 π₯π = π + πβ, π = 0,1, β¦ , π + 1 ...... (46) Karena Persamaan (44) dibangun oleh dua parameter π₯ dan π‘, maka selanjutnya parameter π‘ juga dapat dituliskan sebagai π‘π = π + πβ, π = 0,1, β¦ , π + 1 .... (47) dimana β dapat pula dideskripsikan sebagai Ξπ₯ dan Ξπ‘ (β = Ξπ₯ = Ξπ‘) dengan menyesuaikan parameter yang dimaksud. Selanjutnya Persamaan (44) akan diubah ke dalam bentuk formula metode beda hingga. Pencarian solusi persamaan diferensial melalui pendekatan numerik dilakukan dengan memanfaatkan polinomial Taylor yang dapat dituliskan πΉ π₯π + βπ₯ = πΉ π₯π + βπ₯ βπ₯ 2 π 2 πΉ 2 ππ₯ 2
π₯π
+
ππΉ
+
ππ₯ π₯ π βπ₯ 3 π 3 πΉ
3! ππ₯ 3 π₯ π
+β―
................................................................... (48)
Gambar 5. Kurva suatu fungsi f(x) yang dibagi sama besar berjarak h.20
9
sedangkan polinomial Taylor untuk πΉπ+1 dan πΉπβ1 dapat dinyatakan πΉπ+1 = πΉ π₯π + βπ₯ βπ₯ 3 π 3 πΉ
ππΉ ππ₯ π₯ π
+
βπ₯ 2 π 2 πΉ 2 ππ₯ 2 π₯ π
+
Kemudian jika dijumlahkan dengan maka
πΉπ+1 + πΉπβ1 = 2πΉ π₯π + (βπ₯)2
+ β― ....................... (49)
3! ππ₯ 3 π₯ π
2(βπ₯)4 π 4 πΉ
βπ₯ 3 π 3 πΉ
ππΉ ππ₯ π₯ π
+
βπ₯ 2 π 2 πΉ 2 ππ₯ 2 π₯ π
+
+ β― ....................... (50)
3! ππ₯ 3 π₯ π
π2πΉ ππ₯ 2
π₯π
atau
ππΉ
π2πΉ
πΉπ+1 β πΉπβ1 = 2βπ₯
ππ₯ π₯ π 2(βπ₯)3 π 3 πΉ 3!
+
ππ₯ 3 π₯ π
ππ₯ 2 π₯ π
+ β― ...... (51)
Turunan πΉ π₯π terhadap π₯ adalah ππ₯ π₯ π
πΉπ +1 βπΉπ β1
=
2βπ₯
β
(βπ₯)3 π 3 πΉ 3!
ππ₯ 3 π₯ π
+β―
=
=
πΉπ +1 β2πΉπ +πΉπ β1 βπ₯ 2
πΉπ +1 β2πΉπ +πΉπ β1 βπ₯ 2
ππΉ
atau
ππ‘ π‘ π +1 ,π₯ π
ππ₯ π₯ π
β
πΉπ +1 βπΉπ β1 2βπ₯
β
(βπ₯)3 π 3 πΉ 3!
ππ₯ 3 π₯ π
+β―
.......................................................... (52b) Selanjutnya ruas kanan pada Persamaan (52b) dapat dikatakan sebagai error dari Persamaan (51) yang artinya nilainya dapat diabaikan dengan memberikan notasi O. Sehingga Persamaan (51) dapat dituliskan kembali menjadi ππΉ ππ₯ π₯ π
π
=
πΉπ +1 βπΉπ β1 2βπ₯
+ π(βπ₯) ... (53)
πΉ π₯ π ,π‘ π +1 β πΉ(π₯ π ,π‘ π β1 ) 2Ξt
+π
=
+
terhadap π₯
(βπ₯)2 π 4 πΉ 12
ππ₯ 4 π₯ π
(55a)
+ π(βπ₯ 2 ) ... (55b)
πΉππ +1 +πΉππ 2βπ‘
(βπ‘) ... (56)
Sementara aproksimasi untuk turunan π2πΉ
kedua terhadap jarak ( 2 ) pada ππ₯ Persamaan (44) berdasarkan Persamaan (55b) dalam iterasi waktu ke-π dapat dituliskan persamaannya π2πΉ ππ₯ 2
π₯π
=
πΉππβ1 +2πΉππ βπΉππβ1 βπ₯ 2
+ π(βπ₯ 2 ) ... (57)
Kemudian Persamaan (56) dan (57) disubstitusikan ke Persamaan (44) maka
πΉ π₯ π +1 ,π‘ π β 2πΉ π₯ π ,π‘ π + πΉ(π₯ π β1 ,π‘ π ) Ξt
ππ₯ 2 π₯ π
Selanjutnya dilakukan aproksimasi ππΉ turunan waktu ( ) pada Persamaan (44) ππ‘ berdasarkan Persamaan (53). Dimana dapat dituliskan persamaannya
........................................................... (52a) ππΉ
+
π2πΉ
+ β― ...... (54)
Turunan kedua dari πΉ π₯π adalah
Jika Persamaan (49) dikurangi dengan Persamaan (50) maka
ππΉ
ππ₯ 4 π₯ π
4!
πΉπβ1 = πΉ π₯π β βπ₯
Persamaan (49) Persamaan (50),
Ξx 2
+ π πΉ 2πΉ = 0
πΉ π₯π , π‘π +1 β πΉ π₯π , π‘π β1 = 2π 2 π(πΉ π₯π+1 , π‘π β 2πΉ π₯π , π‘π + πΉ π₯π β1 , π‘π ) Ξx + π πΉ(π₯π , π‘π ) 2 πΉ(π₯π , π‘π ) π₯π‘ πΉ π₯π , π‘π +1 = 2π 2 π(πΉ π₯π +1 , π‘π β 2πΉ π₯π , π‘π + πΉ π₯π β1 , π‘π ) + π₯π₯ π πΉ(π₯π , π‘π ) 2 πΉ(π₯π , π‘π ) + πΉ π₯π , π‘π β1 ............ (58) Persamaan (58) inilah yang selanjutnya diubah ke dalam bahasa pemrograman MATLAB.
3.3.2.2 Interpolasi Lagrange Pada metode finite-difference sebelumnya telah diuraikan bahwa untuk menentukan solusi numerik, cara pertama dengan memilih angka integer sembarang yaitu N dimana N > 0 dan membagi
10
interval [a, b]. Sebelum menentukan nilai a dan b, harus diketahui terlebih dahulu nilai titik-titik sebelum a dan nilai titiktitik setelah b, agar mempermudah dalam menentukan nilai a dan b yang dimaksud. Untuk itu digunakanlah interpolasi Lagrange dalam menentukan titik-titik tersebut.
π π₯ = πΏ0 π₯ π¦0 + πΏ1 (π₯)π¦1 ... (60) Substitusi Persamaan (59a) dan (59b) ke Persamaan (60), maka akan didapat
dan πΏ1 π₯ =
π₯βπ₯ 1 π₯ 0 βπ₯ 1
π₯βπ₯ 0 π₯ 1 βπ₯ 0
π₯ 0 βπ₯ 1
π¦0 +
π₯βπ₯ π π¦ π₯ 1 βπ₯ 0 1
... (61)
dan ketika π₯=π₯0 π π₯0 =
Interpolasi Lagrange diterapkan untuk mendapatkan fungsi polinomial P(x) berderajat tertentu yang melewati sejumlah titik data. Misalnya, untuk mendapatkan fungsi polinomial berderajat satu yang melewati dua buah titik yaitu (x0, y0) dan (x1, y1).21 Langkah pertama yang dilakukan adalah mendefinisikan fungsi berikut 10 πΏ0 π₯ =
π₯βπ₯ 1
π π₯ =
π₯ 0 βπ₯ 1 π₯ 0 βπ₯ 1
π¦0 +
π₯ 0 βπ₯ π π₯ 1 βπ₯ 0
π¦1 = π¦0 ... (62a)
dan pada saat π₯=π₯1 π π₯1 =
π₯ 1 βπ₯ 1 π₯ 0 βπ₯ 1
π¦0 +
π₯ 1 βπ₯ π π₯ 1 βπ₯ 0
π¦1 = π¦1 ...(62b)
Dari persamaan tersebut dapat disimpulkan bahwa Persamaan (61) benarbenar melewati titik (π₯π , π¦π ) dan (π₯1 , π¦1 ).22
........... (59a)
Persamaan (61) dinamakan interpolasi Lagrange derajat 1. Nama interpolasi ini diambil dari nama penemunya, yaitu Joseph Louis Lagrange yang berkebangsaan Perancis. Bentuk umum interpolasi Lagrange derajat β€ n untuk (n+1) titik berbeda adalah
........... (59b)
kemudian definisikan fungsi polinomial sebagai berikut π
π¦π πΏπ π₯ = π¦0 πΏ0 π₯ + π¦1 πΏ1 π₯ + β― + π¦π πΏπ β¦ β¦ β¦ β¦ β¦ β¦ β¦ β¦ β¦ β¦ β¦ β¦ . (63)
π π₯ = π=0
dengan i = 0, 1, 2, ...., n dan π π₯ β π₯π π₯ β π₯0 π₯ β π₯1 β¦ π₯ β π₯πβ1 π₯ β π₯π+1 β¦ π₯ β π₯π πΏπ π₯ = = β¦ β¦ (64) π₯π β π₯π (π₯π β π₯) π₯π β π₯π β¦ π₯π β π₯πβ1 π₯π β π₯π+1 β¦ (π₯π β π₯π ) π =0 π β π
Mudah dibuktikan bahwa: πΏπ π₯π =
π₯
1 ,π = π 0 ,π β π
Jika terdapat N data yang terdiri dari titik-titik π₯0 , π₯1 , π₯2 , π₯3 , dan seterusnya,
π π₯ =
π₯βπ₯ 1 (π₯βπ₯ 2 ) π¦ π₯ 0 βπ₯ 1 (π₯ 0 βπ₯ 2 ) 0
+
β2β (β3β)
ββ (β3β)
ββ (β2β)
π¦0 +
π₯0
β
π₯1
β
π₯2
β
π₯3
β
π₯4
π
dan jarak antara titik satu dengan lainnya adalah h, maka Persamaan (63) dapat ditulis untuk tiga titik terdekat (π₯0 , π₯1 , π₯2 , π₯3 )
dan polinom interpolasi P(x) melalui setiap titik data.23
π π₯ =
β
π₯βπ₯ 0 (π₯βπ₯ 2 ) π¦ π₯ 1 βπ₯ 0 (π₯ 1 βπ₯ 2 ) 1
β(ββ)
π¦1 +
+
π₯βπ₯ 0 (π₯βπ₯ 1 ) π¦ π₯ 2 βπ₯ 0 (π₯ 2 βπ₯ 1 ) 2
ββ (β2β) 2β (β)
π¦2
π π₯ = 3π¦0 β 3π¦1 +π¦2 ............................................................................ (65)
11
Persamaan (65) kemudian dibuat dalam bentuk algoritma dengan menggunakan bahasa MATLAB. 3.3.3 Skala Ulang Parameter Dalam pengerjaannya, nilai variabelvariabel yang digunakan pada persamaan NLS ini masih mengalami masalah dalam perhitungan numeriknya, dimana nilai yang digunakan masih terlalu besar untuk cakupan numerik pada software MATLAB sehingga error yang dihasilkan juga menjadi cukup besar. Untuk mengatasi masalah ini, maka dilakukanlah skala ulang terhadap nilai parameter pendukung. Skala ulang parameter ini tidak mengubah persamaan, tetapi dilakukan dengan cara memasukkan variabel pengali lain yang nilainya kecil yang dapat ditentukan sendiri, ke dalam parameter t, x, dan F, sehingga nilai ketiga parameter tersebut menjadi lebih kecil (penurunan lengkap dapat dilihat pada Lampiran B dan C) π‘β
π‘ π
;
π₯β
π₯
π π
π
;
dan
πΉ β ππΉ
π π
............................................................ (66) dimana π dan π merupakan variabel skala ulang yang dimaksud. Nilai kedua variabel ini harus disesuaikan dengan parameter t, x, dan F agar diperoleh hasil yang tepat. Dalam hal ini nilai variabel π dan π berturut-turut digunakan nilai 0,0005 dan 0,03 untuk gangguan pertama dan kedua, serta 0,001 dan 0,03 untuk gangguan ketiga. Nilai ini diperoleh dari hasil trial and error yang dilakukan selama penelitian. Parameter t, x, dan F yang telah diskala ulang, kemudian disubstitusi ke Persamaan (58), sehingga persamaannya menjadi: Ξt πΉ π₯π , π‘π +1 = 2ππ 2 (πΉ π₯π +1 , π‘π β Ξx 2 2πΉ π₯π , π‘π + πΉ π₯πβ1 , π‘π ) + π 2 Ξt π πΉ1 2 πΉ1 + πΉ π₯π , π‘π ................... (67)
Selain persamaan utama, persamaan anzats yang digunakan juga mengalami perubahan. Dengan cara yang sama variabel t, x, dan F yang telah diskala ulang dan disubstitusi ke Persamaan (40) [halaman 7], sehingga diperoleh hasil πΉ=
2 π
sech
π₯ π
............ (68)
Dalam penelitian ini, solusi numerik dari persamaan NLS stabil diberikan tiga kasus gangguan. Hal ini dengan cara memberikan variabel tambahan pada persamaan ansatz, sehingga dapat diamati perubahan dari masing-masing keadaan setelah diberi gangguan. Ketiga persamaan tersebut adalah 2
πΉ π₯π , π‘π =
π
πΉ π₯π , π‘π =
π
2
π₯
π ππβ
1 + π ... (69a)
π
π ππβ
1+π
π₯ π
1+π
.......................................................... (69b) πΉ π₯π , π‘π =
2 π
(π ππβ
π ππβ
π₯ π +π₯ 0
π₯ π βπ₯ 0 π
π
+
π ππ ).......... (69c)
Persamaan (69a), (69b), dan (69c) selanjutnya digunakan untuk menggantikan Persamaan ansatz (68) solusi stabil. Sama seperti Persamaan (68) persamaan-persamaan tersebut digunakan untuk menentukan solusi baru dari persamaan NLS soliton DNA yang digunakan. Dari ketiga persamaan ini akan dihasilkan solusi baru yang berbeda dengan solusi stabilnya. Persamaan-persamaan yang telah diubah ke dalam bentuk persamaan metode beda hingga dan interpolasi Lagrange, serta telah dilakukan skala ulang kemudian disusun dalam bahasa pemrograman MATLAB [Lampiran D]. Dari program simulasi ini akan ditampilkan output berupa grafik tiga dimensi yang merupakan hasil solusi numerik dari soliton DNA model PBD ini. Solusi numerik yang telah diperoleh ini selanjutnya diamati dan dianalisa, sehingga dapat dijelaskan fenomena yang terjadi.
BAB 4 HASIL DAN PEMBAHASAN
Dari Gambar 6 tampak bahwa dari solusi yang diperoleh membentuk profil soliton dengan amplitudo yang cukup stabil yaitu 1,189 pm. Dari hasil ini, seperti pada penelitian Hermanudin13, maka dapat dikatakan persamaan ansatz (tebakan) yang digunakan sudah tepat. Dari gambar juga tampak gambaran umum dari proses replikasi (denaturasi) DNA yang bergerak dominan ke arah un.
yn (pm)
Melalui penerapan metode bedahingga dengan interpolasi Lagrange sebagai syarat batas terkait, maka solusi numerik dari dinamika dan interaksi soliton DNA model PBD dapat dicari dan disajikan dalam bentuk grafik tiga dimensi. Hasil perhitungan dibagi menjadi empat keadaan, yaitu keadaan tanpa gangguan kecil, dengan gangguan kecil, serta interaksi antara dua buah solusi soliton. Untuk solusi numerik pada kondisi stabil (tanpa gangguan), dihasilkan solusi soliton yang memiliki
karakteristik stabil sejak waktu awal (T = 0) hingga waktu akhirnya. Hasil solusi numerik pada kondisi stabil ini ditunjukkan oleh Gambar 6.
T (s)
nl (pm)
yn (pm)
(a)
nl (pm) (b) Gambar 6. Karekteristik solusi persamaan NLS soliton DNA model PBD hingga orde-3 stabil. (a) profil soliton DNA dalam tiga dimensi. (b) plot hubungan yn (pm) terhadap nl (pm), dimana grafik berwarna merah menunjukkan grafik pada saat Tawal, dan grafik biru menunjukkan grafik pada saat Takhir.
13
4.1. Simulasi Perambatan Soliton Akibat Gangguan pada Amplitudo Setelah diperoleh solusi untuk persamaan NLS kubik soliton DNA model PBD yang stabil, selanjutnya, solusi stabil tersebut diberikan gangguan kecil terhadap amplitudo yaitu dengan cara mengalikan Persamaan (68) yang stabil dengan suatu nilai 1+π [lihat Persamaan (69a) pada halaman 11]. Dalam hal ini digunakan nilai Ξ΅ = 0,25.
yn (pm)
Pada kasus pertama ini, tampak terjadi perubahan dari solusi stabil (Gambar 6) sebelumnya setelah diberikan gangguan. Soliton DNA yang terbentuk mengalami undulasi seperti yang ditunjukkan pada Gambar 7(a). Pada saat
undulasi terjadi penyempitan yang diiringi dengan kenaikan amplitudo, yang pada waktu awal amplitudo sebesar 1,495 pm, sedangkan pada puncak undulasinya amplitudo mencapai 2,87 pm. Pada Gambar 7(b) juga ditunjukkan bahwa amplitudo untuk solusi gangguan pertama ini lebih tinggi daripada solusi stabilnya. Hal ini menunjukkan bahwa gangguan yang diberikan pada ansatz akan mempengaruhi amplitudo dari soliton. Dengan memberikan gangguan pada solusi stabilnya, artinya akan mengubah profil soliton sendiri. Selain mengalami perubahan amplitudo, tampak pula pada Gambar 7(b), soliton mengalami dispersi yang lebih besar daripada solusi stabilnya (Gambar 6(b)).
(a)
T (s)
yn (pm)
nl (pm)
nl (pm) (b) Gambar 7. Karekteristik solusi persamaan NLS soliton DNA model PBD hingga orde-3 Gangguan I. (a) profil soliton DNA dalam tiga dimensi. (b) plot hubungan yn (pm) terhadap nl (pm), dimana grafik berwarna merah menunjukkan grafik pada saat Tawal, grafik biru menunjukkan grafik pada saat Takhir, dan grafik hijau menunjukkan saat terjadinya undulasi.
14
Artinya, gangguan yang diberikan juga mempengaruhi hubungan dispersi pada persamaan Hamiltoniannya.
yn (pm)
Selanjutnya, gangguan kecil kembali diberikan pada persamaan ansatz untuk solusi stabil dengan cara mengubah Persamaan (68) dengan Persamaan (69b) [lihat halaman 11]. Untuk gangguan kedua ini digunakan Ξ΅ = 0,5. Profil yang ditunjukkan pada solusi gangguan kedua ini hampir sama dengan gangguan kecil sebelumnya, tetapi seperti yang ditunjukkan pada Gambar 8(a), undulasi yang terjadi tampak lebih lebar dibandingkan dengan gangguan Ξ΅ = 0,25, namun amplitudo undulasi yang terbentuk
lebih kecil, yaitu 2,437 pm. Dari Gambar 8(b) juga terlihat bahwa amplitudo soliton ini lebih tinggi daripada solusi stabilnya. Dispersi yang terjadi pada kasus dengan gangguan Ξ΅ = 0,5 lebih besar dan sama seperti kasus Ξ΅ = 0,25 dimana undulasi pada kasus dengan Ξ΅ = 0,5 tersebut juga mengakibatkan pengurangan jumlah nukleotida dalam proses denaturasi, hanya saja pada kasus ini, nukleotida yang berkurang lebih sedikit daripada kasus dengan Ξ΅ = 0,25. Dari hasil simulasi gangguan dengan dua nilai Ξ΅ berbeda tersebut tampak bahwa pada keduanya terjadi peristiwa undulasi.
T (s) (a)
yn (pm)
nl (pm)
nl (pm) (b) Gambar 8. Karekteristik solusi persamaan NLS soliton DNA model PBD hingga orde-3 Gangguan II. (a) profil soliton DNA dalam tiga dimensi. (b) plot hubungan yn (pm) terhadap nl (pm), dimana grafik berwarna merah menunjukkan grafik pada saat Tawal, dan grafik biru menunjukkan grafik pada saat Takhir.
15
Ketika terjadi penyempitan pada soliton atau dalam kata lain terjadi undulasi, efek nonlinier mengalami ketidakstabilan yang lebih dominan daripada efek dispersinya. Setiap terjadi undulasi menunjukkan terjadinya pengurangan jumlah eksitasi nukleotida yang terlibat dalam proses denaturasi, dimana nukleotida yang pada awalnya dapat meregang terhalangi oleh efek nonlinier ini.
4.2. Simulasi Interaksi Dua Buah Soliton
yn (pm)
Pada simulasi ini, ditinjau propagasi dua soliton yang memiliki jarak x diantara
keduanya. Persamaan modifikasi yang digunakan adalah Persamaan (69c), [lihat halaman 11] dengan membuat variasi pada π nilai ΞΈ, yaitu ΞΈ = 0, ΞΈ= , dan ΞΈ = Ο. Hasil 2 simulasi untuk ΞΈ = 0 ditunjukkan oleh Gambar 9. Gambar 9(a) menunjukkan bahwa untuk kasus dengan ΞΈ = 0 ini, akan terbentuk dua buah soliton yang kemudian mengalami superposisi dengan memberikan kondisi x = 0,03 pm dan ΞΈ = 0. Artinya diberikan jarak antara soliton satu dengan yang lain adalah 0,03 pm dengan beda fase 0. Pada gangguan ini tampak pada gambar, kedua soliton yang awalnya terpisah dengan jarak 0,03 pm
T (s) nl (pm)
yn (pm)
(a)
nl (pm) (b) Gambar 9. Karekteristik solusi persamaan NLS soliton DNA model PBD hingga orde-3 Gangguan III. (a) profil soliton DNA dalam tiga dimensi untuk ΞΈ = 0. (b) plot hubungan yn (pm) terhadap nl (pm), dimana grafik berwarna merah menunjukkan grafik pada saat Tawal, dan grafik biru menunjukkan grafik pada saat Takhir.
16
bahwa kondisi ini menggambarkan keadaan tarik-menarik antar kedua soliton. Perubahan ini tampak jelas pada Gambar 9(b), dimana pada grafik berwarna merah yang menunjukkan kondisi pada T = 0, masih terdapat dua buah profil soliton yang terpisah, tetapi kemudian pada grafik biru yang menunjukkan kondisi pada T akhir, hanya terdapat satu profil soliton. Dua soliton yang terbentuk menunjukkan bahwa nukleotida yang terlibat dalam proses denaturasi terlokalisasi dalam dua ruang. Artinya, ketika terjadi superposisi, nukleotida-nukleotida yang terlokalisasi dalam dua ruang menyatu kembali.
yn (pm)
menjalar dengan bentuk dan kecepatan yang sama, namun pada T = 1,37 x 10-10 s kedua soliton semakin mendekat dan keduanya berangsur-angsur bergabung menjadi satu soliton (mengalami superposisi). Sebelum terjadi superposisi, masing-masing soliton memiliki amplitudo sebesar 0,6594 pm. Amplitudo terus menurun hingga 0.44 pm pada T = 1,37 x 10-10 s, kemudian seiring bertambahnya waktu terjadi superposisi dan amplitudo berangsur-angsur kembali meningkat seperti amplitudo awalnya. Jika iterasi diteruskan (T diperpanjang) kedua gelombang terpisah kembali dengan bentuk dan kecepatan yang sama dengan sebelum terjadinya superposisi. Jelas
T (s)
nl (pm)
yn (pm)
(a)
nl (pm) (b) Gambar 10. Karekteristik solusi persamaan NLS soliton DNA model PBD hingga orde-3 Gangguan III. π (a) profil soliton DNA dalam tiga dimensi untuk ΞΈ= . 2
(b) plot hubungan yn (pm) terhadap nl (pm), dimana grafik berwarna merah menunjukkan grafik pada saat Tawal, dan grafik biru menunjukkan grafik pada saat Takhir.
17
π
kedua soliton saling tolak menolak dan mengalami perubahan yang asimetrik. Dapat terlihat jelas juga dari grafik hubungan yn (pm) dan nl (pm) pada Gambar 10(b), awalnya profil amplitudo dua soliton sama (grafik merah), namun kemudian di T akhir, kedua soliton memiliki amplitudo yang berbeda dari sebelumnya, selain itu terjadi pula dispersi yang cukup signifikan dari kondisi awalnya (grafik berwarna biru). Untuk ΞΈ = Ο, profil soliton yang terbentuk ditunjukkan oleh Gambar 11. Dari Gambar 11(a), tampak dua soliton yang terbentuk awalnya memiliki profil yang sama dengan amplitudo yang semakin menurun hingga saat T = 1,35 x 10-10 s.
yn (pm)
Sementara itu untuk kondisi ΞΈ= , hasil 2 simulasi ditunjukkan oleh Gambar 10. Pada Gambar 10(a) tampak bahwa, pada awalnya terbentuk dua buah soliton yang terpisah dengan jarak 0,03 pm dengan amplitudo dan kecepatan yang sama, namun kemudian amplitudo berangsurangsur mengecil seiring pertambahan waktu dari keadaan awal amplitudo 0,6237 pm, kemudian pada saat T = 1,35 x 10-10 s hingga waktu akhir amplitudo salah satu soliton kembali meningkat, namun pada soliton yang lain, amplitudo mengalami penurunan. Artinya, terjadi perbedaan amplitudo dari kedua soliton, yang satu meningkat sementara yang lain menurun. Hal ini menunjukkan
T (s) (a)
yn (pm)
nl (pm)
nl (pm) (b) Gambar 11. Karekteristik solusi persamaan NLS soliton DNA model PBD hingga orde-3 Gangguan III. (a) profil soliton DNA dalam tiga dimensi untuk ΞΈ = Ο. (b) plot hubungan yn (pm) terhadap nl (pm), dimana grafik berwarna merah menunjukkan grafik pada saat Tawal, dan grafik biru menunjukkan grafik pada saat Takhir.
Kemudian amplitudo soliton kembali meningkat (meskipun tidak signifikan) seiring dengan pertambahan waktu hingga iterasi akhir. Jarak antara kedua soliton tersebut juga semakin melebar dari kondisi awalnya (Gambar 11(b)). Pada waktu awal, amplitudo masing-masing soliton adalah 0,5642 pm, sedangkan di akhir amplitudo menurun menjadi 0,4524 pm. Berbeda dengan kondisi sebelumnya π (ΞΈ = ), masing-masing soliton mengalami 2 perubahan yang sama, atau dapat dikatakan kedua soliton saling tolak menolak tetapi dengan perubahan yang simetrik. Dari ketiga kasus interaksi dua soliton yang diamati, tampak bahwa dinamika soliton DNA model PBD bergantung pada fase awal yang cukup memberikan dampak yang signifikan pada dinamikanya. Terlihat bahwa soliton mengalami undulasi, kenaikan dan penurunan amplitudo, serta dispersi yang semakin meningkat.
BAB 5 KESIMPULAN Pada penelitian sebelumnya, telah diperoleh solusi analitik untuk persamaan NLS kubik DNA model PBD.12 Selanjutya, berdasarkan solusi analitik yang telah diperoleh, dikembangkan kembali penelitian untuk mencari solusi numerik persamaan NLS kubik DNA model PBD dengan menggunakan metode beda hingga dan interpolasi Lagrange. Dari hasil yang diperoleh menunjukkan bahwa solusi antara solusi analitik dan solusi numerik sudah sesuai. DNA model PBD yang digunakan pada penelitian ini menunjukkan perilaku dinamika DNA pada saat proses denaturasi replikasi DNA. Hasil dari solusi numerik yang dilakukan diperolehlah profil soliton dinamika terkait. DNA mengalami dinamika yang cukup stabil dari proses awal hingga akhirnya. Hasil ini juga menunjukkan bahwa persamaan ansatz yang diperoleh dari solusi analitik sudah tepat.
Dalam penelitian ini, persamaan ansatz yang diperoleh dari solusi analtik persamaan NLS kubik DNA model PBD, diberikan gangguan. Dilakukan sebanyak tiga gangguan pada solusi stabil persamaan NLS DNA model PBD ini. Masing-masing hasil solusi numerik dari beberapa kasus gangguan yang ditinjau menunjukkan hasil yang berbeda-beda. Untuk gangguan yang relatif kecil diperoleh hasil bahwa terjadi peristiwa undulasi pada soliton yang awalnya stabil. Undulasi ini mengakibatkan soliton semakin menyempit dan mengalami kenaikan amplitudo. Hasil ini menunjukkan bahwa pada saat undulasi, terjadi pengurangan jumlah nukleotida yang terlibat dalam proses denaturasi pada rantai DNA. Selain itu, gangguan ini juga menyebabkan efek pada koefisien dispersinya, sehingga tampak terjadi dispersi pada solusi numeriknya. Untuk gangguan yang lebih besar, diperoleh hasil yang hampir sama dengan kasus gangguan sebelumnya. Perbedaannya, pada gangguan yang lebih besar undulasi yang terjadi lebih lebar dan amplitudo puncak undulasinya lebih rendah jika dibandingkan dengan gangguan yang lebih kecil. Artinya, pengurangan jumlah nukleotida yang terlibat dalam proses denaturasi pada rantai DNA lebih sedikit. Selain itu, dispersi yang terjadi tampak lebih besar untuk gangguan yang lebih besar. Peristiwa undulasi yang terjadi pada kedua gangguan tersebut ini menunjukkan bahwa efek nonlinier mengalami ketidakstabilan yang lebih besar daripada efek dispersinya, yang menyebabkan soliton mengalami penyempitan. Efek dispersi lebih dominan ketidakstabilannya pada saat soliton tidak mengalami undulasi (keadaan sebaliknya). Berbeda dengan kasus gangguan kecil, pada interaksi dua soliton dengan jarak x, pada proses denaturasinya, nukleotida terlokalisasi dalam dua ruang. Telah ditinjau kasus tiga kondisi fase awal yaitu dengan membuat variasi pada beda fase antara kedua soliton, pada saat ΞΈ = 0, π , dan Ο. Pada saat ΞΈ = 0, kedua soliton 2 saling tarik menarik, pada awalnya dua
soliton bergerak dengan amplitudo yang sama yang nilainya terus menurun, kemudian saat T = 1,37 x 10-10 s, kedua soliton mengalami superposisi, sehingga terbentuklah satu soliton. Soliton ini kemudian bergerak dengan mengalami peningkatan amplitudo hingga waktu akhir. Dinamika ini akan berlanjut, dimana soliton akan kembali terpisah pada T tertentu kemudian mengalami superposisi kembali, begitu pula π seterusnya. Untuk beda fase ΞΈ = , kedua 2 soliton saling tolak menolak dan mengalami perubahan yang asimetrik. Pada selang waktu T = 0 s hingga T = 1,35 x 10-10 s soliton bergerak dengan amplitudo yang sama yang nilainya terus menurun. Kemudian pada T > 1,35 x 10-10 s hingga akhir iterasi waktu, jarak antar kedua soliton semakin melebar, lalu salah satu soliton bergerak dengan mengalami peningkatan amplitudo, sedangkan soliton yang lain terus bergerak dengan amplitudo yang semakin menurun dari kondisi awalnya. Dan untuk beda fase ΞΈ = Ο, kedua soliton tolak menolak tetapi dengan perubahan yang simetrik. Pada selang waktu T = 0 s hingga T = 1,35 x 10-10 s dua soliton bergerak dengan amplitudo yang sama yang nilainya terus menurun, selanjutnya pada saat T > 1,35 x 10-10 s hingga akhir iterasi waktu, jarak antar kedua soliton semakin melebar, kemudian amplitudo kedua soliton kembali meningkat (walaupun tidak terlalu signifikan) dan terus bergerak stabil hingga T akhir. Dari hasil yang diperpoleh ini, dapat disimpulkan pula bahwa hasil solusi numerik yang diperoleh dari dinamika DNA model PBD yang diberikan efek gangguan akan menunjukkan hasil yang berbeda dengan solusi eksak tanpa gangguan.
SARAN β SARAN Penelitian tentang dinamika permodelan DNA merupakan salah satu hal yang menarik untuk dilakukan. Untuk pengembangan selanjutnya, diharapkan dapat melakukan beberapa hal,
diantaranya dengan menggunakan metode numerik yang memiliki tingkat akurasi yang lebih baik lagi, sehingga dapat dilakukan perhitungan dengan lebih cepat dan tepat. Selain itu, dapat pula dilakukan (baik numerik maupun analitik) ekspansi potensial morse dengan orde yang lebih tinggi. Dengan mengambil ekspansi hingga orde yang lebih tinggi diharapkan akan diperoleh hasil yang lebih akurat seperti keadaan sebenarnya. Pada penelitian ini, permodelan yang digunakan masih bersifat ideal, dalam artian mengabaikan pengaruh-pengaruh luar yang bersifat mengganggu. Selanjutnya mungkin dapat dilakukan penelitian dengan memperhitungkan pengaruhpengaruh luar tersebut. Dalam permodelan DNA sendiri ada banyak sekali model yang terus dikembangkan, yang mungkin dapat diteliti selain daripada model PBD ini.
DAFTAR PUSTAKA 1. Peyrard, M. dan Bishop, A.R. (1989). Statistical mechanics of a nonlinear model for DNA denaturation. Phys. Rev. Lett. 62, 2755-2758. 2. Englander, S.W. et al. (1980). Nature of the open state in long polynucleotide double helices: possibility of soliton excitations. Proc. Natl. Acad. Sci. USA 77, 7222-7226. 3. Yomosa, S. (1984). Solitary excitations in deoxyribonucleic acid (DNA) double helices. Phys. Rev. A 30, 474480. 4. Homma, S. dan Takeno, S. (1984). A coupled base-rotator model for structure and dynamics of DNA, Prog. Theor. Phys. 72, 679-693. 5. ZdravkoviΔ, S., TuszyΕski, J.A. dan SatariΔ, M.V. (2005). Peyrard-BishopDauxois model of DNA dynamics and impact of viscosity. Journal of Computational and Theoretical Nanoscience Vol.2: 1β9. 6. Zhang, Ch.T. (1989). Harmonic and subharmonic resonances of microwave
20
absorption in DNA, Phys. Rev. A-40, 2148-2153. 7. Hadi, M. dan Nurlita, I. βSoliton dan DNAβ. Fisikanet. 2006. Web. 17 Februari 2011.
8. Yakushevich, L.V. (1998). Nonlinear Physics of DNA. Wiley Series in Nonlinear Science, John Wiley, Chichester, 181. 9. Dauxois, T. (1991). Dynamics of breathers modes in a nonlinear βhelicoidalβ model of DNA. Phys. Lett. A-159, 390-395. 10. Dauxois, T. dan Peyrard, M. (1991). Dynamics of breather modes in a nonlinear helicoidal model of DNA. Lecture Notes in Physics 393, Dijon,p.79.
18. Ludmila, V. dan Yakushevich. (2001). Is DNA a nonlinear dynamical system where solitary conformational waves are possible?. J. Biosci., Vol.26, No.3: 305-313. 19. Anonim. βFinite Differenceβ. Wikipedia. 2011. Web. 17 Februari 2011. 20. Recktenwald, G.W. βFiniteDifference Approximations to The Heat Equationβ. 2004. Web. 17 Februari 2011. <www.f.kth.se/~jjalap/numme/FDhea t.pdf> 21. Suparno, S. (2008). Komputasi untuk Sains dan Teknik βMenggunakan Matlab-. Departemen Fisika FMIPA, Univeristas Indonesia,181-182.
11. ZdravkoviΔ, S. dan SatariΔ, M.V. (2001). Impact of viscosity on DNA dynamics. Phys. Scripta 64, 612-615.
22. Supriyanto. (2007). Komputasi untuk Sains dan Teknik. Universitas Indonesia. Jakarta.
12. Hermanudin, D. (2011). Efek osilasi anharmonik pada soliton Deoxyribo Nucleic Acid Peyrard-BishopDauxois. Skripsi. Departemen FisikaFMIPA IPB.
23. Munir, R. (2006). Metode Numerik. Informatika. Bandung
13. Anonim. βDNAβ. Wikipedia. 2011. Web. 02 Maret 2011. 14. Hamdani, S. βDNA sebagai Kode Genetikβ. 2012. Web. 03 Mei 2012. 15. Anonim. βDNA-RNA-Proteinβ. 2012. Web. 03 Mei 2012. 16. Anonim. βDNA vs. RNA and replicationβ.2011. Web. 18 Februari 2011. 17. Anonim. βReplikasi DNAβ. Wikipedia. 2010. Web. 18 Februari 2011.
LAMPIRAN
22
Lampiran A. Diagram Alir Penelitian Mulai
Penelusuran literatur
Membuat sintak simulasi MATLAB dengan menggunakan metode beda hingga (finite-difference)
Menganalisis hasil solusi numerik
Mengolah dan menyusun hasil penelitian
23
Lampiran B. Skala Ulang Parameter Persamaan NLS Tahap 1 Persamaan NLS awal π
ππΉ ππ‘
+π
π2πΉ ππ₯ 2
+ π πΉ1 2 πΉ1 = 0 .................................... (26)
dengan persamaan ansatz (tebakan) πΉ =
2π π
π
sππβ π
π
π πππ .......................................... (40)
misalkan: ππ‘ = π ;
π₯=π₯
π π
;
πΉ=πΉ
π π
................................ (B.1)
Substitusi permisalan pada persamaan (64) ke (26) ππ ππ
ππΉ ππ ππΉ ππ
ππ π
ππΉ ππ
+π +π
π ππΉ π ππ
+
π π2πΉ π ππ₯ 2 π2πΉ ππ₯ 2
+π
π2πΉ ππ₯ 2
+ π πΉ 2 πΉ = 0 ..................................................................................... (B.2)
+ π πΉ 2 πΉ = 0 ........................................................................................ (B.3) π π2 πΉ π
ππ₯ 2
+π
π
π
π
π
πΉ
2
πΉ = 0 .................................................................. (B.4)
2
+ πΉ πΉ = 0 ................................................................................................ (B.5)
Substitusi permisalan pada persamaan (B.1) ke (40) π
πΉ=πΉ
πΉ=
π
...................................................................................................................... (B.1)
π
2π
π
π
sech π₯
π
π
π
π
π ππ ................................................................................. (B.6)
πΉ = 2 sech π₯ π ππ ..................................................................................................... (B.7)
Lampiran C. Skala Ulang Parameter Persamaan NLS Tahap 2 Skala ulang tahap 2 dilakukan, karena pada skala ulang tahap 1, skala yang diperoleh masih terlalu besar, dan belum compatible dengan software MATLAB yang digunakan. Hal dilakukan dengan cara membuat skala ulang kembali pada variabel x dan F.
24
Misalkan: π₯=
π₯ π
πΉ = ππΉ ............................................. (C.1)
;
Substitusi persamaan (C.1) ke (B.5), sehingga ππ π
ππΉ ππ‘
ππΉ ππ‘
+ ππ
+π
π2πΉ ππ₯ 2
π2πΉ ππ₯ 2
+ π 3 πΉ 2 πΉ = 0 .................................................................................. (C.2)
+ π 2 πΉ 2 πΉ = 0 ........................................................................................ (C.3)
Substitusi persamaan (C.1) ke (B.7) πΉ= πΉ=
πΉ π 2 π
sech
π₯ π
π ππ ................................................................................................... (C.4)
Jika dilakukan substitusi persamaan (C.1) ke persamaan (B.1), maka variabel t, x, dan F dapat disederhanakan:
ππ‘ = π
; π₯β
π₯
π π
π
;
πΉ β ππΉ
π π
.......................... (66)
25
Lampiran D. Source Code Program dengan software Matlab clear all clc % M N % % %
>> ..:Definisi Kondisi Awal dan Inkremen:.. << = input ('M = '); = input ('N = '); x01 = input ('x01 = '); teta = input ('teta = '); x01 dan teta digunakan pada saat menjalankan simulasi gangguan ketiga; epsilon = input ('epsilon = '); % gunakan epsilon=0.25 untuk gangguan kedua; % gunakan epsilon=0.5 untuk gangguan ketiga; dx = 0.01; dt = 2.5000e-005; x0 = -N/2*dx; % >> ..:Definisi Nilai Parameter Pembangun:.. << D = 0.1*1.6e5; mass=5.1e-13; a=3e-2; k=24; K=8; l=340; lamda=10*l; q=2*pi/lamda; h=4; sigma=1e10; tau=0; e2=0.001; % >> ..:Definisi hubungan dispersi:.. << wg =2*a*sqrt(D/mass); w=sqrt(wg+(2*k*(1-cos(q*l))+2*K*(cos(q*l*h)+1))/mass); % >> ..:Definisi kecepatan group:.. << Vg=(l*(k*sin(l*q)-K*h*sin(l*q*h))/(mass*w)); % >> ..:Definisi koefisien orde-1 sampai orde-3:.. << alfa =(-3*a)/sqrt(2); beta =7*a^2/3; % >> ..:Definisi hubungan F1 dengan F0 dan F2:.. << miu=-2*alfa/(1+4*K/(mass*wg.^2)); delta=(wg.^2)*alfa/(4*w^2+2*k*(cos(2*q*l-1)/mass2*K*(cos(2*h*q*l)+1)... /mass-wg^2)); % >> ..:Definisi koefisien dispersi(P) dan koefisien nonlinear(Q):.. << P=(((l^2)*(k*cos(l*q)-K*(h^2)*cos(l*q*h))/mass)-Vg^2)/(2*w); Q=(-wg.^2*(2*alfa*(miu+delta)+3*beta)/(2*w)); % >> ..:Parameter rescale:.. << phi = 0.0005; % 0.001 untuk gangguan ketiga <<
26
psi = 0.03; % >> ..:Matriks F dan Y:.. << F=zeros(N,M); Y=zeros(N,M); for n=1:N x(n,1)=x0+(n-1)*dx; F(n,1)=((sqrt(2)/psi)*sech(x(n,1)/sqrt(phi)); %.: untuk gangguan pertama :. % F(n,1)=((sqrt(2)/psi)*sech(x(n,1/sqrt(phi)))*(1+epsilon); %.: untuk gangguan kedua :. % F(n,1)=((sqrt(2)/psi)*sech((1+epsilon)*(x(n,1)/sqrt(phi))* % (1+epsilon); %.: untuk gangguan ketiga :. % F(n,1)=(sqrt(2)/psi)*((sech(((x(n,1)+x01)/sqrt(phi))))+... % ((sech((x(n,1)-x01)/sqrt(phi)))*(exp(1i*teta))); end % >> ..:Running Program:.. << for m=1:M for n=2:N-1 % >> ..:Interpolasi Lagrange:.. << if n==2 F(n-1,m)=3*F(n,m)-3*F(n+1,m)+F(n+2,m); end if n==N-1 F(n+1,m)=3*F(n,m)-3*F(n-1,m)+F(n-2,m); end % >> ..:Metode Finite Difference:.. << if m==1 F(n,m+1)=1i*((phi^2*r*(F(n+1,m)-2*(F(n,m))+F(n-1,m))+... (psi^2*dt*(conj(F(n,m))*(F(n,m)^2))+F(n,m); else F(n,m+1)=2*1i*((phi^2*r*(F(n+1,m)-2*(F(n,m))+F(n1,m)))+... (psi^2*dt*(conj(F(n,m))*(F(n,m)^2)))+F(n,m-1); end t(m)=(m-1)*dt; theta(n,m)=(n*q*l)-(w*t(m)/e2); Y(n,m)= ((e2*(psi*F(n,m)/(sqrt(Q/sigma))))*(2*(cos((sigma*tau)+... theta(n,m)))))+((e2^2*((psi*F(n,m)/(sqrt(Q/sigma)))^2))... *(miu+(2*delta*(cos(2*((sigma*tau)+theta(n,m))))))); %Y adalah persamaan gelombang untuk soliton DNA% end m end % >> ..:Grafik 3 dimensi soliton DNA:.. << figure surf(t/sigma,x*sqrt(P/sigma)/sqrt(phi),abs(Y(:,1:M))); view(0,90); colorbar shading interp
27
xlabel ('T'); ylabel ('x'); zlabel ('Y'); % >> ..:Grafik Hubungan Y terhadap x:.. << figure plot((x*sqrt(P/sigma)),abs(Y(:,0))); xlabel ('x(pm)'); ylabel ('Y(pm)'); figure plot((x*sqrt(P/sigma)),abs(Y(:,M))); xlabel ('x(pm)'); ylabel ('Y(pm)');
