kinbot 2.2.2__tar.gz → 2.2.3__tar.gz
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- {kinbot-2.2.2 → kinbot-2.2.3}/KinBot.egg-info/PKG-INFO +31 -32
- {kinbot-2.2.2 → kinbot-2.2.3}/KinBot.egg-info/requires.txt +2 -2
- {kinbot-2.2.2 → kinbot-2.2.3}/PKG-INFO +31 -32
- kinbot-2.2.3/README.md +131 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/find_motif.py +14 -7
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/fragments.py +24 -10
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/modify_geom.py +4 -2
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/optimize.py +4 -2
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/parameters.py +1 -1
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/pes.py +1 -1
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/pp_settings.py +1 -1
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reac_family.py +1 -1
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reaction_generator.py +19 -41
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/stationary_pt.py +7 -1
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_sella_vts.tpl.py +8 -2
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/rotdPy.tpl +4 -3
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/utils.py +37 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/vrc_tst_scan.py +7 -35
- {kinbot-2.2.2 → kinbot-2.2.3}/pyproject.toml +3 -3
- kinbot-2.2.2/README.md +0 -132
- {kinbot-2.2.2 → kinbot-2.2.3}/KinBot.egg-info/SOURCES.txt +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/KinBot.egg-info/dependency_links.txt +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/KinBot.egg-info/entry_points.txt +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/KinBot.egg-info/top_level.txt +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/LICENSE +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/__init__.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/ase_modules/__init__.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/ase_modules/calculators/__init__.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/ase_modules/calculators/gaussian.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/ase_modules/calculators/nn_pes.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/ase_modules/calculators/qchem.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/ase_modules/constraints.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/ase_modules/io/__init__.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/ase_modules/io/formats.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/ase_modules/io/gaussian.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/ase_modules/io/zmatrix.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/bfgs.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/bond_combinations.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/cheminfo.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/config_log.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/conformers.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/constants.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/exceptions.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/frequencies.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/geometry.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/hindered_rotors.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/irc.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/kb.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/license_message.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/mesmer.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/mess.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/molpro.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/orca.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/postprocess.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/pp_tables.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/qc.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reac_General.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reaction_finder.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reaction_finder_bimol.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/__init__.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_12_shift_S_F.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_12_shift_S_R.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_Cyclic_Ether_Formation.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_Diels_alder_addition.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_HO2_Elimination_from_PeroxyRadical.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_Intra_Diels_alder_R.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_Intra_RH_Add_Endocyclic_F.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_Intra_RH_Add_Endocyclic_R.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_Intra_RH_Add_Exocyclic_F.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_Intra_RH_Add_Exocyclic_R.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_Intra_R_Add_Endocyclic_F.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_Intra_R_Add_ExoTetCyclic_F.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_Intra_R_Add_Exocyclic_F.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_Intra_disproportionation_F.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_Intra_disproportionation_R.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_Korcek_step2.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_Korcek_step2_even.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_Korcek_step2_odd.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_R_Addition_COm3_R.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_R_Addition_CSm_R.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_R_Addition_MultipleBond.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_Retro_Ene.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_abstraction.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_barrierless_saddle.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_beta_delta.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_bimol_disproportionation_R.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_birad_recombination_F.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_birad_recombination_R.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_combinatorial.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_cpd_H_migration.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_h2_elim.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_homolytic_scission.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_intra_H_migration.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_intra_H_migration_suprafacial.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_intra_OH_migration.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_intra_OH_migration_Exocyclic_F.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_intra_R_migration.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_ketoenol.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_r12_cycloaddition.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_r12_insertion_R.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_r13_insertion_CO2.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_r13_insertion_ROR.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_r13_insertion_RSR.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_r14_birad_scission.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_r14_cyclic_birad_scission_R.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reactions/reac_r22_cycloaddition.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reader_gauss.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/reader_qchem.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/symmetry.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/thread_kinbot.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/__init__.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_gauss_hir.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_gauss_irc.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_gauss_opt_well.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_gauss_ring_conf.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_gauss_ts_end.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_gauss_ts_search.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_gauss_vts.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_nwchem_freq_well.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_nwchem_irc.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_nwchem_opt_well.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_nwchem_ts_end.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_nwchem_ts_search.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_nwchem_ts_search_ase_constraints.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_qchem_hir.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_qchem_irc.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_qchem_opt_well.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_qchem_ts_end.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_qchem_ts_search.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_sella_hir.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_sella_irc.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_sella_opt_well.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_sella_ring_conf.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_sella_ts_end.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/ase_sella_ts_search.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/local_molpro.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_2tst.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_atom.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_barrier.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_barrier_union.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_barrierless.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_bimol.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_bimol_union.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_core_rr.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_dummy.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_fragment.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_fragment_OH.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_freerotor.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_header.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_hinderedrotor.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_hinderedrotorgeom.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_outerrrho.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_pst.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_pstfragment.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_rrho.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_termol.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_ts.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_tunneling.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_variational.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_well.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/mess_well_union.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/molpro.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/molpro_ts_search.tpl.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/molpro_vts.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/pbs.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/pbs_mesmer.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/pbs_mess.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/pbs_mess_uq.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/pbs_molpro.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/pbs_python.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/pesviewer.inp.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/rotdPy_1d_corr.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/rotdPy_calc.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/rotdPy_frag.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/rotdPy_surf.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/slurm.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/slurm_mesmer.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/slurm_mess.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/slurm_mess_uq.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/slurm_molpro.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/slurm_orca.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/tpl/slurm_python.tpl +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/uncertaintyAnalysis.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/vrc_tst_surfaces.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/kinbot/zmatrix.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/setup.cfg +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/tests/__init__.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/tests/bfgs.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/tests/cheminfo.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/tests/dihedrals.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/tests/find_motif.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/tests/frequencies.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/tests/geometry.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/tests/modify_geom.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/tests/multimolecular.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/tests/resonance.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/tests/symmetry.py +0 -0
- {kinbot-2.2.2 → kinbot-2.2.3}/tests/test_kinbot.py +0 -0
|
@@ -1,6 +1,6 @@
|
|
|
1
1
|
Metadata-Version: 2.1
|
|
2
2
|
Name: kinbot
|
|
3
|
-
Version: 2.2.
|
|
3
|
+
Version: 2.2.3
|
|
4
4
|
Summary: Automated reaction kinetics for gas-phase species
|
|
5
5
|
Author-email: Judit Zádor <jzador@sandia.gov>, Ruben Van de Vijver <ruben.vandevijver@ugent.be>, Carles Martí <cmartia@sandia.gov>, Clément Soulié <csoulie@sandia.gov>, Amanda Dewyer <adewyer@sandia.gov>
|
|
6
6
|
Maintainer-email: Judit Zádor <jzador@sandia.gov>, Clément Soulié <csoulie@sandia.gov>
|
|
@@ -47,8 +47,8 @@ Classifier: Topic :: Scientific/Engineering :: Chemistry
|
|
|
47
47
|
Requires-Python: >=3.10
|
|
48
48
|
Description-Content-Type: text/markdown
|
|
49
49
|
License-File: LICENSE
|
|
50
|
-
Requires-Dist: numpy
|
|
51
|
-
Requires-Dist: ase
|
|
50
|
+
Requires-Dist: numpy>=1.17.0
|
|
51
|
+
Requires-Dist: ase>=3.19
|
|
52
52
|
Requires-Dist: networkx
|
|
53
53
|
Requires-Dist: rmsd>=1.5.1
|
|
54
54
|
Requires-Dist: sella
|
|
@@ -155,37 +155,36 @@ See [list](https://github.com/zadorlab/KinBot/wiki/KinBot-file-structure).
|
|
|
155
155
|
* Clément Soulié (csoulie@sandia.gov)
|
|
156
156
|
|
|
157
157
|
## Papers using KinBot
|
|
158
|
-
|
|
159
|
-
|
|
160
|
-
|
|
161
|
-
|
|
162
|
-
|
|
163
|
-
|
|
164
|
-
|
|
165
|
-
|
|
166
|
-
|
|
167
|
-
|
|
168
|
-
|
|
169
|
-
|
|
170
|
-
|
|
171
|
-
|
|
172
|
-
|
|
173
|
-
|
|
174
|
-
|
|
175
|
-
|
|
176
|
-
|
|
177
|
-
|
|
178
|
-
|
|
179
|
-
|
|
158
|
+
1. Almeida, T. G., Martí, C., Kurtén, T., Zádor, J., Johansen, S. L.: _Theoretical analysis of the OH-Initiated atmospheric oxidation reactions of imidazole_. Phys. Chem. Chem. Phys., **2024** 26 23570-23587. https://doi.org/10.1039/D4CP02103G
|
|
159
|
+
2. Yuan, E. C.-Y., Kumar, A., Guan, X., Hermes, E. D., Rosen, A. S., Zádor, J., Head-Gordon, T., Blau, S. M.: _Analytical ab initio Hessian from a Deep Learning Potential for Transition State Optimization_. Nat. Comm., **2024**
|
|
160
|
+
3. Doner, A. C., Zádor, J., Rotavera, B.: _Stereoisomer-dependent rate coefficients and reaction mechanisms of 2-ethyloxetanylperoxy radicals_. Proc. Combust. Inst., **2024**, 40, 105578. https://doi.org/10.1016/j.proci.2024.105578
|
|
161
|
+
4. Hansen, N. A, Price, T. D., Filardi, L. R., Gurses, S. M., Zhou, W., Hansen, N., Osborn, D. L. Zádor, J., Kronawitter, C. X.: _The photoionization of methoxymethanol: Fingerprinting a reactive C2 oxygenate in a complex reactive mixture_. J. Chem. Phys., **2024**, 160, 124306. https://doi.org/10.1063/5.0197827
|
|
162
|
+
5. Martí, C., Devereux, C., Najm, H. N., Zádor, J.: _Evaluation of rate coefficients in the gas-phase using a machine learned potential_. J. Phys. Chem. A, **2024**, 128, 1958–1971. https://doi.org/10.1021/acs.jpca.3c07872
|
|
163
|
+
6. Lang, J., Foley, C. D., Thawoos, S., Behzadfar, A., Liu, Y., Zádor, J., Suits, A. G.: _Reaction dynamics of S(3P) with 1,3-butadiene and isoprene: Crossed beam scattering, low temperature flow experiments, and high-level electronic structure calculations_. Farad. Discuss., **2024**, 251, 550-572. https://doi.org/10.1039/D4FD00009A
|
|
164
|
+
7. Wang, D., Tian, Z.-Y., Zheng, Z.-H., Li, W., Wu, L.-N., Kuang, J.-J., Yang, J.-Z.: _Experimental and modeling study of the n, n-dimethylformamide pyrolysis at atmospheric pressure_. Combust. Flame, **2024**, 260, 113240. https://doi.org/10.1016/j.combustflame.2023.113240
|
|
165
|
+
8. Doner, A. C., Zádor, J., Rotavera, B.: _Unimolecular reactions of 2,4-dimethyloxetanyl radicals._ J. Phys. Chem A, **2023**, 127, 2591–2600 https://doi.org/10.1021/acs.jpca.2c08290
|
|
166
|
+
9. Li, H., Lang, J., Foley, C. D., Zádor, J., Suits, A. G.: _Sulfur (3P) reaction with conjugated dienes gives cyclization to thiophenes under single collision conditions._ J. Phys. Chem. Letters, **2023**, 14, 7611–7617. https://doi.org/10.1021/acs.jpclett.3c01953
|
|
167
|
+
10. Martí, C., Michelsen, H. A., Najm, H. N., Zádor, J.: _Comprehensive kinetics on the C7H7 potential energy surface under combustion conditions._ J. Phys. Chem. A, **2023**, 127, 1941–1959. https://pubs.acs.org/doi/full/10.1021/acs.jpca.2c08035
|
|
168
|
+
11. Zádor, J, Martí, C., Van de Vijver, R., Johansen, S. L., Yang, Y., Michelsen, H. A., Najm, H. N.: _Automated reaction kinetics of gas-phase organic species over multiwell potential energy surfaces._ J. Phys. Chem. A, **2023**, 127, 565–588. https://doi.org/10.1021/acs.jpca.2c06558
|
|
169
|
+
12. Lockwood, K. S., Ahmed, S. F., Huq, N. A., Stutzman, S. C., Foust, T. D., Labbe, N. J.: _Advances in predictive chemistry enable a multi-scale rational design approach for biofuels with advantaged properties_ Sustainable Energy Fuels, **2022**, 6, 5371-5383. https://doi.org/10.1039/D2SE00773H
|
|
170
|
+
13. Takahashi, L., Yoshida, S., Fujima, J., Oikawa, H., Takahashi, K.: _Unveiling the reaction pathways of hydrocarbons via experiments, computations and data science._ Phys. Chem. Chem. Phys., **2022**, 24, 29841-29849. https://pubs.rsc.org/en/content/articlelanding/2022/CP/D2CP04499D
|
|
171
|
+
14. Doner, A. C., Zádor, J., Rotavera, B.: _Stereoisomer-dependent unimolecular kinetics of 2,4-dimethyloxetane peroxy radicals._ Faraday Discuss., **2022**, 238, 295-319. https://doi.org/10.1039/D2FD00029F
|
|
172
|
+
15. Ramasesha, K., Savee, J. D., Zádor, J., Osborn, D. L.: _A New Pathway for Intersystem Crossing: Unexpected Products in the O(3P) + Cyclopentene Reaction._ J. Phys. Chem. A, **2021**, 125 9785-9801. https://doi.org/10.1021/acs.jpca.1c05817
|
|
173
|
+
16. Rogers, C. O, Lockwood, K. S., Nguyen, Q. L. D., Labbe, N. J.: _Diol isomer revealed as a source of methyl ketene from propionic acid unimolecular decomposition._ Int. J. Chem. Kinet., **2021**, 53, 1272–1284. https://doi.org/10.1002/kin.21532
|
|
174
|
+
17. Lockwood, K. S., Labbe, N. J.: _Insights on keto-hydroperoxide formation from O2 addition to the beta-tetrahydrofuran radical._ Proceedings of the Combustion Institute, **2021**, 38, 1, 533. https://doi.org/10.1016/j.proci.2020.06.357
|
|
175
|
+
18. Sheps, L., Dewyer, A. L., Demireva, M., and Zádor, J.: _Quantitative Detection of Products and Radical Intermediates in Low-Temperature Oxidation of Cyclopentane._ J. Phys. Chem. A **2021**, 125, 20, 4467. https://doi.org/10.1021/acs.jpca.1c02001
|
|
176
|
+
19. Zhang, J., Vermeire, F., Van de Vijver, R., Herbinet, O.; Battin-Leclerc, F., Reyniers, M.-F., Van Geem, K. M.: _Detailed experimental and kinetic modeling study of 3-carene pyrolysis._ Int. J. Chem. Kinet., **2020**, 52, 785-795. https://doi.org/10.1002/kin.21400
|
|
177
|
+
20. Van de Vijver, R., Zádor, J.: _KinBot: Automated stationary point search on potential energy surfaces._ Computer Physics Communications, **2020**, 248, 106947. https://doi.org/10.1016/j.cpc.2019.106947
|
|
178
|
+
21. Joshi, S. P., Seal, P., Pekkanen, T. T., Timonen, R. S., Eskola, A. J.: _Direct Kinetic Measurements and Master Equation Modelling of the Unimolecular Decomposition of Resonantly-Stabilized CH2CHCHC(O)OCH3 Radical and an Upper Limit Determination for CH2CHCHC(O)OCH3+O2 Reaction._ Z. Phys. Chem., **2020**, 234, 1251. https://doi.org/10.1515/zpch-2020-1612
|
|
180
179
|
|
|
181
180
|
Older Version of KinBot:
|
|
182
|
-
|
|
183
|
-
|
|
184
|
-
|
|
185
|
-
|
|
186
|
-
|
|
187
|
-
|
|
188
|
-
|
|
181
|
+
1. Van de Vijver, R., Van Geem, K. M., Marin, G. B., Zádor, J.: _Decomposition and isomerization of 1-pentanol radicals and the pyrolysis of 1-pentanol._ Combustion and Flame, **2018,** 196, 500. https://doi.org/10.1016/j.combustflame.2018.05.011
|
|
182
|
+
2. Grambow, C. A., Jamal, A., Li, Y.-P., Green, W. H., Zádor, J., Suleimanov, Y. V.: _Unimolecular reaction pathways of a g-ketohydroperoxide from combined application of automated reaction discovery methods._ J. Am. Chem. Soc., 2018, 140, 1035. https://doi.org/10.1021/jacs.7b11009
|
|
183
|
+
3. Rotavera, B., Savee, J. D., Antonov, I. O., Caravan, R. L., Sheps, L., Osborn, D. L., Zádor, J., Taatjes, C. A.: _Influence of oxygenation in cyclic hydrocarbons on chain-termination reactions from R + O2: tetrahydropyran and cyclohexane._ Proceedings of the Combustion Institute, **2017,** 36, 597. https://doi.org/10.1016/j.proci.2016.05.020
|
|
184
|
+
4. Antonov, I. O., Zádor, J., Rotavera, B., Papajak, E., Osborn, D. L., Taatjes, C. A., Sheps, L.: _Pressure-Dependent Competition among Reaction Pathways from First- and Second-O2 Additions in the Low-Temperature Oxidation of Tetrahydrofuran._ J. Phys. Chem. A, **2016,** 120 6582. https://doi.org/10.1021/acs.jpca.6b05411
|
|
185
|
+
5. Antonov, I. O., Kwok, J., Zádor, J., Sheps, L.: OH + 2-butene: A combined experimental and theoretical study in the 300-800 K temperature range. J. Phys. Chem. A, **2015,** 119, 7742. https://doi.org/10.1021/acs.jpca.5b01012
|
|
186
|
+
6. Zádor, J., Miller, J.A.: _Adventures on the C3H5O potential energy surface: OH + propyne, OH + allene and related reactions._ Proceedings of the Combustion Institute, **2015,** 35, 181. https://doi.org/10.1016/j.proci.2014.05.103
|
|
187
|
+
7. Rotavera, B., Zádor, J., Welz, O., Sheps, L., Scheer, A.M., Savee, J.D., Ali, M.A., Lee, T.S., Simmons, B.A., Osborn, D.L., Violi, A., Taatjes, C.A.: _Photoionization mass spectrometric measurements of initial reaction pathways in low-temperature oxidation of 2,5-dimethylhexane._ J. Phys. Chem. A, **2014,** 44, 10188. https://doi.org/10.1021/jp507811d
|
|
189
188
|
|
|
190
189
|
## Acknowledgement
|
|
191
190
|
This research was supported by the Exascale Computing Project (ECP), Project Number: 17-SC-20-SC, a collaborative effort of two DOE organizations, the Office of Science and the National Nuclear Security Administration, responsible for the planning and preparation of a capable exascale ecosystem including software, applications, hardware, advanced system engineering, and early test bed platforms to support the nation's exascale computing imperative. RVdV was also supported by the AITSTME project as part of the Predictive Theory and Modeling component of the Materials Genome Initiative. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology and Engineering Solutions of Sandia, LLC., a wholly owned subsidiary of Honeywell International, Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.
|
|
@@ -1,6 +1,6 @@
|
|
|
1
1
|
Metadata-Version: 2.1
|
|
2
2
|
Name: kinbot
|
|
3
|
-
Version: 2.2.
|
|
3
|
+
Version: 2.2.3
|
|
4
4
|
Summary: Automated reaction kinetics for gas-phase species
|
|
5
5
|
Author-email: Judit Zádor <jzador@sandia.gov>, Ruben Van de Vijver <ruben.vandevijver@ugent.be>, Carles Martí <cmartia@sandia.gov>, Clément Soulié <csoulie@sandia.gov>, Amanda Dewyer <adewyer@sandia.gov>
|
|
6
6
|
Maintainer-email: Judit Zádor <jzador@sandia.gov>, Clément Soulié <csoulie@sandia.gov>
|
|
@@ -47,8 +47,8 @@ Classifier: Topic :: Scientific/Engineering :: Chemistry
|
|
|
47
47
|
Requires-Python: >=3.10
|
|
48
48
|
Description-Content-Type: text/markdown
|
|
49
49
|
License-File: LICENSE
|
|
50
|
-
Requires-Dist: numpy
|
|
51
|
-
Requires-Dist: ase
|
|
50
|
+
Requires-Dist: numpy>=1.17.0
|
|
51
|
+
Requires-Dist: ase>=3.19
|
|
52
52
|
Requires-Dist: networkx
|
|
53
53
|
Requires-Dist: rmsd>=1.5.1
|
|
54
54
|
Requires-Dist: sella
|
|
@@ -155,37 +155,36 @@ See [list](https://github.com/zadorlab/KinBot/wiki/KinBot-file-structure).
|
|
|
155
155
|
* Clément Soulié (csoulie@sandia.gov)
|
|
156
156
|
|
|
157
157
|
## Papers using KinBot
|
|
158
|
-
|
|
159
|
-
|
|
160
|
-
|
|
161
|
-
|
|
162
|
-
|
|
163
|
-
|
|
164
|
-
|
|
165
|
-
|
|
166
|
-
|
|
167
|
-
|
|
168
|
-
|
|
169
|
-
|
|
170
|
-
|
|
171
|
-
|
|
172
|
-
|
|
173
|
-
|
|
174
|
-
|
|
175
|
-
|
|
176
|
-
|
|
177
|
-
|
|
178
|
-
|
|
179
|
-
|
|
158
|
+
1. Almeida, T. G., Martí, C., Kurtén, T., Zádor, J., Johansen, S. L.: _Theoretical analysis of the OH-Initiated atmospheric oxidation reactions of imidazole_. Phys. Chem. Chem. Phys., **2024** 26 23570-23587. https://doi.org/10.1039/D4CP02103G
|
|
159
|
+
2. Yuan, E. C.-Y., Kumar, A., Guan, X., Hermes, E. D., Rosen, A. S., Zádor, J., Head-Gordon, T., Blau, S. M.: _Analytical ab initio Hessian from a Deep Learning Potential for Transition State Optimization_. Nat. Comm., **2024**
|
|
160
|
+
3. Doner, A. C., Zádor, J., Rotavera, B.: _Stereoisomer-dependent rate coefficients and reaction mechanisms of 2-ethyloxetanylperoxy radicals_. Proc. Combust. Inst., **2024**, 40, 105578. https://doi.org/10.1016/j.proci.2024.105578
|
|
161
|
+
4. Hansen, N. A, Price, T. D., Filardi, L. R., Gurses, S. M., Zhou, W., Hansen, N., Osborn, D. L. Zádor, J., Kronawitter, C. X.: _The photoionization of methoxymethanol: Fingerprinting a reactive C2 oxygenate in a complex reactive mixture_. J. Chem. Phys., **2024**, 160, 124306. https://doi.org/10.1063/5.0197827
|
|
162
|
+
5. Martí, C., Devereux, C., Najm, H. N., Zádor, J.: _Evaluation of rate coefficients in the gas-phase using a machine learned potential_. J. Phys. Chem. A, **2024**, 128, 1958–1971. https://doi.org/10.1021/acs.jpca.3c07872
|
|
163
|
+
6. Lang, J., Foley, C. D., Thawoos, S., Behzadfar, A., Liu, Y., Zádor, J., Suits, A. G.: _Reaction dynamics of S(3P) with 1,3-butadiene and isoprene: Crossed beam scattering, low temperature flow experiments, and high-level electronic structure calculations_. Farad. Discuss., **2024**, 251, 550-572. https://doi.org/10.1039/D4FD00009A
|
|
164
|
+
7. Wang, D., Tian, Z.-Y., Zheng, Z.-H., Li, W., Wu, L.-N., Kuang, J.-J., Yang, J.-Z.: _Experimental and modeling study of the n, n-dimethylformamide pyrolysis at atmospheric pressure_. Combust. Flame, **2024**, 260, 113240. https://doi.org/10.1016/j.combustflame.2023.113240
|
|
165
|
+
8. Doner, A. C., Zádor, J., Rotavera, B.: _Unimolecular reactions of 2,4-dimethyloxetanyl radicals._ J. Phys. Chem A, **2023**, 127, 2591–2600 https://doi.org/10.1021/acs.jpca.2c08290
|
|
166
|
+
9. Li, H., Lang, J., Foley, C. D., Zádor, J., Suits, A. G.: _Sulfur (3P) reaction with conjugated dienes gives cyclization to thiophenes under single collision conditions._ J. Phys. Chem. Letters, **2023**, 14, 7611–7617. https://doi.org/10.1021/acs.jpclett.3c01953
|
|
167
|
+
10. Martí, C., Michelsen, H. A., Najm, H. N., Zádor, J.: _Comprehensive kinetics on the C7H7 potential energy surface under combustion conditions._ J. Phys. Chem. A, **2023**, 127, 1941–1959. https://pubs.acs.org/doi/full/10.1021/acs.jpca.2c08035
|
|
168
|
+
11. Zádor, J, Martí, C., Van de Vijver, R., Johansen, S. L., Yang, Y., Michelsen, H. A., Najm, H. N.: _Automated reaction kinetics of gas-phase organic species over multiwell potential energy surfaces._ J. Phys. Chem. A, **2023**, 127, 565–588. https://doi.org/10.1021/acs.jpca.2c06558
|
|
169
|
+
12. Lockwood, K. S., Ahmed, S. F., Huq, N. A., Stutzman, S. C., Foust, T. D., Labbe, N. J.: _Advances in predictive chemistry enable a multi-scale rational design approach for biofuels with advantaged properties_ Sustainable Energy Fuels, **2022**, 6, 5371-5383. https://doi.org/10.1039/D2SE00773H
|
|
170
|
+
13. Takahashi, L., Yoshida, S., Fujima, J., Oikawa, H., Takahashi, K.: _Unveiling the reaction pathways of hydrocarbons via experiments, computations and data science._ Phys. Chem. Chem. Phys., **2022**, 24, 29841-29849. https://pubs.rsc.org/en/content/articlelanding/2022/CP/D2CP04499D
|
|
171
|
+
14. Doner, A. C., Zádor, J., Rotavera, B.: _Stereoisomer-dependent unimolecular kinetics of 2,4-dimethyloxetane peroxy radicals._ Faraday Discuss., **2022**, 238, 295-319. https://doi.org/10.1039/D2FD00029F
|
|
172
|
+
15. Ramasesha, K., Savee, J. D., Zádor, J., Osborn, D. L.: _A New Pathway for Intersystem Crossing: Unexpected Products in the O(3P) + Cyclopentene Reaction._ J. Phys. Chem. A, **2021**, 125 9785-9801. https://doi.org/10.1021/acs.jpca.1c05817
|
|
173
|
+
16. Rogers, C. O, Lockwood, K. S., Nguyen, Q. L. D., Labbe, N. J.: _Diol isomer revealed as a source of methyl ketene from propionic acid unimolecular decomposition._ Int. J. Chem. Kinet., **2021**, 53, 1272–1284. https://doi.org/10.1002/kin.21532
|
|
174
|
+
17. Lockwood, K. S., Labbe, N. J.: _Insights on keto-hydroperoxide formation from O2 addition to the beta-tetrahydrofuran radical._ Proceedings of the Combustion Institute, **2021**, 38, 1, 533. https://doi.org/10.1016/j.proci.2020.06.357
|
|
175
|
+
18. Sheps, L., Dewyer, A. L., Demireva, M., and Zádor, J.: _Quantitative Detection of Products and Radical Intermediates in Low-Temperature Oxidation of Cyclopentane._ J. Phys. Chem. A **2021**, 125, 20, 4467. https://doi.org/10.1021/acs.jpca.1c02001
|
|
176
|
+
19. Zhang, J., Vermeire, F., Van de Vijver, R., Herbinet, O.; Battin-Leclerc, F., Reyniers, M.-F., Van Geem, K. M.: _Detailed experimental and kinetic modeling study of 3-carene pyrolysis._ Int. J. Chem. Kinet., **2020**, 52, 785-795. https://doi.org/10.1002/kin.21400
|
|
177
|
+
20. Van de Vijver, R., Zádor, J.: _KinBot: Automated stationary point search on potential energy surfaces._ Computer Physics Communications, **2020**, 248, 106947. https://doi.org/10.1016/j.cpc.2019.106947
|
|
178
|
+
21. Joshi, S. P., Seal, P., Pekkanen, T. T., Timonen, R. S., Eskola, A. J.: _Direct Kinetic Measurements and Master Equation Modelling of the Unimolecular Decomposition of Resonantly-Stabilized CH2CHCHC(O)OCH3 Radical and an Upper Limit Determination for CH2CHCHC(O)OCH3+O2 Reaction._ Z. Phys. Chem., **2020**, 234, 1251. https://doi.org/10.1515/zpch-2020-1612
|
|
180
179
|
|
|
181
180
|
Older Version of KinBot:
|
|
182
|
-
|
|
183
|
-
|
|
184
|
-
|
|
185
|
-
|
|
186
|
-
|
|
187
|
-
|
|
188
|
-
|
|
181
|
+
1. Van de Vijver, R., Van Geem, K. M., Marin, G. B., Zádor, J.: _Decomposition and isomerization of 1-pentanol radicals and the pyrolysis of 1-pentanol._ Combustion and Flame, **2018,** 196, 500. https://doi.org/10.1016/j.combustflame.2018.05.011
|
|
182
|
+
2. Grambow, C. A., Jamal, A., Li, Y.-P., Green, W. H., Zádor, J., Suleimanov, Y. V.: _Unimolecular reaction pathways of a g-ketohydroperoxide from combined application of automated reaction discovery methods._ J. Am. Chem. Soc., 2018, 140, 1035. https://doi.org/10.1021/jacs.7b11009
|
|
183
|
+
3. Rotavera, B., Savee, J. D., Antonov, I. O., Caravan, R. L., Sheps, L., Osborn, D. L., Zádor, J., Taatjes, C. A.: _Influence of oxygenation in cyclic hydrocarbons on chain-termination reactions from R + O2: tetrahydropyran and cyclohexane._ Proceedings of the Combustion Institute, **2017,** 36, 597. https://doi.org/10.1016/j.proci.2016.05.020
|
|
184
|
+
4. Antonov, I. O., Zádor, J., Rotavera, B., Papajak, E., Osborn, D. L., Taatjes, C. A., Sheps, L.: _Pressure-Dependent Competition among Reaction Pathways from First- and Second-O2 Additions in the Low-Temperature Oxidation of Tetrahydrofuran._ J. Phys. Chem. A, **2016,** 120 6582. https://doi.org/10.1021/acs.jpca.6b05411
|
|
185
|
+
5. Antonov, I. O., Kwok, J., Zádor, J., Sheps, L.: OH + 2-butene: A combined experimental and theoretical study in the 300-800 K temperature range. J. Phys. Chem. A, **2015,** 119, 7742. https://doi.org/10.1021/acs.jpca.5b01012
|
|
186
|
+
6. Zádor, J., Miller, J.A.: _Adventures on the C3H5O potential energy surface: OH + propyne, OH + allene and related reactions._ Proceedings of the Combustion Institute, **2015,** 35, 181. https://doi.org/10.1016/j.proci.2014.05.103
|
|
187
|
+
7. Rotavera, B., Zádor, J., Welz, O., Sheps, L., Scheer, A.M., Savee, J.D., Ali, M.A., Lee, T.S., Simmons, B.A., Osborn, D.L., Violi, A., Taatjes, C.A.: _Photoionization mass spectrometric measurements of initial reaction pathways in low-temperature oxidation of 2,5-dimethylhexane._ J. Phys. Chem. A, **2014,** 44, 10188. https://doi.org/10.1021/jp507811d
|
|
189
188
|
|
|
190
189
|
## Acknowledgement
|
|
191
190
|
This research was supported by the Exascale Computing Project (ECP), Project Number: 17-SC-20-SC, a collaborative effort of two DOE organizations, the Office of Science and the National Nuclear Security Administration, responsible for the planning and preparation of a capable exascale ecosystem including software, applications, hardware, advanced system engineering, and early test bed platforms to support the nation's exascale computing imperative. RVdV was also supported by the AITSTME project as part of the Predictive Theory and Modeling component of the Materials Genome Initiative. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology and Engineering Solutions of Sandia, LLC., a wholly owned subsidiary of Honeywell International, Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.
|
kinbot-2.2.3/README.md
ADDED
|
@@ -0,0 +1,131 @@
|
|
|
1
|
+
[](https://gitter.im/zadorlab/KinBot)
|
|
2
|
+
|
|
3
|
+
# KinBot: Automated Reaction Kinetics of Gas-Phase Organic Species over Multiwell Potential Energy Surfaces
|
|
4
|
+
|
|
5
|
+
<p>
|
|
6
|
+
<img src="https://raw.githubusercontent.com/zadorlab/KinBot/master/graphics/kinbot_logo_V2.png" width="220" height="240" />
|
|
7
|
+
</p>
|
|
8
|
+
|
|
9
|
+
## Description
|
|
10
|
+
This repository contains the KinBot code version 2.2.1,
|
|
11
|
+
a tool for automatically searching for reactions on the potential energy surface.
|
|
12
|
+
|
|
13
|
+
If you are using this tool in scientific publications, please reference the following publications:
|
|
14
|
+
|
|
15
|
+
* Ruben Van de Vijver, Judit Zádor: KinBot: _Automated stationary point search on potential energy surfaces_, Comp. Phys. Comm., **2019**, 248, 106947. https://doi.org/10.1016/j.cpc.2019.106947
|
|
16
|
+
```
|
|
17
|
+
@article{Vijver2020,
|
|
18
|
+
author = {Van de Vijver, Ruben and Z\'ador, Judit},
|
|
19
|
+
title = {KinBot: Automated stationary point search on potential energy surfaces},
|
|
20
|
+
journal = {Comput. Phys. Commun.},
|
|
21
|
+
volume = {248},
|
|
22
|
+
pages = {106947},
|
|
23
|
+
year = {2020},
|
|
24
|
+
type = {Journal Article}
|
|
25
|
+
}
|
|
26
|
+
```
|
|
27
|
+
* Judit Zádor, Carles Martí, Ruben Van de Vijver, Sommer L. Johansen, Yoona Yang, Hope A. Michelsen, Habib N. Najm: _Automated reaction kinetics of gas-phase organic species over multiwell potential energy surfaces_, J. Phys. Chem. A, **2023**, 127, 565–588. https://doi.org/10.1021/acs.jpca.2c06558
|
|
28
|
+
|
|
29
|
+
```
|
|
30
|
+
@article{Zador2022,
|
|
31
|
+
author = {Z\'ador, Judit and Mart\'i, Carles and Van de Vijver, Ruben and Johansen, Sommer L. and Yang, Yoona and Michelsen, Hope A. and Najm, Habib N.},
|
|
32
|
+
title = {Automated reaction kinetics of gas-phase organic species over multiwell potential energy surfaces},
|
|
33
|
+
journal = {J. Phys. Chem. A},
|
|
34
|
+
volume = {127},
|
|
35
|
+
pages = {565-588},
|
|
36
|
+
year = {2023},
|
|
37
|
+
type = {Journal Article}
|
|
38
|
+
}
|
|
39
|
+
```
|
|
40
|
+
|
|
41
|
+
We appreciate if you send us the DOI of your published paper that used KinBot, so we can feature it here below.
|
|
42
|
+
|
|
43
|
+
## How to Install
|
|
44
|
+
|
|
45
|
+
KinBot can be installed both in three different ways, from the PyPI index (`pip install`), from the conda-forge repo (`conda install`) or by cloning this github repo and then install it locally.
|
|
46
|
+
|
|
47
|
+
### PyPI
|
|
48
|
+
|
|
49
|
+
pip install kinbot
|
|
50
|
+
|
|
51
|
+
> **Note**
|
|
52
|
+
> KinBot only works with Python >= 3.10.
|
|
53
|
+
|
|
54
|
+
### conda-forge
|
|
55
|
+
|
|
56
|
+
conda install -c conda-forge kinbot
|
|
57
|
+
|
|
58
|
+
### From Github
|
|
59
|
+
|
|
60
|
+
If you want to have the very last version of KinBot without waiting for a
|
|
61
|
+
release or you want to modify KinBot acccording to your needs you can clone the project
|
|
62
|
+
from github:
|
|
63
|
+
|
|
64
|
+
git clone git@github.com:zadorlab/KinBot.git
|
|
65
|
+
|
|
66
|
+
and then, from within the KinBot directory produced after cloning, type:
|
|
67
|
+
|
|
68
|
+
pip install -e .
|
|
69
|
+
|
|
70
|
+
> **Note**
|
|
71
|
+
> If you want to modify KinBot yourself it's better to fork the project
|
|
72
|
+
> into your own repository and then clone it.
|
|
73
|
+
|
|
74
|
+
## How to Run
|
|
75
|
+
To run a single-well exploration of KinBot, make an input file (e.g. input.json) and run:
|
|
76
|
+
|
|
77
|
+
kinbot input.json
|
|
78
|
+
|
|
79
|
+
To run a full PES search, make an input file (e.g. input.json) and run:
|
|
80
|
+
|
|
81
|
+
pes input.json
|
|
82
|
+
|
|
83
|
+
You can find additional command line arguments in the manual.
|
|
84
|
+
|
|
85
|
+
## Documentation
|
|
86
|
+
See the [wiki](https://github.com/zadorlab/KinBot/wiki) for keywords, and our [tutorial](https://hackmd.io/@jzador/kinbot_workshop_2023#/) for a more hands-on introduction to the code.
|
|
87
|
+
|
|
88
|
+
## List of files in this project
|
|
89
|
+
See [list](https://github.com/zadorlab/KinBot/wiki/KinBot-file-structure).
|
|
90
|
+
|
|
91
|
+
## Authors
|
|
92
|
+
* Judit Zádor (jzador@sandia.gov)
|
|
93
|
+
* Ruben Van de Vijver
|
|
94
|
+
* Amanda Dewyer
|
|
95
|
+
* Carles Martí
|
|
96
|
+
* Clément Soulié (csoulie@sandia.gov)
|
|
97
|
+
|
|
98
|
+
## Papers using KinBot
|
|
99
|
+
1. Almeida, T. G., Martí, C., Kurtén, T., Zádor, J., Johansen, S. L.: _Theoretical analysis of the OH-Initiated atmospheric oxidation reactions of imidazole_. Phys. Chem. Chem. Phys., **2024** 26 23570-23587. https://doi.org/10.1039/D4CP02103G
|
|
100
|
+
2. Yuan, E. C.-Y., Kumar, A., Guan, X., Hermes, E. D., Rosen, A. S., Zádor, J., Head-Gordon, T., Blau, S. M.: _Analytical ab initio Hessian from a Deep Learning Potential for Transition State Optimization_. Nat. Comm., **2024**
|
|
101
|
+
3. Doner, A. C., Zádor, J., Rotavera, B.: _Stereoisomer-dependent rate coefficients and reaction mechanisms of 2-ethyloxetanylperoxy radicals_. Proc. Combust. Inst., **2024**, 40, 105578. https://doi.org/10.1016/j.proci.2024.105578
|
|
102
|
+
4. Hansen, N. A, Price, T. D., Filardi, L. R., Gurses, S. M., Zhou, W., Hansen, N., Osborn, D. L. Zádor, J., Kronawitter, C. X.: _The photoionization of methoxymethanol: Fingerprinting a reactive C2 oxygenate in a complex reactive mixture_. J. Chem. Phys., **2024**, 160, 124306. https://doi.org/10.1063/5.0197827
|
|
103
|
+
5. Martí, C., Devereux, C., Najm, H. N., Zádor, J.: _Evaluation of rate coefficients in the gas-phase using a machine learned potential_. J. Phys. Chem. A, **2024**, 128, 1958–1971. https://doi.org/10.1021/acs.jpca.3c07872
|
|
104
|
+
6. Lang, J., Foley, C. D., Thawoos, S., Behzadfar, A., Liu, Y., Zádor, J., Suits, A. G.: _Reaction dynamics of S(3P) with 1,3-butadiene and isoprene: Crossed beam scattering, low temperature flow experiments, and high-level electronic structure calculations_. Farad. Discuss., **2024**, 251, 550-572. https://doi.org/10.1039/D4FD00009A
|
|
105
|
+
7. Wang, D., Tian, Z.-Y., Zheng, Z.-H., Li, W., Wu, L.-N., Kuang, J.-J., Yang, J.-Z.: _Experimental and modeling study of the n, n-dimethylformamide pyrolysis at atmospheric pressure_. Combust. Flame, **2024**, 260, 113240. https://doi.org/10.1016/j.combustflame.2023.113240
|
|
106
|
+
8. Doner, A. C., Zádor, J., Rotavera, B.: _Unimolecular reactions of 2,4-dimethyloxetanyl radicals._ J. Phys. Chem A, **2023**, 127, 2591–2600 https://doi.org/10.1021/acs.jpca.2c08290
|
|
107
|
+
9. Li, H., Lang, J., Foley, C. D., Zádor, J., Suits, A. G.: _Sulfur (3P) reaction with conjugated dienes gives cyclization to thiophenes under single collision conditions._ J. Phys. Chem. Letters, **2023**, 14, 7611–7617. https://doi.org/10.1021/acs.jpclett.3c01953
|
|
108
|
+
10. Martí, C., Michelsen, H. A., Najm, H. N., Zádor, J.: _Comprehensive kinetics on the C7H7 potential energy surface under combustion conditions._ J. Phys. Chem. A, **2023**, 127, 1941–1959. https://pubs.acs.org/doi/full/10.1021/acs.jpca.2c08035
|
|
109
|
+
11. Zádor, J, Martí, C., Van de Vijver, R., Johansen, S. L., Yang, Y., Michelsen, H. A., Najm, H. N.: _Automated reaction kinetics of gas-phase organic species over multiwell potential energy surfaces._ J. Phys. Chem. A, **2023**, 127, 565–588. https://doi.org/10.1021/acs.jpca.2c06558
|
|
110
|
+
12. Lockwood, K. S., Ahmed, S. F., Huq, N. A., Stutzman, S. C., Foust, T. D., Labbe, N. J.: _Advances in predictive chemistry enable a multi-scale rational design approach for biofuels with advantaged properties_ Sustainable Energy Fuels, **2022**, 6, 5371-5383. https://doi.org/10.1039/D2SE00773H
|
|
111
|
+
13. Takahashi, L., Yoshida, S., Fujima, J., Oikawa, H., Takahashi, K.: _Unveiling the reaction pathways of hydrocarbons via experiments, computations and data science._ Phys. Chem. Chem. Phys., **2022**, 24, 29841-29849. https://pubs.rsc.org/en/content/articlelanding/2022/CP/D2CP04499D
|
|
112
|
+
14. Doner, A. C., Zádor, J., Rotavera, B.: _Stereoisomer-dependent unimolecular kinetics of 2,4-dimethyloxetane peroxy radicals._ Faraday Discuss., **2022**, 238, 295-319. https://doi.org/10.1039/D2FD00029F
|
|
113
|
+
15. Ramasesha, K., Savee, J. D., Zádor, J., Osborn, D. L.: _A New Pathway for Intersystem Crossing: Unexpected Products in the O(3P) + Cyclopentene Reaction._ J. Phys. Chem. A, **2021**, 125 9785-9801. https://doi.org/10.1021/acs.jpca.1c05817
|
|
114
|
+
16. Rogers, C. O, Lockwood, K. S., Nguyen, Q. L. D., Labbe, N. J.: _Diol isomer revealed as a source of methyl ketene from propionic acid unimolecular decomposition._ Int. J. Chem. Kinet., **2021**, 53, 1272–1284. https://doi.org/10.1002/kin.21532
|
|
115
|
+
17. Lockwood, K. S., Labbe, N. J.: _Insights on keto-hydroperoxide formation from O2 addition to the beta-tetrahydrofuran radical._ Proceedings of the Combustion Institute, **2021**, 38, 1, 533. https://doi.org/10.1016/j.proci.2020.06.357
|
|
116
|
+
18. Sheps, L., Dewyer, A. L., Demireva, M., and Zádor, J.: _Quantitative Detection of Products and Radical Intermediates in Low-Temperature Oxidation of Cyclopentane._ J. Phys. Chem. A **2021**, 125, 20, 4467. https://doi.org/10.1021/acs.jpca.1c02001
|
|
117
|
+
19. Zhang, J., Vermeire, F., Van de Vijver, R., Herbinet, O.; Battin-Leclerc, F., Reyniers, M.-F., Van Geem, K. M.: _Detailed experimental and kinetic modeling study of 3-carene pyrolysis._ Int. J. Chem. Kinet., **2020**, 52, 785-795. https://doi.org/10.1002/kin.21400
|
|
118
|
+
20. Van de Vijver, R., Zádor, J.: _KinBot: Automated stationary point search on potential energy surfaces._ Computer Physics Communications, **2020**, 248, 106947. https://doi.org/10.1016/j.cpc.2019.106947
|
|
119
|
+
21. Joshi, S. P., Seal, P., Pekkanen, T. T., Timonen, R. S., Eskola, A. J.: _Direct Kinetic Measurements and Master Equation Modelling of the Unimolecular Decomposition of Resonantly-Stabilized CH2CHCHC(O)OCH3 Radical and an Upper Limit Determination for CH2CHCHC(O)OCH3+O2 Reaction._ Z. Phys. Chem., **2020**, 234, 1251. https://doi.org/10.1515/zpch-2020-1612
|
|
120
|
+
|
|
121
|
+
Older Version of KinBot:
|
|
122
|
+
1. Van de Vijver, R., Van Geem, K. M., Marin, G. B., Zádor, J.: _Decomposition and isomerization of 1-pentanol radicals and the pyrolysis of 1-pentanol._ Combustion and Flame, **2018,** 196, 500. https://doi.org/10.1016/j.combustflame.2018.05.011
|
|
123
|
+
2. Grambow, C. A., Jamal, A., Li, Y.-P., Green, W. H., Zádor, J., Suleimanov, Y. V.: _Unimolecular reaction pathways of a g-ketohydroperoxide from combined application of automated reaction discovery methods._ J. Am. Chem. Soc., 2018, 140, 1035. https://doi.org/10.1021/jacs.7b11009
|
|
124
|
+
3. Rotavera, B., Savee, J. D., Antonov, I. O., Caravan, R. L., Sheps, L., Osborn, D. L., Zádor, J., Taatjes, C. A.: _Influence of oxygenation in cyclic hydrocarbons on chain-termination reactions from R + O2: tetrahydropyran and cyclohexane._ Proceedings of the Combustion Institute, **2017,** 36, 597. https://doi.org/10.1016/j.proci.2016.05.020
|
|
125
|
+
4. Antonov, I. O., Zádor, J., Rotavera, B., Papajak, E., Osborn, D. L., Taatjes, C. A., Sheps, L.: _Pressure-Dependent Competition among Reaction Pathways from First- and Second-O2 Additions in the Low-Temperature Oxidation of Tetrahydrofuran._ J. Phys. Chem. A, **2016,** 120 6582. https://doi.org/10.1021/acs.jpca.6b05411
|
|
126
|
+
5. Antonov, I. O., Kwok, J., Zádor, J., Sheps, L.: OH + 2-butene: A combined experimental and theoretical study in the 300-800 K temperature range. J. Phys. Chem. A, **2015,** 119, 7742. https://doi.org/10.1021/acs.jpca.5b01012
|
|
127
|
+
6. Zádor, J., Miller, J.A.: _Adventures on the C3H5O potential energy surface: OH + propyne, OH + allene and related reactions._ Proceedings of the Combustion Institute, **2015,** 35, 181. https://doi.org/10.1016/j.proci.2014.05.103
|
|
128
|
+
7. Rotavera, B., Zádor, J., Welz, O., Sheps, L., Scheer, A.M., Savee, J.D., Ali, M.A., Lee, T.S., Simmons, B.A., Osborn, D.L., Violi, A., Taatjes, C.A.: _Photoionization mass spectrometric measurements of initial reaction pathways in low-temperature oxidation of 2,5-dimethylhexane._ J. Phys. Chem. A, **2014,** 44, 10188. https://doi.org/10.1021/jp507811d
|
|
129
|
+
|
|
130
|
+
## Acknowledgement
|
|
131
|
+
This research was supported by the Exascale Computing Project (ECP), Project Number: 17-SC-20-SC, a collaborative effort of two DOE organizations, the Office of Science and the National Nuclear Security Administration, responsible for the planning and preparation of a capable exascale ecosystem including software, applications, hardware, advanced system engineering, and early test bed platforms to support the nation's exascale computing imperative. RVdV was also supported by the AITSTME project as part of the Predictive Theory and Modeling component of the Materials Genome Initiative. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology and Engineering Solutions of Sandia, LLC., a wholly owned subsidiary of Honeywell International, Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.
|
|
@@ -7,13 +7,14 @@ def start_motif(motif, natom, bond, atom, allover, eqv):
|
|
|
7
7
|
chain = [-999] * natom
|
|
8
8
|
nsteps = -1
|
|
9
9
|
motifset = []
|
|
10
|
+
mask = [] # to mask new patterns that allow eqv atoms in some cases
|
|
10
11
|
find_motif(motif, visit, chain, nsteps, 0, -1,
|
|
11
|
-
motifset, allover, natom, bond, atom, eqv)
|
|
12
|
+
motifset, allover, natom, bond, atom, eqv, mask, True)
|
|
12
13
|
return motifset
|
|
13
14
|
|
|
14
15
|
|
|
15
16
|
def find_motif(motif, visit, chain, nsteps, current,
|
|
16
|
-
previous, motifset, allover, natom, bond, atom, eqv):
|
|
17
|
+
previous, motifset, allover, natom, bond, atom, eqv, mask, mask_current):
|
|
17
18
|
"""
|
|
18
19
|
This recursive function finds a specific motif in the structure.
|
|
19
20
|
FIXIT - the comments here
|
|
@@ -28,6 +29,7 @@ def find_motif(motif, visit, chain, nsteps, current,
|
|
|
28
29
|
allover < 0: all-over path-finding mode
|
|
29
30
|
allover = > 0: only atom 'allover' is used for the search
|
|
30
31
|
eqv: array for equvivalent atoms
|
|
32
|
+
mask: mask irregulara patterns that contain equivalent atoms
|
|
31
33
|
"""
|
|
32
34
|
if nsteps == -1:
|
|
33
35
|
nsteps = 0
|
|
@@ -37,13 +39,14 @@ def find_motif(motif, visit, chain, nsteps, current,
|
|
|
37
39
|
current = i
|
|
38
40
|
visit = [0] * natom
|
|
39
41
|
chain = [-999] * natom
|
|
42
|
+
mask_current = True
|
|
40
43
|
find_motif(motif, visit, chain, nsteps, current, previous,
|
|
41
|
-
motifset, allover, natom, bond, atom, eqv)
|
|
44
|
+
motifset, allover, natom, bond, atom, eqv, mask, mask_current)
|
|
42
45
|
else:
|
|
43
46
|
current = allover
|
|
44
47
|
visit = [0] * natom
|
|
45
48
|
find_motif(motif, visit, chain, nsteps, current, previous,
|
|
46
|
-
motifset, allover, natom, bond, atom, eqv)
|
|
49
|
+
motifset, allover, natom, bond, atom, eqv, mask, mask_current)
|
|
47
50
|
|
|
48
51
|
if nsteps > -1:
|
|
49
52
|
if nsteps > natom:
|
|
@@ -55,9 +58,12 @@ def find_motif(motif, visit, chain, nsteps, current,
|
|
|
55
58
|
if current in mylist:
|
|
56
59
|
eqv_list = mylist[:]
|
|
57
60
|
eqv_list.remove(current)
|
|
58
|
-
for m in motifset:
|
|
61
|
+
for mi, m in enumerate(motifset):
|
|
62
|
+
if mask[mi] == False:
|
|
63
|
+
continue
|
|
59
64
|
if m[nsteps] in eqv_list:
|
|
60
|
-
if any([True for i in eqv_list if i in chain]):
|
|
65
|
+
if any([True for i in eqv_list if i in chain[:nsteps]]):
|
|
66
|
+
mask_current = False
|
|
61
67
|
break
|
|
62
68
|
return 0
|
|
63
69
|
if visit[current] == 1:
|
|
@@ -77,6 +83,7 @@ def find_motif(motif, visit, chain, nsteps, current,
|
|
|
77
83
|
motifset.append(chain[:chain.index(-999)])
|
|
78
84
|
else:
|
|
79
85
|
motifset.append(chain[:])
|
|
86
|
+
mask.append(mask_current)
|
|
80
87
|
return 0
|
|
81
88
|
|
|
82
89
|
visit[current] = 1
|
|
@@ -86,7 +93,7 @@ def find_motif(motif, visit, chain, nsteps, current,
|
|
|
86
93
|
for i in range(natom):
|
|
87
94
|
current = i
|
|
88
95
|
find_motif(motif, visit, chain, nsteps, current, previous,
|
|
89
|
-
motifset, allover, natom, bond, atom, eqv)
|
|
96
|
+
motifset, allover, natom, bond, atom, eqv, mask, mask_current)
|
|
90
97
|
visit[current] = 0
|
|
91
98
|
|
|
92
99
|
if nsteps > 0:
|
|
@@ -1,4 +1,3 @@
|
|
|
1
|
-
import rmsd
|
|
2
1
|
import os
|
|
3
2
|
from typing import Any
|
|
4
3
|
from kinbot import kb_path
|
|
@@ -13,6 +12,7 @@ import logging
|
|
|
13
12
|
import copy
|
|
14
13
|
import math
|
|
15
14
|
from kinbot.constants import BOHRtoANGSTROM
|
|
15
|
+
from kinbot.utils import reorder_coord
|
|
16
16
|
|
|
17
17
|
logger = logging.getLogger('KinBot')
|
|
18
18
|
|
|
@@ -68,12 +68,14 @@ class Fragment(StationaryPoint):
|
|
|
68
68
|
self.geom = np.array(geom)
|
|
69
69
|
|
|
70
70
|
self.recentre()
|
|
71
|
-
|
|
71
|
+
self.frag_name: str
|
|
72
|
+
if len(self.geom) == 1:
|
|
73
|
+
self.frag_type = 'Monoatomic'
|
|
74
|
+
elif abs(np.prod(self.atoms.get_moments_of_inertia())) < 1e-6:
|
|
72
75
|
self.frag_type = 'Linear'
|
|
73
76
|
else:
|
|
74
77
|
self.frag_type = 'Nonlinear'
|
|
75
78
|
|
|
76
|
-
self.frag_name: str
|
|
77
79
|
Fragment.set_fragnames(self)
|
|
78
80
|
|
|
79
81
|
super(Fragment, self).__init__(name=self.frag_name,
|
|
@@ -550,6 +552,13 @@ class Fragment(StationaryPoint):
|
|
|
550
552
|
raise TypeError('HOMO mode selected for pivot point orientation, but cubefile is not available.')
|
|
551
553
|
step: NDArray[float32] = np.zeros(3, dtype=float32)
|
|
552
554
|
dim: NDArray[int16] = np.zeros(3, dtype=int16)
|
|
555
|
+
elements: list[str] = []
|
|
556
|
+
num2str: dict[str, str] = {
|
|
557
|
+
'1': 'H',
|
|
558
|
+
'6': 'C',
|
|
559
|
+
'7': 'N',
|
|
560
|
+
'8': 'O',
|
|
561
|
+
'16': 'S'}
|
|
553
562
|
# Saves info from the header: step size, origin and dimension
|
|
554
563
|
for ln, line in enumerate(cubefile):
|
|
555
564
|
if ln == 2:
|
|
@@ -573,8 +582,9 @@ class Fragment(StationaryPoint):
|
|
|
573
582
|
continue
|
|
574
583
|
# Save the cube geometry
|
|
575
584
|
elif ln > 5 and ln <= 5 + natm:
|
|
576
|
-
x, y, z = line.split()
|
|
585
|
+
num, _, x, y, z = line.split()
|
|
577
586
|
cube_geom[ln-6] = np.array([x, y, z], dtype=float32)
|
|
587
|
+
elements.append(num2str[num])
|
|
578
588
|
elif ln > 5 and len(line.split()) == 2:
|
|
579
589
|
start: int = ln + 1
|
|
580
590
|
break
|
|
@@ -590,16 +600,20 @@ class Fragment(StationaryPoint):
|
|
|
590
600
|
xsize = int(ysize*dim[1])
|
|
591
601
|
|
|
592
602
|
# Change the geom for the one in the cube to ensure same orientation
|
|
593
|
-
|
|
603
|
+
tmp_atm = Atoms(elements, positions=cube_geom)
|
|
604
|
+
tmp_stp = StationaryPoint.from_ase_atoms(tmp_atm)
|
|
605
|
+
tmp_stp.characterize()
|
|
606
|
+
reorder_coord(self, tmp_stp)
|
|
607
|
+
self.geom = tmp_stp.geom
|
|
594
608
|
self.recentre()
|
|
595
609
|
|
|
596
|
-
for eq in self.equiv:
|
|
597
|
-
|
|
598
|
-
|
|
599
|
-
|
|
610
|
+
# for eq in self.equiv:
|
|
611
|
+
# if index in eq:
|
|
612
|
+
# to_integrate: list[int] = eq
|
|
613
|
+
# break
|
|
600
614
|
|
|
601
615
|
# max_integral = -1
|
|
602
|
-
max_var: float = np.inf
|
|
616
|
+
# max_var: float = np.inf
|
|
603
617
|
# all_orient = []
|
|
604
618
|
|
|
605
619
|
|
|
@@ -17,7 +17,6 @@ from ase.calculators.singlepoint import SinglePointCalculator
|
|
|
17
17
|
from ase.calculators.emt import EMT
|
|
18
18
|
from ase.constraints import FixInternals
|
|
19
19
|
from ase.data import atomic_numbers, covalent_radii
|
|
20
|
-
from sella import Sella
|
|
21
20
|
|
|
22
21
|
from kinbot import bfgs
|
|
23
22
|
from kinbot import constants
|
|
@@ -203,7 +202,10 @@ def modify_coordinates(species, name, geom, changes, bond, write_files=0):
|
|
|
203
202
|
gs += '{}, {:.8f}, {:.8f}, {:.8f}, \n'.format(at, x, y, z)
|
|
204
203
|
logger.debug("For the following initial geometry:\n" + gs)
|
|
205
204
|
|
|
206
|
-
|
|
205
|
+
# step_tol=1E-7, grad_tol=1E-7, line_tol=1E-10,
|
|
206
|
+
# inhess=None, max_step=100, max_lin_step=1000,
|
|
207
|
+
# use_grad_tol=1, use_step_tol=1
|
|
208
|
+
opt = bfgs.BFGS(step_tol=1E-12, grad_tol=1E-12)
|
|
207
209
|
x_opt, x_i, g_i = opt.optimize(cost_fct, x0)
|
|
208
210
|
|
|
209
211
|
new_geom = np.reshape(x_opt, (species.natom, 3))
|
|
@@ -416,11 +416,13 @@ class Optimize:
|
|
|
416
416
|
else:
|
|
417
417
|
names.append(self.name + '_well')
|
|
418
418
|
names.append(self.name + '_well_high')
|
|
419
|
-
if self.par['high_level'] == 1
|
|
419
|
+
if self.par['high_level'] == 1 and \
|
|
420
|
+
'prod' not in self.species.name:
|
|
420
421
|
for count in range(self.species.hir.nrotation):
|
|
421
422
|
for rot_num in range(self.par['nrotation']):
|
|
422
423
|
names.append('hir/' + self.name + '_hir_' + str(count) + '_' + str(rot_num).zfill(2))
|
|
423
|
-
if self.par['conformer_search'] == 1
|
|
424
|
+
if self.par['conformer_search'] == 1 and \
|
|
425
|
+
'prod' not in self.species.name:
|
|
424
426
|
for count in range(self.species.confs.conf):
|
|
425
427
|
names.append('conf/' + self.name + '_' + str(count).zfill(zf))
|
|
426
428
|
for count in range(self.species.confs.cyc_conf):
|
|
@@ -206,7 +206,7 @@ class Parameters:
|
|
|
206
206
|
# Basis set to scan bonds in barrierless_saddle family
|
|
207
207
|
'barrierless_saddle_basis_high': '6-31G',
|
|
208
208
|
# for Gaussian, request CalcAll for TS optimization
|
|
209
|
-
'calcall_ts':
|
|
209
|
+
'calcall_ts': 1,
|
|
210
210
|
# Quantum chemistry method to use for high-level L2
|
|
211
211
|
'high_level_method': 'M062X',
|
|
212
212
|
# Basis set to use for high-level
|
|
@@ -1400,7 +1400,7 @@ def create_rotdpy_inputs(par, bless, vdW) -> None:
|
|
|
1400
1400
|
parent=str(reactant),
|
|
1401
1401
|
mult=pp_info['frags_mult'][frag_num]))
|
|
1402
1402
|
|
|
1403
|
-
fragnames = Fragment.get_fragnames()
|
|
1403
|
+
fragnames: list[str] = Fragment.get_fragnames()
|
|
1404
1404
|
|
|
1405
1405
|
# Set the pivot points on each fragments and create the surfaces
|
|
1406
1406
|
surfaces: list[VRC_TST_Surface] = []
|
|
@@ -5,7 +5,6 @@ import numpy as np
|
|
|
5
5
|
from numpy.typing import NDArray
|
|
6
6
|
from ase import Atoms
|
|
7
7
|
|
|
8
|
-
from kinbot import constants
|
|
9
8
|
from kinbot.fragments import Fragment
|
|
10
9
|
from kinbot.vrc_tst_surfaces import VRC_TST_Surface
|
|
11
10
|
|
|
@@ -220,6 +219,7 @@ def create_surface(dist,
|
|
|
220
219
|
info[findex] += 'COM'
|
|
221
220
|
weights[findex].append(1)
|
|
222
221
|
reac_weights[findex].append(1)
|
|
222
|
+
all_pp_dists[findex].append(0.0)
|
|
223
223
|
continue
|
|
224
224
|
elif pps_dists is None:
|
|
225
225
|
info[findex] += 'on atom'
|