cmelab / cmelab/pl-validation

Generation of forcefields

Open
#9 1 comment 0 reactions 2 assignees Claimed by @jbieri View on GitHub
Dominant language
Jupyter Notebook
Stars
0
Forks
0
PR merge metrics
No merged PRs in 30d

Description

When generating a new forcefield from a smiles string using ff_generation.ipynb the resulting mol2 and forcefield do not polymerize. Polymerization works exactly as intended, but when running the code below the following error occurs. Additionally the monomer does not error out when running a hoomd simulation. The full error can be found in the attached .txt file

```
espmol = espaloma_mol(system_file)
molecule = Molecule(num_mols=1, compound=espmol)

molff = FF_from_file(ff_filepath)
system = Pack(molecules=molecule,density=0.01 * u.g/u.cm**3,packing_expand_factor=5)
system.apply_forcefield(r_cut=2.5, force_field=molff, auto_scale=True,remove_charges=True, remove_hydrogens=True)
```
```
ParameterizationError: No parameters found for connection ,
group: dihedral_type, identifiers: [['N2', 'C0', 'C0', 'S1'], ['N2', 'C0', 'C0', 'S1']] in the Forcefield.
```
[ParameterizationError.txt](https://github.com/user-attachments/files/20778806/ParameterizationError.txt)

Contributor guide

No contributing guide indexed for this repository

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.