[CP2K-user] [CP2K:21269] TDDFT periodic with Hybrid functional

oandr...@gmail.com oandresg15 at gmail.com
Fri Mar 7 09:57:21 UTC 2025


Dear. Prof Jutter, 

Thank you for your email , 
If I may , I did a calculation without the description of XC in the 
PROPERTIES section in a system that is a singlet and I compared with the 
RKS and UKS solutions , 
However I have some differences in the first excitations energies. 
I was wondering if this difference can be as a result of the KERNEL ,  ( 
though I expected differences between RKS and UKS , I am not sure If in 
this case is correct_



 R-TDDFPT states of multiplicity 1
 Transition dipoles calculated using velocity formulation

         State    Excitation        Transition dipole (a.u.)       
 Oscillator
         number   energy (eV)       x           y           z     strength 
(a.u.)
        
 ------------------------------------------------------------------------
 TDDFPT|      1       1.28794  -4.0827E-01 -2.5202E-01 -1.7990E-01   
8.28480E-03
 TDDFPT|      2       1.33263   5.2202E-01 -8.8606E-01  3.4546E-03   
3.45302E-02
 TDDFPT|      3       1.54642  -5.2481E-01  8.9282E-01  9.5216E-04   
4.06359E-02
 TDDFPT|      4       1.56965   1.1752E+00  7.0660E-01 -3.0917E-02   
7.23448E-02
 TDDFPT|      5       1.63983  -4.6491E-01  7.8123E-01  7.6389E-04   
3.32030E-02
 TDDFPT|      6       1.74362   8.1594E-01  4.7197E-01  8.3036E-02   
3.82494E-02
 TDDFPT|      7       1.83057  -2.1331E-01  3.3635E-01  2.2744E-04   
7.11434E-03
 TDDFPT|      8       1.87648  -6.5168E-01 -3.7868E-01  2.2878E-02   
2.61403E-02
 TDDFPT|      9       1.88707   3.5067E-01  8.0384E-01 -3.9408E-03   
3.55596E-02
 TDDFPT|     10       1.88755   1.0765E+00  1.9611E-01 -4.4855E-03   
5.53681E-02


 U-TDDFPT states of multiplicity 1
 Transition dipoles calculated using velocity formulation

         State    Excitation        Transition dipole (a.u.)       
 Oscillator
         number   energy (eV)       x           y           z     strength 
(a.u.)
        
 ------------------------------------------------------------------------
 TDDFPT|      1       0.84089   6.5188E-01  1.1321E+00  8.8625E-04   
3.51567E-02
 TDDFPT|      2       0.93313   6.3627E-01  3.6431E-01  1.7204E-01   
1.29659E-02
 TDDFPT|      3       1.16569   1.0576E+00  6.2774E-01  1.3961E-01   
4.37563E-02
 TDDFPT|      4       1.27747   4.3857E-01  7.2817E-01  3.5816E-03   
2.26149E-02
 TDDFPT|      5       1.28776   2.8845E-01  1.7849E-01  1.2716E-01   
4.14021E-03
 TDDFPT|      6       1.33256   3.6916E-01  6.2625E-01  2.5210E-03   
1.72533E-02
 TDDFPT|      7       1.42774   2.1708E-01  3.7900E-01  5.2884E-03   
6.67373E-03
 TDDFPT|      8       1.44971   6.4799E-02  1.8375E-02  1.9259E-01   
1.47852E-03
 TDDFPT|      9       1.54640   3.7002E-01  6.3152E-01  6.5700E-04   
2.02970E-02
 TDDFPT|     10       1.56956   8.3101E-01  4.9878E-01  2.1858E-02   
3.61402E-02

Thank you again for the time, 

Andres Ortega-Guerrreo
nanotech at surface Empa 

El viernes, 7 de marzo de 2025 a la(s) 10:38:59 a.m. UTC+1, Jürg Hutter 
escribió:

> Hi
>
> you don't have to repeat the XC section within TDDFT if you want to
> use the same functional as for the ground state.
>
> The XC option in the TDDFT section allows to change the functional
> for the kernel only. However, for hybrids/ADMM this is restricted.
>
> regards
> JH
>
> ________________________________________
> From: cp... at googlegroups.com <cp... at googlegroups.com> on behalf of 
> oandr... at gmail.com <oandr... at gmail.com>
> Sent: Friday, March 7, 2025 7:43 AM
> To: cp2k
> Subject: [CP2K:21266] TDDFT periodic with Hybrid functional
>
> Dear Cp2k devs and users,
>
> I am trying to do a tddft calculation for a periodic system using a 
> modified hybrid functional where HF is 15%
>
> &AUXILIARY_DENSITY_MATRIX_METHOD
> ADMM_PURIFICATION_METHOD NONE
> METHOD BASIS_PROJECTION
> EXCH_SCALING_MODEL NONE
> EXCH_CORRECTION_FUNC PBEX
> &END AUXILIARY_DENSITY_MATRIX_METHOD
> &QS
> EPS_PGF_ORB 1.0000000000000001E-032
> METHOD GPW
> &END QS
> &MGRID
> NGRIDS 5
> CUTOFF 7.0000000000000000E+002
> REL_CUTOFF 7.0000000000000000E+001
> &END MGRID
> &XC
> DENSITY_CUTOFF 1.0000000000000000E-010
> GRADIENT_CUTOFF 1.0000000000000000E-010
> TAU_CUTOFF 1.0000000000000000E-010
> &XC_FUNCTIONAL NO_SHORTCUT
> &PBE T
> SCALE_X 8.4999999999999998E-001
> SCALE_C 1.0000000000000000E+000
> &END PBE
> !&PBE_HOLE_T_C_LR T
> ! SCALE_X 1.4999999999999999E-001
> ! CUTOFF_RADIUS 8.0000000000000000E+000
> !&END PBE_HOLE_T_C_LR
> &END XC_FUNCTIONAL
> &HF
> FRACTION 1.5000000000000022E-001
> &SCREENING
> EPS_SCHWARZ 9.9999999999999995E-008
> SCREEN_ON_INITIAL_P F
> &END SCREENING
> &INTERACTION_POTENTIAL
> POTENTIAL_TYPE TRUNCATED
> CUTOFF_RADIUS 8.0000000000000000E+000
> T_C_G_DATA "t_c_g.dat"
> &END INTERACTION_POTENTIAL
> &MEMORY
> EPS_STORAGE_SCALING 1.0000000000000001E-001
> MAX_MEMORY 40000
> &END MEMORY
> &END HF
> &END XC
>
>
> I was wondering If I should define the XC_FUNCTIONAL in the TDDFPT section 
> , or if like this is enough
>
> &PROPERTIES
> &TDDFPT
> KERNEL FULL
> NSTATES 10
> MAX_ITER 100
> CONVERGENCE [eV] 1.0e-7
> ADMM_KERNEL_XC_CORRECTION T
> &MGRID
> CUTOFF 700
> REL_CUTOFF 70
> &END MGRID
> &PRINT
> &NTO_ANALYSIS
> CUBE_FILES T
> STRIDE 3 3 3
> &END NTO_ANALYSIS
> &END PRINT
> &END TDDFPT
> &END PROPERTIES
>
> or should i include
>
> &XC
> &XC_FUNCTIONAL NO_SHORTCUT
> &PBE T
> SCALE_X 8.4999999999999998E-001
> SCALE_C 1.0000000000000000E+000
> &END PBE
> !&PBE_HOLE_T_C_LR T
> ! SCALE_X 1.4999999999999999E-001
> ! CUTOFF_RADIUS 8.0000000000000000E+000
> !&END PBE_HOLE_T_C_LR
> &END XC_FUNCTIONAL
> &HF
> FRACTION 1.5000000000000022E-001
> &SCREENING
> EPS_SCHWARZ 9.9999999999999995E-008
> SCREEN_ON_INITIAL_P F
> &END SCREENING
> &INTERACTION_POTENTIAL
> POTENTIAL_TYPE TRUNCATED
> CUTOFF_RADIUS 8.0000000000000000E+000
> T_C_G_DATA "t_c_g.dat"
> &END INTERACTION_POTENTIAL
> &END HF
> &END XC
>
> In the XC definition of TDDFPT ,
>
> Thank you for your suggestions and help
>
> Andres Ortega-Guerrero
> nanotech at surfaces | EMPA
>
> --
> You received this message because you are subscribed to the Google Groups 
> "cp2k" group.
> To unsubscribe from this group and stop receiving emails from it, send an 
> email to cp2k+uns... at googlegroups.com<mailto:cp2k+uns... at googlegroups.com
> >.
> To view this discussion visit 
> https://groups.google.com/d/msgid/cp2k/2057888c-c55f-46f1-877e-88d130612009n%40googlegroups.com
> <
> https://groups.google.com/d/msgid/cp2k/2057888c-c55f-46f1-877e-88d130612009n%40googlegroups.com?utm_medium=email&utm_source=footer
> >.
>

-- 
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/cp2k/50430063-731c-4d68-83d8-f590c3ecdba9n%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20250307/b6b6643a/attachment-0001.htm>


More information about the CP2K-user mailing list