电子材料与器件计算方法基础——学习使用nwchem——homework3_2

已于 2023-05-20 09:52:11 修改
阅读量348 收藏
点赞数 2
分类专栏： 电材作业 文章标签： 学习 bash linux
于 2023-05-20 09:45:58 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_46538207/article/details/130768041
版权

在这里插入图片描述
在这里插入图片描述

（b）

在目录下创建nwchem与tmp文件夹
修改run_job.sh如下：

#!/bin/bash

nohup /opt/nwchem/bin/LINUX64/nwchem Benzene.nw >& 1.out &

修改Benzene.nw文件中路径

permanent_dir /home/dft_user26/home_work3/problem_2/problem_2/nwchem
scratch_dir /home/dft_user26/home_work3/problem_2/problem_2/tmp

运行程序

bash run_job.sh

运行完成后，打开1.out文件查看结果

搜索关键词Root，可查看相应的激发态
得到如下结果：
单重态

  ----------------------------------------------------------------------------
  Root   1 singlet a              0.195708784 a.u.                5.3255 eV
  ----------------------------------------------------------------------------
     Transition Moments    X  0.00000   Y -0.00000   Z -0.00000
     Transition Moments   XX  0.00000  XY  0.00000  XZ  0.00000
     Transition Moments   YY -0.00000  YZ -0.00000  ZZ  0.00000
     Dipole Oscillator Strength                         0.00000

     Occ.   20  a   ---  Virt.   22  a     -0.70653
     Occ.   21  a   ---  Virt.   23  a      0.70653
  ----------------------------------------------------------------------------
  Root   2 singlet a              0.201045907 a.u.                5.4707 eV
  ----------------------------------------------------------------------------
     Transition Moments    X  0.00000   Y  0.00000   Z  0.00000
     Transition Moments   XX  0.00000  XY  0.00000  XZ -0.00000
     Transition Moments   YY -0.00000  YZ  0.00000  ZZ -0.00000
     Dipole Oscillator Strength                         0.00000

     Occ.   16  a   ---  Virt.   36  a      0.06857
     Occ.   17  a   ---  Virt.   32  a      0.05767
     Occ.   18  a   ---  Virt.   31  a     -0.05767
     Occ.   19  a   ---  Virt.   25  a     -0.13947
     Occ.   20  a   ---  Virt.   23  a     -0.68758
     Occ.   21  a   ---  Virt.   22  a     -0.68758
  ----------------------------------------------------------------------------
  Root   3 singlet a              0.271514533 a.u.                7.3883 eV
  ----------------------------------------------------------------------------
     Transition Moments    X  0.00000   Y  2.27580   Z -0.00000
     Transition Moments   XX -0.00000  XY  0.00000  XZ  0.00000
     Transition Moments   YY -0.00000  YZ  0.00000  ZZ -0.00000
     Dipole Oscillator Strength                         0.93750

     Occ.   12  a   ---  Virt.   29  a      0.05483
     Occ.   14  a   ---  Virt.   34  a     -0.06791
     Occ.   15  a   ---  Virt.   33  a      0.06791
     Occ.   16  a   ---  Virt.   28  a     -0.09259
     Occ.   16  a   ---  Virt.   34  a      0.07716
     Occ.   17  a   ---  Virt.   31  a     -0.10474
     Occ.   17  a   ---  Virt.   35  a     -0.07663
     Occ.   18  a   ---  Virt.   32  a      0.10474
     Occ.   20  a   ---  Virt.   22  a      0.67952
     Occ.   21  a   ---  Virt.   23  a      0.67952
  ----------------------------------------------------------------------------
  Root   4 singlet a              0.271514535 a.u.                7.3883 eV
  ----------------------------------------------------------------------------
     Transition Moments    X -2.27580   Y  0.00000   Z  0.00000
     Transition Moments   XX -0.00000  XY  0.00000  XZ  0.00000
     Transition Moments   YY  0.00000  YZ  0.00000  ZZ -0.00000
     Dipole Oscillator Strength                         0.93750

     Occ.   12  a   ---  Virt.   28  a     -0.05483
     Occ.   14  a   ---  Virt.   33  a     -0.06791
     Occ.   15  a   ---  Virt.   34  a     -0.06791
     Occ.   16  a   ---  Virt.   29  a     -0.09259
     Occ.   16  a   ---  Virt.   33  a     -0.07716
     Occ.   17  a   ---  Virt.   32  a     -0.10474
     Occ.   18  a   ---  Virt.   31  a     -0.10474
     Occ.   18  a   ---  Virt.   35  a      0.07663
     Occ.   20  a   ---  Virt.   23  a     -0.67952
     Occ.   21  a   ---  Virt.   22  a      0.67952
  ----------------------------------------------------------------------------
  Root   5 singlet a              0.337087837 a.u.                9.1726 eV
  ----------------------------------------------------------------------------
     Transition Moments    X  0.00000   Y -0.00000   Z  0.00000
     Transition Moments   XX -4.42936  XY -0.00001  XZ  0.00000
     Transition Moments   YY  4.42936  YZ  0.00000  ZZ  0.00000
     Dipole Oscillator Strength                         0.00000

     Occ.   13  a   ---  Virt.   34  a     -0.05712
     Occ.   15  a   ---  Virt.   35  a      0.06288
     Occ.   16  a   ---  Virt.   32  a     -0.13127
     Occ.   17  a   ---  Virt.   33  a     -0.05745
     Occ.   17  a   ---  Virt.   36  a     -0.07859
     Occ.   18  a   ---  Virt.   34  a     -0.05745
     Occ.   19  a   ---  Virt.   22  a      0.88340
     Occ.   21  a   ---  Virt.   25  a      0.39076
  ----------------------------------------------------------------------------
  Root   6 singlet a              0.337087839 a.u.                9.1726 eV
  ----------------------------------------------------------------------------
     Transition Moments    X  0.00000   Y -0.00000   Z  0.00000
     Transition Moments   XX  0.00001  XY -4.42936  XZ -0.00000
     Transition Moments   YY -0.00001  YZ -0.00000  ZZ  0.00000
     Dipole Oscillator Strength                         0.00000

     Occ.   13  a   ---  Virt.   33  a     -0.05712
     Occ.   14  a   ---  Virt.   35  a     -0.06288
     Occ.   16  a   ---  Virt.   31  a      0.13127
     Occ.   17  a   ---  Virt.   34  a     -0.05745
     Occ.   18  a   ---  Virt.   33  a      0.05745
     Occ.   18  a   ---  Virt.   36  a     -0.07859
     Occ.   19  a   ---  Virt.   23  a     -0.88340
     Occ.   20  a   ---  Virt.   25  a     -0.39076
  ----------------------------------------------------------------------------
  Root   7 singlet a              0.354630769 a.u.                9.6500 eV
  ----------------------------------------------------------------------------
     Transition Moments    X  0.00000   Y -0.00000   Z  0.00000
     Transition Moments   XX -0.00000  XY -0.00000  XZ -0.66192
     Transition Moments   YY -0.00000  YZ -0.00000  ZZ -0.00000
     Dipole Oscillator Strength                         0.00000

     Occ.   12  a   ---  Virt.   22  a      0.06165
     Occ.   14  a   ---  Virt.   23  a     -0.12883
     Occ.   15  a   ---  Virt.   22  a      0.12883
     Occ.   16  a   ---  Virt.   23  a     -0.27653
     Occ.   18  a   ---  Virt.   25  a      0.18153
     Occ.   19  a   ---  Virt.   26  a      0.21831
     Occ.   20  a   ---  Virt.   29  a      0.15119
     Occ.   21  a   ---  Virt.   24  a      0.86599
     Occ.   21  a   ---  Virt.   28  a     -0.15119
  ----------------------------------------------------------------------------
  Root   8 singlet a              0.354630771 a.u.                9.6500 eV
  ----------------------------------------------------------------------------
     Transition Moments    X  0.00000   Y -0.00000   Z  0.00000
     Transition Moments   XX  0.00000  XY -0.00000  XZ  0.00000
     Transition Moments   YY -0.00000  YZ -0.66192  ZZ  0.00000
     Dipole Oscillator Strength                         0.00000

     Occ.   12  a   ---  Virt.   23  a      0.06165
     Occ.   14  a   ---  Virt.   22  a     -0.12883
     Occ.   15  a   ---  Virt.   23  a     -0.12883
     Occ.   16  a   ---  Virt.   22  a      0.27653
     Occ.   17  a   ---  Virt.   25  a      0.18153
     Occ.   19  a   ---  Virt.   27  a      0.21831
     Occ.   20  a   ---  Virt.   24  a      0.86599
     Occ.   20  a   ---  Virt.   28  a      0.15119
     Occ.   21  a   ---  Virt.   29  a      0.15119

