<div dir="ltr">Dear CP2K community,<div><br></div><div>I am computing a transition metal oxide for which I am applying DFT+U corrections in a supercell composed of 384 atoms.</div><div>To analyse structural properties at room temperature (295K) and standard pressure (1 atm) I need to carry out MD simulations. Due to the complexity of the system (of hexagonal symmetry), I start with the <b>NVT</b> ensemble and eventually reach equilibration. Electronic convergence runs smoothly and efficiently with the OT scheme. Since I do not know the range of the supercell volumes that corresponds to a pressure of 1atm, I started from the experimental values but, as expected, I get average pressures of the order of 500 to 1000 bars depending on the composition of the system. I could keep trying with different volumes and approach to the desire 1 atm pressure. This task, however, is very demanding and rather cumbersome.</div><div>I then decided to consider the <b>NPT_F</b> ensemble. However, the electronic DFT+U energy starts increasing significantly during the first steps of the simulation and, eventually, convergence problems start to manifest due to the Mulliken population analysis, with a negative contribution to the DFT+U energy. I could also try the LOWDIN method for charges but forces are not available for this scheme. Thus, it seems that my settings for the NPT_F ensemble are not appropriate for this DFT+U simulation. I describe the relevant MD settings of the input file below:</div><div><br></div><div><div><i>&Motion</i><br></div><div><i>  &MD</i></div><div><i>    Timestep 0.5</i></div><div><i>    Steps 2000</i></div><div><i>    Ensemble NPT_F</i></div><div><i>    Temperature 300.0</i></div><div><i>    &Thermostat</i></div><div><i>       REGION massive</i></div><div><i>       TYPE CSVR</i></div><div><i>       &CSVR</i></div><div><i>         TIMECON [fs] 10.</i></div><div><i>       &END CSVR</i></div><div><i>    &End Thermostat</i></div><div><i><br></i></div><div><i>   &Barostat</i></div><div><i>     Pressure 1.0</i></div><div><i>     TIMECON [fs] 10.<br></i></div><div><i>   &End barostat<br></i></div></div><div><i><br></i></div><div><i>&End Motion</i></div><div><i><br></i></div><div> The message I obtain is the following</div><div><i> *** WARNING in dft_plus_u.F:1457 :: DFT+U energy contibution is negative ***</i><br></div><div><div><i> *** possibly due to unphysical Mulliken charges. Check your input, if    ***</i></div><div><i> *** this warning persists or try a different method!                     ***</i></div></div><div><i><br></i></div><div>Again, this problems does not arise with NVT simlation, only with NPT. I am wondering if I need of an extra flag in the motion/thermostat session to "tame" the MD run and the electronic convergence. I have carried out an extensive research in the web but it seems there is no previous record of such a type of simulation with CP2K. </div><div><br></div><div>Any help would be appreciated.</div><div>Thank you very much in advance.</div><div><br></div><div>Ivan</div></div>