Performance of the M11 and M11-L density functionals for calculations of electronic excitation energies by adiabatic time-dependent density functional theory

Adiabatic time-dependent density functional theory is a powerful method for calculating electronic excitation energies of complex systems, but the quality of the results depends on the choice of approximate density functional. In this article we test two promising new density functionals, M11 and M11-L, against databases of 214 diverse electronic excitation energies, and we compare the results to those for 16 other density functionals of various kinds and to time-dependent Hartree-Fock. Charge transfer excitations are well known to be the hardest challenge for TDDFT. M11 is a long-range-corrected hybrid meta-GGA, and it shows better performance for charge transfer excitations than any of the other functionals except M06-HF, which is a specialized functional that does not do well for valence excitations. Several other long-range-corrected hybrid functionals also do well, and we especially recommend M11, ωB97X, and M06-2X for general spectroscopic applications because they do exceptionally well on ground-state properties as well as excitation energies. Local functionals are preferred for many applications to extended systems because of their significant cost advantage for large systems. M11-L is a dual-range local functional and-unlike all previous local functionals-it has good performance for Rydberg states as well as for valence states. Thus it is highly recommended for excitation energy calculations on extended systems.     

Improved description of nuclear magnetic resonance chemical shielding constants using the M06-L meta-generalized-gradient-approximation density functional

The performance of the M06-L density functional has been tested for four databases of NMR isotropic chemical shielding constants. Comparison with the B3LYP, BLYP, HCTH, KT1, KT2, LSDA, OPBE, OLYP, PBE, TPSS, and VSXC functionals shows that M06-L has improved performance for calculating NMR chemical shielding constants, especially for highly correlated systems. We also found that VSXC and M06-L have encouraging accuracy for calculating (13)C chemical shielding constants, and both functionals perform very well for the chemical shielding constants in the o-benzyne molecule.     

Physical basis for constrained lattice density functional theory

To study nucleation phenomena in an open system, a constrained lattice density functional theory (LDFT) method has been developed before to identify the unstable directions of grand potential functional and to stabilize nuclei by imposing a suitable constraint. In this work, we answer several questions about the method on a fundamental level, and give a firmer basis for the constrained LDFT method. First, we demonstrate that the nucleus structure and free energy barrier from a volume constraint method are equivalent to those from a surface constraint method. Then, we show that for the critical nucleus, the constrained LDFT method in fact produces a bias-free solution for both the nucleus structure and nucleation barrier. Finally, we give a physical interpretation of the Lagrange multiplier in the constraint method, which provides the generalized force to stabilize a nucleus in an open system. The Lagrange multiplier is found to consist of two parts: part I of the constraint produces an effective pressure, and part II imposes a constraint to counteract the supersaturation.     

Solid state density functional calculations for the group 11 monohalides

The group 11 monohalides are characterized by a variety of different solid-state structures and modifications. The copper halides crystallize mainly in a cubic zinc blende structure, while silver halides are mostly found in a rock salt modification. Completely different are the gold compounds where relativistic effects change the symmetry from a cubic to a chainlike AuX arrangement (X = F, Cl, Br, I) with short Au-Au internuclear distances. Here we present a systematic study of all solid state group 11 halides by scalar relativistic density functional theory for the experimentally known and observed structures, as well as for other unknown modifications and compare their relative stability.     

Performance of parallel TURBOMOLE for density functional calculations

The parallelization of density functional treatments of molecular electronic energy and first-order gradients is described, and the performance is documented. The quadrature required for exchange correlation terms and the treatment of exact Coulomb interaction scales virtually linearly up to 100 nodes. The RI-J technique to approximate Coulomb interactions (by means of an auxiliary basis set approximation for the electron density) even shows superlinear speedup on distributed memory architectures. The bottleneck is then linear algebra. Demonstrative application examples include molecules with up to 300 atoms and 3000 basis functions that can now be treated in a few hours per geometry optimization cycle in C₁ symmetry. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1746–1757, 1998     

Adsorption energy distribution from the Aranovich-Donohue lattice density functional theory