三重态

  ----------------------------------------------------------------------------
  Root   1 triplet a              0.074348521 a.u.                2.0231 eV
  ----------------------------------------------------------------------------
     Transition Moments                    Spin forbidden
     Oscillator Strength                   Spin forbidden

     Occ.   19  a   ---  Virt.   25  a      0.28613
     Occ.   20  a   ---  Virt.   23  a      0.67056
     Occ.   20  a   ---  Virt.   47  a     -0.06563
     Occ.   21  a   ---  Virt.   22  a      0.67056
     Occ.   21  a   ---  Virt.   48  a     -0.06563
  ----------------------------------------------------------------------------
  Root   2 triplet a              0.148608157 a.u.                4.0438 eV
  ----------------------------------------------------------------------------
     Transition Moments                    Spin forbidden
     Oscillator Strength                   Spin forbidden

     Occ.   20  a   ---  Virt.   23  a      0.70449
     Occ.   21  a   ---  Virt.   22  a     -0.70449
  ----------------------------------------------------------------------------
  Root   3 triplet a              0.148608157 a.u.                4.0438 eV
  ----------------------------------------------------------------------------
     Transition Moments                    Spin forbidden
     Oscillator Strength                   Spin forbidden

     Occ.   20  a   ---  Virt.   22  a     -0.70449
     Occ.   21  a   ---  Virt.   23  a     -0.70449
  ----------------------------------------------------------------------------
  Root   4 triplet a              0.186462043 a.u.                5.0739 eV
  ----------------------------------------------------------------------------
     Transition Moments                    Spin forbidden
     Oscillator Strength                   Spin forbidden

     Occ.   20  a   ---  Virt.   22  a     -0.70689
     Occ.   21  a   ---  Virt.   23  a      0.70689
  ----------------------------------------------------------------------------
  Root   5 triplet a              0.219201583 a.u.                5.9648 eV
  ----------------------------------------------------------------------------
     Transition Moments                    Spin forbidden
     Oscillator Strength                   Spin forbidden

     Occ.   19  a   ---  Virt.   22  a      0.78309
     Occ.   19  a   ---  Virt.   48  a     -0.06800
     Occ.   21  a   ---  Virt.   25  a      0.60402
     Occ.   21  a   ---  Virt.   49  a     -0.06391
  ----------------------------------------------------------------------------
  Root   6 triplet a              0.219201583 a.u.                5.9648 eV
  ----------------------------------------------------------------------------
     Transition Moments                    Spin forbidden
     Oscillator Strength                   Spin forbidden

     Occ.   19  a   ---  Virt.   23  a      0.78309
     Occ.   19  a   ---  Virt.   47  a     -0.06800
     Occ.   20  a   ---  Virt.   25  a      0.60402
     Occ.   20  a   ---  Virt.   49  a     -0.06391
  ----------------------------------------------------------------------------
  Root   7 triplet a              0.337007403 a.u.                9.1704 eV
  ----------------------------------------------------------------------------
     Transition Moments                    Spin forbidden
     Oscillator Strength                   Spin forbidden

     Occ.    8  a   ---  Virt.   22  a     -0.06215
     Occ.    9  a   ---  Virt.   23  a     -0.06215
     Occ.   10  a   ---  Virt.   25  a      0.08580
     Occ.   12  a   ---  Virt.   23  a     -0.06440
     Occ.   14  a   ---  Virt.   22  a      0.28514
     Occ.   15  a   ---  Virt.   23  a      0.28514
     Occ.   16  a   ---  Virt.   22  a     -0.63262
     Occ.   16  a   ---  Virt.   48  a      0.06497
     Occ.   17  a   ---  Virt.   25  a     -0.41525
     Occ.   17  a   ---  Virt.   49  a      0.05947
     Occ.   19  a   ---  Virt.   27  a     -0.08358
     Occ.   19  a   ---  Virt.   32  a     -0.10846
     Occ.   20  a   ---  Virt.   24  a     -0.24496
     Occ.   20  a   ---  Virt.   28  a     -0.25097
     Occ.   20  a   ---  Virt.   34  a      0.10173
     Occ.   21  a   ---  Virt.   29  a     -0.25097
     Occ.   21  a   ---  Virt.   33  a     -0.10173
     Occ.   21  a   ---  Virt.   36  a     -0.05781
  ----------------------------------------------------------------------------
  Root   8 triplet a              0.337007407 a.u.                9.1704 eV
  ----------------------------------------------------------------------------
     Transition Moments                    Spin forbidden
     Oscillator Strength                   Spin forbidden

     Occ.    8  a   ---  Virt.   23  a     -0.06215
     Occ.    9  a   ---  Virt.   22  a      0.06215
     Occ.   11  a   ---  Virt.   25  a      0.08580
     Occ.   12  a   ---  Virt.   22  a     -0.06440
     Occ.   14  a   ---  Virt.   23  a      0.28514
     Occ.   15  a   ---  Virt.   22  a     -0.28514
     Occ.   16  a   ---  Virt.   23  a      0.63262
     Occ.   16  a   ---  Virt.   47  a     -0.06497
     Occ.   18  a   ---  Virt.   25  a     -0.41525
     Occ.   18  a   ---  Virt.   49  a      0.05947
     Occ.   19  a   ---  Virt.   26  a     -0.08358
     Occ.   19  a   ---  Virt.   31  a      0.10846
     Occ.   20  a   ---  Virt.   29  a     -0.25097
     Occ.   20  a   ---  Virt.   33  a     -0.10173
     Occ.   20  a   ---  Virt.   36  a      0.05781
     Occ.   21  a   ---  Virt.   24  a     -0.24496
     Occ.   21  a   ---  Virt.   28  a      0.25097
     Occ.   21  a   ---  Virt.   34  a     -0.10173

