From e9101e2b237bd4c200df008da235955574883372 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Zhanghao=20Zhouyin=20=28=E5=91=A8=E5=AF=85=E5=BC=A0?= =?UTF-8?q?=E7=9A=93=29?= <64253517+floatingCatty@users.noreply.github.com> Date: Sun, 11 Oct 2026 20:14:41 -0400 Subject: [PATCH 01/20] Publish ad-negf original English article at its canonical URL --- _pages/research-ad-negf.md | 16 ++++------------ 1 file changed, 4 insertions(+), 12 deletions(-) diff --git a/_pages/research-ad-negf.md b/_pages/research-ad-negf.md index d993276..756caf2 100644 --- a/_pages/research-ad-negf.md +++ b/_pages/research-ad-negf.md @@ -8,20 +8,12 @@ description_zh: "AD-NEGF 为紧束缚量子输运计算增加梯度,用于灵 paper_year: 2023 paper_url: "https://doi.org/10.1103/PhysRevB.108.195143" code_url: "https://github.com/floatingCatty/ADNEGF" +research_locale: en +translation_url: "/zh/research/ad-negf/" +alternate_en_url: "/research/ad-negf/" +alternate_zh_url: "/zh/research/ad-negf/" --- -
- -我在天津大学做量子输运的时候,常常用 NEGF(非平衡格林函数)方法算电子通过一个器件的透射和电流。给定描述电子能量和耦合的哈密顿量,再加上电极和偏压,程序就能算出结果。但如果想知道某个参数改一点,电流会变多少,往往只能改完再跑一次。参数多的时候,这个办法很慢。 - -我们当时想试试,能不能直接对整个输运计算求导。比如,电流对某个轨道耦合、原子位置或其他输入变量的梯度。有了这些导数,可以分析模型对什么最敏感,也可以尝试反过来调整哈密顿量,让它接近某个目标输运性质。 - -AD-NEGF 用 PyTorch 实现了 NEGF 计算流程,并针对部分求解步骤写了专门的反向传播。NEGF 里有矩阵求逆和需要数值求解的中间量,直接依赖一般的自动微分不一定合适,所以论文中也用了隐式求导的方法。我们验证了梯度计算,还展示了从目标输运性质反过来优化哈密顿量的例子。 - -不过,这里的设计主要还在模型层面。哈密顿量参数怎样对应到实际能够制造的器件结构,散射、缺陷、接触等效应怎么处理,都需要另外建模。能对求解器求导,确实方便了很多原本只能反复试算的工作,但离直接用于真实器件设计,还有不少事情要做。 - -
-
When I was working on quantum transport at Tianjin University, I often used the non-equilibrium Green's function (NEGF) method to calculate electron transmission and current through a device. Given a Hamiltonian describing electronic energies and couplings, along with electrodes and an applied bias, the program produces a result. But if we want to know how much the current changes when a parameter is adjusted, we usually have to run the calculation again. That gets slow when there are many parameters. From 80307810bf05dfdbcf7834d0c21fc7768c1890b7 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Zhanghao=20Zhouyin=20=28=E5=91=A8=E5=AF=85=E5=BC=A0?= =?UTF-8?q?=E7=9A=93=29?= <64253517+floatingCatty@users.noreply.github.com> Date: Sun, 11 Oct 2026 20:14:44 -0400 Subject: [PATCH 02/20] Add full Chinese ad-negf article with reciprocal hreflang --- _pages/research-ad-negf-zh.md | 28 ++++++++++++++++++++++++++++ 1 file changed, 28 insertions(+) create mode 100644 _pages/research-ad-negf-zh.md diff --git a/_pages/research-ad-negf-zh.md b/_pages/research-ad-negf-zh.md new file mode 100644 index 0000000..89f7bd9 --- /dev/null +++ b/_pages/research-ad-negf-zh.md @@ -0,0 +1,28 @@ +--- +layout: research +permalink: /zh/research/ad-negf/ +title: "AD-NEGF:可微的量子输运计算" +title_en: "AD-NEGF: differentiable quantum transport calculations" +description: "AD-NEGF computes gradients through tight-binding quantum-transport calculations for sensitivity analysis and inverse Hamiltonian design." +description_zh: "AD-NEGF 为紧束缚量子输运计算增加梯度,用于灵敏度分析和哈密顿量逆向设计。" +paper_year: 2023 +paper_url: "https://doi.org/10.1103/PhysRevB.108.195143" +code_url: "https://github.com/floatingCatty/ADNEGF" +research_locale: zh +translation_url: "/research/ad-negf/" +translation_of: "/research/ad-negf/" +alternate_en_url: "/research/ad-negf/" +alternate_zh_url: "/zh/research/ad-negf/" +--- + +
+ +我在天津大学做量子输运的时候,常常用 NEGF(非平衡格林函数)方法算电子通过一个器件的透射和电流。给定描述电子能量和耦合的哈密顿量,再加上电极和偏压,程序就能算出结果。但如果想知道某个参数改一点,电流会变多少,往往只能改完再跑一次。参数多的时候,这个办法很慢。 + +我们当时想试试,能不能直接对整个输运计算求导。比如,电流对某个轨道耦合、原子位置或其他输入变量的梯度。有了这些导数,可以分析模型对什么最敏感,也可以尝试反过来调整哈密顿量,让它接近某个目标输运性质。 + +AD-NEGF 用 PyTorch 实现了 NEGF 计算流程,并针对部分求解步骤写了专门的反向传播。NEGF 里有矩阵求逆和需要数值求解的中间量,直接依赖一般的自动微分不一定合适,所以论文中也用了隐式求导的方法。我们验证了梯度计算,还展示了从目标输运性质反过来优化哈密顿量的例子。 + +不过,这里的设计主要还在模型层面。哈密顿量参数怎样对应到实际能够制造的器件结构,散射、缺陷、接触等效应怎么处理,都需要另外建模。能对求解器求导,确实方便了很多原本只能反复试算的工作,但离直接用于真实器件设计,还有不少事情要做。 + +
From b7b1b46d0d3c1d867b5d1f3fc9f104a69ce87db1 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Zhanghao=20Zhouyin=20=28=E5=91=A8=E5=AF=85=E5=BC=A0?= =?UTF-8?q?=E7=9A=93=29?= <64253517+floatingCatty@users.noreply.github.com> Date: Sun, 11 Oct 2026 20:14:48 -0400 Subject: [PATCH 03/20] Publish deeptb-negf original English article at its canonical URL --- _pages/research-deeptb-negf.md | 16 ++++------------ 1 file changed, 4 insertions(+), 12 deletions(-) diff --git a/_pages/research-deeptb-negf.md b/_pages/research-deeptb-negf.md index 477c973..c2d020c 100644 --- a/_pages/research-deeptb-negf.md +++ b/_pages/research-deeptb-negf.md @@ -8,20 +8,12 @@ description_zh: "DeePTB 和 DPNEGF 结合,计算石墨烯、氮化硼和二硫 paper_year: 2025 paper_url: "https://arxiv.org/abs/2512.11291" code_url: "https://github.com/deepmodeling/dpnegf" +research_locale: en +translation_url: "/zh/research/deeptb-negf/" +alternate_en_url: "/research/deeptb-negf/" +alternate_zh_url: "/zh/research/deeptb-negf/" --- -
- -做量子输运计算时,NEGF 通常需要先拿到体系的哈密顿矩阵。用 DFT 直接构造当然比较可靠,但材料一大,或者想换很多个结构反复计算,这一步往往就是主要开销。DeePTB 能从原子结构预测紧束缚哈密顿量,于是一个自然的问题是:能不能把它直接接到量子输运求解器里? - -这篇合作工作尝试了这件事。我们把 DeePTB 预测的哈密顿量接到 DPNEGF,研究石墨烯、六方氮化硼和 MoS₂ 这三种二维材料,并与传统 DFT-NEGF 得到的能带和电子透射谱做比较。预印本报告了较好的一致性,以及在这些测试体系上的计算效率提升。 - -这里实际上有两层需要分开验证的近似。一层是学习到的哈密顿量是否准确,另一层是输运求解器能否正确处理接触、器件几何和边界条件。即使能带看起来相近,透射也不一定自然就会准确。所以我觉得把模型和物性计算接在一起,并做完整的输运对照,比单独展示一个好的哈密顿量误差更有意义。 - -目前公开的是预印本。它说明了这条计算路线在几个典型二维材料上的可行性,但不能据此说所有缺陷、界面或实际器件都已有同等的精度。 - -
-
A quantum transport calculation usually begins with a Hamiltonian. Constructing it directly from DFT can be reliable, but it becomes expensive when the system is large or when we want to compare many different structures. DeePTB can predict a tight-binding Hamiltonian from an atomic structure, which raises a natural question: can we use that Hamiltonian directly in a quantum transport solver? From dae476a363a4509df8b54348c89addf2c4237ed2 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Zhanghao=20Zhouyin=20=28=E5=91=A8=E5=AF=85=E5=BC=A0?= =?UTF-8?q?=E7=9A=93=29?= <64253517+floatingCatty@users.noreply.github.com> Date: Sun, 11 Oct 2026 20:14:51 -0400 Subject: [PATCH 04/20] Add full Chinese deeptb-negf article with reciprocal hreflang --- _pages/research-deeptb-negf-zh.md | 28 ++++++++++++++++++++++++++++ 1 file changed, 28 insertions(+) create mode 100644 _pages/research-deeptb-negf-zh.md diff --git a/_pages/research-deeptb-negf-zh.md b/_pages/research-deeptb-negf-zh.md new file mode 100644 index 0000000..d5e0457 --- /dev/null +++ b/_pages/research-deeptb-negf-zh.md @@ -0,0 +1,28 @@ +--- +layout: research +permalink: /zh/research/deeptb-negf/ +title: "DeePTB-NEGF:从学习到的哈密顿量计算量子输运" +title_en: "DeePTB-NEGF: transport from learned Hamiltonians" +description: "A DeePTB and DPNEGF workflow for quantum transport in graphene, h-BN and MoS2, compared with first-principles DFT-NEGF calculations." +description_zh: "DeePTB 和 DPNEGF 结合,计算石墨烯、氮化硼和二硫化钼的量子输运,并与 DFT-NEGF 对照。" +paper_year: 2025 +paper_url: "https://arxiv.org/abs/2512.11291" +code_url: "https://github.com/deepmodeling/dpnegf" +research_locale: zh +translation_url: "/research/deeptb-negf/" +translation_of: "/research/deeptb-negf/" +alternate_en_url: "/research/deeptb-negf/" +alternate_zh_url: "/zh/research/deeptb-negf/" +--- + +
+ +做量子输运计算时,NEGF 通常需要先拿到体系的哈密顿矩阵。用 DFT 直接构造当然比较可靠,但材料一大,或者想换很多个结构反复计算,这一步往往就是主要开销。DeePTB 能从原子结构预测紧束缚哈密顿量,于是一个自然的问题是:能不能把它直接接到量子输运求解器里? + +这篇合作工作尝试了这件事。我们把 DeePTB 预测的哈密顿量接到 DPNEGF,研究石墨烯、六方氮化硼和 MoS₂ 这三种二维材料,并与传统 DFT-NEGF 得到的能带和电子透射谱做比较。预印本报告了较好的一致性,以及在这些测试体系上的计算效率提升。 + +这里实际上有两层需要分开验证的近似。一层是学习到的哈密顿量是否准确,另一层是输运求解器能否正确处理接触、器件几何和边界条件。即使能带看起来相近,透射也不一定自然就会准确。所以我觉得把模型和物性计算接在一起,并做完整的输运对照,比单独展示一个好的哈密顿量误差更有意义。 + +目前公开的是预印本。它说明了这条计算路线在几个典型二维材料上的可行性,但不能据此说所有缺陷、界面或实际器件都已有同等的精度。 + +
From 75534bc0f5dc89151668241cfdf4c2e7fc102452 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Zhanghao=20Zhouyin=20=28=E5=91=A8=E5=AF=85=E5=BC=A0?= =?UTF-8?q?=E7=9A=93=29?= <64253517+floatingCatty@users.noreply.github.com> Date: Sun, 11 Oct 2026 20:14:55 -0400 Subject: [PATCH 05/20] Publish ferh-magnetism original English article at its canonical URL --- _pages/research-ferh-magnetism.md | 14 ++++---------- 1 file changed, 4 insertions(+), 10 deletions(-) diff --git a/_pages/research-ferh-magnetism.md b/_pages/research-ferh-magnetism.md index 5ce63b0..d4ce280 100644 --- a/_pages/research-ferh-magnetism.md +++ b/_pages/research-ferh-magnetism.md @@ -7,18 +7,12 @@ description: "A note on why density-functional predictions of magnetism in FeRh description_zh: "关于 FeRh 合金的磁性预测为什么依赖晶格常数和 DFT 交换关联泛函。" paper_year: 2025 paper_url: "https://doi.org/10.1016/j.commatsci.2024.113561" +research_locale: en +translation_url: "/zh/research/ferh-magnetism/" +alternate_en_url: "/research/ferh-magnetism/" +alternate_zh_url: "/zh/research/ferh-magnetism/" --- -
- -FeRh 是铁和铑组成的一种合金,磁性对晶格常数特别敏感。原子之间的距离只要有一点变化,不同磁性状态之间的能量关系就可能跟着改变。这给 DFT 计算出了一个难题:算出来的结果,到底有多少来自材料本身,有多少来自我们使用的近似? - -在 DFT 里,电子的交换关联效应需要用一个泛函来近似。实际做计算时,可选的方法有很多,例如 LDA、GGA 和 SCAN。它们在一些材料上可以给出很接近的结果,换一个体系却不一定如此。FeRh 就是一个很适合拿来比较的例子,因为结构、电子和磁性之间联系得很紧。 - -这篇工作比较了几种交换关联泛函对 FeRh 的描述,包括晶格常数和磁性状态等结果。我们发现,没有哪一个泛函能在所有相关性质上都给出同样可靠的描述。因此研究 FeRh 这类体系时,需要关心计算对泛函和结构参数的依赖,不能只报出一次计算得到的磁性结果。 - -
-
FeRh is an alloy of iron and rhodium whose magnetism is unusually sensitive to the lattice constant. A small change in the distance between atoms can change the relative energies of different magnetic states. That makes it a difficult system for DFT: how much of a prediction comes from the material, and how much comes from the approximation we chose? From 3d64b87fb34762d57b2731512021adce46ed7f10 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Zhanghao=20Zhouyin=20=28=E5=91=A8=E5=AF=85=E5=BC=A0?= =?UTF-8?q?=E7=9A=93=29?= <64253517+floatingCatty@users.noreply.github.com> Date: Sun, 11 Oct 2026 20:14:58 -0400 Subject: [PATCH 06/20] Add full Chinese ferh-magnetism article with reciprocal hreflang --- _pages/research-ferh-magnetism-zh.md | 25 +++++++++++++++++++++++++ 1 file changed, 25 insertions(+) create mode 100644 _pages/research-ferh-magnetism-zh.md diff --git a/_pages/research-ferh-magnetism-zh.md b/_pages/research-ferh-magnetism-zh.md new file mode 100644 index 0000000..5484e62 --- /dev/null +++ b/_pages/research-ferh-magnetism-zh.md @@ -0,0 +1,25 @@ +--- +layout: research +permalink: /zh/research/ferh-magnetism/ +title: "FeRh 中交换关联泛函对磁性的影响" +title_en: "Exchange-correlation functionals and magnetism in FeRh" +description: "A note on why density-functional predictions of magnetism in FeRh change with lattice constant and exchange-correlation functional." +description_zh: "关于 FeRh 合金的磁性预测为什么依赖晶格常数和 DFT 交换关联泛函。" +paper_year: 2025 +paper_url: "https://doi.org/10.1016/j.commatsci.2024.113561" +research_locale: zh +translation_url: "/research/ferh-magnetism/" +translation_of: "/research/ferh-magnetism/" +alternate_en_url: "/research/ferh-magnetism/" +alternate_zh_url: "/zh/research/ferh-magnetism/" +--- + +
+ +FeRh 是铁和铑组成的一种合金,磁性对晶格常数特别敏感。原子之间的距离只要有一点变化,不同磁性状态之间的能量关系就可能跟着改变。这给 DFT 计算出了一个难题:算出来的结果,到底有多少来自材料本身,有多少来自我们使用的近似? + +在 DFT 里,电子的交换关联效应需要用一个泛函来近似。实际做计算时,可选的方法有很多,例如 LDA、GGA 和 SCAN。它们在一些材料上可以给出很接近的结果,换一个体系却不一定如此。FeRh 就是一个很适合拿来比较的例子,因为结构、电子和磁性之间联系得很紧。 + +这篇工作比较了几种交换关联泛函对 FeRh 的描述,包括晶格常数和磁性状态等结果。我们发现,没有哪一个泛函能在所有相关性质上都给出同样可靠的描述。因此研究 FeRh 这类体系时,需要关心计算对泛函和结构参数的依赖,不能只报出一次计算得到的磁性结果。 + +
From b314625a20826b0d00377da1f21f0246097cbe95 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Zhanghao=20Zhouyin=20=28=E5=91=A8=E5=AF=85=E5=BC=A0?= =?UTF-8?q?=E7=9A=93=29?= <64253517+floatingCatty@users.noreply.github.com> Date: Sun, 11 Oct 2026 20:15:02 -0400 Subject: [PATCH 07/20] Publish logdet-entropy original English article at its canonical URL --- _pages/research-logdet-entropy.md | 16 ++++------------ 1 file changed, 4 insertions(+), 12 deletions(-) diff --git a/_pages/research-logdet-entropy.md b/_pages/research-logdet-entropy.md index 2e45b33..b589aa2 100644 --- a/_pages/research-logdet-entropy.md +++ b/_pages/research-logdet-entropy.md @@ -12,20 +12,12 @@ paper_doi: "10.1016/j.neucom.2025.131520" paper_authors: - "Zhanghao Zhouyin" - "Ding Liu" +research_locale: en +translation_url: "/zh/research/logdet-entropy/" +alternate_en_url: "/research/logdet-entropy/" +alternate_zh_url: "/zh/research/logdet-entropy/" --- -
- -我本科时做过一段时间神经网络理论,最感兴趣的是网络内部的信息到底发生了什么变化。比如一层层特征传递下去,哪些东西保留下来了,哪些被丢掉了。信息论给我们提供了熵和互信息这些工具,但真要把它们用在高维神经网络特征上,估计本身就不容易。 - -直接估计高维概率分布需要大量样本。一些常见的熵估计器,在维度很高或者数据分布比较特殊时也容易出现偏差。我们于是试了一种从矩阵出发的办法:利用样本构成的矩阵及其对数行列式(LogDet)来估计信息量,省去直接重建完整概率密度的步骤。 - -论文中我们做了几组对比实验,也把估计器用来观察神经网络训练过程中不同层的特征变化,希望借此分析信息瓶颈理论里讨论的现象。这篇工作很早就开始了,后来经过几轮修改,发表在 2025 年的 Neurocomputing。 - -这类指标的意义也要放在具体估计方法里看。它可以帮助比较和观察网络,但不能简单地把一个估计数值当作对网络内部信息的完整解释。我后来转到物理计算,仍然常常碰到类似的问题:先要弄清楚我们用的量到底测量了什么,才谈得上从结果里得出结论。 - -
-
As an undergraduate, I spent some time studying neural-network theory. I was interested in what happens to information inside a network: what survives as features pass through layers, and what gets discarded. Information theory gives us entropy and mutual information, but estimating them for high-dimensional neural-network features is difficult in practice. From 4548bf60b84f02df1e43b6806e54a765308990af Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Zhanghao=20Zhouyin=20=28=E5=91=A8=E5=AF=85=E5=BC=A0?= =?UTF-8?q?=E7=9A=93=29?= <64253517+floatingCatty@users.noreply.github.com> Date: Sun, 11 Oct 2026 20:15:05 -0400 Subject: [PATCH 08/20] Add full Chinese logdet-entropy article with reciprocal hreflang --- _pages/research-logdet-entropy-zh.md | 32 ++++++++++++++++++++++++++++ 1 file changed, 32 insertions(+) create mode 100644 _pages/research-logdet-entropy-zh.md diff --git a/_pages/research-logdet-entropy-zh.md b/_pages/research-logdet-entropy-zh.md new file mode 100644 index 0000000..1b1d25f --- /dev/null +++ b/_pages/research-logdet-entropy-zh.md @@ -0,0 +1,32 @@ +--- +layout: research +permalink: /zh/research/logdet-entropy/ +title: "LogDet:用矩阵估计神经网络里的信息量" +title_en: "LogDet: estimating information in neural networks" +description: "The LogDet entropy estimator measures high-dimensional information using matrix operations and was applied to the analysis of neural networks." +description_zh: "LogDet 熵估计器利用矩阵运算处理高维数据,用于分析神经网络特征中的信息变化。" +paper_year: 2025 +paper_url: "https://doi.org/10.1016/j.neucom.2025.131520" +paper_title: "Understanding neural networks with logarithm determinant entropy estimator" +paper_doi: "10.1016/j.neucom.2025.131520" +paper_authors: + - "Zhanghao Zhouyin" + - "Ding Liu" +research_locale: zh +translation_url: "/research/logdet-entropy/" +translation_of: "/research/logdet-entropy/" +alternate_en_url: "/research/logdet-entropy/" +alternate_zh_url: "/zh/research/logdet-entropy/" +--- + +
+ +我本科时做过一段时间神经网络理论,最感兴趣的是网络内部的信息到底发生了什么变化。比如一层层特征传递下去,哪些东西保留下来了,哪些被丢掉了。信息论给我们提供了熵和互信息这些工具,但真要把它们用在高维神经网络特征上,估计本身就不容易。 + +直接估计高维概率分布需要大量样本。一些常见的熵估计器,在维度很高或者数据分布比较特殊时也容易出现偏差。我们于是试了一种从矩阵出发的办法:利用样本构成的矩阵及其对数行列式(LogDet)来估计信息量,省去直接重建完整概率密度的步骤。 + +论文中我们做了几组对比实验,也把估计器用来观察神经网络训练过程中不同层的特征变化,希望借此分析信息瓶颈理论里讨论的现象。这篇工作很早就开始了,后来经过几轮修改,发表在 2025 年的 Neurocomputing。 + +这类指标的意义也要放在具体估计方法里看。它可以帮助比较和观察网络,但不能简单地把一个估计数值当作对网络内部信息的完整解释。我后来转到物理计算,仍然常常碰到类似的问题:先要弄清楚我们用的量到底测量了什么,才谈得上从结果里得出结论。 + +
From ebdd7468b68d9d41aa82a2389b2b34f3c9061a6b Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Zhanghao=20Zhouyin=20=28=E5=91=A8=E5=AF=85=E5=BC=A0?= =?UTF-8?q?=E7=9A=93=29?= <64253517+floatingCatty@users.noreply.github.com> Date: Sun, 11 Oct 2026 20:15:15 -0400 Subject: [PATCH 09/20] Preserve existing English prose at /research/moire-phonon/ --- _pages/research-moire-phonon.md | 16 ++++------------ 1 file changed, 4 insertions(+), 12 deletions(-) diff --git a/_pages/research-moire-phonon.md b/_pages/research-moire-phonon.md index 152a2cb..48b1a18 100644 --- a/_pages/research-moire-phonon.md +++ b/_pages/research-moire-phonon.md @@ -14,20 +14,12 @@ paper_authors: - "Jyun-jie Jiang" - "Xianghua Kong" - "Hong Guo" +research_locale: en +translation_url: "/zh/research/moire-phonon/" +alternate_en_url: "/research/moire-phonon/" +alternate_zh_url: "/zh/research/moire-phonon/" --- -
- -把两层石墨烯叠在一起,稍稍转过一个角度,便会出现比原本晶格大得多的莫尔图案。角度越小,莫尔周期越长,原子也会跟着重新排列。我最近做的一项工作,就是想弄清楚魔角附近这种结构变化究竟是怎么发生的。 - -我们最开始看到的是两种不同的形变。离魔角稍远一些,两层石墨烯主要表现为层间距离的起伏;接近魔角后,却逐渐变成两层一起弯曲。用一个包含上万个原子的超胞直接看位移,图很复杂,也不容易知道到底是哪些自由度在起作用。后来我们把位移投影到声子模式上,发现事情简单了许多。在 1.08° 的体系里,11,164 个原子的重构位移,超过 99.5% 的谱权重可以由两个 A1 对称性的声子模式描述。 - -沿着这两个模式继续分析,会看到其中的层对称弯曲模式随着转角减小逐渐软化,也就是结构沿这个方向变形所需要的恢复力越来越弱。到了魔角附近,原来稳定的形态便倾向于沿这个模式发生静态形变。我们把它叫作莫尔声子凝聚(moiré phonon condensation)。这里的“凝聚”指软声子模式形成静态位移,不是说材料里发生了玻色-爱因斯坦凝聚。 - -我们还建立了一个比较简单的连续介质模型,把这种变化和莫尔周期增长后应力与弯曲能之间的竞争联系起来。紧束缚计算也显示,这些形变模式会影响电子结构。现在论文还是预印本,关于实验中怎样直接观察这些模式,还需要进一步研究。 - -
-
Put two sheets of graphene on top of each other and rotate one slightly. A moiré pattern appears, much larger than the original atomic lattice. As the twist angle gets smaller, the moiré period grows and the atoms rearrange. In a recent project, I wanted to understand what causes this reconstruction near the magic angle. From be2c6c7004bd5bb8841d52d320fbd2d460536506 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Zhanghao=20Zhouyin=20=28=E5=91=A8=E5=AF=85=E5=BC=A0?= =?UTF-8?q?=E7=9A=93=29?= <64253517+floatingCatty@users.noreply.github.com> Date: Sun, 11 Oct 2026 20:15:19 -0400 Subject: [PATCH 10/20] Add Chinese moire-phonon research article with reciprocal hreflang --- _pages/research-moire-phonon-zh.md | 34 ++++++++++++++++++++++++++++++ 1 file changed, 34 insertions(+) create mode 100644 _pages/research-moire-phonon-zh.md diff --git a/_pages/research-moire-phonon-zh.md b/_pages/research-moire-phonon-zh.md new file mode 100644 index 0000000..3f26b7e --- /dev/null +++ b/_pages/research-moire-phonon-zh.md @@ -0,0 +1,34 @@ +--- +layout: research +permalink: /zh/research/moire-phonon/ +title: "魔角石墨烯里的声子凝聚" +title_en: "Moiré phonon condensation in twisted bilayer graphene" +description: "A note on moiré phonon condensation, soft modes, and structural reconstruction in magic-angle twisted bilayer graphene." +description_zh: "关于魔角双层石墨烯结构重构、软声子模式和莫尔声子凝聚的一些计算与想法。" +paper_year: 2026 +paper_url: "https://arxiv.org/abs/2607.06711" +paper_title: "Moiré Phonon Condensation in Magic-Angle Twisted Bilayer Graphene" +paper_doi: "10.48550/arXiv.2607.06711" +paper_authors: + - "Zhanghao Zhouyin" + - "Jyun-jie Jiang" + - "Xianghua Kong" + - "Hong Guo" +research_locale: zh +translation_url: "/research/moire-phonon/" +translation_of: "/research/moire-phonon/" +alternate_en_url: "/research/moire-phonon/" +alternate_zh_url: "/zh/research/moire-phonon/" +--- + +
+ +把两层石墨烯叠在一起,稍稍转过一个角度,便会出现比原本晶格大得多的莫尔图案。角度越小,莫尔周期越长,原子也会跟着重新排列。我最近做的一项工作,就是想弄清楚魔角附近这种结构变化究竟是怎么发生的。 + +我们最开始看到的是两种不同的形变。离魔角稍远一些,两层石墨烯主要表现为层间距离的起伏;接近魔角后,却逐渐变成两层一起弯曲。用一个包含上万个原子的超胞直接看位移,图很复杂,也不容易知道到底是哪些自由度在起作用。后来我们把位移投影到声子模式上,发现事情简单了许多。在 1.08° 的体系里,11,164 个原子的重构位移,超过 99.5% 的谱权重可以由两个 A1 对称性的声子模式描述。 + +沿着这两个模式继续分析,会看到其中的层对称弯曲模式随着转角减小逐渐软化,也就是结构沿这个方向变形所需要的恢复力越来越弱。到了魔角附近,原来稳定的形态便倾向于沿这个模式发生静态形变。我们把它叫作莫尔声子凝聚(moiré phonon condensation)。这里的“凝聚”指软声子模式形成静态位移,不是说材料里发生了玻色-爱因斯坦凝聚。 + +我们还建立了一个比较简单的连续介质模型,把这种变化和莫尔周期增长后应力与弯曲能之间的竞争联系起来。紧束缚计算也显示,这些形变模式会影响电子结构。现在论文还是预印本,关于实验中怎样直接观察这些模式,还需要进一步研究。 + +
From 1ce9f965f27b6de79bb82d3ee96dc828392546ad Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Zhanghao=20Zhouyin=20=28=E5=91=A8=E5=AF=85=E5=BC=A0?= =?UTF-8?q?=E7=9A=93=29?= <64253517+floatingCatty@users.noreply.github.com> Date: Sun, 11 Oct 2026 20:15:26 -0400 Subject: [PATCH 11/20] Preserve existing English prose at /research/nqs-impurity/ --- _pages/research-nqs-impurity.md | 16 ++++------------ 1 file changed, 4 insertions(+), 12 deletions(-) diff --git a/_pages/research-nqs-impurity.md b/_pages/research-nqs-impurity.md index 941cb04..f643e92 100644 --- a/_pages/research-nqs-impurity.md +++ b/_pages/research-nqs-impurity.md @@ -15,20 +15,12 @@ paper_authors: - "Ao Chen" - "Nicola Lanatà" - "Hong Guo" +research_locale: en +translation_url: "/zh/research/nqs-impurity/" +alternate_en_url: "/research/nqs-impurity/" +alternate_zh_url: "/zh/research/nqs-impurity/" --- -
- -电子之间的相互作用比较强时,计算材料的电子状态会变得很困难。一个直接的办法是把许多电子放在一起求解,但体系稍微大一点,可能的量子态数量就增长得很快。量子嵌入方法换了一种处理方式:把其中最难的一小部分当作杂质问题来算,再让它与周围的有效环境反复交换信息。 - -我做的这篇工作,考虑的是 ghost Gutzwiller approximation(gGA)中的杂质求解器。通常可以用精确对角化作为参照,只是它的计算开销随杂质自由度增加会变得很大。我们尝试用神经量子态(NQS)表示多电子波函数,网络采用图 Transformer,让杂质轨道之间不必拘泥于某一种固定连接方式。 - -麻烦出现在嵌入计算的迭代里。杂质求解器一次得到的结果还要送回外层自洽过程,电荷和其他物理量只要有一点采样噪声,误差就可能一轮轮积累。我们因此加入了误差控制,在 Anderson 晶格模型的测试中,结果可以与使用精确对角化杂质求解器的计算很好地吻合。 - -原先我也以为主要时间会花在训练神经网络上。实际分析计算开销后发现,对这个问题来说,更难的是把需要的物理量采样到足够高的精度。神经量子态本身能表示复杂波函数,并不意味着自洽计算就会很便宜。怎样减少高精度采样的成本,是这项工作留下的问题。 - -
-
When electrons interact strongly, calculating the electronic state of a material becomes difficult. We could try solving for many electrons at once, but the number of possible quantum states grows rapidly with system size. Quantum embedding takes a different approach. It treats a smaller, difficult part as an impurity problem and repeatedly exchanges information between that problem and an effective environment. From b1c92f735d78f31f513970348743701b8b798afb Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Zhanghao=20Zhouyin=20=28=E5=91=A8=E5=AF=85=E5=BC=A0?= =?UTF-8?q?=E7=9A=93=29?= <64253517+floatingCatty@users.noreply.github.com> Date: Sun, 11 Oct 2026 20:15:29 -0400 Subject: [PATCH 12/20] Add Chinese nqs-impurity research article with reciprocal hreflang --- _pages/research-nqs-impurity-zh.md | 35 ++++++++++++++++++++++++++++++ 1 file changed, 35 insertions(+) create mode 100644 _pages/research-nqs-impurity-zh.md diff --git a/_pages/research-nqs-impurity-zh.md b/_pages/research-nqs-impurity-zh.md new file mode 100644 index 0000000..c7bbcb9 --- /dev/null +++ b/_pages/research-nqs-impurity-zh.md @@ -0,0 +1,35 @@ +--- +layout: research +permalink: /zh/research/nqs-impurity/ +title: "用神经网络求解量子嵌入里的杂质问题" +title_en: "Neural quantum states for quantum embedding" +description: "A research note on a neural-quantum-state impurity solver in ghost Gutzwiller quantum embedding and its sampling bottleneck." +description_zh: "用神经量子态求解 ghost Gutzwiller 量子嵌入中的杂质模型,以及高精度采样带来的计算开销。" +paper_year: 2026 +paper_url: "https://doi.org/10.1103/rkd8-q6yl" +paper_title: "Neural-quantum-states impurity solver for quantum embedding problems" +paper_doi: "10.1103/rkd8-q6yl" +paper_authors: + - "Zhanghao Zhouyin" + - "Tsung-Han Lee" + - "Ao Chen" + - "Nicola Lanatà" + - "Hong Guo" +research_locale: zh +translation_url: "/research/nqs-impurity/" +translation_of: "/research/nqs-impurity/" +alternate_en_url: "/research/nqs-impurity/" +alternate_zh_url: "/zh/research/nqs-impurity/" +--- + +
+ +电子之间的相互作用比较强时,计算材料的电子状态会变得很困难。一个直接的办法是把许多电子放在一起求解,但体系稍微大一点,可能的量子态数量就增长得很快。量子嵌入方法换了一种处理方式:把其中最难的一小部分当作杂质问题来算,再让它与周围的有效环境反复交换信息。 + +我做的这篇工作,考虑的是 ghost Gutzwiller approximation(gGA)中的杂质求解器。通常可以用精确对角化作为参照,只是它的计算开销随杂质自由度增加会变得很大。我们尝试用神经量子态(NQS)表示多电子波函数,网络采用图 Transformer,让杂质轨道之间不必拘泥于某一种固定连接方式。 + +麻烦出现在嵌入计算的迭代里。杂质求解器一次得到的结果还要送回外层自洽过程,电荷和其他物理量只要有一点采样噪声,误差就可能一轮轮积累。我们因此加入了误差控制,在 Anderson 晶格模型的测试中,结果可以与使用精确对角化杂质求解器的计算很好地吻合。 + +原先我也以为主要时间会花在训练神经网络上。实际分析计算开销后发现,对这个问题来说,更难的是把需要的物理量采样到足够高的精度。神经量子态本身能表示复杂波函数,并不意味着自洽计算就会很便宜。怎样减少高精度采样的成本,是这项工作留下的问题。 + +
From 85be2ed7a727da7424e1cc3f0d7809811eba9152 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Zhanghao=20Zhouyin=20=28=E5=91=A8=E5=AF=85=E5=BC=A0?= =?UTF-8?q?=E7=9A=93=29?= <64253517+floatingCatty@users.noreply.github.com> Date: Sun, 11 Oct 2026 20:15:35 -0400 Subject: [PATCH 13/20] Preserve existing English prose at /research/slem/ --- _pages/research-slem.md | 18 ++++-------------- 1 file changed, 4 insertions(+), 14 deletions(-) diff --git a/_pages/research-slem.md b/_pages/research-slem.md index c0fd87a..6abac2d 100644 --- a/_pages/research-slem.md +++ b/_pages/research-slem.md @@ -7,22 +7,12 @@ description: "SLEM uses local, rotation-equivariant representations to predict q description_zh: "SLEM 利用局域性和旋转等变性,从原子结构预测哈密顿量、重叠矩阵等量子算符。" paper_year: 2025 paper_url: "https://arxiv.org/abs/2407.06053" +research_locale: en +translation_url: "/zh/research/slem/" +alternate_en_url: "/research/slem/" +alternate_zh_url: "/zh/research/slem/" --- -
- -在电子结构计算里,一张能带图其实只是结果的一部分。哈密顿矩阵记录轨道的能量和耦合,重叠矩阵描述轨道之间的重叠,密度矩阵则和电子的分布有关。它们都是后面继续算物性时要用到的东西。 - -DeePTB 里我们主要学习紧束缚参数。后来做 SLEM 时,想把这件事推进到更一般的原子轨道基组,直接从原子结构预测 DFT 里的这些矩阵。 - -这里有个挺麻烦的地方:原子轨道有方向。比如 p、d 轨道,不能把它们当成原子上附着的普通数字。晶体转一个角度,描述轨道耦合的矩阵也得按对应规则变化;但仅仅换个坐标方向,又不能改变材料本身的能谱。所以我们把旋转等变性(equivariance)的要求直接写进了网络结构里。 - -我们还限制网络只使用附近一定范围内的原子环境。一般的图神经网络随着消息传递层数增加,会把越来越远的原子也带进来。SLEM 则保持严格的局域性,预测大体系时就比较容易并行。为了处理高角动量轨道,我们也用了 SO(2) 卷积来降低计算开销。 - -论文里,我们在几类二维和三维材料上测试了这些算符的预测精度和数据效率。局域性和对称性确实帮了不少忙,但严格的局域近似也有适用范围,遇到明显的长程效应仍然需要小心。另外,预测出哈密顿量不代表已经算出了电导率或光学响应。大矩阵后面的这些计算,是我后来继续研究的问题之一。 - -
-
A band diagram is only part of what comes out of an electronic-structure calculation. The Hamiltonian matrix records orbital energies and couplings, the overlap matrix describes how the orbitals overlap, and the density matrix relates to where electrons are distributed. These are the quantities we often need for later property calculations. From dbb4e41c2c5350aadca85448cc9f744420136b7f Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Zhanghao=20Zhouyin=20=28=E5=91=A8=E5=AF=85=E5=BC=A0?= =?UTF-8?q?=E7=9A=93=29?= <64253517+floatingCatty@users.noreply.github.com> Date: Sun, 11 Oct 2026 20:15:37 -0400 Subject: [PATCH 14/20] Add Chinese slem research article with reciprocal hreflang --- _pages/research-slem-zh.md | 29 +++++++++++++++++++++++++++++ 1 file changed, 29 insertions(+) create mode 100644 _pages/research-slem-zh.md diff --git a/_pages/research-slem-zh.md b/_pages/research-slem-zh.md new file mode 100644 index 0000000..8e5352b --- /dev/null +++ b/_pages/research-slem-zh.md @@ -0,0 +1,29 @@ +--- +layout: research +permalink: /zh/research/slem/ +title: "SLEM:从原子结构预测量子算符" +title_en: "SLEM: predicting quantum operators from atomic structures" +description: "SLEM uses local, rotation-equivariant representations to predict quantum-operator matrices from atomic structures." +description_zh: "SLEM 利用局域性和旋转等变性,从原子结构预测哈密顿量、重叠矩阵等量子算符。" +paper_year: 2025 +paper_url: "https://arxiv.org/abs/2407.06053" +research_locale: zh +translation_url: "/research/slem/" +translation_of: "/research/slem/" +alternate_en_url: "/research/slem/" +alternate_zh_url: "/zh/research/slem/" +--- + +
+ +在电子结构计算里,一张能带图其实只是结果的一部分。哈密顿矩阵记录轨道的能量和耦合,重叠矩阵描述轨道之间的重叠,密度矩阵则和电子的分布有关。它们都是后面继续算物性时要用到的东西。 + +DeePTB 里我们主要学习紧束缚参数。后来做 SLEM 时,想把这件事推进到更一般的原子轨道基组,直接从原子结构预测 DFT 里的这些矩阵。 + +这里有个挺麻烦的地方:原子轨道有方向。比如 p、d 轨道,不能把它们当成原子上附着的普通数字。晶体转一个角度,描述轨道耦合的矩阵也得按对应规则变化;但仅仅换个坐标方向,又不能改变材料本身的能谱。所以我们把旋转等变性(equivariance)的要求直接写进了网络结构里。 + +我们还限制网络只使用附近一定范围内的原子环境。一般的图神经网络随着消息传递层数增加,会把越来越远的原子也带进来。SLEM 则保持严格的局域性,预测大体系时就比较容易并行。为了处理高角动量轨道,我们也用了 SO(2) 卷积来降低计算开销。 + +论文里,我们在几类二维和三维材料上测试了这些算符的预测精度和数据效率。局域性和对称性确实帮了不少忙,但严格的局域近似也有适用范围,遇到明显的长程效应仍然需要小心。另外,预测出哈密顿量不代表已经算出了电导率或光学响应。大矩阵后面的这些计算,是我后来继续研究的问题之一。 + +
From 4515edeb70e3a803c5ca7d6e6e60020813580863 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Zhanghao=20Zhouyin=20=28=E5=91=A8=E5=AF=85=E5=BC=A0?= =?UTF-8?q?=E7=9A=93=29?= <64253517+floatingCatty@users.noreply.github.com> Date: Sun, 11 Oct 2026 20:15:42 -0400 Subject: [PATCH 15/20] Preserve existing English prose at /research/vbcsr/ --- _pages/research-vbcsr.md | 37 ++++--------------------------------- 1 file changed, 4 insertions(+), 33 deletions(-) diff --git a/_pages/research-vbcsr.md b/_pages/research-vbcsr.md index 77c6aef..24fe87e 100644 --- a/_pages/research-vbcsr.md +++ b/_pages/research-vbcsr.md @@ -9,41 +9,12 @@ project_kind: "software" code_url: "https://github.com/DeepElectron/vbcsr" paper_url: "https://arxiv.org/abs/2609.04397" paper_year: 2026 +research_locale: en +translation_url: "/zh/research/vbcsr/" +alternate_en_url: "/research/vbcsr/" +alternate_zh_url: "/zh/research/vbcsr/" --- -
- -做电子结构计算时,我经常要处理很大的矩阵。原子轨道只和附近的轨道有明显的耦合,所以矩阵里很多元素都是零。把这些零也存下来、算一遍,实在是浪费。稀疏矩阵格式能避免这件事,可当矩阵分布在很多台机器上,事情又会复杂不少。 - -VBCSR 是为这类计算写的一套工具。一个原子或者一个局部自由度组,往往对应一小块矩阵,不同的块大小还可能不一样。我们用可变大小的块来存储这些数据,通过 MPI 把矩阵分到不同进程,再对块内运算做优化。底层是 C++,也提供 Python 接口,日常试算法会方便一些。 - -其中一个我比较看重的操作是稀疏矩阵乘法。比如计算密度矩阵时,可能要反复对稀疏矩阵做乘法;如果每次相乘都产生大量很小的非零元素,矩阵很快就会变密,内存和计算量一起涨。VBCSR 支持按阈值过滤乘积中的小元素,让后续运算继续保持稀疏。仓库里放了密度矩阵纯化、态密度和量子动力学等例子。 - -阈值过滤并不是免费的近似。扔掉多小的元素合适,要看具体体系和希望控制的误差;MPI 通信开销也会随矩阵分布方式而改变。这些细节很琐碎,但大规模计算很多时候就是被它们卡住的。 - -### 如何检查实际性能 - -仓库已经有一套 [benchmark 生成脚本](https://github.com/DeepElectron/vbcsr/blob/main/generate_benchmark_report.py),底层调用 [`tests/benchmark_large_scale.py`](https://github.com/DeepElectron/vbcsr/blob/main/tests/benchmark_large_scale.py)。下面是一组特意缩小过的单机测试参数,便于先检查安装和计时流程,不是我们论文使用的正式性能数据。 - -```bash -git clone https://github.com/DeepElectron/vbcsr.git -cd vbcsr -# Install a working C++ / MPI / BLAS environment first (see the project README). -python -m pip install -e . -python -m pip install scipy matplotlib -python generate_benchmark_report.py \ - --output-dir bench-local --family random \ - --min-block 4 --max-block 8 --target-neighbors 8 \ - --spmv-blocks 80,160 --dense-blocks 100 --k-values 4,8 \ - --sparse-blocks 20,40 --min-seconds 0.5 --min-iterations 3 -``` - -运行后可以打开 `bench-local/benchmark_report.md` 和 `bench-local/benchmark_report_data.json`,查看 SpMV、稀疏矩阵乘多个向量和 SpMM 的计时,以及与 SciPy 的误差比较。真正要比较速度,必须同时写清 CPU、BLAS、线程数、块大小、非零结构和重复次数;小规模测试通过后,再逐步提高规模和运行时长。 - -尤其要区分两件事:**运算结果与参照一致**和**在某台机器上更快**。如果输出里的相对差异无法接受,就不该只引用加速比。阈值过滤的误差也需要单独验证。这里给的是可执行的测试入口,我没有在这次网站修改中重新跑这些性能数字。 - -
-
In electronic-structure calculations, I often have to work with very large matrices. Atomic orbitals mainly couple to nearby orbitals, so many matrix entries are zero. Storing and multiplying all those zeros would be wasteful. Sparse matrix formats avoid that, but distributing the calculation across many machines introduces other problems. From 439a615e4ff5a631deaa7a0b93817d87aacf88ad Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Zhanghao=20Zhouyin=20=28=E5=91=A8=E5=AF=85=E5=BC=A0?= =?UTF-8?q?=E7=9A=93=29?= <64253517+floatingCatty@users.noreply.github.com> Date: Sun, 11 Oct 2026 20:15:45 -0400 Subject: [PATCH 16/20] Add Chinese vbcsr research article with reciprocal hreflang --- _pages/research-vbcsr-zh.md | 50 +++++++++++++++++++++++++++++++++++++ 1 file changed, 50 insertions(+) create mode 100644 _pages/research-vbcsr-zh.md diff --git a/_pages/research-vbcsr-zh.md b/_pages/research-vbcsr-zh.md new file mode 100644 index 0000000..3367f1b --- /dev/null +++ b/_pages/research-vbcsr-zh.md @@ -0,0 +1,50 @@ +--- +layout: research +permalink: /zh/research/vbcsr/ +title: "VBCSR:给大规模稀疏矩阵写的一套工具" +title_en: "VBCSR: distributed sparse matrices for scientific computing" +description: "VBCSR is a distributed variable-block sparse matrix library for electronic-structure and scientific computing workloads, with MPI and thresholded matrix products." +description_zh: "VBCSR 是一套支持 MPI 的变块稀疏矩阵工具,面向电子结构计算和其他科学计算任务。" +project_kind: "software" +code_url: "https://github.com/DeepElectron/vbcsr" +paper_url: "https://arxiv.org/abs/2609.04397" +paper_year: 2026 +research_locale: zh +translation_url: "/research/vbcsr/" +translation_of: "/research/vbcsr/" +alternate_en_url: "/research/vbcsr/" +alternate_zh_url: "/zh/research/vbcsr/" +--- + +
+ +做电子结构计算时,我经常要处理很大的矩阵。原子轨道只和附近的轨道有明显的耦合,所以矩阵里很多元素都是零。把这些零也存下来、算一遍,实在是浪费。稀疏矩阵格式能避免这件事,可当矩阵分布在很多台机器上,事情又会复杂不少。 + +VBCSR 是为这类计算写的一套工具。一个原子或者一个局部自由度组,往往对应一小块矩阵,不同的块大小还可能不一样。我们用可变大小的块来存储这些数据,通过 MPI 把矩阵分到不同进程,再对块内运算做优化。底层是 C++,也提供 Python 接口,日常试算法会方便一些。 + +其中一个我比较看重的操作是稀疏矩阵乘法。比如计算密度矩阵时,可能要反复对稀疏矩阵做乘法;如果每次相乘都产生大量很小的非零元素,矩阵很快就会变密,内存和计算量一起涨。VBCSR 支持按阈值过滤乘积中的小元素,让后续运算继续保持稀疏。仓库里放了密度矩阵纯化、态密度和量子动力学等例子。 + +阈值过滤并不是免费的近似。扔掉多小的元素合适,要看具体体系和希望控制的误差;MPI 通信开销也会随矩阵分布方式而改变。这些细节很琐碎,但大规模计算很多时候就是被它们卡住的。 + +### 如何检查实际性能 + +仓库已经有一套 [benchmark 生成脚本](https://github.com/DeepElectron/vbcsr/blob/main/generate_benchmark_report.py),底层调用 [`tests/benchmark_large_scale.py`](https://github.com/DeepElectron/vbcsr/blob/main/tests/benchmark_large_scale.py)。下面是一组特意缩小过的单机测试参数,便于先检查安装和计时流程,不是我们论文使用的正式性能数据。 + +```bash +git clone https://github.com/DeepElectron/vbcsr.git +cd vbcsr +# Install a working C++ / MPI / BLAS environment first (see the project README). +python -m pip install -e . +python -m pip install scipy matplotlib +python generate_benchmark_report.py \ + --output-dir bench-local --family random \ + --min-block 4 --max-block 8 --target-neighbors 8 \ + --spmv-blocks 80,160 --dense-blocks 100 --k-values 4,8 \ + --sparse-blocks 20,40 --min-seconds 0.5 --min-iterations 3 +``` + +运行后可以打开 `bench-local/benchmark_report.md` 和 `bench-local/benchmark_report_data.json`,查看 SpMV、稀疏矩阵乘多个向量和 SpMM 的计时,以及与 SciPy 的误差比较。真正要比较速度,必须同时写清 CPU、BLAS、线程数、块大小、非零结构和重复次数;小规模测试通过后,再逐步提高规模和运行时长。 + +尤其要区分两件事:**运算结果与参照一致**和**在某台机器上更快**。如果输出里的相对差异无法接受,就不该只引用加速比。阈值过滤的误差也需要单独验证。这里给的是可执行的测试入口,我没有在这次网站修改中重新跑这些性能数字。 + +
From 351b372b94f4122f72f6e8d0a441500f90c249c6 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Zhanghao=20Zhouyin=20=28=E5=91=A8=E5=AF=85=E5=BC=A0?= =?UTF-8?q?=E7=9A=93=29?= <64253517+floatingCatty@users.noreply.github.com> Date: Sun, 11 Oct 2026 20:16:13 -0400 Subject: [PATCH 17/20] Make original research catalog fully English and link localized research URLs --- _pages/research.html | 77 +++++++++++++++++++++++--------------------- 1 file changed, 41 insertions(+), 36 deletions(-) diff --git a/_pages/research.html b/_pages/research.html index 8ea465f..b39538c 100644 --- a/_pages/research.html +++ b/_pages/research.html @@ -4,6 +4,10 @@ title: "Research" description: "Research papers and open-source work by Zhanghao Zhouyin on quantum materials, electronic structure, quantum transport and scientific agents." description_zh: "周寅张皓的量子材料、电子结构、量子输运和科学计算软件研究。" +research_locale: en +translation_url: "/zh/research/" +alternate_en_url: "/research/" +alternate_zh_url: "/zh/research/" --- @@ -31,74 +35,75 @@
{% include lamp-nav.html active="research" %}
-

