Single-Atom Substituents in Copper Surfaces May Adsorb Multiple CO Molecules.

Copper is a good CO2 electroreduction catalyst as products beyond CO form, but efficiency and selectivity is low. Experiments have shown that an admixture of other elements can help, and computational screening studies have pointed out various promising candidates based on the adsorption of a single CO molecule as a descriptor. Our calculations of CO adsorption on surfaces where a first row transition metal atom replaces a Cu atom show that multiple CO molecules, not just one, bind to the substitutional atom. For Fe, Co, and Ni atoms, a decrease in binding energy is found, but the reverse trend, namely, increasing bond strength, is found for V, Cr, and Mn and the first three CO molecules. Magnetic moment, charge, and position of the substitutional atom are also strongly affected by the CO adsorption in most cases. Magnetic moment is stepwise reduced to zero, and the outward displacement of the substitutional atom increased.

Evidence of Sharp Transitions between Octahedral and Capped Trigonal Prism States of the Solvation Shell of the Aqueous Fe3+ Ion.

The structure of the solvation shell of the aqueous Fe3+ ion has been a subject of controversy due to discrepancies between experiments and different levels of theory. We address this issue by performing simulations for a wide range of ion concentrations, using various potential energy functions, supplemented by density functional theory calculations of selected configurations. The solvation shell undergoes abrupt transitions between two states: a hexacoordinated octahedral (OH) state and a capped trigonal prism (CTP) state with 7-fold coordination. The lifetime of these states is dependent on concentration. In dilute FeCl3 solutions, the lifetimes of both are similar (≈1 ns). However, the lifetime of the OH state increases with ion concentration, while that of the CTP state decreases slightly. When a uniform negative background charge is used instead of explicit counterions, the lifetime of the OH state is greatly overestimated. These findings underscore the need for further experimental measurements and higher-level simulations.

Variational Density Functional Calculations of Excited States: Conical Intersection and Avoided Crossing in Ethylene Bond Twisting.

Theoretical studies of photochemical processes require a description of the energy surfaces of excited electronic states, especially near degeneracies, where transitions between states are most likely. Systems relevant to photochemical applications are typically too large for high-level multireference methods, and while time-dependent density functional theory (TDDFT) is efficient, it can fail to provide the required accuracy. A variational, time-independent density functional approach is applied to the twisting of the double bond and pyramidal distortion in ethylene, the quintessential model for photochemical studies. By allowing for symmetry breaking, the calculated energy surfaces exhibit the correct topology around the twisted-pyramidalized conical intersection even when using a semilocal functional approximation, and by including explicit self-interaction correction, the torsional energy curves are in close agreement with published multireference results. The findings of the present work point to the possibility of using a single determinant time-independent density functional approach to simulate nonadiabatic dynamics, even for large systems where multireference methods are impractical and TDDFT is often not accurate enough.

Method for Calculating Excited Electronic States Using Density Functionals and Direct Orbital Optimization with Real Space Grid or Plane-Wave Basis Set.

A direct orbital optimization method is presented for density functional calculations of excited electronic states using either a real space grid or a plane-wave basis set. The method is variational, provides atomic forces in the excited states, and can be applied to Kohn-Sham (KS) functionals as well as orbital-density-dependent (ODD) functionals including explicit self-interaction correction. The implementation for KS functionals involves two nested loops: (1) An inner loop for finding a stationary point in a subspace spanned by the occupied and a few virtual orbitals corresponding to the excited state; (2) an outer loop for minimizing the energy in a tangential direction in the space of the orbitals. For ODD functionals, a third loop is used to find the unitary transformation that minimizes the energy functional among occupied orbitals only. Combined with the maximum overlap method, the algorithm converges in challenging cases where conventional self-consistent field algorithms tend to fail. The benchmark tests presented include two charge-transfer excitations in nitrobenzene and an excitation of CO to degenerate π* orbitals where the importance of complex orbitals is illustrated. The application of this method to several metal-to-ligand charge-transfer and metal-centered excited states of an FeII photosensitizer complex is described, and the results are compared to reported experimental estimates. This method is also used to study the effect of the Perdew-Zunger self-interaction correction on valence and Rydberg excited states of several molecules, both singlet and triplet states, and the performance compared to semilocal and hybrid functionals.

Mn Dimer Can Be Described Accurately with Density Functional Calculations When Self-Interaction Correction Is Applied.

Qualitatively incorrect results are obtained for the Mn dimer in density functional theory calculations using the generalized gradient approximation (GGA), and similar results are obtained from local density and meta-GGA functionals. The coupling is predicted to be ferromagnetic rather than antiferromagnetic, and the bond between the atoms is predicted to be an order of magnitude too strong and approximately an Ångstrøm too short. Explicit, self-interaction correction (SIC) applied to a commonly used GGA energy functional, however, provides close agreement with both experimental data and high-level, multireference wave function calculations. These results show that the failure is not due to a strong correlation but rather the single electron self-interaction that is necessarily introduced in estimates of the classical Coulomb and exchange-correlation energy when only the total electron density is used as the input. The corrected functional depends explicitly on the orbital densities and can, therefore, avoid the introduction of a self-Coulomb interaction. The error arises because of an overstabilization of bonding d-states in the minority spin channel resulting from an overestimate of the d-electron self-interaction in the semilocal exchange-correlation functionals. Since the computational effort in the SIC calculations scales with the system size in the same way as for regular semilocal functional calculations, this approach provides a way to calculate properties of Mn nanoclusters as well as biomolecules and extended solids, where Mn dimers and larger cluster are present, while multireference wave function calculations can only be applied to small systems.

