d****z 发帖数: 9 | 1 根本吸引不了聪明的新人投入其中,本身也确实没有什么值得研究的新方向。
现有的中老年研究人员知识缺口很大,普遍统计力学,加声光电磁知识理论缺乏,
这样要解释现象,解决问题就是盲人骑瞎马。
其实上次量化得诺奖的Pople就没有什么理论贡献,数值计算方面也不如 Pulay有创造
力。
物理DFT方面倒还是有点小用。 |
F****I 发帖数: 270 | 2 Garnet Chan还算聪明吧,大概很年轻(~20岁)的时候就改进Lieb-Oxford bound了
您觉得高温超导机制算一个值得研究的方向吗?
有了量子计算机之后DFT的优势会被指数级的削弱吧
【在 d****z 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 根本吸引不了聪明的新人投入其中,本身也确实没有什么值得研究的新方向。 : 现有的中老年研究人员知识缺口很大,普遍统计力学,加声光电磁知识理论缺乏, : 这样要解释现象,解决问题就是盲人骑瞎马。 : 其实上次量化得诺奖的Pople就没有什么理论贡献,数值计算方面也不如 Pulay有创造 : 力。 : 物理DFT方面倒还是有点小用。
|
d****z 发帖数: 9 | 3 我觉得Mazziotti, Chan都不错,但因为方向局限,难有突破。
你觉得Aspuru-Guzik如何?
高温超导机制不能算量化研究反向吧?那是无穷大自由度加电子输运,量化方法根本无
从下手。你要研究你肯定是要用物理方法和物理凝聚态的人交流,何必要自称量化。
等量子计算机发明了,你觉得还会有人拿它来进行现在愚蠢的量化热化学计算码?
【在 F****I 的大作中提到】![](/moin_static193/solenoid/img/up.png) : Garnet Chan还算聪明吧,大概很年轻(~20岁)的时候就改进Lieb-Oxford bound了 : 您觉得高温超导机制算一个值得研究的方向吗? : 有了量子计算机之后DFT的优势会被指数级的削弱吧
|
F****I 发帖数: 270 | 4 Mazziotti和Chan的东西区别很大吧,前者要面对N-representable问题,而这个很可能
不能在多项式时间内解决。DMRG是目前处理静态电子相关和强关联系统最强大的工具,
而这个领域是量子化学最难处理的问题之一,个人态度很积极
不过Chan和Aspuru-Guzik的原始想法其实都来自物理系。我觉得固体同样是量子化学的对象
,用凝聚态物理排除量子化学完全是画地为牢
量化热化学计算? 这是什么,对不起没明白
【在 d****z 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 我觉得Mazziotti, Chan都不错,但因为方向局限,难有突破。 : 你觉得Aspuru-Guzik如何? : 高温超导机制不能算量化研究反向吧?那是无穷大自由度加电子输运,量化方法根本无 : 从下手。你要研究你肯定是要用物理方法和物理凝聚态的人交流,何必要自称量化。 : 等量子计算机发明了,你觉得还会有人拿它来进行现在愚蠢的量化热化学计算码?
|
d****z 发帖数: 9 | 5 我怎么觉得Aspuru-Guzik就是会呼悠呢?最近发文章很猛,养好多物理系来的薄厚,他
本人有本事码?
你要搞固体,你为什么非要把它叫“量子化学”呢。化学这个学科就是伴娘,研究的对
象空间,时间尺度都有限的,不然都叫物理好了。再加化学从炼丹术士而来,根本上不
是搞理论的,所以化学还是以有机化学为主,其他的都可叫可不叫化学,你又为何执着
于化学。
Thermochemistry啦,就是计算总能啦。
的对象
【在 F****I 的大作中提到】![](/moin_static193/solenoid/img/up.png) : Mazziotti和Chan的东西区别很大吧,前者要面对N-representable问题,而这个很可能 : 不能在多项式时间内解决。DMRG是目前处理静态电子相关和强关联系统最强大的工具, : 而这个领域是量子化学最难处理的问题之一,个人态度很积极 : 不过Chan和Aspuru-Guzik的原始想法其实都来自物理系。我觉得固体同样是量子化学的对象 : ,用凝聚态物理排除量子化学完全是画地为牢 : 量化热化学计算? 这是什么,对不起没明白
|
F****I 发帖数: 270 | 6 和Aspuru-Guzik本人没有接触过,不便评价。不过他做的方向都是热门,无论谁解决问
题都会有很大的影响。
炼丹本来炼的就是无机物,有机物有也烧成碳了,而且还是固体。为什么化学就要有有
机化学为主?再说一个学科的历史又不意味着现状和未来,这里的逻辑是什么?
我为什么不能把固体叫量子化学?叫了化学就不能同时叫物理吗? 这里有什么不相容原理
没理由不能用量子计算机计算总能,算得更准更大,肯定会有人去干这种事情 |
h********r 发帖数: 821 | 7 啥叫物理方法?如果有人出来说你们用的都是数学方法,咋反驳?
另外,开会聊的时候有人说量子计算机很扯,跟随者stupid。
个人认为所有攻击性词语都不好,对科学发展都没有益处。不喜欢的,不做就是了。当
然,出于利益关系(比如有funding竞争)的另当别论。
【在 d****z 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 我觉得Mazziotti, Chan都不错,但因为方向局限,难有突破。 : 你觉得Aspuru-Guzik如何? : 高温超导机制不能算量化研究反向吧?那是无穷大自由度加电子输运,量化方法根本无 : 从下手。你要研究你肯定是要用物理方法和物理凝聚态的人交流,何必要自称量化。 : 等量子计算机发明了,你觉得还会有人拿它来进行现在愚蠢的量化热化学计算码?
|
d****z 发帖数: 9 | 8 用量化方法算石墨烯总能会告诉你 "there are two-dimensional gas of massless
Dirac fermions in graphene"?
"用量子计算机计算总能,算得更准更大"绝不是一个聪明人该干的事。
其实前面的什么“N-representable"的问题是另一类聪明人不该干的,没有一个伟大
的数学家或应数学家觉得这是一个有趣的问题,只有一个加拿大的三流数学家在攻这个
问题,可笑的是他还曾经想让费曼对这个问题感兴趣。
容原理
【在 F****I 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 和Aspuru-Guzik本人没有接触过,不便评价。不过他做的方向都是热门,无论谁解决问 : 题都会有很大的影响。 : 炼丹本来炼的就是无机物,有机物有也烧成碳了,而且还是固体。为什么化学就要有有 : 机化学为主?再说一个学科的历史又不意味着现状和未来,这里的逻辑是什么? : 我为什么不能把固体叫量子化学?叫了化学就不能同时叫物理吗? 这里有什么不相容原理 : 没理由不能用量子计算机计算总能,算得更准更大,肯定会有人去干这种事情
|
h********r 发帖数: 821 | 9 基本是外行,讲错了请纠正。
David确实很聪明,但是目前的东西既有理论上的困难,也还没有到实际应用的程度,
还处于摸索阶段。DMRG适用的范围小了点。
的对象
【在 F****I 的大作中提到】![](/moin_static193/solenoid/img/up.png) : Mazziotti和Chan的东西区别很大吧,前者要面对N-representable问题,而这个很可能 : 不能在多项式时间内解决。DMRG是目前处理静态电子相关和强关联系统最强大的工具, : 而这个领域是量子化学最难处理的问题之一,个人态度很积极 : 不过Chan和Aspuru-Guzik的原始想法其实都来自物理系。我觉得固体同样是量子化学的对象 : ,用凝聚态物理排除量子化学完全是画地为牢 : 量化热化学计算? 这是什么,对不起没明白
|
d****z 发帖数: 9 | 10 这个太好反驳了,因为数学方法里没有物理。
【在 h********r 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 啥叫物理方法?如果有人出来说你们用的都是数学方法,咋反驳? : 另外,开会聊的时候有人说量子计算机很扯,跟随者stupid。 : 个人认为所有攻击性词语都不好,对科学发展都没有益处。不喜欢的,不做就是了。当 : 然,出于利益关系(比如有funding竞争)的另当别论。
|
|
|
F****I 发帖数: 270 | 11
>: 用量化方法算石墨烯总能会告诉你 "there are two-dimensional gas of
massless
so what?
你说“你觉得还会有人拿它来进行现在愚蠢的量化热化学计算码”
又没说这个人聪明不聪明
so what?
【在 d****z 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 这个太好反驳了,因为数学方法里没有物理。
|
h********r 发帖数: 821 | 12 还有不少搞化学的人认为石墨烯没意思呢。还是那句话,萝卜白菜各有所爱。聪明不聪
明,水平几流这类offensive的议题没有多大意思。
【在 d****z 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 用量化方法算石墨烯总能会告诉你 "there are two-dimensional gas of massless : Dirac fermions in graphene"? : "用量子计算机计算总能,算得更准更大"绝不是一个聪明人该干的事。 : 其实前面的什么“N-representable"的问题是另一类聪明人不该干的,没有一个伟大 : 的数学家或应数学家觉得这是一个有趣的问题,只有一个加拿大的三流数学家在攻这个 : 问题,可笑的是他还曾经想让费曼对这个问题感兴趣。 : : 容原理
|
d****z 发帖数: 9 | 13 最后劝一句,
你要是有这个时间,精力,脑力,真没别要局限于量子化学里的问题。
个人觉得这些问题都不是好问题。而且现在的里面从业人员烂的居多。
老的一时半会不会退。
爬出这个烂坑,外面广阔的天地在等着你。
【在 F****I 的大作中提到】![](/moin_static193/solenoid/img/up.png) : : >: 用量化方法算石墨烯总能会告诉你 "there are two-dimensional gas of : massless : so what? : 你说“你觉得还会有人拿它来进行现在愚蠢的量化热化学计算码” : 又没说这个人聪明不聪明 : so what?
|
F****I 发帖数: 270 | 14 谢谢建议。您能否推荐一个广阔的天地?除了实验有机化学以外 |
h********r 发帖数: 821 | 15 “物理”方法也没有化学。
【在 d****z 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 这个太好反驳了,因为数学方法里没有物理。
|
d****z 发帖数: 9 | 16 我觉得“不少搞化学的人认为石墨烯没意思”就是把石墨烯当成另一个碳六十,
而根本不了解“there are two-dimensional gas of massless Dirac fermions in
graphene and so on .....".
【在 h********r 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 还有不少搞化学的人认为石墨烯没意思呢。还是那句话,萝卜白菜各有所爱。聪明不聪 : 明,水平几流这类offensive的议题没有多大意思。
|
d****z 发帖数: 9 | 17 个人觉得在现在这个时代搞一点工程应用类比较好。
金融统计网络运筹策略里面都有很多可做有趣的问题。
【在 F****I 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 谢谢建议。您能否推荐一个广阔的天地?除了实验有机化学以外
|
F****I 发帖数: 270 | 18 谢谢您的建议,我觉得我不喜欢实用的东西,这点比较麻烦
【在 d****z 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 个人觉得在现在这个时代搞一点工程应用类比较好。 : 金融统计网络运筹策略里面都有很多可做有趣的问题。
|
S*M 发帖数: 10832 | 19 QM/MM研究biochemical反应也算量子化学吧
你们觉得这个方向怎么样?还有没有搞头 |
f**********e 发帖数: 1994 | 20 老烂坑
【在 S*M 的大作中提到】![](/moin_static193/solenoid/img/up.png) : QM/MM研究biochemical反应也算量子化学吧 : 你们觉得这个方向怎么样?还有没有搞头
|
|
|
s****s 发帖数: 66 | 21 转一个有关理论物理的八卦,我想也适用于量子化学吧:
------------------------
突然想起来个古老的八卦, 听说的。大概是01或者02年的时候。俞允强(学过电动力学
的应该知道他哈)还是李淼(中科院物理所搞超弦的)来着,招收理论物理研究生。面
试。
老师:你家兄弟几个?
考生:就我一个。
老师:你学理论物理,你父母怎么办?
考生:……
------------------------
现在这个急功近利的社会,很难孕育理论界的人才呀!
另外Pople还是有贡献的吧,至少从Pople basis set 来说。
【在 d****z 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 根本吸引不了聪明的新人投入其中,本身也确实没有什么值得研究的新方向。 : 现有的中老年研究人员知识缺口很大,普遍统计力学,加声光电磁知识理论缺乏, : 这样要解释现象,解决问题就是盲人骑瞎马。 : 其实上次量化得诺奖的Pople就没有什么理论贡献,数值计算方面也不如 Pulay有创造 : 力。 : 物理DFT方面倒还是有点小用。
|
m*****e 发帖数: 1506 | 22 Gaussian 98养活了多少人啊
【在 d****z 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 根本吸引不了聪明的新人投入其中,本身也确实没有什么值得研究的新方向。 : 现有的中老年研究人员知识缺口很大,普遍统计力学,加声光电磁知识理论缺乏, : 这样要解释现象,解决问题就是盲人骑瞎马。 : 其实上次量化得诺奖的Pople就没有什么理论贡献,数值计算方面也不如 Pulay有创造 : 力。 : 物理DFT方面倒还是有点小用。
|
S*M 发帖数: 10832 | 23 有没有搞这方面的同学反驳一下
我也好奇,化学家能不能真的算出点儿生物学家关心的东西
别最后就是自娱自乐
【在 f**********e 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 老烂坑
|
s****s 发帖数: 66 | 24 生物学家关心的东西很大,只能用QM/MM做了,没什么不和谐的东西吧?
我是没做过的,难道算出来都是错的?
【在 S*M 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 有没有搞这方面的同学反驳一下 : 我也好奇,化学家能不能真的算出点儿生物学家关心的东西 : 别最后就是自娱自乐
|
F****I 发帖数: 270 | 25 想起一个学长说要发起海洛因超弦基金,用贩毒来支持理论物理。一段时间之后说如果弦论被实验验证了,他就提头来见。一直没明白这个过程怎么实现
【在 s****s 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 转一个有关理论物理的八卦,我想也适用于量子化学吧: : ------------------------ : 突然想起来个古老的八卦, 听说的。大概是01或者02年的时候。俞允强(学过电动力学 : 的应该知道他哈)还是李淼(中科院物理所搞超弦的)来着,招收理论物理研究生。面 : 试。 : 老师:你家兄弟几个? : 考生:就我一个。 : 老师:你学理论物理,你父母怎么办? : 考生:…… : ------------------------
|
a**********e 发帖数: 157 | 26 量化的任务的是解决化学中的实际问题,比较IQ/晒书袋/搞fancy math 都是无意义
的。 |
F****I 发帖数: 270 | 27 这种工作还是挺多的,比如一些新出的方法还没有解析梯度,或者在某些程序包里没有
。把解析梯度写进去可以大大提高计算效率。真是一片广阔的天地....
【在 a**********e 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 量化的任务的是解决化学中的实际问题,比较IQ/晒书袋/搞fancy math 都是无意义 : 的。
|
c*****m 发帖数: 2191 | 28 对的,量化的落脚点最终是化学
你要去解决实际的化学问题,比如反应的机理,分子间的作用,分子的结构
所有的计算方法,软件只是工具,归根到底解决的是实际的化学问题
如果不去解决化学问题,那就真是数学和物理了
【在 a**********e 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 量化的任务的是解决化学中的实际问题,比较IQ/晒书袋/搞fancy math 都是无意义 : 的。
|
c*****m 发帖数: 2191 | 29 找一个热门的反应,发个JACS也没问题
如果是个普通的反应,就没镐头了
【在 S*M 的大作中提到】![](/moin_static193/solenoid/img/up.png) : QM/MM研究biochemical反应也算量子化学吧 : 你们觉得这个方向怎么样?还有没有搞头
|
t******t 发帖数: 3045 | 30 哈哈
【在 s****s 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 转一个有关理论物理的八卦,我想也适用于量子化学吧: : ------------------------ : 突然想起来个古老的八卦, 听说的。大概是01或者02年的时候。俞允强(学过电动力学 : 的应该知道他哈)还是李淼(中科院物理所搞超弦的)来着,招收理论物理研究生。面 : 试。 : 老师:你家兄弟几个? : 考生:就我一个。 : 老师:你学理论物理,你父母怎么办? : 考生:…… : ------------------------
|
|
|
h********r 发帖数: 821 | 31 好奇地问一句:解析能量梯度应该是提高精度,而不一定提高计算的效率吧?只要有总
能量程序就可以做数值梯度,虽然很不准,但是要算的量还比较少。
Correct me if I'm wrong.
【在 F****I 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 这种工作还是挺多的,比如一些新出的方法还没有解析梯度,或者在某些程序包里没有 : 。把解析梯度写进去可以大大提高计算效率。真是一片广阔的天地....
|
a**********e 发帖数: 157 | 32 精度会高一些,计算量也小一些吧(至少二阶梯度吧),因为有些common
intermediates不依赖于构型(如果做gradient的话)。
【在 h********r 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 好奇地问一句:解析能量梯度应该是提高精度,而不一定提高计算的效率吧?只要有总 : 能量程序就可以做数值梯度,虽然很不准,但是要算的量还比较少。 : Correct me if I'm wrong.
|
s****s 发帖数: 66 | 33 当你在对高对称性分子做数值解,由于做数值取样,某些单点可能能收敛到不是你想要
的电子态上。这时你会由衷的感谢那些研究一阶甚至二阶能量解析解的化学家。
【在 h********r 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 好奇地问一句:解析能量梯度应该是提高精度,而不一定提高计算的效率吧?只要有总 : 能量程序就可以做数值梯度,虽然很不准,但是要算的量还比较少。 : Correct me if I'm wrong.
|
h********r 发帖数: 821 | 34 这是另外一个问题。数值梯度做不对很容易理解,而且大部分分子没有高对称性。
【在 s****s 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 当你在对高对称性分子做数值解,由于做数值取样,某些单点可能能收敛到不是你想要 : 的电子态上。这时你会由衷的感谢那些研究一阶甚至二阶能量解析解的化学家。
|
h********r 发帖数: 821 | 35 没有深入钻研过,但是算二阶梯度本身的scaling比算总能量的要高吧,这样的话数值
梯度的计算量scaling和总能量的基本一样,所以感觉计算量会小一些。没搞懂你说的
common intermediates是个啥玩意儿,请指教。
【在 a**********e 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 精度会高一些,计算量也小一些吧(至少二阶梯度吧),因为有些common : intermediates不依赖于构型(如果做gradient的话)。
|
C**R 发帖数: 1047 | |
C**R 发帖数: 1047 | |
C**R 发帖数: 1047 | 38 我学物理的,不懂你们在说什么。我们做物理的最关心的是如何把exchang和
correlation搞对。比如研究quantum monte carlo啊,研究GW啊,研究bethe-salpeter
equation啊,研究hybrid functional啊。剩下算的东西就是software engineering级
别的工作。
【在 s****s 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 当你在对高对称性分子做数值解,由于做数值取样,某些单点可能能收敛到不是你想要 : 的电子态上。这时你会由衷的感谢那些研究一阶甚至二阶能量解析解的化学家。
|
F****I 发帖数: 270 | 39 我试过的方法里面解析梯度都比数值梯度要快,越大的分子越明显。对这方面不是很熟悉,我觉得大概是这样的。比如最简单的HartreeFock,用Hellmann-Feynman定理得到解析梯度的表达,实际是收敛的分子轨道再乘以一些新的积分,相对数值梯度不需要重新做SCF,所以比较快。非变分方法的解析梯度会更复杂一些,比如CPHF技术,但也是用单点能计算的结果再此之上做一些计算。数值梯度至多可以用上一个计算收敛的结果作为initial guess, (SCF, coupled-cluster),这样会慢很多。不过更复杂的方法也可能解析梯度太复杂了,molcas里的CASPT2好像现在都没有解析梯度...
【在 h********r 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 好奇地问一句:解析能量梯度应该是提高精度,而不一定提高计算的效率吧?只要有总 : 能量程序就可以做数值梯度,虽然很不准,但是要算的量还比较少。 : Correct me if I'm wrong.
|
F****I 发帖数: 270 | 40 你们对那些model Hamiltonian怎么看呢?比如Anderson model
salpeter
【在 C**R 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 我学物理的,不懂你们在说什么。我们做物理的最关心的是如何把exchang和 : correlation搞对。比如研究quantum monte carlo啊,研究GW啊,研究bethe-salpeter : equation啊,研究hybrid functional啊。剩下算的东西就是software engineering级 : 别的工作。
|
|
|
c******n 发帖数: 16403 | 41 还是做量子计算的人聪明, 讨论的帖子有真才实料, 搞有机的多半都是灌水唧唧歪歪
一些无聊的八卦。
【在 F****I 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 你们对那些model Hamiltonian怎么看呢?比如Anderson model : : salpeter
|
F****I 发帖数: 270 | 42 我一直很佩服做有机实验的人啊。自己做的时候产率经常为零。开学的时候交了玻璃仪
器的押金,学期末都成了负数
【在 c******n 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 还是做量子计算的人聪明, 讨论的帖子有真才实料, 搞有机的多半都是灌水唧唧歪歪 : 一些无聊的八卦。
|
C**R 发帖数: 1047 | 43 好东西啊。真正学理论物理的人,不大看得起ab inito的东东。所谓的物理直觉好,就
是可以直接猜到hamiltonian长什么样子。当然缺点就是model都有adjustable
parameters的。好的物理学家就是尽量减少fitting parameters,混饭吃的就是小修小
补,加一项小修正,发个paper之类的。
不过,想解释高温超导,再等个50年吧。
【在 F****I 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 你们对那些model Hamiltonian怎么看呢?比如Anderson model : : salpeter
|
C**R 发帖数: 1047 | 44 唉,怪不得,量子化学从30年代就和做物理的分家了。原来还停留在薛定谔第一篇文章
算氢原子的水平。
熟悉,我觉得大概是这样的。比如最简单的HartreeFock,用Hellmann-Feynman定理得
到解析梯度的表达,实际是收敛的分子轨道再乘以一些新的积分,相对数值梯度不需要
重新做SCF,所以比较快。非变分方法的解析梯度会更复杂一些,比如CPHF技术,但也
是用单点能计算的结果再此之上做一些计算。数值梯度至多可以用上一个计算收敛的结
果作为initial guess, (SCF, coupled-cluster),这样会慢很多。不过更复杂的方法
也可能解析梯度太复杂了,molcas里的CASPT2好像现在: 都没有解析梯度...
【在 F****I 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 我试过的方法里面解析梯度都比数值梯度要快,越大的分子越明显。对这方面不是很熟悉,我觉得大概是这样的。比如最简单的HartreeFock,用Hellmann-Feynman定理得到解析梯度的表达,实际是收敛的分子轨道再乘以一些新的积分,相对数值梯度不需要重新做SCF,所以比较快。非变分方法的解析梯度会更复杂一些,比如CPHF技术,但也是用单点能计算的结果再此之上做一些计算。数值梯度至多可以用上一个计算收敛的结果作为initial guess, (SCF, coupled-cluster),这样会慢很多。不过更复杂的方法也可能解析梯度太复杂了,molcas里的CASPT2好像现在都没有解析梯度...
|
F****I 发帖数: 270 | 45 thx.
【在 C**R 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 好东西啊。真正学理论物理的人,不大看得起ab inito的东东。所谓的物理直觉好,就 : 是可以直接猜到hamiltonian长什么样子。当然缺点就是model都有adjustable : parameters的。好的物理学家就是尽量减少fitting parameters,混饭吃的就是小修小 : 补,加一项小修正,发个paper之类的。 : 不过,想解释高温超导,再等个50年吧。
|
F****I 发帖数: 270 | 46 啥样子叫超越了氢原子的水平呢?
【在 C**R 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 唉,怪不得,量子化学从30年代就和做物理的分家了。原来还停留在薛定谔第一篇文章 : 算氢原子的水平。 : : 熟悉,我觉得大概是这样的。比如最简单的HartreeFock,用Hellmann-Feynman定理得 : 到解析梯度的表达,实际是收敛的分子轨道再乘以一些新的积分,相对数值梯度不需要 : 重新做SCF,所以比较快。非变分方法的解析梯度会更复杂一些,比如CPHF技术,但也 : 是用单点能计算的结果再此之上做一些计算。数值梯度至多可以用上一个计算收敛的结 : 果作为initial guess, (SCF, coupled-cluster),这样会慢很多。不过更复杂的方法 : 也可能解析梯度太复杂了,molcas里的CASPT2好像现在: 都没有解析梯度...
|
C**R 发帖数: 1047 | 47 我真是不懂你们在做什么,如果说每个学校都有一个做量子化学的,你们真是命好啊。
population
【在 F****I 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 啥样子叫超越了氢原子的水平呢?
|
C**R 发帖数: 1047 | 48 当然是多体啦,quasiparticle啦。多是reciprocal space thinking啦。不过呢,因为
强关联多是d orbital的local bonding effect,近些年在物理界real space thinking
很有回潮。
不过呢,做物理的真的很不喜欢molecules。你去问物理的人,他们基本都不知道LCAO
是什么。但是没有不懂tight binding的,虽然都是一个东西。
【在 F****I 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 啥样子叫超越了氢原子的水平呢?
|
C**R 发帖数: 1047 | 49 btw,我说的多体可不是HF那么简单的近似。
thinking
LCAO
【在 C**R 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 当然是多体啦,quasiparticle啦。多是reciprocal space thinking啦。不过呢,因为 : 强关联多是d orbital的local bonding effect,近些年在物理界real space thinking : 很有回潮。 : 不过呢,做物理的真的很不喜欢molecules。你去问物理的人,他们基本都不知道LCAO : 是什么。但是没有不懂tight binding的,虽然都是一个东西。
|
F****I 发帖数: 270 | 50 那为啥K Wilson说“Ab initio quantum chemistry is an emerging computational area that is fifty years ahead of lattice gauge theory"呢
Anderson的Resonance Valence Bond也是从化学系类比来呢?
thinking
LCAO
【在 C**R 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 当然是多体啦,quasiparticle啦。多是reciprocal space thinking啦。不过呢,因为 : 强关联多是d orbital的local bonding effect,近些年在物理界real space thinking : 很有回潮。 : 不过呢,做物理的真的很不喜欢molecules。你去问物理的人,他们基本都不知道LCAO : 是什么。但是没有不懂tight binding的,虽然都是一个东西。
|
|
|
C**R 发帖数: 1047 | 51 wilson也就懂group renormalization那些统计物理的东西,不不明白他凭什么把ab
initio的东西和lattice gauge theory这种抽象的model system比。也不知道他这个50
年是怎么算的。
学物理的从化学那边借点东西也很正常啊。
我的一个committee member是做dft的(物理的,mit出身)。他说他们这些做dft的人
,最大的一个梦想就是将来可以精确的predict化学反应速率。想要把化学反应速率算
对,就要把各个反应阶段的free energy算对。这个就必须要把exchange-correlation
搞对。这个是最根本的东西。我还发现有个做化学的就是用pbe算了算反应函,和实验
比了比,发现比较接近,就得了nobel奖。真是shame啊。你想想,就是water kohn也没
有得物理奖,对不?
当然我不是贬低做量化的人,任何一个学科存在都是有其必要性的。而且现在强关联体
系的东东似乎用量化的一些方法做可以得到更transparent的理解。物理的做dft的大方
向是要把能量算对。但是现在也没有什么理论进展,所以算的快也是绝大多数做dft人
主攻的方向。这样就和你们做量化的殊途同归了。
area that is fifty years ahead of lattice gauge theory"呢
【在 F****I 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 那为啥K Wilson说“Ab initio quantum chemistry is an emerging computational area that is fifty years ahead of lattice gauge theory"呢 : Anderson的Resonance Valence Bond也是从化学系类比来呢? : : thinking : LCAO
|
a**********e 发帖数: 157 | 52 譬如用coupled cluster算nuclear gradient,只需要把一个Lambda vector(与构型变
化无关,只在那一个计算的构型求解)求出来,那么算不同dimension(x1, y1, x2, .
..)上的gradient,不同的只有integral gradient(和轨道系数变化?),其他的量不
变。
【在 h********r 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 没有深入钻研过,但是算二阶梯度本身的scaling比算总能量的要高吧,这样的话数值 : 梯度的计算量scaling和总能量的基本一样,所以感觉计算量会小一些。没搞懂你说的 : common intermediates是个啥玩意儿,请指教。
|
a**********e 发帖数: 157 | 53 这个算算反应han就得物理奖的是谁啊?
50
correlation
【在 C**R 的大作中提到】![](/moin_static193/solenoid/img/up.png) : wilson也就懂group renormalization那些统计物理的东西,不不明白他凭什么把ab : initio的东西和lattice gauge theory这种抽象的model system比。也不知道他这个50 : 年是怎么算的。 : 学物理的从化学那边借点东西也很正常啊。 : 我的一个committee member是做dft的(物理的,mit出身)。他说他们这些做dft的人 : ,最大的一个梦想就是将来可以精确的predict化学反应速率。想要把化学反应速率算 : 对,就要把各个反应阶段的free energy算对。这个就必须要把exchange-correlation : 搞对。这个是最根本的东西。我还发现有个做化学的就是用pbe算了算反应函,和实验 : 比了比,发现比较接近,就得了nobel奖。真是shame啊。你想想,就是water kohn也没 : 有得物理奖,对不?
|
a****o 发帖数: 6612 | 54 你对量子化学很有误解。 Kohn 和 Pople 一起分享了1998年的诺贝尔化学奖。
PBE算反应焓是决计得不了诺贝尔奖的。得诺贝尔奖的倒是有可能用PBE算算东西。
50
correlation
【在 C**R 的大作中提到】![](/moin_static193/solenoid/img/up.png) : wilson也就懂group renormalization那些统计物理的东西,不不明白他凭什么把ab : initio的东西和lattice gauge theory这种抽象的model system比。也不知道他这个50 : 年是怎么算的。 : 学物理的从化学那边借点东西也很正常啊。 : 我的一个committee member是做dft的(物理的,mit出身)。他说他们这些做dft的人 : ,最大的一个梦想就是将来可以精确的predict化学反应速率。想要把化学反应速率算 : 对,就要把各个反应阶段的free energy算对。这个就必须要把exchange-correlation : 搞对。这个是最根本的东西。我还发现有个做化学的就是用pbe算了算反应函,和实验 : 比了比,发现比较接近,就得了nobel奖。真是shame啊。你想想,就是water kohn也没 : 有得物理奖,对不?
|
a****o 发帖数: 6612 | 55 没有关系。好多搞工程的,还停留在牛顿力学的水平呢。
【在 C**R 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 唉,怪不得,量子化学从30年代就和做物理的分家了。原来还停留在薛定谔第一篇文章 : 算氢原子的水平。 : : 熟悉,我觉得大概是这样的。比如最简单的HartreeFock,用Hellmann-Feynman定理得 : 到解析梯度的表达,实际是收敛的分子轨道再乘以一些新的积分,相对数值梯度不需要 : 重新做SCF,所以比较快。非变分方法的解析梯度会更复杂一些,比如CPHF技术,但也 : 是用单点能计算的结果再此之上做一些计算。数值梯度至多可以用上一个计算收敛的结 : 果作为initial guess, (SCF, coupled-cluster),这样会慢很多。不过更复杂的方法 : 也可能解析梯度太复杂了,molcas里的CASPT2好像现在: 都没有解析梯度...
|
a**********e 发帖数: 157 | 56 关键是能够解决重要的问题,make contribution。譬如分子激发态的问题,很复杂,
又化学上很重要,从头算提供了各种方法,尽管有这样那样的缺点;批评容易,提出更
好的建设性的方法才是困难的。 |
s****s 发帖数: 66 | 57 HF挺好的嘛哈哈,物理的越做越小,化学家想弄大点的分子。你看玩更大分子的的,HF
的近似都算不动了,只好用QM/MM了。
【在 C**R 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 唉,怪不得,量子化学从30年代就和做物理的分家了。原来还停留在薛定谔第一篇文章 : 算氢原子的水平。 : : 熟悉,我觉得大概是这样的。比如最简单的HartreeFock,用Hellmann-Feynman定理得 : 到解析梯度的表达,实际是收敛的分子轨道再乘以一些新的积分,相对数值梯度不需要 : 重新做SCF,所以比较快。非变分方法的解析梯度会更复杂一些,比如CPHF技术,但也 : 是用单点能计算的结果再此之上做一些计算。数值梯度至多可以用上一个计算收敛的结 : 果作为initial guess, (SCF, coupled-cluster),这样会慢很多。不过更复杂的方法 : 也可能解析梯度太复杂了,molcas里的CASPT2好像现在: 都没有解析梯度...
|
h*****y 发帖数: 298 | 58 DFT已经是凝聚态物理里研究材料的标准工具了。model hamitonian最后不回到数值解
的有几个?现在计算机的计算能力越来越强,不管凝聚态还是化学,ab initio是个趋
势。强关联问题现在DFT还不能完全处理,这反过来也说明这个领域大有可为。DMFT, G
-W, S-B这些都是可以继续完善的方案。得到精确ec也不是完全没有办法,LDA方案中的
ec参数就可以用QMC算出来。回头说说传统理论物理,现在的趋势是数学越来越fancy,
但是在问题的实质上并没有新的突破。一个例子就是 ferroelectricity,多少年了都是
Ginzburg–Landau theory,DFT一下子就把认识推进到了微观层次。不管物理还是化学
,理论还是计算,只要有新见解,就是好的研究,只是打着纯理论的幌子,即使方法再
高深,数学推导再繁琐,也是不可取的。
50
correlation
【在 C**R 的大作中提到】![](/moin_static193/solenoid/img/up.png) : wilson也就懂group renormalization那些统计物理的东西,不不明白他凭什么把ab : initio的东西和lattice gauge theory这种抽象的model system比。也不知道他这个50 : 年是怎么算的。 : 学物理的从化学那边借点东西也很正常啊。 : 我的一个committee member是做dft的(物理的,mit出身)。他说他们这些做dft的人 : ,最大的一个梦想就是将来可以精确的predict化学反应速率。想要把化学反应速率算 : 对,就要把各个反应阶段的free energy算对。这个就必须要把exchange-correlation : 搞对。这个是最根本的东西。我还发现有个做化学的就是用pbe算了算反应函,和实验 : 比了比,发现比较接近,就得了nobel奖。真是shame啊。你想想,就是water kohn也没 : 有得物理奖,对不?
|
C**R 发帖数: 1047 | 59 “得到精确ec也不是完全没有办法,LDA方案中的
ec参数就可以用QMC算出来。”
一看这句就知道你是半壶水。
G
【在 h*****y 的大作中提到】![](/moin_static193/solenoid/img/up.png) : DFT已经是凝聚态物理里研究材料的标准工具了。model hamitonian最后不回到数值解 : 的有几个?现在计算机的计算能力越来越强,不管凝聚态还是化学,ab initio是个趋 : 势。强关联问题现在DFT还不能完全处理,这反过来也说明这个领域大有可为。DMFT, G : -W, S-B这些都是可以继续完善的方案。得到精确ec也不是完全没有办法,LDA方案中的 : ec参数就可以用QMC算出来。回头说说传统理论物理,现在的趋势是数学越来越fancy, : 但是在问题的实质上并没有新的突破。一个例子就是 ferroelectricity,多少年了都是 : Ginzburg–Landau theory,DFT一下子就把认识推进到了微观层次。不管物理还是化学 : ,理论还是计算,只要有新见解,就是好的研究,只是打着纯理论的幌子,即使方法再 : 高深,数学推导再繁琐,也是不可取的。 :
|
h********r 发帖数: 821 | 60 LCAO和tight binding根本不是一个层面上的东西,理论化学中和tight binding等价的
是Hukel theory。
thinking
LCAO
【在 C**R 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 当然是多体啦,quasiparticle啦。多是reciprocal space thinking啦。不过呢,因为 : 强关联多是d orbital的local bonding effect,近些年在物理界real space thinking : 很有回潮。 : 不过呢,做物理的真的很不喜欢molecules。你去问物理的人,他们基本都不知道LCAO : 是什么。但是没有不懂tight binding的,虽然都是一个东西。
|
|
|
h********r 发帖数: 821 | 61 一阶的解析梯度我相信比较快,但是二阶的应该scaling比总能量高(比如算振动频率
比算总能量慢)。解析能量梯度推导太麻烦了,所以从来没有认真钻研过Pulay的那一
套,认识比较模糊。
熟悉,我觉得大概是这样的。比如最简单的HartreeFock,用Hellmann-Feynman定理得
到解析梯度的表达,实际是收敛的分子轨道再乘以一些新的积分,相对数值梯度不需要
重新做SCF,所以比较快。非变分方法的解析梯度会更复杂一些,比如CPHF技术,但也
是用单点能计算的结果再此之上做一些计算。数值梯度至多可以用上一个计算收敛的结
果作为initial guess, (SCF, coupled-cluster),这样会慢很多。不过更复杂的方法
也可能解析梯度太复杂了,molcas里的CASPT2好像现在: 都没有解析梯度...
【在 F****I 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 我试过的方法里面解析梯度都比数值梯度要快,越大的分子越明显。对这方面不是很熟悉,我觉得大概是这样的。比如最简单的HartreeFock,用Hellmann-Feynman定理得到解析梯度的表达,实际是收敛的分子轨道再乘以一些新的积分,相对数值梯度不需要重新做SCF,所以比较快。非变分方法的解析梯度会更复杂一些,比如CPHF技术,但也是用单点能计算的结果再此之上做一些计算。数值梯度至多可以用上一个计算收敛的结果作为initial guess, (SCF, coupled-cluster),这样会慢很多。不过更复杂的方法也可能解析梯度太复杂了,molcas里的CASPT2好像现在都没有解析梯度...
|
F****I 发帖数: 270 | 62 所以他们是engineer...
【在 a****o 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 没有关系。好多搞工程的,还停留在牛顿力学的水平呢。
|
F****I 发帖数: 270 | 63 如果是单个体系固定基组的计算,快慢也可能是prefactor导致的。我觉得analytical frequency,gradient和单点能的scaling是一样的,因为没有额外的积分变换和更高阶的excitation引入,当然我的感觉可能是错的
【在 h********r 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 一阶的解析梯度我相信比较快,但是二阶的应该scaling比总能量高(比如算振动频率 : 比算总能量慢)。解析能量梯度推导太麻烦了,所以从来没有认真钻研过Pulay的那一 : 套,认识比较模糊。 : : 熟悉,我觉得大概是这样的。比如最简单的HartreeFock,用Hellmann-Feynman定理得 : 到解析梯度的表达,实际是收敛的分子轨道再乘以一些新的积分,相对数值梯度不需要 : 重新做SCF,所以比较快。非变分方法的解析梯度会更复杂一些,比如CPHF技术,但也 : 是用单点能计算的结果再此之上做一些计算。数值梯度至多可以用上一个计算收敛的结 : 果作为initial guess, (SCF, coupled-cluster),这样会慢很多。不过更复杂的方法 : 也可能解析梯度太复杂了,molcas里的CASPT2好像现在: 都没有解析梯度...
|
h*****y 发帖数: 298 | 64 如果我哪里说错了,你完全可以解释清楚。
【在 C**R 的大作中提到】![](/moin_static193/solenoid/img/up.png) : “得到精确ec也不是完全没有办法,LDA方案中的 : ec参数就可以用QMC算出来。” : 一看这句就知道你是半壶水。 : : G
|
a****o 发帖数: 6612 | 65 单点能的计算使用固定的分子几何构象。分子中原子的坐标是固定的。计算成本基本完
全是为考虑电子间相互作用。但是能量梯度不同,与分子中原子的数目以及对称性有关
。即使不考虑原子数目(或者说独立内坐标的数目),能量解析梯度也要远比单点能费
时间。根据不同理论,多费的时间多少,差别甚大。
如果你算过构型优化和频率,你会发现频率分析中的Hessian(能量的二级梯度)非常
耗费时间。 二级能量梯度基本上都是数值梯度,原因是解析梯度形式复杂,过于耗时。
计算能量梯度需要的积分,可以由原来的原子轨道积分获得,但是角动量要加一。 计
算能量梯队需要比单点能计算更多的积分,因而计算量也随之增大。(相比单点能的积
分,数目上要多出最高角动量的积分。比如你用DZ的基组,最高是d轨道。那么能量梯
度就需要类似f轨道的积分。要知道原子d轨道是五重简并,到f轨道就是七重了。依此
类推。)同理,如果是二级解析梯度,角动量要加二。DZ基组的主族元素,就需要类似
f轨道的积分,从基组的角度上,相当于QZ的计算了
analytical frequency,gradient和单点能的scaling是一样的,因为没有额外的积分变
换和更高阶的excitation引入,当然我的感觉可能是错的
【在 F****I 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 如果是单个体系固定基组的计算,快慢也可能是prefactor导致的。我觉得analytical frequency,gradient和单点能的scaling是一样的,因为没有额外的积分变换和更高阶的excitation引入,当然我的感觉可能是错的
|
F****I 发帖数: 270 | 66 我不否定解析梯度和频率比单个单点耗时,您针对帖子是这种增加来自于prefactor还是scaling。那么,对于HF的解析梯度和hessian,它的scaling是否超过了N^4?
>二级能量梯度基本上都是数值梯度,原因是解析梯度形式复杂,过于耗时。
在molpro上测试了一下HF/cc-pVDZ对于丙烷的解析和数值Hessian
解析的
PROGRAMS * TOTAL FREQ OPTG HF INT
CPU TIMES * 18.29 14.39 3.62 0.02 0.17
REAL TIME * 19.56 SEC
DISK USED * 64.30 MB
SF USED * 0.65 MB
GA USED * 0.00 MB (max) 0.00 MB (current)
**********************************************************************************************************************************
数值的
PROGRAMS * TOTAL FREQ FREQ OPTG HF INT
CPU TIMES * 131.42 113.12 14.39 3.62 0.02 0.17
REAL TIME * 209.82 SEC
DISK USED * 155.30 MB
SF USED * 0.65 MB
GA USED * 0.00 MB (max) 0.00 MB (current)
**********************************************************************************************************************************
这里没有CFOUR/MRCC支持coupled-cluster的解析Hessian,以后有机会再比较一下吧
时。
【在 a****o 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 单点能的计算使用固定的分子几何构象。分子中原子的坐标是固定的。计算成本基本完 : 全是为考虑电子间相互作用。但是能量梯度不同,与分子中原子的数目以及对称性有关 : 。即使不考虑原子数目(或者说独立内坐标的数目),能量解析梯度也要远比单点能费 : 时间。根据不同理论,多费的时间多少,差别甚大。 : 如果你算过构型优化和频率,你会发现频率分析中的Hessian(能量的二级梯度)非常 : 耗费时间。 二级能量梯度基本上都是数值梯度,原因是解析梯度形式复杂,过于耗时。 : 计算能量梯度需要的积分,可以由原来的原子轨道积分获得,但是角动量要加一。 计 : 算能量梯队需要比单点能计算更多的积分,因而计算量也随之增大。(相比单点能的积 : 分,数目上要多出最高角动量的积分。比如你用DZ的基组,最高是d轨道。那么能量梯 : 度就需要类似f轨道的积分。要知道原子d轨道是五重简并,到f轨道就是七重了。依此
|
C**R 发帖数: 1047 | 67 你丫就是想说tight binding是a very simple linear combination of atomic
orbitals molecular orbitals?
The Hückel method or Hückel molecular orbital method (HMO) proposed by
Erich Hückel in 1930, is a very simple linear combination of atomic
orbitals molecular orbitals (LCAO MO) method for the determination of
energies of molecular orbitals of pi electrons in conjugated hydrocarbon
systems, such as ethene, benzene and butadiene.[1][2] It is the theoretical
basis for the Hückel's rule; the extended Hückel method developed by Roald
Hoffmann is the basis of the Woodward-Hoffmann rules.[3] It was later
extended to conjugated molecules such as pyridine, pyrrole and furan that
contain atoms other than carbon, known in this context as heteroatoms.[4]
【在 h********r 的大作中提到】![](/moin_static193/solenoid/img/up.png) : LCAO和tight binding根本不是一个层面上的东西,理论化学中和tight binding等价的 : 是Hukel theory。 : : thinking : LCAO
|
a****o 发帖数: 6612 | 68 数值解中你利用对称性了吗?个别例子可能有差别。
还是scaling。那么,对于HF的解析梯度和hessian,它的scaling是否超过了N^4?
【在 F****I 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 我不否定解析梯度和频率比单个单点耗时,您针对帖子是这种增加来自于prefactor还是scaling。那么,对于HF的解析梯度和hessian,它的scaling是否超过了N^4? : >二级能量梯度基本上都是数值梯度,原因是解析梯度形式复杂,过于耗时。 : 在molpro上测试了一下HF/cc-pVDZ对于丙烷的解析和数值Hessian : 解析的 : PROGRAMS * TOTAL FREQ OPTG HF INT : CPU TIMES * 18.29 14.39 3.62 0.02 0.17 : REAL TIME * 19.56 SEC : DISK USED * 64.30 MB : SF USED * 0.65 MB : GA USED * 0.00 MB (max) 0.00 MB (current)
|
C**R 发帖数: 1047 | 69 首先
Accurate quantum Monte Carlo (QMC) simulations for the energy of the HEG
have been performed for several intermediate values of the density, in turn
providing accurate values of the correlation energy density.[3]
Vosko-Wilk-Nusair (VWN) [4]
Perdew-Zunger (PZ81) [5]
Cole-Perdew (CP) [6]
Perdew-Wang (PW92) [7]
HEG就是homogeneous electron gas,现在绝大多数functional都是基于这个
approximation的,问题就是这个HEG approximation很难进一步有突破(因为qmc实在
是太慢了),所以你说的那个什么精确Eex,就是你不懂。
在问你个我知道答案问题,为什么QMC算Cuprates高温超导算的很好。但是一碰
到铁基超导就歇菜了呢?
发信人: CCMR (花和尚), 信区: Chemistry
标 题: Re: 量子化学已经接近完全死亡了
发信站: BBS 未名空间站 (Thu Jul 28 23:10:03 2011, 美东)
“得到精确ec也不是完全没有办法,LDA方案中的
ec参数就可以用QMC算出来。”
一看这句就知道你是半壶水。
【在 h*****y 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 如果我哪里说错了,你完全可以解释清楚。
|
h*****y 发帖数: 298 | 70 HEG的基态能是解析的,QMC做了精确的参数化,这就是为什么大家都用HEG作为半经验
ec项。带来的问题就是引入了平均场近似。但是这不是QMC的问题,而是DFT计算的问题
。QMC能得到精确ec这有什么问题? 关于铁基超导的问题你可以描述一下到底什么叫歇
菜?然后我才可以考虑。
turn
【在 C**R 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 首先 : Accurate quantum Monte Carlo (QMC) simulations for the energy of the HEG : have been performed for several intermediate values of the density, in turn : providing accurate values of the correlation energy density.[3] : Vosko-Wilk-Nusair (VWN) [4] : Perdew-Zunger (PZ81) [5] : Cole-Perdew (CP) [6] : Perdew-Wang (PW92) [7] : HEG就是homogeneous electron gas,现在绝大多数functional都是基于这个 : approximation的,问题就是这个HEG approximation很难进一步有突破(因为qmc实在
|
|
|
F****I 发帖数: 270 | 71 用了。再给一个没有空间对称性的,乙醇,HF/cc-pVDZ opt + freq, Gaussian09
解析频率
Job cpu time: 0 days 0 hours 0 minutes 17.0 seconds.
数值频率
Job cpu time: 0 days 0 hours 2 minutes 56.0 seconds.
【在 a****o 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 数值解中你利用对称性了吗?个别例子可能有差别。 : : 还是scaling。那么,对于HF的解析梯度和hessian,它的scaling是否超过了N^4?
|
a****o 发帖数: 6612 | 72 有意思。我没有比较过。不过能有解析解,就不会用近似的数值解了。其实我是说,用
了对称性,数值解会更快一些。 但是公平的比较是都不用对称性。
很好奇大基组会怎样,还有MP2, CC的又会怎样。HF的单点能下面能量是最低的,轨道
也是最优的。这两个条件可能会是解析Hessian的计算比较便利一些。如果到相关能计
算的方法,情形可能完全不同了。
另外数值解是近似方法,其计算时间与取多少点(三点,五点)也有关系。
【在 F****I 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 用了。再给一个没有空间对称性的,乙醇,HF/cc-pVDZ opt + freq, Gaussian09 : 解析频率 : Job cpu time: 0 days 0 hours 0 minutes 17.0 seconds. : 数值频率 : Job cpu time: 0 days 0 hours 2 minutes 56.0 seconds.
|
F****I 发帖数: 270 | 73 好吧好吧
乙醇,MP2/cc-pVDZ, Gaussian09
解析频率
Job cpu time: 0 days 0 hours 1 minutes 56.0 seconds.
数值频率
Job cpu time: 0 days 0 hours 8 minutes 25.0 seconds.
更大基组自行测试吧
【在 a****o 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 有意思。我没有比较过。不过能有解析解,就不会用近似的数值解了。其实我是说,用 : 了对称性,数值解会更快一些。 但是公平的比较是都不用对称性。 : 很好奇大基组会怎样,还有MP2, CC的又会怎样。HF的单点能下面能量是最低的,轨道 : 也是最优的。这两个条件可能会是解析Hessian的计算比较便利一些。如果到相关能计 : 算的方法,情形可能完全不同了。 : 另外数值解是近似方法,其计算时间与取多少点(三点,五点)也有关系。
|
h********r 发帖数: 821 | 74 讨论问题没必要说脏话。
theoretical
Roald
【在 C**R 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 你丫就是想说tight binding是a very simple linear combination of atomic : orbitals molecular orbitals? : The Hückel method or Hückel molecular orbital method (HMO) proposed by : Erich Hückel in 1930, is a very simple linear combination of atomic : orbitals molecular orbitals (LCAO MO) method for the determination of : energies of molecular orbitals of pi electrons in conjugated hydrocarbon : systems, such as ethene, benzene and butadiene.[1][2] It is the theoretical : basis for the Hückel's rule; the extended Hückel method developed by Roald : Hoffmann is the basis of the Woodward-Hoffmann rules.[3] It was later : extended to conjugated molecules such as pyridine, pyrrole and furan that
|
h********r 发帖数: 821 | 75 这里计算的是opt+freq,还是freq本身?数值梯度不准可能导致opt多做几圈。
【在 F****I 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 用了。再给一个没有空间对称性的,乙醇,HF/cc-pVDZ opt + freq, Gaussian09 : 解析频率 : Job cpu time: 0 days 0 hours 0 minutes 17.0 seconds. : 数值频率 : Job cpu time: 0 days 0 hours 2 minutes 56.0 seconds.
|
h********r 发帖数: 821 | 76 我猜小分子的时候prefactor重要,分子变大scaling就重要了,结果也许会颠倒(待验
证)。
多谢你的辛勤计算!
【在 F****I 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 好吧好吧 : 乙醇,MP2/cc-pVDZ, Gaussian09 : 解析频率 : Job cpu time: 0 days 0 hours 1 minutes 56.0 seconds. : 数值频率 : Job cpu time: 0 days 0 hours 8 minutes 25.0 seconds. : 更大基组自行测试吧
|
F****I 发帖数: 270 | 77 opt+freq. opt用的解析梯度,freq用解析或者数值
【在 h********r 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 这里计算的是opt+freq,还是freq本身?数值梯度不准可能导致opt多做几圈。
|
t****a 发帖数: 3544 | 78 赞这句~~~
【在 a****o 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 没有关系。好多搞工程的,还停留在牛顿力学的水平呢。
|
t****a 发帖数: 3544 | 79 赞~~ 过于数学化是数学物理了
G
【在 h*****y 的大作中提到】![](/moin_static193/solenoid/img/up.png) : DFT已经是凝聚态物理里研究材料的标准工具了。model hamitonian最后不回到数值解 : 的有几个?现在计算机的计算能力越来越强,不管凝聚态还是化学,ab initio是个趋 : 势。强关联问题现在DFT还不能完全处理,这反过来也说明这个领域大有可为。DMFT, G : -W, S-B这些都是可以继续完善的方案。得到精确ec也不是完全没有办法,LDA方案中的 : ec参数就可以用QMC算出来。回头说说传统理论物理,现在的趋势是数学越来越fancy, : 但是在问题的实质上并没有新的突破。一个例子就是 ferroelectricity,多少年了都是 : Ginzburg–Landau theory,DFT一下子就把认识推进到了微观层次。不管物理还是化学 : ,理论还是计算,只要有新见解,就是好的研究,只是打着纯理论的幌子,即使方法再 : 高深,数学推导再繁琐,也是不可取的。 :
|
F****I 发帖数: 270 | 80 插一句,如果不对请指正。
> HEG的基态能是解析的,QMC做了精确的参数化,
好像HEG的相关能(作为基态能的一部分)只有低密度和高密度极限的解析行为是已知的,一般情况下只能猜一个表达式,并且保证正确的渐进行为。VWN,PW92都是这样。因而Exc的精度被猜出的表达式的泛函形式限制。或者您所说的解析指复变函数里的解析函数的那种解析? 还是说能够写成一个closed-form就叫解析(这个意义下好像也没有证明)
>带来的问题就是引入了平均场近似。
只要有库仑相关就超过了平均场近似的水平
【。 在 hahoney (hahoney) 的大作中提到: 】 |
|
|
a*****9 发帖数: 1884 | 81
严重握手
【在 F****I 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 我一直很佩服做有机实验的人啊。自己做的时候产率经常为零。开学的时候交了玻璃仪 : 器的押金,学期末都成了负数
|