neb_utils.tcl

::pwtk::manual_neb

::pwtk::manual_neb

SYNOPSIS

proc ::pwtk::manual_neb {args} {

USAGE

   manual_neb  ?-noCI_thr THRESHOLD?  ?-CI_thr THRESHOLD?  ?-CI CI_LIST?  ?-onlyCI?  head

PURPOSE

Perform no-CI + manual-CI NEB calculations.

First, the "no-CI" NEB is performed up to a reduced threshold to stabilize the minimum-energy path and then CI-NEB is restarted with CI_scheme = 'manual' using a full-precision threshold.

OPTIONS

-noCI_thr THRESHOLD ... path threshold for the preliminary no-CI NEB calculation; if not specified a default preliminary-threshold is used.
-CI_thr THRESHOLD ... path threshold for the CI-NEB calculation; if not specified, the current path_thr is used, but if path_thr is not set, the neb.x default is used.
-CI CI_LIST ... list of climbing-images for the CLIMBING_IMAGES card
-onlyCI ... skip the 'no-CI' calculation, i.e., perform only CI calculation

ARGUMENTS

head ... the rootname for the I/O files

SIDE EFFECTS

"prefix" and outdir_postfix are set to "head", whereas all other modifications are poped away after completion because both calculations are performed with input_pushpop.

SOURCE

    variable noCI_path_thr
    variable neb_options

    set narg 1
    set options [concat $neb_options {
        {CI.arg   {} "list of climbing images for the CLIMBING_IMAGES card"}
    }]
    set usage "?-noCI_thr THRESHOLD?  ?-CI_thr THRESHOLD?  ?-CI CI_LIST?  ?-onlyCI?  head"
    ::pwtk::parseOpt_
    ::pwtk::checkOType_ -noCI_thr $opt(noCI_thr) {number posreal} "positive real number"
    ::pwtk::checkOType_ -CI_thr   $opt(CI_thr)   {number posreal} "positive real number"
    ::pwtk::checkOType_ -CI       $opt(CI)    {numberlist posint} "list of positive integer numbers"
    set head $args

    ifset opt(noCI_thr) $noCI_path_thr

    outdir_postfix $head
    CONTROL " prefix = '$head' "

    ::pwtk::input::pushpop {
        # 1. run no-CI NEB
        if { ! $opt(onlyCI) } {
            ::pwtk::input::pushpop {
                PATH "
                   CI_scheme = 'no-CI'
                   path_thr  = $opt(noCI_thr)
                "
                runNEB neb.$head.no-CI
            }
            PATH {
                restart_mode = 'restart'
            }
        }

        # 2. run manual CI-NEB
        PATH {
            CI_scheme    = 'manual'
        }
        ifnotempty opt(CI_thr) {
            PATH " path_thr = $opt(CI_thr) "
        }
        ifnotempty opt(CI) {
            CLIMBING_IMAGES $opt(CI)
        }
        runNEB neb.$head.manual
    }    
}

TABLE OF CONTENTS

::pwtk::manual_neb