welcome: please sign in
location: Diff for "bdfopt"
Differences between revisions 3 and 38 (spanning 35 versions)
Revision 3 as of 2014-11-14 02:27:51
Size: 1372
Editor: 124
Comment:
Revision 38 as of 2018-10-17 22:59:47
Size: 2333
Editor: leiyibo
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. This BDFOPT module can be used to find minimal points, saddle points and also conical intersection points. For details of the interface, see '''bdfopt/''', '''module/bdfopt_mod.F90,dlfind_module.F90''', '''sys_util/bdf_dlfind_util.F90''' for interface.

Two optimizers can be used by specifying the '''solver''' keyword:

1. [solver=1] Original optimizer in BDF developed by Dr. Yong Zhang using '''redundant internal coordinates''' (more efficient!).

2. [solver=0] Interface to the package - 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

The optimization information is stored in filename.pes1 [ground state/excited state] or filename.pes2 [conical intersection].

'''Currently, the following methods are available:'''

1. HF/MCSCF with GRAD module

2. HF/DFT/TD-DFT with RESP modules

''' 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].


== Notice ==
=== MCSCF OPT ===
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.
}}}
   For MCSCF geometry optimization, SCF module should use keyword maxiter and set maxiter = 2.
}}}
Line 19: Line 48:
=== Ground-state opt === 1. [[Ground-state opt]]
Line 21: Line 50:
{{{ 2. [[Excited state opt]]
Line 23: Line 52:
$COMPASS
Title
 CH2 Molecule test run, cc-pvqz
Basis
 3-21g
# cc-pvdz
Geometry
O 0.00000000 -0.00000000 0.00174418
H -0.78045951 0.00000000 -0.10962269
H 0.78045951 0.00000000 -0.10962269
End geometry
Check
Skeleton
$END

$bdfopt
$end

$XUANYUAN
$END

$SCF
RHF
$END

$resp
geom
norder
1
method
1
$end

}}}

=== Excited state opt ===

{{{

$COMPASS
Title
 CH2 Molecule test run, cc-pvqz
Basis
 3-21g
Geometry
 O 0.000000 0.000000 0.243991
 H 0.000000 1.178727 -0.45964
 H 0.000000 -1.178727 -0.45964
End geometry
Check
Skeleton
$END

$bdfopt
$end

$XUANYUAN
$END

$SCF
RHF
$END

$tddft
imethod
1
itda
0
idiag
1
istore
1
iexit
1
crit_e
1.d-8
crit_vec
1.d-6
$end

$resp
geom
norder
1
method
2
states
1
1 4 1
nfiles
1
$end

}}}
3. [[Conical intersection]]
Line 119: Line 55:
=== Maxcycle ===
{{{#!wiki
xxx
}}}
=== Method ===
{{{#!wiki
xxx
}}}
=== 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

bdfopt

Introduction

Geometry optimiser of BDF package. This BDFOPT module can be used to find minimal points, saddle points and also conical intersection points. For details of the interface, see bdfopt/, module/bdfopt_mod.F90,dlfind_module.F90, sys_util/bdf_dlfind_util.F90 for interface.

Two optimizers can be used by specifying the solver keyword:

1. [solver=1] Original optimizer in BDF developed by Dr. Yong Zhang using redundant internal coordinates (more efficient!).

2. [solver=0] Interface to the package - 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

The optimization information is stored in filename.pes1 [ground state/excited state] or filename.pes2 [conical intersection].

Currently, the following methods are available:

1. HF/MCSCF with GRAD module

2. HF/DFT/TD-DFT with RESP modules

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].

Notice

MCSCF OPT

   For MCSCF geometry optimization, SCF module should use keyword maxiter and set maxiter = 2.

Examples

1. Ground-state opt

2. Excited state opt

3. Conical intersection

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

bdfopt (last edited 2022-11-15 06:59:28 by wzou)