<div dir="ltr">Hi Vladimir,<div> <div>Thank you so much! I've noticed the Mulliken's shortage, so I'm wondering whetther I can get the right information by qualitative estimate based on the Mulliken. </div><div><br></div><div>Here are two figures. In first one, the green line is the react, the purple line is the product. And the 142O 143O are the focused atoms(figure 2). It shows, after dissociation, the 142O 143O have a spin moment peak, though the value is small, but the contrast is obvious. Can it prove the radical thing?</div><div><br><div><br>在 2020年4月6日星期一 UTC+8下午9:39:04,Vladimir Rybkin写道:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear Du,<div><br></div><div>in principle, population analysis will help, although Mulliken should be used with care. There are other options in CP2K as well, e.g. Hirshfeld, RESP. It may also be useful to print the spin density cube and analyse spin density distribution visually. </div><div><br></div><div>It's not clear to me what the dissociation products are and which atoms in the output belong to them, but as far as I see, spin moments on all atoms are very small, i.e. there's no considerable spin polarisation. So, I would say you don't have radicals.</div><div><br></div><div>Yours, </div><div><br></div><div>Vladimir<br><br>понедельник, 6 апреля 2020 г., 3:16:59 UTC+2 пользователь Du написал:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear developers and users,<div><blockquote style="margin:0 0 0 40px;border:none;padding:0px"><div>Hello! I have a problem in my work, and I hope someone can give me some guidance. My system is 40H2O + HSO5- + H+ + NGP, about 180 atoms. and I focus on the HSO5- dissociation. Now the point is that I don't know whether the dissociation product species are radicals or not. So I cut a piece of the MD simulation when HSO5- gets dissociated, to do a single point energy calculation in cp2k, with the mulliken polarization analysis. So my question is, Can I figure out the species are radicals or not based on the information in Mulliken polarization analysis? If it can, I attach the Mulliken output profile, and I wonder how to use it? If not, What other ways can I use to do it? </div></blockquote></div><blockquote style="margin:0 0 0 40px;border:none;padding:0px"><div>Thank you!</div></blockquote><div style="text-align:right">Du </div><blockquote style="margin:0 0 0 40px;border:none;padding:0px"><div> </div></blockquote><blockquote style="margin:0 0 0 40px;border:none;padding:0px"><div><div> </div></div></blockquote></div></blockquote></div></div></blockquote></div></div></div></div>