24小时热门版块排行榜    

查看: 1561  |  回复: 9

zhangying09

铜虫 (小有名气)

[交流] 【求助】优化散射区和中间区域的原子 已有4人参与

各位前辈,我想问一下在ATK里怎么优化散射区和中间区域的原子,让他们之间的接触距离更加合理
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kaypu

铁杆木虫 (小有名气)



小木虫(金币+0.5):给个红包,谢谢回帖交流
上次我也问这个问题了,没有回复的 同求这个问题 我就是里面的几何限制不知道该如何设置
2楼2010-06-02 07:38:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+2):多谢专家指导 2010-06-02 08:25:27
还没有一个非常合适的方法。

本人曾在此版回答过类似的问题。

另一种方法,将它转换为equlivalent bulk system of a two-probe system,然后采用只有优化c轴的方法来优化结构。
3楼2010-06-02 08:24:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangying09

铜虫 (小有名气)

引用回帖:
Originally posted by valenhou001 at 2010-06-02 08:24:25:
还没有一个非常合适的方法。

本人曾在此版回答过类似的问题。

另一种方法,将它转换为equlivalent bulk system of a two-probe system,然后采用只有优化c轴的方法来优化结构。

那怎样在转化到TWO—prob电极系统呢??
还有就是看一篇文章上是用material studio 优化的,我想问一下版主可知道怎么优化呢
4楼2010-06-02 09:02:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
gavinliu7390(金币+1):谢谢交流! 2010-06-02 16:34:41
Virtual nanolab是可以直接输出equlivalent bulk system of a two-probe system。
pwscf里面是可以直接只优化c轴方向的。对castep还不是非常清楚怎么操作。
5楼2010-06-02 12:50:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangying09

铜虫 (小有名气)

引用回帖:
Originally posted by valenhou001 at 2010-06-02 12:50:23:
Virtual nanolab是可以直接输出equlivalent bulk system of a two-probe system。
pwscf里面是可以直接只优化c轴方向的。对castep还不是非常清楚怎么操作。

我直接指定一个距离,不优化可以吗
6楼2010-06-02 14:34:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
gavinliu7390(金币+1):谢谢交流! 2010-06-02 16:34:55
可以,但是要做很多点。比较麻烦。
7楼2010-06-02 16:33:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangying09

铜虫 (小有名气)

做很多点什么意思啊,是不是一步步改变散射区和中间区域的距离,但是最优点的判据是什么呢???还有这事我做的优化文件贴出来大家看看

from ATK.TwoProbe import *
import ATK
ATK.setVerbosityLevel(2)
# Restore old two-probe configuration
vnl_file = VNLFile("carbon6.vnl"
configurations = vnl_file.readAtomicConfigurations()
ini_two_probe_conf = configurations["carbon6"]
# Create parameters for electrodes
electrode_bz = brillouinZoneIntegrationParameters( (5,5,100) )
electrode_occ = eigenstateOccupationParameters(
temperature = 1300*Kelvin
)
electrode_params = ElectrodeParameters(
eigenstate_occupation_parameters = electrode_occ,
brillouin_zone_integration_parameters = electrode_bz
)
# Collect parameters into a two-probe calculation method
two_probe_method=TwoProbeMethod(
electrode_parameters = (electrode_params, electrode_params),
exchange_correlation_type = LDA.PZ
)
# Define parameters for geometry optimization
opt_params = geometricOptimizationParameters(
force_tolerance = 0.05*eV/Ang,
max_steps = 200
)

# Optimize geometry for the two-probe
opt_two_probe_conf = calculateOptimizedAtomicGeometry(
atomic_configuration = ini_two_probe_conf,
method = two_probe_method,
optimization_parameters = opt_params,
)
# Print and export two-probe configuration to a VNL file
print '\n','='*50,'\n Two-probe optimized configuration:\n','='*50
nlPrint(opt_two_probe_conf)
vnl_file = VNLFile("carbon6opt"
vnl_file.addToSample(opt_two_probe_conf, "carbon6opt

但是有一个问题是,我的中间区域在material studio 里已经优化好了,这样再优化不知道是不是合理
8楼2010-06-02 16:47:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangying09

铜虫 (小有名气)

引用回帖:
Originally posted by valenhou001 at 2010-06-02 16:33:04:
可以,但是要做很多点。比较麻烦。

做很多点什么意思啊,是不是一步步改变散射区和中间区域的距离,但是最优点的判据是什么呢???还有这事我做的优化文件贴出来大家看看

from ATK.TwoProbe import *
import ATK
ATK.setVerbosityLevel(2)
# Restore old two-probe configuration
vnl_file = VNLFile("carbon6.vnl"
configurations = vnl_file.readAtomicConfigurations()
ini_two_probe_conf = configurations["carbon6"]
# Create parameters for electrodes
electrode_bz = brillouinZoneIntegrationParameters( (5,5,100) )
electrode_occ = eigenstateOccupationParameters(
temperature = 1300*Kelvin
)
electrode_params = ElectrodeParameters(
eigenstate_occupation_parameters = electrode_occ,
brillouin_zone_integration_parameters = electrode_bz
)
# Collect parameters into a two-probe calculation method
two_probe_method=TwoProbeMethod(
electrode_parameters = (electrode_params, electrode_params),
exchange_correlation_type = LDA.PZ
)
# Define parameters for geometry optimization
opt_params = geometricOptimizationParameters(
force_tolerance = 0.05*eV/Ang,
max_steps = 200
)

# Optimize geometry for the two-probe
opt_two_probe_conf = calculateOptimizedAtomicGeometry(
atomic_configuration = ini_two_probe_conf,
method = two_probe_method,
optimization_parameters = opt_params,
)
# Print and export two-probe configuration to a VNL file
print '\n','='*50,'\n Two-probe optimized configuration:\n','='*50
nlPrint(opt_two_probe_conf)
vnl_file = VNLFile("carbon6opt"
vnl_file.addToSample(opt_two_probe_conf, "carbon6opt

但是有一个问题是,我的中间区域在material studio 里已经优化好了,这样再优化不知道是不是合理
9楼2010-06-02 16:48:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

youngjfly

银虫 (初入文坛)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
aylayl08(金币+2):谢谢指点 2010-09-29 19:05:14
引用回帖:
Originally posted by zhangying09 at 2010-06-02 16:48:41:



做很多点什么意思啊,是不是一步步改变散射区和中间区域的距离,但是最优点的判据是什么呢???还有这事我做的优化文件贴出来大家看看

from ATK.TwoProbe import *
import ATK
ATK.setVerbosityLevel ...

因为优化好中心分子了,所以ATK的自动优化有可能会破坏你以前的结构,教程里面的两个H原子就跑偏了。所以还是手动调节距离比较稳妥,只要计算total energy,取最小的能量就行。本人愚见
真的很想把你搞好。
10楼2010-09-29 16:26:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhangying09 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:33 by ZPa0EcMwuECS
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:49 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:47 by ZPa0EcMwuECS
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 01:36 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 01:03 by ZPa0EcMwuECS
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-30 00:59 by winsaint
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 00:52 by ZPa0EcMwuECS
[基金申请] 中青基了要发朋友圈吗? +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[基金申请] 为什么到现在没收到通知? +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见