适用于DNA电荷输运研究的紧束缚模型方法

来源 :山东大学 | 被引量 : 0次 | 上传用户:chenshuae9o
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
作为生命遗传信息的载体,DNA是一类非常重要的生命大分子,其对绝大多数生命体的发育和机体运作起指导作用。近年来,由于在基因损伤修复等许多基本生命过程中的关键作用和在纳米电子学中作为分子器件的重要应用前景,DNA分子中电荷输运问题已成为多个学科比如化学、物理、生物前沿领域中的研究热点。在研究初期围绕DNA分子是否具有导电能力这一问题研究者们进行了诸多争论,不同实验中DNA表现出导体、半导体、绝缘体甚至是超导体等相差悬殊的导电行为。经过十多年来的研究,尽管对具体机制仍存在争议,研究者们已经基本认同DNA分子具有输运电荷的能力,但其程度受DNA分子的序列、碱基堆叠的完整性、构型涨落、溶剂环境以及电极(或电荷给体受体)与碱基的耦合等因素的影响很大,其中某些因素对DNA导电性的影响是颠覆性的。这些因素由分子微观结构决定,在实验上不容易控制,这就迫切要求理论计算学家通过模拟计算从微观角度研究DNA分子中的电荷输运行为,解释电荷输运机制。对DNA这类生命大分子体系的电子结构的计算一直是理论计算领域中一个具有挑战性的课题,主要是因为生命大分子体系所含原子数目太多且不具备周期性,利用常规量子化学从头算方法所需的计算条件当前计算机能力尚不能承受。在DNA电荷输运的计算模拟研究中,紧束缚模型方法因其简单明了、物理图像清晰等优点而被广泛采用。该方法把DNA分子抽象为一维格点体系,每个格点提供一个在位轨道,电子在紧邻格点轨道间能发生跃迁,其状态通过在位能和迁移积分描述。通过设定这两个参数的值,紧束缚模型方法可以计算DNA中电荷输运所涉及的电子状态、空穴状态或前线电子能带结构。尽管在DNA的电荷输运研究中紧束缚模型方法被广泛应用,但对其方法本身的研究还很少,这导致模型中参数取值缺少标准,各种文献中所采用的参数不一致,严重影响了该方法的准确性,不利于DNA电荷输运的研究。针对这一问题,我们从最基本的理论公式出发,根据DNA的分子结构和电子结构逐级做近似,推导出参数计算公式,建立系统化的参数计算方案,完善适用于DNA分子的紧束缚模型方法,使其能够更加合理地应用于DNA电荷输运研究。我们提出以下参数计算方案:第一,对孤立格点体系作量子化学计算以得到格点在位轨道;第二,对格点及其紧邻格点组成的子体系作量子化学计算,得到该子体系的单电子哈密顿作为计算参数专用的有效DNA哈密顿,总体系中其他原子作为点电荷作用于该有效哈密顿或直接忽略;第三,计算有效DNA哈密顿基于格点在位轨道的矩阵元,将矩阵元带入相应公式得到在位能和迁移积分。根据所提出的参数计算方案,我们对DNA的紧束缚模型方法做了以下研究:通过计算各种序列组合下空穴和电子的紧束缚模型参数,研究了DNA分子的一级结构(碱基序列)对参数的影响。计算结果表明,格点在位能主要取决于格点代表的碱基种类,四种碱基的空穴在位能的顺序是G<A<T<C,电子在位能的顺序恰好相反。在位能在其紧邻格点极化作用下发生能级分裂。对空穴来说,紧邻格点的极化作用能导致格点在位能的降低,其中3′侧紧邻格点的极化作用大于5′侧格点,四种碱基中G碱基的极化能力最强。理想B-型构象下格点间迁移积分取值在0.02-0.12eV,主要取决于相邻格点轨道之间的重叠。我们给出的这些参数可以用来构建任意序列DNA分子的紧束缚模型哈密顿,描述其中的空穴或电子状态,进而研究序列对DNA分子电荷输运行为的影响。以DNA分子中碱基对间距离和扭转角为代表,研究了DNA分子的二级结构(构象)对紧束缚模型参数的影响。计算结果表明,格点间迁移积分受构象尤其是扭转角的影响非常大。当格点间距离取标准值而扭转角为0°时,GG格点间和AA格点间的迁移积分均能达到近0.8eV,而当扭转角取标准值36°时,两类迁移积分均不足0.1eV。当固定扭转角而增加格点间距离时,格点间迁移积分随着距离的增加呈指数函数形式衰减。计算过程中我们还考虑了紧邻格点极化作用对迁移积分的影响,当格点二体体系明显不对称时,格点间的相互极化差异较大,导致不能利用前线轨道的能级分裂来计算迁移积分。另外,我们还分别计算了从大量晶体数据统计分析得出的平均A-型构象和平均B-型构象中所有序列组合对应在位能和迁移积分,比较了其差异。我们把计算出的参数应用于SSH模型,研究了DNA分子的势阱势垒结构对其中空穴极化子的影响。研究表明单阱结构中空穴以极化子的形式存在于GC碱基对组成的势阱中;当势阱跨度较窄时,空穴电荷在势垒上有一定分布;势阱宽度为3个GC碱基对时,极化子最稳定。在周期性阱垒结构中,势阱和势垒的宽度对极化子形状都有显著影响。通过拟合格点间迁移积分同格点间距的关系,我们得出DNA链polyA-polyT和polyG-polyC中电子-声子耦合因子,然后将其应用于SSH模型以研究这两种链中空穴极化子的静态性质及其在电场下的动态行为。计算结果表明polyG-polyC中的空穴极化子的定域程度、稳定性均强于polyA-polyT中极化子,但在电场作用下polyA-polyT中极化子的迁移率较高,利于长程电荷迁移现象的发生。除了对DNA分子的紧束缚模型方法研究外,我们还对之前提出的利用点群对称性简化矩阵元计算的生成子方法进行了一定研究。按照初始基函数在对称操作下的变换关系,可将初始基函数空间划分为对称性恒定子空间。在此基础上我们利用群表示理论和投影算符的性质从理论上证明了生成子方法,并将其推广到点群的高维不可约表示。作为演示,我们把该方法应用于苯分子的紧束缚模型计算和Heisenberg半经验价键模型计算中。本论文的创新点主要包括:1.从最根本的量子力学方程出发,通过逐级近似,得到了紧束缚模型参数的计算公式,据此建立了系统化的参数计算方案,完善了适用于DNA分子的紧束缚模型方法。2.根据所提出的参数计算方案,研究了DNA分子紧束缚模型中在位能和迁移积分随DNA分子的序列、构象的变化关系。3.根据在位能计算结果构建DNA势阱势垒结构,并利用SSH模型研究了其中空穴极化子的性质;根据迁移积分同格点间距离的关系得出SSH模型中电子-声子耦合因子。4.基于对称性恒定子空间的特性,利用群表示理论证明了之前提出的生成子方法,并将其推广到点群的高维不可约表示,使这一方法得以完善。
其他文献
目的探讨生态免疫肠内营养支持对溃疡性结肠炎(UC)患者临床症状、营养状况、免疫功能的影响。方法采用随机、对照研究,将106例溃疡性结肠炎患者随机分为对照组(常规药物+少渣
目的探讨通过采取预防措施降低老年髋部骨折患者压疮的发生率。方法对60例老年髋部骨折患者入院后行Braden评分,分析原因,采取针对性护理措施。结果 60例患者中1例出现Ⅰ度压
2007年全球金融危机以来,中央4万亿刺激政策和地方政府不断涌现的各种投资平台使国民经济取得惊人增长,但房地产泡沫和地方债务负担不断加剧。党的第十八届三中全会提出全面
考虑全样本场合三重两参数威布尔分布混合模型,提出了一个优选准则,从而得到参数的逆矩估计和区间估计;并通过实例说明该方法的简便、易行性。
以黄芪为试材,用微波消解仪提取黄酮,分光光度法测定黄酮提取率,采用单因子试验及正交实验研究黄芪黄酮的最佳提取条件。结果表明:乙醇浓度95%,提取时间20min,料液比1∶15g/m
小组合作在美术教学中起到积极的作用,是当前教学中新兴的一种教育模式,能够促进学生形成良好的学习习惯,有利于课堂氛围的提升。美术是对学生进行美育的重要课程,小组合作能够更
目的探讨平阳霉素(PYM)联合透明质酸钠(HA)对体外培养的人淋巴管畸形内皮细胞(HLMECs)的作用。方法对新鲜婴幼儿淋巴管畸形组织进行HLMECs原代培养。利用淋巴管内皮特异性标
1936年10月,中国共产党领导的那场震惊中外的红军长征,是中国革命史上的一座丰碑。红军长征的胜利,不仅为中国革命的胜利奠定了坚实的基础,而且对世界也产生了深远的影响。长征中所表现
在人-车-环境组成的闭环驾驶系统中,驾驶员由于驾驶技能良莠不齐、生理与心理状态、性格与体力情况等诸多方面各不相同,被认为是驾驶系统中不确定因素最多,最不可靠的环节。而智能汽车通过通信设备,传感器和智能算法替代了人在车辆行驶中的作用,被认为可以极大的解决交通安全问题。而路径跟随控制作为智能汽车领域的关键技术之一,在车辆智能化的发展下,已经成为了研究的热点。当前大部分关于智能汽车路径跟随控制的研究方法