A density functional especially designed for hydrogen-only systems

Comparative Density Functional Theory Study of Magnetic Exchange Couplings in Dinuclear Transition-Metal Complexes

Multicenter transition-metal complexes (MCTMs) with magnetically interacting ions have been proposed as components for information-processing devices and storage units. For any practical application of MCTMs as magnetic units, it is crucial to characterize their magnetic behavior, and in particular, the isotropic magnetic exchange coupling, J, between its magnetic centers. Due to the large size of typical MCTMs, density functional theory is the only practical electronic structure method for evaluating the J coupling. Here, we assess the accuracy of different density functional approximations for predicting the magnetic couplings of eight dinuclear transition-metal complexes, including five dimanganese, two dicopper, and one divanadium with known reliable experimental J couplings spanning from ferromagnetic to strong antiferromagnetic. The density functionals considered include global hybrid functionals which mix semilocal density functional approximations and exact exchange with a fixed admixing parameter, six local hybrid functionals where the admixing parameters are extended to be spatially dependent, the SCAN and r2SCAN meta-generalized gradient approximations (GGAs), and two widely used GGAs. We found that global hybrids tested in this work have a tendency to over-correct the error in magnetic coupling parameters from the Perdew-Burke-Ernzerhof (PBE) GGA as seen for manganese complexes. The performance of local hybrid density functionals shows no improvement in terms of bias and is scattered without a clear trend, suggesting that more efforts are needed for the extension from global to local hybrid density functionals for this particular property. The SCAN and r2SCAN meta-GGAs are found to perform as well as benchmark global hybrids on most tested complexes. We further analyze the charge density redistribution of meta-GGAs as well as global and local hybrid density functionals with respect to that of PBE, in connection to the self-interaction error or delocalization error.

The smallest proton-bound dimer H5+: theoretical progress

The protonated hydrogen dimer, H5+, is the smallest system including proton transfer, and has been of long-standing interest since its first laboratory observation in 1962. H5+ and its isotopologues are the intermediate complexes in deuterium fractionation reactions, and are of central importance in molecular astrophysics. The recently recorded infrared spectra of both H5+ and D5+ reveal a rich vibrational dynamics of the cations, which presents a challenge for standard theoretical approaches. Although H5+ is a four-electron ion, which makes highly accurate electronic structure calculations tractable, the construction of ab initio-based potential energy and dipole moment surfaces has proved a hard task. In the same vein, the difficulties in treating the nuclear motion could also become cumbersome due to their high dimensionality, floppiness and/or symmetry. These systems are prototypical examples for studying large-amplitude motions, as they are highly delocalized, interconverting between equivalent minima through internal rotation and proton transfer motions requiring state-of-the-art treatments. Recent advances in the computational vibrational spectroscopy of the H5+ cation and its isotopologues are reported from full quantum spectral simulations, providing important information in a rigorous manner, and open perspectives for further future investigations. This article is part of a discussion meeting issue 'Advances in hydrogen molecular ions: H3+, H5+ and beyond'.

Coupled-Perturbed Scheme for the Calculation of Electronic g-Tensors with Local Hybrid Functionals

A coupled-perturbed Kohn-Sham (CPKS) scheme for calculating second-order magnetic properties has been developed for the case of general occupied-orbital-dependent (OOD) exchange-correlation functionals involving the exact-exchange energy density. The origin of the coupling terms in the functional derivatives of OOD functionals with respect to the orbitals has been thoroughly analyzed, and general expressions for the resulting coupling terms have been obtained. The generalized CPKS scheme thus obtained has been implemented within the MAG-ReSpect code and tested in calculations of electronic g-tensors with local hybrid functionals. Compared to previously tested global hybrids, like B3LYP, thermochemically optimized local hybrids provide only little to moderate improvement for test sets of main-group radicals and paramagnetic transition-metal complexes. Closer analyses point to possible areas in which the fundamentally more flexible local hybrids may be improved for the property at hand.

First-principles simulations of vibrational states and spectra for H5(+) and D5(+) clusters using multiconfiguration time-dependent Hartree approach

Simulations of the infrared (IR) spectra of the H5(+) and D5(+) clusters are carried out in the whole energy range, using a recent, reliable "on the fly" DFT-based potential energy surface, and its corresponding dipole moment surface. For the present study we adopted a recently proposed four-dimensional quantum model to describe the proton transfer motion between the two vibrating H2 or D2 units. Time-dependent and time-independent approaches within the multiconfiguration time-dependent Hartree method are employed for investigating the vibrational dynamics of the complexes. The obtained spectra are compared with recent experimental data available for energies up to 4500 and 3500 cm(-1) for the H5(+) and D5(+), respectively. Even though the present results are based on a reduced dimensional model, the infrared spectra are shown to be in good qualitative accord with those observed experimentally. Also as the reported data are subject to the potential energy surface, comparisons with previous theoretical calculations based on an analytical ab initio parameterized surface are also presented. The differences on the topology of the potentials are discussed in connection with their effect on the spectral features. We found that the main characteristics of the experimentally observed spectra are reproduced by both surfaces, evaluating in this way the sensitivity of such computations on the quality of the underlying potential. This finding serves to connect aspects of the potential surface of these systems to their spectral complexity, and could be indicative to calibrate intrinsic errors in their calculation for future studies.

