Ab initio benchmark study for the oxidative addition of CH4 to Pd: importance of basis-set flexibility and polarization

To obtain a state-of-the-art benchmark potential energy surface (PES) for the archetypal oxidative addition of the methane C-H bond to the palladium atom, we have explored this PES using a hierarchical series of ab initio methods (Hartree-Fock, second-order Møller-Plesset perturbation theory, fourth-order Møller-Plesset perturbation theory with single, double and quadruple excitations, coupled cluster theory with single and double excitations (CCSD), and with triple excitations treated perturbatively [CCSD(T)]) and hybrid density functional theory using the B3LYP functional, in combination with a hierarchical series of ten Gaussian-type basis sets, up to g polarization. Relativistic effects are taken into account either through a relativistic effective core potential for palladium or through a full four-component all-electron approach. Counterpoise corrected relative energies of stationary points are converged to within 0.1-0.2 kcal/mol as a function of the basis-set size. Our best estimate of kinetic and thermodynamic parameters is -8.1 (-8.3) kcal/mol for the formation of the reactant complex, 5.8 (3.1) kcal/mol for the activation energy relative to the separate reactants, and 0.8 (-1.2) kcal/mol for the reaction energy (zero-point vibrational energy-corrected values in parentheses). This agrees well with available experimental data. Our work highlights the importance of sufficient higher angular momentum polarization functions, f and g, for correctly describing metal-d-electron correlation and, thus, for obtaining reliable relative energies. We show that standard basis sets, such as LANL2DZ+ 1f for palladium, are not sufficiently polarized for this purpose and lead to erroneous CCSD(T) results. B3LYP is associated with smaller basis set superposition errors and shows faster convergence with basis-set size but yields relative energies (in particular, a reaction barrier) that are ca. 3.5 kcal/mol higher than the corresponding CCSD(T) values

Docking to heme proteins.

In silico screening has become a valuable tool in drug design, but some drug targets represent real challenges for docking algorithms. This is especially true for metalloproteins, whose interactions with ligands are difficult to parametrize. Our docking algorithm, EADock, is based on the CHARMM force field, which assures a physically sound scoring function and a good transferability to a wide range of systems, but also exhibits difficulties in case of some metalloproteins. Here, we consider the therapeutically important case of heme proteins featuring an iron core at the active site. Using a standard docking protocol, where the iron-ligand interaction is underestimated, we obtained a success rate of 28% for a test set of 50 heme-containing complexes with iron-ligand contact. By introducing Morse-like metal binding potentials (MMBP), which are fitted to reproduce density functional theory calculations, we are able to increase the success rate to 62%. The remaining failures are mainly due to specific ligand-water interactions in the X-ray structures. Testing of the MMBP on a second data set of non iron binders (14 cases) demonstrates that they do not introduce a spurious bias towards metal binding, which suggests that they may reliably be used also for cross-docking studies.

Pi*-pi* bonding interactions generated by halogen oxidation of zirconium(IV) redox-active ligand complexes.

The new complex, [Zr(pda)2]n (1, pda2- = N,N'-bis(neo-pentyl)-ortho-phenylenediamide, n = 1 or 2), prepared by the reaction of 2 equiv of pdaLi2 with ZrCl4, reacts rapidly with halogen oxidants to afford the new product ZrX2(disq)2 (3, X = Cl, Br, I; disq- = N,N'-bis(neo-pentyl)-ortho-diiminosemiquinonate) in which each redox-active ligand has been oxidized by one electron. The oxidation products 3a-c have been structurally characterized and display an unusual parallel stacked arrangement of the disq- ligands in the solid state, with a separation of approximately 3 A. Density functional calculations show a bonding-type interaction between the SOMOs of the disq- ligands to form a unique HOMO while the antibonding linear combination forms a unique LUMO. This orbital configuration leads to a closed-shell-singlet ground-state electron configuration (S = 0). Temperature-dependent magnetism measurements indicate a low-lying triplet excited state at approximately 750 cm-1. In solution, 3a-c show strong disq--based absorption bands that are invariant across the halide series. Taken together these spectroscopic measurements provide experimental values for the one- and two-electron energies that characterize the pi-stacked bonding interaction between the two disq- ligands.

Basis set effects on the energy and hardness profiles of the hydrogen fluoride dimer

In earlier work, the present authors have shown that hardness profiles are less dependent on the level of calculation than energy profiles for potential energy surfaces (PESs) having pathological behaviors. At variance with energy profiles, hardness profiles always show the correct number of stationary points. This characteristic has been used to indicate the existence of spurious stationary points on the PESs. In the present work, we apply this methodology to the hydrogen fluoride dimer, a classical difficult case for the density functional theory methods