研究Research

-

这里放一些我做过的研究和当时的想法。欢迎交流。A few projects I've worked on, and some thoughts about them. Feel free to get in touch.

+

Research

+

A few projects I've worked on, and some thoughts about them. Feel free to get in touch.

+

中文研究目录 ↗

2026 · arXiv -

魔角石墨烯里的声子凝聚Moiré phonon condensation in twisted bilayer graphene

-

把上万个原子的结构变化投影到声子模式上,发现主要由两个集体模式描述。Projecting a reconstruction of more than ten thousand atoms onto phonon modes reveals that two collective modes account for most of it.

- 阅读Read ↗ +

Moiré phonon condensation in twisted bilayer graphene

+

Projecting a reconstruction of more than ten thousand atoms onto phonon modes reveals that two collective modes account for most of it.

+ Read ↗
2026 · Physical Review B -

用神经网络求解量子嵌入里的杂质问题Neural quantum states for quantum embedding

-

用神经量子态处理强关联电子的杂质模型,也研究了自洽迭代里高精度采样的开销。Using neural quantum states for an impurity model of interacting electrons, including the cost of accurate sampling in self-consistent calculations.

- 阅读Read ↗ +

Neural quantum states for quantum embedding

+

Using neural quantum states for an impurity model of interacting electrons, including the cost of accurate sampling in self-consistent calculations.

