<div>Hi 
Nimish,</div><div><br></div><div>I think the e density cube can only be printed when using quickstep, i.e. non-pw dft. I'm not sure how difficult it would be to extend these printing functions for the pw density too.</div><div></div><div><br></div><div>Cheers,</div><div>Fabian<br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Thursday, 4 February 2021 at 18:48:35 UTC+1 nimi...@gmail.com wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Hi,<br><br>I want to obtain the electron density in the Gaussian cube format of a pure Plane Wave DFT calculation using the 'SIRIUS' library within CP2K-7.1. I have modified an example from the CP2K website (<a href="https://www.cp2k.org/howto:running_qe_computation" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://www.cp2k.org/howto:running_qe_computation&source=gmail&ust=1612552475659000&usg=AFQjCNE6AYTTjFGezryJ4bpqob1O4qLBTw">https://www.cp2k.org/howto:running_qe_computation</a>) by including the section controlling the printing of electron density cube file. The calculation terminates normally but does not produce any Gaussian cube file. <br><br>I have attached the input script (Si7Ge_cube_input.txt) that I used and the outputs from the CP2K program (Si7Ge_cube_output.txt & sirius_output.txt) with this email. <br><br>Any input that helps me resolve the issue will be much appreciated.<br><img alt=""><br>Regards,<br>Nimish Dwarkanath,<br>Ph.D. student,<br>JNCASR, Bengaluru (INDIA)</blockquote></div>