GEO_OPT and CELL_OPT difference

174 views
Skip to first unread message

sumit agrawal

unread,
May 13, 2021, 4:04:49 AM5/13/21
to cp...@googlegroups.com
Hi everyone,

What is the difference between RUN TYPE: GEO_OPT and CELL_OPT ? Since both are doing optimization. 

Can we do GEO_OPT for crystal lattice ?


Thanks,
Sumit

Krack Matthias (PSI)

unread,
May 13, 2021, 5:00:21 AM5/13/21
to cp...@googlegroups.com

Yes. GEO_OPT relaxes only the atomic positions while the cell parameters are kept fixed. CELL_OPT relaxes both concurrently.

 

M.

 

--
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...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/CALSDoYaWDjwMjEaWUQ3tM58RD%2B78guPXzV3VCStt5_AK5DrDRg%40mail.gmail.com.

sumit agrawal

unread,
May 13, 2021, 5:11:43 AM5/13/21
to cp...@googlegroups.com
Thanks Matthias,

One more query: If i want to calculate binding energy for "A molecule inside the crystal lattice and a molecule on the crystal surface". Then I need to do both the operation (CELL_OPT first then GEO_OPT) or only GEO_OPT/CELL_OPT is enough?

Please let me know.

THanks ,
Sumit

Krack Matthias (PSI)

unread,
May 13, 2021, 5:33:18 AM5/13/21
to cp...@googlegroups.com

Hi Sumit

 

CELL_OPT includes GEO_OPT. It is up to you if you want to relax the cell of your surface or bulk system in addition to all atomic positions.

 

M.

 

sumit agrawal

unread,
May 14, 2021, 12:06:47 AM5/14/21
to cp...@googlegroups.com
Thanks Matthias for the information. 

sumit agrawal

unread,
May 14, 2021, 12:15:02 AM5/14/21
to cp...@googlegroups.com
Hii,

Now I have done cell_opt for my 2x2x2 lattice system. I prepared an argon lattice. There will be 32 argon atoms on my lattice. Now I replace one argon atom with a CO molecule (Carbon dioxide).  I again do the cell_opt for this system.

After optimization I want to do vibrational analysis for this system. My calculation stopped when SCF started.  I don't know what I am doing wrong. Can you please suggest something?

The error message is:

SCF PARAMETERS         Density guess:                                   RESTART
                        --------------------------------------------------------
                        max_scf:                                             100
                        max_scf_history:                                       0
                        max_diis:                                              4
                        --------------------------------------------------------
                        eps_scf:                                        1.00E-06
                        eps_scf_history:                                0.00E+00
                        eps_diis:                                       1.00E-01
                        eps_eigval:                                     1.00E-05
                        --------------------------------------------------------
                        level_shift [a.u.]:                                 0.00
                        --------------------------------------------------------
                        Mixing method:                           DIRECT_P_MIXING
                        --------------------------------------------------------
                        Outer loop SCF in use
                        No variables optimised in outer loop
                        eps_scf                                         1.00E-06
                        max_scf                                               15
                        No outer loop optimization
                        step_size                                       5.00E-01
"ar31-co.out" 120L, 8381C                                                            


Here I attached my input and out files.

ar31-co.inp
ar31-co.out
Reply all
Reply to author
Forward
Message has been deleted
0 new messages