24小时热门版块排行榜    

查看: 1007  |  回复: 2

identation

金虫 (著名写手)

[交流] 【求助】请教siesta可否计算磁性体系的介电极化呢?

It has to be stressed that the macroscopic polarization can only be meaningfully calculated
using this approach for insulators. Therefore, the presence of an energy gap is necessary,
and no band can cross the Fermi level. The program performs a simple check of this
condition, just by counting the electrons in the unit cell ( the number must be even for
a non-magnetic system, and the total spin polarization must have an integer value for
spin polarized systems), however is the responsability of the user to check that the system
under study is actually an insulator (for both spin components if spin polarized).

这段话,请教siesta可不可以计算考虑自旋的介电极化呢?
做的是磁性材料的介电性能计算,要考虑自旋,pwscf没法计算自旋体系的电极化,因为磁性体系一般要用smearing和nbnd多于电子数的一半。
因此想问一下siesta可不可以实现这个功能。,
以上英文是从手册上找到的,但不敢确认,想请教大家一下!
另外,可否发一下siesta计算多铁性材料的例子?
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

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

linghuntuzi

金虫 (著名写手)


identation(金币+2):谢谢参与
本人认为可以,不过也不确定
交流无止境
2楼2010-09-26 08:14:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

★ ★ ★ ★ ★
identation(金币+2):谢谢参与
aylayl08(金币+2):感谢提供信息 2010-10-21 14:58:30
zzy870720z(金币+2):谢谢转载 2010-10-23 14:24:09
引用回帖:
Originally posted by identation at 2010-09-25 16:53:42:
It has to be stressed that the macroscopic polarization can only be meaningfully calculated
using this approach for insulators. Therefore, the presence of an energy gap is necessary,
and no band  ...

这个是我从一个blog上看到的,不知道对你有没有用。网址:http://blog.chinaunix.net/u2/85704/showart_2022856.html
========================================
Regarding real part of dielectric function



Regarding real part of dielectric function
Systemlabel.EPSIMG  generated by SIESTA, contains the imaginary part of the dielectric function in the energy range specified in the input file of SIESTA. the program The programs contained in the directory (siesta-2.0/Util/Optical) calculate different optical properties of a solid out of the file.

There are two small programs:

input.f

and

optical.f


Both can be compile independently:

     f77( f95, f90,....) input.f -o inpu
     
     f77( f95, f90,....)  optical.f -o optical


* The first program creates the input file for optical (e2.dat), containing the
imaginary part of the dielectric function, from the file Systemlabel.EPSIMG

* The second reads a file called 'e2.dat' containing
the imaginary part of the dieletric function, and creates the
files :
'epsilon_real.out'  containing the real part of the dielectric function
'epsilon_img.out'  containing the imaginary part of the dielectric function
'refrac_index.out'   containing the refraction index
'absorp_index.out'  containing the extintion coefficient
'absorp_coef.out'    containing the absorption coefficient in cm**-1
'reflectance.out'     containing the reflectance
'conductivity.out'     containing the optical conductivity in (ohm*m)**-1

To do so, you just need to type:

    ./ input < Systemlabel.EPSIMG
  
    ./optical
Karim
=========================================
弘德明志博学笃行
3楼2010-10-21 14:55:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 identation 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 uFOAtb1CDngh 2026-09-02 9/450 2026-09-03 20:20 by jyNB2vrbdUpR
[考博] 找导师 +7 yuanjiabao 2026-08-29 8/400 2026-09-03 20:15 by yinnelly
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 YLHlRHNwYkce 2026-09-02 8/400 2026-09-03 18:50 by jyNB2vrbdUpR
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 YLHlRHNwYkce 2026-09-02 7/350 2026-09-03 18:50 by jyNB2vrbdUpR
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 YLHlRHNwYkce 2026-09-02 11/550 2026-09-03 18:42 by jyNB2vrbdUpR
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 rFEsUBKXRll0 2026-09-02 6/300 2026-09-03 16:49 by T0rGB46095mJ
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ero8OE6tv9cu 2026-09-02 8/400 2026-09-03 16:38 by T0rGB46095mJ
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ero8OE6tv9cu 2026-09-02 5/250 2026-09-03 16:38 by T0rGB46095mJ
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 ero8OE6tv9cu 2026-09-02 10/500 2026-09-03 16:27 by T0rGB46095mJ
[基金申请] 国社科又开始会评了,不知道这次命运如何 +10 雨打竹帘 2026-08-30 14/700 2026-09-03 13:17 by qsd10086
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 YLHlRHNwYkce 2026-09-02 8/400 2026-09-03 12:03 by T0rGB46095mJ
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +7 瞬息宇宙 2026-09-01 14/700 2026-09-02 17:33 by ma0526
[基金申请] 学科评审组评审是指会评吗? +5 瞬息宇宙 2026-08-31 5/250 2026-09-02 10:13 by 雨冰共舞
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
信息提示
请填处理意见