过渡金属类立方烷原子簇化合物电子结构和二阶非线性光学性质的密度泛函理论研究

采用基于第一原理的含时密度泛函理论(TDDFT)对一系列具有类立方烷簇芯结构的过渡金属簇合物二阶非线性光学性质进行了研究。结果证明, 由于簇芯结构的对称性的影响, 这一类簇分子的二阶非线性光学系数的数值要小于三核欠完整类立方烷体系。通过对电子结构的分析, 发现二阶非线性光学性质主要是由簇芯内电荷的迁移造成的。轨道分析显示, -S原子对于电荷的迁移起主要的传递作用。定域化轨道分析证明簇分子中存在的多中心键有利于簇芯内电荷的迁移。分子模拟的研究表明:虽然类立方烷结构簇分子的值比较小, 但是通过合理的配体设计, 获得具有较大值非线性光学晶体是可能的。

三核欠完整立方烷金属原子簇二阶非线性光学性质的理论研究

研究了具有欠完整立方烷构型的过渡金属原子簇分子的二阶非线性光学性质。利用TDDFT方法计算了选取的簇分子及相应模拟构造分子的静态和动态的一阶非线性光学超级化率(ijk);并计算了不同金属、桥原子和配体以及簇芯对该类化合物一阶超级化率的影响。选取其中的一个簇分子为基本模型,分析了该分子的电子结构和分子轨道,在微观水平上阐述了其非线性光学性质的可能起源。认为由过渡金属和硫原子组成的簇芯和与桥原子相连的配体对该类簇合物的二阶非线性光学性质的起决定性的作用。

某些过渡金属钌配合物非线性光学性能的从头算和密度泛函理论研究

实验发现金属钌的某些配合物具有大的非线性光学性能 .我们的理论研究表明金属钌配合物[Ru(NH3 ) 4LDLA]n + (n =2 ,3;LD,LA=吡啶衍生物配体 )的非线性光学性能取决于推电子基团LD 的给电子能力以及拉电子基团LA 的受电子能力 ,LD 和LA 的推拉电子能力越强越有利于提高配合物的二阶非线性光学系数 β,因此带正电荷的LA 能大幅度提高 β值 .虽然增加共轭体系的长度有利于提高 β值 ,但在Ru的配体中 ,吡啶环间或吡啶环与苯环间不一定要保持共面也会有大的β值 .DFT和abinitio方法的计算结果对比表明 ,对于含过渡金属Ru的化合物 ,在HF水平上难以得到满意的结论 ,由从头算有限场方法计算得到的 β值偏小 ,而用TDDFT方法能得到可与实验值符合较好的结果 .

一水甲酸锂晶体三阶非线性光学性质理论研究

用含时耦合微扰Hartree Fork (CPHF)方法 ,有限场MP2方法以及有向气体方法 ,首次分别在分子和晶体水平上计算了一水甲酸锂晶体三次谐波极化率 .表明超分子γ( -3ω ;ω ,ω ,ω)数值与超分子链长有线性关系 ;晶体χ( -3ω ;ω ,ω ,ω) 计算值主要受到电子相关影响 ,其次受到基组和局域场影响 ,在非共振条件下 ,频率色散的影响甚小 .最后 。

一水甲酸锂晶体二阶非线性光学系数的理论计算

采用有限场方法和含时耦合微扰方法,在6-31+G*基组水平上并考虑了电子相关效应和色散效应的影响,计算了一水甲酸锂晶体的宏观二阶非线性光学极化率和二阶非线性光学系数。其中非线性光学系数d31和d32与Roberts报道的相符,而另一个系数d33则比文献值大。计算结果还表明,在非共振情况下,电子相关效应对非线性光学极化率的影响远远超过了色散效应,同时表明Roberts报道的非线性光学系数比Singh等人报道的更为合理。

过渡金属偶氮苯衍生物二阶非线性光学性质的密度泛函理论研究

