Hallo <div><br /></div><div>Please have a look at the examples in https://manual.cp2k.org/trunk/methods/dft/constrained.html</div><div>One of the problem of your input seems to be that no optimiser has been indicated in the outer_scf section. </div><div>Regards</div><div>Marcella<br /><br /></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, July 29, 2025 at 9:31:01 AM UTC+2 vishnu...@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 all,<div><br></div><div>I am trying to do a Constrained DFT geometry relaxation for a semiconductor-dopant pair. The aim is to put constrain a positive charge on the semiconductor (+1 for e.g) and an equivalent negative charge on the dopant (-1) while keeping the overall system neutral. </div><div>However, when the code reaches the CDFT part it crashes without specifying the error.</div><br> CDFT EXTERNAL SCF WAVEFUNCTION OPTIMIZATION <br><br>  ---------------------------------- CDFT -------------------------------------- <br><br>  Optimizing a density constraint in an external SCF loop  <br><br>  Type of constraint:     Hirshfeld <br><br>  Number of constraints:          2 <br><br>  Using fragment densities:       F <br><br>  Calculating atomic CDFT charges <br><br><br> ******************************************************************************* <br><br> *   ___                                                                       * <br><br> *  /   \                                                                      * <br><br> * [ABORT]                                                                     * <br><br> *  \___/                                                                      * <br><br> *    |                                                                        * <br><br> *  O/|                                                                        * <br><br> * /| |                                                                        * <br><br> * / \                                                     qs_scf_output.F:976 * <br><br> *******************************************************************************<div><br></div><div>Typically CP2k gives a very clear error message, but here there is nothing. I have attached my input file and the full output till the error. The system has 78 atoms, with the first 58 atoms being the semiconductor, and the remaining 20 being dopant. My CP2K version is 8.2, built using GCC 11.2.0 on the Purdue Anvil supercomputer.</div><div><br></div><div>I am new to constrained DFT and CP2k, so apologies if there are some very basic mistakes that I've missed. </div><div><br></div><div>Thank you, I am immensely grateful for any help. </div><div><br></div><div>Regards,</div><div>Vishnu Raghuraman</div><div>Postdoctoral Research Associate, Chemistry Department</div><div>University of Illinois at Urbana-Champaign</div></blockquote></div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/5dbfdca0-829f-432d-b519-5b3c517184f8n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/5dbfdca0-829f-432d-b519-5b3c517184f8n%40googlegroups.com</a>.<br />