快捷下单入口 关于 合作 招聘 新人手册 会员中心

    热线:400-152-6858

    测试狗科研服务

    预存 免费试测 登录
    Document
    当前位置:文库百科文章详情
    ASE原子模拟环境入门之纳米颗粒
    来源: 时间:2022-12-16 16:42:25 浏览:2037次

    ASE提供模块 ase.cluster,用于建立具有普通晶体形式的金属纳米粒子



    构建和优化纳米颗粒


    ase.cluster.Octahedron(),除了可生成严格的八面体纳米粒子外,它还提供了一个cutoff关键字来切割八面体的角。这就产生了“截短的八面体”,这是纳米粒子中一个众所周知的结构基序。此外,晶格将与银的大块FCC结构一致。

    练习:尝试使用ase.cluster.Octahedron()来制造不同形态的截短八面体。构建一个含有55个原子的银的14面体颗粒。像往常一样,用ASE GUI来将它可视化。

    ASE提供了一个基于有效介质理论的力场代码,ase.calculators.emt.EMT主攻FCC金属(Cu,  Ag, Au, Pt等)。这比DFT快得多,所以我们用它来优化我们的立方体八面体。




    基态


    金属纳米粒子最有趣的问题之一是它们的电子结构和其他性质如何依赖于尺寸。一个小的纳米粒子就像一个分子,只有几个离散能级。一个大的纳米颗粒则像块体材料,具有连续的态密度。我们计算一下纳米粒子的Kohn-Sham能谱(以及态密度)

    我们设置了一些参数来节省时间,因为这仅是范例计算。我们想要一个更小的基集,也想要一个比正常情况下电子更少的PAW数据集。我们还想使用费米弥散函数,因为在费米能级附近可能有多个电子态:

    
    

    from gpaw import GPAW, FermiDirac

    calc = GPAW(mode='lcao', basis='sz(dzp)', setups={'Ag': '11'},occupations=FermiDirac(0,1))

    这里有一些特定于gpaw的关键字——如果使用其他代码,这些变量将有其他名称。

    练习:使用GPAW运行一个优化过的Ag55 颗粒的单点计算。

    计算之后,将基态存储到文件中:

    calc.write('groundstate.gpw')




    态密度


    一旦我们保存了.gpw文件,我们可以用它写一个新的脚本来加载它并获得DOS:

    import matplotlib.pyplot as plt

    from gpaw import GPAW

    calc = GPAW('groundstate.gpw')

    energies, dos = calc.get_dos(npts=500, width=0.1)

    efermi = calc.get_fermi_level()

    在这个例子中,我们对DOS进行采样的时候使用了0.1 eV的高斯展宽。你需要在图中标记费米能级。一个好的方法是画一条垂直线:plt.axvline(efermi)

    练习:使用matplotlib画出DOS,并标出费米能级。

    练习:观察这个DOS图,它应该被理解为离散的还是连续的?

    这张图应该告诉我们,已经有55个原子,大量的d电子已经形成一个连续的带(回想一下,我们使用的是0.1 eV高斯展宽)。与此同时,少数s电子的能量在更大范围内离散分布,我们清楚地看到孤立的峰值:s态仍然清晰地量子化并有明显的间隙。贵金属Cu, Ag和Au的特征是,它们的d能带被完全占据,所以费米能级位于这些s态之间。电子数不同的团簇可能具有更高或更低的费米能级,这强烈地影响了它们的反应性。我们可以推测,在55个原子时,独立Ag纳米粒子的性质可能强烈地依赖于大小。

    上述分析是推测性的。为了验证这个分析,我们需要计算s、p和d轨道的投影DOS,看看我们的假设是否正确。




    练习答案


    优化立方八面体:

    from ase.cluster import Octahedron

    from ase.calculators.emt import EMT

    from ase.optimize import BFGS

    atoms = Octahedron('Ag', 5, cutoff=2)

    atoms.calc = EMT()

    opt = BFGS(atoms, trajectory='opt.traj')

    opt.run(fmax=0.01)


    计算基态:

    from gpaw import GPAW, FermiDirac

    from ase.io import read

    atoms = read('opt.traj')

    calc = GPAW(mode='lcao', basis='sz(dzp)', txt='gpaw.txt',

               occupations=FermiDirac(0.1),

               setups={'Ag': '11'})

    atoms.calc = calc            

    atoms.center(vacuum=4.0)

    atoms.get_potential_energy()

    atoms.calc.write('groundstate.gpw')


    画DOS图:

    import matplotlib.pyplot as plt

    from gpaw import GPAW

    from ase.dft.dos import DOS

    calc = GPAW('groundstate.gpw')

    dos = DOS(calc, npts=800, width=0.1)

    energies = dos.get_energies()

    weights = dos.get_dos()

    ax = plt.gca()

    ax.plot(energies, weights)

    ax.set_xlabel(r'$E - E_{\mathrm{Fermi}}$ [eV]')

    ax.set_ylabel('DOS [1/eV]')

    plt.savefig('dos.png')

    plt.show()


    计算狗·模拟计算

    评论 / 文明上网理性发言
    12条评论
    全部评论 / 我的评论
    最热 /  最新
    全部 3小时前 四川
    文字是人类用符号记录表达信息以传之久远的方式和工具。现代文字大多是记录语言的工具。人类往往先有口头的语言后产生书面文字,很多小语种,有语言但没有文字。文字的不同体现了国家和民族的书面表达的方式和思维不同。文字使人类进入有历史记录的文明社会。
    点赞12
    回复
    全部
    查看更多评论
    相关文章

    【科研干货】电化学表征:循环伏安法详解(上)

    2019-10-25

    晶体结构可视化软件 VESTA使用教程(下篇)

    2021-01-22

    手把手教你用ChemDraw 画化学结构式:基础篇

    2021-06-19

    【科研干货】电化学表征:循环伏安法详解(下)

    2019-10-25

    基础理论丨一文了解XPS(概念、定性定量分析、分析方法、谱线结构)

    2020-05-03

    Zeta电位的基本理论、测试方法和应用

    2020-08-24

    项目推荐/Project
    第一性原理-锂硫路径

    第一性原理-锂硫路径

    第一性原理-形成能

    第一性原理-形成能

    热门文章/popular

    【科研干货】电化学表征:循环伏安法详解(上)

    晶体结构可视化软件 VESTA使用教程(下篇)

    手把手教你用ChemDraw 画化学结构式:基础篇

    电化学实验基础之电化学工作站篇 (二)三电极和两电极体系的搭建 和测试

    【科研干货】电化学表征:循环伏安法详解(下)

    基础理论丨一文了解XPS(概念、定性定量分析、分析方法、谱线结构)

    微信扫码分享文章

    意见反馈

    有奖举报

    商务合作

    ...

    更多

    公众号

    关注我们 了解更多

    小程序

    随时预约 掌握进度

    举报有奖

    TEL: 191-3608-6524

    如:在网络上恶意使用“测试狗”等相关关键词误导用户点击、恶意盗用测试狗商标、冒称官方工作人员等情形,请您向我们举报,经查实后,我们将给予您奖励。

    举报内容:

    200

    上传附件:
    文件格式不正确,请重新上传文件格式不正确,请重新上传文件格式不正确,请重新上传
    文件格式:jpg、jpeg、png、gif、tif、doc、docx、ppt、pptx、xls、xlsx、pdf、zip、rar
    联系方式
    姓名
    电话
    提交意见

    意见反馈

    Suggestions

    您可以在此留下您宝贵的意见,您的意见或问题反馈将会成为我们不断改进的动力。

    意见类型
    测试服务
    网站功能
    财务报账
    其他类型
    意见内容

    200

    联系方式
    姓名
    电话
    提交意见

    收起

    01

    专属信用额度,先测后付0元下单

    02

    下单享高额积分,万千好礼免费兑

    200
    200元无门槛优惠券
    立即激活 立即下单
    已使用
    已作废

    全流程在线可视化,便捷高效触手可及

    如下单过程中有任何疑问或需要帮助,请随时咨询专属顾问~
    9
    9折无门槛优惠券

    支付一笔订单后可领取

    立即领取 立即下单
    已使用
    已作废

    免费测+惊喜盲盒+高额福利,多重福利大放送

    新人免费测

    双双开盲盒(100%中奖)

    邀请人专享

    受邀人专享

    创建/加入团队,解锁定制化权益

    01

    1500元团队专属优惠券

    02

    万元大额信用额度,享先测后付

    03

    团队成员统一开票报销;

    04

    专业工程师课题专属服务

    领取成功,请下单
    请您支付一笔订单后才可以领取优惠券
    Document
    关于我们 新手帮助 测试干货 商务合作 基金查询 相关资质 模拟计算 现场测试 服务项目 科研绘图 同步辐射 电池行业

    联系方式/contact

    400-152-6858

    工作时间/work time

    09:00-18:00

    测试狗公众号

    关注我们 了解更多

    测试狗小程序

    随时预约 掌握进度

    蜀公网安备51010602000648号

    蜀ICP备17005822号-1

    成都世纪美扬科技有限公司

    Copyright@测试狗·科研服务