+ Read ↗
2025 · arXiv -

DeePTB-NEGF:从学习到的哈密顿量计算量子输运DeePTB-NEGF: transport from learned Hamiltonians

-

将学习到的紧束缚哈密顿量接到 NEGF 求解器,与三种二维材料的 DFT-NEGF 结果对照。Using learned tight-binding Hamiltonians for NEGF transport, compared with first-principles results in three 2D materials.

- 阅读Read ↗ +

DeePTB-NEGF: transport from learned Hamiltonians

+

Using learned tight-binding Hamiltonians for NEGF transport, compared with first-principles results in three 2D materials.

+ Read ↗
2025 · Computational Materials Science -

FeRh 中交换关联泛函对磁性的影响Exchange-correlation functionals and magnetism in FeRh

-

比较几种 DFT 泛函对 FeRh 结构和磁性的预测,看看结果为什么会不同。Comparing how several DFT functionals describe the structure and magnetism of FeRh.

- 阅读Read ↗ +

Exchange-correlation functionals and magnetism in FeRh

+

Comparing how several DFT functionals describe the structure and magnetism of FeRh.

+ Read ↗
2025 · ICLR -

SLEM:从原子结构预测量子算符SLEM: predicting quantum operators from atomic structures

-

在网络里加入局域性和旋转对称性的约束,预测哈密顿量、重叠矩阵等量子算符。Predicting Hamiltonian, overlap and other quantum-operator matrices using a model built around locality and rotational symmetry.