Theory and Applications of Generalized Pipek-Mezey Wannier Functions.

The theory for the generation of Wannier functions within the generalized Pipek-Mezey approach (Lehtola, S.; Jónsson, H. J. Chem. Theory Comput. 2014, 10, 642) is presented and an implementation thereof is described. Results are shown for systems with periodicity in one, two, and three dimensions as well as isolated molecules. The generalized Pipek-Mezey Wannier functions (PMWF) are highly localized orbitals consistent with chemical intuition where a distinction is maintained between σ- and π-orbitals. The PMWF method is compared with the so-called maximally localized Wannier functions (MLWFs) that are frequently used for the analysis of condensed matter calculations. Whereas PMWFs maximize the localization criterion of Pipek and Mezey, MLWFs maximize that of Foster and Boys and have the disadvantage of mixing σ- and π-orbitals in many cases. The PMWF orbitals turn out to be as localized as the MLWF orbitals as evidenced by cross-comparison of the values of the PMWF and MLWF objective functions for the two types of orbitals. Our implementation in the atomic simulation environment (ASE) is compatible with various representations of the wave function, including real-space grids, plane waves, and linear combinations of atomic orbitals. The projector-augmented wave formalism for the representation of atomic core electrons is also supported. Results of calculations with the GPAW software are described here, but our implementation can also use output from other electronic structure software such as ABINIT, NWChem, and VASP.

Implementation of Constrained DFT for Computing Charge Transfer Rates within the Projector Augmented Wave Method.

Combining constrained density function theory (cDFT) with Marcus theory is an efficient and promising way to address charge transfer reactions. Here, we present a general and robust implementation of cDFT within the projector augmented wave (PAW) framework. PAW pseudopotentials offer a reliable frozen-core electron description across the whole periodic table, with good transferability, as well as facilitate the extraction of all-electron quantities. The present implementation is applicable to two different wave function representations, atomic-centered basis sets (LCAO) and the finite-difference (FD) approximation utilizing real-space grids. LCAO can be used for large systems, molecular dynamics, or quick initialization, while more accurate calculations are achieved with the FD basis. Furthermore, the calculations can be performed with flexible boundary conditions, ranging from isolated molecules to periodic systems in one-, two-, or three-dimensions. As such, this implementation is relevant for a wide variety of applications. We also present how to extract the electronic coupling element and reorganization energy from the resulting diabatic cDFT-PAW wave functions for the parametrization of Marcus theory. Here, the combined method is applied to important test cases where practical implementations of DFT fail due to the self-interaction error, such as the dissociation of the helium dimer cation, and it is compared to other established cDFT codes. Moreover, for charge localization in a diamine cation, where it was recently shown that the commonly used generalized gradient and hybrid functionals of DFT failed to produce the localized state, cDFT produces qualitatively and quantitatively accurate results when benchmarked against self-interaction corrected DFT and high-level CCSD(T) calculations at a fraction of the computational cost.

Effect of Complex-Valued Optimal Orbitals on Atomization Energies with the Perdew-Zunger Self-Interaction Correction to Density Functional Theory.

The spurious interaction of an electron with itself-self-interaction error-is one of the biggest problems in modern density-functional theory. Some of its most glaring effects, such as qualitatively incorrect charge separation upon dissociation, can be removed by an approximate self-interaction correction due to Perdew and Zunger (PZ) (Perdew, J.; Zunger, A. Phys. Rev. B 1981, 23, 5048). However, the correction introduces an explicit dependence on the occupied orbital densities, which makes proper minimization of the functional difficult. Previous work (Vydrov et al., J. Chem. Phys. 2006, 124, 094108) has suggested that the application of the PZ correction results in worse atomization energies than those obtained with the uncorrected parent functional. But, it has only recently been found that the correct minimization of the PZ energy functional requires complex-valued orbitals, which have not been used in previous work on atomization energies. Here, we study the effect of the proper use of complex-valued orbitals on the atomization energies of molecules in the W4-11 data set (Karton, A.; Daon, S.; Martin, J. M. Chem. Phys. Lett. 2001, 510, 165). We find that the correction has a tendency to weaken the binding of molecules. The correction using complex-valued orbitals is invariably found to yield better atomization energies than the correction with real-valued orbitals. The correction applied to the PBEsol exchange-correlation functional results in a functional that is more accurate than the (uncorrected) PBE functional.

Direct measurement and modulation of single-molecule coordinative bonding forces in a transition metal complex.

Coordination chemistry has been a consistently active branch of chemistry since Werner's seminal theory of coordination compounds inaugurated in 1893, with the central focus on transition metal complexes. However, control and measurement of metal-ligand interactions at the single-molecule level remain a daunting challenge. Here we demonstrate an interdisciplinary and systematic approach that enables measurement and modulation of the coordinative bonding forces in a transition metal complex. Terpyridine is derived with a thiol linker, facilitating covalent attachment of this ligand on both gold substrate surfaces and gold-coated atomic force microscopy tips. The coordination and bond breaking between terpyridine and osmium are followed in situ by electrochemically controlled atomic force microscopy at the single-molecule level. The redox state of the central metal atom is found to have a significant impact on the metal-ligand interactions. The present approach represents a major advancement in unravelling the nature of metal-ligand interactions and could have broad implications in coordination chemistry.