Size: 911
Comment:
|
Size: 2759
Comment:
|
Deletions are marked like this. | Additions are marked like this. |
Line 7: | Line 7: |
Line 9: | Line 8: |
Line 13: | Line 11: |
== Introduction == Geometry optimiser of BDF package. Two optimizers can be used: 1. Interface to DL-FIND. For details, see: http://ccpforge.cse.rl.ac.uk/gf/project/dl-find/ "DL-FIND: An Open-Source Geometry Optimizer for Atomistic Simulations", Johannes Kästner, Joanne M. Carr, Thomas W. Keal, Walter Thiel, Adrian Wander and Paul Sherwood, J. Phys. Chem. A, 2009, 113 (43), 11856-11865. DOI: 10.1021/jp9028968 2. Original optimizer in BDF developed by Dr. Yong Zhang. This BDFOPT module can be used to find minimal points, saddle points and also conical intersection points, when combined with TD-DFT and RESP modules. For details of the interface, see '''bdfopt/''', '''module/bdfopt_mod.F90,dlfind_module.F90''', '''sys_util/bdf_dlfind_util.F90''' for interface. == Current limitations == This preliminary version has several limitations: 1. Change of symmetry during optimization is not correctly handled for excited state optimization, because the correct input file needs to be prepared. 2. Conical intersection optimization does not work properly! 3. Since RESP only support nonrelativistic calculations currently, spin-free X2C based approach can not be used in optimization [to be changed in future]. == Examples == 1. [[Ground-state opt]] 2. [[ Excited state opt ]] === Conical intersection with pp-TDA === |
|
Line 14: | Line 43: |
Geometry optimiser of BDF package. This module can be used to find minimal points, saddle points and also conical intersection points. | $COMPASS Title CH2 Molecule test run, cc-pvqz Basis sto-3g # 6-31GP Geometry C 0.00000000 -0.67759997 -0.00000000 H 0.92414474 -1.21655197 0.00000000 H -0.92414474 -1.21655197 0.00000000 C -0.00000000 0.67759997 -0.00000000 H -0.92414474 1.21655197 0.00000000 H 0.92414474 1.21655197 -0.00000000 End geometry Check Skeleton $END $bdfopt solver 0 imulti 1 icoord 0 $end $XUANYUAN $END $SCF RHF charge 2 spin 1 $END $tddft imethod 4 itda 1 idiag 1 istore 1 iexit 2 crit_e 1.d-8 crit_vec 1.d-6 $end $resp geom norder 1 method 2 states 2 1 1 1 1 1 2 nfiles 1 $end $resp quad fnac double norder 1 method 2 pairs 1 1 1 1 1 1 2 nfiles 1 $end |
Line 17: | Line 130: |
== General keywords == === IPRT === |
|
Line 18: | Line 133: |
== General keywords == === Keyword1 === {{{#!wiki xxx }}} === Keyword2 === {{{#!wiki xxx }}} == Expert keywords == === Keyword3 === {{{#!wiki xxx }}} === Keyword4 === {{{#!wiki xxx }}} |
=== SOLVER === =0, DLFIND; =1, BDF optimizer |
Line 37: | Line 136: |
= Depend Files = | === MAXCYCLE === |
Line 39: | Line 138: |
|| Filename || Description || Format || || || || || |
=== TOLGRAD === |
Line 42: | Line 140: |
= Examples = | === TOLENE === |
Line 44: | Line 142: |
= Insert latex = You can insert latex formula as {{{#!latex \begin{equation} \rho_{\mu}=\sum_i w_i*\chi_{\mu}(r_i)\chi_{\mu}(r_i) \end{equation} }}} |
=== IOPT === Algorithm for optimization, default L-BFGS |
Line 52: | Line 145: |
You can also insert a latex formular inline The basis is <<latex($\chi_{\mu}$)>> |
=== ICOORD === Coordinates === IMULTI === multi-state optimization for conical intersection === ILINE === line search |
bdfopt
Contents
Introduction
Geometry optimiser of BDF package. Two optimizers can be used:
1. Interface to DL-FIND. For details, see: http://ccpforge.cse.rl.ac.uk/gf/project/dl-find/
"DL-FIND: An Open-Source Geometry Optimizer for Atomistic Simulations", Johannes Kästner, Joanne M. Carr, Thomas W. Keal, Walter Thiel, Adrian Wander and Paul Sherwood, J. Phys. Chem. A, 2009, 113 (43), 11856-11865. DOI: 10.1021/jp9028968
2. Original optimizer in BDF developed by Dr. Yong Zhang.
This BDFOPT module can be used to find minimal points, saddle points and also conical intersection points, when combined with TD-DFT and RESP modules. For details of the interface, see bdfopt/, module/bdfopt_mod.F90,dlfind_module.F90, sys_util/bdf_dlfind_util.F90 for interface.
Current limitations
This preliminary version has several limitations:
1. Change of symmetry during optimization is not correctly handled for excited state optimization, because the correct input file needs to be prepared.
2. Conical intersection optimization does not work properly!
3. Since RESP only support nonrelativistic calculations currently, spin-free X2C based approach can not be used in optimization [to be changed in future].
Examples
Conical intersection with pp-TDA
$COMPASS Title CH2 Molecule test run, cc-pvqz Basis sto-3g # 6-31GP Geometry C 0.00000000 -0.67759997 -0.00000000 H 0.92414474 -1.21655197 0.00000000 H -0.92414474 -1.21655197 0.00000000 C -0.00000000 0.67759997 -0.00000000 H -0.92414474 1.21655197 0.00000000 H 0.92414474 1.21655197 -0.00000000 End geometry Check Skeleton $END $bdfopt solver 0 imulti 1 icoord 0 $end $XUANYUAN $END $SCF RHF charge 2 spin 1 $END $tddft imethod 4 itda 1 idiag 1 istore 1 iexit 2 crit_e 1.d-8 crit_vec 1.d-6 $end $resp geom norder 1 method 2 states 2 1 1 1 1 1 2 nfiles 1 $end $resp quad fnac double norder 1 method 2 pairs 1 1 1 1 1 1 2 nfiles 1 $end
General keywords
IPRT
SOLVER
=0, DLFIND; =1, BDF optimizer
MAXCYCLE
TOLGRAD
TOLENE
IOPT
Algorithm for optimization, default L-BFGS
ICOORD
Coordinates
IMULTI
multi-state optimization for conical intersection
ILINE
line search