Hi all,<div><br /></div><div>I'm trying to run QM/MM simulations on a solvated system with 40 atoms in the QM region using the M06L functional. I've run the same calculation using BLYP with no error, but when I specify M06L in my input file, the SCF calculations fail to converge after 21 iterations. The same convergence error occurs when I use revM06L.<br /><br />I've also tried implementing a finer grid with the additional parameters below, but get an error saying "tau with finer grids not implemented"<br />&XC_GRID<br />        USE_FINER_GRID<br />      &END XC_GRID<br /><br />Since I'm new to cp2k, I'm guessing there is an error in my input. Based on examples I've seen on here, I've tried both of these formats for specifying the functional (the latter says it doesn't recognize LIBXC as a subsection of XC_FUNCTIONAL).<br />&XC<br />      DENSITY_CUTOFF     1.0E-12<br />      GRADIENT_CUTOFF    1.0E-12<br />      TAU_CUTOFF         1.0E-12<br />      &XC_FUNCTIONAL<br />        &MGGA_X_M06_L<br />        &END MGGA_X_M06_L<br />        &MGGA_C_M06_L<br />        &END MGGA_C_M06_L<br />      &END XC_FUNCTIONAL<br />    &END XC<br /><br />--------------------------------------------------------------<br /><br />&XC<br />      DENSITY_CUTOFF     1.0E-12<br />      GRADIENT_CUTOFF    1.0E-12<br />      TAU_CUTOFF         1.0E-12<br />      &XC_FUNCTIONAL<br />        &LIBXC<br />         FUNCTIONAL MGGA_X_M06_L<br />        &END LIBXC<br />        &LIBXC<br />         FUNCTIONAL MGGA_C_M06_L<br />        &END LIBXC<br />      &END XC_FUNCTIONAL<br />    &END XC<br /><br />I also consistently get a warning that my restart file doesn't exist even though it does seem to be output during my calculation. I've attached examples of my input, output, md.log, and mdp files and have the commands I'm running below. Any tips would be greatly appreciated! <br /><br />I'm running gromacs 2022.5, cp2k 2024.1, and plumed 2.9<br /><br />Thanks,</div><div>Emma<br /><br />gmx_mpi_d grompp -f rxn_QMMM_DEF.mdp -c confout.gro -t state.cpt -n index.ndx -p topol.top -qmi rxn_QMMM_DEF_cp2k.inp -o rxn_QMMM_DEF.tpr<br />srun --mpi=pmi2 -n $SLURM_TASKS_PER_NODE gmx_mpi_d mdrun -ntomp $SLURM_CPUS_PER_TASK -s rxn_QMMM_DEF.tpr -plumed plumed.dat<br /><br /></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 on the web visit <a href="https://groups.google.com/d/msgid/cp2k/450eeb0b-3226-473d-9a87-55873c87ddd2n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/450eeb0b-3226-473d-9a87-55873c87ddd2n%40googlegroups.com</a>.<br />