- 阅读Read ↗ +

SLEM: predicting quantum operators from atomic structures

+

Predicting Hamiltonian, overlap and other quantum-operator matrices using a model built around locality and rotational symmetry.

+ Read ↗
2025 · Neurocomputing -

LogDet:用矩阵估计神经网络里的信息量LogDet: estimating information in neural networks

-

用对数行列式构造高维熵估计器,观察神经网络训练时的信息变化。A matrix-based entropy estimator for high-dimensional data, used to study information changes during neural-network training.

- 阅读Read ↗ +

LogDet: estimating information in neural networks

+

A matrix-based entropy estimator for high-dimensional data, used to study information changes during neural-network training.

+ Read ↗
- + 2024 · Nature Communications -

DeePTB:用机器学习做大规模电子结构计算DeePTB: machine learning for large-scale electronic structure

-

从第一性原理数据学习随原子环境变化的紧束缚参数,用来计算更大材料的电子结构。Learning environment-dependent tight-binding parameters from first-principles data to calculate electronic structure in larger materials.

- 阅读Read ↗ +

DeePTB: machine learning for large-scale electronic structure

+

Learning environment-dependent tight-binding parameters from first-principles data to calculate electronic structure in larger materials.

+ Read ↗
2023 · Physical Review B -

AD-NEGF:可微的量子输运计算AD-NEGF: differentiable quantum transport calculations