We propose a new methodology projected for the estimation of the adsorption energy distribution from the monolayer part of a single nitrogen adsorption isotherm determined at 77 K based on the lattice density functional theory (DFT) via the Aranovich-Donohue formalism. At first sight, the presented approach is computationally more difficult than a classical one. However, it is more flexible and comprehensible. Next, we developed a numerical program and used it for the estimation of the adsorption energy distribution from the experimental data on carbon black samples. The main nitrogen molecule-carbon black surface interaction energy can be estimated as approximately 7-8 kJ/mol, but the heterogeneity of the investigated materials differs significantly. Furthermore, we compare the results obtained from the lattice DFT via the Aranovich-Donohue formalism with the solution of the integral equation with the kernel represented by the classical monolayer localized Fowler-Guggenheim isotherm equation. The similarity between these two independent approaches is observed. The proposed methodology can be used for the investigation of the energetic heterogeneity of not only the carbonaceous materials but also the other "flat-surfaced" solids.     

A density functional model for the binary crystal of hard spheres with vacancies

We study the stability of a binary mixture of hard spheres in the crystalline state in which a small fraction of lattice sites in the solid structure are vacant. The optimum vacancy concentration is obtained by minimizing the free energy of the inhomogeneous solid state. We use the modified weighted density functional approximation to compute the free energy. The necessary input for the theory is the thermodynamic properties of the homogeneous state of the mixture and is obtained from the solutions of the corresponding Percus-Yevick integral equations for the binary system. We compute the free energy for the crystal having two kinds of ordered structures in which (i) both the species lie in a disordered manner on a single face-centered-cubic lattice and (ii) each of the two species lie on a separate cubic lattice. Our theoretical model obtains equilibrium vacancy fraction of O(10(-5)) near the freezing point in both cases. The vacancy concentration decreases with the increase of density in both cases.     

Performance of the Harris functional for extended basis sets at the Hartree-Fock and density functional levels

The Harris functional is a noniterative variational procedure that uses an input charge density to produce an energy that is surprisingly accurate compared to the converged Kohn-Sham self-consistent result. We adapted and generalized this functional for the Hartree-Fock closed- and open-shell cases as well as examined its use for hybrid density functional methods such as B3LYP. Analysis of extended basis set calculations shows that at the B3LYP level an input density formed from a double zeta + polarization orbital basis is accurate enough to reproduce the energy of triple zeta + double polarization + diffuse orbital basis. For large molecules this translates into a computational speed that can be an order of magnitude faster. In the case of Hartree-Fock calculations a "bootstrapping technique" that employs successive applications of the Harris functional can further reduce computational times while retaining sufficient accuracy.     

On the use of the dual process Langmuir model for predicting unary and binary isosteric heats of adsorption

Analytic expressions for unary and binary isosteric heats of adsorption as a function of the adsorbed phase loading were derived from the dual process Langmuir (DPL) model using the Clausius-Clapeyron equation. Unary isosteric heats of adsorption predicted from these expressions for several adsorbate-adsorbent systems were compared to values in the literature predicted from the well-accepted graphical approach using Toth and unilan models (Adsorption Equilibrium Data Handbook; Prentice Hall: NJ, 1989). Predictions from the DPL model were also compared to rare experimental unary and binary isosteric heats of adsorption in the literature for another adsorbate-adsorbent system. In all cases, very good agreement was obtained, showing that the DPL model can be used in adsorption process modeling for accurately predicting not only ideal and nonideal mixed-gas adsorption equilibria (Langmuir 2011, 27, 4700), but also unary and even binary isosteric heats of adsorption.     

Thermodynamics of symmetric dimers: lattice density functional theory predictions and simulations

A new lattice density functional theory (DFT) approach is proposed for symmetric dimers taking into account all possible configurations for molecules adjacent to a central dimer. Comparison with Monte Carlo simulations shows significant improvement of the proposed model compared to previously developed version of lattice DFT for dimers. It is shown that the new model gives accurate analytical solutions over a wide range of densities and temperatures. Phase transitions in dimers are analyzed and fundamental differences between dimers and monomers are discussed.     

Performance of dispersion-corrected density functional theory for the interactions in ionic liquids

