<div>Well that was easy !</div><div>Thanks a lot for your help !<br /></div><br /><div class="gmail_quote"><div dir="auto" class="gmail_attr">Le vendredi 26 janvier 2024 à 14:31:16 UTC+1, Jürg Hutter a écrit :<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Hi
<br>
<br>you cannot use k-points and/or SMEARING with TDDFT.
<br>
<br>regards
<br>JH
<br>
<br>________________________________________
<br>From: <a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> on behalf of Maxime Hodée <<a href data-email-masked rel="nofollow">maxime...@gmail.com</a>>
<br>Sent: Friday, January 26, 2024 2:24 PM
<br>To: cp2k
<br>Subject: [CP2K:19818] Periodit TD-DFT calculation failing : CPASSERT Failed
<br>
<br>Dear all,
<br>I am a new user of CP2K, and I am trying to learn by testing various calculation
<br>Right now, I want to calculate the UV-Vis absorption of a perodic anthracene crystal.
<br>In a prior calculation I successfully optimized the cell parameters, but I am having trouble doing the TD-DFT calculation.
<br>
<br>I followed the indications given on the CP2K website, however after the SCF optimization, the calculations crashes instantly with this error :
<br>
<br> *******************************************************************************
<br> *   ___                                                                       *
<br> *  /   \                                                                      *
<br> * [ABORT]                                                                     *
<br> *  \___/                             CPASSERT failed                          *
<br> *    |                                                                        *
<br> *  O/|                                                                        *
<br> * /| |                                                                        *
<br> * / \                                                qs_tddfpt2_methods.F:619 *
<br> *******************************************************************************
<br>
<br>
<br> ===== Routine Calling Stack =====
<br>
<br>            4 tddfpt
<br>            3 qs_energies_properties
<br>            2 qs_energies
<br>            1 CP2K
<br>
<br>My input file is (sorry for the long text) :
<br>
<br>&GLOBAL
<br>  PRINT_LEVEL medium
<br>  WALLTIME 259000
<br>  EXTENDED_FFT_LENGTHS .TRUE.
<br>  RUN_TYPE ENERGY
<br>  PREFERRED_DIAG_LIBRARY SL
<br>&END GLOBAL
<br>&FORCE_EVAL
<br>  METHOD Quickstep
<br>&PROPERTIES
<br> &TDDFPT
<br>  KERNEL FULL
<br>  NSTATES 5
<br>  MAX_ITER 100
<br>  CONVERGENCE [eV] 1.0e-7
<br>  RKS_TRIPLETS F
<br> &DIPOLE_MOMENTS
<br>  DIPOLE_FORM BERRY
<br> &END DIPOLE_MOMENTS
<br>&END TDDFPT
<br>&END PROPERTIES
<br>  STRESS_TENSOR ANALYTICAL
<br>  &PRINT
<br>    &STRESS_TENSOR ON
<br>    &END STRESS_TENSOR
<br>    &FORCES ON
<br>    &END FORCES
<br>  &END PRINT
<br>  &SUBSYS
<br>    &CELL
<br>      ABC 8.531 5.946 11.164
<br>      ALPHA_BETA_GAMMA 90.01 124.98 90.00
<br>      PERIODIC XYZ
<br>      SYMMETRY NONE
<br>      MULTIPLE_UNIT_CELL 1 1 1
<br>    &END CELL
<br>   &TOPOLOGY
<br>      MULTIPLE_UNIT_CELL 1 1 1
<br>      COORD_FILE_FORMAT PDB
<br>      COORD_FILE_NAME input.pdb
<br>   &GENERATE
<br>      REORDER
<br>   &END GENERATE
<br>   &END TOPOLOGY
<br>    &KIND H
<br>      POTENTIAL GTH-PBE
<br>      BASIS_SET ORB DZVP-MOLOPT-PBE-GTH-q1
<br>      BASIS_SET AUX_FIT admm-tzp-q1
<br>    &END KIND
<br>    &KIND C
<br>      POTENTIAL GTH-PBE
<br>      BASIS_SET ORB DZVP-MOLOPT-PBE-GTH-q4
<br>      BASIS_SET AUX_FIT admm-tzp-q4
<br>    &END KIND
<br>  &END SUBSYS
<br>  &DFT
<br>    POTENTIAL_FILE_NAME GTH_POTENTIALS
<br>    CHARGE 0
<br>    BASIS_SET_FILE_NAME GTH_BASIS_SETS
<br>    BASIS_SET_FILE_NAME BASIS_MOLOPT_UZH
<br>    BASIS_SET_FILE_NAME BASIS_ADMM_UZH
<br>   &MGRID
<br>      REL_CUTOFF 30
<br>      NGRIDS 4
<br>      CUTOFF 400
<br>    &END MGRID
<br>    &QS
<br>      EXTRAPOLATION_ORDER 3
<br>      EXTRAPOLATION ASPC
<br>      EPS_DEFAULT 1e-10
<br>    &END QS
<br>    &XC
<br>      &VDW_POTENTIAL
<br>        POTENTIAL_TYPE PAIR_POTENTIAL
<br>        &PAIR_POTENTIAL
<br>          REFERENCE_FUNCTIONAL PBE
<br>          SCALING 0.0
<br>          R_CUTOFF 10.0
<br>          PARAMETER_FILE_NAME ./dftd3.dat
<br>          TYPE DFTD3
<br>        &END PAIR_POTENTIAL
<br>      &END VDW_POTENTIAL
<br>      &XC_FUNCTIONAL PBE
<br>        &PBE
<br>          PARAMETRIZATION ORIG
<br>        &END PBE
<br>      &END XC_FUNCTIONAL
<br>      &XC_GRID
<br>        USE_FINER_GRID FALSE
<br>        XC_DERIV PW
<br>      &END XC_GRID
<br>    &END XC
<br>    &KPOINTS
<br>      SYMMETRY FALSE
<br>      FULL_GRID .True.
<br>      SCHEME MONKHORST-PACK 4 8 4
<br>    &END KPOINTS
<br>    &POISSON
<br>      PERIODIC XYZ
<br>      POISSON_SOLVER PERIODIC
<br>    &END POISSON
<br>    &SCF
<br>      ADDED_MOS 14
<br>      EPS_SCF 1e-06
<br>      MAX_SCF 400
<br>      SCF_GUESS ATOMIC
<br>      &DIAGONALIZATION .TRUE.
<br>        ALGORITHM STANDARD
<br>      &END DIAGONALIZATION
<br>      &PRINT
<br>        &RESTART
<br>          BACKUP_COPIES 0
<br>        &END RESTART
<br>      &END PRINT
<br>      &MIXING .TRUE.
<br>        METHOD BROYDEN_MIXING
<br>        BETA 1.25
<br>        ALPHA 0.2
<br>        NBUFFER 24
<br>      &END MIXING
<br>      &SMEAR TRUE
<br>        METHOD FERMI_DIRAC
<br>        ELECTRONIC_TEMPERATURE [K] 300
<br>      &END SMEAR
<br>    &END SCF
<br>  &END DFT
<br>&END FORCE_EVAL
<br>
<br>I would be really greatful if I could have help regarding this issue
<br>
<br>Thanks in advance for your help,
<br>Maxime Hodée
<br>
<br>
<br>--
<br>You received this message because you are subscribed to the Google Groups "cp2k" group.
<br>To unsubscribe from this group and stop receiving emails from it, send an email to <a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a><mailto:<a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>>.
<br>To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/822372b9-bba8-472e-a19b-83c06a6e0b38n%40googlegroups.com" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=fr&q=https://groups.google.com/d/msgid/cp2k/822372b9-bba8-472e-a19b-83c06a6e0b38n%2540googlegroups.com&source=gmail&ust=1706363283437000&usg=AOvVaw36wvyobCPcbcdzltGQi7B3">https://groups.google.com/d/msgid/cp2k/822372b9-bba8-472e-a19b-83c06a6e0b38n%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/822372b9-bba8-472e-a19b-83c06a6e0b38n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=fr&q=https://groups.google.com/d/msgid/cp2k/822372b9-bba8-472e-a19b-83c06a6e0b38n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1706363283437000&usg=AOvVaw3_ORBsPQc90X1DpY5FjzuB">https://groups.google.com/d/msgid/cp2k/822372b9-bba8-472e-a19b-83c06a6e0b38n%40googlegroups.com?utm_medium=email&utm_source=footer</a>>.
<br></blockquote></div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/1b2b0c6e-9832-4234-b54a-a67f81544c38n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/1b2b0c6e-9832-4234-b54a-a67f81544c38n%40googlegroups.com</a>.<br />