-

让 NEGF 的输运结果可以对模型参数求导,用于灵敏度分析和哈密顿量逆向设计。Computing derivatives of NEGF transport results with respect to model parameters for sensitivity analysis and Hamiltonian design.

- 阅读Read ↗ +

AD-NEGF: differentiable quantum transport calculations

+

Computing derivatives of NEGF transport results with respect to model parameters for sensitivity analysis and Hamiltonian design.

+ Read ↗
-

软件Software

- +

Software

+
Open source -

flower:把 Agent 跑过的计算整理清楚flower: keeping scientific-agent calculations reproducible

-

记录 Agent 的计算步骤、参数和结果,把有用的流程导出成可以重复运行的 protocol。Recording an agent's steps, parameters and results, then exporting useful runs as reproducible protocols.

- 阅读Read ↗ +

flower: keeping scientific-agent calculations reproducible

+

Recording an agent's steps, parameters and results, then exporting useful runs as reproducible protocols.

+ Read ↗
Open source -

VBCSR:给大规模稀疏矩阵写的一套工具VBCSR: distributed sparse matrices for scientific computing

-

支持 MPI 的变块稀疏矩阵与乘法工具,主要用于电子结构等科学计算。Variable-block sparse matrices and multiplication across MPI processes, mainly for electronic structure and related numerical work.

