如何用g09如何获得Bq的受力

alystone

请问各位老师，    我可以用g03的massage关键词，得到我结构中Bq的受力。但是在G09中删除了massage的此类功能，我想用punch+iop(7/32)来输出，但是也输出不出来。    请教有哪位老师知道怎么在g09中得到Bq的受力情况。g03的输入输出文件见附件，以供参考。

谢谢！

sobereva

g09试了，可以跑，结果和你传的g03一样
massage在g09里只是变成了deprecated，但还能用

alystone

sobereva 发表于 2018-1-17 18:14
g09试了，可以跑，结果和你传的g03一样
massage在g09里只是变成了deprecated，但还能用

请问sob老师用的哪个g09版本，输入文件和我的这个输入文件一模一样吗？

sobereva

alystone 发表于 2018-1-17 18:30
请问sob老师用的哪个g09版本，输入文件和我的这个输入文件一模一样吗？

一模一样
G09W D.01

alystone

sobereva 发表于 2018-1-17 19:04
一模一样
G09W D.01

谢谢sob老师。我想是因为刚才那个体系太简单了，我的体系比较大，然后计算的时候出现这样的错误，应该怎么解决。谢谢！ 具体错误见附件g09计算出的log文件。 不加massage关键词。能够算下去，但是结果是错的。加上massage关键词，就这样报错了。

sobereva

alystone 发表于 2018-1-17 20:30
谢谢sob老师。我想是因为刚才那个体系太简单了，我的体系比较大，然后计算的时候出现这样的错误，应该怎 ...

你的赝势没有被读入，搜索 Pseudopotential Parameters，会看到Fe并没有赋予赝势
把背景电荷和赝势字段交换一下看看
massage我自己倒不用。如果之后还有问题，把背景电荷只加一小点，然后逐渐增加，反复尝试找原因。

alystone

sobereva 发表于 2018-1-17 21:06
你的赝势没有被读入，搜索 Pseudopotential Parameters，会看到Fe并没有赋予赝势
把背景电荷和赝势字 ...

谢谢sob老师。的确如你所说，就是位置不对。但是不管怎么换，好像都不对。
情况1：
   坐标
   空
   基组
   空
   赝势
   空
   Bq
这种情况下，Bq读不进去。

情况2：
    坐标
   空
   基组
   空
   Bq
   空
   赝势
这种情况下，赝势基组读不进去。

请问sob老师，gaussian可以把基组和赝势信息写在外边，用引用文件的方法来避免这种情况吗。或者还有其他什么方法。再次麻烦您了。谢谢！

sobereva

alystone 发表于 2018-1-17 21:40
谢谢sob老师。的确如你所说，就是位置不对。但是不管怎么换，好像都不对。
情况1：
   坐标

引用文件并不解决问题，在计算时只不过会把引用的文件进行展开而已。
关于读取顺序，在这里有详细说明
http://sobereva.com/g09/m_input.htm

alystone

sobereva 发表于 2018-1-17 23:57
引用文件并不解决问题，在计算时只不过会把引用的文件进行展开而已。
关于读取顺序，在这里有详细说明
...

谢谢sob老师。

alystone

sobereva 发表于 2018-1-17 23:57
引用文件并不解决问题，在计算时只不过会把引用的文件进行展开而已。
关于读取顺序，在这里有详细说明
...

sob老师您好，如果我用charge代替massage用来添加背景电荷，是否有方法来输出背景电荷的受力。谢谢

sobereva

alystone 发表于 2018-1-18 08:10
sob老师您好，如果我用charge代替massage用来添加背景电荷，是否有方法来输出背景电荷的受力。谢谢

应该没办法