; position restraints for Na+ of THF Na+ I- [ position_restraints ] ; i funct fcx fcy fcz 1 1 1000 1000 1000