所以对于第一问有：
得到了8个激发态，其激发能如下：
单重态：

  1. 5.3255 eV
  2. 5.4707 eV
  3. 7.3883 eV
  4. 7.3883 eV
  5. 9.1726 eV
  6. 9.1726 eV
  7. 9.6500 eV
  8. 9.6500 eV

三重态：
9. 2.0231 eV
10. 4.0438 eV
11. 4.0438 eV
12. 5.0739 eV
13. 5.9648 eV
14. 5.9648 eV
15. 9.1704 eV
16. 9.1704 eV

（c）

在对称性相同的激发态之间，电偶极矩不会发生变化，因此它们之间的跃迁不会导致偶极跃迁。

在本题中，单重态的 1、2、5、6、7、8 号激发态的激发能与基态有较大的能量差异，因此它们的跃迁对应的波长可能在紫外光区域。而三重态的 2、4、5、6、7、8 号激发态的激发能与基态有较小的能量差异，因此它们的跃迁对应的波长可能在可见光区域。

因为只有在电偶极矩发生变化的激发态之间才会发生光激发导致的电偶极跃迁，因此在单重态中只有第 3 和第 4 号激发态及在三重态中只有第 2 号激发态满足电偶极跃迁选择定则。这是因为这些激发态的电偶极矩与基态不同，可能发生光激发导致的偶极跃迁。

