Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision		 Previous revision
Last revisionBoth sides next revision
kohn-sham_random-phase_approximation [2023/02/12 17:09] – Documentation for RIRPA/URIRPA is introduced: separate subsection and two sentences in the first paragraph. chemieegortrushikohn-sham_random-phase_approximation [2023/08/13 20:11] – [RIRPA program] chemieegortrushi
Line 319: Line 319:
  
 Example input file for spin-restricted calculations for the CO molecule: Example input file for spin-restricted calculations for the CO molecule:
-<code> +<code - examples/co_rirpa.inp>
-memory,8000,m ! memory specification+
 gthresh,twoint=1d-20,energy=1d-10,orbital=1d-8 ! tighter thresholds are recommended gthresh,twoint=1d-20,energy=1d-10,orbital=1d-8 ! tighter thresholds are recommended
 gdirect ! integral-direct mode gdirect ! integral-direct mode
Line 345: Line 344:
 </code> </code>
 As well as an example of spin-unrestricted calculation for the NH<sub>2</sub> molecule: As well as an example of spin-unrestricted calculation for the NH<sub>2</sub> molecule:
-<code> +<code - examples/nh2_urirpa.inp>
-memory,8000,m ! memory specification+
 gthresh,twoint=1d-20,energy=1d-10,orbital=1d-8 ! tighter thresholds are recommended gthresh,twoint=1d-20,energy=1d-10,orbital=1d-8 ! tighter thresholds are recommended
 gdirect ! integral-direct mode gdirect ! integral-direct mode
Line 374: Line 372:
 </code> </code>
 The following options are available for the RIRPA and URIRPA programs:\\ The following options are available for the RIRPA and URIRPA programs:\\
-* **orb** record number containing the orbital coefficients, eigenvalues, etc. from the preceding DFT calculation (default: ‘2100.2’ and ‘2200.2’ for RIRPA and URIRPA, respectively)\\ +  * **orb** record number containing the orbital coefficients, eigenvalues, etc. from the preceding DFT calculation (default: ‘2100.2’ and ‘2200.2’ for RIRPA and URIRPA, respectively)\\ 
-* **dfit** logical flag to enable density fitting during the reference energy calculation (default: ’1’)\\ +  * **dfit** logical flag to enable density fitting during the reference energy calculation (default: ’1’)\\ 
-* **sigma** logical flag to enable σ-functional calculation (default: ’1’)\\ +  * **sigma** logical flag to enable σ-functional calculation (default: ’1’)\\ 
-* **sigma_param** string containing a name for the parametrization used. Choose 'PBE_S2' [6], 'PBE0_S2' [6], 'TPSS_W' 1[5], 'B3LYP_W1' [5] parameterisation in combination with a preceding DFT calculation with PBE, PBE0, TPSS or B3LYP exchange correlation functional, respectively (default: ‘PBE_S2’)\\ +  * **sigma_param** string containing a name for the parametrization used. Choose 'PBE_S2' [6], 'PBE0_S2' [6], 'TPSS_W' 1[5], 'B3LYP_W1' [5] parameterisation in combination with a preceding DFT calculation with PBE, PBE0, TPSS or B3LYP exchange correlation functional, respectively (default: ‘PBE_S2’)\\ 
-* **write_sigma** logical flag to enable writing of sigma.dat file with reference energy, frequency integration weights and σ-values (default: ’0’)\\ +  * **write_sigma** logical flag to enable writing of sigma.dat file with reference energy, frequency integration weights and σ-values (default: ’0’)\\ 
-* **thr_overlap_ri** threshold for processing RI basis according to Section IIB2 in Ref. [7] (default: ‘1d-99’)\\ +  * **thr_overlap_ri** threshold for processing RI basis according to Section IIB2 in Ref. [7] (default: ‘1d-99’)\\ 
-* **thr_fai_ri** threshold for processing RI basis according to Section IIB5 in Ref. [7] (default: ‘1d-14’)\\ +  * **thr_fai_ri** threshold for processing RI basis according to Section IIB5 in Ref. [7] (default: ‘1d-14’)\\ 
-* **thr_rpa** threshold for throwing out contributions corresponding to small eigenvalue differences during construction of the response function (default: ‘1d-6’)\\ +  * **thr_rpa** threshold for throwing out contributions corresponding to small eigenvalue differences during construction of the response function (default: ‘1d-6’)\\ 
-* **nquadint** number of logarithmically spaced intervals for frequency integration (default ‘1’)\\ +  * **nquadint** number of logarithmically spaced intervals for frequency integration (default ‘1’)\\ 
-* **nquad** number of points per interval for frequency integration (default '50')\\ +  * **nquad** number of points per interval for frequency integration (default '50')\\ 
-* **w0** caling factor for rational the function mapping the Gauss–Legendre quadrature for the interval [−1, 1] to the interval [0, ∞], see Eqs. 37-38 in Ref. [4] for details (default: ‘2.5’)\\ +  * **w0** scaling factor for the rational function mapping the Gauss–Legendre quadrature for the interval [−1, 1] to the interval [0, ∞], see Eqs. 37-38 in Ref. [4] for details (default: ‘2.5’)\\ 
-* **vc_scal** scaling factor for the Coulomb kernel, which can be used to mimic the effect of the inclusion of the exact-exchange kernel. In the special case of non-spin-polarized two-electron systems, the RPA calculation with a Coulomb kernel scaled by 1/2 is equivalent to including of the exact-exchange kernel. Implemented only in RIRPA (default: ‘1d0’)\\ +  * **vc_scal** scaling factor for the Coulomb kernel, which can be used to mimic the effect of the inclusion of the exact-exchange kernel. In the special case of non-spin-polarized two-electron systems, the RPA calculation with a Coulomb kernel scaled by 1/2 is equivalent to including of the exact-exchange kernel. Implemented only in RIRPA (default: ‘1d0’)\\ 
-* **verb** determines the level of verbosity in the output file, integer values of 0, 1, 3 provide different levels of verbosity (default ’0’)+  * **verb** determines the level of verbosity in the output file, integer values of 0, 1, 3 provide different levels of verbosity (default ’0’)