Comments (6)
Ok one need to add the line of parmed_structure.residues[0].name = "Lig"
from openff-bespokefit.
Hi @xiki-tempula glad you managed to make a bespoke force field! To export the parameters for a different simulation engine it is recommended that you make an OpenMM system using the OpenFF-Toolkit and then use OpenFF-Interchange or parmed as you have above. Here are some notebook examples using interchange or parmed, I hope this helps let m know if you have any issues.
from openff-bespokefit.
@jthorton Thanks. I think this is essentially what I'm doing in this code. I wonder if you mind help me on this?
I have attached the file that I'm using here. Thank you.
Archive.zip
from openff-bespokefit.
@xiki-tempula not quite you need to use the OpenFF-Toolkit to read offxml force field files rather than OpenMM directly, see the code attached which should work on your example. Let me know if there is anything else I can help with!
from openff.toolkit.typing.engines.smirnoff import ForceField
from openff.toolkit.topology import Molecule
import parmed
# load the molecule and the force field using the toolkit
molecule = Molecule.from_file("Archive/lig.sdf")
ff = ForceField("Archive/lig.offxml")
# create an openmm system
openmm_system = ff.create_openmm_system(molecule.to_topology())
# Convert OpenMM System to a ParmEd structure.
parmed_structure = parmed.openmm.load_topology(
molecule.to_topology().to_openmm(), openmm_system, molecule.conformers[0]
)
parmed_structure.save("system.top", overwrite=True)
parmed_structure.save("system.gro", overwrite=True)
from openff-bespokefit.
Thanks.
from openff-bespokefit.
Thanks. I wonder if it is possible to add the molecule name?
If you check the gromacs output, it is
[ moleculetype ]
; Name nrexcl
3
Which lacks the name of moleculetype.
This should normally be something like
[ moleculetype ]
; Name nrexcl
lig 3
from openff-bespokefit.
Related Issues (20)
- Error with linear angle during QC generation within a torsion drive HOT 2
- Optionally skip individual torsions when QC step fails
- About the way of adding more cpmds to an existing set
- Update installation docs HOT 2
- bespokefit gives the torsion a worse fit than the OpenFF HOT 19
- Automate running Psi4-based fits using default installation instructions HOT 5
- Use OpenFF 2.1.0 as default initial force field HOT 1
- Using psi4 HOT 11
- Test environments broken with OpenFF Toolkit 0.15.0 HOT 1
- SMILES shown in output is wrong HOT 11
- Provide a simple path to retrieving cached datasets, or package them with BespokeFit HOT 1
- BespokeFit appears to ignore `--qc-compute-n-cores` with XTB HOT 1
- Document environment variables HOT 4
- Python 3.12 support
- QCFractal 0.54 breaks tests HOT 1
- BespokeFit-client HOT 1
- Make docker container HOT 6
- Add API token support for HTTP endpoint? HOT 7
- LinearTorsionError for a nitrile group HOT 2
- Fragmentation failed
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from openff-bespokefit.