Thermal nucleation of cavities in liquid helium at negative pressures

We have investigated the nucleation rate at which cavities are formed in 4He and 3He at negative pressures due to thermal fluctuations. To this end, we have used a density functional that reproduces the He liquid-gas interface along the coexistence line. The inclusion of thermal effects in the calculation of the barrier against nucleation results in a sizable decrease of the absolute value of the tensile strength above 1.5 K.

Response of doped 4He droplets

In the framework of a finite-range density-functional theory, we compute the response of 4HeN clusters doped with a rare-gas molecule. For this purpose, the mean field for the 4He atoms, their wave functions and effective quasiparticle interaction, are self-consistently calculated for a variety of particle numbers in the cluster. The response function is then evaluated for several multipolarities in each drop and the collective states are consequently located from the peaks of the strength function. The spectra of pure droplets approach those previously extracted with a similar algorithm resorting to a zero-range density functional. The spectra of doped clusters are sensitive to the presence of the impurity and are worth a future systematic investigation.

Multipole surface-plasmon modes on simple metals

The average multipole surface-plasmon energy for simple metals, as well as that of ordinary surface and bulk plasmons, is obtained using energy-weighted moments of the electronic response to sufficiently general external perturbations. A local approximation of exchange and correlation effects is used within a jellium model. Band-structure effects are incorporated through an effective electronic mass. Taking advantage of the transparency of the method, we analyze under what circumstances such modes might be observable. It is shown that due to an interplay between Coulomb and kinetic energies, the multipole modes become unobservable for increasing values of the transferred momentum (q) parallel to the surface. The value of q at which the multipole mode becomes unobservable is much smaller than the cutoff value for Landau damping. The effect of the electronic surface diffuseness is also analyzed. We compare our results with previous density-functional calculations and with recent experimental data for Na, K, and Cs.

Structure and energetics of mixed 4He-3He drops

Using a finite-range density functional, we have investigated the energetics and structural features of mixed helium clusters. The possibility of doping the cluster with a molecule of sulfur hexafluoride is also considered. It is seen that the repulsion introduced by the impurity strongly modifies the properties of the smallest drops. Although only a qualitative comparison is possible, the gross features displayed by our calculations are in agreement with recent experimental findings.

Structure and far-infrared edge modes of quantum antidots at zero magnetic field

We have investigated edge modes of different multipolarity sustained by quantum antidots at zero magnetic field. The ground state of the antidot is described within a local-density-functional formalism. Two sum rules, which are exact within this formalism, have been derived and used to evaluate the energy of edge collective modes as a function of the surface density and the size of the antidot.

Transverse dipole spin modes in quantum dots

We have carried out a systematic analysis of the transverse dipole spin response of a large-size quantum dot within time-dependent current density functional theory. Results for magnetic fields corresponding to integer filling factors are reported, as well as a comparison with the longitudinal dipole spin response. As in the two-dimensional electron gas, the spin response at high-spin magnetization is dominated by a low-energy transverse mode.

Ground-state self-consistent calculation of quantum dots under magnetic fields: Addition spectrum

A density-functional self-consistent calculation of the ground-state electronic density of quantum dots under an arbitrary magnetic field is performed. We consider a parabolic lateral confining potential. The addition energy, E(N+1)-E(N), where N is the number of electrons, is compared with experimental data and the different contributions to the energy are analyzed. The Hamiltonian is modeled by a density functional, which includes the exchange and correlation interactions and the local formation of Landau levels for different equilibrium spin populations. We obtain an analytical expression for the critical density under which spontaneous polarization, induced by the exchange interaction, takes place.

From nonwetting to prewetting: The asymptotic behavior of 4He drops on alkali substrates

We investigate the spreading of 4He droplets on alkali-metal surfaces at zero temperature, within the frame of finite range density-functional theory. The equilibrium configurations of several 4HeN clusters and their asymptotic trend with increasing particle number N, which can be traced to the wetting behavior of the quantum fluid, are examined for nanoscopic droplets. We discuss the size effects inferring that the asymptotic properties of large droplets correspond to those of the prewetting film.

Integer filling facfor phases and isospin in vertical diatomic artificial molecules

Integer filling factor phases of many-electron vertically coupled diatomic artificial quantum dot molecules are investigated for different values of the interdot coupling. The experimental results are analyzed within local-spin density functional theory for which we have determined a simple lateral confining potential law that can be scaled for the different coupling regimes, and Hartree-Fock theory. Maximum density droplets composed of electrons in both bonding and antibonding or just bonding states are revealed, and interesting isospin-flip physics appears for weak interdot coupling when the systematic depopulation of antibonding states leads to changes in isospin.