24小时热门版块排行榜    

查看: 1564  |  回复: 3
【悬赏金币】回答本帖问题,作者清清清行将赠送您 20 个金币

清清清行

新虫 (初入文坛)

[求助] FAC程序包计算磁子截面 已有1人参与

在用FAC计算Ni,Co,Rh中高Z原子的磁子截面的时候,最后结果报错显示“too many free orbitals, enlarge norbmap2: 5000 >= 5000”,有同学知道是哪里有问题吗,如何修改。

发自小木虫Android客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

清清清行

新虫 (初入文坛)

还有一个问题是:计算直接电离磁截面时,根据输出磁截面求出的总截面与FAC输出的总截面之间误差很大,有遇到这种情况的同学吗?

发自小木虫Android客户端
2楼2024-10-24 22:43:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rlafite

至尊木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
# calculate the photoionization and
# radiative recombination cross sections
InitializeMPI(8)
SetAtom('Ni')

# specify the configurations for both recombining
# and recombined ions.
Closed('1s')
Config('2s2 2p5 3s2 3p6 3d8 4s2', group = 'n1')
# Config('2s1 2*1', group = 'n2')
Config('2s2 2p6 3s2 3p6 3d8 4s2', group = 'rn2')

ConfigEnergy(0)
OptimizeRadial(['rn2'])
ConfigEnergy(1)

# configuration interaction between n=1 and n=2
# complexes are included for the recombining ion.
Structure('Ni0p.lev.b', ['n1'])
Structure('Ni0p.lev.b', ['rn2'])
MemENTable('Ni0p.lev.b')
PrintTable('Ni0p.lev.b', 'Ni0p.lev', 1)

RRTable('Ni0p.rr.b', ['rn2'], ['n1'])
PrintTable('Ni0p.rr.b', 'Ni0p.rr', 1)
FinalizeMPI()
3楼2024-10-25 00:43:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rlafite

至尊木虫 (正式写手)

【答案】应助回帖

# calculate the photoionization and
# radiative recombination cross sections
InitializeMPI(8)
SetAtom('Co')

# specify the configurations for both recombining
# and recombined ions.
Closed('1s')
Config('2s1 2p6 3s2 3p6 3d7 4s2', group = 'n1')
# Config('2s1 2*1', group = 'n2')
Config('2s2 2p6 3s2 3p6 3d7 4s2', group = 'rn2')

ConfigEnergy(0)
OptimizeRadial(['rn2'])
ConfigEnergy(1)

# configuration interaction between n=1 and n=2
# complexes are included for the recombining ion.
Structure('Co0s.lev.b', ['n1'])
Structure('Co0s.lev.b', ['rn2'])
MemENTable('Co0s.lev.b')
PrintTable('Co0s.lev.b', 'Co0s.lev', 1)

RRTable('Co0s.rr.b', ['rn2'], ['n1'])
PrintTable('Co0s.rr.b', 'Co0s.rr', 1)
FinalizeMPI()
4楼2024-10-25 00:44:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 清清清行 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] E0414, 我的本子有没有希望? +7 布布和一二 2026-06-17 9/450 2026-06-21 19:29 by 布布和一二
[考研] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 19:18 by q1p9ixelzp
[教师之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 18:58 by q1p9ixelzp
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 18:00 by q1p9ixelzp
[考研] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 4/200 2026-06-21 17:45 by q1p9ixelzp
[论文投稿] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 17:40 by q1p9ixelzp
[论文投稿] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 17:37 by q1p9ixelzp
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 4/200 2026-06-21 17:25 by ky8v3skvft
[考博] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 5/250 2026-06-21 17:20 by ky8v3skvft
[博后之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 6/300 2026-06-21 17:05 by ky8v3skvft
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 5/250 2026-06-21 17:00 by ky8v3skvft
[教师之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +5 vlfdw50o0u 2026-06-20 7/350 2026-06-21 16:57 by ky8v3skvft
[论文投稿] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 vlfdw50o0u 2026-06-20 5/250 2026-06-21 16:45 by ky8v3skvft
[教师之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 7/350 2026-06-21 16:40 by ky8v3skvft
[论文投稿] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 9/450 2026-06-21 16:37 by ky8v3skvft
[考博] 【全奖博士/科研助理/博后招生】新加坡南洋理工大学机械与航空航天学院 +3 robohaha 2026-06-15 4/200 2026-06-21 00:36 by lanbing1
[硕博家园] 一篇论文同时出现在两个期刊,一模一样,这算不算学术不端,请各位老师斧正。 +11 mahl005 2026-06-14 12/600 2026-06-18 17:01 by 770904134
[基金申请] 希望面上有个好结果 +7 碧水00 2026-06-16 7/350 2026-06-18 12:18 by wuke100666
[论文投稿] 三区计算机方向期刊推荐 +5 1457340941 2026-06-15 5/250 2026-06-17 13:16 by 会议编辑
[论文投稿] sci论文二审求助 +5 潘倍倍 2026-06-15 5/250 2026-06-16 10:15 by xs74101122
信息提示
请填处理意见