 -- Examples of i-PI working with LAMMPS --
 
 -- Example with Tersoff graphene potential --

   Uses the Tersoff parameters of Lindsay and Broido PRB 81, 205441 (2010)
   to run a short simulation of graphene with the LAMMPS MD code.
 
 -- Example with the q-TIP4P-F water potential --

 * This gives an example of water with the q-TIP4P-F potential
   of Scott Habershon, Thomas E. Markland and David E. Manolopoulos,
   J. Chem. Phys., 131, 024501, (2009).

 * State point: (N, V, T) = (216 water molecules, 43737 a_0**3, 298 K)

 * This demonstrates the convergence of the kinetic and potential energy
   using the ring polymer contraction (RPC) method of Thomas E. Markland and
   David E. Manolopoulos, J. Chem. Phys. 129, 024105, (2008), and
   the PIGLET method of Michele Ceriotti and 
   David Manolopoulos, Phys. Rev. Lett., 109, 100604, (2012).


** Run the examples automatically:

 * First, it is necessary to patch and compile LAMMPS, and create a make.in
   file containing the path to the lammps executable e.g.

```
LAMMPS:=~/bin/lmp_serial
```

 * Make sure you source the env.sh file in the root directory of i-pi to
   set all the necessary enviroment variables.
 
 * The runs can be done automatically using the Makefile provided. The make 
   targets are self-explanatory. To run the RPC example, for instance, just type:

$ make h2o-pimd.1

 * To clean up output files from all examples:

$ make clean


** Run the examples manually:

 * Go back to the example directory and run

$ i-pi input.xml
 
   the wrapper will start and sit waiting on the UDS /tmp/ipi.
 
 * Open a separate terminal and run the LAMMPS driver code using:
 
$ path/lammps/src/lmp_serial < in.lmp

   For the runs h2o-pimd+mts.4, h2o-pimd+rpc, h2o-pimd+sc+rpc instead use:

$ path/lammps/src/lmp_serial < in_long.lmp
$ path/lammps/src/lmp_serial < in_short.lmp

   You can run multiple instances of the code; it is so fast that parallel 
   scaling won't be appreciable.

 * If your system does not support Unix domain sockets, just set in input.xml
   <socket mode="inet"> <port> port_no </port>

   To make the client socket connect, the in.lmp file should 
   be changed so that it has the line:

   fix 1 all driver hostname port_no

   where hostname is the address the socket is binding to.