采用含时密度泛函理论方法(TDDFT)计算了IB族过渡金属及第5周期的金属杂化偶氮苯生色团的二阶非线性光学极化率。研究了金属杂化偶氮苯生色团的电子激发跃迁。结果表明, 与金属离子络合后的杂化偶氮苯生色团, 二阶非线性光学极化率明显改变, 是金属的推拉电子效应的结果。IB族金属的在原有机偶氮苯共轭体系的电荷转移方向的强拉电子效应, 增大了电荷转移范围, 使二阶非线性光学极化率增大。与Nb络合的杂化生色团, 沿原电荷转移相反的方向拉动电荷, 电荷转移跃迁局限于金属离子附近, 未能产生大的激发跃迁偶极矩, 体系的二阶非线性光学响应因而降低。激发跃迁能量是另一个主要影响因素, 同一类跃迁中, 金属杂化生色团的跃迁能量越低, 其二阶非线性光学响应就越大。络合Rh的生色团与Nb杂化生色团类似, 电荷转移跃迁范围都比较小, 同时由于Rh与羧基结合时几乎垂直羧基平面, 电荷同时沿X, Y方向转移跃迁, 属于二维电荷转移类型。

无机-有机杂化非线性光学晶体材料LAP和d-LAP晶体红外吸收边计算

根据双原子谐振子模型近似,提出了估算分子晶体材料红外吸收边的理论方法,然后利用超分子量子化学从头算,计算了L-精氨酸一水磷酸盐(LAP)晶体和氘化L-精氨酸一水磷酸盐(d-LAP)晶体单元超分子的红外振动光谱,其中LAP晶体超分子计算值与晶体红外光谱实验值吻合,表明超分子计算能有效地模拟无机-有机杂化非线性光学晶体的红外振动光谱。在归属了超分子重要基频线后,分析其红外强度,并估算了泛频频率。根据我们提出的方法,估算了这两种晶体的红外吸收边,结果与实验值较吻合。表明我们建议的理论方法是合理的。此外,通过计算表明,如果非线性光学晶体材料的红外吸收边是由与活泼H有关的伸缩振动泛频频率决定,那么活泼H的氘化是一种改善红外吸收边的有效途径。

尿素分子团簇线性和一阶非线性光学极化率的理论计算

检验分子晶体非线性光学极化率的可加性，对于确认有向气体模型(Oriented-Gas Model) 具有重要的科学意义，对探索新型实用的分子非线性光学材料也有实际的意义。本文对尿素晶体中分子簇非线性光学极化率的可加性做了详尽的理论研究。首次做到了尿素分子最近邻排布团簇的非线性光学极化率的从头算，从而扩展了前人的工作。研究结果有效地确认了有向气体模型。同时，我们还提出了壳层结构模型 (Shell-Structure Model)，用以解决有向气体模型中的缺陷，可实现更为有效可靠地计算晶体宏观非线性光学系数。

LAP晶体二阶NLO响应电子起源理论研究

The properties of excited states and second-order nonlinear optics for supramolecule [(H2N)2+C- NH(CH2)3CH(NH3)+COO-H2PO4-H2O] from L-arginine phosphate monohydrate (LAP) crystals were calculated by using the CIS and finite-field RHF methods, respectively. The results revealed that the second-order nonlinear optical (NLO) response for this crystal originates largely from the contributions of intramolecular charge transfer of arginine, and the intermolecular charge transfer between arginine and H2PO4- plays an important role as well. The conclusion can help to understand the origin of NLO response in the inorganic-organic hybrid materials.

Ir-C_(60)配合物衍生结构三阶非线性光学性质的DFT计算研究

基于密度泛函理论(DFT),利用高精度量子化学计算方法对配位铱原子的富勒烯衍生物分子Ir(η5 C60Me5)(CO)2的非线性光学性质进行了计算研究,结果显示具有较高对称性的Ir(η5 C60Me5)(CO)2衍生物分子,具有较高的三阶非线性光学响应.为了揭示Ir-C60的非线性光学响应的起源,研究了Ir(η5 C60Me5)(CO)2的电子结构性质,并且基于含时密度泛函理论(TDDFT)研究了Ir(η5 C60Me5)(CO)2配合物的激发态性质,探讨了富d电子及d空轨道的重过渡金属铱与富π轨道及离域电子的富勒烯相互系统作用对体系三阶非线性光学性质的影响机制.