- 阅读Read ↗ +

VBCSR: distributed sparse matrices for scientific computing

+

Variable-block sparse matrices and multiplication across MPI processes, mainly for electronic structure and related numerical work.

+ Read ↗
-

完整论文见The papers are also listed under 论文Publications ↗.

+

The papers are also listed under Publications ↗.

{% include lamp-footer.html %}
From 618ab3450e3f77d236d4b9b11733da8e56632b0c Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Zhanghao=20Zhouyin=20=28=E5=91=A8=E5=AF=85=E5=BC=A0?= =?UTF-8?q?=E7=9A=93=29?= <64253517+floatingCatty@users.noreply.github.com> Date: Sun, 11 Oct 2026 20:16:17 -0400 Subject: [PATCH 18/20] Publish complete Chinese research catalog with static links --- _pages/research-zh.html | 112 ++++++++++++++++++++++++++++++++++++++++ 1 file changed, 112 insertions(+) create mode 100644 _pages/research-zh.html diff --git a/_pages/research-zh.html b/_pages/research-zh.html new file mode 100644 index 0000000..89cb691 --- /dev/null +++ b/_pages/research-zh.html @@ -0,0 +1,112 @@ +--- +layout: raw +permalink: /zh/research/ +title: "研究" +description: "周寅张皓的研究项目与科研软件:机器学习电子结构、量子材料、量子输运和科学计算 Agent。" +description_zh: "周寅张皓的量子材料、电子结构、量子输运和科学计算软件研究。" +research_locale: zh +translation_url: "/research/" +alternate_en_url: "/research/" +alternate_zh_url: "/zh/research/" +--- + + + + + +研究 | 周寅张皓 + +{% include lamp-head.html %} + + + + +
+ {% include lamp-nav.html active="research" %} +
+

