Skip to content
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

ONIOM #763

Merged
merged 35 commits into from
Feb 1, 2023
Merged

ONIOM #763

Show file tree
Hide file tree
Changes from 29 commits
Commits
Show all changes
35 commits
Select commit Hold shift + click to select a range
df4e62d
small changes to .gitignore
Nov 7, 2022
46549a5
Allocation of solvation
Nov 7, 2022
4b5a71b
Default methods (gfn2:gfnff) for the ONIOM calculation
Nov 7, 2022
5c84053
.gitignore
Nov 8, 2022
ac5e428
license update for man/xtb.1.adoc
Nov 22, 2022
c4ce55a
ONIOM Jacobian matrix
Bugfixe Nov 24, 2022
e096d11
Changes related to the ONIOM
Bugfixe Jan 26, 2023
a349831
Turbomole settings for ONIOM
Bugfixe Jan 26, 2023
90b855f
general changes
Bugfixe Jan 26, 2023
ad4851a
Comments
Bugfixe Jan 26, 2023
5e64b84
inner:outer charges
Bugfixe Jan 26, 2023
6cf9986
To change setparam for ONIOM
Bugfixe Jan 26, 2023
f7d39c6
allocate wfn for oniom in main
Bugfixe Jan 26, 2023
d0743a2
Merge branch 'main' of github.com:Albkat/xtb into mult
Bugfixe Jan 26, 2023
ebcc66c
oniom_settings data type to control ONIOM
Bugfixe Jan 27, 2023
9e775a4
ONIOM: option to print inner region logs
Bugfixe Jan 28, 2023
6dae69b
ONIOM: to change k factor from fixed to derive
Bugfixe Jan 28, 2023
e525dcf
add ONIOM specs in man
Bugfixe Jan 29, 2023
86cc025
Add --ceasfiles to delete all external files
Bugfixe Jan 29, 2023
aaba942
add 'mute' to redirect ORCA and TURBOMOLE output
Bugfixe Jan 29, 2023
c20eb8d
add protectCoord subroutine:
Bugfixe Jan 29, 2023
56b04c3
add silent into xcontrol
Bugfixe Jan 29, 2023
27df558
Pretty headers
Bugfixe Jan 29, 2023
cbd066d
fixed rdrm reference
Bugfixe Jan 29, 2023
571bf65
Cleaning the source
Bugfixe Jan 30, 2023
51d8f35
fixed class default statement
Bugfixe Jan 30, 2023
6c43a11
Update bias_path.f90
cplett Jan 30, 2023
43bd990
Update optimizer.f90
cplett Jan 30, 2023
689493a
Update scf_module.F90
cplett Jan 30, 2023
41dd08d
Update anc.f90
Albkat Jan 30, 2023
0fcd044
Update main.F90
Albkat Jan 30, 2023
130a842
tests added
Bugfixe Feb 1, 2023
eea3011
new test for iniom
Bugfixe Feb 1, 2023
27b2941
Update xtb.1.adoc
Albkat Feb 1, 2023
e76c559
Update orca.f90
Albkat Feb 1, 2023
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
13 changes: 13 additions & 0 deletions man/xcontrol.7.adoc
Original file line number Diff line number Diff line change
Expand Up @@ -326,6 +326,19 @@ $modef
*mode*='int'::
can set by --modef via cmdline

$oniom
~~~~~~
*inner logs*='bool'::
to print optimization log files for model region geometry (`high.inner_region.log` and `low.inner_region.log`)

*derived k*='bool'::
to calculate prefactor *k* and create jacobian dynamically (see more )

*silent*='bool'::
to hide the execution runs of external software



