ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - Nguyễn Thị Trang NGHIÊN CỨU MÔ PHỎNG CỦA SỐ BÀI TOÁN VỀ TĨNH ĐIỆN TRONG SINH HỌC LUẬN VĂN THẠC SĨ KHOA HỌC Hà Nội - 2015 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - Nguyễn Thị Trang NGHIÊN CỨU MÔ PHỎNG CỦA SỐ BÀI TOÁN VỀ TĨNH ĐIỆN TRONG SINH HỌC Chuyên ngành: Vật lí lí thuyết vật lí toán Mã số: 60 44 01 03 LUẬN VĂN THẠC SĨ KHOA HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC: TS NGUYỄN THẾ TOÀN Hà Nội - 2015 LỜI CẢM ƠN Để hoàn thành luận văn "Nghiên cứu mô số toán tĩnh điện sinh học", xin bày tỏ lòng biết ơn sâu sắc tới TS Nguyễn Thế Toàn - Giảng viên Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội, người tận tình hướng dẫn, trang bị kiến thức động viên suốt trình thực luận văn Tôi xin gửi lời cảm ơn chân thành tới thày cô giáo môn Vật lý lý thuyết - Khoa Vật lý - Trường Đại học Khoa học Tự nhiên trang bị kiến thức chuyên môn cần thiết tạo điều kiện thuận lợi để hoàn thành luận văn Xin chân thành cảm ơn quan tâm, giúp đỡ, tạo điều kiện Ban chủ nhiệm khoa Vật Lý, phòng Sau Đại học trường Đại Học Khoa Học Tự Nhiên - Đại Học Quốc Gia Hà Nội Tôi xin gửi lời cảm ơn tới gia đình, quan bạn bè, đồng nghiệp sát cánh, giúp đỡ động viên suốt trình học tập hoàn thành luận văn Cuối cùng, xin gửi lời cảm ơn quỹ NAFOSTED đề tài số 103.02-2012.75 hỗ trợ trình thực đề tài nghiên cứu Tôi xin chân thành cảm ơn! Hà Nội, tháng năm 2015 Học viên thực Nguyễn Thị Trang MỤC LỤC MỞ ĐẦU 1 Lý chọn đề tài Phương pháp nghiên cứu 3 Bố cục luận văn CHƯƠNG TỔNG QUAN VỀ CẤU TRÚC VÀ MỘT SỐ TÍNH CHẤT VẬT LÝ CỦA PHÂN TỬ ADN .5 1.1 Cấu trúc hóa học phân tử ADN 1.2 Cấu trúc không gian phân tử ADN 1.2.1 Cấu trúc sơ cấp 1.2.2 Hiện tượng cuộn xoắn cấu trúc bậc cao ADN 1.3 ADN xem polymer tích điện, tính chất đàn hồi ADN 11 1.3.1 Định nghĩa độ dài quán tính dây polymer 12 1.3.2 Đóng góp hiệu ứng tĩnh điện vào độ dài quán tính ADN 14 CHƯƠNG PHƯƠNG TRÌNH POISSON - BOLTZMANN VÀ LÝ THUYẾT TĨNH ĐIỆN DEBYE - HUCKEL 19 2.1 Phương trình Poisson - Boltzmann 19 2.1.1 Phương trình Poisson - Boltzmann; Phương trình trường trung bình tự hợp 19 2.1.2 Tuyến tính hóa phương trình Poison - Boltzmann (PB); Phương trình Debye Huckel (DH) 21 2.2 Phương trình Poisson - Boltzmann cho quanh hình trụ tích điện….24 2.2.1 Nghiệm phương trình Debye - Huckel cho quanh hình trụ tích điện 25 2.2.2 Hiện tượng cô đọng "Manning" .26 CHƯƠNG PHƯƠNG PHÁP MÔ PHỎNG MONTE CARLO 30 3.1 Phương pháp mô Monte Carlo 30 3.2 Điều kiện cân chi tiết thuật toán Metropolis 33 3.2.1 Điều kiện cân chi tiết .33 3.2.2 Thuật toán Metropolis 34 CHƯƠNG NHỮNG KẾT QUẢ MÔ PHỎNG THU ĐƯỢC CHO DÂY POLYMER TÍCH ĐIỆN 37 4.1 Mô hình dây polymer tích điện 37 4.2 Những kết mô thu 41 4.2.1 Sự phụ thuộc cách tỷ lệ độ dài đầu - cuối vào bán kính chắn 42 4.2.2 Những kết thu giới hạn bán kính chắn lớn 46 4.2.3 Hàm tương quan góc liên kết 51 KẾT LUẬN 54 TÀI LIỆU THAM KHẢO 56 PHỤ LỤC 59 DANH MỤC HÌNH VẼ Hình 1.1 Cấu trúc loại base phân tử ADN…………………………………6 Hình 1.2 Cấu trúc hóa học phân tử ADN……………………………………7 Hình 1.3 Cấu trúc xoắn kép ADN………………………………………… Hình 1.4 Một số dạng cấu trúc hình học ADN………………………………9 Hình 1.5 Giản đồ sơ lược bước xoắn ADN hạt nhân tế bào…….10 Hình 1.6 Cấu trúc nucleosome…………………………………………………11 Hình 1.7 Mô hình chuỗi "blobs" tĩnh điện…………………………………… 16 Hình 2.1 Các ion linh động phân bố quanh hình trụ tích điện phân tử ADN………………………………………………………………………… 26 Hình 2.2 Phản ion phân bố que cứng tích điện với tham số Manning ξeff…………………………………………………………………………… 29 10 Hình 3.1 Giản đồ phương pháp mô Monte Carlo………………………32 11 Hình 4.1 Dây polymer tích điện mô hình hóa chuỗi gồm N cầu ……………………………………………………………………………… 37 12 Hình 4.2 Phép dịch chuyển pivot thay đổi cấu hình polymer cách lấy phần chuỗi polymer từ monomer chọn ngẫu nhiên đến cuối chuỗi quay theo góc ngẫu nhiên quanh trục ngẫu nhiên………………38 13 Hình 4.3 Phép dịch chuyển flip thay đổi cấu hình polymer cách chọn ngẫu nhiên monomer thực phép quay monomer quanh trục nối hai monomer lân cận…………………………………………………………… 39 14 Hình 4.4 Quá trình đạt trạng thái cân hệ với cấu hình ban đầu chọn dạng khối cầu Gaussian …………………………………………………….40 15 Hình 4.5 Các cấu hình ban đầu mô phỏng: (a) Dạng que, (b) Dạng khối cầu Gaussian……………………………………………………………………….40 16 Hình 4.6 Sự phụ thuộc bình phương độ dài đầu cuối PE (�2 ) vào bán kính �� chắn rs với số monomer khác nhau: N=64, N=128, N=256, N=512, N=2048, N=4096 Mũi tên bên phải hình vẽ cho giá trị �2 ⁄�� mô Coulomb �� không chắn (�� → ∞)………………………………………………………….41 17 Hình 4.7 Đồ thị lý thuyết phụ thuộc � = ���[�2 ]/����� vào �� theo lý �� thuyết OSF với ��~��2 (đường nét liền) ��~�� (đường nét đứt)…………… 44 18 Hình 4.8 Kết mô phụ thuộc � vào �� với giá trị N khác nhau: N=64, N=128, N=256, N=512, N=1024, N=2048, N=4096…………….45 19 Hình 4.9 Mật độ điện tích tuyến tính η thể hàm bán kính chắn với giá trị khác N: N=64, N=128, N=256, N=512, N=1024, N=2048 Khi N tăng mật độ điện tích tuyến tính η ��⁄�� ≫ tiến tới đường liền nét hình vẽ (không phụ thuộc N)…………………………………… 48 20 Hình 4.10 Đồ thị phụ thuộc ��/�� vào �� tính toán theo phương trình (4.6), (4.7), (4.8) sử dụng liệu R2 thu từ mô với mật độ điện ee tích tuyến tính � = ��…………………………………………………………50 21 Hình 4.11 Logarit hàm tương quan góc �(�) phụ thuộc vào khoảng cách x (trong đơn vị ��) dọc theo chuỗi PE cho trường hợp N = 512, rs = 50lB…….52 BẢNG KÝ HIỆU CÁC CHỮ VIẾT TẮT STT Viết tắt ADN BACF DH KK ĐDQT MC OSF PB PE Cụm từ viết tắt Axit deoxyribonucleic Hàm tương quan góc (The bond angle correlation function) Debye - Huckel Khokhlov - Khachaturian Độ dài quán tính Monte Carlo Odijk, Skolnick Fixman Poisson - Bolztmann Polymer tích điện (Polyelectrolyte) MỞ ĐẦU Lý chọn đề tài Trải qua lịch sử nghiên cứu lâu dài, axit deoxyribonucleic (ADN) biết đến phân tử axit nucleic mang thông tin di truyền mã hóa cho hoạt động sinh trưởng phát triển dạng sống bao gồm số virus Trong thể sống ADN tìm thấy nhân tế bào tế bào chất (đối với sinh vật nhân sơ) Tại đó, ADN tham gia vào trình tổng hợp, điều tiết số lượng protein, trình trao đổi chất, sinh trưởng, phát triển hoạt động di truyền qua hệ… Chính ADN đóng vai trò quan trọng hoạt động thể sống hay gọi "phân tử sống" Với tầm quan trọng vậy, nghiên cứu vấn đề liên quan đến ADN trở thành đề tài nhận quan tâm nhiều nhà khoa học hàng đầu nước quốc tế Chúng ta dễ dàng tìm thấy tạp chí khoa học có số ảnh hưởng (IF) cao, tạp chí hàng đầu như: Nature, Science… nhiều công trình nghiên cứu ADN Ngoài vai trò quan trọng ngành khoa học sống (sinh học, công nghệ sinh học, bệnh học, di truyền, …), ADN biết phân tử có cấu trúc trật tự tinh thể chiều, bền (có thể tồn lâu kể sau vật chủ mất) Do vậy, ADN ứng cử viên tiềm cho nghiên cứu liên quan tới khoa học vật liệu công nghệ vật liệu mới, mô vật liệu tự nhiên, vật liệu sinh học…của kỷ 21 Trong dung môi nước phân tử ADN trở nên tích điện âm xem dây polymer tích điện với mật độ điện tích tuyến tính khoảng −1e/1.7A� (e điện tích proton) Và ta coi ADN hình trụ tích điện với bán kính 1nm, mật độ điện tích bề mặt 1e/1n�2 Đây mật độ điện tích lớn quan sát hệ sinh học Do hiệu ứng tĩnh điện đóng vai trò quan trọng cấu trúc hoạt tính hệ ADN Chẳng hạn việc đóng gói ADN hạt nhân tế bào theo cấu trúc đa cấp khác nhau, việc gói ADN virus, việc đọc gene từ ADN… có mặt protein tương tác với ADN dùng hiệu ứng tĩnh điện Trong thí dụ nêu trên, tính chất đàn hồi ADN tham số có đóng góp quan trọng vào lượng tự cấu trúc đóng gói Các nghiên cứu toán đàn hồi ADN cho thấy, tương tác tĩnh điện Coulomb monomer có ảnh hưởng lớn tới độ đàn hồi phân tử Một số kết nghiên cứu cho kết tương tác tĩnh điện làm thay đổi độ đàn hồi phân tử ADN Mặc dù ý tưởng tương tác tĩnh điện làm tăng độ cứng dây polymer chấp nhận mặt định tính kiểm chứng số kết thực nghiệm, song phụ thuộc độ dài tĩnh điện �� vào bán kính chắn �� dung môi đề tài gây nhiều tranh cãi Lý thuyết OSF cho kết độ dài tĩnh điện �� tỉ lệ với bình phương bán kính chắn ��, số lý thuyết tính toán máy tính khác lại cho kết �� tỉ lệ với �� mũ , �� tỉ lệ với �� với số mũ nhỏ Như toán phụ thuộc �� vào �� dù phát biểu rõ ràng trực quan, lý thuyết thống Với mong muốn nghiên cứu ảnh hưởng số hiệu ứng tĩnh điện lên phân tử ADN đặc biệt, hi vọng tìm kết rõ ràng cho toán phụ thuộc độ dài tĩnh điện, khoảng cách đầu cuối dây polymer vào bán kính chắn dung môi Bên cạnh tìm hiểu nguyên nhân dẫn tới không thống kết luận phụ thuộc �� vào �� nêu lý tác giả lựa chọn đề tài "Nghiên cứu mô số toán tĩnh điện sinh học" đề tài cho luận văn Trong luận văn này, sử dụng mô Monte - Carlo máy tính để xem xét phụ thuộc độ dài tĩnh điện le khoảng cách đầu - cuối Ree dây polymer vào bán kính chắn rs Để thực mô phỏng, sử dụng thư viện SimEngine viết TS Nguyễn Thế Toàn Đây thư viện mô vật lý modular dùng ngôn ngữ C++ mở rộng OpenMP, OpenCL để song song hóa tính toán CPU nhiều lõi GPGPU (card đồ họa tính toán) Điều Chúng ta thấy đường ��/�� → �� bão hòa giá trị �2�2⁄72� ��� �� ≫ �, nhiên theo phương trình (4.9) việc giảm xảy �� lớn (��/� ≈ 0.25) Độ lệch khỏi ��~��2 �� lớn thể hình (4.10) điều kiện �� ≤ � bị vi phạm � giảm nhẹ theo �� 4.2.3 Hàm tương quan góc liên kết Một phương pháp chuẩn sử dụng để xác định độ dài quán tính PE tính " độ dài tắt dần" hàm tương quan góc liên kết (The bond angle correlation function - BACF) dọc theo trục polymer Nếu giả thiết hàm giảm theo hàm mũ dọc theo trục của dây theo công thức: (|� ′ − �|) = 〈cos|< (�� , ��′)|〉~ exp (− � |�′ − �| � ) (4.15) Trong �⃗⃗� �⃗⃗′� hai vecto liên kết thứ s s' < (��, �′�) góc ⃗⃗ ⃗⃗ hai vecto ký hiệu 〈… 〉 giá trị trung bình theo tất cấu hình polymer Để có sai số thống kê tốt cặp số s s' lấy trung bình dọc theo dây cách tịnh tiến số s s' giữ nguyên |�′ − �| Chúng chứng minh định nghĩa độ dài quán tính có phạm vi ứng dụng hạn chế PE Nếu �� nhỏ lớn kết độ dài quán tính thu từ BACF không đáng tin cậy Ngược lại phạm vi mà định nghĩa áp dụng chứng minh độ dài quán tính thu phù hợp với kết thu mục 4.2.1 4.2.2 trình bày Với �� nhỏ, hiệu ứng thể tích loại trừ đóng vai trò quan trọng việc giảm theo hàm mũ không dây polymer không dây Gauss Cụ thể là, chưa có phương pháp tính toán sử dụng để tách riêng hiệu ứng thể tích loại trừ thu �� số tắt dần Theo tài liệu tham khảo số [43], việc tắt dần dây polyemer có tương tác thể tích loại trừ hàm mũ trường hợp Việc xác định độ dài quán tính sử dụng BACF trở nên không đáng tin cậy �� lớn tương tác Coulomb không bị chắn kéo dài dây 51 A B -1 C -2 � 200 � x 400 Hình 4.11 Logarit hàm tương quan góc �(�) phụ thuộc vào khoảng cách x (trong đơn vị ��) dọc theo chuỗi PE cho trường hợp � = 512, �� = 50�� Trên hình (4.11) vẽ ln hàm BACF, hàm độ dài dọc theo trục PE cho trường hợp N=512 �� = 50��, với giá trị chọn hiệu ứng thể tích loại trừ xem tối thiểu Kết cho thấy có miền A, B C Chúng ta quan sát miền hình vẽ, miền A khoảng cách ngắn dọc theo trục PE monomer nằm sộ v ád ấ nà n hi "blob" tĩnh điện tương tác Coulomb nhỏ t cắ át cd BACF miền tắt dần độ dài liên ầ theo trục PE (vùng B) việc tắt dần theo hàm mũ độ dài tắt dần số rõ ràng Cuối khoảng cách lớn xấp xỉ độ dài nhanh (vùng C) Hiệu ứng ứng suất giảm 0, liên kết đầu cuối không tương quan Như độ dài quán tính quan tâm để so g n cuố i) trở ề nê n n ữ a lớn , s t t r ê p ýc h ủ vù h ta h B u A y C ế F t tr t ĩo t ng ì s B n kh h nn i hg � PE lần BACF lại giảm a ln kết �� Ở khoảng cách lớn dọc đ ứn i n g đủ k lớn h để � in ệ g n B l n M l n ôn h v đù có mộ t ( ) v độ ù dài n tắt g dần hằ C ng số c Có h ộ thể i đt v qu ế l bán kính chắn 10% độ dài dải Thậm chí có điều kiện chặt độ lớn �� phương pháp BACF không đáng tin cậy Cụ thể �� > � hàm BACF ��(�(�)) vùng B nhỏ Nghĩa hàm mũ không định nghĩa rõ ràng Do tất hạn chế này, nên dùng phương pháp BACF để tính �� dải giới hạn �� hiệu ứng thể tích loại trừ nhỏ �� không lớn so với L (việc giảm độ dài tắt dần vùng B nhỏ 0.1) Các giá trị độ dài quán tính dùng phương pháp liệt kê bảng (4.1) Bảng 4.1 So sánh tính toán độ dài ����� sử dụng phương pháp BACF tính toán tỉ số ���/�� mục 4.2 (trong đơn vị ��) N � � ��� 2048 1024 512 ���/� � � 100 4590 � 3682 150 9000 7484 80 2535 2180 100 3733 3111 50 1083 809 Trong bảng này, liệt kê kết ���⁄�� thu mục trước sử dụng công thức Bresler - Frenkel Hệ số �⁄�� cần thiết độ dài BACF ����� đo dọc theo trục polymer tự nhiên �� thu mục trước đo dọc theo trục dây polymer hiệu dụng Chúng ta nhận thấy phương pháp xác định độ dài quán tính khoảng chênh lệch từ 20 - 25% Đây phù hợp tốt cho thấy phương pháp sử dụng tương thích với Một lần kết phương pháp BACF nhấn mạnh kết thu mục trước Đó phương pháp cho kết phù hợp với lý thuyết OSF để mô tả PE dẻo 53 KẾT LUẬN Trong luận văn này, sử dụng phương pháp mô MC để nghiên cứu phụ thuộc độ dài quán tính tĩnh điện PE hàm bán kính chắn dung môi Chúng không mô polymer dài polymer nghiên cứu trước để thu dải tỉ lệ mà kết hợp việc sử dụng nhiều phương pháp phân tích kết khác có tính đến thăng giáng địa phương monomer Các phương pháp cho kết tương thích với phù hợp tốt với lý thuyết OSF Bên cạnh việc giới thiệu chi tiết lý thuyết tĩnh điện dung môi, PE (lý thuyết Poisson - Boltzmann, lý thuyết Debye - Huckel) trình bày kết lý thuyết OSF - KK nghiên cứu phụ thuộc độ dài tĩnh điện dây PE vào bán kính chắn, luận văn áp dụng phương pháp mô MC để mô chuỗi PE có thay đổi cường độ chắn dung dịch thu số kết mô đáng ý Cụ thể là: - Khảo sát phụ thuộc khoảng cách đầu cuối dây PE vào bán kính chắn cho dây PE có kích thước khác nhau: Từ dây có số monomer nhỏ (N =64) đến dây polymer có kích thước lớn (N=4096) - Tiến hành mô hệ lớn cho phép quan sát chuyển pha dải tỷ lệ, đặc biệt giải tỷ lệ ���~��6/5, kết quan trọng mà mô trước chưa thể rõ - Khớp kết mô cho trường hợp chắn yếu với lý thuyết biết để thu mật độ điện tích hiệu dụng độ dài quán tính cho trường hợp Đồng thời khảo sát chi tiết công trình nghiên cứu trước việc có tính đến giả thiết tranh "blobs" tích điện thăng giáng nhiệt độ monomer Kết cho thấy, bỏ qua thăng giáng độ dài quán tính tỉ lệ với bán kính chắn với số mũ nhỏ - Chúng sử dụng phương pháp BACF để khớp tính độ dài quán tính cho dây PE (mặc dù phương pháp thường áp dụng dây polymer 54 trung hòa) Chúng rằng, phương pháp áp dụng cho PE tích điện dải hạn chế bán kính chắn Trong dải bán kính chắn hạn chế phụ thuộc �� vào �� dung môi khác nhau, lần cho kết phù hợp với lý thuyết OSF - Dựa kết mô thu giải thích lý mô trước không cho kết phù hợp với lý thuyết OSF Như kết luận cuối sau kết nghiên cứu cho thấy, lý thuyết OSF - KK phù hợp để mô tả polymer tích điện Trong tương lai, trang bị hệ thống máy tính đại tiếp tục mô hệ có số monomer lớn hơn, hệ có monomer trung hòa monomer tích điện…để mô tả hoàn thiện tranh "blobs" tĩnh điện hi vọng tìm giải tỉ lệ với � = 55 TÀI LIỆU THAM KHẢO Tiếng Việt Trịnh Xuân Hoàng, Bài giảng môn mô mô hình hóa, Viện Vật lý, Viện KHVN Tiếng Anh Andy Wing - Chi Lau (2000), "Fluctuation ADN Correlation Effects in Electrostatics of Highly-Charged Surfaces", Physic, Stanta Barbara University of California Barrat J L., Boyer D (1993), J Phys II, 3, pp 343 Barrat J.-L., Joanny J.-F (1993), Europhys Lett, 24, pp 333 Barrat J.-L., Joanny J.-F (1996), Adv Chem Phys., 94, pp Bratko D., Dawson K A (1993), J Chem Phys., 24, pp 5352 Chapman.D.L (1913), Philos Mag., 25, pp 475 Debye P.W., H¨uckel E (1923), Z Phys., 24, pp 185 Fuoss R.M., Katachalsky A., Lifson S (1997), Proc Natl Acad Sci USA, 37, pp 275 10 Gouy G (1917), Ann Phys., 7, pp 129 11 Ha B.-Y., Thirumalai D (1995), Macromolecules, 28, pp 577 12 Ha B.-Y., Thirumalai D (1999), J Chem Phys., 110, pp 7533 13 Hao Li, Witten T A (1995), Macromolecules, 28, pp 5921 14 Hubbard J., (1959), Phys Rev Lett., 3, pp 77 15 Hsu H.-P., Paul W., Binder K (2010), "Standard Definition of Persistence Length Do Not Describe the Local Intrinsic Stiffness of Real Polymer Chains", Macromolecules, Vol.43, pp 3094-3102 16 Israelachvili J.N (1992), "Intermolecular ADN Surface Forces", Academic Press Inc., San Diego 17 IUPAC Compendium of Chemical Terminology (1997), the Gold Book, 2nd edition 18 Jo¨nsson.B., Peterson C., So¨derberg B (1995), J Phys Chem 99, pp 1251 56 19 Khokhlov A R., Khachaturian K A (1982), Polymer, 23, pp 1793 20 Lal M (1969), Mol Phys, 17, pp 57 21 Landau L D., Lifshitz E M (1996), Statistical Physics, Butterworth Heinemann Oxford 22 Lee Sell, Nguyen T T (2008), "Radial distribution of RNA genome packaged inside spherical viruses", Phys Rev Lett., 100, pp 198102 23 Lee Sell, Le Tung T., Nguyen T T (2010), "Reentrant Behavior of DivalentCounterion-Mediated DNA-DNA Electrostatic Interaction", Phys Rev Lett., 105, pp 248101 24 Lee Sell, Tran C V., Nguyen T T (2011), "Inhibition of DNA ejection from bacteriophages by Mg+2 counterion", J Chem Phys., 134, pp 125104 25 Madras N., Sokal A D., Stat.J (1988), Phys, 50, pp 109 26 Micka U., Kremer K (1996), Phys Rev E, 54, pp 2653 27 Netz R R., Orland H (1999), Eur Phys J B, 8, pp 81 28 Nguyen T T., Grosberg A Yu., Shklovskii B I (2000), "Lateral correlation of multivalent counterion is the universal mechanism of charge inversion", Electrostatic effects in Soft Matter and Biophysics., NATO SCIENCE SERIES: II: Mathematics, Physics and Chemistry, Volume 46, Kluwer Academic Publishers, Dordrecht 29 Nguyen T T., Grosberg A Yu., Shklovskii B.I (2001), "Screening and giant charge inversion in electrolytes", in More is different., edited by N P Ong and R N Bhatt, Princeton University Press, Princeton and Oxford, pp 285 30 Nguyen T T., Shklovskii B I (2002), "Persistence length of a polyelctrolyte in salty water: Monte Carlo study", Physical review E, 66, 021801 31 Nguyen T.T (2013), "Strongly correlated electrostatics of viral genome packaging", J Biol Phys., 39, pp 247 32 Odijk T., Sci J Polym (1977), Polym Phys Ed., 15, pp 477 57 33 Petra Bacova (2000), "Linear semiflexible polyelectrolytes in solution", Physical, Macromoleculer chemistry, Faculty of science, Charles University in Pargue 34 Reed C E., Reed W F (1991), J Chem Phys., 94, pp 8479 35 Rubinstein M., Colby R H (2003), Polymer Physics, Oxford University Press 36 Safran S.A (1994), "Statistical Thermodynamics of Surfaces, Interfaces, ADN membranes", Reading, Addison-Wesley Publishing Com 37 Seidel C., Bunsenges Ber (1996), Phys Chem., 100, pp 757 38 Skolnick J., Fixman M (1977), Macromolecules, 10, pp 944 39 Sierra-Martin B., Lietor-Santos J J., Fernandez-Barbero A., Nguyen T T., Fernandez-Nieves A (2011), "Swelling thermodynamics of microgel particles", In: Microgels Suspension - Fundamentals and Application, Editors: A FernandezNieves, H M Wyss, J Mattsson, D A Weitz 40 Stratonovitch R L (1957), Dokl Akad Nauk USSR, 115, pp 1907 41 Ullner M., Jo¨nsson B., So¨derberg B., Peterson C (1996), J Chem Phys, 104, pp 3048 42 Ullner M., Jo¨nsson B., Peterson C., Sommelius O., So¨derberg B (1997), J Chem Phys., 107, pp 1279 43 Ullner M., Woodward C E (2002), Macromolecules, 35, pp 1437 44 Verwey E., Overbeek J.Th.G (1948), "The Theory of the Stability of Lyophobic Colloids", Elsevier, Amsterdam ADN New York 45 Vu Van Quyen, Hoang Thuy Nga, Nguyen The Toan (2015), "Condensation of a polyelectrolyte by an oppositely charged spherical macro - ion", Student Scientific conference, Faculty of Physic, Vietnam National University 58 PHỤ LỤC Trong mô phỏng, thực code tính toán dài phức tạp Phần phụ lục, giới thiệu thuật toán pivot flip - Những thuật toán quan trọng để thay đổi cấu hình polymer Thuật toán (phép dịch chuyển) Flip void flip_move ( const int site ) { double x , y , z , t [ ] , deltaEC , deltaER , deltaE , temp , scalar ; double alpha1 , alpha2 , dx1 , dx2 , dy1 , dy2 , dz1 , dz2 ; double a11 , a21 , a31 , a12 , a22 , a32 , a13 , a23 , a33 ; double cpst , spst , cost , gamma, cosg , sing , omcg ; int i , j , idx , isite , cw ; if ((site )&&(site != (Nbead −1) ) ) { cpst = pos [ (site +1)∗3+0] − pos [ (site −1) ∗ 3+0]; spst = pos [ (site +1)∗3+1] − pos [ (site −1) ∗ 3+1]; cost = pos [ ( sit e +1) ∗3+2] − pos [ (site −1) ∗ 3+2]; temp = sqrt ( sqr (cost )+sqr (spst )+sqr ( cpst ) ) ; cpst /= temp ; spst /= temp ; cost /= temp ; gamma = DGMAX∗ ( two ∗ arand ( ) −one ) ; cossin (gamma, cosg , sing ) ; omcg = one −cosg ; temp = cost ∗ omcg ; a11 = cpst ∗ cpst ∗ omcg+cosg ; a22 = spst ∗ spst ∗ omcg+cosg ; a33 = cost ∗ temp+cosg ; t [ ] = spst ∗ cpst ∗ omcg ; t[ ] = cos t ∗ sing ; a12 = t [0] − t [ ] ; a21 = t [0]+ t [ ] ; t [ ] = cpst ∗ temp ; t [ ] = spst ∗ sing ; a31 = t [0] − t [ ] ; a13 = t [0]+ t [ ] ; t [ ] = spst ∗ temp ; t [ ] = cpst ∗ s i n g ; a23 = t [0] − t [ ] ; a32 = t [0]+ t [ ] ; t [ ] = pos [ site ∗ 3+0] − pos [ (site −1) ∗ 3+0]; 59 t [ ] = pos [ site ∗ 3+1] − pos [ (site −1) ∗ 3+1]; t [ ] = pos [ site ∗ 3+2] − pos [ ( sit e −1) ∗ 3+2]; x = pos [ ( site −1) ∗3+0]+ t [ ] ∗ a11+t [ ] ∗ a21+t [ ] ∗ a31 ; y = pos [ ( site −1) ∗3+1]+ t [ ] ∗ a12+t [ ] ∗ a22+t [ ] ∗ a32 ; z = pos [ ( site −1) ∗3+2]+ t [ ] ∗ a13+t [ ] ∗ a23+t [ ] ∗ a33 ; } else i f ( site ==0) { rand_unitvector (&cpst ,& spst ,& cost ) ; x = pos [1 ∗ 3+0] − R0 ∗ cpst ; y = pos [1 ∗ 3+1] − R0 ∗ spst ; z = pos [1 ∗ 3+2] − R0 ∗ cost ; } else { rand_unitvector (&cpst ,& spst ,& c o s t ) ; x = pos [ (site −1) ∗3+0]+R0 ∗ cpst ; y = pos [ ( site −1) ∗3+1]+R0 ∗ spst ; z = pos [ (sit e −1) ∗3+2]+R0 ∗ cost ; } deltaEC = zero ; #pragma omp parallel for private ( idx ) reduction (+: deltaEC ) for ( j =0; j site ) idx = t_idx [ site ] + j − site − 1; else continue ; EC_t_temp [ j ] = ECoulomb ( x−pos [ j ∗ 3+0] , \ y−pos [ j ∗ 3+1] , \ z −pos [ j ∗ 3+2] , \ ZZbead , \ RbeadMin2 ) ; deltaEC += EC_t_temp [ j ] −EC_t [ idx ] ; } if ( deltaEC > Infinity _ ) { / ∗ over la pping ∗ / return ; } deltaE = deltaEC ∗ KCoulomb ; if ( deltaE >0) { if ( log ( arand ( ) ) > −Tinv ∗ deltaE ) { / ∗ reject ∗ / return ; 60 } } / ∗ accept the move ∗ / pos [ site ∗ 3+0] = x ; pos [ site ∗ 3+1] = y ; pos [ site ∗ 3+2] = z ; EC+=deltaEC ; PE_accept++; flip _ accept ++; for ( i =0; i 1)) { x = pos_temp[(site+2*direction)*3+0]- \ pos_temp[(site+ direction)*3+0]; y = pos_temp[(site+2*direction)*3+1]- \ pos_temp[(site+ direction)*3+1]; z = pos_temp[(site+2*direction)*3+2]- \ pos_temp[(site+ direction)*3+2]; alphatemp2 = (-x2*x-y2*y-z2*z)/R0_2; if(fabs(alphatemp2)one?zero:M_PI; deltaER + = sqr(alphatemp2-M_PI)- \ sqr(alpha_t[site+direction-1]-M_PI);}} #endif deltaE=deltaEC*KCoulomb; if(deltaE>0) { if(log(arand()) > -Tinv*deltaE) { /* reject move */ return; }} /* accept move */ EC += deltaEC; ree = Ree(); PE_accept++; if (direction-1) { /* Update new position */ Memcpy(pos, pos_temp,site*3); } else { /* Update new position */ Memcpy(pos+(site+1)*3, pos_temp+(site+1)*3,(Nbead-(site+1))*3); }} 64 ...ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - Nguyễn Thị Trang NGHIÊN CỨU MÔ PHỎNG CỦA SỐ BÀI TOÁN VỀ TĨNH ĐIỆN TRONG SINH HỌC Chuyên ngành: Vật lí lí thuyết vật lí toán. .. lí toán Mã số: 60 44 01 03 LUẬN VĂN THẠC SĨ KHOA HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC: TS NGUYỄN THẾ TOÀN Hà Nội - 2 015 LỜI CẢM ƠN Để hoàn thành luận văn "Nghiên cứu mô số toán tĩnh điện sinh học" , xin... giả lựa chọn đề tài "Nghiên cứu mô số toán tĩnh điện sinh học" đề tài cho luận văn Trong luận văn này, sử dụng mô Monte - Carlo máy tính để xem xét phụ thuộc độ dài tĩnh điện le khoảng cách đầu