研究

+

这里放一些我做过的研究和当时的想法。欢迎交流。

+

English research index ↗

+
+
+
+ + 2026 · arXiv +

魔角石墨烯里的声子凝聚

+

把上万个原子的结构变化投影到声子模式上,发现主要由两个集体模式描述。

+ 阅读 ↗ +
+ + 2026 · Physical Review B +

用神经网络求解量子嵌入里的杂质问题

+

用神经量子态处理强关联电子的杂质模型,也研究了自洽迭代里高精度采样的开销。

+ 阅读 ↗ +
+ + 2025 · arXiv +

DeePTB-NEGF:从学习到的哈密顿量计算量子输运

+

将学习到的紧束缚哈密顿量接到 NEGF 求解器,与三种二维材料的 DFT-NEGF 结果对照。

+ 阅读 ↗ +
+ + 2025 · Computational Materials Science +

FeRh 中交换关联泛函对磁性的影响

+

比较几种 DFT 泛函对 FeRh 结构和磁性的预测,看看结果为什么会不同。

+ 阅读 ↗ +
+ + 2025 · ICLR +

SLEM:从原子结构预测量子算符

+

在网络里加入局域性和旋转对称性的约束,预测哈密顿量、重叠矩阵等量子算符。

+ 阅读 ↗ +
+ + 2025 · Neurocomputing +

