Transition Temperature for Spin-Crossover Materials with the Mean Value Ensemble Hubbard-U Correction

Calculation of transition temperatures T 1/2 for thermally driven spin-crossover in condensed phases is challenging, even with sophisticated state-of-the-art density functional approximations. The first issue is the accuracy of the adiabatic crossover energy difference ?E-HL between the low- and high-spin states of the bistable metal-organic complexes. The other is the proper inclusion of entropic contributions to the Gibbs free energy from the electronic and vibrational degrees of freedom. We discuss the effects of treatments of both contributions upon the calculation of thermochemical properties for a set of 20 spin-crossover materials using a Hubbard-U correction obtained from a reference ensemble spin-state. The U values obtained from a simplest bimolecular representation may overcorrect, somewhat, the ?E-HL values, hence giving somewhat excessive reduction of the T-1/2 results with respect to their U = 0 values in the crystalline phase. We discuss the origins of the discrepancies by analyzing different sources of uncertainties. By use of a first-coordination-sphere approximation and the assumption that vibrational contributions from the outermost atoms in a metal-organic complex are similar in both low- and high-spin states, we achieve T-1/2 results with the low-cost, widely used PBE generalized gradient density functional approximation comparable to those from the more costly, more sophisticated r2SCAN meta-generalized gradient approximation. The procedure is promising for use in high-throughput materials screening, because it combines rather low computational effort requirements with freedom from user manipulation of parameters.