| 查看: 680 | 回复: 3 | ||
yourne木虫 (小有名气)
|
[求助]
用ATK初步计算报错
|
|
使用new calculator,，选ATK-DFT,版本11.2.2，每次都立即出现如下报错： ** Back Engine Exception 541 : Diagonalization error, overlap matrix not positive definite, try to increase interaction_max_range ** Location : mathutils.cpp:1062 之前都是按照VNLtutorials来的，材料是TiB2，求助如何解决？ |
回复此楼» 猜你喜欢
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
假如你的研究生提出不合理要求
已经有5人回复
-
所感
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有28人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
-
实验室接单子
已经有3人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
guohuazhong
至尊木虫 (职业作家)
- 1ST强帖: 1
- 应助: 52 (初中生)
- 贵宾: 0.032
- 金币: 35131.5
- 散金: 5
- 红花: 19
- 沙发: 1
- 帖子: 3924
- 在线: 885.6小时
- 虫号: 550098
- 注册: 2008-04-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2012-11-06 19:06:15
yourne
木虫 (小有名气)
- 应助: 11 (小学生)
- 金币: 1979.5
- 散金: 24
- 红花: 1
- 帖子: 237
- 在线: 169.8小时
- 虫号: 1670096
- 注册: 2012-03-06
- 性别: GG
- 专业: 金属材料表面科学与工程
|
new calculator文件： # ------------------------------------------------------------- # Bulk configuration # ------------------------------------------------------------- # Set up lattice lattice = Hexagonal(3.0292*Angstrom, 3.2284*Angstrom) # Define elements elements = [Titanium, Boron, Boron, Titanium, Titanium, Titanium, Titanium, Titanium, Titanium, Titanium] # Define coordinates cartesian_coordinates = [[ 0. , 0. , 0. ], [ 1.5146 , 0.87462961, 1.6142 ], [ 1.5146 , -0.87462961, 1.6142 ], [ 1.5146 , -2.62336415, 0. ], [ 1.5146 , 2.62336415, 0. ], [ 3.0292 , 0. , 0. ], [ 3.0292 , 0. , 3.2284 ], [ 1.5146 , -2.62336415, 3.2284 ], [ 0. , 0. , 3.2284 ], [ 1.5146 , 2.62336415, 3.2284 ]]*Angstrom # Set up configuration bulk_configuration = BulkConfiguration( bravais_lattice=lattice, elements=elements, cartesian_coordinates=cartesian_coordinates ) # ------------------------------------------------------------- # Calculator # ------------------------------------------------------------- #---------------------------------------- # Basis Set #---------------------------------------- basis_set = GGABasis.DoubleZeta #---------------------------------------- # Exchange-Correlation #---------------------------------------- exchange_correlation = GGA.PBE numerical_accuracy_parameters = NumericalAccuracyParameters( grid_mesh_cutoff=300.0*eV, k_point_sampling=(3, 3, 3), ) calculator = LCAOCalculator( basis_set=basis_set, exchange_correlation=exchange_correlation, numerical_accuracy_parameters=numerical_accuracy_parameters, ) bulk_configuration.setCalculator(calculator) bulk_configuration.update() nlsave('F:/ykc/TiB2.nc', bulk_configuration) # ------------------------------------------------------------- # Bandstructure # ------------------------------------------------------------- bandstructure = Bandstructure( configuration=bulk_configuration, route=['G', 'M', 'L', 'A', 'G', 'K', 'H', 'A'], points_per_segment=50, bands_above_fermi_level=All ) nlsave('F:/ykc/TiB2.nc', bandstructure) 接触ATK时间很短，请勿见笑。。 |
3楼2012-11-06 20:40:19
kaypu 
铁杆木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 399
- 散金: 130
- 帖子: 292
- 在线: 113.4小时
- 虫号: 864391
- 注册: 2009-10-07
- 专业: 凝聚态物性 II ：电子结构
4楼2012-12-04 21:15:30