Dear cp2k users,<div><br></div><div>I'm new to cp2k and NEGF calculations. Recently I found model inputs in <a href="https://github.com/schulkov/cp2k/tree/negf_samples/tests/QS/regtest-negf-fft">https://github.com/schulkov/cp2k/tree/negf_samples/tests/QS/regtest-negf-fft</a>. I'm confused by the *kpoints.inp. Why doesn't Force_EVAL 0 (the entrie system) section include kpoints section? Is it calculated on gamma point?</div><div><br></div><div>Here is the inputs</div><div><br>@SET project      au111_c6h4s2_kpoints_0.50V<br><br>@SET fermi_level  0.38253510<br><br>@SET v_bias_volt  0.50<br><br>@SET temperature  293.0<br><br>
<br><br>@SET basis_set_path BASIS_MOLOPT<br><br>@SET potential_path GTH_POTENTIALS<br><br>
<br><br>@SET xc_functional  pade<br><br>
<br><br>@SET cutoff         500<br><br>@SET added_mos      10000<br><br>@SET max_scf        300<br><br>@SET eps_scf        1e-7<br><br>
<br><br>&GLOBAL<br><br>   PROJECT              ${project}<br><br>   RUN_TYPE             negf<br><br>   PRINT_LEVEL          low<br><br>   EXTENDED_FFT_LENGTHS .true.<br><br>&END GLOBAL<br><br>
<br><br>&MULTIPLE_FORCE_EVALS<br><br>   # Order in which the input file parser should store force environments.<br><br>   #<br><br>   # 'FORCE_EVAL_ORDER a b c' means that the a-th force environment --<br><br>   # which holds the entire system -- will have the internal index '0',<br><br>   # while the b-th and c-th force environments -- which correspond to the<br><br>   # 1st and 2nd semi-infinite electrodes respectively -- will have<br><br>   # the internal indices '1' and '2'. Note that internal indices start from 0.<br><br>   #<br><br>   FORCE_EVAL_ORDER 1 2 3<br><br>
<br><br>   # Different force environments contain non-identical molecular structures.<br><br>   # By default, the input parser ignores 'subsys' sections from all but the force<br><br>   # environment with the internal index '0'.<br><br>   MULTIPLE_SUBSYS .true.<br><br>&END<br><br>
<br><br>&NEGF<br><br>   &CONTACT<br><br>      FORCE_EVAL_SECTION 1<br><br>
<br><br>      &BULK_REGION<br><br>         MOLNAME L2 L3<br><br>      &END<br><br>      &SCREENING_REGION<br><br>         MOLNAME L1<br><br>      &END<br><br>
<br><br>      # comment the next line to recompute Fermi level for this (bulk) contact from scratch<br><br>      FERMI_LEVEL  ${fermi_level}<br><br>      TEMPERATURE [K] ${temperature}<br><br>      ELECTRIC_POTENTIAL [eV] ${v_bias_volt}/2<br><br>   &END CONTACT<br><br>
<br><br>   &CONTACT<br><br>      FORCE_EVAL_SECTION 2<br><br>
<br><br>      &BULK_REGION<br><br>         MOLNAME R2 R3<br><br>      &END<br><br>      &SCREENING_REGION<br><br>         MOLNAME R1<br><br>      &END<br><br>
<br><br>      # comment the next line to recompute Fermi level for this (bulk) contact from scratch.<br><br>      # Alternatively the keyword REFINE_FERMI_LEVEL can be given, so the next line will be<br><br>      # interpreted as an initial guess<br><br>      FERMI_LEVEL  ${fermi_level}<br><br>      TEMPERATURE [K] ${temperature}<br><br>      ELECTRIC_POTENTIAL [eV] -${v_bias_volt}/2<br><br>   &END CONTACT<br><br>
<br><br>   &SCATTERING_REGION<br><br>      MOLNAME L0 S R0<br><br>   &END SCATTERING_REGION<br><br>
<br><br>   &MIXING<br><br>      ALPHA   0.05<br><br>      METHOD  broyden_mixing<br><br>      NBUFFER 8<br><br>   &END MIXING<br><br>
<br><br>   INTEGRATION_MIN_POINTS 64<br><br>   INTEGRATION_MAX_POINTS 768<br><br>
<br><br>   ENERGY_LBOUND -3.0<br><br>   EPS_DENSITY   1e-4<br><br>   EPS_SCF       1e-3<br><br>   MAX_SCF         24<br><br>
<br><br>   #V_SHIFT           0.07621247<br><br>   V_SHIFT_MAX_ITERS   30<br><br>
<br><br>   # split available MPI processes into a number of groups with 12 processes per group;<br><br>   # each group will then compute lesser Green's functions at its own set of energy points<br><br>   #NPROC_POINT 12<br><br>
<br><br>   # recompute retarded surface Green's function at every NEGF self-consistent iteration.<br><br>   # In case of NEGF calculations with applied bias, it may be worth to disable caching,<br><br>   # so adaptive integration routines will prefer points that minimised overall error<br><br>   # in best possible way to ones that have precomputed Green's functions.<br><br>   DISABLE_CACHE true<br><br>&END NEGF<br><br>
<br><br>&FORCE_EVAL<br><br>   METHOD Quickstep<br><br>
<br><br>   &DFT<br><br>      BASIS_SET_FILE_NAME ${basis_set_path}<br><br>      POTENTIAL_FILE_NAME ${potential_path}<br><br>      WFN_RESTART_FILE_NAME au111_c6h4s2_0.50V.wfn<br><br>
<br><br>      &MGRID<br><br>         CUTOFF ${cutoff}<br><br>      &END MGRID<br><br>
<br><br>      &SCF<br><br>         MAX_SCF   ${max_scf}<br><br>         EPS_SCF   ${eps_scf}<br><br>         SCF_GUESS restart<br><br>         ADDED_MOS ${added_mos}<br><br>
<br><br>         &SMEAR<br><br>            ELECTRONIC_TEMPERATURE [K] ${temperature}<br><br>            METHOD fermi_dirac<br><br>         &END SMEAR<br><br>
<br><br>         &MIXING<br><br>            ALPHA   0.3<br><br>            METHOD  broyden_mixing<br><br>            NBUFFER 8<br><br>         &END MIXING<br><br>      &END SCF<br><br>
<br><br>      &POISSON<br><br>         PERIODIC xyz<br><br>      &END POISSON<br><br>
<br><br>      &XC<br><br>         &XC_FUNCTIONAL ${xc_functional}<br><br>         &END XC_FUNCTIONAL<br><br>      &END XC<br><br>   &END DFT<br><br>
<br><br>   &SUBSYS<br><br>      &CELL<br><br>         # Lattice constant = 2.920 Angstroms<br><br>         A 8.76000000000       0.0000000000        0.0000000000<br><br>         B 4.3800000000        7.5863825040        0.0000000000<br><br>         C 0.0000000000        0.0000000000       64.5128794640<br><br>         PERIODIC xyz<br><br>      &END CELL<br><br>
<br><br>      &COORD<br><br> Au         0.0000000000        0.0000000000        0.0000000000 L3<br><br> Au         2.9200000000        0.0000000000        0.0000000000 L3<br><br> Au         5.8400000000        0.0000000000        0.0000000000 L3<br><br> Au         1.4600000000        2.5287941680        0.0000000000 L3<br><br> Au         4.3800000000        2.5287941680        0.0000000000 L3<br><br> Au         7.3000000000        2.5287941680        0.0000000000 L3<br><br> Au         2.9200000000        5.0575883360        0.0000000000 L3<br><br> Au         5.8400000000        5.0575883360        0.0000000000 L3<br><br> Au         8.7600000000        5.0575883360        0.0000000000 L3<br><br> Au         1.4600000000        0.8429312920        2.3841700720 L3<br><br> Au         4.3800000000        0.8429312920        2.3841700720 L3<br><br> Au         7.3000000000        0.8429312920        2.3841700720 L3<br><br> Au         2.9200000000        3.3717254600        2.3841700720 L3<br><br> Au         5.8400000000        3.3717254600        2.3841700720 L3<br><br> Au         8.7600000000        3.3717254600        2.3841700720 L3<br><br> Au         4.3800000000        5.9005196280        2.3841700720 L3<br><br> Au         7.3000000000        5.9005196280        2.3841700720 L3<br><br> Au        10.2200000000        5.9005196280        2.3841700720 L3<br><br> Au         0.0000000000        1.6858628760        4.7683401440 L3<br><br> Au         2.9200000000        1.6858628760        4.7683401440 L3<br><br> Au         5.8400000000        1.6858628760        4.7683401440 L3<br><br> Au         1.4600000000        4.2146570440        4.7683401440 L3<br><br> Au         4.3800000000        4.2146570440        4.7683401440 L3<br><br> Au         7.3000000000        4.2146570440        4.7683401440 L3<br><br> Au         2.9200000000        6.7434512120        4.7683401440 L3<br><br> Au         5.8400000000        6.7434512120        4.7683401440 L3<br><br> Au         8.7600000000        6.7434512120        4.7683401440 L3<br><br> Au         0.0000000000        0.0000000000        7.1525099240 L2<br><br> Au         2.9200000000        0.0000000000        7.1525099240 L2<br><br> Au         5.8400000000        0.0000000000        7.1525099240 L2<br><br> Au         1.4600000000        2.5287941680        7.1525099240 L2<br><br> Au         4.3800000000        2.5287941680        7.1525099240 L2<br><br> Au         7.3000000000        2.5287941680        7.1525099240 L2<br><br> Au         2.9200000000        5.0575883360        7.1525099240 L2<br><br> Au         5.8400000000        5.0575883360        7.1525099240 L2<br><br> Au         8.7600000000        5.0575883360        7.1525099240 L2<br><br> Au         1.4600000000        0.8429312920        9.5366799960 L2<br><br> Au         4.3800000000        0.8429312920        9.5366799960 L2<br><br> Au         7.3000000000        0.8429312920        9.5366799960 L2<br><br> Au         2.9200000000        3.3717254600        9.5366799960 L2<br><br> Au         5.8400000000        3.3717254600        9.5366799960 L2<br><br> Au         8.7600000000        3.3717254600        9.5366799960 L2<br><br> Au         4.3800000000        5.9005196280        9.5366799960 L2<br><br> Au         7.3000000000        5.9005196280        9.5366799960 L2<br><br> Au        10.2200000000        5.9005196280        9.5366799960 L2<br><br> Au         0.0000000000        1.6858628760       11.9208500680 L2<br><br> Au         2.9200000000        1.6858628760       11.9208500680 L2<br><br> Au         5.8400000000        1.6858628760       11.9208500680 L2<br><br> Au         1.4600000000        4.2146570440       11.9208500680 L2<br><br> Au         4.3800000000        4.2146570440       11.9208500680 L2<br><br> Au         7.3000000000        4.2146570440       11.9208500680 L2<br><br> Au         2.9200000000        6.7434512120       11.9208500680 L2<br><br> Au         5.8400000000        6.7434512120       11.9208500680 L2<br><br> Au         8.7600000000        6.7434512120       11.9208500680 L2<br><br> Au         0.0000000000        0.0000000000       14.3050198480 L1<br><br> Au         2.9200000000        0.0000000000       14.3050198480 L1<br><br> Au         5.8400000000        0.0000000000       14.3050198480 L1<br><br> Au         1.4600000000        2.5287941680       14.3050198480 L1<br><br> Au         4.3800000000        2.5287941680       14.3050198480 L1<br><br> Au         7.3000000000        2.5287941680       14.3050198480 L1<br><br> Au         2.9200000000        5.0575883360       14.3050198480 L1<br><br> Au         5.8400000000        5.0575883360       14.3050198480 L1<br><br> Au         8.7600000000        5.0575883360       14.3050198480 L1<br><br> Au         1.4600000000        0.8429312920       16.6891899200 L1<br><br> Au         4.3800000000        0.8429312920       16.6891899200 L1<br><br> Au         7.3000000000        0.8429312920       16.6891899200 L1<br><br> Au         2.9200000000        3.3717254600       16.6891899200 L1<br><br> Au         5.8400000000        3.3717254600       16.6891899200 L1<br><br> Au         8.7600000000        3.3717254600       16.6891899200 L1<br><br> Au         4.3800000000        5.9005196280       16.6891899200 L1<br><br> Au         7.3000000000        5.9005196280       16.6891899200 L1<br><br> Au        10.2200000000        5.9005196280       16.6891899200 L1<br><br> Au         0.0000000000        1.6858628760       19.0733599920 L1<br><br> Au         2.9200000000        1.6858628760       19.0733599920 L1<br><br> Au         5.8400000000        1.6858628760       19.0733599920 L1<br><br> Au         1.4600000000        4.2146570440       19.0733599920 L1<br><br> Au         4.3800000000        4.2146570440       19.0733599920 L1<br><br> Au         7.3000000000        4.2146570440       19.0733599920 L1<br><br> Au         2.9200000000        6.7434512120       19.0733599920 L1<br><br> Au         5.8400000000        6.7434512120       19.0733599920 L1<br><br> Au         8.7600000000        6.7434512120       19.0733599920 L1<br><br> Au         0.0000000000        0.0000000000       21.4575297720 L0<br><br> Au         2.9200000000        0.0000000000       21.4575297720 L0<br><br> Au         5.8400000000        0.0000000000       21.4575297720 L0<br><br> Au         1.4600000000        2.5287941680       21.4575297720 L0<br><br> Au         4.3800000000        2.5287941680       21.4575297720 L0<br><br> Au         7.3000000000        2.5287941680       21.4575297720 L0<br><br> Au         2.9200000000        5.0575883360       21.4575297720 L0<br><br> Au         5.8400000000        5.0575883360       21.4575297720 L0<br><br> Au         8.7600000000        5.0575883360       21.4575297720 L0<br><br> Au         1.4600000000        0.8429312920       23.8416998440 L0<br><br> Au         4.3800000000        0.8429312920       23.8416998440 L0<br><br> Au         7.3000000000        0.8429312920       23.8416998440 L0<br><br> Au         2.9200000000        3.3717254600       23.8416998440 L0<br><br> Au         5.8400000000        3.3717254600       23.8416998440 L0<br><br> Au         8.7600000000        3.3717254600       23.8416998440 L0<br><br> Au         4.3800000000        5.9005196280       23.8416998440 L0<br><br> Au         7.3000000000        5.9005196280       23.8416998440 L0<br><br> Au        10.2200000000        5.9005196280       23.8416998440 L0<br><br> Au         0.0000000000        1.6858628760       26.2258699160 L0<br><br> Au         2.9200000000        1.6858628760       26.2258699160 L0<br><br> Au         5.8400000000        1.6858628760       26.2258699160 L0<br><br> Au         1.4600000000        4.2146570440       26.2258699160 L0<br><br> Au         4.3800000000        4.2146570440       26.2258699160 L0<br><br> Au         7.3000000000        4.2146570440       26.2258699160 L0<br><br> Au         2.9200000000        6.7434512120       26.2258699160 L0<br><br> Au         5.8400000000        6.7434512120       26.2258699160 L0<br><br> Au         8.7600000000        6.7434512120       26.2258699160 L0<br><br>  S         4.3897374445        2.5067369638       34.1662712986 S<br><br>  S         4.3871636646        2.5054981547       27.9597162296 S<br><br>  C         4.4049049466        2.4402425171       32.4352453855 S<br><br>  C         4.4016468646        2.4387739728       29.6952874054 S<br><br>  C         3.1883608421        2.4306559799       31.7528215477 S<br><br>  C         3.1865224158        2.4293893824       30.3790798006 S<br><br>  C         5.6207135759        2.4311099196       31.7510402036 S<br><br>  C         5.6188095460        2.4303813928       30.3755718748 S<br><br>  H         2.2525370715        2.4222730120       32.3198101873 S<br><br>  H         2.2491024319        2.4197512514       29.8148804049 S<br><br>  H         6.5589168926        2.4243905511       32.3145243256 S<br><br>  H         6.5551277581        2.4213380461       29.8088004207 S<br><br> Au         8.7600000000        6.7434512120       35.9028394760 R0<br><br> Au         5.8400000000        6.7434512120       35.9028394760 R0<br><br> Au         2.9200000000        6.7434512120       35.9028394760 R0<br><br> Au         7.3000000000        4.2146570440       35.9028394760 R0<br><br> Au         4.3800000000        4.2146570440       35.9028394760 R0<br><br> Au         1.4600000000        4.2146570440       35.9028394760 R0<br><br> Au         5.8400000000        1.6858628760       35.9028394760 R0<br><br> Au         2.9200000000        1.6858628760       35.9028394760 R0<br><br> Au         0.0000000000        1.6858628760       35.9028394760 R0<br><br> Au        10.2200000000        5.9005196280       38.2870095480 R0<br><br> Au         7.3000000000        5.9005196280       38.2870095480 R0<br><br> Au         4.3800000000        5.9005196280       38.2870095480 R0<br><br> Au         8.7600000000        3.3717254600       38.2870095480 R0<br><br> Au         5.8400000000        3.3717254600       38.2870095480 R0<br><br> Au         2.9200000000        3.3717254600       38.2870095480 R0<br><br> Au         7.3000000000        0.8429312920       38.2870095480 R0<br><br> Au         4.3800000000        0.8429312920       38.2870095480 R0<br><br> Au         1.4600000000        0.8429312920       38.2870095480 R0<br><br> Au         8.7600000000        5.0575883360       40.6711796200 R0<br><br> Au         5.8400000000        5.0575883360       40.6711796200 R0<br><br> Au         2.9200000000        5.0575883360       40.6711796200 R0<br><br> Au         7.3000000000        2.5287941680       40.6711796200 R0<br><br> Au         4.3800000000        2.5287941680       40.6711796200 R0<br><br> Au         1.4600000000        2.5287941680       40.6711796200 R0<br><br> Au         5.8400000000        0.0000000000       40.6711796200 R0<br><br> Au         2.9200000000        0.0000000000       40.6711796200 R0<br><br> Au         0.0000000000        0.0000000000       40.6711796200 R0<br><br> Au         8.7600000000        6.7434512120       43.0553494000 R1<br><br> Au         5.8400000000        6.7434512120       43.0553494000 R1<br><br> Au         2.9200000000        6.7434512120       43.0553494000 R1<br><br> Au         7.3000000000        4.2146570440       43.0553494000 R1<br><br> Au         4.3800000000        4.2146570440       43.0553494000 R1<br><br> Au         1.4600000000        4.2146570440       43.0553494000 R1<br><br> Au         5.8400000000        1.6858628760       43.0553494000 R1<br><br> Au         2.9200000000        1.6858628760       43.0553494000 R1<br><br> Au         0.0000000000        1.6858628760       43.0553494000 R1<br><br> Au        10.2200000000        5.9005196280       45.4395194720 R1<br><br> Au         7.3000000000        5.9005196280       45.4395194720 R1<br><br> Au         4.3800000000        5.9005196280       45.4395194720 R1<br><br> Au         8.7600000000        3.3717254600       45.4395194720 R1<br><br> Au         5.8400000000        3.3717254600       45.4395194720 R1<br><br> Au         2.9200000000        3.3717254600       45.4395194720 R1<br><br> Au         7.3000000000        0.8429312920       45.4395194720 R1<br><br> Au         4.3800000000        0.8429312920       45.4395194720 R1<br><br> Au         1.4600000000        0.8429312920       45.4395194720 R1<br><br> Au         8.7600000000        5.0575883360       47.8236895440 R1<br><br> Au         5.8400000000        5.0575883360       47.8236895440 R1<br><br> Au         2.9200000000        5.0575883360       47.8236895440 R1<br><br> Au         7.3000000000        2.5287941680       47.8236895440 R1<br><br> Au         4.3800000000        2.5287941680       47.8236895440 R1<br><br> Au         1.4600000000        2.5287941680       47.8236895440 R1<br><br> Au         5.8400000000        0.0000000000       47.8236895440 R1<br><br> Au         2.9200000000        0.0000000000       47.8236895440 R1<br><br> Au         0.0000000000        0.0000000000       47.8236895440 R1<br><br> Au         8.7600000000        6.7434512120       50.2078593240 R2<br><br> Au         5.8400000000        6.7434512120       50.2078593240 R2<br><br> Au         2.9200000000        6.7434512120       50.2078593240 R2<br><br> Au         7.3000000000        4.2146570440       50.2078593240 R2<br><br> Au         4.3800000000        4.2146570440       50.2078593240 R2<br><br> Au         1.4600000000        4.2146570440       50.2078593240 R2<br><br> Au         5.8400000000        1.6858628760       50.2078593240 R2<br><br> Au         2.9200000000        1.6858628760       50.2078593240 R2<br><br> Au         0.0000000000        1.6858628760       50.2078593240 R2<br><br> Au        10.2200000000        5.9005196280       52.5920293960 R2<br><br> Au         7.3000000000        5.9005196280       52.5920293960 R2<br><br> Au         4.3800000000        5.9005196280       52.5920293960 R2<br><br> Au         8.7600000000        3.3717254600       52.5920293960 R2<br><br> Au         5.8400000000        3.3717254600       52.5920293960 R2<br><br> Au         2.9200000000        3.3717254600       52.5920293960 R2<br><br> Au         7.3000000000        0.8429312920       52.5920293960 R2<br><br> Au         4.3800000000        0.8429312920       52.5920293960 R2<br><br> Au         1.4600000000        0.8429312920       52.5920293960 R2<br><br> Au         8.7600000000        5.0575883360       54.9761994680 R2<br><br> Au         5.8400000000        5.0575883360       54.9761994680 R2<br><br> Au         2.9200000000        5.0575883360       54.9761994680 R2<br><br> Au         7.3000000000        2.5287941680       54.9761994680 R2<br><br> Au         4.3800000000        2.5287941680       54.9761994680 R2<br><br> Au         1.4600000000        2.5287941680       54.9761994680 R2<br><br> Au         5.8400000000        0.0000000000       54.9761994680 R2<br><br> Au         2.9200000000        0.0000000000       54.9761994680 R2<br><br> Au         0.0000000000        0.0000000000       54.9761994680 R2<br><br> Au         8.7600000000        6.7434512120       57.3603692480 R3<br><br> Au         5.8400000000        6.7434512120       57.3603692480 R3<br><br> Au         2.9200000000        6.7434512120       57.3603692480 R3<br><br> Au         7.3000000000        4.2146570440       57.3603692480 R3<br><br> Au         4.3800000000        4.2146570440       57.3603692480 R3<br><br> Au         1.4600000000        4.2146570440       57.3603692480 R3<br><br> Au         5.8400000000        1.6858628760       57.3603692480 R3<br><br> Au         2.9200000000        1.6858628760       57.3603692480 R3<br><br> Au         0.0000000000        1.6858628760       57.3603692480 R3<br><br> Au        10.2200000000        5.9005196280       59.7445393200 R3<br><br> Au         7.3000000000        5.9005196280       59.7445393200 R3<br><br> Au         4.3800000000        5.9005196280       59.7445393200 R3<br><br> Au         8.7600000000        3.3717254600       59.7445393200 R3<br><br> Au         5.8400000000        3.3717254600       59.7445393200 R3<br><br> Au         2.9200000000        3.3717254600       59.7445393200 R3<br><br> Au         7.3000000000        0.8429312920       59.7445393200 R3<br><br> Au         4.3800000000        0.8429312920       59.7445393200 R3<br><br> Au         1.4600000000        0.8429312920       59.7445393200 R3<br><br> Au         8.7600000000        5.0575883360       62.1287093920 R3<br><br> Au         5.8400000000        5.0575883360       62.1287093920 R3<br><br> Au         2.9200000000        5.0575883360       62.1287093920 R3<br><br> Au         7.3000000000        2.5287941680       62.1287093920 R3<br><br> Au         4.3800000000        2.5287941680       62.1287093920 R3<br><br> Au         1.4600000000        2.5287941680       62.1287093920 R3<br><br> Au         5.8400000000        0.0000000000       62.1287093920 R3<br><br> Au         2.9200000000        0.0000000000       62.1287093920 R3<br><br> Au         0.0000000000        0.0000000000       62.1287093920 R3<br><br>     &END COORD<br><br>
<br><br>     &KIND Au<br><br>        BASIS_SET SZV-MOLOPT-SR-GTH<br><br>        POTENTIAL GTH-PBE-q11<br><br>     &END KIND<br><br>     &KIND C<br><br>        BASIS_SET TZVP-MOLOPT-GTH<br><br>        POTENTIAL GTH-PADE-q4<br><br>     &END KIND<br><br>     &KIND S<br><br>        BASIS_SET TZVP-MOLOPT-GTH<br><br>        POTENTIAL GTH-PADE-q6<br><br>     &END KIND<br><br>     &KIND H<br><br>        BASIS_SET TZVP-MOLOPT-GTH<br><br>        POTENTIAL GTH-PADE-q1<br><br>     &END KIND<br><br>  &END SUBSYS<br><br>&END FORCE_EVAL<br><br>
<br><br>&FORCE_EVAL<br><br>   METHOD Quickstep<br><br>
<br><br>   &DFT<br><br>      BASIS_SET_FILE_NAME ${basis_set_path}<br><br>      POTENTIAL_FILE_NAME ${potential_path}<br><br>      WFN_RESTART_FILE_NAME au111_c6h4s2_0.50V-lead1.kp<br><br>
<br><br>      &MGRID<br><br>         CUTOFF ${cutoff}<br><br>      &END MGRID<br><br>
<br><br>      &SCF<br><br>         MAX_SCF   ${max_scf}<br><br>         EPS_SCF   ${eps_scf}<br><br>         SCF_GUESS restart<br><br>         ADDED_MOS ${added_mos}<br><br>
<br><br>         &SMEAR<br><br>            ELECTRONIC_TEMPERATURE [K] ${temperature}<br><br>            METHOD fermi_dirac<br><br>         &END SMEAR<br><br>
<br><br>         &MIXING<br><br>            ALPHA   0.3<br><br>            METHOD  broyden_mixing<br><br>            NBUFFER 8<br><br>         &END MIXING<br><br>
<br><br>         &PRINT<br><br>            &RESTART<br><br>               FILENAME lead1-RESTART<br><br>            &END RESTART<br><br>         &END PRINT<br><br>      &END SCF<br><br>
<br><br>      &POISSON<br><br>         PERIODIC xyz<br><br>      &END POISSON<br><br>
<br><br>      &XC<br><br>        &XC_FUNCTIONAL ${xc_functional}<br><br>        &END XC_FUNCTIONAL<br><br>      &END XC<br><br>
<br><br>      &KPOINTS<br><br>         SCHEME  MONKHORST-PACK  8 8 8<br><br>         SYMMETRY OFF<br><br>         EPS_GEO 1.e-8<br><br>         FULL_GRID ON<br><br>      &END KPOINTS<br><br>   &END DFT<br><br>
<br><br>   &SUBSYS<br><br>      &CELL<br><br>         # Lattice constant = 2.920 Angstroms<br><br>         A 8.76000000000       0.0000000000        0.0000000000<br><br>         B 4.3800000000        7.5863825040        0.0000000000<br><br>         C 0.0000000000        0.0000000000        7.1525099240<br><br>         PERIODIC xyz<br><br>      &END CELL<br><br>
<br><br>      &COORD<br><br> Au         0.0000000000        0.0000000000        7.1525099240 L2<br><br> Au         2.9200000000        0.0000000000        7.1525099240 L2<br><br> Au         5.8400000000        0.0000000000        7.1525099240 L2<br><br> Au         1.4600000000        2.5287941680        7.1525099240 L2<br><br> Au         4.3800000000        2.5287941680        7.1525099240 L2<br><br> Au         7.3000000000        2.5287941680        7.1525099240 L2<br><br> Au         2.9200000000        5.0575883360        7.1525099240 L2<br><br> Au         5.8400000000        5.0575883360        7.1525099240 L2<br><br> Au         8.7600000000        5.0575883360        7.1525099240 L2<br><br> Au         1.4600000000        0.8429312920        9.5366799960 L2<br><br> Au         4.3800000000        0.8429312920        9.5366799960 L2<br><br> Au         7.3000000000        0.8429312920        9.5366799960 L2<br><br> Au         2.9200000000        3.3717254600        9.5366799960 L2<br><br> Au         5.8400000000        3.3717254600        9.5366799960 L2<br><br> Au         8.7600000000        3.3717254600        9.5366799960 L2<br><br> Au         4.3800000000        5.9005196280        9.5366799960 L2<br><br> Au         7.3000000000        5.9005196280        9.5366799960 L2<br><br> Au        10.2200000000        5.9005196280        9.5366799960 L2<br><br> Au         0.0000000000        1.6858628760       11.9208500680 L2<br><br> Au         2.9200000000        1.6858628760       11.9208500680 L2<br><br> Au         5.8400000000        1.6858628760       11.9208500680 L2<br><br> Au         1.4600000000        4.2146570440       11.9208500680 L2<br><br> Au         4.3800000000        4.2146570440       11.9208500680 L2<br><br> Au         7.3000000000        4.2146570440       11.9208500680 L2<br><br> Au         2.9200000000        6.7434512120       11.9208500680 L2<br><br> Au         5.8400000000        6.7434512120       11.9208500680 L2<br><br> Au         8.7600000000        6.7434512120       11.9208500680 L2<br><br>      &END COORD<br><br>
<br><br>      &KIND Au<br><br>         BASIS_SET SZV-MOLOPT-SR-GTH<br><br>         POTENTIAL GTH-PBE-q11<br><br>      &END KIND<br><br>   &END SUBSYS<br><br>&END FORCE_EVAL<br><br>
<br><br>&FORCE_EVAL<br><br>   METHOD Quickstep<br><br>
<br><br>   &DFT<br><br>      BASIS_SET_FILE_NAME ${basis_set_path}<br><br>      POTENTIAL_FILE_NAME ${potential_path}<br><br>      WFN_RESTART_FILE_NAME au111_c6h4s2_0.50V-lead2.kp<br><br>
<br><br>      &MGRID<br><br>         CUTOFF ${cutoff}<br><br>      &END MGRID<br><br>
<br><br>      &SCF<br><br>         MAX_SCF   ${max_scf}<br><br>         EPS_SCF   ${eps_scf}<br><br>         SCF_GUESS restart<br><br>         ADDED_MOS ${added_mos}<br><br>
<br><br>         &SMEAR<br><br>            ELECTRONIC_TEMPERATURE [K] ${temperature}<br><br>            METHOD fermi_dirac<br><br>         &END SMEAR<br><br>
<br><br>         &MIXING<br><br>            ALPHA   0.3<br><br>            METHOD  broyden_mixing<br><br>            NBUFFER 8<br><br>         &END MIXING<br><br>
<br><br>         &PRINT<br><br>            &RESTART<br><br>               FILENAME lead2-RESTART<br><br>            &END RESTART<br><br>         &END PRINT<br><br>      &END SCF<br><br>
<br><br>      &POISSON<br><br>         PERIODIC xyz<br><br>      &END POISSON<br><br>
<br><br>      &XC<br><br>         &XC_FUNCTIONAL ${xc_functional}<br><br>         &END XC_FUNCTIONAL<br><br>      &END XC<br><br>
<br><br>      &KPOINTS<br><br>         SCHEME  MONKHORST-PACK  8 8 8<br><br>         SYMMETRY OFF<br><br>         EPS_GEO 1.e-8<br><br>         FULL_GRID ON<br><br>      &END KPOINTS<br><br>   &END DFT<br><br>
<br><br>   &SUBSYS<br><br>      &CELL<br><br>         # Lattice constant = 2.920 Angstroms<br><br>         A 8.76000000000       0.0000000000        0.0000000000<br><br>         B 4.3800000000        7.5863825040        0.0000000000<br><br>         C 0.0000000000        0.0000000000        7.1525099240<br><br>         PERIODIC xyz<br><br>      &END CELL<br><br>
<br><br>      &COORD<br><br> Au         8.7600000000        6.7434512120       50.2078593240 R2<br><br> Au         5.8400000000        6.7434512120       50.2078593240 R2<br><br> Au         2.9200000000        6.7434512120       50.2078593240 R2<br><br> Au         7.3000000000        4.2146570440       50.2078593240 R2<br><br> Au         4.3800000000        4.2146570440       50.2078593240 R2<br><br> Au         1.4600000000        4.2146570440       50.2078593240 R2<br><br> Au         5.8400000000        1.6858628760       50.2078593240 R2<br><br> Au         2.9200000000        1.6858628760       50.2078593240 R2<br><br> Au         0.0000000000        1.6858628760       50.2078593240 R2<br><br> Au        10.2200000000        5.9005196280       52.5920293960 R2<br><br> Au         7.3000000000        5.9005196280       52.5920293960 R2<br><br> Au         4.3800000000        5.9005196280       52.5920293960 R2<br><br> Au         8.7600000000        3.3717254600       52.5920293960 R2<br><br> Au         5.8400000000        3.3717254600       52.5920293960 R2<br><br> Au         2.9200000000        3.3717254600       52.5920293960 R2<br><br> Au         7.3000000000        0.8429312920       52.5920293960 R2<br><br> Au         4.3800000000        0.8429312920       52.5920293960 R2<br><br> Au         1.4600000000        0.8429312920       52.5920293960 R2<br><br> Au         8.7600000000        5.0575883360       54.9761994680 R2<br><br> Au         5.8400000000        5.0575883360       54.9761994680 R2<br><br> Au         2.9200000000        5.0575883360       54.9761994680 R2<br><br> Au         7.3000000000        2.5287941680       54.9761994680 R2<br><br> Au         4.3800000000        2.5287941680       54.9761994680 R2<br><br> Au         1.4600000000        2.5287941680       54.9761994680 R2<br><br> Au         5.8400000000        0.0000000000       54.9761994680 R2<br><br> Au         2.9200000000        0.0000000000       54.9761994680 R2<br><br> Au         0.0000000000        0.0000000000       54.9761994680 R2<br><br>      &END COORD<br><br>
<br><br>      &KIND Au<br><br>         BASIS_SET SZV-MOLOPT-SR-GTH<br><br>         POTENTIAL GTH-PBE-q11<br><br>      &END KIND<br><br>   &END SUBSYS<br><br>&END FORCE_EVAL  <br></div><div><br></div><div>Thank you very much!</div><div><br></div><div>Best regards,</div><div>Mengxuan</div><div><p><br></p></div><p>
</p>