非公企业党建工作汇报材料(高斯贝尔公司)
苏仙区高斯贝尔公司党支部基层组织建设年情况汇报
(201*年11月30日)
高斯贝尔数码科技股份有限公司是一家集贸、工、科于一体的高科技民营企业。现有员工3000多人。主要生产经营数字电视设备、卫星电视接收设备、微波介质陶瓷、安防监控等产品。公司于201*年9月成立了党支部。支部委员会有书记1名,委员2名,现有正式和预备党员22名,入党积极分子20余名,这些同志广泛分布在公司管理层和生产、经营、技术、检验等第一线,是企业的骨干力量,为公司的和谐发展发挥了重要作用。
今年来,公司党支部以基层组织建设年活动为平台,以“同步上,共发展”为目标,加强党建工作,坚持从企业实际出发,充分发挥党组织的核心作用和党员的先锋模范作用,为公司可持续发展提供坚强的组织保障,多次被评为市、区“先进基层党组织”,实现了党建工作和企业发展良性互动。我们的主要做法是:
一、实施“基础工程”,强力推进支部自身建设
公司党支部以制度建设为保障,着力把党支部自身建设作为“车头工程”来抓,推进了党建工作顺利开展。
1、优选支部班子。公司党支部坚持把建设好领导班子摆在各项工作的突出位置,把党性强、作风正、有奉献精神、党员和群众拥护的同志选进支部班子。支部委员会现有3名成员,平均年龄28岁,全部具有本科以上学历,均为公司中层以上干部,确保了支部班子年轻化、知识化、富有活力。同时,实行党支部和经营组织交叉任职,将党的工作和生产经营相结合,无论生产、销售、行政、财务、外贸等部门,都活跃着党员的身
影,确保了党的工作覆盖到企业经营活动的所有领域,不断扩大党组织的影响力和向心力,党支部成了业主和员工“政治生活的知己、经济运作的参谋、法律保护的依托。”
2、激活党员管理。为彻底改变部分党员“角色”错位、“个人主义思想”较强、“党员意识”弱化等问题,党支部把人性化管理贯穿企业管理之中,切实把党员的思想作风建设同党的组织建设相结合、同企业经营活动相结合、同企业文化建设相结合、同提高党员队伍素质相结合,在党团中深入开展爱国主义、社会主义、集体主义教育,引导他们把热爱社会和热爱企业有机结合起来,努力为企业作贡献。同时,把建设学习型企业、培育学习型员工作为党建工作的拓展和延伸,每年投入30万多元对党员干部、骨干员工进行专业培训,并采取电视、网络、多媒体等教学手段,开展党员的理论知识学习和业务能力培训。同时,开展“双培养”活动,把基层一线员工培养成党员,把党员培养成岗位技术能手。近三年来,公司党支部先后组织党员学习27次,举办知识讲座12场次,培训党员和入党积极分子170余人次。“平常时刻看得出,关键时刻站得出,危难时刻豁得出”,党员以模范行动影响、带动了其他员工。
3、创新活动载体。公司党支部充分发挥党的政治优势,将思想政治工作同岗位技能培训相结合,在提高员工生产技能中,增强他们爱岗敬业意识。公司先后成立了文艺宣传队,建立了文化长廊和读书会等,每逢重大节日庆典,公司党支部积极配合区工会、共青团、妇联等群团组织,及时将党员、团员组织起来,举办文艺晚会、卡拉OK、书画创作、体育比赛、劳动技能竞赛等,并积极开展“为企业发展献策”,“岗位能手竞赛”,“人才传帮带”等活动,营造党员争当技术新能手,争做先锋模范的氛围。拓展了思想政治工作载体,使员工在潜移默化中培养了健康向上的情操,大大增强了党组织的凝聚力和感召力。今年“七一”前夕,党支部与工会联
手组织全体党员参观桂东第一军规广场,并重温入党誓词,增强党员的自豪感和责任感。
二、实施“形象工程”,着力彰显共产党员先进性
公司党支部以开展活动为抓手,着力体现党员先进性,引导党员在建树形象中彰显先进性。
1、建言献策聚合力。如何使党建工作“为企业所需要、为业主所理解、为职工所拥护、为党员所欢迎”?这是公司党支部一直探索的问题。党支部紧扣生产经营中心,积极寻找党的方针政策、法规要求与生产经营实际的结合点,建言献策,参与决策,发挥帮助促进的作用,增强了党组织在企业决策中的影响力。每逢重大决策出台前,党支部书记参加经理办公会,在企业经营、人事、工资、保障等方面直接关系到员工利益的问题上都有发言权,积极为公司生产经营出谋划策。为培养一支能打会拼的企业队伍,党支部提出“以德治企、以人为本”的管理新理念,积极倡导德治,为公司董事会采纳,促进了公司的凝聚人心;为支持企业全力以赴抓生产,使产品尽快通过ISO9000质量论证,公司党支部多次组织党员技术骨干进行业务培训,并建议邀请商检部门的专家授课,确保了质量论证顺利通过。今年来,党员提出的10条合理化建议都被董事会采纳,公司完成产品销售额比去年大幅增加,成为郴州市民营企业的中的“标兵和明星”。
2、责任示范当表率。根据党员特长和岗位性质,支部开展“佩戴党徽,亮出党员身份”活动,使党员时刻不忘自己身份,时刻保持责任感和使命感。在此基础上,为表现优秀的党员设置了党员先锋岗,发挥其表率作用。同时,把每个党员所在部门或班组车间,作为党员责任区,要求责任区做到“三有三无”,即困难面前有党员、重点岗位有党员、突击攻关有党员;党员身边无事故、无次品、无违纪。原公司党支部书记肖嗣彪同
志在“为企业建言献策”活动中提出建立“公司科技创新成果展示厅”的建议被董事会采纳,展厅建成后已成为展示公司科技实力,宣传公司最新产品的主要窗口。
3、扶贫助困献爱心。公司党支部建立健全了“委员带党员、党员带群众”等联系帮扶制度,规定每个支委成员联系3-4名党员,每个党员联系1户贫困户,切实开展联系帮扶活动,目前已建立联系帮扶示范点2个,党员共联系贫困户20户,为群众办实事30余件,帮助解决实际问题17个。公司职工贺丽兰家里遭受火灾,损失惨重,支部领导知道情况后,一方面上门看望慰问,另一方面积极组织公司党员捐款捐物,共计2.4万元,使他们一家走出了生活的困境,重新树立了生活的信心。
三、实施“动力工程”,全力促进企业发展壮大
公司党支部本着“有为才有位”的原则,把党建工作同企业发展紧密结合,做到党建工作与企业生产经营目标同向、思想同心、工作同步。
1、发挥部门推动力。党支部充分发挥桥梁纽带作用,通过加强联系,对下及时掌握情况,对上及时反映情况,帮助企业主协调社会关系,为企业发展壮大创造良好的外部环境。在开发新产品过程中,公司党支部积极向上级党组织汇报企业动态和面临的困难,得到了市区领导的大力支持,从项目论证,到立项,到筹集资金,市区领导都亲自挂帅,积极改善外部环境。同时,党支部主动与公安、检察、纪检、劳动、工商等执纪执法部门协调,并切实发挥监督作用,促进企业依法经营,依法纳税,文明经商,既为企业创造了良好的经营环境,又树立了公司的良好形象。
2、聚合外在牵引力。党支部充分利用自身资源优势,多侧面、多角度与所在的农村党组织开展“村企同行、携手共建”活动,建立村企党支部共建联络制度,通过共筑交流平台、共建互助通道,大力开展思想结对、办实事结对和服务结对,取得了良好成效。今年来,公司党支部与观山洞
街道观山洞村党支部进行结对共建“资源共享,优势互补”的原则,把支部的自身建设、经济发展与结对共建工作有机结合起来,努力构建基层党建新格局。村党支部利用熟悉村务、了解情况等优势,帮助企业协调关系,及时解决企业生产中一些涉及群众利益的问题;企业党支部充分发挥自身优势,积极参与村里的公益事业建设,帮助村里转移劳动力,村企双方关系日益融洽,村里和企业相互支持共谋发展,优势互补,实现共赢。
3、开发自身原动力。党支部从企业发展和企业员工利益出发,定期召开企业主与员工见面会和企业员工座谈会,主动架设劳资双方互相沟通的桥梁。每月初支部召开企业主和员工见面会,企业主主动解答员工提出的问题,并有针对性地提出解决办法;每月末召开员工座谈会,员工们简要回顾当月的生产、生活情况,存在和需要解决的问题等,归纳梳理,及时向企业主通报,促进了企业主和员工的内部协调。同时,建立“联系户”,规定每个党员联系1-3名企业员工,1-3名员负责一个车间或班组,促使党支部及时掌握每个班组、车间的工作情况和每个员工的思想动态,从而为加强内部协调提供了客观、可靠的第一手材料,确保党的工作覆盖到企业每个角落,党支部成了关心、维护、代表员工利益的“主心骨”。
扩展阅读:经典高斯总结
一般地,优化所得驻点的性质(极小点还是过渡态)要靠频率来确定;而对过渡态,要确定反应路径(即到底是哪个反应的过渡态)必需要做IRC了,不然靠不住的(往往用QST找到的过渡态并不一定就是连接输入反应物和产物的过渡态)。
在我们用QST2或QST3来优化过渡态时,需输入反应物和产物,实际上反应物和产物的输入顺序是没有关系的。就是说,先输反应物后输产物和先输产物后输反应物得到的是同样的过渡态。这也好理解,QST2里对过渡态的初始猜测实际上是程序自动将输入的反应物和产物的各变量取个平均,所以输入顺序是没有关系的。对QST3和TS,需人为指定过渡态的初始猜测。
上面说的反应物和产物的输入顺序没有关系,有个前提条件,就是反应物和产物的自旋多重度一致。对QST3,因为是人为指定过渡态的初始猜测,所以没有影响;而对QST2,过渡态的自旋多重度默认和后面输入的一致。如果反应物和产物的自旋多重度一致,那随便先输哪个都没关系;而如果反应物和产物的自旋多重度不一致,这时该怎么办???
当反应物基态多重度为1,产物基态多重度为3时,到底将过渡态的多重度定为1还是3?还是两个都试,哪个能量低取哪个?假如取1,在下来做IRC验证时,Reverse还好办,Forward却是按多重度为1做的,这样怎么能和多重度为3的产物连接起来?
要将其关联起来,需有一个为激发态,这样才能在自旋多重度上保持一致。可这种处理对吗?还有更好的思想吗?另外,我觉得IRC难用得很,很难控制,也不知自己钻进了哪条死胡同!!
一直对此很迷茫,希望各位大侠援助啊!
对于这点确实很迷茫,我觉得好像得用cas解决,有看过类似的文献,上面有用cas得出过反应物基态经激发态再回到基态得出产物的反映路径。我没有试过,对于cas我很是头疼,理论和实践都一无所知。作IRC能说明哪些问题?
irc做完,就说明我们完成了这个反应路径地计算,
在SUMMARYOFREACTIONPATHFOLLOWING:,我们可以看到分子过渡态的键长键角与能量随着反应坐标的变化而变化,如果你将反应坐标与能量作图,就可以得到一条过渡态曲线,由于irc计算是结构沿着反应路径的方向,在每个点进行优化的,所以如果你找的过渡态是正确且优化是成功的话,TS确实是连接两个minimum的。
一般IRC的两个终点并不恰好就是反应物和生成物,除非不断加大Maxpoints。然而此时IRC又往往会出错,或者并不能达到指定的Maxpoints,例如,我想正反两个方向各算20各点,IRC计算正常完成,但结果中并没有41个构型。如果这时作能量-反应坐标曲线,如何确定反应物和生成物的反应坐标?如果你的过渡态确实做得漂亮,那么就一定可以成功地由过渡态找到与之对应的反应物、生成物。可这只是理论上的东东,现实中,由于反应的势能面实在太复杂,你千辛万苦所得到的过渡态是很难天衣无缝地连接反应物、生成物的。不断加大Maxpoints简直是Killtime。
一般地,我们可以将IRC最后得到的两个结构的所有数据分别与你做QST时所用的两个结构进行对比,只要结构数据的差别、能量差别比较小,就可以基本认定(因为反应的势能面太复杂)你IRC是成功的了。
说得好!取点越多越好,但我发现往往作不到这一点,所以我也是将IRC最后得到的两个结构的数据分别生产物态和产物态的两个结构进行对比。
而且我发现,作IRC时频率计算给出的讯息很重要,呵呵!过渡态寻找小结
刚刚做了一段时间的过渡态,期间碰到了许多的困难。寻找过渡态不是一件容易的事(对于我和大多数刚涉及量化的人来说),因此我希望通过写这个经验小结能对大家有些帮助。
1.首先遇到的问题是,用哪种方法来寻找过渡态?
GAUSSIAN提供的方法是QSTN和TSN方法。两种方法各有优点和缺点。QSTN方法特别QST3方法要求输入反应物,过渡态的猜测结构,产物这三者的结构。特别麻烦。但很管用,一般不会出现不收敛的情况。对于TSN(对应关键词为OPT=TS)方法,只要求输入过渡态的初始结构,但这个初始结构非常的关键,如果结构不好,则很容易出现不收敛的情况。所以我建议,如果是刚开始做过度态的话,用QSTN方法是好的选择,等有了“感觉”之后,再用TSN方法。
2.怎么解决经常出现的错误?
在找过度态的时候,经常碰到的一些问题就是不收敛(1).,有一个错误的本征值(错误信息为:thereisawrongsigneigenvalueinhessianmatrix.....)(2),和LINK9999错误导致退出。(3)
对于不收敛的情况,可以分为两类,比如提示信息里的CONVERGENCEFAILER提醒收敛到了10(-5),而此时你设定的SCF循环次数也仅仅是64步,那么完全有希望通过加大SCF循环次数来达到收敛的目的。倘若只收敛到10(-3)或10(-2),此时加大循环次数可能就没用了。结果还是CONVERGEFAILER。此时可采用SCF=QC,来达到强制收敛的目的。因为SCF=QC(LINK508)的计算量比默认的L502要大,所以不到万不得以就不用它了。出现第二个错误可以直接用关键词OPT=NOEIGEN来实现。
LINK9999出错是因为已经走完了默认的步数,但还未完成。系统会自动跳出。出现这种情况大多数就是因为优化步数和SCF步数超过了默认值。可用OPT(MAXCYCLE=100)和SCF(MAXCYCLE=300)来改错。
3.怎么样控制过渡态的优化,使得过渡态不至于收敛到其他的分子结构中去?我用GAUSSVIEW可以解决这个问题,当刚开始运行GAUSSIAN时,你用GVIEW去打开输出文件时,你可以看到你的过渡态的初始输入结构,当一个循环过后(从上一个LINK502到下一个LINK502),你再打开输出文件,你就可以清晰地看到优化一步后分子的构型,这样就可以随时监控过度态分子的结构,倘若已经有收敛到其他分子构型的趋势时,你就可以把它给KILL了,而不至于需要等全部工作结束后,打开输出文件才知道已经不是想要的过渡态了。如果收敛到其它的构型上去,可以考虑缩小OPT的步长.iop(1/8=2或3)即可。4.还需要加其它的关键词吗?
建议在OPT中加入CALCFC。这样可以加大找到过渡态的几率。本人深有体会!
先写这么多了,难免有错误和不恰当之处,还希望大家来指点和补充!我找过渡态的经验是:找过渡态的关键点,是如何根据反应体系的立体结构和反应过程化学键的断裂和形成的轨道的空间形状,提出一个合理初始的过渡态的结构,再进行优化。提出的初始过渡态结构越接近真实结构,就越容易找到,否则花一年时间也可能找不到一个过渡态。而找过渡态的QSTN方法、TSN方法的差别不是太大。
这几个帖子都很好,我也找了十几个简单的过渡态,不是直接找到的,而是先用目标过渡态的四元环结构做,然后对四元环结构逐步“修饰”直到那个四元环就是我需要的四元环为止,这个方法来得一般比较简单,但也有无法“修饰”到目标过渡态的情况,其他方法都是一样的。10(-5)表示十的-5次方。
每个循环都会列出SCF计算的结果,和收敛值啊。可以找的到的。iop(1/8)=2
相当于在找过度态的时候,以2A(A为基本单位长度)为单位来寻找过渡态。就象有一百个抽屉里只放了一个苹果,当然是一个一个抽屉打开,找到苹果的把握最大了。此时意味着iop(1/8)=1。
有虚频说明结构处于不稳定状态。过渡态的初始结构猜测应该根据已有的知识和经验啊,哪种构型或构象最可能是过渡态,注意有时候猜测可能与计算结果完全相反。反应坐标
势能面是研究化学反应历程的基础。对一个势能面来说,我们比较关心其关键点。势能面上的关键点是指势能的极值点,包括反应物(R)、产物(P)、中间体(I)与过渡态(TS)或鞍点(saddlepoint)等。连接反应物、过渡态和产物的反应途径,是一条能量最低的路径(minimumenergypath),称为反应坐标(reactioncoordinate)。各关键点是通过反应坐标联结起来的,过渡态是其上面的一级鞍点。
*不同类型的反应坐标
我们常说的反应坐标、最小能量路径等用语,严格地说,是有一定区别的。(1)最速下降路径:也称最小陡降路径。从过渡态出发,体系会沿着最陡的斜坡向反应物深谷和产物深谷移动,移动路径与等能量面正交,方向与势能梯度相反。这条路径的任何地方,其垂直方向有一极小点,就是垂直切割最速下降路径的交点。(2)内禀反应坐标:内禀反应坐标就是连接反应物、过渡态和产物的最速下降路径。
(3)最小能量路径(MEP):通常选定一个反应坐标,使能量相对于其它坐标为最小。这样的反应路径,常常会引起弯曲甚至不连续,特别地与最速下降路径偏离甚远,而后者是真正的最小能量路径。举例来说,RelaxedScan的结果就是一条MEP。
(4)经典轨迹:上述路径在势能面上对应着一组原子无限缓慢地移动(零动量),而轨迹是通过在包括核最初动量的合适初始条件下求解反应物的运动方程得到的。动力学效应,例如离心力,将使分子沿着不同于最速下降的路径移动。随着分子象台球一样从势能山脊上弹起又滑下山坡,轨迹可能十分繁长。irc输入和结果%Chk=irc
#phf/sto-3girc(forward,calcfc)guess=readgeom=allcheck
******************************************%nproc=2%Chk=irc2
#phf/sto-3girc(reverse,calcfc)
guess=readgeom=allcheckscf(maxcyc=200)(results:)
SUMMARYOFREACTIONPATHFOLLOWING:与scan看结果一样是看OPti上的那个orient。
如果是单独作irc(forward)和irc(reverse),取出最后opt结果上面的那个orient来看irc到底是个什么样的过程。第一个是过渡态的orient了。关于虚频的一个简单的理解首先,什么是频率。
中学的时候我们学过简谐振动,对应的回复力是f=-kx,对应的能量曲线,是一个开口向上的二次函数E=kx^2/2.这样的振动,对应的x=0的点是能量极小值点(简单情况下也就是最小值点)。这时的振动频率我们也会求:ω=2πsqrt(k/m)。显然它是一个正的频率,也就是通常意义下的振动频率。
那么,一维情况下,如果能量曲线是一个开口向下的二次曲线呢?首先,从能量上看,这是个不稳定的点,中学的物理书上称为“不稳平衡”。用现在的观点看,就是这一点导数是零(受力为0),且是能量极大值。如果套用上面的公式,“回复力”f=-k"x(实际上已经不是回复,而是让x越来越远了),这里k"是个负数,ω=2πsqrt(k"/m)显然就是一个虚数了,即所谓的虚频。Gaussian里面给出一个负的频率,就是对应这个虚频的。
实际情况下,分子的能量是一个高维的势能面,构型优化的时候,有时得到了极小值点,这样这个点的任意方向上,都可以近似为开口向上的二次函数,这样这里对应的振动频率就都是正的。对于极大值点,在每个方向都是开口向下的二次函数,那么频率就会都是负的当然一般优化很少会遇到这样的情况。对于频率有正有负的情况,说明找到的点在某些方向上是极大值,有些方向上是极小值。如果要得到稳定的能量最低构型,显然需要通过微调分子的构型,消去所有的虚频。如何微调?要看虚频的振动方向。想象着虚频对应的就是开口向下的二次函数,显然,把分子坐标按照振动的方向移动一点点,分子应该就可以顺着势能面找到新的稳定点,但是也不能太小。而所谓的过渡态,则是连接反应物和产物之间的最低能量路径上的能量极大值。好比山谷中的A,B两点,它们之间的一个小土丘,就是过渡态,从A到B的反应,需要越过的是这个小土丘,而不是两边的高山。这样,过渡态就是在一个方向上是极大值,而在其它方向上都是极小值的点。因此,过渡态只有一个虚频。
频率分析只能在势能面的稳定点进行,这样,频率分析就必须在已经优化好的结构上进行。频率分析的另外一个用处是判断稳定点的本质。稳定点表述的是在势能面上力为零的点,它即可能是极小值,也可能是鞍点。极小值在势能面的各个方向都是极小的。而鞍点则是在某些方向上是极小的,但在某一个方向上是极大的,因为鞍点是连接两个极小值的点。有关鞍点的信息:1.负的频率2.频率相应简正振动的模式
当一个结构产生负的振动频率时,可以表明在该振动方向可能存在着能量更低的结构。判断所得鞍点是不是需要的鞍点的方法,就是察看它的简正振动模式,分析是不是可以导向所需要的产物或反应物。进一步的,更好的办法是通过IRC计算来判断反应物,产物与得到的鞍点是否有关系。
向大家推荐一本书,翻译成中文的(原文我没有见过),可能是《化学反应中的电子》。
在寻找过渡态计算中opt=ts,需要计算频率,如果有一个虚频就说明优化所得的构型是个过渡态。有时在做部分限制优化opt=minimal时,也能得到一个虚频(一般发生在frozen部分的化学键上)。我想问一下,虚频的数值大小本身具有什么物理意义。后一种情况是否也可以认定为过渡态。我的理解:虚频就是在这个振动方向上力常数是负的
过渡态之所以有一个虚频,是因为它是一个鞍点,鞍点上有一个负的梯度的方向分子沿这个方向振动时将转化为反应物或产物。就象从山上掉下来,受到的力可以认为是负的。但并不是有一个虚频就是过渡态,还要在这个虚频的振动方向上分别指向反应物和产物
Gaussian程序频率计算中的几个问题
在Gaussian计算中,为了确定优化得到的几何结构是势能面上的局域极小点还是鞍点,或者要得到相关的热力学性质,经常需要对优化后的几何结构进行振动分析。这里我们将讨论几个频率计算中常见的一些问题。希望能对初学Gaussian的人有所帮助。首先,原则上说,振动频率分析只对稳定结构有意义。这里所说的稳定结构包括是势能面上的局域极小点和鞍点。如下图1所示是一维自由度上的势能面,A和B处在势能面的局域极小点,而处在势能面的鞍点上。他们在都处在平衡位(原子核受力为零),不同的是,A和B来说离开平衡位会受到指向平衡位处的力,而C离开平衡位会受到远离平衡位的力。因此A和B处在稳定平衡点,C处在不稳定平衡点。实际上,一个分子可以有很多的自由度,如果在所有自由度上分子都处在稳定平衡,就是稳定的分子。频率分析得结果是所有频率都是正的,表明这是一个局域的极小点。如果分子只在一个自由度上处于不稳定平衡位,其他自由度上都处在稳定平衡位,说明该结构是一阶鞍点。分子在稳定自由度方向上的振动才是真实的振动,在不稳定自由度方向上的实际上是不会有振动的。不过我们可以对不稳定方向上的运动也按振动来做数学处理,会的到负的振动频率,我们称它为虚频。虚频的出现表明该结构为鞍点。图1势能面上的局域极小点和鞍点
第二,Gaussian计算中,频率的计算一定要在和分子结构优化相同的方法,基组下进行,否则计算的结果是没有意义的。我们知道,任何理论水平下的计算,都是在一定的近似下进行的,不同的理论水平的近似程度是不同的。在一种理论水平A下优化的稳定结构Geom_A会和另一种理论水平B下优化的稳定结构Geom_B有差别,也就是说Geom_A不会是理论水平B下的稳定结构。根据前面我们所讨论的,在理论水平B下对一个不稳定的结构进行频率分析是没有意义的。图2示意说明了不同理论水平下稳定点结构的不同。图2不同理论水平下优化的稳定结构是不同的
第三,频率计算中可以考虑同位素效应(Freq=ReadIsotopes)。在波恩-奥本海默近似下,对于同一种元素采用不同的同位素对几何优化和电子结构计算没有影响,频率计算所需的力常数矩阵(Hessian矩阵)也不会变化,变化的只是约化质量。容易理解,重的同位素会导致低的振动频率。实际上,原子序数大的元素的同位素效应非常不明显,一般只需考虑H原子的同位素效应。第四,各种方法计算的频率和实验结果之间存在系统误差,需要乘以一个约化因子来进行校正(Scale=f)。一般来说,理论计算的频率值会比实验结果大。下面是一些理论水平下的约化因子。注意频率和零点能的约化因子是可以不同的。更多水平下的约化因子需要查文献获得。
方法:约化因子约化因子(频率)(ZPE)HF/3-21G0.90850.9409HF/6-31G(d)0.89290.9135MP2(Full)/6-31G(d)0.94270.9646MP2(FC)/6-31G(d)0.94340.9676BLYP/6-31G(d)0.99401.0119B3LYP/6-31G(d)0.96130.9804SVWN/6-31G(d)0.98331.00
第五,Gaussian的频率计算有时会遇到下面的警告:
Warning--explicitconsiderationof3degreesoffreedomasvibrationsmaycausesignificanterror
这时一般有两种可能:一种可能是,优化的几何结构不够精确,还没有达到稳定点。对于这种情况,需要考虑用OPT=tight或OPT=VeryTight,结合Int=fine或者Int=VeryFine进行更加精确的优化。第二种可能是,在计算的振动模式中,有部分的低频模式对应于内转动模式,在热力学分析中应该按自由转子或受阻转子模型处理,如果按谐振子模型处理会造成较大的误差。采用受阻内转子模型,可以用Freq=Hindered-Rotor计算或者自己手动计算这些振动模式对热力学性质的贡献。
第六,对于计算出的力常数较小的振动模式,其势能面形状可能不满足谐振子模型的要求,力常数函数需要考虑非谐性的成分,这时可以用Freq=Anharmonic来进行非谐性处理。
第七,在频率输出前,会有这样的输出:
Lowfrequencies----0.00080.00030.001340.627559.380866.4408Lowfrequencies---1799.18923809.46043943.3536
上面行的六个模式实际对应于平动(前三个)和转动(后三个),理论上,这六个数都应该是零。如果是对过渡态进行的频率计算,第一行会先输出虚频,然后才给出平动转动的数值。一般来说,平动的三个数都是接近于零的。如果转动的三个数不接近于零(一般解析频率10个波数以内,数值频率50个波数以内),说明需要进行OPT=Tight或OPT=VeryTight计算。第二行的几个实频率应该和随后输出中的相应频率进行对比,如果基本一样,说明频率计算结果很好;如果差别较大,说明这些频率受到的平动和转动的污染较大。
第八,频率计算后程序会进行一步优化计算,正常会得到四个判据都是YES的结果,这说明优化的结构很好。但是有时也会遇到四个判据不全是YES的情况,这时需要仔细分析,一般有三种情况:1、Force和RMSForce都是收敛的,Displacement和RMSDisplacement虽然是NO,但是都比较接近收敛判据;2、如果Force和RMSForce都是收敛的,Displacement和RMSDisplacement不收敛且数值远大于收敛判据;3、Force和RMSForce不收敛。对于第一种情况,我们可以不管它,仍然可以认为计算结果是可靠的。对于第二和第三种情况,一般说明优化过程中估算的Hessian是不准确的,优化的结构可能还没有达到稳定点,需要重新优化。如果反复优化仍然无法解决这个问题,建议在优化过程中使用OPT=CalcAll。
第九,优化一个局域极小点,收敛后从力常数本征值看应该是局域极小点(没有负的本征值),但是频率计算中会出现虚频和负的本征值。这种情况下虚频一般是由转动模式造成的,说明分子中两个基团之间的相互位不是很合适,需要绕转动的键相对转动到一个合适的位,重新优化。绕某个键转动两个基团,有时可以很方便地用修改二面角的方法实现:OPT=Modredundant结合分子描述后输入:*mn*[+=]value。其中m,n为连接两个基团的键的顶端原子。有时,用OPT=CalcAll也可以解决这个问题。第十,Gaussian程序中,频率计算是和温度、压力无关的。因为Gaussian所考虑的是原子核在一定核构型和电子运动状态构成的势场内振动,一般来说,温度和压力的变化不会影响分子的结构和电子运动状态,所以振动的势能函数是与温度、压力无关的。Freq=Isotopes关键词要求输入的温度和压力影响的只是平动,转动和振动的平衡分布和它们对热力学性质的贡献。
分子激发态计算的理论方法一1.激发态计算方法的分类
按照化学模型,目前计算激发态的方法可以分为:____DFT(TDDFT,MC-DFT,DFT+CI...)Chemical____|
Model|____Wavefunction____|~~~~Abinitio(MR-CI,CIS...)|____Semiempirical(Zindo/1,INDO/S,CNDO/S...)
这与基态计算方法的分类是一样的,从中看不出激发态计算的特点。另一种分类是按照获得激发能(或激发态能量)的方式,可以分为两类:i)电子组态方法。选择不同的电子组态的线性组合进行迭代,得到激发态的绝对能量(有些程序经过简单的后期处理也能给出相对能量)。分为多参考(MR)方法和单参考(SR)方法两种。常见的多参考方法有:MCSCF(有CASSCF,RASSCF等多种形式),MR-CISD,CASPT2等。标准的
MR的计算过程一般是三步:HF〉MCSCF(常用CASSCF)〉MRCISD或CASPT2...R是MR的特例,不需要中间的MCSCF步骤,计算中只需要对MR程序模块指定一个参考组态即可,计算量大大减少,但是精度也会下降。对于一些无法用SR很好描述的体系(如镧系化合物),计算是失败的,这时就必须用MR方法。ii)线性响应理论(LRT)。基态用普通的理论方法求得,然后通过求解特殊的线性响应方程,获得激发态相对于基态的垂直激发能(相对能量)。有的方法求解方程需要做迭代,有的则不需要。这一方法要求作为参考态的基态必须被很好地描述。常见的有TDDFT,CIS,CCSD-LRT等,共同点是基态都使用SR模型。基态用MR模型的有MCLR(MCSCF-LRT),MR-AQCC-LRT等方法,用于计算那些无法用SR模型很好描述基态的体系。2.一些方法的特点及软件
i)Zindo/1,INDO/S,CNDO/S都属于半经验方法,可用于研究原子数在100左右的大分子体系。软件:Zindo/1Gaussian,Zindo。Zindo结合ZINDO-MN模块可以用Zindo/1方法计算溶剂对激发态的影响。INDO/S,CNDO/S商业MOPAC的MOS-F模块。免费的MOPAC6.0/7.0无此功能
ii)CIS基态能量用Hartree-Fock方法得到,激发态精度也是HF级别的,得到半定量的结果。适于研究原子数在50左右的中等分子体系。此外还有精度更高一些的CIS(D),计算量当然更大。Gaussian03的CIS还能考虑溶剂环境对激发态影响。程序化的这些方法一般只能计算单重态和三重态。扩展CIS(XCIS)方法还能计算双重态和四重态。
软件:CISGaussian,Q-Chem,CaChe,商业的MOPAC,TURBOMOLE,最新版MWChemCIS(D)Q-Chem,TURBOMOLEXCISQ-Chem
iii)TD-DFT和TD-HF目前比较流行的激发态计算方法,适用于十几个原子构成的中小分子,
小分子计算也能得到令人满意的结果。TD-DFT目前的缺陷是无法研究多电子激发和Rydberg态。TD-HF早期也叫RPA,是TD-DFT的特殊形式,精度较差。Gaussian03中TDHF/TDDFT的激发态计算还能考虑溶剂环境影响。程序化的TD-DFT一般只提供对单重态和三重态计算的正式支持。最新版
NMChem的TD-DFT除了计算闭壳层分子的单重态和三重态,还能计算开壳层分子的双重态和四重态。软件:
TD-HFGaussian,ADF,NWChem,Q-Chem,DeMon,GAMESS-US,ACESII,TURBOMOLETD-DFT
Gaussian,
ADF,
NWChem,
Q-Chem。DeMon,
TURBOMOLE,CADPAC,BDF。
一些计算晶体的程序也可以做分子的TD-DFT计算,但是由于使用平面波基组,结果不太好。
iv)CC-LRTCC-LRT包括CCS-LRT(等价于CIS),CC2-LRT(对CCSD-LRT的近似),CCSD-LRT,CC3-LRT(对CCSDT-LRT的近似)等不同激发级别的形式。最常用的是CCSD-LRT,有很多种形式,如EOM-CCSD,VOOD-CCSD等。适用于中小分子和小分子的计算,结果与TD-DFT相当,可能还略好一些。软件:
CC-LRTDALTON。CCSD-LRT的其他形式可在ACESII,Q-Chem,PSI,MOLPRO中找到。TURBOMOLE可做CC2-LRT计算。v)SAC-CI。
精度略高于TD-DFT的方法,适于研究十几个原子构成的有机分子,还可以包含多电子激发的影响。软件:SAC-CIGaussian03vi)DFT+CI混合方法。
结合了DFT和CI的优点,适用于中小分子激发态计算。软件:DFT+CITURBOMOLE
vii)MCSCF(有CASSCF,RASSCF等)。
做(电子组态方法)多参考计算的基础。适用于中小分子激发态计算。常用电子组态形式的CASSCF。线性响应理论形式(MCLR)用的较少。软件:MCSCF
GAMESS-US,
GAMESS-UK,
Gaussian,
DALTON,
MOLPRO,MOLCAS,COLUMBUS,NWChemMCLRDALTON,GAMESS-UK
viii)MR-CI及其变体MR-CISD是计算(几个原子构成的)小分子激发态的标准方法。由于这一方法存在着大小一致性问题,因而又发展了很多变体。比较简单的是对MR-CISD做事后修正(+Q),主要有Davidson修正和Pople修正两种。复杂一些的有MR-ACPF、MR-AQCC、MR-CEPA等,其中表现最好的是MR-AQCC,以及后来发展的虚轨道修正MR-AQCC(MR-AQCC-V)。有的程序还允许对MR-CI做激发级别的修正,得到MR-CIS,MR-CISD(T)等。软件:MR-CISDGAMESS-US,GAMESS-UK,DALTON,MOLPRO,MOLCAS,DIRAC,MOLFDIR,COLUMBUS
变体:MOLPRO,MOLCAS,COLUMBUS。COLUMBUS还能做MR-AQCC-LRT计算,此方法综合了MR-AQCC与线性响应理论的特点。MOLPRO对MR-CI及其变体使用了内部收缩,速度更快。
ix)MR-MPn多参考微扰理论。目前主要有二阶微扰和三阶微扰。在程序作了优化的前提下,它们比相同级别的MR-CI计算速度快(例如在MOLPRO中,CASPT2比MR-CISD快),但是精度略低。目前应用的难题是入侵态问题。有些程序通过与MR-CI混合可以部分地解决这一问题。
软件:CASPT2MOLPRO,MOLCAS,COLUMBUS,GAMESS-UK。GAMESS-US的MC-QDPT方法和Gaussian的CASSCF+MP2也属于MR-MP2。
CASPT3MOLPRO,GAMESS-UK
x)MR-CC和FCIMR-CI的理论极限是完全组态相互作用(FCI),常用的MR-CISD是把FCI在CI空间的展开截断到二级的形式。
鉴于耦合簇CC能获得比CI更好的结果,而且不会有大小一致性的问题,如果把FCI展开到CC空间并在一定的级别截断,就得到MR-CC。
由于FCI计算量相当大,所以目前也只能计算较轻原子构成的双原子分子。从理论上讲,MR-CCSD是小分子激发态计算最实用、最精确的方法。但目前还有入侵态的问题没有解决。公开发布的程序中目前还没有此方法。xi)MC-DFT,CAS-DFT:多参考的DFT,尚在研究中。MOLCAS声称将在6.0版包含此方法。
xii)MR-Gn多参考的Gaussian-n理论,可想而知是相当昂贵的方法。目前只有GAMESS-US小组的一篇文章研究过这一理论。将来能否用于激发态计算并获得高精度结果还未知。3.方法的选择构成分子的原子数:
2.....10.......20........50..........100..........
|_____________________||______________|半经验CIS|_____________________|MCSCF(注1)|______________|TD-DFT(注2),CCSD-LRT|_____|
MR-CI及其变体,MR-MPn(注2)|FCI
可见,越往下精度越高,适用分子越小。
注1:CASSCF的计算量和活性空间大小2A和活性电子数E有关。C(2A,E)=(2A)!/[(2A-E)!*E!]。C(2A,E)越大则计算量越大。RASSCF由于对组态进行了筛选,能适当减少计算量。对于小分子,MCSCF一般只作为高级MR计算(MR-CI,MR-MPn)的中间步骤。
注2:对于含有重元素的小分子激发态计算,还要考虑自旋轨道耦合(SO)造成的电子态分裂。例如AtBr的3Pi会分裂为0-,0+,1,2四个态。能够在激发态计算中包含SO的程序有MR-CISD及其变体+SO:COLUMBUS,MOLPRO,GAMESS-USMR-MPn+SO:MOLPRO,GAMESS-US
MCSCF+SO:COLUMBUS,MOLPRO,GAMESS-US,MOLCAS,TD-DFT+SO:ADF201*,尚在测试中,目前只能对原子获得较好的激发能完全相对论MR-CISD:DIRAC,MOLFDIR完全相对论TD-DFT:BDF,尚在测试中
如果是用Gaussian或者GAMESS-US,我的理解是:1.先用AM1进行几何优化。2.用优化的结构进行Hartree-Fock计算。3.用第二步的分子轨道进行CIS计算(Gaussian使用Guess=Read读入第二步的checkpoint文件;GAMESS使用guess=moread,并把第二步的MO复制到$VEC中)。Gaussian的CIS自动给出激发能和谐振强度,GAMESS只给出激发能,计算谐振强度需要另外的计算步骤。4.Gaussian中的Zindo和CIS使用方法相同,速度更快,但适用的原子有限,如果计算量不是太大的话,还可以用TD,结果更好。另外Gaussian的CIS只能算单态、三态,其它的像双态、四重态、五重态可能也会给出结果,但是分析很麻烦。用GAMESS进行CIS计算可以定义任意多重度,比较直观。
友情提示:本文中关于《非公企业党建工作汇报材料(高斯贝尔公司)》给出的范例仅供您参考拓展思维使用,非公企业党建工作汇报材料(高斯贝尔公司):该篇文章建议您自主创作。
来源:网络整理 免责声明:本文仅限学习分享,如产生版权问题,请联系我们及时删除。