how to reduce the memory usage for HFX
Matt W
mattwa... at gmail.com
Thu Sep 28 03:07:59 UTC 2017
- SCREEN_ON_INITIAL_P
<https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT/XC/HF/SCREENING.html#list_SCREEN_ON_INITIAL_P>
should be the most important. You need to read in a converged GGA
calculation.
Also EPS_PGF_ORB doesn't need to be so small.
On Wednesday, September 27, 2017 at 10:56:06 PM UTC+8, Xiaoming Wang wrote:
>
> Hi all,
>
> Is there any way to reduce the memory usage for HFX calculations? I am
> dealing with a large scale HFX calculation with the box size of 47*47*47 A
> and 22560 atoms. After testing, I find that 30Gb mem is required for each
> MPI task. I am using the MPI/Openmp version. Following is my input, any
> suggestion on the improvement of the parameters is appreciated.
>
> &GLOBAL
> PROJECT_NAME mysys
> RUN_TYPE GEO_OPT
> PRINT_LEVEL MEDIUM
> &END GLOBAL
> &FORCE_EVAL
> METHOD QS
> &DFT
> CHARGE -1
> LSD
> BASIS_SET_FILE_NAME BASIS_MOLOPT
> BASIS_SET_FILE_NAME BASIS_ADMM_MOLOPT
> POTENTIAL_FILE_NAME GTH_POTENTIALS
> WFN_RESTART_FILE_NAME pbe.wfn
> &QS
> EPS_PGF_ORB 1.0e-32
> &END
> &MGRID
> CUTOFF 250
> REL_CUTOFF 50
> &END MGRID
> &XC
> &XC_FUNCTIONAL
> &PBE
> SCALE_X 0.75
> SCALE_C 1.0
> &END PBE
> &PBE_HOLE_T_C_LR
> CUTOFF_RADIUS 10.0
> SCALE_X 0.25
> &END PBE_HOLE_T_C_LR
> &END XC_FUNCTIONAL
> &HF
> FRACTION 0.25
> &SCREENING
> EPS_SCHWARZ 1.0e-6
> &END SCREENING
> &INTERACTION_POTENTIAL
> POTENTIAL_TYPE TRUNCATED
> CUTOFF_RADIUS 10.0
> T_C_G_DATA ./t_c_g.dat
> &END INTERACTION_POTENTIAL
> &MEMORY
> MAX_MEMORY 32000
> EPS_STORAGE_SCALING 0.1
> &END MEMORY
> &END HF
> &END XC
> &SCF
> MAX_SCF 100
> EPS_SCF 1.0e-6
> CHOLESKY INVERSE
> SCF_GUESS RESTART
> &OT
> ROTATION
> PRECONDITIONER FULL_KINETIC
> ENERGY_GAP 0.001
> &END OT
> &OUTER_SCF
> EPS_SCF 1.0e-6
> MAX_SCF 10
> &END OUTER_SCF
> &END SCF
> &PRINT
> &MULLIKEN OFF
> &END MULLIKEN
> &MO_CUBES
> WRITE_CUBE F
> NHOMO 2
> &END MO_CUBES
> &END PRINT
> &AUXILIARY_DENSITY_MATRIX_METHOD
> METHOD BASIS_PROJECTION
> ADMM_PURIFICATION_METHOD MO_DIAG
> &END AUXILIARY_DENSITY_MATRIX_METHOD
> &END DFT
> &SUBSYS
> &CELL
> ABC [angstrom] 47.216 47.216 47.216
> ALPHA_BETA_GAMMA [deg] 90 90 90
> PERIODIC XYZ
> SYMMETRY CUBIC
> &END CELL
> &TOPOLOGY
> COORD_FILE_FORMAT XYZ
> COORD_FILE_NAME el.xyz
> &END TOPOLOGY
> &KIND Cs
> ELEMENT Cs
> BASIS_SET DZVP-MOLOPT-SR-GTH
> AUX_FIT_BASIS_SET cFIT7
> POTENTIAL GTH-PBE-q9
> &END KIND
> &KIND Pb
> ELEMENT Pb
> BASIS_SET DZVP-MOLOPT-SR-GTH
> AUX_FIT_BASIS_SET cFIT6
> POTENTIAL GTH-PBE-q4
> &END KIND
> &KIND Br
> ELEMENT Br
> BASIS_SET DZVP-MOLOPT-SR-GTH
> AUX_FIT_BASIS_SET cFIT6
> POTENTIAL GTH-PBE-q7
> &END KIND
> &END SUBSYS
> &END FORCE_EVAL
> &MOTION
> &GEO_OPT
> MAX_ITER 500
> OPTIMIZER LBFGS
> &END GEO_OPT
> &END MOTION
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20170927/f48cc4ae/attachment.htm>
More information about the CP2K-user
mailing list