A method implemented in parameterize is inspired by GAAMP 1. A detailed description of our method is be published 2.

General procedure

  1. Get initial parameters (from GAFF or CGenFF)

  2. Minimize the geometry of molecule with QM

  3. Get atomic charges:
    1. Compute ESP with QM

    2. Fit the atomic charges to reproduce the QM results

  4. Get dihedral angle parameters:
    1. Detect rotatable dihedral angles

    2. Scan the dihedral angles with QM

    3. Fit the parameters to reproduce the QM results

  5. Check parameterization quality:
    1. Plot rotamer QM and MM energies


L. Huang and B. Roux, Automated Force Field Parameterization for Nonpolarizable and Polarizable Atomic Models Based on Ab Initio Target Data, J. Chem. Theory Comput., 2013, 9, 8, 3543-3556. DOI: 10.1021/ct4003477


R. Galvelis, S. Doerr, J. M. Damas, M. J. Harvey, and G. De Fabritiis, A Scalable Molecular Force Field Parameterization Method Based on Density Functional Theory and Quantum-Level Machine Learning, J. Chem. Inf. Model. 2019, 59, 8, 3485-3493. DOI: 10.1021/acs.jcim.9b00439 arXiv: 1907.06952