$opt
~~~~
*engine*='method'::
Expand Down
18 changes: 12 additions & 6 deletions man/xtb.1.adoc
Original file line number Diff line number Diff line change
Expand Up @@ -24,8 +24,8 @@ geometries, frequencies and non-covalent interactions (GFN)
as well as with the ionisation potential and
electron affinity (IPEA) parametrisation of the GFN1 Hamiltonian.
The generalized born (GB) model with solvent accessable surface area (SASA)
is also available available in this version.
Ground state calculations for the simplified Tamm-Danceoff approximation (sTDA)
is also available in this version.
Ground state calculations for the simplified Tamm-Dancoff approximation (sTDA)
with the vTB model are currently not implemented.

GEOMETRY INPUT
Expand Down Expand Up @@ -66,6 +66,9 @@ OPTIONS
*-c, --chrg* 'INT'::
specify molecular charge as 'INT', overrides `.CHRG` file and `xcontrol` option

*-c, --chrg* 'INT:INT'::
specify charges for 'inner region:outer region' for `oniom` calculation, overrides `.CHRG` file and `xcontrol` option

*-u, --uhf* 'INT'::
specify number of unpaired electrons as 'INT', overrides `.UHF` file and `xcontrol` option

Expand All @@ -79,9 +82,9 @@ OPTIONS
use tblite library as implementation of for xTB

*--oniom* 'METHOD' 'LIST'::
use subtractive embedding via ONIOM method. 'METHOD' is given as `inner:outer`
where `inner` can be 'orca', 'turbomole', 'gfn2', 'gfn1', or 'gfnff' and
`outer` can be 'gfn2', 'gfn1', or 'gfnff'.
use subtractive embedding via ONIOM method. 'METHOD' is given as `high:low`
where `high` can be 'orca', 'turbomole', 'gfn2', 'gfn1', or 'gfnff' and
`low` can be 'gfn2', 'gfn1', or 'gfnff'.
The inner region is given as a comma separated indices directly in the commandline
or in a file with each index on a separate line.

Expand Down Expand Up @@ -263,6 +266,9 @@ GENERAL
*-h, --help*::
show help page

*--cut*::
create inner region for `oniom` calculation without performing any calcultion

ENVIRONMENT VARIABLES
---------------------
`xtb(1)` accesses a path-like variable to determine the location of its
Expand Down Expand Up @@ -421,7 +427,7 @@ After `xtb(1)` has evaluated the all input sources it immediately enters
the production mode. Severe errors will lead to an abnormal termination
which is signalled by the printout to STDERR and a non-zero return value
(usually 128). All non-fatal errors are summerized in the end of the calculation
in one block, right bevor the timing analysis.
in one block, right before the timing analysis.

To aid the user to fix the problems generating these warnings a brief
summary of each warning with its respective string representation in the
Expand Down
2 changes: 1 addition & 1 deletion src/bias_path.f90
Original file line number Diff line number Diff line change
Expand Up @@ -282,7 +282,7 @@ subroutine bias_path(env, mol, chk, calc, egap, et, maxiter, epot, grd, sigma)
enddo
factor2= factor2 + alp_change
enddo bias_loop

!! ------------------------------------------------------------------------
! output and find path yielding product
!! ------------------------------------------------------------------------
Expand Down
4 changes: 2 additions & 2 deletions src/extern/driver.f90
Original file line number Diff line number Diff line change
Expand Up @@ -117,7 +117,7 @@ subroutine singlepoint(self, env, mol, chk, printlevel, restart, &
!$omp critical (turbo_lock)
inquire (file='gradient', exist=exist)
if (exist) then
call rdtm(mol%n, .true., energy, gradient, xyz_cached)
call rdtm(env,mol%n, .true., energy, gradient, xyz_cached)
cache = all(abs(xyz_cached - mol%xyz) < 1.e-10_wp)
end if
if (.not. cache) then
Expand All @@ -135,7 +135,7 @@ subroutine singlepoint(self, env, mol, chk, printlevel, restart, &
end if
write (env%unit, '(72("="))')

call rdtm(mol%n, .true., energy, gradient, xyz_cached)
call rdtm(env,mol%n, .true., energy, gradient, xyz_cached)
end if
!$omp end critical (turbo_lock)

Expand Down
Loading