The First Principles Investigations of the Thermoelectric Properties of GaN with p-and n-Type Doping

We investigate the thermoelectric properties of GaN with p- and n-type doping by the first principles calculation and the semi-classical Boltzmann theory. We find that the power factors （Sacr） of p-type GaN （-3500 μW/mK2） is about twice that of the n-type （-1750 μW/mK2）, which indicates the thermoelectric properties of p-type GaN would be better. Thermal conductivity of GaN crystal decreases rapidly as the temperature increases, but it is still too large for thermoelectric applications. The figure of merit （ZT） estimated at 1500 K is 0.134 for p-type GaN crystal and 0.062 for the n-type.

Noble-metal-free plasmonic MoO_(3-x)-based S-scheme heterojunction for photocatalytic dehydrogenation of benzyl alcohol to storable H2 fuel and benzaldehyde

Simultaneous generation of H_(2) fuel and value-added chemicals has attracted increasing attention since the photogenerated electrons and holes can be both employed to convert solar light into chemical energy.Herein,for realizing UV-visible-NIR light driven dehydrogenation of benzyl alcohol(BA)into benzaldehydes(BAD)and H_(2),a novel localized surface plasmon resonance(LSPR)enhanced S-scheme heterojunction was designed by combining noble-metal-free plasmon MoO_(3-x) as oxidation semiconductor and Zn_(0.1)Cd_(0.9)S as reduction semiconductor.The photoredox system of Zn_(0.1)Cd_(0.9)S/MoO_(3-x) displayed an unconventional reaction model,in which the BA served as both electron donor and acceptor.The S-scheme charge transfer mechanism induced by the formed internal electric field enhanced the redox ability of charge carriers thermodynamically and boosted charge separation kinetically.Moreover,due to the LSPR effect of MoO_(3-x) nanosheets,Zn_(0.1)Cd_(0.9)S/MoO_(3-x) photocatalysts exhibited strong absorption in the region of full solar spectrum.Therefore,the Zn_(0.1)Cd_(0.9)S/MoO_(3-x) composite generated H_(2) and BAD simultaneously via selective oxidation of BA with high production(34.38 and 33.83 mmol×g^(–1) for H_(2) and BAD,respectively)upon full solar illumination.Even under NIR light irradiation,the H_(2) production rate could up to 94.5 mmol×g^(–1)×h^(–1).In addition,the Zn_(0.1)Cd_(0.9)S/MoO_(3-x) composite displayed effective photocatalytic H_(2) evolution rate up to 149.2 mmol×g^(–1)×h^(–1) from water,which was approximate 6 times that of pure Zn_(0.1)Cd_(0.9)S.This work provides a reference for rational design of plasmonic S-scheme heterojunction photocatalysts for coproduction of high-value chemicals and solar fuel production.

Theoretical Investigations of Nonlinear Optical Properties of Transition Metal Cluster Anions

In the framework of density functional theory （DFT）, the electronic excitations and nonlinear optical （NLO） properties of six binuclear transition metal cluster anions with the formula of [Ch2M-（μ-Ch）2-M＇CN]^2- （M = Mo, W; Ch = S, Se; M＇ = Cu, Ag） have been systemically investigated at both cases of gas phase and DMF solution. The obtained electronic absorption spectra reveal that the element replacements of metals M and ligands Ch have significant influence on the absorptions, especially on the low-lying ones. In addition, the transitions of μ-Ch→M are dominant for the low-lying excitations, whereas the transitions of M＇→M as well as Ch→M are mainly responsible for the higher excitations. The calculated molecular first and second hyperpolarizabilities present the remarkable element substitution and solvent effects. The analyses show that the transitions involving μ-Ch→M charge transfer make the critical contributions to the first hyperpolarizability t, and that the charge transfers from the moieties of MCh4 to M＇CN as well as those of μ-Ch→M and M＇→M are responsible for the second hyperpolarizability y. Moreover, the introduction of solvent leads to the results that the transitions within the moieties of MCh4 and M＇CN make larger contributions to the hyperpolarizability, especially to γ.

Theoretical Studies on the First Hyperpolarizabilities of One-dimensional Donor-bridge-acceptor Chromophores and New Applications of BLA in Determining Molecular First Hyperpolarizabilities

