24小时热门版块排行榜    

查看: 2097  |  回复: 7

fangyongxinxi

新虫 (正式写手)

[交流] 【求助】如何利用软件,数值求解“一维无限深方势井”问题。 已有2人参与

想数值求解一些可以精确求解的量子系统,
选方势井作为例子,
一般的解析求解方法:写出H,及薛定谔方程,解析推导求解,利用边界条件得到最终的phi= sin() 形式。平时看到的电子密度分布,其仅仅是最终phi的图像形式。
我现在的想法是:如何直接通过求解偏微分薛定谔方程,数值得到phi(x),(如何直接避开量子数 n)。
同样的问题出现中氢原子中:得到薛定谔方程后,能不能直接求解二阶偏微分方程来得到数值解?而不需要分解变量,利用特殊函数求解。
谢谢。
回复此楼

» 猜你喜欢

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

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

coolrainbow

木虫 (著名写手)

未来国家冻凉

当然能了,这早就研究烂了的东西了~~具体讲就是有限查分或有限元解微分方程

[ Last edited by coolrainbow on 2010-10-16 at 16:36 ]
技术博客:http://hi.baidu.com/coolrainbow/blog
2楼2010-10-16 16:14:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fangyongxinxi

新虫 (正式写手)

你说的是求解微分方程的方法。不是我问题想要的结果。
引用回帖:
Originally posted by coolrainbow at 2010-10-16 16:14:56:
当然能了,这早就研究烂了的东西了~~具体讲就是有限查分或有限元解微分方程

[ Last edited by coolrainbow on 2010-10-16 at 16:36 ]

3楼2010-10-17 14:52:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fangyongxinxi

新虫 (正式写手)

因为这个方程解比较复杂
4楼2010-10-17 14:53:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

coolrainbow

木虫 (著名写手)

未来国家冻凉

fangyongxinxi(金币+1): 2010-10-20 11:05:50
那我就不明白你到底想要什么

现在做量子力学的,即使算数值解也要对S方程做一些变换才能做,不可能完全在原始的3N元变量的S方程求解,因为需要的格点数目太多了
技术博客:http://hi.baidu.com/coolrainbow/blog
5楼2010-10-17 15:16:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫


★ ★
zhou2009(金币+2): 2010-10-18 06:16:03
fangyongxinxi(金币+1): 2010-10-20 11:05:46
Levine的量子化学里给出了Numerov数值求解一维问题的例子,甚至给出了在excel里操作的方法,可以参考。将方法推广到三维求解氢原子这样单粒子体系也是容易的。
如上所说,数值精确求解多粒子体系很困难。这种所谓basis-set free方法用在量子化学里一般是将多粒子关联体系转化到独立粒子模型求解,比如数值的KS-DFT已经有不少人研究并投入应用了,有兴趣可以看看Numerical methods for electronic structure calculations of materials一文。
6楼2010-10-17 22:02:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gzgzgz

铜虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
请教多位专家,类似有限元格点能动态划分格点大小,或者能自适应调整基组及步进改善计算分辨率的实现是否是有希望的,能否超过纯解析及辅助数值求解?
引用回帖:
Originally posted by sobereva at 2010-10-17 22:02:16:
Levine的量子化学里给出了Numerov数值求解一维问题的例子,甚至给出了在excel里操作的方法,可以参考。将方法推广到三维求解氢原子这样单粒子体系也是容易的。
如上所说,数值精确求解多粒子体系很困难。这种所谓 ...

[ Last edited by gzgzgz on 2010-11-12 at 11:44 ]
7楼2010-11-12 11:32:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

coolrainbow

木虫 (著名写手)

未来国家冻凉


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by gzgzgz at 2010-11-12 11:32:00:
请教多位专家,类似有限元格点能动态划分格点大小,或者能自适应调整基组及步进改善计算分辨率的实现是否是有希望的,能否超过纯解析及辅助数值求解?


[ Last edited by gzgzgz on 2010-11-12 at 11:44 ]

有限元解SCF的文献浩如烟海,仅仅2010年我就见过两篇了
技术博客:http://hi.baidu.com/coolrainbow/blog
8楼2010-11-12 17:34:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fangyongxinxi 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 2026国自然函评费到账 +20 羊腰板 2026-08-21 23/1150 2026-08-25 16:34 by zsna
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 +6 ziyangfang 2026-08-25 8/400 2026-08-25 16:28 by 宝贝的猪
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 18/900 2026-08-25 14:50 by Weed114114
[基金申请] 什么时候开奖? +12 CrisMessi 2026-08-18 13/650 2026-08-25 14:40 by ZJTJZ
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[基金申请] 今天基金会出结果吗?20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[基金申请] 范进中举一文的中心思想 +7 炎黄贵胄 2026-08-22 8/400 2026-08-25 08:48 by ZJTJZ
[基金申请] 只有每年这种时候来逛逛小木虫 +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见