搜索

文章查询

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

TiAl电子态结构的ab initio计算

张树东 王传航 唐伟 孙阳 孙宁泽 孙召玉 徐慧

TiAl电子态结构的ab initio计算

张树东, 王传航, 唐伟, 孙阳, 孙宁泽, 孙召玉, 徐慧
PDF
HTML
导出引用
导出核心图
  • 应用完全活动基自洽场方法, 结合N电子价态微扰近似(NEVPT2), 对TiAl金属二聚体的基态和若干最低电子激发态的势能曲线进行了计算. 完全活动空间由Al的3个价电子(3s23p1)轨道和Ti的4个价电子(3d24s2)轨道构成, 计算基组选用Karlsruhe group的价分裂全电子基组def2-nZVPP(n = T, Q). 在确认TiAl的基态为四重态的基础上, 在核间距R = 0.200—0.500 nm范围内, 扫描获得了TiAl基态和最低二个激发态的完整势能曲线, 并对电子态进行了标识, 发现在0.255 nm附近存在电子态结构的“突变”. 在R > 0.255 nm区域, 基态和两个激发态分别为X4Δ, A4Π和B4Γ; 在R < 0.255 nm区域, 基态仍为X4Δ, 但两个激发态变为A'4Φ和B'4Π, 且存在激发态简并解除的现象. 基于NEVPT2修正后的势能曲线, 获得了TiAl电子态的平衡核间距、束缚能、激发能、跃迁偶极矩等特征参数, 并解释了实验上观测不到TiAl电子跃迁光谱的原因. 电子激发态存在“突变”的结构特征, 可为分析理解TiAl合金在室温下的脆性问题提供参考.
      通信作者: 张树东, zhangsd2@126.com
    • 基金项目: 国家自然科学基金(批准号: 11705101)资助的课题
    [1]

    Lewandowski J J, Seifi M 2016 Annu. Rev. Mater. Res. 46 151

    [2]

    Hug G, Loiseau A, Lasalmonie A 1986 Philos. Mag. A 54 47

    [3]

    Hong T, Watsonyang T, Guo X, Freeman A, Oguchi T, Xu J 1990 Phys. Rev. B 41 12462

    [4]

    Hall E L, Huang S C 1989 J. Mater. Res. 4 595

    [5]

    Hussain A, Hayat S S, Choudhry M A 2011 Physica B 406 1961

    [6]

    Yuan X, Yin S, Lian Y, Yan P Y, Xu H F, Yan B 2019 Chin. Phys. B 28 043101

    [7]

    Zhao H Y, Ma H M, Wang J, Liu Y 2016 Chin. Phys. Lett. 33 108105

    [8]

    Tang F D, Du Q H, Petrovic C, Zhang W, He M Q, Zhang L Y 2019 Chin. Phys. B 28 037104

    [9]

    Wan M J, Jin C G, Yu Y, Huang D H, Shao J X 2017 Chin. Phys. B 26 033101

    [10]

    Chen G, Peng Y, Zheng G, Qi Z, Wand M, Yu H, Dong C, Liu C T 2016 Nat. Mater. 15 876

    [11]

    李兴华, 杨绍利 2011 材料导报 25 94

    Li X H, Yang S L 2011 Mater. Reports 25 94

    [12]

    宋成粉, 樊沁娜, 李蔚, 刘永利, 张林 2011 物理学报 60 063104

    Song C F, Fan X N, Li W, Liu Y L, Zhang L 2011 Acta Phys. Sin. 60 063104

    [13]

    闫蕴琪, 张振祺 2000 材料导报 14 31

    Yan Y Q, Zhang Z Q 2000 Mater. Reports 14 31

    [14]

    Behm J M, Morse M D 1994 J. Chem. Phys. 101 6500

    [15]

    Behm J M, Brugh D J, Morse M D 1994 J. Chem. Phys. 101 6487

    [16]

    Behm J M, Arrington C A, Langenberg J D, Morse M D 1993 J. Chem. Phys. 99 6394

    [17]

    宋庆功, 秦国顺, 杨宝宝, 蒋清杰, 胡雪兰 2016 物理学报 65 046102

    Song Q G, Qin G S, Yang B B, Jiang Q J, Hu X L 2016 Acta Phys. Sin. 65 046102

    [18]

    Zope R R, Mishin Y 2003 Phys. Rev. B 68 024102

    [19]

    Hu H, Ren Y S, Wu X Z, Liu W G, Luo J J 2019 Int. J. Mod. Phys. B 33 1950097

    [20]

    Jeong B, Kim J, Lee T, Kim S W, Ryu S 2018 Sci. Rep. 8 15200

    [21]

    Ouyang Y, Wang J, Liu F, Liu Y, Du Y, He Y 2009 J. Mol. Struc. Theochem 905 106

    [22]

    Neese F 2012 Comput. Mol. Sci. 2 73

    [23]

    Gerard H, Davidson E R, Eisenstein O 2002 Mol. Phys. 100 533

    [24]

    Larsson S 2011 Int. J. Quantum Chem. 111 3424

    [25]

    Zhang S D, Wand M X, Wand Z F, Xu K, Dong J P 2017 J. Phys. Soc. Jpn. 86 074301

    [26]

    Weigend F, Ahlrichs R 2005 Phys. Chem. Chem. Phys. 7 3297

    [27]

    Barysz M 2016 J. Chem. Theory Comput. 12 1614

    [28]

    NIST Standard Reference Database 78 https://dx.doi.org/10.18434/T4W30F [2019-9-4]

    [29]

    Smith E B 2014 Basic Physical Chemistry: The Route to Understanding (London: Imperial College Press)

  • 图 1  SA-CAS (7, 10)对TiAl最低6条电子态势能曲线的扫描计算, 其中采用的基组分别是(a) def2-TZVP; (b) def2-TZVPP, (c) def2-QZVPP

    Fig. 1.  The lowest 6 potential energy curves of TiAl calculated by SA-CAS (7, 10) with basis set of (a) def2-TZVP; (b) def2-TZVPP; (c) def2-QZVPP.

    图 2  活动基轨道能量随核间距的变化

    Fig. 2.  Variation of CAS orbital energy with nuclear distance.

    图 3  活动基轨道上电子占居数随核间距的变化

    Fig. 3.  Variation of occupied electrons in the CAS orbital with nuclear distance.

    图 4  TiAl最低四重电子态的标识(SA-CAS计算结果)

    Fig. 4.  Identification of the lowest quadruple electronic state of TiAl (SA-CAS calculation results).

    图 5  TiAl电子态势能曲线的SC-NEVPT2动态相关修正计算

    Fig. 5.  Dynamic correlation correction calculation of SC-NEVPT2 for TiAl electronic potential energy curve.

    表 1  CAS (7, 10)/def2-TZVP计算的活动基分子轨道(MO14−MO23)系数(Eh = 2625.5 kJ/mol)

    Table 1.  Coefficients of the CAS orbital (MO14−MO23) calculated by CAS (7, 10)/def2-TZVP.

    MO No.14151617181920212223
    Energy/Eh–0.4035–0.18380.0029–0.0159–0.01590.07430.07430.04390.17140.1714
    Number of occupied electron1.9571.7720.6690.5950.5950.4940.4940.2220.0990.099
    Symbolσσσππδδσππ
    Tisσ12.147.57.7000013.400
    Tipzσ7.52.61.2000038.200
    Tipxπ0004.9000000
    Tipy0001.8000000
    Tidz2σ71.786.1000010.400
    Tidxzπ00031.411.400034.631.2
    Tidyz00011.431.400031.234.6
    Tidx2y2δ0000082.617.2000
    Tidxy0000017.282.6000
    Alsσ72.97.81.20000000
    Alpzσ0.539.23.2000036.600
    Alpxπ00034.412.500016.715.1
    Alpy00012.534.400015.116.7
    下载: 导出CSV

    表 2  两组${\text{π}}$轨道的组成分析

    Table 2.  Composition analysis of two π orbits

    OrbitalR = 0.200 nmR = 0.240 nmR = 0.280 nmR = 0.490 nm
    ${\rm{(1{\text{π}})(2{\text{π}})}}$${aligned}& {\rm{Ti(}}3{{\rm{p}}_x}, {\rm{ }}3{{\rm{p}}_y}{\rm{) }}7{\rm{\% }} \\& {\rm{Ti(}}3{{\rm{d}}_{xz}}, 3{{\rm{d}}_{yz}}{\rm{) }}60{\rm{\% }} \\ &{\rm{Al(}}3{{\rm{p}}_x}, 3{{\rm{p}}_y}{\rm{) }}28{\rm{\% }}{aligned} $${aligned}& {\rm{Ti(3}}{{\rm{p}}_x}, {\rm{ 3}}{{\rm{p}}_y}{\rm{) }}7{\rm{\% }} \\& {\rm{Ti(3}}{{\rm{d}}_{xz}}, 3{{\rm{d}}_{yz}}{\rm{) }}5{\rm{7\% }} \\& {\rm{Al(3}}{{\rm{p}}_x}, 3{{\rm{p}}_y}{\rm{) }}3{\rm{2\% }}{aligned} $${aligned}& {\rm{Ti(3}}{{\rm{p}}_x}, {\rm{ }}3{{\rm{p}}_y}{\rm{) 3\% }} \\ &{\rm{Ti(3}}{{\rm{d}}_{xz}}, 3{{\rm{d}}_{yz}}{\rm{) 73\% }} \\& {\rm{Al(3}}{{\rm{p}}_x}, 3{{\rm{p}}_y}{\rm{) 21\% }}{aligned} $${\rm{Ti}}(3{{\rm{d}}_{xz}}, 3{{\rm{d}}_{yz}}){\rm{ }}1{\rm{00}}\% $
    ${\rm{(3{\text{π}})(4{\text{π}})}}$${aligned} &{\rm{Ti(3}}{{\rm{d}}_{xz}}{\rm{, 3}}{{\rm{d}}_{yz}}{\rm{) 52\% }} \\ &{\rm{Al(3}}{{\rm{p}}_x}{\rm{, 3}}{{\rm{p}}_y}{\rm{) 36\% }}{aligned} $${aligned} &{\rm{Ti(3}}{{\rm{d}}_{xz}}, 3{{\rm{d}}_{yz}}{\rm{) 52\% }} \\& {\rm{Al(3}}{{\rm{p}}_x}, 3{{\rm{p}}_y}{\rm{) 40\% }}{aligned} $${aligned}& {\rm{Ti(3}}{{\rm{p}}_x}, {\rm{ }}3{{\rm{p}}_y}{\rm{) 12\% }} \\ &{\rm{Ti(3}}{{\rm{d}}_{xz}}, 3{{\rm{d}}_{yz}}{\rm{) 34\% }} \\ &{\rm{Al(3}}{{\rm{p}}_x}, 3{{\rm{p}}_y}{\rm{) 52\% }}{aligned} $${\rm{Al(3}}{{\rm{p}}_x}, 3{{\rm{p}}_y}{\rm{) 99\% }}$
    下载: 导出CSV

    表 3  R = 0.490 nm处活动基分子轨道MO14-MO23组成

    Table 3.  Composition of CAS orbitals MO14-MO23 at R = 0.490 nm.

    Orbital No121314151617
    Energy/Hartree–1.79778–1.7976–0.37486–0.215130.020010.03822
    Occupied electron2.000002.00001.919761.899000.624630.55980
    Ti s002.894.500
    Ti pz099.80.5000
    Ti px55.70.10000
    Ti py44.300000
    Ti dxz000055.843.6
    Ti dyz000044.155.2
    Al s0095.82.800
    Orbital No181920212223
    Energy/Hartree–0.00791–0.006820.08850.089790.056510.1227
    Occupied electron0.517860.517500.410580.406450.101590.04283
    Ti pz0000.091.45.3
    Ti dx2y2001.498.600
    Ti dxy0098.41.400
    Al pz0000792.7
    Al px55.343.30000
    Al py43.754.80000
    下载: 导出CSV

    表 4  基态及最低激发态的组态及跃迁偶极矩

    Table 4.  Configuration and transition dipole moment of the ground state and the lowest excited state

    R/nmstateMain configurationExcitation energy/cm–1Transition dipole moment T2/Debye2Possible quartet stateIdetified state
    0.285Ground state${{\rm{\text{σ} }}^{\rm{2}}}{{\rm{\text{σ} }}^{\rm{2}}}{{\text{π}}^{\rm{2}}}{{\rm{\text{δ} }}^{\rm{1}}}{{\text{π}}^{\rm{0}}}$04ΔX4Δ
    1st excited state${{\rm{\text{σ} }}^{\rm{2}}}{{\rm{\text{σ} }}^{\rm{2}}}{{\text{π}}^{\rm{2}}}{{\rm{\text{δ} }}^{\rm{0}}}{{\text{π}}^{\rm{1}}}$32120.0344ΠA4Π
    2nd excited state${{\rm{\text{σ} }}^{\rm{2}}}{{\rm{\text{σ} }}^{\rm{2}}}{{\text{π}}^{\rm{1}}}{{\rm{\text{δ} }}^{\rm{1}}}{{\text{π}}^{\rm{1}}}$346204Σ, 4Δ(2), 4ΓB4Γ
    0.240Ground state${{\rm{\text{σ} }}^{\rm{2}}}{{\rm{\text{σ} }}^{\rm{2}}}{{\text{π}}^{\rm{2}}}{{\rm{\text{δ} }}^{\rm{1}}}$04ΔX4Δ
    1st excited state${{\rm{\text{σ} }}^{\rm{2}}}{{\rm{\text{σ} }}^{\rm{1}}}{{\text{π}}^{\rm{3}}}{{\rm{\text{δ} }}^{\rm{1}}}$41400.008244Π, 4ΦA'4Φ
    2nd excited state${{\rm{\text{σ} }}^{\rm{2}}}{{\rm{\text{σ} }}^{\rm{1}}}{{\text{π}}^{\rm{3}}}{{\rm{\text{δ} }}^{\rm{1}}}$47270.008694Π, 4ΦB'4Π
    3rd excited state${{\rm{\text{σ} }}^{\rm{2}}}{{\rm{\text{σ} }}^{\rm{1}}}{{\text{π}}^{\rm{3}}}{{\rm{\text{δ} }}^{\rm{1}}}$50740.005514Π, 4ΦB'4Π
    下载: 导出CSV

    表 5  TiAl最低3个四重态的结构参数

    Table 5.  Structural parameters of the lowest three quadruple states of TiAl.

    StateRe/nmDe/cm–1
    CASNEVPT2CASNEVPT2
    X4Δ0.2880.26630168151
    A4Π0.320$\left\{{aligned}& {0.248} \\ & {0.296} {aligned} \right.$796$\left\{{aligned}& {3845} \\ & {3406} {aligned} \right.$
    B4Γ0.324$\left\{{aligned}& {0.248} \\ & {0.306} {aligned} \right.$711$\left\{{aligned}& {2884} \\ & {3406} {aligned} \right.$
    下载: 导出CSV
  • [1]

    Lewandowski J J, Seifi M 2016 Annu. Rev. Mater. Res. 46 151

    [2]

    Hug G, Loiseau A, Lasalmonie A 1986 Philos. Mag. A 54 47

    [3]

    Hong T, Watsonyang T, Guo X, Freeman A, Oguchi T, Xu J 1990 Phys. Rev. B 41 12462

    [4]

    Hall E L, Huang S C 1989 J. Mater. Res. 4 595

    [5]

    Hussain A, Hayat S S, Choudhry M A 2011 Physica B 406 1961

    [6]

    Yuan X, Yin S, Lian Y, Yan P Y, Xu H F, Yan B 2019 Chin. Phys. B 28 043101

    [7]

    Zhao H Y, Ma H M, Wang J, Liu Y 2016 Chin. Phys. Lett. 33 108105

    [8]

    Tang F D, Du Q H, Petrovic C, Zhang W, He M Q, Zhang L Y 2019 Chin. Phys. B 28 037104

    [9]

    Wan M J, Jin C G, Yu Y, Huang D H, Shao J X 2017 Chin. Phys. B 26 033101

    [10]

    Chen G, Peng Y, Zheng G, Qi Z, Wand M, Yu H, Dong C, Liu C T 2016 Nat. Mater. 15 876

    [11]

    李兴华, 杨绍利 2011 材料导报 25 94

    Li X H, Yang S L 2011 Mater. Reports 25 94

    [12]

    宋成粉, 樊沁娜, 李蔚, 刘永利, 张林 2011 物理学报 60 063104

    Song C F, Fan X N, Li W, Liu Y L, Zhang L 2011 Acta Phys. Sin. 60 063104

    [13]

    闫蕴琪, 张振祺 2000 材料导报 14 31

    Yan Y Q, Zhang Z Q 2000 Mater. Reports 14 31

    [14]

    Behm J M, Morse M D 1994 J. Chem. Phys. 101 6500

    [15]

    Behm J M, Brugh D J, Morse M D 1994 J. Chem. Phys. 101 6487

    [16]

    Behm J M, Arrington C A, Langenberg J D, Morse M D 1993 J. Chem. Phys. 99 6394

    [17]

    宋庆功, 秦国顺, 杨宝宝, 蒋清杰, 胡雪兰 2016 物理学报 65 046102

    Song Q G, Qin G S, Yang B B, Jiang Q J, Hu X L 2016 Acta Phys. Sin. 65 046102

    [18]

    Zope R R, Mishin Y 2003 Phys. Rev. B 68 024102

    [19]

    Hu H, Ren Y S, Wu X Z, Liu W G, Luo J J 2019 Int. J. Mod. Phys. B 33 1950097

    [20]

    Jeong B, Kim J, Lee T, Kim S W, Ryu S 2018 Sci. Rep. 8 15200

    [21]

    Ouyang Y, Wang J, Liu F, Liu Y, Du Y, He Y 2009 J. Mol. Struc. Theochem 905 106

    [22]

    Neese F 2012 Comput. Mol. Sci. 2 73

    [23]

    Gerard H, Davidson E R, Eisenstein O 2002 Mol. Phys. 100 533

    [24]

    Larsson S 2011 Int. J. Quantum Chem. 111 3424

    [25]

    Zhang S D, Wand M X, Wand Z F, Xu K, Dong J P 2017 J. Phys. Soc. Jpn. 86 074301

    [26]

    Weigend F, Ahlrichs R 2005 Phys. Chem. Chem. Phys. 7 3297

    [27]

    Barysz M 2016 J. Chem. Theory Comput. 12 1614

    [28]

    NIST Standard Reference Database 78 https://dx.doi.org/10.18434/T4W30F [2019-9-4]

    [29]

    Smith E B 2014 Basic Physical Chemistry: The Route to Understanding (London: Imperial College Press)

  • [1] 刘慧, 邢伟, 施德恒, 孙金锋, 朱遵略. PS自由基X2Π态的势能曲线和光谱性质. 物理学报, 2013, 62(20): 203104. doi: 10.7498/aps.62.203104
    [2] 于坤, 张晓美, 刘玉芳. 从头计算研究BCl+基态和激发态的势能曲线和光谱性质. 物理学报, 2013, 62(6): 063301. doi: 10.7498/aps.62.063301
    [3] 郭雨薇, 张晓美, 刘彦磊, 刘玉芳. BP+基态和激发态的势能曲线和光谱性质的研究. 物理学报, 2013, 62(19): 193301. doi: 10.7498/aps.62.193301
    [4] 黄多辉, 万明杰, 王藩侯, 杨俊升, 曹启龙, 王金花. GeS分子基态和低激发态的势能曲线与光谱性质. 物理学报, 2016, 65(6): 063102. doi: 10.7498/aps.65.063102
    [5] 罗华锋, 万明杰, 黄多辉. BH+离子基态及激发态的势能曲线和跃迁性质的研究. 物理学报, 2018, 67(4): 043101. doi: 10.7498/aps.67.20172409
    [6] 高 峰, 杨传路, 张晓燕. 多参考组态相互作用方法研究ZnHg低激发态(1∏,3∏)的势能曲线和解析势能函数. 物理学报, 2007, 56(5): 2547-2552. doi: 10.7498/aps.56.2547
    [7] 黄多辉, 王藩侯, 杨俊升, 万明杰, 曹启龙, 杨明超. SnO分子的X1Σ+, a3Π和A1Π态的势能曲线与光谱性质. 物理学报, 2014, 63(8): 083102. doi: 10.7498/aps.63.083102
    [8] 李瑞, 张晓美, 李奇楠, 罗旺, 金明星, 徐海峰, 闫冰. SiS低激发态势能曲线和光谱性质的全电子组态相互作用方法研究. 物理学报, 2014, 63(11): 113102. doi: 10.7498/aps.63.113102
    [9] 杨子元, 王参军, 许启明, 魏 群. 轴对称晶场中d3离子激发态对4A2基态自旋哈密顿参量的影响. 物理学报, 2007, 56(1): 507-511. doi: 10.7498/aps.56.507
    [10] 王新强, 杨传路, 苏涛, 王美山. BH分子基态和激发态解析势能函数和光谱性质. 物理学报, 2009, 58(10): 6873-6878. doi: 10.7498/aps.58.6873
    [11] 周 勋, 令狐荣锋, 杨向东, 朱正和, 吕 兵. MgH分子X2Σ+,A2Π和B2Σ+电子态的势能函数. 物理学报, 2008, 57(2): 816-821. doi: 10.7498/aps.57.816
    [12] 高双红, 任兆玉, 郭平, 郑继明, 杜恭贺, 万丽娟, 郑琳琳. 石墨烯量子点的磁性及激发态性质. 物理学报, 2011, 60(4): 047105. doi: 10.7498/aps.60.047105
    [13] 张锦芳, 任雅娜, 王军民, 杨保东. 铯原子激发态双色偏振光谱. 物理学报, 2019, 68(11): 113201. doi: 10.7498/aps.68.20181872
    [14] 陈恒杰. LiAl分子基态、激发态势能曲线和振动能级. 物理学报, 2013, 62(8): 083301. doi: 10.7498/aps.62.083301
    [15] 刘晓军, 苗凤娟, 李瑞, 张存华, 李奇楠, 闫冰. GeO分子激发态的电子结构和跃迁性质的组态相互作用方法研究. 物理学报, 2015, 64(12): 123101. doi: 10.7498/aps.64.123101
    [16] 蔡绍洪, 周业宏. 氯乙烯在外电场下的激发态结构研究. 物理学报, 2010, 59(11): 7749-7755. doi: 10.7498/aps.59.7749
    [17] 田原野, 郭福明, 曾思良, 杨玉军. 原子激发态在高频强激光作用下的光电离研究. 物理学报, 2013, 62(11): 113201. doi: 10.7498/aps.62.113201
    [18] 赵翠兰, 王丽丽, 赵丽丽. 有限深抛物势量子盘中极化子的激发态性质. 物理学报, 2015, 64(18): 186301. doi: 10.7498/aps.64.186301
    [19] 施德恒, 牛相宏, 孙金锋, 朱遵略. BF自由基X1+和a3态光谱常数和分子常数研究. 物理学报, 2012, 61(9): 093105. doi: 10.7498/aps.61.093105
    [20] 韩晓琴, 肖夏杰, 刘玉芳. HNO(1A’)自由基的从头算势能曲线. 物理学报, 2013, 62(19): 193101. doi: 10.7498/aps.62.193101
  • 引用本文:
    Citation:
计量
  • 文章访问数:  178
  • PDF下载量:  6
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-09-05
  • 修回日期:  2019-10-23
  • 上网日期:  2019-11-27
  • 刊出日期:  2019-12-01

TiAl电子态结构的ab initio计算

  • 曲阜师范大学物理工程学院, 曲阜 273165
  • 通信作者: 张树东, zhangsd2@126.com
    基金项目: 国家自然科学基金(批准号: 11705101)资助的课题

摘要: 应用完全活动基自洽场方法, 结合N电子价态微扰近似(NEVPT2), 对TiAl金属二聚体的基态和若干最低电子激发态的势能曲线进行了计算. 完全活动空间由Al的3个价电子(3s23p1)轨道和Ti的4个价电子(3d24s2)轨道构成, 计算基组选用Karlsruhe group的价分裂全电子基组def2-nZVPP(n = T, Q). 在确认TiAl的基态为四重态的基础上, 在核间距R = 0.200—0.500 nm范围内, 扫描获得了TiAl基态和最低二个激发态的完整势能曲线, 并对电子态进行了标识, 发现在0.255 nm附近存在电子态结构的“突变”. 在R > 0.255 nm区域, 基态和两个激发态分别为X4Δ, A4Π和B4Γ; 在R < 0.255 nm区域, 基态仍为X4Δ, 但两个激发态变为A'4Φ和B'4Π, 且存在激发态简并解除的现象. 基于NEVPT2修正后的势能曲线, 获得了TiAl电子态的平衡核间距、束缚能、激发能、跃迁偶极矩等特征参数, 并解释了实验上观测不到TiAl电子跃迁光谱的原因. 电子激发态存在“突变”的结构特征, 可为分析理解TiAl合金在室温下的脆性问题提供参考.

English Abstract

    • 过渡金属铝化物合金属于金属间合金, 是开发在恶劣环境中仍具有高强度和优异抗氧化性材料的重要候选物之一, 对其性质的研究正受到实验和理论计算的广泛关注[1-9]. 作为钛铝合金材料的基本单元, 对TiAl双原子体系几何结构和电子态结构的认识, 为了解该合金的性质提供有意义的重要参考. 钛铝合金在室温下的脆性问题一直困扰着人们, 也限制了其作为结构材料的使用[10-13]. Behm等[14-16]通过共振双光子电离光谱技术对第一周期过渡金属铝化物进行了实验研究, 但唯独没有观测到ScAl, TiAl和FeAl的光谱. 就是TiAl的基态, 也仅仅是参照MnAl的实验结果加以推断, 得出TiAl的基态可能是4Δ态. 理论研究方面, 目前多采用密度泛函方法计算合金的性质[17-20], 但对TiAl双原子体系的研究较少. Ouyang等[21]通过B3LYP密度泛函理论研究了3d金属铝化物的基态, 在6-311+G(2d)(Al)+Lan12DZ (过渡金属原子)基组下计算了键长、谐振动频率和解离能D0. 其中对TiAl的计算中, 在默认TiAl的基态为4Δ态下获得该态的键长、谐振动频率和解离能D0分别是0.277 nm, 234.4 cm–1和1.42 eV. 尚未见对TiAl激发态计算研究的报道.

      本文基于ORCA软件包提供的完全活动基自洽场(CASSCF)计算模块, 尝试对TiAl基态及其最低若干激发态的势能曲线开展了扫描计算, 并在考虑N电子价态微扰近似(NEVPT2)的基础上, 获得了TiAl最低3个四重电子态的特征参数, 分析讨论实验上观测不到TiAl电子光谱的原因.

    • 所有计算是使用ORCA程序包(版本4.1.1)在计算服务器上完成的[22]. 该程序包的突出优点之一就是提供了丰富的激发态计算手段. 使用该程序包, 在完全活动基自洽场(CASSCF)计算[23,24]的基础上, 通过非收缩的多参考组态相互作用(uncontracted MRCI)方法, 完成了对ZnAl最低几个电子态的势能曲线扫描计算[25]. 由于过渡金属Zn中的d轨道已被电子占满, 故活动空间的选取并未包含d轨道. 但最近在进行ScAl的计算时, 仅仅因为d轨道加入活动空间, 就因为遇到分子轨道收敛困难使得MRCI计算无法顺利进行. 这里报道对TiAl的计算, 也遇到类似的问题, 即用MRCI方法进行动态相关修正计算时, 存在收敛困难或者存在势能曲线不光滑. 为此, 尝试采用强收缩的N电子价态微扰理论(strong contracted N-electron valence state perturbation theory), 即SC-NEVPT2 (简称NEVPT2), 对CAS的计算结果进行动态相关修正计算, 获得了较为理想的势能曲线图像. 同时也发现, 采用Weigend 和Ahlrichs[26]提出的Karlsruhe全电子分裂价基组, 例如def2-TZVP, def2-TZVPP和def2-QZVPP, 均能完成势能曲线的扫描计算, 但def2-TZVP的计算精度明显不够, 能够用来分析讨论的结果至少需要在def2-TZVPP及以上的基组.

      CAS计算的核心问题是如何选取活动基空间轨道. 对于势能曲线扫描计算, 由于核间距的变化, 活动空间轨道的组成和次序也会发生变化[27]. 本文的计算中, 首先用密度泛函方法对基态构型进行了优化, 即采用B3LYP/def2-TZVPP优化基态构型, 获得基态的平衡核间距Re. 在保持该核间距不变的条件下, 用MP2方法进行单点能计算并对分子轨道进行分析, 生成自然轨道. 将该自然轨道作为CAS计算中初始轨道猜测输入. 分析CAS计算的活动基轨道组成, 最终确定选用了主要由Al和Ti的10个价电子所在轨道构成的活动空间, 即CAS (7, 10). 在该空间下, 用态平均尝试计算最低的四重态能级, 根据能级分布特征, 选择了最低3个具有二重简并特征的电子态作为势能曲线扫描计算对象. 在顺利完成CAS扫描计算的基础上, 进而选择NEVPT2修正计算.

    • 由原子数据库数据[28]可知, Al原子的基电子组态和基电子态分别是KL3s23p12Pu, Ti原子的基电子组态及基电子态为KL3s23p63d24s2和a3Fg. 按照分子理论[29], 原子群表示分解为异核双原子${C_{\infty \upsilon }}$点群表示的对应关系, Pu$\Sigma _{\rm{u}}^ + + {\Pi _{\rm{u}}}$, Fg$\Sigma _{\rm g}^ - \!+ {\Pi _{\rm{g}}} \!+ {\Delta _{\rm{g}}} \!+ {\Phi _{\rm{g}}}$, 通过$\Sigma _{\rm{u}}^{\rm{ + }} + {\Pi _{\rm{u}}}$$\Sigma _{\rm{g}}^ - + {\Pi _{\rm{g}}} + $${\Delta _{\rm{g}}} + {\Phi _{\rm{g}}} $的直积可推导出对应最低离解极限Al(2Pu)+Ti(3Fg), TiAl可能的电子态有${}^{2, 4}\Sigma _{\rm{u}}^{\rm{ + }} + $$ {}^{2, 4}\Sigma _{\rm{u}}^ - (2) + {}^{2, 4}{\Pi _{\rm{u}}}(3) \!+\! {}^{2, 4}{\Delta _{\rm{u}}}(2) + {}^{2, 4}{\Phi _{\rm{u}}}(2) \!+\! {}^{2, 4}{\Gamma _{\rm{u}}}$, 共11个二重态和11个四重态, 考虑态的简并性, 计算中需要各计算19个二重态和四重态的根. 显然, 要完整地获得如此众多的电子态的信息是困难的, 为此围绕基态和最低的若干激发态开展分析计算.

      采用密度泛函B3LYP/def2-TZVP和B3LYP/def2-QZVPP对二重态和四重态进行构型优化计算, 通过比较稳定构型下的单点能, 两种基组下得出的结论是一致的, 即TiAl的基态为四重态, 而非二重态. 两种方法得到的基态平衡核间距分别为0.27206 nm和0.27136 nm, 这与Ouyang等[21]用6-311+G(2d)(Al)+Lan12DZ(Ti)基组的计算值0.277 nm基本一致. 下面的计算围绕四重态进行.

    • Al原子有3个价电子(3s23p1), Ti原子有4个价电子(3d24s2), 这7个价电子共涉及10个原子轨道. 由群论分析可知, 在${C_{\infty \upsilon }}$点群下, Ti的4s, $(3{{\rm{d}}_{xz}}, 3{{\rm{d}}_{yz}})$, ($3{{\rm{d}}_{{x^2} - {y^2}}}, 3{{\rm{d}}_{xy}}$), $(3{{\rm{d}}_{{z^2}}})$ 分别构成${C_{\infty \upsilon }}$点群的${\Sigma ^ + }, \;\Pi,\;\Delta, \;{\Sigma ^ + }$不可约表示基函数, 即将Ti原子的4s, 3d轨道对称性匹配为$\text{σ}$(4s), ${\text{π}} (3{{\rm{d}}_{xz}}, $ $3{{\rm{d}}_{yz}}), \text{δ} (3{{\rm{d}}_{{x^2} - {y^2}}}, 3{{\rm{d}}_{xy}}), \sigma (3{{\rm{d}}_{{z^2}}})$轨道; 同样, Al的3s和3p轨道分别对称性匹配为${\rm{\text{σ} (3 s), {\text{π}} (3}}{{\rm{p}}_x}, {\rm{3}}{{\rm{p}}_y}), $ $\text{σ} (3{{\rm{p}}_z}) $轨道. 由对称性相同的轨道进行组合而得到TiAl的分子轨道. 由这10个价电子原子轨道构成的10个分子轨道为${\rm{\text{σ} (4) + {\text{π}} (2) + \text{δ} }}$, 其中${\rm{{\text{π}}, \text{δ} }}$分别为二重简并轨道.

      以B3LYP/def2-TZVP优化后的TiAl分子轨道为初始轨道猜测, 对TiAl基态进行CAS (7, 10)的尝试计算. 因为暂时不能确定基态的对称性, 即尚未对基态进行指认, 故计算时未选择对称性限制, 只设置为四重态而已. 在核间距R = 0.27206 nm下, 计算得到的活动基轨道分布如表1所列.

      MO No.14151617181920212223
      Energy/Eh–0.4035–0.18380.0029–0.0159–0.01590.07430.07430.04390.17140.1714
      Number of occupied electron1.9571.7720.6690.5950.5950.4940.4940.2220.0990.099
      Symbolσσσππδδσππ
      Tisσ12.147.57.7000013.400
      Tipzσ7.52.61.2000038.200
      Tipxπ0004.9000000
      Tipy0001.8000000
      Tidz2σ71.786.1000010.400
      Tidxzπ00031.411.400034.631.2
      Tidyz00011.431.400031.234.6
      Tidx2y2δ0000082.617.2000
      Tidxy0000017.282.6000
      Alsσ72.97.81.20000000
      Alpzσ0.539.23.2000036.600
      Alpxπ00034.412.500016.715.1
      Alpy00012.534.400015.116.7

      表 1  CAS (7, 10)/def2-TZVP计算的活动基分子轨道(MO14−MO23)系数(Eh = 2625.5 kJ/mol)

      Table 1.  Coefficients of the CAS orbital (MO14−MO23) calculated by CAS (7, 10)/def2-TZVP.

      表1可见, 这10个分子活动基轨道编号为MO14-MO23. 注意轨道的次序是按照轨道上电子占居数的多少排列的, 而不是按照轨道的能量高低. 这10个分子轨道恰好均主要由Al和Ti原子的10个价电子对应的原子轨道构成. 按照轨道能量由低到高的次序, 这10个分子轨道依次为σσππσσδδππ, 为方便讨论, 将其依次编号为(1σ)(2σ)(1π)(2π)(3σ)(4σ)(1δ)(2δ)(3π)(4π). 但需要指出的是: 1)${\rm{(1{\text{π}} )}}$${\rm{(2{\text{π}} )}}$为简并轨道, 类似的简并轨道还有(1δ)和(2δ), ${\rm{(3{\text{π}} )}}$${\rm{(4{\text{π}} )}}$; 2)在势能曲线的扫描计算过程中, 随着核间距的变化, 这10个分子轨道的能级次序可能发生变化, 需要注意观察每个核间距下分子轨道的排序.

      基于上面的在R = 0.27206 nm下CAS (7, 10)的波函数作为TiAl基态及激发态势能曲线扫描计算的初始值, 进行势能曲线的扫描计算. 扫描区间设定为R = 0.200—0.500 nm. 首先进行态平均自洽场计算, 即SA-CAS (7, 10)的扫描计算, 在此基础上进行强收缩-N电子价态微扰(SC-NEVPT2)的动态相关的修正计算和非收缩-多参考组态相互作用的修正计算. 但MRCI修正计算的势能曲线存在明显的跳变现象, 这里仅给出SC-NEVPT2修正计算的结果.

    • 经多次尝试, 分别在基组为def2-TZVP, def2-TZVPP, def2-QZVPP下完成了对TiAl最低6条四重态势能曲线的态平均CASSCF扫描计算, 即SA-CAS (7, 10)计算, 结果如图1所示, 结果可归结如下.

      图  1  SA-CAS (7, 10)对TiAl最低6条电子态势能曲线的扫描计算, 其中采用的基组分别是(a) def2-TZVP; (b) def2-TZVPP, (c) def2-QZVPP

      Figure 1.  The lowest 6 potential energy curves of TiAl calculated by SA-CAS (7, 10) with basis set of (a) def2-TZVP; (b) def2-TZVPP; (c) def2-QZVPP.

      1) 随着计算基组的增大, 势能曲线更为光滑. 在def2-TZVP基组下, 曲线在R = 0.37 nm附近出现明显的不正常起伏, 这是波函数不稳定造成的. 但基组增大到def2-TZVPP及以后, 函数的稳定性已能满足计算要求. Karlsruhe group提供的这套基组非常适合非相对论条件下对含过渡金属原子分子体系的计算, 而且def2-QZVPP给出最精确的能量计算值. 在本计算中, 基组def2-TZVP, def2-TZVPP, def2-QZVPP所用的收缩后的基函数数目分别为82, 106, 163. 为了得到相对可靠的计算结果, 后续的计算均采用def2-TZVPP基组或更大的def2-QZVPP基组.

      2) 势能曲线在R = 0.255 nm附近存在跳变. 在def2-TZVP基组下这种跳变极为明显, 随着基组的加大, 这种跳变仍无法消除. 怀疑在该核间距处存在电子轨道的跳变, 后面对活动基轨道能量及其轨道电子占居数随核间距的变化进行分析, 也说明这一推断.

      3) TiAl最低的3个四重态均为二重简并态. 计算的6条势能曲线均表现为每2条曲线重合的特征, 即二重简并的特征, 故其对应的电子态应为$\Pi, \Delta, \Phi $等, 而不会是$\Sigma $态. 具体电子态的指认通过电子组态及跃迁偶极矩加以分析.

      4) 6条势能曲线趋于同一离解极限. 尽管在$ R\!= $ 0.500 nm处各势能曲线并未完全重合, 但能量差仅为270 cm–1 (def2-TZVP), 220 cm–1 (def2-TZVPP), 200 cm–1 (def2-QZVPP).

      5) TiAl最低的三个四重态似乎均为束缚态, 但除了基态外, 另外2个激发态的束缚能极小. 以def2-QZVPP计算的结果进行分析, 三个态的束缚能分别大约为3000, 700, 600 cm–1. 更精确的值应由后面考虑动态相关计算后加以确定.

      6) 基态的平衡核间距均大于DFT优化给出的结果. 前面通过用密度泛函B3LYP/def2-QZVPP对基态构型优化得出平衡核间为0.271 nm, 而SA-CAS (7, 10)/def2-TZVPP和SA-CAS (7, 10)/def2-QZVPP势能曲线扫描计算给出的Re分别为0.285 nm和0.286 nm.

    • 为了分析图1中势能曲线在R = 0.255 nm附近的“突变”, 下面以图1(b)中def2-TZVPP基组下的计算结果为例, 从SA-CAS (7, 10)计算输出文件中逐点提取出10个活动基轨道在各核间距下的轨道能量分布次序以及各轨道上的电子占居数, 结果如图2图3所示.

      图  2  活动基轨道能量随核间距的变化

      Figure 2.  Variation of CAS orbital energy with nuclear distance.

      图  3  活动基轨道上电子占居数随核间距的变化

      Figure 3.  Variation of occupied electrons in the CAS orbital with nuclear distance.

      图2图3所示的结果可以看出, 在R = 0.255 nm前后, 轨道能级以及轨道上的电子占居数均发生“剧烈”变化.

      1) 在R = 0.200—0.255 nm之间: 除了能量接近的轨道${\rm{(1\text{σ} )}}$${\rm{(1{\text{π}} )(2{\text{π}} )}}$R = 0.215 nm处出现交叉外, 其余轨道能级次序均保持不变, 且电子占居数在该核间距范围内变化不大, 如在R = 0.210 nm和R = 0.240 nm处的电子占居数分布为:

      $ \begin{split} R = \,& 0.210\;{\rm{ nm}},{{\left( {1\text{σ}} \right)}^{1.93}}{{\left( {1{\text{π}}} \right)}^{1.24}}{{\left( {2{\text{π}}} \right)}^{1.24}}{{\left( {2\text{σ}} \right)}^{1.23}}{{\left( {1\text{δ}} \right)}^{0.50}}\\&{{\left( {2\text{δ}} \right)}^{0.50}}{{\left( {3\text{σ}} \right)}^{0.14}}{{\left( {3{\text{π}}} \right)}^{0.10}}{{\left( {4{\text{π}}} \right)}^{0.10}}{{\left( {4\text{σ}} \right)}^{0.02}};\\ R = \,& 0.240\;{\rm{ nm}},{{\left( {1\text{σ}} \right)}^{1.93}}{{\left( {2\text{σ}} \right)}^{1.23}}{{\left( {1{\text{π}}} \right)}^{1.18}}{{\left( {2{\text{π}}} \right)}^{1.18}}{{\left( {1\text{δ}} \right)}^{0.50}}\\&{{\left( {2\text{δ}} \right)}^{0.50}}{{\left( {3\text{σ}} \right)}^{0.14}}{{\left( {3{\text{π}}} \right)}^{0.15}}{{\left( {4{\text{π}}} \right)}^{0.15}}{{\left( {4\text{σ}} \right)}^{0.02}}. \end{split}$

      比较可发现, 电子占居数变化和轨道能级变化主要发生在轨道${\rm{(1{\text{π}} )(2{\text{π}} )}}$${\rm{(3{\text{π}} )(4{\text{π}} )}}$之间. ${\rm{(1{\text{π}} )(2{\text{π}} )}}$的轨道能随核间距增大在线性增加, 而${\rm{(3{\text{π}} )(4{\text{π}} )}}$的轨道在线性降低; 轨道${\rm{(1{\text{π}} )(2{\text{π}} )}}$上的电子占居数从1.24减小为1.18, 而轨道${\rm{(3{\text{π}} )(4{\text{π}} )}}$上的电子数却从0.10增加为0.15. 考虑到${\rm{(1{\text{π}} )(2{\text{π}} )}}$等为简并轨道, 上面的电子占居分布可简化或等效为:

      $R = {\rm{ }}0.210\;{\rm{ nm}}, {\left( \text{σ} \right)^{1.93}}{\left( {\text{π}} \right)^{2.5}}{\left(\text{σ} \right)^{1.23}}{\left(\text{δ} \right)^1}$, 近似为${\left(\text{σ} \right)^2}{\left( {\text{π}} \right)^3}{\left( \text{σ} \right)^1}{\left(\text{δ} \right)^1}$;

      $R = {\rm{ }}0.240\;{\rm{ nm}}, {\left( \text{σ}\right)^{1.93}}{\left( \text{σ} \right)^{1.23}}{\left( {\text{π}} \right)^{2.36}}{\left( \text{δ} \right)^1}$, 近似为${\left(\text{σ}\right)^2}{\left(\text{σ} \right)^1}{\left( {\text{π}} \right)^3}{\left(\text{δ}\right)^1}$.

      可见, 这里仅仅发生了第2个σ轨道与第1个π轨道的能级次序调换.

      为了更清晰地了解这两组π轨道随核间距的变化, 围绕“特殊”点R = 0.255 nm前后选取了4个特征点, 即R = 0.200, 0.240, 0.280, 0.490 nm, 从计算输出文件中提取出这4个点处这两组${\text{π}}$轨道的主要原子轨道构成, 如表2所列, 与前面理论分析一样, 这两组π 轨道均主要由${\rm{Ti(3}}{{\rm{d}}_{xz}}, 3{{\rm{d}}_{yz}}{\rm{)}}$${\rm{Al}}(3{{\rm{p}}_x}, 3{{\rm{p}}_y})$构成. 但随着核间距的变化, 原子轨道所占比例发生明显变化. 比如${\rm{(1{\text{π}} )(2{\text{π}} )}}$轨道, 在R = 0.200和0.240 nm处, ${\rm{Ti}}(3{{\rm{d}}_{xz}}, 3{{\rm{d}}_{yz}})$的占比近似为${\rm{Al}}(3{{\rm{p}}_x}, 3{{\rm{p}}_y})$的2倍. 而且无论${\rm{(1{\text{π}} )(2{\text{π}} )}}$还是${\rm{(3{\text{π}} )(4{\text{π}} )}}$, 轨道占比以${\rm{Ti}}(3{{\rm{d}}_{xz}}, 3{{\rm{d}}_{yz}})$为主.

      OrbitalR = 0.200 nmR = 0.240 nmR = 0.280 nmR = 0.490 nm
      ${\rm{(1{\text{π}})(2{\text{π}})}}$${aligned}& {\rm{Ti(}}3{{\rm{p}}_x}, {\rm{ }}3{{\rm{p}}_y}{\rm{) }}7{\rm{\% }} \\& {\rm{Ti(}}3{{\rm{d}}_{xz}}, 3{{\rm{d}}_{yz}}{\rm{) }}60{\rm{\% }} \\ &{\rm{Al(}}3{{\rm{p}}_x}, 3{{\rm{p}}_y}{\rm{) }}28{\rm{\% }}{aligned} $${aligned}& {\rm{Ti(3}}{{\rm{p}}_x}, {\rm{ 3}}{{\rm{p}}_y}{\rm{) }}7{\rm{\% }} \\& {\rm{Ti(3}}{{\rm{d}}_{xz}}, 3{{\rm{d}}_{yz}}{\rm{) }}5{\rm{7\% }} \\& {\rm{Al(3}}{{\rm{p}}_x}, 3{{\rm{p}}_y}{\rm{) }}3{\rm{2\% }}{aligned} $${aligned}& {\rm{Ti(3}}{{\rm{p}}_x}, {\rm{ }}3{{\rm{p}}_y}{\rm{) 3\% }} \\ &{\rm{Ti(3}}{{\rm{d}}_{xz}}, 3{{\rm{d}}_{yz}}{\rm{) 73\% }} \\& {\rm{Al(3}}{{\rm{p}}_x}, 3{{\rm{p}}_y}{\rm{) 21\% }}{aligned} $${\rm{Ti}}(3{{\rm{d}}_{xz}}, 3{{\rm{d}}_{yz}}){\rm{ }}1{\rm{00}}\% $
      ${\rm{(3{\text{π}})(4{\text{π}})}}$${aligned} &{\rm{Ti(3}}{{\rm{d}}_{xz}}{\rm{, 3}}{{\rm{d}}_{yz}}{\rm{) 52\% }} \\ &{\rm{Al(3}}{{\rm{p}}_x}{\rm{, 3}}{{\rm{p}}_y}{\rm{) 36\% }}{aligned} $${aligned} &{\rm{Ti(3}}{{\rm{d}}_{xz}}, 3{{\rm{d}}_{yz}}{\rm{) 52\% }} \\& {\rm{Al(3}}{{\rm{p}}_x}, 3{{\rm{p}}_y}{\rm{) 40\% }}{aligned} $${aligned}& {\rm{Ti(3}}{{\rm{p}}_x}, {\rm{ }}3{{\rm{p}}_y}{\rm{) 12\% }} \\ &{\rm{Ti(3}}{{\rm{d}}_{xz}}, 3{{\rm{d}}_{yz}}{\rm{) 34\% }} \\ &{\rm{Al(3}}{{\rm{p}}_x}, 3{{\rm{p}}_y}{\rm{) 52\% }}{aligned} $${\rm{Al(3}}{{\rm{p}}_x}, 3{{\rm{p}}_y}{\rm{) 99\% }}$

      表 2  两组${\text{π}}$轨道的组成分析

      Table 2.  Composition analysis of two π orbits

      2) 在R = 0.255—0.500 nm之间: 由图2可见, 轨道能级次序有多处交叉, 首先是0.255 nm处的${\rm{(3{\text{π}} )(4{\text{π}} )}}$${\rm{(1\text{δ})(2\text{δ} )}}$交叉, 然后是0.300 nm附近的${\rm{(3{\text{π}} )(4{\text{π}} )}}$分别与${\rm{(1{\text{π}} )(2{\text{π}} )}}$${\rm{(3\text{σ} )}}$交叉. 从图3的电子占居数分析, ${\rm{(2\text{σ} )}}$轨道占居数明显增加, 且逐渐接近电子双占居. 而${\rm{(1{\text{π}} )(2{\text{π}} )}}$电子占居数在R = 0.255 nm前后发生了急剧变化, 电子占居数几乎减半. 同时, ${\rm{(3{\text{π}} )(4{\text{π}} )}}$${\rm{(1{\text{π}} )(2{\text{π}} )}}$恰好相反, 电子占居数急剧增加, 相当于轨道${\rm{(1{\text{π}} )(2{\text{π}} )}}$上的电子有一半转移至${\rm{(3{\text{π}} )(4{\text{π}} )}}$轨道上. 总之, 在该核间距段, ${\rm{(1{\text{π}} )(2{\text{π}} )}}$, ${\rm{(3{\text{π}} )(4{\text{π}} )}}$, ${\rm{(1\text{δ})(2\text{δ} )}}$上的电子占居数差距变小, 均趋向于占居数0.5. 为便于观察比较, 下面列出R = 0.280 nm和R = 0.490 nm处的活动基分子轨道次序以及电子占居分布:

      $ \begin{split} R =\,&0.280\;{\rm{nm}},{\left( {1\text{σ} } \right)^{1.95}}{\left( {2\text{σ} } \right)^{1.68}}{\left( {1{\text{π}} } \right)^{0.69}}{\left( {2{\text{π}} } \right)^{0.66}}{\left( {3\text{σ} } \right)^{0.32}}\\&{\left( {3{\text{π}} } \right)^{0.41}}{\left( {4{\text{π}} } \right)^{0.39}}{\left( {1\text{δ} } \right)^{0.44}}{\left( {\rm 2d} \right)^{0.44}}{\left( {4\text{σ} } \right)^{0.02}};\\ R =\, &0.490\;{\rm{ nm}},{\left( {1\text{σ} } \right)^{1.92}}{\left( {2\text{σ} } \right)^{1.90}}{\left( {3{\text{π}} } \right)^{0.52}}{\left( {4{\text{π}} } \right)^{0.52}}{\left( {1{\text{π}} } \right)^{0.62}}\\&{\left( {2{\text{π}} } \right)^{0.56}}{\left( {3\text{σ} } \right)^{0.10}}{\left( {1\text{δ} } \right)^{0.41}}{\left( {2\text{δ} } \right)^{0.41}}{\left( {4\text{σ} } \right)^{0.04}}. \end{split} $

      同样, 考虑到$(1{\text{π}} )(2{\text{π}} )$等为简并轨道, 上面的电子占居分布可简化或等效为:

      $R \!=\!0.280\;{\rm{ nm}}, {\left( \text{σ} \right)^{1.95}}{\left( \text{σ} \right)^{1.68}}{\left( {\text{π}} \right)^{1.35}}{\left( \text{σ} \right)^{0.32}}{\left( {\text{π}} \right)^{0.80}}{\left( \text{δ} \right)^{0.88}}$, 近似为${\left( \text{σ} \right)^2}{\left( \text{σ} \right)^2}{\left( {\text{π}} \right)^1}{\left( {\text{π}} \right)^1}{\left( \text{δ}\right)^1};$

      $R \!=\!0.490\;{\rm{ nm}}, {\left( \text{σ} \right)^{1.92}}{\left( \text{σ} \right)^{1.90}}{\left( {\text{π}} \right)^{1.0}}{\left( {\text{π}} \right)^{1.18}}{\left( {1\text{δ} } \right)^{0.82}}$, 近似为${\left( \text{σ} \right)^2}{\left( \text{σ} \right)^2}{\left( {\text{π}} \right)^1}{\left( {\text{π}} \right)^1}{\left( \text{δ} \right)^1}$.

      分析R = 0.490 nm处的活动基轨道组成, 见表3, 其中编号14—23为活动基分子轨道, 主要电子占居分布可近似表示为(3sAl)1.92 (4sTi)1.90 (3dTi)0.62(3dTi)0.56 (3pAl)0.52 (3pAl)0.52 (3dTi)0.41 (3dTi)0.41, 即(3sAl)1.92 (3pAl)0.93, (4sTi)1.90 (3dTi)2.00, 近似为(3sAl)2(3pAl)1, (4sTi)2 (3dTi)2, 可见, 此时TiAl基本分解为两个中性原子, 这也恰好符合前面对TiAl最低离解的分析.

      Orbital No121314151617
      Energy/Hartree–1.79778–1.7976–0.37486–0.215130.020010.03822
      Occupied electron2.000002.00001.919761.899000.624630.55980
      Ti s002.894.500
      Ti pz099.80.5000
      Ti px55.70.10000
      Ti py44.300000
      Ti dxz000055.843.6
      Ti dyz000044.155.2
      Al s0095.82.800
      Orbital No181920212223
      Energy/Hartree–0.00791–0.006820.08850.089790.056510.1227
      Occupied electron0.517860.517500.410580.406450.101590.04283
      Ti pz0000.091.45.3
      Ti dx2y2001.498.600
      Ti dxy0098.41.400
      Al pz0000792.7
      Al px55.343.30000
      Al py43.754.80000

      表 3  R = 0.490 nm处活动基分子轨道MO14-MO23组成

      Table 3.  Composition of CAS orbitals MO14-MO23 at R = 0.490 nm.

      综上可知, 发生在R = 0.255 nm附近的“突变”现象, 可大致归结为电子占居轨道的突变, 即由${\left( \text{σ} \right)^2}{\left( {\text{π}} \right)^3}{\left( \text{σ} \right)^1}{\left( \text{δ} \right)^1}$跳变为${\left( \text{σ} \right)^2}{\left( \text{σ} \right)^2}{\left( {\text{π}} \right)^1}{\left( {\text{π}} \right)^1}{\left( \text{δ} \right)^1}$, 因而R = 0.255 nm两边的能态是不同的, 下面对能态进行分析标识.

    • 图1(b)图1(c)的势能曲线可得出, 在“突变”点R = 0.255 nm前后的电子态是不同的, 为此, 以def2-TZVPP基组下的计算结果为例, 即从图1(b)所示的计算输出文件中分别提取出R = 0.240 nm和R = 0.285 nm时各电子态的相应的组态以及由基态向激发态的跃迁偶极矩, 结果如表4所列.

      R/nmstateMain configurationExcitation energy/cm–1Transition dipole moment T2/Debye2Possible quartet stateIdetified state
      0.285Ground state${{\rm{\text{σ} }}^{\rm{2}}}{{\rm{\text{σ} }}^{\rm{2}}}{{\text{π}}^{\rm{2}}}{{\rm{\text{δ} }}^{\rm{1}}}{{\text{π}}^{\rm{0}}}$04ΔX4Δ
      1st excited state${{\rm{\text{σ} }}^{\rm{2}}}{{\rm{\text{σ} }}^{\rm{2}}}{{\text{π}}^{\rm{2}}}{{\rm{\text{δ} }}^{\rm{0}}}{{\text{π}}^{\rm{1}}}$32120.0344ΠA4Π
      2nd excited state${{\rm{\text{σ} }}^{\rm{2}}}{{\rm{\text{σ} }}^{\rm{2}}}{{\text{π}}^{\rm{1}}}{{\rm{\text{δ} }}^{\rm{1}}}{{\text{π}}^{\rm{1}}}$346204Σ, 4Δ(2), 4ΓB4Γ
      0.240Ground state${{\rm{\text{σ} }}^{\rm{2}}}{{\rm{\text{σ} }}^{\rm{2}}}{{\text{π}}^{\rm{2}}}{{\rm{\text{δ} }}^{\rm{1}}}$04ΔX4Δ
      1st excited state${{\rm{\text{σ} }}^{\rm{2}}}{{\rm{\text{σ} }}^{\rm{1}}}{{\text{π}}^{\rm{3}}}{{\rm{\text{δ} }}^{\rm{1}}}$41400.008244Π, 4ΦA'4Φ
      2nd excited state${{\rm{\text{σ} }}^{\rm{2}}}{{\rm{\text{σ} }}^{\rm{1}}}{{\text{π}}^{\rm{3}}}{{\rm{\text{δ} }}^{\rm{1}}}$47270.008694Π, 4ΦB'4Π
      3rd excited state${{\rm{\text{σ} }}^{\rm{2}}}{{\rm{\text{σ} }}^{\rm{1}}}{{\text{π}}^{\rm{3}}}{{\rm{\text{δ} }}^{\rm{1}}}$50740.005514Π, 4ΦB'4Π

      表 4  基态及最低激发态的组态及跃迁偶极矩

      Table 4.  Configuration and transition dipole moment of the ground state and the lowest excited state

      1) R = 0.285 nm

      由基电子组态${{\rm{\text{σ} }}^{\rm{2}}}{{\rm{\text{σ} }}^{\rm{2}}}{{\text{π}}^{\rm{2}}}{{\rm{\text{δ} }}^{\rm{1}}}{{\text{π}}^{\rm{0}}}$可以推算出可能的四重态只有一个, 即4Δ, 即基态为X4Δ. 对于第一激发态的组态${{\rm{\text{σ} }}^{\rm{2}}}{{\rm{\text{σ} }}^{\rm{2}}}{{\text{π}}^{\rm{2}}}{{\rm{\text{δ} }}^{\rm{0}}}{{\text{π}}^{\rm{1}}}$, 可能的四重态也只有一个, 即4Π, 将其标记为A4Π. 对于第二激发态的组态${{\rm{\text{σ} }}^{\rm{2}}}{{\rm{\text{σ} }}^{\rm{2}}}{{\text{π}}^{\rm{1}}}{{\rm{\text{δ} }}^{\rm{1}}}{{\text{π}}^{\rm{1}}}$, 可能的四重态有四个, 即4Σ, 4Δ(2), 4Γ, 考虑到势能曲线已经表明该态为二重简并态, 按照Hund定则, 4Γ的能量应该最低, 故第二激发态应该为4Γ态, 且基态X4Δ到4Γ态的跃迁不符合跃迁规则, 这也符合表4中该跃迁偶极矩为0的计算结果, 将其标识为B 4Γ. 图4给出了SA-CAS (7, 10)计算下的标识.

      图  4  TiAl最低四重电子态的标识(SA-CAS计算结果)

      Figure 4.  Identification of the lowest quadruple electronic state of TiAl (SA-CAS calculation results).

      2) R = 0.240 nm

      基态的电子组态和R = 0.285 nm时一样, 故基态为X4Δ. 对于简并的第一激发态, 其主要组态为${{\rm{\text{σ} }}^{\rm{2}}}{{\rm{\text{σ} }}^{\rm{1}}}{{\text{π}}^{\rm{3}}}{{\rm{\text{δ} }}^{\rm{1}}}$, 可能的四重态只有2个, 即4Π, 4Φ. 按照Hund定则, 4Φ的能量应该最低, 故该激发态应该为4Φ态, 且基态X4Δ到4Φ态的跃迁符合跃迁规则, 这也符合表4中该跃迁偶极矩不为0的计算结果, 将其标识为A'4Φ. 对于第二和第三激发态, 虽然势能曲线没有能够重合, 但它们的主要电子组态和第一激发态的一样, 均为${{\rm{\text{σ} }}^{\rm{2}}}{{\rm{\text{σ} }}^{\rm{1}}}{{\text{π}}^{\rm{3}}}{{\rm{\text{δ} }}^{\rm{1}}}{\rm A}$, 可能的四重态也只能为4Π和4Φ, 且基态X4Δ到它们的跃迁偶极矩也不为0, 将其标识为B'4Π. 同时注意到, 基态X4Δ到A'4Φ和B'4Π的跃迁偶极矩远小于X4Δ到A4Π的值.

    • 基于上面的CAS计算及其态的标识, 可将TiAl最低的三个四重态明确标识出来, 并在CAS(7, 10)计算的基础上, 进行了SC-NEVPT2动态相关修正计算, 图5给出了在Karlsruhe提供的最大基组def2-QZVPP下的计算结果. 为方便显示, 取基态势能曲线的最小值为能量参考点, 势能的相对能量用波数表示.

      图  5  TiAl电子态势能曲线的SC-NEVPT2动态相关修正计算

      Figure 5.  Dynamic correlation correction calculation of SC-NEVPT2 for TiAl electronic potential energy curve.

      SC-NEVPT2修正后的曲线与SA-CAS (7, 10)计算相比, 有以下4点不同.

      1) 激发态出现明显双势阱. 对于激发态A4Π和B4Γ而言, CAS计算显示只在“突变”点R = 0.255 nm右方, 即R = 0.26—0.50 nm之间出现势阱, 而SC-NEVPT2修正后, 在该点前后均存在极小值, 即出现双势阱, 两势阱之间的势垒高度大约为2700 cm–1.

      2) 平衡核间距明显缩小. CAS计算和NEVPT2修正计算给出基态X4Δ的平衡核间距分别为0.288 nm和0.266 nm, 即NEVPT2修正后核间距明显减小了. 激发态的情形也类似, 如突变点右方的势阱, A4Π的Re由0.320 nm (CAS)减小为0.296 nm (NEVPT2), B4Γ态的相应变化为0.324 nm (CAS)和0.306 nm (NEVPT2).

      3) 束缚能增大. 以R = 0.500 nm处的能量作为参照, 基态X4Δ的束缚能由3016 cm–1 (CAS)增大为8151 cm–1 (NEVPT2). 激发态A4Π和B4Γ态的势阱深度也相应增加, 具体数值见表5.

      StateRe/nmDe/cm–1
      CASNEVPT2CASNEVPT2
      X4Δ0.2880.26630168151
      A4Π0.320$\left\{{aligned}& {0.248} \\ & {0.296} {aligned} \right.$796$\left\{{aligned}& {3845} \\ & {3406} {aligned} \right.$
      B4Γ0.324$\left\{{aligned}& {0.248} \\ & {0.306} {aligned} \right.$711$\left\{{aligned}& {2884} \\ & {3406} {aligned} \right.$

      表 5  TiAl最低3个四重态的结构参数

      Table 5.  Structural parameters of the lowest three quadruple states of TiAl.

      4) 在R < 0.255 nm区域, A'4Φ和B'4Π态的简并特征趋向于退简并. 这可能是由于核间距的缩小, 电子自旋与电子轨道的耦合得以加强, 使得简并趋于消除.

    • 如果图5所示的势能曲线最为接近真实情况, 那么可以根据这一计算结果来解释为什么实验上没有能够观测到TiAl的光谱. TiAl基态X4Δ的平衡核间距为Re = 0.266 nm (考虑NEVPT2修正), 而结构的“突变”点在0.255 nm附近, 非常靠近Re. 因而TiAl由基态的光跃迁可能存在不稳定, 或者吸收系数很小. 同时, 基态X4Δ与共振激发态A4Π的平衡核间距差距也比较大, 即分别为0.266 nm和0.296 nm, 因而Franck-Conden因子也应该偏小. 这些因素最终决定了实验上很难观测到TiAl的光谱. 同时从电子态的束缚能分析(仅依据NEVPT2修正后的计算结果), 如表5所列, TiAl基态的束缚能De 仅为8151 cm–1 (约1.01 eV). TiAl由基态X4Δ (Re = 0.266 nm)到共振激发态A4Π的垂直激发能和绝热激发能分别为6001 cm–1 (0.74 eV)和4908 cm–1 (0.61 eV), 而A4Π的束缚能仅为3406 cm–1 (0.42 eV). 如果TiAl被垂直激发到激发态A4Π, 其过剩的动能很容易跨越激发态A4Π的浅势阱, 从而导致TiAl离解. TiAl电子结构上的这一特殊情形, 也可为理解TiAl合金表现出的脆性提供参考.

    • 采用CASSCF+NEVPT2方法, 在CAS (7, 10)空间下完成了TiAl基态及最低两个激发态的势能曲线扫描计算. 这几个态均为束缚态, 且完全趋向于中性原子Ti+Al的离解极限. 根据电子组态特征、跃迁偶极矩大小等信息, 对这几个电子态进行了标识, 得出基态为4Δ态, 并发现R = 0.255 nm附近可能存在一个电子态结构的“突变”点. 由于电子在分子轨道中的布局在该点处发生跳变, 导致该点前后的激发态不同. 在R > 0.255 nm区域, 基态和两个激发态分别为X4Δ, A4Π和B4Γ; 在R < 0.255 nm区域, 三个电子态变为X4Δ, A'4Φ和B'4Π, 且存在激发态简并解除的现象. 通过对电子态平衡核间距、激发能、束缚能等特征参数的分析, 得出实验上观测不到TiAl电子跃迁光谱可能与激发态结构存在“突变”有关, 且由于Franck-Conden因子偏小、激发态势阱太浅, 均容易导致TiAl被激发后容易发生解离.

参考文献 (29)

目录

    /

    返回文章
    返回