In the attachment I added a debug report for this bug.  I have circumvented the problem, and works on my own computer,
    although my hack still needs testing on our computing clusters. 
    Providing a in-depth full solution is beyond me since I'm not a CP2K
    developer, but I hope it is of use for someone else.<br><br>best regards<br>marc<br><br>On Friday, February 8, 2013 10:22:41 AM UTC+1, marc wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;">Can this problem be related to this thread somehow?<br><a href="https://groups.google.com/forum/?hl=en&fromgroups=#!searchin/cp2k/Marc$20Van$20Houteghem/cp2k/t5UYWtZt3HI/DlfmK2W-a28J" target="_blank">https://groups.google.com/<wbr>forum/?hl=en&fromgroups=#!<wbr>searchin/cp2k/Marc$20Van$<wbr>20Houteghem/cp2k/t5UYWtZt3HI/<wbr>DlfmK2W-a28J</a><br><br>No succes yet on this matter. I added a big molecular methanol system (128 molecules) and a small one (2 molecules) that reproduce the problem.  Important to mention is that the computation only fails when run in parallel mode.  It does not occur with a NVT ensemble and/or in serial mode (only the combination NPT ensemble and parallel computation fails).  <br>Probably this is due to the stress tensor implementation.<br><br>marc<br><br>On Tuesday, November 27, 2012 3:35:51 PM UTC+1, Noam Bernstein wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex">Hi - I think I've found a small bug in the alchemy code, and I'm not sure how<div>best to fix it. For example for</div><div>     FORCE_EVAL%MIXED%MIXING_TYPE LINEAR_COMBINATION,</div><div>there are two calls to mixed_map_forces(), the first with overwrite=.true., the</div>
<div>second with overwrite=.false.  However, if the first FORCE_ENV has some of</div><div>the atoms missing (e.g. alchemy), the atoms that are missing never have</div><div>their forces overwritten.  Presumably those forces should be set to 0 before</div>
<div>the forces from the second FORCE_ENV are added.  </div><div><br></div><div>Was there some particular reason that the mixed forces aren't all just initialized to 0, </div><div>and instead the overwrite argument is used?  Because I don't see a reliable way of</div>
<div>fixing this problem without either some tedious bookkeeping to find</div><div>all atoms that map to nothing in the overwrite=.true. case, or just setting</div><div>all the forces to 0 before the calls to mixed_map_forces().  I'd be happy</div>
<div>to submit a patch, but if there's some reason it's not a good idea to set</div><div>all the forces to 0, I'd like to understand that.</div><div><br></div><div>           Noam</div>
</blockquote></blockquote>