1
import "spce.lt" # <- defines the "SPCE" (water) molecule type.
2
import "methane.lt" # <- defines the "Methane" molecule type (uses OPLSAA)
4
# Periodic boundary conditions:
5
write_once("Data Boundary") {
11
# The next command generates a (rather dense) cubic lattice with
12
# spacing 3.45 Angstroms. (The pressure must be equilibrated later.)
14
waters = new SPCE [12].move(0.00, 0.00, 3.45)
15
[12].move(0.00, 3.45, 0.00)
16
[12].move(3.45, 0.00, 0.00)
18
# Now add methane molecules:
20
methanes = new Methane [4].move(0, 0, 10.35)
24
# Move the methane molecules slightly to reduce overlap with the water
25
methanes[*][*][*].move(1.725, 1.725, 1.725)