GaffMolecularMechanicsCalculator.h

Go to the documentation of this file.

#ifndef GAFFMOLECULARMECHANICSCALCULATOR_H
#define GAFFMOLECULARMECHANICSCALCULATOR_H

#include "../Interactions/AnglesEvaluator.h"
#include "../Interactions/BondsEvaluator.h"
#include "../Interactions/DihedralsEvaluator.h"
#include "../Interactions/ElectrostaticEvaluator.h"
#include "../Interactions/LennardJonesEvaluator.h"
#include "../MolecularMechanicsCalculator.h"
#include "../Topology/IndexedStructuralTopology.h"
#include "GaffAtomTypeIdentifier.h"
#include "GaffParameters.h"

namespace Scine {

namespace Qmmm {
class InteractionTermEliminator;
} // namespace Qmmm

namespace MolecularMechanics {
class GaffParameters;
class AtomTypesHolder;
struct BondType;
struct AngleType;
struct DihedralType;
struct ImproperDihedralType;
class GaffMolecularMechanicsCalculator final
  : public Utils::CloneInterface<GaffMolecularMechanicsCalculator, MolecularMechanicsCalculator, Core::Calculator> {
 public:
  static constexpr const char* model = "GAFF";
  GaffMolecularMechanicsCalculator();
  ~GaffMolecularMechanicsCalculator() override = default;
  GaffMolecularMechanicsCalculator(const GaffMolecularMechanicsCalculator& rhs);
  void setStructure(const Utils::AtomCollection& structure) override;
  const Utils::Results& calculate(std::string description) override;
  std::string name() const override;

 private:
  // friend class declarations
  friend class Qmmm::InteractionTermEliminator;
  // Setter for the internal copy of the parameters.
  void setParameters(GaffParameters parameters);
  /*
   * @brief Implementation of a calculation.
   */
  const Utils::Results& calculateImpl(std::string description);
  /*
   * @brief Apply settings.
   */
  void applySettings();

  // Other private members:
  std::unique_ptr<BondsEvaluator> bondsEvaluator_;
  std::unique_ptr<AnglesEvaluator> anglesEvaluator_;
  std::unique_ptr<DihedralsEvaluator> dihedralsEvaluator_;
  std::unique_ptr<DihedralsEvaluator> improperDihedralsEvaluator_;
  std::unique_ptr<ElectrostaticEvaluator> electrostaticEvaluator_;
  std::unique_ptr<LennardJonesEvaluator> lennardJonesEvaluator_;
  bool printContributionsMolecularMechanics_;
  std::string atomicChargesFile_;
  std::string atomTypesFile_;
  bool onlyCalculateBondedContribution_;
  bool detectBondsWithCovalentRadii_;
  /*
   * Decides whether the potential terms that correspond to non-bonded interaction beyond the cutoff radius,
   * are ignored during the initialization stage of the calculator. This can result in a great speed-up, however,
   * the calculator may has to be re-initialized when the structure changes significantly.
   */
  bool applyCutoffDuringInitialization_;
  double nonCovalentCutoffRadius_;
  std::string connectivityFilePath_;
  std::string parameterFilePath_;
  /*
   * Tracks whether the parameter file path has been recently changed.
   * This boolean is set to true when the member 'parameterFilePath_' is updated
   * and it is set back to false when these parameters are parsed and stored in the
   * 'parameters_' member.
   */
  bool parameterFilePathHasBeenChanged_ = false;
  // Whether this calculator is currently calculating a Hessian -> no detailed output printing
  bool hessianMode_ = false;
  // Copy of the parameters as a member that is also cloned when the MM calculator is cloned.
  GaffParameters parameters_;
  bool parametersHaveBeenSetInternally_ = false;
  // Private methods:
  void initialize();
  void generatePotentialTerms(const GaffParameters& parameters, const IndexedStructuralTopology& topology,
                              const AtomTypesHolder& atomTypes);
  void generatePotentialTerms(const std::string& parameterPath);
};

} // namespace MolecularMechanics
} // namespace Scine

#endif // GAFFMOLECULARMECHANICSCALCULATOR_H