铭....
关注
  • 2
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
专栏目录
nwhess2freq:从nwchem生成的粗麻布计算频率/模式-开源
04-29
为此使用单独的脚本的好处是，您可以使用相同的粗麻布来计算不同同位素取代的频率，即，您只需运行一次nwchem计算。 从功能上讲，脚本以其当前形式完成。 随时进行修改，翻译，复制，窃取等。有关说明，请参见代码...
探索高效计算化学：深入理解NWCHEM项目
最新发布
gitblog_00057的博客
04-19 423
探索高效计算化学：深入理解NWCHEM项目 项目地址:https://gitcode.com/nwchemgit/nwchem NWCHEM是一个开源的分子模拟软件包，专为高性能计算而设计。它提供了广泛的量子化学和分子力学方法，使得研究人员能够以精确的方式理解和预测分子行为。本文将深入解析其技术特性、应用范围及优势，旨在鼓励更多的科学家和工程师利用NWCHEM进行研究工作。 技术分析 NWCHEM...
nwchem 教程
资源备忘
07-27 2633
https://wiki.alcf.anl.gov/images/b/bc/NWChem_tutorial.pdf
电子材料与器件计算方法基础——学习使用nwchem——homework3_1
weixin_46538207的博客
05-19 634
主要为朋友写的，可能用处没有那么大，聊为参考
电子材料与器件计算方法基础——学习使用nwchem——homework3_3
weixin_46538207的博客
05-31 150
为朋友写的，对其他人可能作用不大，仅为参考
nwchem （ROCM版）编译 -最终目标
jslove1997的博客
02-23 5245
1.编译ucx版本 参考网址 https://github.com/openucx/ucx/wiki/Build-and-run-ROCM-UCX-OpenMPI
nwbas2ecce:nwbas2ecce将nwchem基础集文件转换为ECCE格式-开源
04-29
此外，nwchem附带的许多较新的基础集也不包括在内，例如def2-系列。 nwbas2ecce是一个python脚本，旨在尽可能自动地将nwchem基础集文件转换为ECCE以便将其包含在基础集工具中的文件中。 该过程不是完全自动的，但我...
aaaaaa.rar_计算化学
09-23
例如，使用诸如Gaussian或NWChem等量子化学软件生成的输入和输出数据可以通过MATLAB进行后处理。 2. **分子动力学模拟**：MATLAB可以构建和运行简单的分子动力学模型，跟踪分子系统随时间的演化，以研究热运动、...
nwchempy:使用python解析nwchem输出
05-26
nwchem是一款强大的分子模拟软件，广泛应用于量子化学计算，它提供了多种方法来解决化学和材料科学中的复杂问题。nwchempy则为用户提供了方便的接口，使科学家和工程师能够更轻松地分析nwchem的计算结果。 nwchempy...
Nwchem-6.6.revision27746-src.2015-10-20.tar.gz
01-20
Nwchem (Linux)
chemPackage：创建用于操作NWChem输入和输出文件的Python包
02-28
化学包装 创建Python包以处理NWChem输入和输出文件 依存关系 截至目前，像numpy这样的常规包之外的唯一依赖项是natsort 。 可以通过pip或conda找到并安装此软件包 附加条款 该软件包是高度实验性的。 甚至在我从python2快速攻入python3之前，我都会建议任何人验证它给出的结果是否正确。
DFT的matlab源代码-NWCHEM-examples:您身处扭曲的小福特堡垒之迷，一切都一样
05-26
DFT的matlab源代码NWCHEM-示例 您身处扭曲的小福特堡垒迷宫，所有这些都一样。 NWCHEM（）是功能强大的开源电子结构包。 它的独特卖点是 Tensor压缩引擎（TCE）具有广泛的后HF方法（包括时髦的近似方法，例如CC2），可以使用全局数组（分布式内存）来扩展到超级计算机上一些非常疯狂的大小的计算； 在同一代码中进行平面波和高斯基集计算的能力。 另外，它在开源软件包中具有最大范围的异常DFT功能，包括约束DFT（CDFT）和范围分隔的混合动力（如BNL）。 但是，文档很简洁，在互联网上很少能找到示例。 在这里，我将整理多年来使用过的（.nw）脚本示例，以及vim语法高亮显示等内容，这些内容我已经从扔掉的Gist和电子邮件中删除了。 如果您在这里用Google搜索，则可能会遇到NWCHEM问题。 祝你好运！ 有用的参考 -2012年，非常有用的例子 -2010，此处有很多详细信息； 有关TCE，内存管理等的警告
NWChem Planning
资源备忘
12-24 762
NWChem Planning NWChem Future Strategy Planning Below are the pictures of the planning board of our first meeting on May 2nd. Activities that were going to be pursued: Edo: Internode commu
NWChem Usage Instructions
资源备忘
08-10 753
NWChem Usage Instructions Overview NWChem is a computational chemistry package that is designed to run on high-performance parallel supercomputers as well as conventional workstation clusters. It
NWChem
资源备忘
08-18 3467
NWChem Jump to: navigation, search Contents  [hide]  1 Overview 1.1 Official Documentation1.2 Support1.3 Tutorials 2 NWChem on Blue Gene/Q 2.1 Submission 3 NWChem on
赝势平面波计算机软件,NWChem
weixin_36140358的博客
07-09 512
NWChem是运行在高性能并行超级计算机和通常工作站集群上的计算化学软件，可以用在大多数计算平台上。NWChem使用标准量子力学描述电子波函或密度，计算分子和周期性系统的特性，还可以进行经典分子动力学和自由能模拟。外文名NWChem类别计算化学软件运行平台高性能并行超级计算机功能经典分子动力学和自由能模拟NWChem分子的电子结构编辑语音1) 下面的方法计算原子坐标的能量，解析一阶...
监控系统zabbix报错常见错误处理
gzytester的博客
03-02 1242
zabbix server is nor running 1.启动后无进程： 9月 29 16:24:47 localhost.localdomain systemd[1]: Starting SYSV: Starts and stops Zabbix Server using chkconfig... 9月 29 16:24:47 localhost.localdomain zabbix_server[12859]: Starting zabbix_server: /usr/local/zabbix
在组内集群编译nwchem
jslove1997的博客
02-20 1862
1. 编译openblas 这是nwchem编译需要使用的库，从官网上下载源代码。 #编译 make #安装到指定路径 make install PREFIX=/home/jrf/tools/openblas 2.编译ucx 要使用的OpenMPI为4.0.1版本 所以使用ucx连接。 最开始下载和编译ucx通讯模块。 此处参考 https://github.com/openucx/ucx/wi...
NWCHEM 6 make
资源备忘
11-23 1795
最近要算多组态，可是Gaussian的代码太烂了，Molpro又买不起，决定用NWChem。网上似乎没有什么编译的详细步骤，我摸索了一下，写了这个脚本，根据你自己的集群，修改一些参数，直接在src目录下运行就可以了。[code]#!/bin/csh# Author: coolra
3. 建立络合物的电子结构模型
02-07
为了建立络合物的电子结构模型，需要使用计算机软件或理论计算的方法来解决 Schrödinger 方程，以确定络合物中电子的能量和电子密度分布。常用的软件包括 Gaussian、NWChem 和 Molpro 等。 建立络合物的电子结构...
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值