<br>Thanks for your help. I'll try . <br><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, December 16, 2020 at 11:57:37 PM UTC+8 Maxime Van den Bossche 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;"><div>I forgot to add that, if you are content with Mulliken charges,</div><div>then these should be more readily available (presumably with</div><div>CP2K_INPUT / FORCE_EVAL / DFT / PRINT / MULLIKEN ).</div><div>Maxime<br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">Op woensdag 16 december 2020 om 14:44:43 UTC+1 schreef Maxime Van den Bossche:<br></div></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>Dear Evelyn,</div><div><br></div><div>I'm not a CP2K developer, but I suspect that the following is at play here:</div><div><br>In a DFTB implementation, the basis functions themselves do not need to be known,</div><div>as the required information for e.g. building the Hamiltonian and overlap matrices</div><div>is provided via the (one- and two-center) integrals in the Slater-Koster files.</div><div><br></div><div>So, from only those files, the DFTB module in e.g. CP2K would not know how to</div><div>construct the electron density (or molecular orbitals, for that matter) from e.g.</div><div>the wavefunction coefficients or the density matrix.</div><div><br></div><div>In the DFTB+ code, this can be done (in a post-processing step) by describing</div><div>the basis functions using Slater-type orbitals, see e.g.</div><div><br></div><div><a href="https://dftbplus-recipes.readthedocs.io/en/latest/basics/waveplot.html" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://dftbplus-recipes.readthedocs.io/en/latest/basics/waveplot.html&source=gmail&ust=1608474576785000&usg=AFQjCNFr0eJnzsqZGjO0uYW9KloBBu87DA">https://dftbplus-recipes.readthedocs.io/en/latest/basics/waveplot.html</a></div><div><a href="https://dftbplus.org/fileadmin/DFTBPLUS/public/dftbplus/latest/manual.pdf#chapter.7" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://dftbplus.org/fileadmin/DFTBPLUS/public/dftbplus/latest/manual.pdf%23chapter.7&source=gmail&ust=1608474576786000&usg=AFQjCNEKqRK2Hc1Lgxx_zL8eTP0Guy6X1w">https://dftbplus.org/fileadmin/DFTBPLUS/public/dftbplus/latest/manual.pdf#chapter.7</a></div><div><br></div><div>At least, I think this explains why this is not (yet) implemented in CP2K.</div><div><br></div><div>If you can get your hands on e.g. the wave function coefficients, and have the required</div><div>information on the basis functions, you may of course compute the electron density</div><div>yourself.<br></div><div><br></div><div>Best,</div><div>Maxime</div>
<br><div class="gmail_quote"><div dir="auto" class="gmail_attr">Op zaterdag 12 december 2020 om 14:29:13 UTC+1 schreef <a rel="nofollow">ea...@gmail.com</a>:<br></div></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">hello<br>I'm a beginner, and I want to try running a monolayer of HBN using the DFTB method, and then output the charge of the atom.But there are always mistakes 。I want to know what the problem is.The display error is Electronic density cube file not implemented for DFTB method。<div><div>#CPQA INCLUDE DFTB/nonscc/nonscc_parameter</div><div>#CPQA INCLUDE uff_table</div><div>#CPQA INCLUDE DFTB/nonscc/bn</div><div>#CPQA INCLUDE DFTB/nonscc/nb</div><div>&FORCE_EVAL</div><div> METHOD Quickstep</div><div> &DFT</div><div> &QS</div><div> METHOD DFTB</div><div> &DFTB</div><div> SELF_CONSISTENT F</div><div> &PARAMETER</div><div> PARAM_FILE_PATH DFTB/nonscc</div><div> PARAM_FILE_NAME nonscc_parameter</div><div> &END PARAMETER</div><div> &END DFTB</div><div> &END QS</div><div> &SCF</div><div> SCF_GUESS NONE</div><div> &MIXING</div><div> METHOD DIRECT_P_MIXING</div><div> ALPHA 1.</div><div> &END</div><div> MAX_SCF 5</div><div> &END SCF</div><div> &PRINT</div><div> </div><div> &E_DENSITY_CUBE</div><div> FILENAME cube</div><div> STRIDE 1 1 1</div><div> &END E_DENSITY_CUBE</div><div> </div><div> &END PRINT</div><div> &END DFT</div><div> &SUBSYS</div><div> &CELL</div><div> ABC 64 63 2</div><div> PERIODIC NONE</div><div> &END CELL</div><div> &TOPOLOGY</div><div> COORD_FILE_FORMAT PDB</div><div> COORD_FILE_NAME HBN1.pdb</div><div> &END TOPOLOGY</div><div> &END SUBSYS</div><div>&END FORCE_EVAL</div><div>&GLOBAL</div><div> PROJECT hbn</div><div> RUN_TYPE MD</div><div> PRINT_LEVEL LOW</div><div>&END GLOBAL</div><div>&MOTION</div><div> &MD</div><div> ENSEMBLE NVT </div><div> STEPS 10</div><div> TIMESTEP 1</div><div> TEMPERATURE 300.0</div><div> &THERMOSTAT</div><div> &NOSE</div><div> LENGTH 3</div><div> YOSHIDA 3</div><div> TIMECON [wavenumber_t] 2000</div><div> MTS 2</div><div> &END NOSE</div><div> &END</div><div> &END MD</div><div>&END MOTION</div><div><br></div>Thank you in advance for your help! <br><br>Best regards，<br>Evelyn <br></div></blockquote></div></blockquote></div></blockquote></div>