quemb.molbe.autofrag.autogen

quemb.molbe.autofrag.autogen(mol, frozen_core=True, n_BE=2, write_geom=False, iao_valence_basis=None, iao_valence_only=False, print_frags=True)

Automatic molecular partitioning

Partitions a molecule into overlapping fragments as defined in BE atom-based fragmentations. It automatically detects branched chemical chains and ring systems and partitions accordingly. For efficiency, it only checks two atoms for connectivity (chemical bond) if they are within 3.5 Angstrom. This value is hardcoded as normdist. Two atoms are defined as bonded if they are within 1.8 Angstrom (1.2 for Hydrogen atom). This is also hardcoded as bond & hbond.

Parameters:

  • mol (Mole) – This is required for the options, ‘autogen’, and ‘chain’ as frag_type.

  • frozen_core (bool, optional) – Whether to invoke frozen core approximation. Defaults to True.

  • n_BE (int, optional) – Specifies the order of bootstrap calculation in the atom-based fragmentation, i.e. BE(n). Supported values are 1, 2, 3, and 4 Defaults to 2.

  • write_geom (bool, optional) – Whether to write a ‘fragment.xyz’ file which contains all the fragments in Cartesian coordinates. Defaults to False.

  • iao_valence_basis (str, optional) – Name of minimal basis set for IAO scheme. ‘sto-3g’ is sufficient for most cases. Defaults to None.

  • iao_valence_only (bool, optional) – If True, all calculations will be performed in the valence basis in the IAO partitioning. This is an experimental feature. Defaults to False.

  • print_frags (bool, optional) – Whether to print out the list of resulting fragments. Defaults to True.