Post
``` from openff.toolkit import Molecule, ForceField from tqdm import tqdm # for progress bars mols = Molecule.from_file("inputs.sdf") # loads 50 different molecules ff = ForceField("openff_unconstrained-2.3.0-rc2.offxml") ...
0 reposts0 likes
... for idx, mol in enumerate(tqdm(mols, desc="Processing molecules")): interchange = ff.create_interchange(mol.to_topology()) interchange.minimize() mol.conformers[0] = interchange.positions mol.to_file(f"minimized_mol_{idx}.sdf", file_format="sdf") ```