含氢小分子是一种在自然界中广泛存在的简单物质,对它们的研究可以提高人们对物质科学的基本认识。同时,由于含氢小分子在地球内部和其他星体内部大量存在,其在压力下的结构和性质的研究也具有很重要的实际意义。通过压力,含氢小分子在常压下的结构、性质可以被大大的改变,从而带来很多新奇的现象和性质。含氢小分子在压力的影响下会发生复杂的结构转变,呈现出一系列性质迥异的新相。含氢小分子在压力下存在着的一些相似的变化规律,比如压致分子解离,更是引起了人们的广泛关注。其中最著名的就是水的压致解离。大量的高压实验和理论研究证明,水的压致解离和水中的氢原子在其压致解离曲线周围的三个高压固体相中的动力学行为息息相关。在水解离前的两个高压相VII相和VIII相中,由于氢原子处在特殊的双势阱中运动,从而表现为在低温相VIII时,氢位于其中一个势阱中,表现为有序相;在常温和高温时,由于热运动,氢原子在两个势阱间等几率分布,而其运动就表现为氢原子在两个势阱点来回运动。随着压力的升高,由于氧氧间距被压缩,这种特殊的双势阱势垒渐渐变成单势阱,从而约束了氢原子的运动,使得氢原子处在两个氧原子中间点的位置,形成了对称氢键,从而完整了分子解离。在水之后,人们发现其姐妹分子硫化氢——在分子水平上和水分子相似——在高压下却变现出迥异的性质。首先水在常压下具有很强的氢键,而硫化氢却没有;其次,硫化氢的高压相图远没有水的高压相图复杂;最重要的是,硫化氢在高压下将发生分解,从而形成硫、氢和硫化氢的混合物,而不会生成具有对称氢键的分子解离相。前人的实验虽然定出了硫化氢的解离曲线,但是由于混相的缘故,无法得到硫化氢在解离后的晶体相结构,也没有像水那样给出解离的微观机制。这给予了我们研究动机及研究思路,即:通过从头算分子动力学方法,从硫化氢未解离前的结构出发,先得到稳定的相结构,再通过改变温度压力条件,模拟其解离过程。我们首先从一个实验建议的结构出发(其具有I41/acd空间对称群),在15GPa,100K下进行动力学模拟。计算发现实验建议的结构并不稳定,在200fs的动力学模拟时间内,其氢键键连网络发生了变化。新得到的结构具有Ibca的对称性。我们得到的结构的XRD谱和实验符合的非常好。实验给出了错误的IV相结构的原因是无法测定氢原子的位置,而是在解出了硫原子的位置后手动加入了氢原子,从而得到了不稳定的氢键键连网络结构。我们接着从Ibca结构出发,通过升温的方式穿过实验建议的解离边界。之所以这样做是因为:1.硫化氢的解离压力随温度变化非常大。前人实验发现在室温时(300K),硫化氢在27GPa解离;而在150K时,解离压高达42GPa。这充分说明温度对硫化氢解离的重要性;2.实验建议,在压力超过25GPa时,由于硫硫原子间距被压缩到硫原子轨道相接触的程度,会促使硫化氢分解而不是解离,正像实验中所观察的那样。因此我们选择在15GPa时,通过升温动力学穿过扩展的解离区线,从而期望能得到硫化氢的V相。我们的Parrinello-Raham变胞动力学表明,在15GPa,350K时,动力学超胞和原子均方根位移同时发生了突变,得到了一个新的稳定相。同时,通过对分布函数(RDF)和角关联(ACF)函数的计算,我们发现,在15GPa,350K时,S-S RDF的第一主峰和第二主峰合并为一个峰,显示了硫晶格的一个突变;S-H RDF的第一主峰和第二主峰在300K以下时有明显的零值分界,因为它们分别代表S-H共价键和S…H氢键,但在350K时却有明显的交叠,这说明硫化氢的共价键和氢键有相互转化的迹象,也就是说,氢原子在两个S原子之间跳跃。而ACF函数显示,在350K时,硫化氢分子表现出明显的转动无序。由于存在明显的氢无序,于是我们先求出平均硫原子晶格,得到其晶格对称群为P63/mmc。接着,我们将9个ps的平衡动力学轨迹中超胞内的氢原子通过P63/mmc的周期性边界条件移动到其单胞内,从而得到一个时间和空间平均的氢原子空间分布。我们发现,氢原子空间密度分布在每个硫原子周围都有12个密度极大值。而且这些极大值都位于P63/mmc对称群的Wyckoff位置上。综上结果我们认为,硫化氢在15GPa350K时由IV相转变到V相。我们得到的新结构是一个纯硫化氢晶体相,具有P63/mmc对称群。氢原子在V相中是无序的,并且以分数占据的形式分布在硫原子周围。计算得到的V相的XRD谱也和实验符合。在V相中,氢无序主要表现在两方面,一是H原子在两个S原子之间的跳跃,而这和水的VII相没有不同;二是硫化氢分子的转动无序,而这在水的VII相中不曾发现。我们认为新发现的这个具有氢无序的相结构是分析和理解硫化氢的压致分解的关键。随着压力的升高,和水不一样,氢原子不会被固定在两个硫原子中间,从而不会随着硫原子间距的减小而被固定在两个硫原子之间而形成对称氢键。相反的,由于转动无序的存在,氢原子会被带离两个硫原子之间,甚至是形成两个硫原子直接接触的情况。再加上由于压力带来的硫硫轨道交叠,硫硫共价键终将形成。因此,我们发现的硫化氢的V相具有水的VII、VIII相不具有的动力学性质,从而导致硫化氢在高压下分解而不是形成对称氢键的分子解离相。含氢小分子引起广泛关注的另一个原因涉及到能源问题。现在的能源结构是建立在化石能源(包括煤、石油、天然气等)基础上的,其缺陷就是污染环境,并且储量有限。而氢能源由于其燃烧生成的是水,无任何污染,因此是未来的能源材料。应用氢能源的途径主要是用其做载体,因为地球上的氢主要是以稳定化合物的形式存在,不能直接利用。用氢做能源载体需要满足无污染和方便简单。特别是像汽车等小型机械设备,对以氢为载体的能源应用来说是再合适不过。因此,寻找合适的材料去储存氢是迫切之需。现在储存氢的方式主要集中在将氢固化方面,因为无论是气态还是液态形式的氢,都不利于实际应用。现在对氢的固化储存主要集中在三个方面,一个是利用金属和氢化合生成金属氢化物。金属氢化物的优点就是,通过加热可以很快的放出氢气,并且可以反复充氢,达到循环利用。缺点就是储氢比例不高,且金属成本较高,质量较大不便运输。第二个是通过各种纳米材料做成的孔筛式载体,使得氢吸附在纳米材料表面。但此方法同样面临着材料合成方面的难度。第三个方法是化学储氢,就是通过轻元素和氢的化合物,例如含氢小分子,达到储氢的目的。化学储氢的优点就是材料来源广泛,并且由于轻质元素(例如硼、碳、氮等)较低的质量数,能使氢的储量百分比达到一个很高的值。化学储氢的缺点就是放氢速率不高,放氢温度较高,还有就是重复利用较困难。综合来看,现在还没有一种方式可以在所有的方面达到最合理的要求。对于化学储氢来说,其将来的研究重点就是如何提高储放氢效率,从而达到重复利用的目的。化学储氢里最具前景的材料之一就是氨硼烷(NH3BH3),它的含氢体积百分比和质量百分比都非常高。已有的实验研究主要集中在改进催化剂和应用纳米材料改性等方面。但直到目前为止,还未能达到实际应用的要求。同时,对于氨硼烷本身的基础研究也大量开展起来。包括从实验上测定其结构,观察结构随温度和压力的变化,拉曼光谱的温度和压力效应等。人们发现在氨硼烷中,最重要的相互作用是双氢键,它是由带正电的、和氮相连的+δH与带负电的、和硼相连的Hδ-之间相互吸引造成的。而固态氨硼烷和气态氨硼烷最重要的区别就是在气态时,氨硼烷分子的旋转是整个分子作为一个整体而来的;在固态时,其氨基的旋转和硼基的旋转却是分离的。氨硼烷在常压下有两个固态相。其中低温相具有Pmn21对称群,常温相具有I4mm对称群。在更高的温度下,氨硼烷将分解,释放氢气。氨硼烷的高压研究也很重要,因为压力带来的普遍的结构和性质变化可能对其储氢性能有所提高。更重要的是,通过研究此类储氢材料的压力效应,可以提高人们对储氢材料的理解,进而对储氢材料的设计提供重要参考。人们发现在常温下加压时,氨硼烷会从I4mm相转变到一个具有Cmc21对称群的高压相,转变压大概在1.1GPa~1.4GPa,误差主要来自于不同实验小组用的样品的纯度、性状的不同。在更高的压力下,拉曼实验认为氨硼烷将在5GPa,12GPa时分别发生相变,而XRD实验认为直到12.1GPa,氨硼烷只会发生在~1.1GPa时的一个相变。在更高的压力下,没有相应的研究。为了解释实验上的差异和寻求氨硼烷在更高压力下的结构和动力学性质,我们应用从头算分子动力学,首先从氨硼烷已知的低温常压相出发,通过缓慢加压,发现在3GPa、100K时,体系发生了结构性转变。得到的新的稳定结构具有Cmc21对称性,和常温XRD实验发现的高压相一致。接着,我们在常温条件下对Cmc21结构进行加压。动力学模拟显示,Cmc21相在15GPa之前都保持稳定。而在15GPa时,动力学计算显示,体系的超胞在8个皮秒的模拟时间内没有收敛,而是表现为a、b两个晶轴交替振荡。我们发现这种交替振荡正好对应着Cmc21相单位晶格的一个剪切形变。于是我们又进行了弹性常数的计算,发现体系在~15.3GPa时发生C66的软化,而且C66正好对应着动力学模拟中发现的超胞的形变。当压力继续升高到20GPa,动力学显示超胞迅速收敛到稳定值,得到了一个具有P21对称群的新结构,由于该结构单包中有两个分子,我们将其记为P21(Z=2)。随着压力继续升高,在50GPa、300K时体系又发生了一次结构性转变,体系的超晶胞和均方根位移都发生了突变,得到了一个新相,同样具有P21对称群。由于其单包中有4个分子,记为P21(Z=4)结构。我们的分子动力学一直进行到60GPa,没有发现任何分子解离或分解的迹象。我们还进行了零温焓差曲线的计算。计算结果和动力学结果是大致符合的。通过零温结构分析,我们发现,在Cmc21相中,体系的双氢键数目并不是如最近一篇XRD实验研究文章所宣称的那样一直保持在12个每分子上,而是随着压力的升高而增多。我们利用一篇CSD结构库搜索研究文献提出的双氢键的作用范围和作用角度为判据,得出在6.4GPa时,体系中每个分子参与的双氢键数目从12增大到14。我们认为因为这种双氢键数目的变化而带来的结构变化,对应着Lin等人在拉曼光谱实验中发现的光谱突变,并且这种变化并没有造成体系相结构的变化。当体系转变到P21(Z=4)结构时,我们发现氨硼烷分子的交叉排列构象被破坏。虽然这种分子稳定构象的破坏会带来能量的升高,但相变同时双氢键网络的变化却大大降低了体系的能量。通过对Pmn21、Cmc21、P21(Z=2)和P21(Z=4)结构的动力学模拟轨迹的分析我们发现,在所有相结构中,双氢键键角的NHH部分随压力变化不大,而BHH部分随着压力先增大,在P21(Z=4)相中,由于分子构象的改变,BHH角又大大减小。我们还发现在所有相结构中,双氢键键角的分布都符合前人提出的双氢键键和规律,即:NHH角较直(接近160度),BHH角较弯(接近120度)。最后,我们的从头算分子动力学研究表明,氨硼烷的三个高压相中都存在分子的120度跳跃旋转。由于氨硼烷分子的C3v对称性,这种旋转并不会破坏氨硼烷的有序结构。动力学显示,在常压相下存在的NH3和BH3基团旋转的分离在高压下被保留。并且保持着NH3旋转优先于BH3旋转和分子整体旋转的格局。通过分子旋转势能面的计算,我们发现这种格局可以被它们相差很大的旋转势垒所解释。