From 9465061aa059fa4fb1f40a892efcdb6c4f092985 Mon Sep 17 00:00:00 2001 From: "k.ido" <13773226+k-ido@users.noreply.github.com> Date: Fri, 2 Aug 2024 10:24:35 +0900 Subject: [PATCH 1/3] mod ls_out for lanczos --- doc/en/source/output.rst | 2 +- doc/ja/source/output.rst | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/doc/en/source/output.rst b/doc/en/source/output.rst index 743a277d..bd7e2570 100644 --- a/doc/en/source/output.rst +++ b/doc/en/source/output.rst @@ -328,7 +328,7 @@ xxx\_ls\_out\_yyy.dat ~~~~~~~~~~~~~~~~~~~~~~ This file is the outputted files for :math:`\langle H \rangle`, -:math:`\langle H^2\rangle`, and the optimized parameter :math:`\alpha` +:math:`(\langle H^2\rangle - \langle H \rangle^2)/\langle H \rangle^2`, and the optimized parameter :math:`\alpha` obtained by Power Lanczos method. This file is outputted when ``NVMCCalMode`` = 1, ``NLanczosmode`` = 1 or 2 are set in ``ModPara`` file. Here, xxx is the header indicated by ``CDataFileHead`` in diff --git a/doc/ja/source/output.rst b/doc/ja/source/output.rst index 3d31acf5..68a1102c 100644 --- a/doc/ja/source/output.rst +++ b/doc/ja/source/output.rst @@ -327,7 +327,7 @@ xxx\_ls\_out\_yyy.dat ~~~~~~~~~~~~~~~~~~~~~~ Power Lanczos法により求めた :math:`\langle H \rangle`, -:math:`\langle H^2\rangle` および最適化パラメータ :math:`\alpha` の順に出力されます。 +:math:`(\langle H^2\rangle - \langle H \rangle^2)/\langle H \rangle^2` および最適化パラメータ :math:`\alpha` の順に出力されます。 ``ModPara`` 指定ファイルで ``NVMCCalMode`` =1, ``NLanczosmode`` =1, 2に設定することで計算されます。 xxxには ``CDataFileHead`` で指定されるヘッダが、yyyには ``ModPara`` ファイルの ``NDataIdxStart``, From 92b8485eb42b0cba8ab49eebbb2dbc581fe53bbf Mon Sep 17 00:00:00 2001 From: "k.ido" <13773226+k-ido@users.noreply.github.com> Date: Tue, 6 Aug 2024 16:18:49 +0900 Subject: [PATCH 2/3] add SRinfo --- doc/en/source/output.rst | 24 +++++++++++++++++++++++- 1 file changed, 23 insertions(+), 1 deletion(-) diff --git a/doc/en/source/output.rst b/doc/en/source/output.rst index bd7e2570..9da47986 100644 --- a/doc/en/source/output.rst +++ b/doc/en/source/output.rst @@ -37,11 +37,13 @@ where both ``NDataIdxStart`` and ``NDataQtySmp`` are defined in +------------------------------------------+------------------------------------------------------------------------------------------------------------------------+ | xxx\_out\_yyy.dat | Energy and deviation. | +------------------------------------------+------------------------------------------------------------------------------------------------------------------------+ +| xxx\_SRinfo.dat | Information related to SR optimization method. | ++------------------------------------------+------------------------------------------------------------------------------------------------------------------------+ | xxx\_var\_yyy.dat | Progress information for optimizing variational parameters. | +------------------------------------------+------------------------------------------------------------------------------------------------------------------------+ | xxx\_CalcTimer.dat | Computation time for each processes. | +------------------------------------------+------------------------------------------------------------------------------------------------------------------------+ -| xxx\_time\_zzz.dat | Progress information for Monte Carlo samplings. | +| xxx\_time\_zzz.dat | Progress information for Monte Carlo samplings. | +------------------------------------------+------------------------------------------------------------------------------------------------------------------------+ | xxx\_cisajs\_yyy.dat | One body Green’s functions. | +------------------------------------------+------------------------------------------------------------------------------------------------------------------------+ @@ -148,6 +150,26 @@ where both ``NDataIdxStart`` and ``NDataQtySmp`` are defined in 2.072610167083121837e+02 2.029162857569105916e-01 ... +xxx\_SRinfo.dat.dat +~~~~~~~~~~~~~~~~~~~~~ + +This file provides the optimization information by the SR method at each step. +From the left column, the number of variational parameters as complex numbers ``Npara``, the dimension of the S metric matrix ``Msize``, the number of variational parameters that are not optimized in the input files for trial wavefunctions such as ``Orbital``, the number of variational parameters cut by the value of ``DSROptRedCut`` in the ``ModPara`` file ``diagCut``, the maximum and minimum values of the diagonal components of the S matrix ``sDiagMax`` and ``sDiagMin``, the maximum value of the change in the variational parameters ``absRmax`` and its parameter index ``imax`` are outputted. +The header specified by ``CDataFileHead`` in the ``ModPara`` file is described in xxx. +An example of the file format is as follows. +Note that the number of imaginary parameters is counted in ``optCut`` when the variational parameters are treated as real variables (``ComplexType=0``). + +:: + + #Npara Msize optCut diagCut sDiagMax sDiagMin absRmax imax + 4 4 4 0 4.17626e-02 0.00000e+00 -1.60883e-01 4 + 4 4 4 0 3.53941e-02 0.00000e+00 1.63056e-01 0 + 4 4 4 0 3.28032e-02 0.00000e+00 1.69939e-01 0 + 4 4 4 0 3.31451e-02 0.00000e+00 1.92363e-01 0 + … + + + xxx\_CalcTimer.dat ~~~~~~~~~~~~~~~~~~~ From b68692c142e5ded0377e91865ebf15c5b9956a11 Mon Sep 17 00:00:00 2001 From: "k.ido" <13773226+k-ido@users.noreply.github.com> Date: Tue, 6 Aug 2024 16:50:31 +0900 Subject: [PATCH 3/3] mod qptrans and namelist --- doc/en/source/expert.rst | 115 ++++++++++++++++++++++----------------- doc/ja/source/expert.rst | 95 +++++++++++++++++--------------- 2 files changed, 116 insertions(+), 94 deletions(-) diff --git a/doc/en/source/expert.rst b/doc/en/source/expert.rst index 96fb2d99..30789d42 100644 --- a/doc/en/source/expert.rst +++ b/doc/en/source/expert.rst @@ -243,47 +243,52 @@ User rules - When the head of line is ``#``, the line is skipped. -===================== =================================================================================== -Keywords Details for corresponding files -===================== =================================================================================== -ModPara :math:`^*` Parameters for calculation. -LocSpin :math:`^*` Configurations of the local spins for Hamiltonian. -Trans Transfer and chemical potential for Hamiltonian. -InterAll Two-body interactions for Hamiltonian. -CoulombIntra CoulombIntra interactions. -CoulombInter CoulombInter interactions. -Hund Hund couplings. -PairHop Pair hopping couplings. -Exchange Exchange couplings. -Gutzwiller Gutzwiller factors. -Jastrow Charge Jastrow factors. -DH2 2-site doublon-holon correlation factors. -DH4 4-site doublon-holon correlation factors. -GeneralRBM_PhysLayer A part of the general RBM correlation factor including variational parameters in a physical layer :math:`\alpha_{i\sigma}`. - - GeneralRBM_HiddenLayer A part of the general RBM correlation factor including variational parameters in a hidden layer :math:`\h_{k}`. - GeneralRBM_PhysHidden A part of the general RBM correlation factor including variational parameters which connect physical and hidden layers :math:`\W_{i\sigma k}`. -Orbital :math:`^*` Pair orbital factors with anti-parallel spins :math:`f_{i\uparrow j\downarrow}`. -OrbitalAntiParallel Pair orbital factors with anti-parallel spins :math:`f_{i\uparrow j\downarrow}`. -OrbitalParallel Pair orbital factors with parallel spins :math:`f_{i\sigma j\sigma}`. -OrbitalGeneral Pair orbital factors :math:`f_{i\sigma_1 j\sigma_2}`. -TransSym :math:`^*` Translational and lattice symmetry operation. -InGutzwiller Initial values of Gutzwiller factors. -InJastrow Initial values of charge Jastrow factors. -InDH2 Initial values of 2-site doublon-holon correlation factors. -InDH4 Initial values of 4-site doublon-holon correlation factors. - InGeneralRBM_PhysLayer Initial values of variational parameters of the general RBM correlation factors in a physical layer :math:`\alpha_{i\sigma}`. - InGeneralRBM_HiddenLayer Initial values of variational parameters of the general RBM correlation factors in a hidden layer :math:`\h_{k}`. - InGeneralRBM_PhysHidden Initial values of variational parameters of the general RBM correlation factors which connect physical and hidden layers :math:`\W_{i\sigma k}`. -InOrbital Initial values of pair orbital factors :math:`f_{i\uparrow j\downarrow}`. -InOrbitalAntiParallel Initial values of pair orbital factors :math:`f_{i\uparrow j\downarrow}`. -InOrbitalParallel Initial values of pair orbital factors :math:`f_{i\sigma j\sigma}`. -InOrbitalGeneral Initial values of pair orbital factors :math:`f_{i\sigma_1 j\sigma_2}`. -OneBodyG Output components for Green functions - :math:`\langle c_{i\sigma}^{\dagger}c_{j\sigma}\rangle` -TwoBodyG Output components for Correlation functions - :math:`\langle c_{i\sigma}^{\dagger}c_{j\sigma}c_{k\tau}^{\dagger}c_{l\tau}\rangle` -===================== =================================================================================== +======================== =================================================================================== +Keywords Details for corresponding files +======================== =================================================================================== +ModPara :math:`^*` Parameters for calculation. +LocSpin :math:`^*` Configurations of the local spins for Hamiltonian. +Trans Transfer and chemical potential for Hamiltonian. +InterAll Two-body interactions for Hamiltonian. +CoulombIntra CoulombIntra interactions. +CoulombInter CoulombInter interactions. +Hund Hund couplings. +PairHop Pair hopping couplings. +Exchange Exchange couplings. +Gutzwiller Gutzwiller factors. +Jastrow Charge Jastrow factors. +DH2 2-site doublon-holon correlation factors. +DH4 4-site doublon-holon correlation factors. +GeneralRBM_PhysLayer A part of the general RBM correlation factor including variational parameters + in a physical layer :math:`\alpha_{i\sigma}`. +GeneralRBM_HiddenLayer A part of the general RBM correlation factor including variational parameters + in a hidden layer :math:`h_{k}`. +GeneralRBM_PhysHidden A part of the general RBM correlation factor including variational parameters + which connect physical and hidden layers :math:`W_{i\sigma k}`. +Orbital :math:`^*` Pair orbital factors with anti-parallel spins :math:`f_{i\uparrow j\downarrow}`. +OrbitalAntiParallel Pair orbital factors with anti-parallel spins :math:`f_{i\uparrow j\downarrow}`. +OrbitalParallel Pair orbital factors with parallel spins :math:`f_{i\sigma j\sigma}`. +OrbitalGeneral Pair orbital factors :math:`f_{i\sigma_1 j\sigma_2}`. +TransSym :math:`^*` Translational and lattice symmetry operation. +InGutzwiller Initial values of Gutzwiller factors. +InJastrow Initial values of charge Jastrow factors. +InDH2 Initial values of 2-site doublon-holon correlation factors. +InDH4 Initial values of 4-site doublon-holon correlation factors. +InGeneralRBM_PhysLayer Initial values of variational parameters of the general RBM correlation factors + in a physical layer :math:`\alpha_{i\sigma}`. +InGeneralRBM_HiddenLayer Initial values of variational parameters of the general RBM correlation factors + in a hidden layer :math:`h_{k}`. +InGeneralRBM_PhysHidden Initial values of variational parameters of the general RBM correlation factors + which connect physical and hidden layers :math:`W_{i\sigma k}`. +InOrbital Initial values of pair orbital factors :math:`f_{i\uparrow j\downarrow}`. +InOrbitalAntiParallel Initial values of pair orbital factors :math:`f_{i\uparrow j\downarrow}`. +InOrbitalParallel Initial values of pair orbital factors :math:`f_{i\sigma j\sigma}`. +InOrbitalGeneral Initial values of pair orbital factors :math:`f_{i\sigma_1 j\sigma_2}`. +OneBodyG Output components for Green functions + :math:`\langle c_{i\sigma}^{\dagger}c_{j\sigma}\rangle` +TwoBodyG Output components for Correlation functions + :math:`\langle c_{i\sigma}^{\dagger}c_{j\sigma}c_{k\tau}^{\dagger}c_{l\tau}\rangle` +======================== =================================================================================== ModPara file (modpara.def) ~~~~~~~~~~~~~~~~~~~~~~~~~~ @@ -1850,7 +1855,7 @@ User rules GeneralRBM_PhysLayer file -~~~~~~~~ +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ In the general RBM correlation factors, .. math:: @@ -1982,7 +1987,7 @@ User rules GeneralRBM_HiddenLayer file -~~~~~~~~ +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ In the general RBM correlation factors, .. math:: @@ -2108,7 +2113,7 @@ User rules GeneralRBM_PhysHidden file -~~~~~~~~ +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ In the general RBM correlation factors, .. math:: @@ -2259,6 +2264,7 @@ This file sets the calculation conditions of pair orbitals |\phi_{\rm pair} \rangle = \left[\sum_{i, j=1}^{N_s} f_{ij}c_{i\uparrow}^{\dagger}c_{j\downarrow}^{\dagger} \right]^{N/2}|0 \rangle. + Site numbers :math:`i, j` and the variational parameters :math:`f_{ij}` are indicated. An example of the file format is shown as follows. @@ -2689,10 +2695,10 @@ not done. An example of the file format is shown as follows. ==================================== == TrIdx_TrWeight_and_TrIdx_i_xi == ==================================== - 0 1.000000 - 1 1.000000 - 2 1.000000 - 3 1.000000 + 0 1.000000 0.000000 + 1 1.000000 0.000000 + 2 1.000000 0.000000 + 3 1.000000 0.000000 0 0 0 (continue...) 3 12 1 @@ -2710,7 +2716,7 @@ projection patterns as :math:`N_{\rm TS}`, respectively. - Lines 3 - 5: Header -- Lines 6 - (5+ :math:`N_{\rm TS})`: [int02] [double01] +- Lines 6 - (5+ :math:`N_{\rm TS})`: [int02] [double01] [double02] - Lines (6+ :math:`N_{\rm TS}`) - (5+ :math:`(N_s+1) \times N_{\rm TS}`): [int03] [int04] [int05] [int06] @@ -2745,9 +2751,18 @@ Parameters **Type :** double-type (blank parameter not allowed) **Description :** The weight - :math:`p_{\alpha}\cos ({\boldsymbol K}\cdot {\boldsymbol R})` of the projection + Real part of :math:`p_{\alpha}\cos ({\boldsymbol K}\cdot {\boldsymbol R})` of the projection pattern :math:`(\alpha, {\boldsymbol R})`. +- [ double02 ] + + **Type :** double-type (blank parameter not allowed) + + **Description :** The weight + Imaginary part of :math:`p_{\alpha}\cos ({\boldsymbol K}\cdot {\boldsymbol R})` of the projection + pattern :math:`(\alpha, {\boldsymbol R})`. + + - [ int03 ] **Type :** int-type (blank parameter not allowed) diff --git a/doc/ja/source/expert.rst b/doc/ja/source/expert.rst index 253f0783..fdfa9b85 100644 --- a/doc/ja/source/expert.rst +++ b/doc/ja/source/expert.rst @@ -243,44 +243,44 @@ - :math:`\#` で始まる行は読み飛ばされます。 -======================== ======================================================================= - Keywords 対応するファイルの概要 -======================== ======================================================================= - ModPara :math:`^*` 計算用のパラメータを指定します。 - LocSpin :math:`^*` 局在・遍歴スピンを指定します。 - Trans 一般的な一体相互作用を指定します。 - InterAll 一般的な二体相互作用を指定します。 - CoulombIntra 内部クーロン相互作用を指定します。 - CoulombInter サイト間クーロン相互作用を指定します。 - Hund フント結合を指定します。 - PairHop ペアホッピング相互作用を指定します。 - Exchange 交換相互作用を指定します。 - Gutzwiller 最適化するGutzwiller因子を設定します。 - Jastrow 最適化する電荷Jastrow因子を指定します。 - DH2 最適化する2サイトダブロン・ホロン相関因子を指定します。 - DH4 最適化する4サイトダブロン・ホロン相関因子を指定します。 - GeneralRBM_PhysLayer 一般的なRBM相関因子のうち、最適化する物理層での変分パラメータを指定します。 - GeneralRBM_HiddenLayer 一般的なRBM相関因子のうち、最適化する隠れ層での変分パラメータを指定します。 - GeneralRBM_PhysHidden 一般的なRBM相関因子のうち、最適化する物理層と隠れ層を繋ぐ変分パラメータを指定します。 - Orbital :math:`^*` 反平行のスピンを持つペア軌道因子を指定します。 - OrbitalAntiParallel 反平行のスピンを持つペア軌道因子を指定します。 - OrbitalParallel 平行のスピンを持つペア軌道因子を指定します。 - OrbitalGeneral ペア軌道因子を指定します。 - TransSym :math:`^*` 並進・格子対称演算子を設定します。 - InGutzwiller Gutzwiller因子の初期値を設定します。 - InJastrow 電荷Jastrow因子の初期値を設定します。 - InDH2 2サイトダブロン・ホロン相関因子の初期値を設定します。 - InDH4 4サイトダブロン・ホロン相関因子の初期値を設定します。 +========================== =========================================================================================== + Keywords 対応するファイルの概要 +========================== =========================================================================================== + ModPara :math:`^*` 計算用のパラメータを指定します。 + LocSpin :math:`^*` 局在・遍歴スピンを指定します。 + Trans 一般的な一体相互作用を指定します。 + InterAll 一般的な二体相互作用を指定します。 + CoulombIntra 内部クーロン相互作用を指定します。 + CoulombInter サイト間クーロン相互作用を指定します。 + Hund フント結合を指定します。 + PairHop ペアホッピング相互作用を指定します。 + Exchange 交換相互作用を指定します。 + Gutzwiller 最適化するGutzwiller因子を設定します。 + Jastrow 最適化する電荷Jastrow因子を指定します。 + DH2 最適化する2サイトダブロン・ホロン相関因子を指定します。 + DH4 最適化する4サイトダブロン・ホロン相関因子を指定します。 + GeneralRBM_PhysLayer 一般的なRBM相関因子のうち、最適化する物理層での変分パラメータを指定します。 + GeneralRBM_HiddenLayer 一般的なRBM相関因子のうち、最適化する隠れ層での変分パラメータを指定します。 + GeneralRBM_PhysHidden 一般的なRBM相関因子のうち、最適化する物理層と隠れ層を繋ぐ変分パラメータを指定します。 + Orbital :math:`^*` 反平行のスピンを持つペア軌道因子を指定します。 + OrbitalAntiParallel 反平行のスピンを持つペア軌道因子を指定します。 + OrbitalParallel 平行のスピンを持つペア軌道因子を指定します。 + OrbitalGeneral ペア軌道因子を指定します。 + TransSym :math:`^*` 並進・格子対称演算子を設定します。 + InGutzwiller Gutzwiller因子の初期値を設定します。 + InJastrow 電荷Jastrow因子の初期値を設定します。 + InDH2 2サイトダブロン・ホロン相関因子の初期値を設定します。 + InDH4 4サイトダブロン・ホロン相関因子の初期値を設定します。 InGeneralRBM_PhysLayer 一般的なRBM相関因子のうち、最適化する物理層での変分パラメータの初期値を設定します。 InGeneralRBM_HiddenLayer 一般的なRBM相関因子のうち、最適化する隠れ層での変分パラメータの初期値を設定します。 InGeneralRBM_PhysHidden 一般的なRBM相関因子のうち、最適化する物理層と隠れ層を繋ぐ変分パラメータの初期値を設定します。 - InOrbital ペア軌道因子 :math:`f_{i\uparrow j\downarrow}` の初期値を設定します。 - InOrbitalAntiParallel ペア軌道因子 :math:`f_{i\uparrow j\downarrow}` の初期値を設定します。 - InOrbitalParallel ペア軌道因子 :math:`f_{i\sigma j\sigma}` の初期値を設定します。 - InOrbitalGeneral ペア軌道因子 :math:`f_{i\sigma j\sigma'}` の初期値を設定します。 - OneBodyG 出力する一体グリーン関数を指定します。 - TwoBodyG 出力する二体グリーン関数を指定します。 -======================== ======================================================================= + InOrbital ペア軌道因子 :math:`f_{i\uparrow j\downarrow}` の初期値を設定します。 + InOrbitalAntiParallel ペア軌道因子 :math:`f_{i\uparrow j\downarrow}` の初期値を設定します。 + InOrbitalParallel ペア軌道因子 :math:`f_{i\sigma j\sigma}` の初期値を設定します。 + InOrbitalGeneral ペア軌道因子 :math:`f_{i\sigma j\sigma'}` の初期値を設定します。 + OneBodyG 出力する一体グリーン関数を指定します。 + TwoBodyG 出力する二体グリーン関数を指定します。 +========================== =========================================================================================== ModParaファイル (modpara.def) ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ @@ -582,7 +582,7 @@ ModParaファイル (modpara.def) **形式 :** int型 (デフォルト値=0) - **説明 :** RBMの隠れ層にあるニューロン数 :math:`N_{\rm General RBM}`を指定する整数。 + **説明 :** RBMの隠れ層にあるニューロン数 :math:`N_{\rm General RBM}` を指定する整数。 LocSpin指定ファイル(locspn.def) ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ @@ -1908,7 +1908,7 @@ RBM因子 GeneralRBM_HiddenLayer指定ファイル -~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ RBM因子 @@ -2026,7 +2026,7 @@ RBM因子 - [ int02 ]-[ int06 ] を指定する際、範囲外の整数を指定した場合はエラー終了します。 GeneralRBM_PhysHidden指定ファイル -~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ RBM因子 @@ -2542,10 +2542,10 @@ TransSym指定ファイル(qptransidx.def) ==================================== == TrIdx_TrWeight_and_TrIdx_i_xi == ==================================== - 0 1.000000 - 1 1.000000 - 2 1.000000 - 3 1.000000 + 0 1.000000 0.000000 + 1 1.000000 0.000000 + 2 1.000000 0.000000 + 3 1.000000 0.000000 0 0 0 (continue...) 3 12 1 @@ -2562,7 +2562,7 @@ TransSym指定ファイル(qptransidx.def) - 3-5行: ヘッダ(何が書かれても問題ありません)。 -- 6 - (5+ :math:`N_{\rm TS}`)行: [int02] [double01] +- 6 - (5+ :math:`N_{\rm TS}`)行: [int02] [double01] [double02] - (6+ :math:`N_{\rm TS}`) - (5+ :math:`(N_s+1) \times N_{\rm TS}`)行：[int03] [int04] [int05] [int06] @@ -2595,7 +2595,14 @@ TransSym指定ファイル(qptransidx.def) **形式 :** double型 (空白不可) **説明 :** - 射影パターン :math:`(\alpha, {\boldsymbol R})` の重み :math:`p_{\alpha}\cos ({\boldsymbol K}\cdot {\boldsymbol R})` を指定します。 + 射影パターン :math:`(\alpha, {\boldsymbol R})` の重み :math:`p_{\alpha}\cos ({\boldsymbol K}\cdot {\boldsymbol R})` の実部を指定します。 + +- [ double02 ] + + **形式 :** double型 + + **説明 :** + 射影パターン :math:`(\alpha, {\boldsymbol R})` の重み :math:`p_{\alpha}\exp (i{\boldsymbol K}\cdot {\boldsymbol R})` の虚部を指定します。 - [ int03 ]