User Tools

Site Tools


exercises:2015_ethz_mmm:ls_scf

Differences

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

Link to this comparison view

Both sides previous revision Previous revision
Last revision Both sides next revision
exercises:2015_ethz_mmm:ls_scf [2015/03/05 13:33]
vondele
exercises:2015_ethz_mmm:ls_scf [2015/03/05 13:42]
vondele [Input File]
Line 37: Line 37:
 ===== Input File ===== ===== Input File =====
 <code - 2d_polymer.inp>​ <code - 2d_polymer.inp>​
 +! Number of unit cell replications,​ increase this to create a large polymer.
 +@SET NREP 1
 +
 &GLOBAL &GLOBAL
   PROJECT 2d_polymer   PROJECT 2d_polymer
Line 44: Line 47:
 &​FORCE_EVAL &​FORCE_EVAL
   &DFT   &DFT
 +    ! linear scaling SCF
     &LS_SCF     &LS_SCF
-      ​EPS_FILTER 1E-5  ​! threshold used to determine sparsity +      ​! TRS4, does not need an estimate for the chemical potential 
-      ​MU         -0.15 ! chemical potential, an energy from within in the HOMO-LUMO gap+      PURIFICATION_METHOD TRS4 
 +      ​! threshold used to determine sparsity ​and thus speed and accuracy 
 +      ​EPS_FILTER 1E-5   
 +      ! convergence for the SCF 
 +      EPS_SCF ​   1E-5 
 +      ​! chemical potential, an energy from within in the HOMO-LUMO gap 
 +      MU         -0.15  
 +      S_PRECONDITIONER ATOMIC 
 +    &END 
 + 
 +    ! traditional SCF, diagonalization by default 
 +    &SCF
     &END     &END
  
     &QS     &QS
-      LS_SCF ​ ! Turns on linear scaling DFT. Comment this line to run with normal ​SCF.+      LS_SCF ​ ! Turns on linear scaling DFT. Comment this line to run with standard ​SCF.
  
-! ============================================================================================= +      ​! ============================================================================================= 
-! === Description of Hamiltonian (Density Functional Tight Binding + Particle Meshed Ewald) ===+      ! === Description of Hamiltonian (Density Functional Tight Binding + Smooth ​Particle Meshed Ewald) ===
       METHOD DFTB       METHOD DFTB
       &DFTB       &DFTB
Line 61: Line 76:
         DISPERSION ​        T         DISPERSION ​        T
         &​PARAMETER         &​PARAMETER
-          PARAM_FILE_PATH ​ ./dftb_params/scc+          PARAM_FILE_PATH ​ ./DFTB/scc
           PARAM_FILE_NAME ​ scc_parameter           PARAM_FILE_NAME ​ scc_parameter
           UFF_FORCE_FIELD ​ uff_table           UFF_FORCE_FIELD ​ uff_table
Line 67: Line 82:
       &END DFTB       &END DFTB
     &END QS     &END QS
 +
     &​POISSON     &​POISSON
       &EWALD       &EWALD
        ​EWALD_TYPE SPME        ​EWALD_TYPE SPME
-       ​GMAX ​160+       ! scale the number of points with the system size, roughly 1pt/A. 
 +       GMAX 40 40*${NREP} 40
        ​O_SPLINE 5        ​O_SPLINE 5
       &END EWALD       &END EWALD
     &END POISSON     &END POISSON
-! ==============================End of Hamiltonian Description ================================ +   ! ==============================End of Hamiltonian Description ================================ 
-! =============================================================================================+   ​! =============================================================================================
   &END DFT   &END DFT
  
-! Number of unit cell replications,​ increase this to create a large polymer. 
-@SET NREP 1 
   ​   ​
   &SUBSYS   &SUBSYS
     &CELL     &CELL
 +      ! a 2D system, z-direction is non-periodic
       ABC 23.577999 40.8383 30       ABC 23.577999 40.8383 30
       MULTIPLE_UNIT_CELL 1 ${NREP} 1       MULTIPLE_UNIT_CELL 1 ${NREP} 1
Line 94: Line 110:
   &END SUBSYS   &END SUBSYS
 &END FORCE_EVAL &END FORCE_EVAL
 +
 </​code>​ </​code>​
  
exercises/2015_ethz_mmm/ls_scf.txt · Last modified: 2015/03/05 13:54 by vondele