Potential energy curves for the dissociation of cation-anion associates representing the building units of ionic liquids have been computed with dispersion corrected DFT methods. Non-local van der Waals density functionals (DFT-NL) for the first time as well as an atom pair-wise correction method (DFT-D3) have been tested. Reference data have been computed at the extrapolated MP2/CBS and estimated CCSD(T)/CBS levels of theory. The investigated systems are combined from two cations (1-butyl-3-methylimidazolium and tributyl(methyl)posphonium) and three anions (chloride, dicyanamide, acetate). We find substantial stabilization from London dispersion energy near equilibrium of 5-7 kcal mol(-1) (about 5-6% of the interaction energy). Equilibrium distances are shortened by 0.03-0.09 ? and fundamental (inter-fragment) vibrational frequencies (which are in the range 140-180 cm(-1)) are increased by typically 10 cm(-1) when dispersion corrections are made. The dispersion-corrected hybrid functional potentials are in general in excellent agreement with the corresponding CCSD(T) reference data (typical deviations of about 1-2%). The DFT-D3 method performs unexpectedly well presumably because of cancellation of errors between the dispersion coefficients of the cations and anions. Due to self-interaction error, semi-local density functionals exhibit severe SCF convergence problems, and provide artificial charge-transfer and inaccurate interaction energies for larger inter-fragment distances. Although these problems may be alleviated in condensed phase simulations by effective Coulomb screening, only dispersion-corrected hybrid functionals with larger amounts of Fock-exchange can in general be recommended for such ionic systems.     

Calculation of Verdet constants with time-dependent density functional theory: implementation and results for small molecules

We report the implementation of a method by which to calculate Verdet constants for molecules. The method is based on gauge-including atomic orbitals (GIAOs) and density functional theory. Calculations based on this method afford magneto-optical rotations of the right magnitude for the molecules H2, N2, CO, HF, CH4, C2H2, H2O, and CS2. The results are in satisfactory agreement with experiment. We investigate the dependency of the results on the gauge origin if GIAOs are not chosen, the convergence of the results with the size of the basis set for AOs and GIAOs, and for H2O and CS2 a comparison of gas-phase and liquid phase values. For the small molecules studied here, large polarized basis sets with diffuse functions are required to obtain well converged results. The use of an asymptotically correct Kohn-Sham potential is advantageous.     

Prediction of spectroscopic constants for diatomic molecules in the ground and excited states using time-dependent density functional theory

Spectroscopic constants of the ground and next seven low-lying excited states of diatomic molecules CO, N2, P2, and ScF were computed using the density functional theory SAOP/ATZP model, in conjunction with time-dependent density functional theory (TD-DFT) and a recently developed Slater type basis set, ATZP. Spectroscopic constants, including the equilibrium distances r(e), harmonic vibrational frequency omega(e), vibrational anharmonicity omega(e)x(e), rotational constant B(e), centrifugal distortion constant D(e), the vibration-rotation interaction constant alpha(e), and the vibrational zero-point energy E(n)0 were generated in an effort to establish a reliable database for electron spectroscopy. By comparison with experimental values and a similar model with an established larger Slater-type basis set, et-QZ3P-xD, it was found that this model provides reliably accurate results at reduced computational costs, for both the ground and excited states of the molecules. The over all errors of all eight lowest lying electronic states of the molecules under study using the effective basis set are r(e)(+/-4%), omega(e)(+/-5% mostly without exceeding +/-20%), omega(e)x(e)(+/-5% mostly without exceeding 20%, much more accurate than a previous study on this constant of +/-30%), B(e)(+/-8%), D(e)(+/-10%), alpha(e)(+/-10%), and E(n)0(+/-10%). The accuracy obtained using the ATZP basis set is very competitive to the larger et-QZ3P-xD basis set in particular in the ground electronic states. The overall errors in r(e), omega(e)x(e), and alpha(e) in the ground states were given by +/-0.7, +/-10.1, and +/-8.4%, respectively, using the efficient ATZP basis set, which is competitive to the errors of +/-0.5, +/-9.2, and +/-9.1%, respectively for those constants using the larger et-QZ3P-xD basis set. The latter basis set, however, needs approximately four times of the CPU time on the National Supercomputing Facilities (Australia). Due to the efficiency of the model (TD-DFT, SAOP and ATZP), it will be readily applied to study larger molecular systems.     

On the relation between orbital-localization and self-interaction errors in the density functional theory treatment of organic semiconductors

It is commonly argued that the self-interaction error (SIE) inherent in semilocal density functionals is related to the degree of the electronic localization. Yet at the same time there exists a latent ambiguity in the definitions of the terms "localization" and "self-interaction," which ultimately prevents a clear and readily accessible quantification of this relationship. This problem is particularly pressing for organic semiconductor molecules, in which delocalized molecular orbitals typically alternate with localized ones, thus leading to major distortions in the eigenvalue spectra. This paper discusses the relation between localization and SIEs in organic semiconductors in detail. Its findings provide further insights into the SIE in the orbital energies and yield a new perspective on the failure of self-interaction corrections that identify delocalized orbital densities with electrons.     

