MOMAP对azulene的S1→S0的辐射速率和内转换速率的计算结构优化
我们使用bdf进行azulene的s1激发态结构优化和频率计算。选中simulator → bdf → bdf，界面中设置参数。在basic settings界面中的calculation type选择tddft-opt+freq，方法采用默认的gb3lyp泛函，基组在basis中的all electron类型中，选择6-31g(d,p)。scf settings和tddft settings面板中将use mpec+cosx acceleraton的默认勾选去掉。basic settings、scf settings、tddft settings界面中的其它参数以及opt settings、freq settings等面板的参数使用推荐的默认值，不需要做修改。之后点击 generate files 即可生成对应计算的输入文件。
选中生成的bdf.inp文件，右击选择open with，打开该文件，如下所示：
$compass
title
c10h8
geometry
c 0.79273796 0.49102542 -0.00003307
c -0.70229649 0.61186591 0.00000000
c 1.30022932 1.80163337 -0.00006272
c -0.99262499 1.98726800 0.00007812
c 1.54415132 -0.67887418 0.00000000
c -1.63318173 -0.42094563 -0.00002837
c 0.21877157 2.69859813 0.00000000
c 1.10656346 -2.00562676 0.00005788
c -1.41619168 -1.80093044 -0.00004814
c -0.20258112 -2.49333483 0.00000000
h 2.35092512 2.06249889 -0.00009828
h -1.98777600 2.41348149 0.00017650
h 2.62424717 -0.53731745 0.00001117
h -2.67585843 -0.10561277 -0.00001521
h 0.30641472 3.77916915 0.00002386
h 1.88966566 -2.75951313 0.00017581
h -2.31053950 -2.41870505 -0.00009019
h -0.29054446 -3.57807510 0.00000000
end geometry
basis
6-31g(d,p)
skeleton
group
c(1)
$end
$bdfopt
solver
1
maxcycle
108
iopt
3
hess
final
$end
$xuanyuan
direct
$end
$scf
rks
charge
0
spinmulti
1
dft
gb3lyp
d3
molden
$end
$tddft
imethod
1
isf
0
idiag
1
iroot
6
istore
1
$end
$resp
geom
method
2
nfiles
1
iroot
1
$end
选中bdf.inp文件，右击选择run提交作业，任务结束后bdf.out，bdf.out.tmp，bdf.scf.molden三个结果文件会出现在project中。 选择bdf.out，右击show view，在optimization对话框中，显示结构已经达到收敛标准。
在frequency对话框中，检查频率，若不存在虚频证明结构已经优化到极小点。
选择bdf.out.tmp，右击open with containing folder，打开bdf.out.tmp，在文件开始向下查找第一个tddft计算模块‘module tddft’。根据tddft模块的ground to excited state transition electric dipole moments (au) 中的state 1的跃迁电偶极矩，得到momap需要的参数edma，计算方法为： 需要将单位a.u.转换为debye，因此edma= 0.3646*2.5417=0.9267 debye。
在文件末尾向上查找第一个tddft计算模块‘module tddft’。根据tddft模块的ground to excited state transition electric dipole moments (au) 中的state 1的跃迁电偶极矩，得到momap需要的参数edme，计算方法为： 需要将单位a.u.转换为debye，因此edma= 0.3646*2.5417=0.6507 debye。
在该tddft模块的statistics for [dvdson_rpa_block]:中读取no. 1态的能量为 -385.8030478000 a.u.，即为s1态的单点能。 将s1态的单点能与s0态的单点能相减，即得momap需要的两态的能量差参数ead=0.07502 a.u.。
基于基态结构，做s0和s1之间的非绝热耦合计算。点击azulene_s0.hzw，右击点击copy，设置new file name为nacme，在project中出现nacme.hzw。双击nacme.hzw，使用bdf进行azulene的非绝热耦合计算。选中simulator → bdf → bdf，界面中设置参数。在basic settings界面中的calculation type选择tddft-nac，方法采用默认的gb3lyp泛函，基组在basis中的all electron类型中，选择6-31g(d,p)。scf settings和tddft settings面板中将use mpec+cosx acceleraton的默认勾选去掉。在tddft settings面板中的non-adiabatic coupling内容框中，在默认的coupling between ground and excited-state下，点击‘+’号，basic settings、scf settings、tddft settings界面中的其它参数以及opt settings、freq settings等面板的参数使用推荐的默认值，不需要做修改。之后点击 generate files 即可生成对应计算的输入文件。
选中生成的bdf.inp文件，右击选择open with，打开该文件，如下所示：
$compass
title
c10h8
geometry
c 0.79273796 0.49102542 -0.00003306
c -0.70229648 0.61186591 0.00000000
c 1.30022931 1.80163336 -0.00006271
c -0.99262499 1.98726799 0.00007812
c 1.54415131 -0.67887417 0.00000000
c -1.63318173 -0.42094562 -0.00002837
c 0.21877157 2.69859812 0.00000000
c 1.10656346 -2.00562675 0.00005788
c -1.41619168 -1.80093044 -0.00004813
c -0.20258112 -2.49333482 0.00000000
h 2.35092512 2.06249888 -0.00009827
h -1.98777599 2.41348149 0.00017649
h 2.62424717 -0.53731744 0.00001117
h -2.67585843 -0.10561277 -0.00001520
h 0.30641472 3.77916915 0.00002385
h 1.88966565 -2.75951312 0.00017580
h -2.31053950 -2.41870504 -0.00009018
h -0.29054446 -3.57807510 0.00000000
end geometry
basis
6-31g(d,p)
skeleton
group
c(1)
$end
$xuanyuan
direct
$end
$scf
rks
charge
0
spinmulti
1
dft
gb3lyp
d3
mpec+cosx
molden
$end
$tddft
imethod
1
isf
0
idiag
1
iroot
6
istore
1
$end
$resp
quad
fnac
norder
1
method
2
nfiles
1
single
states
1
1 1 1
$end
选中bdf.inp文件，右击选择run提交作业，任务结束后bdf.out，bdf.scf.molden三个结果文件会出现在project中。 至此，momap对azulene的s1→s0的辐射速率和内转换速率的计算需要的bdf量化软件的结构优化频率结果文件、非绝热耦合结果文件和参数部分都已完成。 下一期会接着介绍使用momap进行azulene的s1→s0的辐射速率和内转换速率的计算。
巴克莱银行试水英国无线支付领域
亚洲如何主导工业机器人的应用
物联网给海上油气检测带来了什么改变
初识 Python 正则表达式
2020年最佳供应链企业榜单出炉，阿里巴巴成入选的唯一中企
MOMAP对azulene的S1→S0的辐射速率和内转换速率的计算结构优化
iPhone8究极曝光，竟然这么丑？
润欣科技超低功耗MCU+BLE方案分享
中芯国际“NAND闪存器件及其形成方法”专利获授权
移动支付市场百花齐放 谁来为“安全”保驾护航？
中国电信移动用户规模已达到了3.36亿户4G用户渗透率达到了83.8%
两款CD4069构成的简易振荡电路详解
融合云计算技术的楼宇对讲系统实现智能互联
单片机控制的多参数实时采集与处理系统设计
什么是真正的全能旗舰？Find X3 Pro告诉你答案
自动浇水机的制作方法
我国首座兆瓦级太阳能塔式热发电站发电成功
cx20106a引脚图及功能
真无线蓝牙耳机有哪些推荐
一号双终端开通城市有哪些
选中生成的bdf.inp文件，右击选择open with，打开该文件，如下所示：
$compass
title
c10h8
geometry
c 0.79273796 0.49102542 -0.00003307
c -0.70229649 0.61186591 0.00000000
c 1.30022932 1.80163337 -0.00006272
c -0.99262499 1.98726800 0.00007812
c 1.54415132 -0.67887418 0.00000000
c -1.63318173 -0.42094563 -0.00002837
c 0.21877157 2.69859813 0.00000000
c 1.10656346 -2.00562676 0.00005788
c -1.41619168 -1.80093044 -0.00004814
c -0.20258112 -2.49333483 0.00000000
h 2.35092512 2.06249889 -0.00009828
h -1.98777600 2.41348149 0.00017650
h 2.62424717 -0.53731745 0.00001117
h -2.67585843 -0.10561277 -0.00001521
h 0.30641472 3.77916915 0.00002386
h 1.88966566 -2.75951313 0.00017581
h -2.31053950 -2.41870505 -0.00009019
h -0.29054446 -3.57807510 0.00000000
end geometry
basis
6-31g(d,p)
skeleton
group
c(1)
$end
$bdfopt
solver
1
maxcycle
108
iopt
3
hess
final
$end
$xuanyuan
direct
$end
$scf
rks
charge
0
spinmulti
1
dft
gb3lyp
d3
molden
$end
$tddft
imethod
1
isf
0
idiag
1
iroot
6
istore
1
$end
$resp
geom
method
2
nfiles
1
iroot
1
$end
选中bdf.inp文件，右击选择run提交作业，任务结束后bdf.out，bdf.out.tmp，bdf.scf.molden三个结果文件会出现在project中。 选择bdf.out，右击show view，在optimization对话框中，显示结构已经达到收敛标准。
在frequency对话框中，检查频率，若不存在虚频证明结构已经优化到极小点。
选择bdf.out.tmp，右击open with containing folder，打开bdf.out.tmp，在文件开始向下查找第一个tddft计算模块‘module tddft’。根据tddft模块的ground to excited state transition electric dipole moments (au) 中的state 1的跃迁电偶极矩，得到momap需要的参数edma，计算方法为： 需要将单位a.u.转换为debye，因此edma= 0.3646*2.5417=0.9267 debye。
在文件末尾向上查找第一个tddft计算模块‘module tddft’。根据tddft模块的ground to excited state transition electric dipole moments (au) 中的state 1的跃迁电偶极矩，得到momap需要的参数edme，计算方法为： 需要将单位a.u.转换为debye，因此edma= 0.3646*2.5417=0.6507 debye。
在该tddft模块的statistics for [dvdson_rpa_block]:中读取no. 1态的能量为 -385.8030478000 a.u.，即为s1态的单点能。 将s1态的单点能与s0态的单点能相减，即得momap需要的两态的能量差参数ead=0.07502 a.u.。
基于基态结构，做s0和s1之间的非绝热耦合计算。点击azulene_s0.hzw，右击点击copy，设置new file name为nacme，在project中出现nacme.hzw。双击nacme.hzw，使用bdf进行azulene的非绝热耦合计算。选中simulator → bdf → bdf，界面中设置参数。在basic settings界面中的calculation type选择tddft-nac，方法采用默认的gb3lyp泛函，基组在basis中的all electron类型中，选择6-31g(d,p)。scf settings和tddft settings面板中将use mpec+cosx acceleraton的默认勾选去掉。在tddft settings面板中的non-adiabatic coupling内容框中，在默认的coupling between ground and excited-state下，点击‘+’号，basic settings、scf settings、tddft settings界面中的其它参数以及opt settings、freq settings等面板的参数使用推荐的默认值，不需要做修改。之后点击 generate files 即可生成对应计算的输入文件。
选中生成的bdf.inp文件，右击选择open with，打开该文件，如下所示：
$compass
title
c10h8
geometry
c 0.79273796 0.49102542 -0.00003306
c -0.70229648 0.61186591 0.00000000
c 1.30022931 1.80163336 -0.00006271
c -0.99262499 1.98726799 0.00007812
c 1.54415131 -0.67887417 0.00000000
c -1.63318173 -0.42094562 -0.00002837
c 0.21877157 2.69859812 0.00000000
c 1.10656346 -2.00562675 0.00005788
c -1.41619168 -1.80093044 -0.00004813
c -0.20258112 -2.49333482 0.00000000
h 2.35092512 2.06249888 -0.00009827
h -1.98777599 2.41348149 0.00017649
h 2.62424717 -0.53731744 0.00001117
h -2.67585843 -0.10561277 -0.00001520
h 0.30641472 3.77916915 0.00002385
h 1.88966565 -2.75951312 0.00017580
h -2.31053950 -2.41870504 -0.00009018
h -0.29054446 -3.57807510 0.00000000
end geometry
basis
6-31g(d,p)
skeleton
group
c(1)
$end
$xuanyuan
direct
$end
$scf
rks
charge
0
spinmulti
1
dft
gb3lyp
d3
mpec+cosx
molden
$end
$tddft
imethod
1
isf
0
idiag
1
iroot
6
istore
1
$end
$resp
quad
fnac
norder
1
method
2
nfiles
1
single
states
1
1 1 1
$end
选中bdf.inp文件，右击选择run提交作业，任务结束后bdf.out，bdf.scf.molden三个结果文件会出现在project中。 至此，momap对azulene的s1→s0的辐射速率和内转换速率的计算需要的bdf量化软件的结构优化频率结果文件、非绝热耦合结果文件和参数部分都已完成。 下一期会接着介绍使用momap进行azulene的s1→s0的辐射速率和内转换速率的计算。
巴克莱银行试水英国无线支付领域
亚洲如何主导工业机器人的应用
物联网给海上油气检测带来了什么改变
初识 Python 正则表达式
2020年最佳供应链企业榜单出炉，阿里巴巴成入选的唯一中企
MOMAP对azulene的S1→S0的辐射速率和内转换速率的计算结构优化
iPhone8究极曝光，竟然这么丑？
润欣科技超低功耗MCU+BLE方案分享
中芯国际“NAND闪存器件及其形成方法”专利获授权
移动支付市场百花齐放 谁来为“安全”保驾护航？
中国电信移动用户规模已达到了3.36亿户4G用户渗透率达到了83.8%
两款CD4069构成的简易振荡电路详解
融合云计算技术的楼宇对讲系统实现智能互联
单片机控制的多参数实时采集与处理系统设计
什么是真正的全能旗舰？Find X3 Pro告诉你答案
自动浇水机的制作方法
我国首座兆瓦级太阳能塔式热发电站发电成功
cx20106a引脚图及功能
真无线蓝牙耳机有哪些推荐
一号双终端开通城市有哪些