<?xml version="1.0" ?><simulation verbosity="medium" threading="true">   
   <output prefix="simulation">    
      <properties stride="10"> [ step, time{picosecond}, conserved, temperature{kelvin}, pot_component_raw(0), pot_component_raw(1),pot_component_raw(2),pot_component_raw(3),pot_component_raw(4), potential, ensemble_temperature{kelvin} ] </properties>    
     <trajectory stride="10" filename="pos" cell_units="angstrom"> positions{angstrom} </trajectory>    
     <trajectory stride="10" filename="force">forces</trajectory>    
     <checkpoint stride="1000" overwrite="False"/>    
     </output>  
    <total_steps>1000000</total_steps>  
    <prng>    
        <seed>31865</seed>    
  </prng>  
       <ffsocket mode="inet" name="driver-lammps1">    
          <address>192.168.100.1</address>    
          <port>XXX0</port>    
          <slots>1</slots>    
          <timeout>600</timeout>    
       </ffsocket>  
       <ffsocket mode="inet" name="driver-lammps2">    
          <address>192.168.100.1</address>    
          <port>XXX1</port>    
          <slots>1</slots>    
          <timeout>600</timeout>    
       </ffsocket>  
       <ffsocket mode="inet" name="driver-lammps3">    
          <address>192.168.100.1</address>    
          <port>XXX2</port>    
          <slots>1</slots>    
          <timeout>600</timeout>    
       </ffsocket>  
     <ffsocket mode="inet" name="driver-dftb">    
            <address>192.168.100.1</address>    
            <port>XXX3</port>    
            <slots>1</slots>    
            <timeout>600</timeout>    
          </ffsocket>  
    <ffsocket mode="inet" name="driver-lammpsd1">    
            <address>192.168.100.1</address>    
            <port>XXX4</port>    
            <slots>1</slots>    
            <timeout>600</timeout>    
          </ffsocket>   
    <ffplumed name="plumed-1">    
             <init_file mode="xyz">./last_structure.xyz</init_file>    
             <plumeddat> plumed-1/plumed.in </plumeddat>    
    </ffplumed>  
    <system>    
        <initialize nbeads="1">      
              <file mode="xyz" units="angstrom"> last_structure.xyz </file>      
              <cell mode="abc" units="angstrom"> [18.428, 18.428, 18.428] </cell>      
            <velocities mode="thermal" units="kelvin"> 298.15 </velocities>      
        </initialize>    
        <forces>      
            <force forcefield="driver-lammps1" weight="0.334">        
              <mts_weights>[1.0,0]</mts_weights>               </force>      
            <force forcefield="driver-lammps2" weight="0.333">        
              <mts_weights>[1.0,0]</mts_weights>               </force>      
            <force forcefield="driver-lammps3" weight="0.333">        
              <mts_weights>[1.0,0]</mts_weights>               </force>      
            <force forcefield="driver-dftb">        
              <mts_weights>[1.0,0]</mts_weights>               </force>      
            <force forcefield="driver-lammpsd1">        
              <mts_weights>[-1.0,1.0]</mts_weights>      </force>      
          </forces>    
          <motion mode="dynamics">      
            <dynamics mode="nvt" splitting="baoab">        
              <thermostat mode="multi">          
                <thermostat mode="svr">                         <tau units="femtosecond"> 10.0 </tau>                       </thermostat>          
                <thermostat mode="gle">            
                  <A shape="(3,3)">
        [   3.023605375000e-3,    8.921193875950e-3,    5.150653685772e-3,
            8.921193875950e-3,    2.632422483298e-2,    1.519829829353e-2,
           -5.150653685772e-3,   -1.519829829353e-2,    0.000000000000e+0 ]
         </A>            
                 </thermostat>          
              </thermostat>        
              <nmts> [1,6] </nmts>        
               <timestep units="femtosecond"> 3.0 </timestep>        
             </dynamics>      
          </motion>    
        <ensemble>      
            <temperature units="kelvin"> 298.15 </temperature>      
          <bias>        
                    <force forcefield="plumed-1" nbeads="1"/>        
          </bias>      
        </ensemble>    
    </system>  
      <smotion mode="metad">    
             <metad>             <metaff> [ plumed-1 ] </metaff>           </metad>    
        </smotion>  
</simulation>