Solubility of hydrogen in PdAg and PdAu binary alloys using density functional theory

The present work deals with the study of palladium-silver (PdAg) and palladium-gold (PdAu) binary alloys over a broad range of temperatures and alloy compositions using density functional theory (DFT) to find possible conditions where the solubility of hydrogen (H) is significantly higher than that of pure palladium (Pd). Several alloy structures, such as Pd(100-x)Ag(x) with x = 14.81, 25.93, 37.04, and 48.51, Pd(100-x)Aux with x = 14.81, 25.93, and 37.04, and Pd(100-x)Cu(x) with x = 25.93 and 48.51 were considered. The lattice constants of these structures were optimized using DFT, and relaxed structures were used for the estimation of binding energy. It was found that the solubility of H in PdAg is higher than pure Pd with a maximum at approximately 30% Ag at 456 K. Also, the solubility of PdAu alloys was higher than pure Pd with a maximum at about 20% Au with a solubility 12 times higher than that of pure Pd. It was found that for a 3.7% H concentration in a PdAg alloy, a cell expansion of 0.15-0.2% occurs, which if ignored may affect the individual binding energy of the O-site by approximately 3.56% and may affect the predicted solubility by approximately 11.8%.     

Fundamental measure density functional theory studies on the freezing of binary hard-sphere and Lennard-Jones mixtures

Free energies and correlation functions of liquid and solid hard-sphere (HS) mixtures are calculated using the fundamental measure density functional theory. Using the thermodynamic perturbation theory the free energies of solid and liquid Lennard-Jones (LJ) mixtures are obtained from correlation functions of HS systems within a single theoretical approach. The resulting azeotrope- and spindle-type solid-liquid phase diagrams of HS and LJ binary mixtures are in good agreement with the corresponding ones from computer simulations.     

An improved and broadly accurate local approximation to the exchange-correlation density functional: The MN12-L functional for electronic structure calculations in chemistry and physics

We report a new local exchange-correlation energy functional that has significantly improved across-the-board performance, including main-group and transition metal chemistry and solid-state physics, especially atomization energies, ionization potentials, barrier heights, noncovalent interactions, isomerization energies of large moleucles, and solid-state lattice constants and cohesive energies.     

Performance of density functional theory on homogeneous gold catalysis

The performance of 32 density functionals to describe the homogeneous gold catalysis of propargyl esters has been tested. These catalytic reactions are very commonly rather complex, numerous intermediates can be found along the reaction profile and the individual reaction steps are often associated to very small barrier heights. In this scenario, the experimental kinetic study of the catalytic mechanisms is very challenging. A computational approach to this problem provides invaluable help to gain insight into these mechanisms. However the high accuracy needed to describe such highly branched paths with low energy transition states poses many practical problems for cost-efficient DFT methods. The lack of accurate experimental or high-level computational data to employ as validation sets for these methods is also an added difficulty. High-level computational data needed to validate these functionals have been obtained at the CCSD/def2-TZVPP//CCSD/def2-SVP level and, with such data, we aim to help discern the best density functional recipe for homogeneous gold catalysis     

New approach to the kinetics of heterogeneous unary nucleation on liquid aerosols of a binary solution

The formation of a droplet on a hygroscopic center may occur either in a barrierless way via Kohler activation or via nucleation by overcoming a free energy barrier. Unlike the former, the latter mechanism of this process has been studied very little and only in the framework of the classical nucleation theory based on the capillarity approximation whereby a nucleating droplet behaves like a bulk liquid. In this paper the authors apply another approach to the kinetics of heterogeneous nucleation on liquid binary aerosols, based on a first passage time analysis which avoids the concept of surface tension for tiny droplets involved in nucleation. Liquid aerosols of a binary solution containing a nonvolatile solute are considered. In addition to modeling aerosols formed through the deliquescence of solid soluble particles, the considered aerosols constitute a rough model of "processed" marine aerosols. The theoretical results are illustrated by numerical calculations for the condensation of water vapor on binary aqueous aerosols with nonvolatile nondissociating solute molecules using Lennard-Jones potentials for the molecular interactions.