J*******3
发帖数: 1651 | 1 本人用SIESTA. 一个general的问题。
通常情况下，设计一个体系，假设选用了Atomic orbitals 作basis，用cutoff radius
和orbital撒数来控制计算精度。这时，那些控制planwaves 和realspace grid 的参
数们都不起作用了吗？换句话说，他们对计算没什么贡献了，对吗？ |
e******g
发帖数: 29 | 2 不对
首先，这个问题跟你用不用DFT无关，即便用HF或post-HF这个问题的表现也是一样
的。
其次，siesta的具体实现里面是没有planwave的。
然后，siesta的计算精度不仅仅是由basis orbital的cutoff radii和数目确定
的，其他的参数也对精度有影响，包括realspace grid density
radius
grid 的参
【在 J*******3 的大作中提到】
: 本人用SIESTA. 一个general的问题。
: 通常情况下，设计一个体系，假设选用了Atomic orbitals 作basis，用cutoff radius
: 和orbital撒数来控制计算精度。这时，那些控制planwaves 和realspace grid 的参
: 数们都不起作用了吗？换句话说，他们对计算没什么贡献了，对吗？
|
J*******3
发帖数: 1651 | 3 为什么？ basis 的选择 通常有三个options：1） atomic orbital 2）plane wave
3） real space without basis。 他们不是平行的吗？ 我的意思是一旦用其中一个
，例如说 1， 那就是说我们用 atomic orbital 同展开 KS equation,另外两个以及相
关的参数在计算过程中就被neglected? 不对吗？ |
e******g
发帖数: 29 | 4 哦，大概明白你问的什么了。。。
在siesta中并没有关于平面波展开KS轨道的功能。也没有相关参数，所以你的这个问题
貌似不是个问题。。
所谓的real space without basis，我并不熟悉，orbital free DFT应该算这一块吧？
不过这一点也跟siesta没有关系
此外，我印象中还没见过某个软件既可以用平面波基组，也可以用原子轨道基组。如果
你知道的话请告诉我，呵呵，谢谢 |
h********r
发帖数: 821 | 5 real space without basis只是计算实现方法，orbital free dft是理论框架，这个是
两回事。
【在 e******g 的大作中提到】
: 哦，大概明白你问的什么了。。。
: 在siesta中并没有关于平面波展开KS轨道的功能。也没有相关参数，所以你的这个问题
: 貌似不是个问题。。
: 所谓的real space without basis，我并不熟悉，orbital free DFT应该算这一块吧？
: 不过这一点也跟siesta没有关系
: 此外，我印象中还没见过某个软件既可以用平面波基组，也可以用原子轨道基组。如果
: 你知道的话请告诉我，呵呵，谢谢
|
J*******3
发帖数: 1651 | 6 很困惑的问：
可是一些DFT books 中叙述说： 如果是solids（ 也就是说有周期性的块材） 通常是
用planewave 作基的？难道不是么？
如果SIESTa，中没有涉及到 planewave，也就是说不管是什么体系，SIESTa 通用
atomic orbitals 来作基的？
那样的话，其中的energycutoff 有什么用呢？ |
e******g
发帖数: 29 | 7
架，这个是
啊哈，了解了一下，不是一回事：）
【在 h********r 的大作中提到】
: real space without basis只是计算实现方法，orbital free dft是理论框架，这个是
: 两回事。
|
e******g
发帖数: 29 | 8
常是
用
在一个计算中采用什么基组跟你的体系无关，只跟你选用的软件，或者计算方法有关。
对bulk体系，平面波基组有天生的优势，但对低维材料来说，平面波基组相对于原子轨
道基组来说计算量比较大。
siesta是基于原子轨道基组的，它的参数貌似跟平面波cutoff之间没有很明显的联
系。这也是siesta缺少一种系统化提高精度的方法的原因之一。
不清楚siesta里哪个参数叫“energycutoff”
【在 J*******3 的大作中提到】
: 很困惑的问：
: 可是一些DFT books 中叙述说： 如果是solids（ 也就是说有周期性的块材） 通常是
: 用planewave 作基的？难道不是么？
: 如果SIESTa，中没有涉及到 planewave，也就是说不管是什么体系，SIESTa 通用
: atomic orbitals 来作基的？
: 那样的话，其中的energycutoff 有什么用呢？
|
J*******3
发帖数: 1651 | 9 sorry, 它不叫energycutoff, 而是meshcutoff， 定义说是 planewave cutoff for
the grid?
既然不是用planewave, 为啥还有这个参数？
另一个问题是：在SIESTA 中用 super cell 不是为了使非周期性的东东具有周期性而
利用 planewave 吗？ |
h********r
发帖数: 821 | 10 推荐看martin的书，siesta也有足够多的参考资料。
【在 J*******3 的大作中提到】
: sorry, 它不叫energycutoff, 而是meshcutoff， 定义说是 planewave cutoff for
: the grid?
: 既然不是用planewave, 为啥还有这个参数？
: 另一个问题是：在SIESTA 中用 super cell 不是为了使非周期性的东东具有周期性而
: 利用 planewave 吗？
|
e******g
发帖数: 29 | 11
cutoff for
性而
关于meshcutoff，这个问题貌似在别的帖子里回答过了，简单说就是siesta的
Hamiltonian可以分成两坨，其中双中心积分是在k空间做的，其他的在real space
做比较容易，但最后要把两块加起来。这里面需要一个Fourier变化倒来倒去。而
meshcutoff指定的是做FFT时的实空间格点密度。这个参数跟基组没有关系
siesta中会把所有的体系都当做周期性体系来算。但按你设置super cell的不同，你
的体系实际上不一定bulk，低维体系可以通过设定真空层来模拟。
程序用到了super cell 跟它用什么基组没有关系
【在 J*******3 的大作中提到】
: sorry, 它不叫energycutoff, 而是meshcutoff， 定义说是 planewave cutoff for
: the grid?
: 既然不是用planewave, 为啥还有这个参数？
: 另一个问题是：在SIESTA 中用 super cell 不是为了使非周期性的东东具有周期性而
: 利用 planewave 吗？
|
J*******3
发帖数: 1651 | |
J*******3
发帖数: 1651 | 13 推荐看martin的书? 能不能把书的全名给一下，以便查找？感激涕凌！ |
h********r
发帖数: 821 | |
