<simulation mode="md" verbosity="medium">
   <output prefix='simulation'>
      <properties stride='10' filename='out'> [ step, time{picosecond}, conserved{electronvolt}, temperature{kelvin}, pressure_cv{megapascal} ] </properties>
      <trajectory filename='pos' stride='10'> positions{angstrom} </trajectory>
      <trajectory filename='xc' stride='10'> x_centroid{angstrom} </trajectory>
      <trajectory filename='vc' stride='10'> v_centroid </trajectory>
      <checkpoint stride='4000'/>
   </output>
   <total_steps> 40 </total_steps>
   <!-- <total_time> 840000000 </total_time> -->
   <prng><seed> 31415 </seed></prng>
   <ffsocket name="ase" mode="unix">
      <address> driver-irpmd-16 </address>
   </ffsocket>
   <system>
      <initialize nbeads='8'>
         <file mode='xyz' units='angstrom'> firstframe.xyz </file>
         <velocities mode="thermal" units='kelvin'> 300 </velocities>
         <cell units='angstrom' mode='abc'> [ 12.412, 12.412, 12.412 ] </cell>
      </initialize>
      <forces><force forcefield="ase"> </force></forces>
      <motion mode='dynamics'>
         <dynamics mode='nvt'>
           <timestep units='femtosecond'> 0.25 </timestep>
           <thermostat mode='pile_g'>
             <tau units='femtosecond'> 100 </tau>
             <pile_lambda> 0.5 </pile_lambda>
           </thermostat> 
         </dynamics>
      </motion>
      <ensemble>
         <temperature units='kelvin'> 300 </temperature>
      </ensemble>
   </system>
</simulation>
