1
import "hexadecane.lt" # <- defines the "Hexadecane" molecule type.
3
# Periodic boundary conditions:
4
write_once("Data Boundary") {
10
molecules = new Hexadecane [12].move(0, 0, 5.2)
15
# NOTE: The spacing between molecules is large. There should be extra room to
16
# move during the initial stages of equilibration. However, you will have to
17
# run the simulation at NPT conditions later to compress the system to a
18
# more realistic density.