7
$\begingroup$

I am studying Pd(II) metal-organic systems and have used GFN2-xTB (with implicit solvation) for high-throughput optimisation of candidate structures. I use the relative energies (from xtb) of different structural isomers as a screening metric based on empirical results. For the sake of validation, I aim to recreate those relative isomer energy calculations using DFT (optimisation, also with implicit solvation), which has previously been shown to agree with the xtb relative energies, for the top candidates.

I am aware that GFN2-xTB is parameterised with a focus on geometries, frequencies and NCIs - not energies, so validation (by DFT) of its use in my workflow is important.

I have a single case so far (although there may be more) where the relative isomer energies from DFT (PBE0/def2-svp/D3BJ) and xTB are opposing and in this case, I expect the difference to be due to a specific chemical functionality.

My question here is about the approaches one might use to understand this type of result. I do not have a lot of expertise in exploring the specific components of DFT/semiempirical calculations. So far, I have tried to understand this effect with a small model system of the different isomers that isolates the effect of the troublesome functionality. However, I am curious about the best way to analyse the results. I.e. is comparing the energy difference of all of the components (electronic energy, Gsolv, dispersion correction) appropriate?

If anyone could point me in the direction of a good resource on this type of issue, that would be appreciated. I am currently in the process of going back to some basic DFT reading.

$\endgroup$
  • 2
    $\begingroup$ +1. Welcome to our site, and thank you for asking your question here! We hope to see much more of you in the future !!! $\endgroup$ – Nike Dattani Sep 2 at 16:07
  • $\begingroup$ So there is one anomalous energy difference with DFT, and you want to know why? How "wrong" is this anomalous DFT result? How many kcal/mol? $\endgroup$ – Nike Dattani Sep 3 at 2:11
  • 1
    $\begingroup$ Differences of a few kcal/mol are already extremely small; the error in your DFT calculations can easily be of that magnitude. $\endgroup$ – Susi Lehtola Sep 3 at 17:36
  • 1
    $\begingroup$ Sorry, the unit associated with the number 10 in my comment should be kJ/mol. My understanding is that DFT error is somewhere around 5 kJ/mol (although that is something I have picked up somewhere and I have no memory as to where...). And my relative energies here are seeing a change in about 20 kJ/mol altogether (going from 10 more stable to 10 less stable is a relative change in 20). I would not expect such a change to be due to error, but if it was, how would you test or show that? $\endgroup$ – Andrew Tarzia Sep 4 at 8:07
  • 1
    $\begingroup$ Andrew Tarzia, next time, please tag me with @NikeDattani so that I get a notification when you've replied. I didn't see that you'd replied yesterday, and would have answered if I knew the answer to my question. I agree with Susi, a DFT energy changing by 5 kcal/mol is no surprise at all. The error in DFT is usually bigger than 5 kcal/mol and your basis set is too small to be so concerned about 5 kcal/mol differences in DFT. $\endgroup$ – Nike Dattani Sep 4 at 23:48
5
$\begingroup$

When modeling metal-organic systems, it is not uncommon for DFT to be wrong by far more than 5 kJ/mol. They are notoriously tricky to model. On the DFT side of things, I can say you're using a pretty small basis set of def2-SVP. That's going to be one factor limiting the quality of your DFT results. Also, Pd is pretty heavy, and relativistic corrections might be important. None of this is to mention the inherent limitations of PBE0.

On the semi-empirical side, GFN2-xTB is not meant for reliable energies of systems, and I'd expect it to do particularly poorly for transition metal complexes. It is going to be very difficult to understand why the two give very different results, seeing how highly empirical GFN2-xTB is by design. Without a thorough benchmarking comparing the energies of many different systems at the PBE0-D3(BJ)/def2-SVP level of theory and GFN2-xTB level of theory, it's going to be hard to make any claims. It should also be mentioned that there are a few differences here regardless. For instance, GFN2-xTB uses the D4 dispersion correction scheme, whereas you are using D3(BJ) at the DFT level. Furthermore, you are using an implicit solvation model with GFN2-xTB but not with DFT.

Since GFN2-xTB is supposed to be alright for vibrational frequencies, you can compare the vibrational contribution to the enthalpy or Gibbs free energy between PBE0-D3(BJ)/def2-SVP and GFN2-xTB to start. That will give you an idea if it's related to poor geometries and/or frequencies. Beyond that, you can compare the electronic energies. The only way to identify any systematic effect here is to compare the energies for many different systems to see if you can identify trends, and there simply may not be any. Overall, I personally would be very hesitant to compare the relative energies of any inorganic or metal-organic systems using any empirical method, including GFN2-xTB.

| cite | improve this answer | |
$\endgroup$
  • 1
    $\begingroup$ +1. I would have wrote the same thing :) I hadn't noticed that the user had answered my comment since I wasn't pinked, haha. $\endgroup$ – Nike Dattani Sep 5 at 0:32
  • $\begingroup$ Hey Andrew, thank you for your response! I agree with what you are saying and am aware of many of those limitations. I appreciate your idea of how to compare the contributions. Overall, I would tend to not do energetic comparisons (especially with xtb) between these systems because I would also be hesitant to trust them. This project has been interesting in the sense that the gfn2 energetics offer a rough guideline in screening, and then experiments and DFT have so far validated them but only for systems where GFNs approximations will lead to a systematic error (in my understanding). $\endgroup$ – Andrew Tarzia Sep 5 at 13:23
  • $\begingroup$ Thank you, @NikeDattani and Andrew for your help! This problem is pushing beyond my expertise! $\endgroup$ – Andrew Tarzia Sep 5 at 13:25

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.