<div dir="ltr">Dear all,<div><br></div><div>I performed a single-point DFT-PBE simulation of an O2 molecule in a triplet state.</div><div>I noticed that the "DDAP SPIN DENSITY charges" on the O atoms have an opposite </div><div>sign with respect to Mulliken or Hirshfeld spin moments (see attached input/output files).</div><div>Is there any reason why the sign is reversed for the DDAPC spin moments, or am I</div><div>doing something wrong?</div><div><br></div><div>On the same topic, I'd be interested in using the charge or spin on selected atoms </div><div>as a reaction coordinate to model charge transfer processes, along the lines of this work: </div><div>M. Sulpizi at al. J. Theor. Comp. Chem 4, 985 (2006) </div><div>http://www.worldscientific.com/doi/abs/10.1142/S0219633605001957</div><div>I think one could use the existing restraining potential provided by the DDAPC_RESTRAINT</div><div>feature of CP2K and perform an umbrella sampling simulation. Is there any plan to </div><div>include the DDAPC charge/spin as a collective variable to perform metadynamics too?</div><div><br></div><div>Thanks a lot!</div><div>Simone</div><div><div>---</div><div>Simone Piccinin</div><div>CNR-IOM c/o SISSA</div><div>Via Bonomea 265</div><div>34136 Trieste</div><div>Tel: +39-040-3787317</div><div>email: picc...@iom.cnr.it</div></div><div> </div><div> </div><div><br></div></div>