Full-dimensional quantum calculations of the dissociation energy, zero-point, and 10 K properties of H7+/D7+ clusters using an ab initio potential energy surface

Diffusion Monte Carlo (DMC) and path-integral Monte Carlo computations of the vibrational ground state and 10 K equilibrium state properties of the H7 (+)/D7 (+) cations are presented, using an ab initio full-dimensional potential energy surface. The DMC zero-point energies of dissociated fragments H5 (+)(D5 (+))+H2(D2) are also calculated and from these results and the electronic dissociation energy, dissociation energies, D0, of 752 ± 15 and 980 ± 14 cm(-1) are reported for H7 (+) and D7 (+), respectively. Due to the known error in the electronic dissociation energy of the potential surface, these quantities are underestimated by roughly 65 cm(-1). These values are rigorously determined for first time, and compared with previous theoretical estimates from electronic structure calculations using standard harmonic analysis, and available experimental measurements. Probability density distributions are also computed for the ground vibrational and 10 K state of H7 (+) and D7 (+). These are qualitatively described as a central H3 (+)/D3 (+) core surrounded by "solvent" H2/D2 molecules that nearly freely rotate.

Vibrational dynamics of the H5(+) and its isotopologues from multiconfiguration time-dependent Hartree calculations

Full-dimensional multiconfiguration time-dependent Hartree (MCTDH) computations are reported for the vibrational states of the H(5)(+) and its H(4)D(+), H(3)D(2)(+), H(2)D(3)(+), HD(4)(+), D(5)(+) isotopologues employing two recent analytical potential energy surfaces of Xie et al. [J. Chem. Phys. 122, 224307 (2005)] and Aguado et al. [J. Chem. Phys. 133, 024306 (2010)]. The potential energy operators are constructed using the n-mode representation adapted to a four-combined mode cluster expansion, including up to seven-dimensional grids, chosen adequately to take advantage in representing the MCTDH wavefunction. An error analysis is performed to quantify the convergence of the potential expansion to reproduce the reference surfaces at the energies of interest. An extensive analysis of the vibrational ground state properties of these isotopes and comparison with the reference diffusion Monte Carlo results by Acioli et al. [J. Chem. Phys. 128, 104318 (2008)] are presented. It is found that these systems are highly delocalized, interconverting between equivalent minima through rotation and internal proton transfer motions even at their vibrational ground state. Isotopic substitution affects the zero-point energy and structure, showing preference in the arrangements of the H and D within the mixed clusters, and the most stable conformers of each isotopomer are the ones with the H in the central position. Vibrational excited states are also computed and by comparing the energies and structures predicted from the two surfaces, the effect of the potential topology on them is discussed.

Full-dimensional (15-dimensional) ab initio analytical potential energy surface for the H7+ cluster

Full-dimensional ab initio potential energy surface is constructed for the H(7)(+) cluster. The surface is a fit to roughly 160,000 interaction energies obtained with second-order MöllerPlesset perturbation theory and the cc-pVQZ basis set, using the invariant polynomial method [B. J. Braams and J. M. Bowman, Int. Rev. Phys. Chem. 28, 577 (2009)]. We employ permutationally invariant basis functions in Morse-type variables for all the internuclear distances to incorporate permutational symmetry with respect to interchange of H atoms into the representation of the surface. We describe how different configurations are selected in order to create the database of the interaction energies for the linear least squares fitting procedure. The root-mean-square error of the fit is 170 cm(-1) for the entire data set. The surface dissociates correctly to the H(5)(+) + H(2) fragments. A detailed analysis of its topology, as well as comparison with additional ab initio calculations, including harmonic frequencies, verify the quality and accuracy of the parameterized potential. This is the first attempt to present an analytical representation of the 15-dimensional surface of the H(7)(+) cluster for carrying out dynamics studies.

DFT and Proton Transfer Reactions: A Benchmark Study on Structure and Kinetics

A significant number of different exchange correlation functionals, ranging from generalized gradient approximations to double hybrids, has been tested on a difficult playground represented by proton transfer reactions. In order to have a complete picture of their performances, both energetics and structural features have been compared and the obtained ranking compared with those issued from the standard test for kinetics (i.e., the DBH24/08 set). Among all of the functionals, the ωB97X, BMK, B1LYP, and PBE0-DH approaches are those providing a good error balance on all four trials. Beyond these figures, the obtained results allow for some general considerations, such as those on the role of Hartree-Fock exchange in reaction barriers or the relation between structure and energetics.

A quantum reaction dynamics study of the translational, vibrational, and rotational motion effects on the HD + H3+ reaction

