计算化学公社

标题: 在siesta中怎样输入zn和co 的lda+u ？ [打印本页]
作者
Author: ZHANGZY    时间: 2016-6-23 20:48
标题: 在siesta中怎样输入zn和co 的lda+u ？
LDAU.FirstIteration .true.
LDAU.PopTol 5.0d-4
LDAU.ThresholdTol 1.0d-2
LDAU.ProjectorGenerationMethod 2
%block LDAU.proj
Co 1 # number of shells of projectors
n=3 2 # n, l
3.00 0.0000 # U(eV), J(eV)
0.000 0.0000 # rc, \omega
%endblock LDAU.proj

co原子的lda+u输入，对吗？
作者
Author: 卡开发发    时间: 2016-6-23 22:51
https://www.mail-archive.com/siesta-l@uam.es/msg08831.html这里提到了这件事：
the parameters to run with LDA+U are fully described in the Docs/siesta.tex file, Section 8.

以及这里提到的https://www.mail-archive.com/siesta-l@uam.es/msg07291.html输出文件记录了一些
LDAU.FirstIteration .false.

LDAU.PopTol 5.0d-4

LDAU.ThresholdTol   1.0d-2

LDAU.ProjectorGenerationMethod 2

%block LDAU.proj

Fe   1           # number of shells of projectors

n=3    2         #  n, l

   2.00 0.0000    # U(eV), J(eV)

   0.000 0.0000  # rc, \omega

%endblock LDAU.proj

trunk的版本我没用过，我不确定你的是不是trunk308的版本，要是是的话我顺带求一份好了，这个现在很难下载。
作者
Author: 卡开发发    时间: 2016-6-24 13:00
从朋友@wei那边要来了一份siesta-ldau的，我把siesta.tex编译然后转换成pdf发上来，里面写的比较清楚，你可以研究一下（说白了也就那些参数）。

等有时间我看看能不能把线性响应确定U的方法整进去，前段时间正好对这块做了些了解。
作者
Author: ZHANGZY    时间: 2016-6-26 06:57
多谢，我这里有siesta-4-b2版本，也包含lda+U功能。运行可以。他的说明书里面没有lda+u，但是功能包含。
我编译siesta-ldau-308.tgz没有成功。不知道错在哪里
作者
Author: 卡开发发    时间: 2016-6-26 14:32
本帖最后由 卡开发发 于 2016-6-26 16:38 编辑
ZHANGZY 发表于 2016-6-26 06:57
多谢，我这里有siesta-4-b2版本，也包含lda+U功能。运行可以。他的说明书里面没有lda+u，但是功能包含。
...

按官网的说法，只有这个308ldau也就是2.67的trunck版才支持。我看4.0b2的手册，确实没有ldau的描述。308这个版本我在OSX下用GNU编译器确实会出错（相同环境下编译过2.02、3.2、4.0b），但是换到RHEL5.7上用ifort11编译就没有问题，你可以试试用intel编译器，mpi用intel或者mpich2试试，还不行的话就make >> make.log，然后把这个make.log上传上来看看。
作者
Author: yasuo    时间: 2017-4-29 20:43
你好  你说4.0-b2有lda+u功能  怎么加u了和没加效果一样啊   4.1-b2编译成功了不能计算
作者
Author: yasuo    时间: 2017-4-29 20:49
ZHANGZY 发表于 2016-6-26 06:57
多谢，我这里有siesta-4-b2版本，也包含lda+U功能。运行可以。他的说明书里面没有lda+u，但是功能包含。
...

你好 你有siesta-4-b2版本加u的算例吗？求发一个  1056118571@qq.com
作者
Author: 卡开发发    时间: 2017-4-30 00:01
本帖最后由 卡开发发 于 2017-4-30 00:03 编辑
yasuo 发表于 2017-4-29 20:49
你好 你有siesta-4-b2版本加u的算例吗？求发一个

现在siesta-4.1b2就已经支持+U了，在这个版本的test中的mno的例子就是dft+U的。编译的时候换intel2015的compiler 和openmpi试试。
作者
Author: yasuo    时间: 2017-5-2 01:39
卡开发发 发表于 2017-4-30 00:01
现在siesta-4.1b2就已经支持+U了，在这个版本的test中的mno的例子就是dft+U的。编译的时候换intel2015的c ...

你好，siesta4.1-b2确实有lda+u的功能，我下载了，按照4.0及以前的版本用gfortran和mpich2编译也能成功，但是计算的时候会报错，是不是一定要用ifort编译，能不能给我看看您的arch.make啊，我的邮箱1056118571@qq.com。谢谢了
作者
Author: ZHANGZY    时间: 2017-5-2 15:19
yasuo，你好，附件有两个arch.make 一个是用gfortran一个使用ifort。
如果你用gfortran，一定要按照手册把相关的动态库全部装一遍。我用的是ubuntu系统，晕乎乎的装上了，试过例子可以。
最好用ifort，据说快不少。
我也是学siesta时间不长，可以互相讨论。
作者
Author: yasuo    时间: 2017-5-2 15:26
好的，谢谢了，我看看，有问题再请教大神
作者
Author: yasuo    时间: 2017-5-2 15:49
siesta 4.1-b2编译出来没有版本号。有的能算，有的不能算。
                           ***********************      
                           *  WELCOME TO SIESTA  *      
                           ***********************      

reinit: Reading from standard input
reinit: Dumped input in INPUT_TMP.74629
************************** Dump of input data file ****************************
作者
Author: 卡开发发    时间: 2017-5-2 17:13
yasuo 发表于 2017-5-2 15:49
siesta 4.1-b2编译出来没有版本号。有的能算，有的不能算。
                           *************** ...

看下test能否正常进行，也有可能是fdf写的有问题



欢迎光临 计算化学公社 (http://ccc.keinsci.com/) Powered by Discuz! X3.3