<div dir="ltr"><div>Hi,</div><div>Thanks for the prompt reply. Then, the spin density for a given molecular orbital cannot be printed, can it be? I understand that the electronic density can be easily by postprocessing the wavefunctions written in the files generated with the command MO_CUBES, but what about the spin?</div><div><br></div><div> - Fernan<br></div><br>On Thursday, May 2, 2019 at 12:29:30 PM UTC+2, Matt W wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hi,<div><br></div><div>yes, if you print with E_DENSITY_CUBE for an unrestricted calculation you should get both total density and a spin density cube files.</div><div><br></div><div>Matt<br><br>On Thursday, May 2, 2019 at 11:26:39 AM UTC+1, Fernan Saiz wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Hello all,</div><div>I would like to ask how CP2k can print the spin densities of the HOMO and/or LUMO in a cube format of a spin-polarised system. I am not sure if the command E_DENSITY_CUBE can do it in combination with other ones.<br></div><div><br></div><div>Best regards,</div><div> - Fernan Saiz, PhD<br></div><div> Institut de Ciència de Materials de Barcelona (ICMAB-CSIC)</div><div> Bellaterra, Barcelona, 08193, Spain<br></div></div></blockquote></div></div></blockquote></div>