creating a core hole with MOM

Matthias matthi... at gmail.com
Mon Jun 11 18:25:47 UTC 2018


Hi,

I am struggling to figure out, how to correctly use the MOM method to 
specify that the carbon atom in methane should have a core hole, i.e. to 
remove one core electron from the system. When I use the following input, 
the program crashes with a segmentation fault. However, since I ran the 
test on the UK national supercomputer Archer using the default version of 
CP2K, I would be surprised if it were compiled incorrectly, + I have been 
able to run other calculations with no problems on the same system. The 
ground state calculation with the same basis set runs smoothly in no time 
at all.

Any help would be greatly appreciated.

Best Regards,

Matthias

P.S. I am aware that a similar question has been asked before, but since it 
didn't receive any responses I thought I should start a new one. 
https://groups.google.com/forum/#!searchin/cp2k/mom%7Csort:date/cp2k/JVKa3FqbUS0/_O-xtwgZAgAJ


&GLOBAL
  PROJECT methane
  RUN_TYPE ENERGY
  PRINT_LEVEL MEDIUM
&END GLOBAL

&FORCE_EVAL
  METHOD Quickstep
  &DFT
    &POISSON
      PERIODIC NONE
      PSOLVER  WAVELET
    &END POISSON
    &QS
      METHOD GAPW
    &END QS

    CHARGE 1
    MULTIPLICITY 2
    LSD

    &SCF
      MAX_ITER_LUMOS 10000
      EPS_SCF 1.0E-6
      SCF_GUESS ATOMIC
      MAX_SCF 60
      EPS_LUMOS  0.000001
      &OUTER_SCF
        EPS_SCF 1.0E-6
        MAX_SCF 6
      &END
      &MOM
        DEOCC_BETA 1
        PROJ_FORMULA SUM
      &END MOM
    &END SCF

    &XC
      &XC_FUNCTIONAL PBE
      &END XC_FUNCTIONAL
    &END XC
  &END DFT

  &SUBSYS
    &CELL
      ABC 10 10 10
      PERIODIC NONE
    &END CELL
    &TOPOLOGY
      &CENTER_COORDINATES
      &END
    &END
    &COORD
      C     0.00000010    -0.00000174     0.00000062
      H     0.97867541     0.00000198     0.49253941
      H     0.00648797     0.72060111    -0.82527843
      H    -0.77313050     0.28064054     0.72381230
      H    -0.21202870    -1.00124727    -0.39108047
    &END COORD
    &KIND H
     &BASIS
  6
  1  0  0  4  1
      0.754732D+02           0.427821D-02
      0.113575D+02           0.328672D-01
      0.260081D+01           0.159578D+00
      0.735503D+00           0.500000D+00
  1  0  0  1  1
      0.231761D+00           1.0000000
  1  0  0  1  1
      0.741675D-01           1.0000000
  1  1  1  1  1
      0.160000D+01           1.0000000
  1  1  1  1  1
      0.450000D+00           1.0000000
  1  2  2  1  1
      0.125000D+01           1.0000000
     &END
     POTENTIAL ALL
     &POTENTIAL
     1    0    0
     0.20000000    0
     &END
    &END KIND
    &KIND C
     &BASIS
  17
  1  0  0  7  1
              0.272316E+08           0.979800E-06
              0.108927E+07           0.441346E-04
              0.871412E+05           0.784056E-03
              0.130520E+05           0.639707E-02
              0.297061E+04           0.333869E-01
              0.841460E+03           0.140202E+00
              0.274589E+03           0.500000E+00
  1  0  0  1  1
              0.992150E+02           1.0000000
  1  0  0  1  1
              0.386940E+02           1.0000000
  1  0  0  1  1
              0.159109E+02           1.0000000
  1  0  0  1  1
              0.676677E+01           1.0000000
  1  0  0  1  1
              0.289837E+01           1.0000000
  1  0  0  1  1
              0.822020E+00           1.0000000
  1  1  1  3  1
              0.963557E+03           0.290977E-02
              0.148240E+03           0.568381E-01
              0.349117E+02           0.500000E+00
  1  1  1  1  1
              0.108936E+02           1.0000000
  1  1  1  1  1
              0.419336E+01           1.0000000
  1  1  1  1  1
              0.171343E+01           1.0000000
  1  1  1  1  1
              0.735662E+00           1.0000000
  1  0  1  1  1  1
              0.323000E+00           1.0000000        1.0000000
  1  0  1  1  1  1
              0.123000E+00           1.0000000        1.0000000
  1  2  2  1  1
              0.520000E+01           1.0000000
  1  2  2  1  1
              0.135000E+01           1.0000000
  1  2  2  1  1
              0.520000E+00           1.0000000
     &END
     POTENTIAL ALL
     &POTENTIAL
     4    2    0
     0.34883045    0
     &END
     &BS
       &BETA
         L 0
         N 1
         NEL -2
       &END BETA
     &END BS
    &END KIND
  &END SUBSYS
&END FORCE_EVAL





-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20180611/29842f7d/attachment.htm>


More information about the CP2K-user mailing list