Not sure, sorry. It was a suggestion 'in principle'.<div><br></div><div>I think you are better off printing all the cubes and deleting unwanted ones. If it is really important to get particular orbitals printed out you'd have to hack the code.<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Saturday, March 6, 2021 at 5:30:14 PM UTC yk...@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;">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><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" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://manual.cp2k.org/cp2k-8_1-branch/CP2K_INPUT/FORCE_EVAL/DFT/PRINT/PDOS/R_LDOS.html&source=gmail&ust=1615488074892000&usg=AFQjCNHt5tpyrMxKPFnow0V2iktnbFIC1A">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 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></blockquote></div>