Spurious Oscillations Caused by Density Functional Approximations: Who is to Blame? Exchange or Correlation?

We analyze the varying susceptibilities of different density functional approximations (DFAs) to present spurious oscillations on the profiles of several vibrational properties. Among other problems, these spurious oscillations cause significant errors in harmonic and anharmonic IR and Raman frequencies and intensities. This work hinges on a judicious strategy to dissect the exchange and correlation components of DFAs and pinpoint the origins of these oscillations. We identify spurious oscillations in derivatives of all energy components with respect to nuclear displacements, including those energy terms that do not involve numerical integrations. These indirect spurious oscillations are attributed to suboptimal electron densities resulting from a self-consistent field procedure using a DFA that exhibits direct spurious oscillations. Direct oscillations stem from inaccurate numerical integration of the exchange and correlation energy density functionals. A thorough analysis of direct spurious oscillations reveals that only a handful of exchange and correlation components are insensitive to spurious oscillations, giving rise to three families of functionals, BH&H, LSDA, and BLYP. Among the functionals in these families, we encounter four widespread DFAs: BLYP, B3LYP, LC-BLYP, and CAM-B3LYP. Certain DFAs like PBE appear less sensitive to spurious oscillations due to compensatory cancellations between their energy components. Additionally, we found non-negligible but small oscillations in PBE and TPSS, which could be safely employed provided a sufficiently large integration grid is used in the calculations. These findings hint at the key components of current approximations to be improved and emphasize the necessity to develop accurate DFAs suitable for studying molecular spectroscopies.

In-situ formation of one-dimensional coordination polymers in molecular junctions.

We demonstrate the bottom-up in-situ formation of organometallic oligomer chains at the single-molecule level. The chains are formed using the mechanically controllable break junction technique operated in a liquid environment, and consist of alternating isocyano-terminated benzene monomers coordinated to gold atoms. We show that the chaining process is critically determined by the surface density of molecules. In particular, we demonstrate that by reducing the local supply of molecules within the junction, either by lowering the molecular concentration or by adding side groups, the oligomerization process can be suppressed. Our experimental results are supported by ab-initio simulations, confirming that the isocyano terminating groups display a high tendency to form molecular chains, as a result of their high affinity for gold. Our findings open the road for the controlled formation of one-dimensional, single coordination-polymer chains as promising model systems of organometallic frameworks.