We present a quantum-chemical analysis of the relationship between the bond length alteration （BLA） and the static first hyperpolarizability of a series of one-dimensional （1D） chromophores with donor-bridge-acceptor （D-B-A） structures. The calculated results show that the parameter BLA can be considered as an indicator to evaluate the molecular first hyper- polarizability. Along the direction of molecular ground-state dipole moments, the evolutions of BLA can be classified into three categories： the first is a non-monotonic line, which represents most chromophores; the second is monotonic increasing; and the third, contrarily, is monotonic decreasing. On the whole, the first hyperpolarizabilities of these studied chromophores are the monotonic functions of BLA along the direction of dipole moments. Therefore, the first hyperpolarizability of these 1D chromophores can be preliminarily evaluated in terms of the development of BLA without a rigorous computation. In other words, one can roughly estimate the relative magnitude of the first hyperpolarizability according to the optimized geometry.

Theoretical Studies on the First Hyperpolarizabilities of One-dimensional Donor-bridge-acceptor Chromophores

The ground-state dipole moments and second-order nonlinear optical （NLO） properties of a series of one-dimensional （1D） chromophores with donor-bridge-acceptor （D-B-A） structures have been investigated by using the second-order MФller-Plesset （MP2） and density functional theory （DFT） methods with the basis set of 6-31＋G（d）. According to the calculated results, the relationship between the molecular static first hyperpolarizability （βμ） and the directions of electron transition has been summarized. In terms of the sign of βμ, these 1D organic chromophores were classified into two categories： type Ⅰ with negative βμ and type Ⅱ bearing positive βμ. The analyses show that the remarkable difference of the first hyperpolarizabilities between Ⅰ and Ⅱ chromophores is associated mainly with the electrostatic interaction between terminal groups and the transport electrons in excited states. Moreover, different from the popular viewpoint, the obtained results also show that most of this series of 1D D-B-A molecules are more charge-separated in the ground states than in the excited states. As a whole, this theoretical investigation, to some extent, can be considered as a useful reference in designing the NLO chromophores with large first hyperpolarizabilities.

Theoretical Studies on the Binding of Cd^(2+) to Uracil and Its Thio-derivatives

The interaction of Cd2＋ with uracil,2-thiouracil,4-thiouracil and 2,4-dithiouracil have been investigated by the density functional theory（DFT） calculations.For uracil and 2,4-dithiouracil,where the two basic sites are the same,Cd2＋ attachment to the heteroatom at position 4 is preferred.However,for the systems where both types of basic centers,a carbonyl and a thiocarbonyl groups,are present,Cd2＋ association with sulfur is favorable.The enhanced stability of these enolic and thiol forms comes from Cd2＋ interaction with two basic sites simultaneously,which thereby triggers a significant aromatization of the ring.More significantly,the Cd2＋ binding energy with uracil and its thio-derivatives is larger than the tautomerization barriers connecting the diketo-like forms with the corresponding enolic-like tautomers.Consequently,when associated with Cd2＋,all tautomers are energetically accessible and should be observed in the gas phase.

Theoretical Study of Hydrated Cd^(2+) Interactions with Guanine

Theoretical study was performed to investigate how the hydration of cadmium ca-tion influences the structure and properties of guanine.The aqueous environment was simulated by both explicit solvent（1-5 water molecules） model and implicit solvent model.For complexes in which Cd2＋ attached to the N（7） and O（6） sites of guanine,energy analysis together with the Natural Bonding Orbital（NBO） analysis were performed to elucidate the bonding characteristics in detail.The most stable structures are penta-coordinate complexes without aqua ligand located at the guanine site.Higher number of water ligands corresponds to higher stabilization energies.Average bonding energies of G-Cd increase with the number of water molecules.Bonding energies of water ligands depend on its position in the complexes.The charge distribution of guanine changed with increasing the number of water ligands,which may also influence the base-pairing pattern of guanine.There is positive charge transfer from guanine to aqua ligand as the number of the hydration waters increases.IEFPCM optimization has results comparable to the [CdG（H2O）5]2＋ structure 5a.

