search

maryprimary / mypydmet

dmet算法测试
0 stars 2 forks source link

来请教觉主了 #1

Open bsonghao opened 2 years ago

bsonghao commented 2 years ago

觉主，我对这个DMET的repo比较感兴趣，可否有空给我讲讲？最近在看Garnet Chan 的论文

maryprimary commented 2 years ago

这个代码当时没有写完就被搁置了，主要是求解mu的时候功能没有实现，我觉得dmet这个算法主要的部分就是构建两个轨道，把这部分写完就差不多了。你可以看一下dmet/init.py里面的update_umat，这一部分写完其实就差不多了，里面涉及BathOrbital。 我是根据这篇文章写的： A Practical Guide to Density Matrix Embedding Theory in Quantum Chemistry. JCTC. Sebastian Wouters

70akaline commented 2 years ago

这个代码当时没有写完就被搁置了，主要是求解mu的时候功能没有实现，我觉得dmet这个算法主要的部分就是构建两个轨道，把这部分写完就差不多了。你可以看一下dmet/init.py里面的update_umat，这一部分写完其实就差不多了，里面涉及BathOrbital。 我是根据这篇文章写的： A Practical Guide to Density Matrix Embedding Theory in Quantum Chemistry. JCTC. Sebastian Wouters

跟楼一下问问 tianyu zhu最近的gw+dmft的文章中提到用dmrg来作为impurity solver 这类dmrg算法可以支持完整的four index Coulomb tersor么？

bsonghao commented 2 years ago

image @maryprimary 我发现你这里的impurity list是你自己定义的？这和原文描述的的不太相符？原文是根据overlap matrix确定fragment的？

maryprimary commented 2 years ago

这个imp_list相当于公式里面的Ax, \sum{A_x} 就相当于在imp_list里面求和 @bsonghao

maryprimary commented 2 years ago

@70akaline 理论上讲是可以的，DMRG会保留所有的算符，只是缩Hilbert空间的大小

bsonghao commented 2 years ago

我知道，但是A和B是你自己定义的？