-
Notifications
You must be signed in to change notification settings - Fork 148
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
aISS docking update for next QCG version #956
Conversation
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks good
src/docking/set_module.f90
Outdated
@@ -245,11 +245,15 @@ subroutine set_directed(env,key,val,nat,at,idMap,xyz) | |||
integer, allocatable :: list(:) | |||
real(wp) :: ddum | |||
logical :: ldum | |||
real(wp), parameter ::au = 627.509541d0 |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
real(wp), parameter ::au = 627.509541d0 | |
use xtb_mctc_convert, only : autokcal |
src/docking/set_module.f90
Outdated
case('expo') | ||
if (getValue(env,val,ddum).and.set2) directedset%expo(1) = ddum | ||
set2 = .false. | ||
case('prefac') | ||
!Prefactor is given in kcal, but xtb energies are in hartree | ||
if (getValue(env,val,ddum).and.set3) directedset%val(1) = ddum / au | ||
set3 = .false. | ||
case('midpoint') | ||
if (getValue(env,val,ddum).and.set4) directedset%expo(2) = ddum | ||
set4 = .false. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Since these instructions aren't documented in xcontrol.7.adoc, it might be beneficial to add them to the xtb-docs
Signed-off-by: cplett <[email protected]>
* Update aISS docking for next generation QCG Signed-off-by: cplett <[email protected]> * Cleanup Signed-off-by: cplett <[email protected]> --------- Signed-off-by: cplett <[email protected]> Signed-off-by: Johannes Gorges <[email protected]>
Introduction of an optional attractive potential during the position screening of the
aISS
algorithm. This will be used to make the solute more attractive for QCG runs. It has the shape of an error function including three parameters that can be used to determine the strength, the width, and the turning point. They are provided via an input file and stored in thedirectedset
object.