Dear Matt, <div><br><div>Thank you for your answer! </div><div>I tried to find the function to print cube files called by &DFT/&PRINT/&MO_CUBE but I lost in such a large number of codes of cp2k... </div><div>Actually at first I tried &R_LDOS by using:</div><div>&PDOS</div><div>  &R_LDOS</div><div>    ERANGE [eV] 0 2</div><div>    LIST 1..117</div><div>    XRANGE 0 $XMAX</div><div>    YRANGE 0 $YMAX</div><div>    ZRANGE 0 $ZMAX</div><div>  &END R_LDOS</div><div>  &EACH</div><div>    MD 1</div><div>  &END EACH</div><div>&END PDOS</div><div>but there's no cube file generated but only a file *. list, there's one line of data in it and that line has the similar format with those lines printed by &PDOS. That line just told me that there's one MO in the range I selected. </div><div>I went back to see if I misunderstood the function of R_LDOS. As manual said that "printing of local PDOS, projected on 3D volume in real space", however, I did not get any result of projection on 3d volume! </div><div>Did I make something wrong? </div><div><br></div><div>Thanks again,</div><div><br></div><div>Yike</div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">在2021年3月5日星期五 UTC+8 下午9:03:06<Matt W> 写道：<br></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">I don't think this is possible without hacking the code.<div><br></div><div>You could probably make</div><div><a href="https://manual.cp2k.org/cp2k-8_1-branch/CP2K_INPUT/FORCE_EVAL/DFT/PRINT/PDOS/R_LDOS.html" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=zh-CN&q=https://manual.cp2k.org/cp2k-8_1-branch/CP2K_INPUT/FORCE_EVAL/DFT/PRINT/PDOS/R_LDOS.html&source=gmail&ust=1615136671424000&usg=AFQjCNHJiL4XIemeOWMOmct2klshSjCfqQ">https://manual.cp2k.org/cp2k-8_1-branch/CP2K_INPUT/FORCE_EVAL/DFT/PRINT/PDOS/R_LDOS.html</a></div><div>select a particular state by choosng the correct energy range.</div><div><br></div><div>It is almost certainly easier to just generate the cubes and then delete the unwanted ones.</div><div><br></div><div>Matt</div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, March 5, 2021 at 11:46:03 AM UTC <a href="" data-email-masked="" rel="nofollow">yk...@gmail.com</a> 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">Dear cp2k users and developers, <div><br><div>I want to investigate some unoccupied molecular orbitals of particular interest. However, my cell has more than 120 atoms so there will be a quite large number of cube files generated if I simply use &FORCE_EVAL/&DFT/&PRINT/&MO_CUBES/NLUMO.</div></div><div>So I want to ask is it possible to print only MO cube files which are selected manually? </div><div><br></div><div>Thanks in advance:) </div><div><br></div><div>Yike Huang, PhD candidate</div><div>Dalian Institute of Chemical Physics, CAS, China. </div></blockquote></div></blockquote></div></div>