|
6 | 6 |
|
7 | 7 | - Geometry optimization for arbitrary state local minimum, MECI, MEP with Scipy library |
8 | 8 | - Geometry optimization for arbitrary state local minimum, MECI, TS with DL-FIND library |
| 9 | +- Geometry optimization with geomeTRIC for state-specific optimization, MECI, MECP, TS, constrained optimizations, and IRC paths |
| 10 | +- Native PySCF-backed initial guesses: `pyscf`, `sad`, and `sap` |
9 | 11 | - Energy gradient calculations interface for nonadiabatic molecular dynamics |
10 | 12 |
|
11 | 13 | ## Prerequisite |
@@ -81,6 +83,18 @@ Results, including log files and `test_report.txt`, will be stored in the curren |
81 | 83 | save_mol=False |
82 | 84 | continue_geom=False |
83 | 85 |
|
| 86 | +[geometric] |
| 87 | + coordsys=tric |
| 88 | + trust=0.1 |
| 89 | + tmax=0.3 |
| 90 | + convergence_set=GAU |
| 91 | + prefix=geometric |
| 92 | + hessian=never |
| 93 | + irc_direction=forward |
| 94 | + constraints_file='' |
| 95 | + enforce=0.0 |
| 96 | + conmethod=0 |
| 97 | + |
84 | 98 | [scf] |
85 | 99 | type=rhf |
86 | 100 | maxit=30 |
@@ -188,6 +202,7 @@ input section handle the basic information of molecular system |
188 | 202 | meci minimum energy conical intersection optimization |
189 | 203 | mep minimum energy path calculation |
190 | 204 | ts transition state optimization |
| 205 | + irc intrinsic reaction coordinate with [optimize]lib=geometric |
191 | 206 | neb nudge elasted band calculation (not avaialble yet) |
192 | 207 |
|
193 | 208 | - system // specify molecular structure or xyz file |
@@ -217,6 +232,14 @@ guess section handle the guess orbitals |
217 | 232 | hcore hcore guess |
218 | 233 | model read orbital from molden |
219 | 234 | json load data from json |
| 235 | + auto load json if the requested file exists; otherwise use huckel |
| 236 | + pyscf run a PySCF SCF calculation and import the converged orbitals |
| 237 | + sad build a PySCF superposition-of-atomic-density initial guess |
| 238 | + sap build a PySCF superposition-of-atomic-potential initial guess |
| 239 | + |
| 240 | + pyscf, sad, and sap use the native OpenQP JSON restart exporter. |
| 241 | + PySCF is installed with OpenQP's Python dependencies. MOKIT is not |
| 242 | + required for these initial guesses. |
220 | 243 |
|
221 | 244 | - file // set the guess orbital or data file |
222 | 245 |
|
@@ -384,6 +407,7 @@ optimize section handle the geometry optimization |
384 | 407 |
|
385 | 408 | scipy use scipy.optimize library (default) |
386 | 409 | dlfind use DL-FIND library |
| 410 | + geometric use geomeTRIC. Supports runtype=optimize, meci, mecp, ts, and irc |
387 | 411 |
|
388 | 412 | - optimizer // choose the scipy optimizer |
389 | 413 |
|
@@ -463,6 +487,52 @@ optimize section handle the geometry optimization |
463 | 487 | True do initial SCF iterations in every optimization step |
464 | 488 | False do not do initial SCF iterations after the first optimization step |
465 | 489 |
|
| 490 | +### [geometric] |
| 491 | + |
| 492 | +geometric section controls the geomeTRIC optimizer backend when [optimize]lib=geometric |
| 493 | + |
| 494 | +- coordsys // coordinate system passed to geomeTRIC |
| 495 | + |
| 496 | + tric translation-rotation internal coordinates (default) |
| 497 | + |
| 498 | +- trust // initial trust radius in Angstrom |
| 499 | + |
| 500 | + 0.1 (default) |
| 501 | + |
| 502 | +- tmax // maximum trust radius in Angstrom |
| 503 | + |
| 504 | + 0.3 (default) |
| 505 | + |
| 506 | +- convergence_set // geomeTRIC convergence preset |
| 507 | + |
| 508 | + GAU (default) |
| 509 | + |
| 510 | +- prefix // prefix for geomeTRIC output files |
| 511 | + |
| 512 | + geometric (default) |
| 513 | + |
| 514 | +- hessian // initial Hessian option passed to geomeTRIC |
| 515 | + |
| 516 | + never (default for minima and constrained optimization) |
| 517 | + first (recommended/automatically used for ts and irc if hessian=never) |
| 518 | + |
| 519 | +- irc_direction // IRC branch direction for runtype=irc |
| 520 | + |
| 521 | + forward (default) |
| 522 | + backward |
| 523 | + |
| 524 | +- constraints_file // geomeTRIC constraints file for constrained optimization |
| 525 | + |
| 526 | + relative paths are resolved relative to the OpenQP input file |
| 527 | + |
| 528 | +- enforce // geomeTRIC constraint enforcement value |
| 529 | + |
| 530 | + 0.0 (default) |
| 531 | + |
| 532 | +- conmethod // geomeTRIC constraint method |
| 533 | + |
| 534 | + 0 (default) |
| 535 | + |
466 | 536 | ### [dlfind] |
467 | 537 |
|
468 | 538 | dlfind section handle the DL-FIND library for geometry optimization |
|
0 commit comments