Time-dependent, quantum reaction dynamics wavepacket approach is employed to investigate the impacts of the translational, vibrational, and rotational motion on the HD+H(3)(+) → H(2)D(+) + H(2) reaction using the Xie-Braams-Bowman potential energy surface [Z. Xie, B. J. Braams, and J. M. Bowman, J. Chem. Phys. 122, 224307 (2005)]. We treat this five atom reaction with a seven-degree-of-freedom model by fixing one Jacobi and one torsion angle related to H(3) (+) at the lowest saddle point geometry of the potential energy surface. The initial state selected reaction probabilities show that the rotational excitations of H(+)-H(2) greatly enhance the reactivity with the reaction probabilities increased double at high rotational states compared to the ground state. However, the vibrational excitations of H(3) (+) hinder the reactivity. The ground state reaction probability shows no reaction threshold for this exoergic reaction, and as the translational energy increases, the reaction probability decreases. Furthermore, reactive resonances and zero point energy play very important roles on the reaction dynamics. The obtained integral cross section has the character of an exoergic reaction without a threshold: it decreases with the translational energy increasing. The calculated thermal rate constants using this seven-degree-of-freedom model are in agreement with a later experiment measurement.

Internal proton transfer and H2 rotations in the H5(+) cluster: a marked influence on its thermal equilibrium state

Classical and path integral Monte Carlo (CMC, PIMC) "on the fly" calculations are carried out to investigate anharmonic quantum effects on the thermal equilibrium structure of the H5(+) cluster. The idea to follow in our computations is based on using a combination of the above-mentioned nuclear classical and quantum statistical methods, and first-principles density functional (DFT) electronic structure calculations. The interaction energies are computed within the DFT framework using the B3(H) hybrid functional, specially designed for hydrogen-only systems. The global minimum of the potential is predicted to be a nonplanar configuration of C(2v) symmetry, while the next three low-lying stationary points on the surface correspond to extremely low-energy barriers for the internal proton transfer and to the rotation of the H2 molecules, around the C2 axis of H5(+), connecting the symmetric C(2v) minima in the planar and nonplanar orientations. On the basis of full-dimensional converged PIMC calculations, results on the quantum vibrational zero-point energy (ZPE) and state of H5(+) are reported at a low temperature of 10 K, and the influence of the above-mentioned topological features of the surface on its probability distributions is clearly demonstrated.

Toward a realistic density functional theory potential energy surface for the H5+ cluster

The potential energy surface of H(5)(+) is characterized using density functional theory. The hypersurface is evaluated at selected configurations employing different functionals, and compared with results obtained from ab initio CCSD(T) calculations. The lowest ten stationary points (minima and saddle-points) on the surface are located, and the features of the short-, intermediate-, and long-range intermolecular interactions are also investigated. A detailed analysis of the surface's topology, and comparisons with extensive CCSD(T) results, as well as a recent ab initio analytical surface, shows that density functional theory calculations using the B3(H) functional represent very well all aspects studied on the H(5)(+) potential. These include the tiny energy difference between the minimum at 1-C(2v) configuration and the 2-D(2d) one corresponding to the transition state for the proton transfer between the two equivalent C(2v) minima, and also the correct asymptotic behavior of the long-range interactions. The calculated binding energy and dissociation enthalpies compare very well with previous benchmark coupled-cluster ab initio data, and with experimental data available. Based on these results the use of such approach to perform first-principles molecular dynamics simulations could provide reliable information regarding the dynamics of protonated hydrogen clusters.

A thermochemically competitive local hybrid functional without gradient corrections

Following the suggestion of local hybrid functionals with position-dependent exact-exchange admixture [J. Jaramillo, G. E. Scuseria, and M. Ernzerhof, J. Chem. Phys. 118, 1068 (2003)], a functional that mixes only local and exact exchange plus local correlation has been constructed. With a simple local mixing function for the position dependence, this Lh-SVWN functional provides atomization energies for the G2-1 set that are competitive with currently available state-of-the-art functionals like, e.g., B3LYP. This is achieved without generalized gradient approximations for exchange or correlation.

Benchmarking approximate density functional theory. I.s/d excitation energies in 3d transition metal cations

The performance of a number of different implementations of density functional theory (DFT) for predicting the s/d interconfigurational energies of the 3d transition metal cations is investigated. Systematic comparisons of computed results with experimental data indicate that gradient corrected correlation functionals, like the LYP GGA, efficiently correct the flaws of the LDA, but reveal shortcomings in the treatment of exchange by currently available GGAs. The admixture of exact exchange in hybrid functionals eventually leads to largely reduced errors. Several basis sets available for the 3d elements are tested in combination with the B3LYP functional. Finally, the influence of variations of the admixture of exact exchange is systematically tested. The results reveal that computed s/d excitation energies obtained for the individual ions depend in markedly different ways on the amount of exact exchange admixture and that there is no single optimal and transferable exchange parameter to create a hybrid functional that yields improved results for all ions alike. Several of the recently devised functionals perform as good as or slightly better than the B3LYP functional in the present study. But given the fact that the B3LYP functional has been identified as the most successful DFT method in an overwhelming number of systematic investigations in very many areas of chemical research, there is no persuasive motivation to recommend its replacement by one of the other functionals, as much less is known about their robustness.