24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1069  |  回复: 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个字符以上)
信息提示
请填处理意见