1. 简介

   基于微扰理论的态求和法（sum-over-states），我们编写了适用于分子（最高5阶）和固体（最高3阶）的非线性光学极化率和吸收系数（最多3光子吸收）的计算程序。

2. 下载（请先登陆） sos_dist240615.rar

3. 使用

   使用之前请熟读相应的参考文献：

   （1）分子：JCP_1994v100p6535_Explicit nondivergent formulas for atomic and molecular dynamic hyperpolarizabilities。

   准备文件sos.inp 和 *.mnt 输入文件，详见example目录H2O.mnt例子

   H2O.mnt可借助multiwfn程序获得。

   （2）固体：PRB_1995v52p14636_Nonlinear optical susceptibilities of semiconductors。

   PRB_2005v72p045223_Scissors implementation within length-gauge formulations of NLO。

   PRB_1993v48p11705_Nonlinear optical response of semiconductors in the IPA。

   JOSAB_1985v2p294_Review of multiphoton absorption in crystalline solids。

          a) 支持yambo5.1.1: 准备文件sos.inp即可，yambo程序修改后自动会生成必需文件。

          b) 支持MS的castep: 准备文件sos.inp 和castep光学计算结果中的*_Optics.bands，*_Optics.cell和*_Optics.cst_ome文件，详见example目录中的例子。

   （3）sos.inp 文件的内容请按里面的提示输入。准备好所有的文件后，放在同一目录，运行sos.exe文件即可。建议打开ms-dos窗口运行sos.exe，不要直接鼠标双击运行，否则如有错误提示，会一闪而过。linux版本的运行类似。

4. 致谢

   （1）输入文件的设计参考elk程序http://elk.sourceforge.net/ 

   （2）分子SOS计算设计（并行部分和输入参数格式）参考multiwfn程序http://sobereva.com/multiwfn/ 

   （3）光学矩阵元素的计算依赖Material Studio 程序 http://accelrys.com/products/collaborative-science/biovia-materials-studio/ 或 Yambo程序 https://www.yambo-code.eu/ 

   （4）对称性设计参考http://www.abinit.org/ 

5. 免责申明

   本程序完全免费，因使用本程序导致的一切问题，请自行承担，暂不公布源码。为避免重复回答问题，请有问题的用户在本贴评论中提问，其它无关的评论一概删除，谢谢支持！

6. 引用

   使用本程序发布结果，建议引用参考文献：

   （1） SHG: https://wulixb.iphy.ac.cn/article/doi/10.7498/aps.67.20181337 

   （2） THG: http://dx.doi.org/10.1016/j.cocom.2016.07.002  

7. 其它说明

    (1) 由于本人时间有限，只提供软件使用的技术支持，无法进行合作。

    (2) 工作太忙了，写程序累又花时间，请予理解，谢谢！

    (3)“支持自旋轨道耦合和自旋极化计算”是指支持QE+Yambo或CASTEP的“自旋轨道耦合和自旋极化计算”结果输出作为输入。

    (4) 建议大家使用Yambo，这个软件的技术支持“及时”和“详细”，同时在计算光学矩阵元的时候它也考虑了[r,V_nl]，理论上比MS的准确。

    (5) 你只要对非线性光学理论有初步的了解（找本书看下，建议《非线性光学物理》[叶佩弦]），就能看懂sos.inp里的参数，初步使用SOS是没问题的。

    (6) 本人目前主要研究原生和含缺陷的二维过渡金属硫属化合物和二维氮化硼体系的线性和非线性光学性能。欢迎交流。

    请不要在此评论，因为现在博客48小时后就不能修改了，我是留作增加附加信息用的。谢谢支持！