| 查看: 846 | 回复: 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，求助如何解决？ |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有7人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
明天应该可查了！？
已经有5人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有7人回复
-
没有任何消息-是不是就凉了
已经有8人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
人气不行了
已经有10人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
» 本主题相关商家推荐: (我也要在这里推广)
kaypu 
铁杆木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 399
- 散金: 130
- 帖子: 292
- 在线: 113.4小时
- 虫号: 864391
- 注册: 2009-10-07
- 专业: 凝聚态物性 II ：电子结构
4楼2012-12-04 21:15:30
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
20