LogDet:用矩阵估计神经网络里的信息量

+

用对数行列式构造高维熵估计器,观察神经网络训练时的信息变化。

+ 阅读 ↗ +
+ + 2024 · Nature Communications +

DeePTB:用机器学习做大规模电子结构计算

+

从第一性原理数据学习随原子环境变化的紧束缚参数,用来计算更大材料的电子结构。

+ 阅读 ↗ +
+ + 2023 · Physical Review B +

AD-NEGF:可微的量子输运计算

+

让 NEGF 的输运结果可以对模型参数求导,用于灵敏度分析和哈密顿量逆向设计。

+ 阅读 ↗ +
+

软件

+ + Open source +

flower:把 Agent 跑过的计算整理清楚

+

记录 Agent 的计算步骤、参数和结果,把有用的流程导出成可以重复运行的 protocol。

+ 阅读 ↗ +
+ + Open source +

VBCSR:给大规模稀疏矩阵写的一套工具

+

支持 MPI 的变块稀疏矩阵与乘法工具,主要用于电子结构等科学计算。

+ 阅读 ↗ +
+
+

完整论文见 论文 ↗.

+
+ {% include lamp-footer.html %} +
+{% include lamp-embers.html %} + + From fd6b077be38c7d36976edfb62ceb9c91c16e5633 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Zhanghao=20Zhouyin=20=28=E5=91=A8=E5=AF=85=E5=BC=A0?= =?UTF-8?q?=E7=9A=93=29?= <64253517+floatingCatty@users.noreply.github.com> Date: Sun, 11 Oct 2026 20:17:32 -0400 Subject: [PATCH 19/20] Route research and About navigation to matching language pages when available --- _includes/lamp-nav.html | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/_includes/lamp-nav.html b/_includes/lamp-nav.html index 05c0154..c63d82f 100644 --- a/_includes/lamp-nav.html +++ b/_includes/lamp-nav.html @@ -6,9 +6,9 @@ From 54bd94ae8f29268edc38da1208508b209dd9d1c3 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Zhanghao=20Zhouyin=20=28=E5=91=A8=E5=AF=85=E5=BC=A0?= =?UTF-8?q?=E7=9A=93=29?= <64253517+floatingCatty@users.noreply.github.com> Date: Sun, 11 Oct 2026 20:18:41 -0400 Subject: [PATCH 20/20] Send Chinese homepage visitors directly to corresponding indexed Chinese research URLs --- index.html | 20 ++++++++++++-------- 1 file changed, 12 insertions(+), 8 deletions(-) diff --git a/index.html b/index.html index 5145898..4a0f16e 100644 --- a/index.html +++ b/index.html @@ -340,7 +340,7 @@

It began with a small flame.

A structural change involving more than ten thousand atoms turned out to be governed by just two collective motions.

Moiré Phonon Condensation in Magic-Angle Twisted Bilayer Graphene

2026preprint
- +
@@ -350,7 +350,7 @@

It began with a small flame.

A neural network can read only a material’s local structure and predict its quantum operators directly — reaching near first-principles accuracy while staying local enough to scale toward whole devices.

Learning local equivariant representations for quantum operators (SLEM)

2025ICLR 2025
- +
@@ -360,7 +360,7 @@

It began with a small flame.

Neural quantum states can act as the impurity solver inside quantum embedding, matching exact diagonalization on correlated models — where the real bottleneck is how precisely you sample observables, not training the network.

Neural-Quantum-States Impurity Solver for Quantum Embedding Problems

2026Phys. Rev. B
- +
@@ -370,7 +370,7 @@

It began with a small flame.

Trained only on first-principles eigenvalues, a deep-learning tight-binding model predicts the electronic structure of structures it has never seen — making million-atom, finite-temperature electronic simulations feasible.

Deep learning tight-binding approach for large-scale electronic simulations (DeePTB)

2024Nature Communications
- +
@@ -380,13 +380,13 @@

It began with a small flame.

Making the non-equilibrium Green’s function method differentiable end to end turns a quantum transport simulator into something you can optimize through — providing gradients for sensitivity analysis and inverse design.

Automatic differentiable nonequilibrium Green's function formalism (AD-NEGF)

2023Phys. Rev. B
- +
All publications → - 研究Research → + 研究Research →
@@ -406,13 +406,13 @@

It began with a small flame.

VBCSR

open source · early

Infrastructure for representing and computing with very large sparse operators in electronic-structure and quantum-transport work — the quiet plumbing that lets bigger simulations run on smaller machines.

github.com/DeepElectron/vbcsr ↗ - 关于 VBCSRAbout VBCSR ↗ + 关于 VBCSRAbout VBCSR ↗

And a few things still forming — AI-native DFT, scientific agents, and the beginnings of new project around.

@@ -568,6 +568,10 @@

Some questions are better explored together.

document.getElementById('lang-toggle').textContent = t.langLabel; document.documentElement.lang = lang === 'zh' ? 'zh-Hans' : 'en'; document.documentElement.setAttribute('data-lang', lang); + document.querySelectorAll('a[data-href-en][data-href-zh]').forEach(function(a) { + var destination = a.getAttribute('data-href-' + lang); + if (destination) a.setAttribute('href', destination); + }); try { localStorage.setItem('lang', lang); } catch (e) {} currentLang = lang; }