Ligand-controlled First Hyperpolarizabilities of a Series of Tetrahedral Iridium Clusters Ir_4(CO)_9L: a TDDFT Study

A series of tetrahedral iridium carbonyl clusters coordinated by systematically varied series of ligands have been studied by TDDFT method focusing on their electronic and non- linear optical properties. The clusters of Ir4（CO）12 （1）, Ir4（μ-CO）3（CO）9 （2）, Ir4（μ-L）（CO）10 （L = dppm 3, dppe 4, （Ph2P）2CHMe 5, Ph2P（CH2）3PPh2 6） and Ir4（CO）10（phen） （phen = 1,10-phen- anthroline） （7） exhibit the first static hyperpolarizabilities of medium magnitude （βtot-10×10^-30 esu）. The second order nonlinear optical response of the seven clusters increase from 0 to 23 ×10^-30 esu; the high symmetric cluster Ir4（CO）12 debases its symmetry and presents the second order nonlinear optical behavior as the coordination style of some carbonyls changes to bridge style, and then the response increases regularly with the systematical variation of the ligands. The origination of the first hyperpolarizability is discussed by the expanded orbital decomposition scheme. The results suggest the d-d electron transition from the apical iridium atom to the other three Ir atoms inside the metal skeleton, and d-πelectron transitions from metals to carbonyls are responsible for the first hyperpolarizabilities. Particularly, for cluster 7, the charge transfer from d orbitals of iridium to π＊ orbirals of phenanthroline originates the first hyperpolarizabilities.

TDDFT Study on the Electronic Excitations and First Hyperpolarizabilities of Mixed-metal Carbonyl Clusters W_2Ir_2(μ-L)(CO)_8(η～5-C_5H_4Me)_2(L=dppe/dppf)

Mixed-metal carbonyl clusters of W2Ir2（CO）10（η^5-C5H4Me）2 1 and W2Ir2（μ-L）（CO）8（η^5-C5H4Me）2 （L = dppe 2, dppf 3） have been studied by TDDFT method focusing on their electronic and nonlinear optical properties. These three clusters exhibit the first static hyperpolarizabilities of medium magnitude （βtot-10×10^-30 esu）. The origin of β is discussed by the new proposed orbital-pair decomposition scheme by Barandes et al. The result suggests that the β values of the two clusters are mainly originated from d-d electron transition within the metal skeleton, and d-p （π＊） electron transition from metals to carbonyls and phenyl. The additional coordination by the electron donor group, ferrocene, makes cluster 3 own much larger β values, and the relatively longer range charge transfer from d orbitals of ferrocene to d orbirals of Ir and W is responsible for the enhanced β values.

Ultra-incompressible Ternary Diborides Re_(0.5)Ir_(0.5)B_2,Re_(0.5)Tc_(0.5)B_2,Os_(0.5)W_(0.5)B_2 and Os_(0.5)Ru_(0.5)B_2:Density Functional Calculation

In this paper, density functional computations have been applied to the structural, elastic and electronic properties of ternary transition metal diborides Re0.5Ir0.5B2, Re0.5Tc0.5B2, Os0.5W0.5B2 and Os0.5Ru0.5B2 in hexagonal （P63/mmc） and orthorhombic （Pmmn） structures with both local density approximation and generalized gradient approximation. LDA gives smaller lattice parameters and larger elastic moduli than GGA. Both results show that the hexagonal ones are more stable than orthorhombic ones except Os0.5Ru0.5B2. Moreover, the hexagonal structure has superior elastic property than orthorhombic one. Generally speaking, the calculated elastic moduli of Re0.5Ir0.5B2 and Os0.5Ru0.5B2 are smaller than those values of Re0.5Tc0.5B2 and Os0.5W0.5B2 within the same structure because of the filling of antibonding states. The relativistic effects result in weaker bonds of Tc-B （Ru-B） than those of Re-B （Os-B）. All the diborides are ultra-incompressible. Re0.5Tc0.5B2 has the largest shear modulus and it is a promising superhard diboride like Os0.5W0.5B2. The elastic properties are in high correlation with the bond strength. The shear moduli are more sensitive than the bulk moduli to the bond strength.