mirror of
https://github.com/openmm/openmm
synced 2026-06-03 06:39:48 +09:00
16 lines
718 B
Python
16 lines
718 B
Python
from openmm.app import *
|
|
from openmm import *
|
|
from openmm.unit import *
|
|
from sys import stdout
|
|
|
|
pdb = PDBFile('input.pdb')
|
|
forcefield = ForceField('amber19-all.xml', 'amber19/tip3pfb.xml')
|
|
system = forcefield.createSystem(pdb.topology, nonbondedMethod=PME, nonbondedCutoff=1*nanometer, constraints=HBonds)
|
|
integrator = LangevinMiddleIntegrator(300*kelvin, 1/picosecond, 0.004*picoseconds)
|
|
simulation = Simulation(pdb.topology, system, integrator)
|
|
simulation.context.setPositions(pdb.positions)
|
|
simulation.minimizeEnergy()
|
|
simulation.reporters.append(DCDReporter('output.dcd', 1000))
|
|
simulation.reporters.append(StateDataReporter(stdout, 1000, step=True, potentialEnergy=True, temperature=True))
|
|
simulation.step(10000)
|