Usage examples
Calculate fast DFT (B97-D3(0)/def2-SV(P)+gcp) single-point energies on GFNn-xTB input geometries:
$ censo -inp ensemble.xyz -part0 on -chrg 1 -solvent h2o > censo.out &
______________________________________________________________
| |
| |
| CENSO - Commandline ENSO |
| v 1.0.3 |
| energetic sorting of CREST Conformer Rotamer Ensembles |
| University of Bonn, MCTC |
| Feb 2021 |
| based on ENSO version 2.0.1 |
| F. Bohle and S. Grimme |
| |
|______________________________________________________________|
This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
----------------------------------------------------------------------------------------------------
PARAMETERS
----------------------------------------------------------------------------------------------------
The configuration file .censorc is read from /home/bohle/1projects/from_tmp1/CENSO/documentation-calcs/1-part0/.censorc.
Reading conformer rotamer ensemble from: /home/bohle/1projects/from_tmp1/CENSO/documentation-calcs/1-part0/ensemble.xyz.
Reading file: censo_solvents.json
--------------------------------------------------
CRE SORTING SETTINGS
--------------------------------------------------
number of atoms in system: 25
number of considered conformers: 22
number of all conformers from input: 22
charge: 1
unpaired: 0
solvent: h2o
temperature: 298.15
evaluate at different temperatures: on
temperature range: 273.15, 278.15, 283.15, 288.15, ...
calculate mRRHO contribution: on
consider symmetry for mRRHO contribution: off
cautious checking for error and failed calculations: on
checking the DFT-ensemble using CREST: off
maxthreads: 2
omp: 2
--------------------------------------------------
CRE CHEAP-PRESCREENING - PART0
--------------------------------------------------
part0: on
starting number of considered conformers: 22
program for part0: tm
functional for fast single-point: b97-d
basis set for fast single-point: def2-SV(P)
threshold g_thr(0) for sorting in part0: 4.0
Solvent model used with xTB: alpb
short-notation:
b97-d-D3/def2-SV(P) // GFNn-xTB (Input geometry)
END of parameters
------------------------------------------------------------
PATHS of external QM programs
------------------------------------------------------------
The following program paths are used:
xTB: /home/abt-grimme/AK-bin/xtb
TURBOMOLE: /home/abt-grimme/TURBOMOLE.7.5.1/bin/em64t-unknown-linux-gnu_smp
Using cefine from /home/bohle/bin/cefine
PARNODES for TM or COSMO-RS calculation was set to 2
----------------------------------------------------------------------------------------------------
Processing data from previous run (enso.json)
----------------------------------------------------------------------------------------------------
INFORMATION: No restart information exists and is created during this run!
----------------------------------------------------------------------------------------------------
CRE CHEAP-PRESCREENING - PART0
----------------------------------------------------------------------------------------------------
program: tm
functional for part0: b97-d
basis set for part0: def2-SV(P)
threshold g_thr(0): 4.0
starting number of considered conformers: 22
Calculating efficient gas-phase single-point energies:
The efficient gas-phase single-point is calculated for:
CONF1, CONF2, CONF3, CONF4, CONF5, CONF6, CONF7, CONF8, CONF9, CONF10, CONF11
CONF12, CONF13, CONF14, CONF15, CONF16, CONF17, CONF18, CONF19, CONF20, CONF21, CONF22
Constructed folders!
Starting 22 ALPB-Gsolv calculations
Running single-point in CONF1/part0_sp
Running single-point in CONF2/part0_sp
Running ALPB_GSOLV calculation in 1-part0/CONF1/part0_sp
Running ALPB_GSOLV calculation in 1-part0/CONF2/part0_sp
Running single-point in CONF3/part0_sp
Running single-point in CONF4/part0_sp
Running ALPB_GSOLV calculation in 1-part0/CONF3/part0_sp
Running ALPB_GSOLV calculation in 1-part0/CONF4/part0_sp
Running single-point in CONF5/part0_sp
Running single-point in CONF6/part0_sp
Running ALPB_GSOLV calculation in 1-part0/CONF5/part0_sp
Running ALPB_GSOLV calculation in 1-part0/CONF6/part0_sp
Running single-point in CONF7/part0_sp
Running single-point in CONF8/part0_sp
Running ALPB_GSOLV calculation in 1-part0/CONF7/part0_sp
Running ALPB_GSOLV calculation in 1-part0/CONF8/part0_sp
Running single-point in CONF9/part0_sp
Running single-point in CONF10/part0_sp
Running ALPB_GSOLV calculation in 1-part0/CONF9/part0_sp
Running ALPB_GSOLV calculation in 1-part0/CONF10/part0_sp
Running single-point in CONF11/part0_sp
Running single-point in CONF12/part0_sp
Running ALPB_GSOLV calculation in 1-part0/CONF11/part0_sp
Running ALPB_GSOLV calculation in 1-part0/CONF12/part0_sp
Running single-point in CONF13/part0_sp
Running single-point in CONF14/part0_sp
Running ALPB_GSOLV calculation in 1-part0/CONF13/part0_sp
Running ALPB_GSOLV calculation in 1-part0/CONF14/part0_sp
Running single-point in CONF15/part0_sp
Running single-point in CONF16/part0_sp
Running ALPB_GSOLV calculation in 1-part0/CONF15/part0_sp
Running ALPB_GSOLV calculation in 1-part0/CONF16/part0_sp
Running single-point in CONF17/part0_sp
Running single-point in CONF18/part0_sp
Running ALPB_GSOLV calculation in 1-part0/CONF17/part0_sp
Running ALPB_GSOLV calculation in 1-part0/CONF18/part0_sp
Running single-point in CONF19/part0_sp
Running single-point in CONF20/part0_sp
Running ALPB_GSOLV calculation in 1-part0/CONF19/part0_sp
Running ALPB_GSOLV calculation in 1-part0/CONF20/part0_sp
Running single-point in CONF21/part0_sp
Running single-point in CONF22/part0_sp
Running ALPB_GSOLV calculation in 1-part0/CONF21/part0_sp
Running ALPB_GSOLV calculation in 1-part0/CONF22/part0_sp
Tasks completed!
The efficient gas-phase single-point was successful for CONF1/part0_sp: E(DFT) = -496.55110270 Gsolv = -0.12279245
The efficient gas-phase single-point was successful for CONF2/part0_sp: E(DFT) = -496.55473079 Gsolv = -0.12241044
The efficient gas-phase single-point was successful for CONF3/part0_sp: E(DFT) = -496.55551554 Gsolv = -0.12099969
The efficient gas-phase single-point was successful for CONF4/part0_sp: E(DFT) = -496.55384665 Gsolv = -0.12156811
The efficient gas-phase single-point was successful for CONF5/part0_sp: E(DFT) = -496.55341065 Gsolv = -0.12251050
The efficient gas-phase single-point was successful for CONF6/part0_sp: E(DFT) = -496.55368828 Gsolv = -0.12175678
The efficient gas-phase single-point was successful for CONF7/part0_sp: E(DFT) = -496.54759208 Gsolv = -0.12365913
The efficient gas-phase single-point was successful for CONF8/part0_sp: E(DFT) = -496.55118920 Gsolv = -0.12200372
The efficient gas-phase single-point was successful for CONF9/part0_sp: E(DFT) = -496.54970422 Gsolv = -0.12167549
The efficient gas-phase single-point was successful for CONF10/part0_sp: E(DFT) = -496.55157272 Gsolv = -0.12001541
The efficient gas-phase single-point was successful for CONF11/part0_sp: E(DFT) = -496.54991969 Gsolv = -0.12296934
The efficient gas-phase single-point was successful for CONF12/part0_sp: E(DFT) = -496.55212770 Gsolv = -0.11751864
The efficient gas-phase single-point was successful for CONF13/part0_sp: E(DFT) = -496.55077718 Gsolv = -0.11854291
The efficient gas-phase single-point was successful for CONF14/part0_sp: E(DFT) = -496.55028792 Gsolv = -0.12090071
The efficient gas-phase single-point was successful for CONF15/part0_sp: E(DFT) = -496.54987374 Gsolv = -0.12311773
The efficient gas-phase single-point was successful for CONF16/part0_sp: E(DFT) = -496.55360018 Gsolv = -0.11274054
The efficient gas-phase single-point was successful for CONF17/part0_sp: E(DFT) = -496.55289531 Gsolv = -0.11378052
The efficient gas-phase single-point was successful for CONF18/part0_sp: E(DFT) = -496.50423172 Gsolv = -0.16468285
The efficient gas-phase single-point was successful for CONF19/part0_sp: E(DFT) = -496.52384365 Gsolv = -0.14390846
The efficient gas-phase single-point was successful for CONF20/part0_sp: E(DFT) = -496.50481476 Gsolv = -0.16519450
The efficient gas-phase single-point was successful for CONF21/part0_sp: E(DFT) = -496.53860751 Gsolv = -0.13212346
The efficient gas-phase single-point was successful for CONF22/part0_sp: E(DFT) = -496.50938650 Gsolv = -0.15826546
----------------------------------------------------------------------------------------------------
Removing high lying conformers by improved energy description
----------------------------------------------------------------------------------------------------
CONF# G [Eh] ΔG [kcal/mol] E [Eh] Gsolv [Eh] Gtot ΔE(DFT) ΔGsolv ΔGtot
GFN2-xTB GFN2-xTB b97-d-D3(0)/def2-SV(P) alpb [Eh] [kcal/mol] [kcal/mol] [kcal/mol]
[ALPB] [ALPB] [gfn2]
CONF1 -34.1746819 0.00 -496.5511027 -0.1227924 -496.6738951 2.28 -0.24 2.04
CONF2 -34.1743917 0.18 -496.5547308 -0.1224104 -496.6771412 0.00 0.00 0.00 <------
CONF3 -34.1742792 0.25 -496.5555155 -0.1209997 -496.6765152 -0.49 0.89 0.39
CONF4 -34.1739821 0.44 -496.5538467 -0.1215681 -496.6754148 0.55 0.53 1.08
CONF5 -34.1740224 0.41 -496.5534106 -0.1225105 -496.6759212 0.83 -0.06 0.77
CONF6 -34.1736309 0.66 -496.5536883 -0.1217568 -496.6754451 0.65 0.41 1.06
CONF7 -34.1725555 1.33 -496.5475921 -0.1236591 -496.6712512 4.48 -0.78 3.70
CONF8 -34.1721876 1.57 -496.5511892 -0.1220037 -496.6731929 2.22 0.26 2.48
CONF9 -34.1719439 1.72 -496.5497042 -0.1216755 -496.6713797 3.15 0.46 3.62
CONF10 -34.1714765 2.01 -496.5515727 -0.1200154 -496.6715881 1.98 1.50 3.48
CONF11 -34.1712638 2.14 -496.5499197 -0.1229693 -496.6728890 3.02 -0.35 2.67
CONF12 -34.1704841 2.63 -496.5521277 -0.1175186 -496.6696463 1.63 3.07 4.70
CONF13 -34.1703687 2.71 -496.5507772 -0.1185429 -496.6693201 2.48 2.43 4.91
CONF14 -34.1694191 3.30 -496.5502879 -0.1209007 -496.6711886 2.79 0.95 3.74
CONF15 -34.1691431 3.48 -496.5498737 -0.1231177 -496.6729915 3.05 -0.44 2.60
CONF16 -34.1664418 5.17 -496.5536002 -0.1127405 -496.6663407 0.71 6.07 6.78
CONF17 -34.1661652 5.34 -496.5528953 -0.1137805 -496.6666758 1.15 5.42 6.57
CONF18 -34.1660003 5.45 -496.5042317 -0.1646828 -496.6689146 31.69 -26.53 5.16
CONF19 -34.1664801 5.15 -496.5238437 -0.1439085 -496.6677521 19.38 -13.49 5.89
CONF20 -34.1658387 5.55 -496.5048148 -0.1651945 -496.6700093 31.32 -26.85 4.48
CONF21 -34.1670164 4.81 -496.5386075 -0.1321235 -496.6707310 10.12 -6.10 4.02
CONF22 -34.1656180 5.69 -496.5093865 -0.1582655 -496.6676520 28.45 -22.50 5.95
----------------------------------------------------------------------------------------------------
--------------------------------------------------
Conformers considered further
--------------------------------------------------
These conformers are below the 4.000 kcal/mol g_thr(0) threshold.
CONF1, CONF2, CONF3, CONF4, CONF5, CONF6, CONF7, CONF8, CONF9, CONF10, CONF11
CONF14, CONF15
Calculating Boltzmann averaged (free) energy of ensemble on input geometries (not DFT optimized)!
temperature /K: avE(T) /a.u. avG(T) /a.u.
----------------------------------------------------------------------------------------------------
298.15 -496.5544580 -496.6764445 <<==part0==
----------------------------------------------------------------------------------------------------
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>END of Part0<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
Ran part0 in 32.8727 seconds
Part : #conf time
--------------------------------------------------
Input : 22 -
Part0_all : 13 32.87s
--------------------------------------------------
All parts : 32.87s
CENSO all done!Calculate free energies in solution phase (CHCl3) on GFNn-xTB geometries:
$ censo -inp ensemble.xyz -part1 on -chrg 1 -solvent chcl3 -smgsolv1 cosmors -func r2scan-3c -basis automatic -P 4 -O 2 > censo.out &
______________________________________________________________
| |
| |
| CENSO - Commandline ENSO |
| v 1.0.3 |
| energetic sorting of CREST Conformer Rotamer Ensembles |
| University of Bonn, MCTC |
| Feb 2021 |
| based on ENSO version 2.0.1 |
| F. Bohle and S. Grimme |
| |
|______________________________________________________________|
This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
----------------------------------------------------------------------------------------------------
PARAMETERS
----------------------------------------------------------------------------------------------------
The configuration file .censorc is read from /home/bohle/1projects/from_tmp1/CENSO/documentation-calcs/2-part1/.censorc.
Reading conformer rotamer ensemble from: /home/bohle/1projects/from_tmp1/CENSO/documentation-calcs/2-part1/ensemble.xyz.
Reading file: censo_solvents.json
--------------------------------------------------
CRE SORTING SETTINGS
--------------------------------------------------
number of atoms in system: 25
number of considered conformers: 22
number of all conformers from input: 22
charge: 1
unpaired: 0
solvent: chcl3
temperature: 298.15
evaluate at different temperatures: on
temperature range: 273.15, 278.15, 283.15, 288.15, ...
calculate mRRHO contribution: on
consider symmetry for mRRHO contribution: off
cautious checking for error and failed calculations: on
checking the DFT-ensemble using CREST: off
maxthreads: 4
omp: 2
--------------------------------------------------
CRE PRESCREENING - PART1
--------------------------------------------------
part1: on
starting number of considered conformers: 22
program for part1: tm
functional for initial evaluation: r2scan-3c
basis set for initial evaluation: def2-mTZVPP
calculate mRRHO contribution: on
program for mRRHO contribution: xtb
GFN version for mRRHO and/or GBSA_Gsolv: gfn2
Apply constraint to input geometry during mRRHO calculation: on
solvent model applied with xTB: alpb
evaluate at different temperatures: off
threshold g_thr(1) and G_thr(1) for sorting in part1: 3.5
solvent model for Gsolv contribution of part1: cosmors
short-notation:
r2scan-3c + COSMORS[chcl3] + GmRRHO(GFN2[alpb]-bhess) // GFNn-xTB (Input geometry)
END of parameters
------------------------------------------------------------
PATHS of external QM programs
------------------------------------------------------------
The following program paths are used:
xTB: /home/abt-grimme/AK-bin/xtb
TURBOMOLE: /home/abt-grimme/TURBOMOLE.7.5.1/bin/em64t-unknown-linux-gnu_smp
Setup of COSMO-RS:
ctd = BP_TZVPD_FINE_C30_1601.ctd
cdir = "/home/bohle/COSMOlogic/COSMOthermX19/COSMOtherm/CTDATA-FILES"
ldir = "/home/bohle/COSMOlogic/COSMOthermX19/COSMOtherm/CTDATA-FILES"
Using /home/bohle/COSMOlogic/COSMOthermX19/COSMOtherm/DATABASE-COSMO/BP-TZVP-COSMO
as path to the COSMO-RS NORMAL DATABASE.
Using cefine from /home/bohle/bin/cefine
PARNODES for TM or COSMO-RS calculation was set to 2
Using COSMOtherm from /home/bohle/COSMOlogic/COSMOthermX19/COSMOtherm/BIN-LINUX/cosmotherm
----------------------------------------------------------------------------------------------------
Processing data from previous run (enso.json)
----------------------------------------------------------------------------------------------------
INFORMATION: No restart information exists and is created during this run!
----------------------------------------------------------------------------------------------------
CRE PRESCREENING - PART1
----------------------------------------------------------------------------------------------------
program: tm
functional for part1 and 2: r2scan-3c
basis set for part1 and 2: def2-mTZVPP
Solvent: chcl3
solvent model for Gsolv contribution: cosmors
threshold g_thr(1) and G_thr(1): 3.5
starting number of considered conformers: 22
calculate mRRHO contribution: on
program for mRRHO contribution: xtb
GFN version for mRRHO and/or GBSA_Gsolv: gfn2
Apply constraint to input geometry during mRRHO calculation: on
temperature: 298.15
Calculating single-point energies and solvation contribution (G_solv):
The prescreening COSMO-RS is calculated for:
CONF1, CONF2, CONF3, CONF4, CONF5, CONF6, CONF7, CONF8, CONF9, CONF10, CONF11
CONF12, CONF13, CONF14, CONF15, CONF16, CONF17, CONF18, CONF19, CONF20, CONF21, CONF22
Constructed folders!
Constructed folders!
Starting 22 COSMO-RS-Gsolv calculations.
Running COSMO-RS calculation in CONF1/r2scan-3c/COSMO
Running COSMO-RS calculation in CONF2/r2scan-3c/COSMO
Running COSMO-RS calculation in CONF3/r2scan-3c/COSMO
Running COSMO-RS calculation in CONF4/r2scan-3c/COSMO
Running COSMO-RS calculation in CONF5/r2scan-3c/COSMO
Running COSMO-RS calculation in CONF6/r2scan-3c/COSMO
Running COSMO-RS calculation in CONF7/r2scan-3c/COSMO
Running COSMO-RS calculation in CONF8/r2scan-3c/COSMO
Running COSMO-RS calculation in CONF9/r2scan-3c/COSMO
Running COSMO-RS calculation in CONF10/r2scan-3c/COSMO
Running COSMO-RS calculation in CONF11/r2scan-3c/COSMO
Running COSMO-RS calculation in CONF12/r2scan-3c/COSMO
Running COSMO-RS calculation in CONF13/r2scan-3c/COSMO
Running COSMO-RS calculation in CONF14/r2scan-3c/COSMO
Running COSMO-RS calculation in CONF15/r2scan-3c/COSMO
Running COSMO-RS calculation in CONF16/r2scan-3c/COSMO
Running COSMO-RS calculation in CONF17/r2scan-3c/COSMO
Running COSMO-RS calculation in CONF18/r2scan-3c/COSMO
Running COSMO-RS calculation in CONF19/r2scan-3c/COSMO
Running COSMO-RS calculation in CONF20/r2scan-3c/COSMO
Running COSMO-RS calculation in CONF21/r2scan-3c/COSMO
Running COSMO-RS calculation in CONF22/r2scan-3c/COSMO
Tasks completed!
prescreening COSMO-RS calculation was successful for CONF1/r2scan-3c/COSMO: -0.09391084
prescreening COSMO-RS calculation was successful for CONF2/r2scan-3c/COSMO: -0.09344716
prescreening COSMO-RS calculation was successful for CONF3/r2scan-3c/COSMO: -0.09189166
prescreening COSMO-RS calculation was successful for CONF4/r2scan-3c/COSMO: -0.09269361
prescreening COSMO-RS calculation was successful for CONF5/r2scan-3c/COSMO: -0.09328873
prescreening COSMO-RS calculation was successful for CONF6/r2scan-3c/COSMO: -0.09278491
prescreening COSMO-RS calculation was successful for CONF7/r2scan-3c/COSMO: -0.09457155
prescreening COSMO-RS calculation was successful for CONF8/r2scan-3c/COSMO: -0.09275254
prescreening COSMO-RS calculation was successful for CONF9/r2scan-3c/COSMO: -0.09316162
prescreening COSMO-RS calculation was successful for CONF10/r2scan-3c/COSMO: -0.09116989
prescreening COSMO-RS calculation was successful for CONF11/r2scan-3c/COSMO: -0.09348580
prescreening COSMO-RS calculation was successful for CONF12/r2scan-3c/COSMO: -0.08900263
prescreening COSMO-RS calculation was successful for CONF13/r2scan-3c/COSMO: -0.08998048
prescreening COSMO-RS calculation was successful for CONF14/r2scan-3c/COSMO: -0.09122770
prescreening COSMO-RS calculation was successful for CONF15/r2scan-3c/COSMO: -0.09272288
prescreening COSMO-RS calculation was successful for CONF16/r2scan-3c/COSMO: -0.08509833
prescreening COSMO-RS calculation was successful for CONF17/r2scan-3c/COSMO: -0.08598277
prescreening COSMO-RS calculation was successful for CONF18/r2scan-3c/COSMO: -0.13224220
prescreening COSMO-RS calculation was successful for CONF19/r2scan-3c/COSMO: -0.11057302
prescreening COSMO-RS calculation was successful for CONF20/r2scan-3c/COSMO: -0.13223689
prescreening COSMO-RS calculation was successful for CONF21/r2scan-3c/COSMO: -0.09915029
prescreening COSMO-RS calculation was successful for CONF22/r2scan-3c/COSMO: -0.12536341
--------------------------------------------------
Removing high lying conformers
--------------------------------------------------
CONF# E(GFNn-xTB) ΔE(GFNn-xTB) E [Eh] Gsolv [Eh] Gtot ΔGtot
[a.u.] [kcal/mol] r2scan-3c COSMO-RS-normal [Eh] [kcal/mol]
[r2scan-3c]
CONF1 -34.1599548 0.00 -497.3021467 -0.0939108 -497.3960575 1.82
CONF2 -34.1594484 0.32 -497.3049956 -0.0934472 -497.3984428 0.33
CONF3 -34.1592530 0.44 -497.3064930 -0.0918917 -497.3983847 0.36
CONF4 -34.1590954 0.54 -497.3049058 -0.0926936 -497.3975994 0.86
CONF5 -34.1589250 0.65 -497.3056747 -0.0932887 -497.3989634 0.00 <------
CONF6 -34.1587465 0.76 -497.3049670 -0.0927849 -497.3977519 0.76
CONF7 -34.1578065 1.35 -497.2983644 -0.0945716 -497.3929360 3.78
CONF8 -34.1571241 1.78 -497.3029570 -0.0927525 -497.3957096 2.04
CONF9 -34.1570372 1.83 -497.3008248 -0.0931616 -497.3939864 3.12
CONF10 -34.1563217 2.28 -497.3038986 -0.0911699 -497.3950685 2.44
CONF11 -34.1559251 2.53 -497.3008250 -0.0934858 -497.3943108 2.92
CONF12 -34.1554750 2.81 -497.3045828 -0.0890026 -497.3935854 3.37
CONF13 -34.1554585 2.82 -497.3036749 -0.0899805 -497.3936554 3.33
CONF14 -34.1537285 3.91 -497.3022176 -0.0912277 -497.3934453 3.46
CONF15 -34.1529637 4.39 -497.3009381 -0.0927229 -497.3936610 3.33
CONF16 -34.1514138 5.36 -497.3063638 -0.0850983 -497.3914621 4.71
CONF17 -34.1512255 5.48 -497.3047147 -0.0859828 -497.3906975 5.19
CONF18 -34.1510249 5.60 -497.2564102 -0.1322422 -497.3886524 6.47
CONF19 -34.1509627 5.64 -497.2746319 -0.1105730 -497.3852050 8.63
CONF20 -34.1508506 5.71 -497.2550494 -0.1322369 -497.3872863 7.33
CONF21 -34.1504892 5.94 -497.2896633 -0.0991503 -497.3888135 6.37
CONF22 -34.1504331 5.98 -497.2596139 -0.1253634 -497.3849773 8.78
--------------------------------------------------
Conformers considered further
--------------------------------------------------
Below the g_thr(1) threshold of 3.5 kcal/mol.
CONF1, CONF2, CONF3, CONF4, CONF5, CONF6, CONF8, CONF9, CONF10, CONF11, CONF12
CONF13, CONF14, CONF15
--------------------------------------------------
Calculating prescreening G_mRRHO with implicit solvation!
The prescreening G_mRRHO calculation is now performed for:
CONF1, CONF2, CONF3, CONF4, CONF5, CONF6, CONF8, CONF9, CONF10, CONF11, CONF12
CONF13, CONF14, CONF15
Constructed folders!
Starting 14 G_RRHO calculations.
Running GFN2-xTB mRRHO in CONF1/rrho_part1
Running GFN2-xTB mRRHO in CONF2/rrho_part1
Running GFN2-xTB mRRHO in CONF3/rrho_part1
Running GFN2-xTB mRRHO in CONF4/rrho_part1
Running GFN2-xTB mRRHO in CONF5/rrho_part1
Running GFN2-xTB mRRHO in CONF6/rrho_part1
Running GFN2-xTB mRRHO in CONF8/rrho_part1
Running GFN2-xTB mRRHO in CONF9/rrho_part1
Running GFN2-xTB mRRHO in CONF10/rrho_part1
Running GFN2-xTB mRRHO in CONF11/rrho_part1
Running GFN2-xTB mRRHO in CONF12/rrho_part1
Running GFN2-xTB mRRHO in CONF13/rrho_part1
Running GFN2-xTB mRRHO in CONF14/rrho_part1
Running GFN2-xTB mRRHO in CONF15/rrho_part1
Tasks completed!
The prescreening G_mRRHO calculation @ c1 was successful for CONF1/rrho_part1: 0.18236493
The prescreening G_mRRHO calculation @ c1 was successful for CONF2/rrho_part1: 0.18246714
The prescreening G_mRRHO calculation @ c1 was successful for CONF3/rrho_part1: 0.18235782
The prescreening G_mRRHO calculation @ c1 was successful for CONF4/rrho_part1: 0.18239049
The prescreening G_mRRHO calculation @ c1 was successful for CONF5/rrho_part1: 0.18241907
The prescreening G_mRRHO calculation @ c1 was successful for CONF6/rrho_part1: 0.18225430
The prescreening G_mRRHO calculation @ c1 was successful for CONF8/rrho_part1: 0.18193879
The prescreening G_mRRHO calculation @ c1 was successful for CONF9/rrho_part1: 0.18214342
The prescreening G_mRRHO calculation @ c1 was successful for CONF10/rrho_part1: 0.18217367
The prescreening G_mRRHO calculation @ c1 was successful for CONF11/rrho_part1: 0.18267361
The prescreening G_mRRHO calculation @ c1 was successful for CONF12/rrho_part1: 0.18167483
The prescreening G_mRRHO calculation @ c1 was successful for CONF13/rrho_part1: 0.18145764
The prescreening G_mRRHO calculation @ c1 was successful for CONF14/rrho_part1: 0.18268845
The prescreening G_mRRHO calculation @ c1 was successful for CONF15/rrho_part1: 0.18320604
--------------------------------------------------
* Gibbs free energies of part1 *
--------------------------------------------------
CONF# G(GFNn-xTB) ΔG(GFNn-xTB) E [Eh] Gsolv [Eh] GmRRHO [Eh] Gtot ΔGtot
[a.u.] [kcal/mol] r2scan-3c COSMO-RS-normal GFN2 [Eh] [kcal/mol]
[r2scan-3c] [alpb]-bhess
CONF1 -33.9775899 0.00 -497.3021467 -0.0939108 0.1823649 -497.2136926 1.79
CONF2 -33.9769813 0.38 -497.3049956 -0.0934472 0.1824671 -497.2159757 0.36
CONF3 -33.9768952 0.44 -497.3064930 -0.0918917 0.1823578 -497.2160268 0.32
CONF4 -33.9767049 0.56 -497.3049058 -0.0926936 0.1823905 -497.2152089 0.84
CONF5 -33.9765060 0.68 -497.3056747 -0.0932887 0.1824191 -497.2165444 0.00 <------
CONF6 -33.9764922 0.69 -497.3049670 -0.0927849 0.1822543 -497.2154976 0.66
CONF8 -33.9751853 1.51 -497.3029570 -0.0927525 0.1819388 -497.2137708 1.74
CONF9 -33.9748938 1.69 -497.3008248 -0.0931616 0.1821434 -497.2118430 2.95
CONF10 -33.9741480 2.16 -497.3038986 -0.0911699 0.1821737 -497.2128948 2.29
CONF11 -33.9732515 2.72 -497.3008250 -0.0934858 0.1826736 -497.2116372 3.08
CONF12 -33.9738002 2.38 -497.3045828 -0.0890026 0.1816748 -497.2119106 2.91
CONF13 -33.9740008 2.25 -497.3036749 -0.0899805 0.1814576 -497.2121977 2.73
CONF14 -33.9710400 4.11 -497.3022176 -0.0912277 0.1826885 -497.2107568 3.63
CONF15 -33.9697576 4.91 -497.3009381 -0.0927229 0.1832060 -497.2104549 3.82
Additional global 'fuzzy-threshold' based on the standard deviation of (G_mRRHO):
Std_dev(G_mRRHO) = 0.272 kcal/mol
Fuzzythreshold = 0.309 kcal/mol
Final sorting threshold G_thr(1) = 3.500 + 0.309 = 3.809 kcal/mol
Spearman correlation coefficient between (E + Solv) and (E + Solv + mRRHO) = 0.908
--------------------------------------------------
Conformers considered further
--------------------------------------------------
Considered CONF14 because of increased fuzzythr.
These conformers are below the 3.809 kcal/mol G_thr(1) threshold.
CONF1, CONF2, CONF3, CONF4, CONF5, CONF6, CONF8, CONF9, CONF10, CONF11, CONF12
CONF13, CONF14
Calculating Boltzmann averaged free energy of ensemble on input geometries (not DFT optimized)!
temperature /K: avE(T) /a.u. avGmRRHO(T) /a.u. avGsolv(T) /a.u. avG(T) /a.u.
----------------------------------------------------------------------------------------------------
298.15 -497.3054054 0.1823775 -0.0928881 -497.2159160 <<==part1==
----------------------------------------------------------------------------------------------------
Calculating unbiased GFNn-xTB energy
Constructed folders!
Starting 13 xTB - single-point calculations.
gfn2-xTB energy for CONF1/GFN_unbiased = -34.1515395
gfn2-xTB energy for CONF2/GFN_unbiased = -34.1509941
gfn2-xTB energy for CONF3/GFN_unbiased = -34.1515326
gfn2-xTB energy for CONF4/GFN_unbiased = -34.1508491
gfn2-xTB energy for CONF5/GFN_unbiased = -34.1514366
gfn2-xTB energy for CONF6/GFN_unbiased = -34.1512893
gfn2-xTB energy for CONF8/GFN_unbiased = -34.1493072
gfn2-xTB energy for CONF9/GFN_unbiased = -34.1491025
gfn2-xTB energy for CONF10/GFN_unbiased = -34.1491972
gfn2-xTB energy for CONF11/GFN_unbiased = -34.1487085
gfn2-xTB energy for CONF12/GFN_unbiased = -34.1484725
gfn2-xTB energy for CONF13/GFN_unbiased = -34.1480733
gfn2-xTB energy for CONF14/GFN_unbiased = -34.1461266
Tasks completed!
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>END of Part1<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
Ran part1 in 277.0194 seconds
Part : #conf time
--------------------------------------------------
Input : 22 -
Part1_initial_sort : 14 -
Part1_all : 14 277.02s
--------------------------------------------------
All parts : 277.02s
CENSO all done!
Calculate free energies on populated, DFT optimized conformers
Restarting calculations
Last updated