Search PubMed⌕ Search

Biomedical subjects

David J Tozer

Publications and source records attributed to David J Tozer.

9 recordsLinked to original sources

Influence of Coulomb-attenuation on exchange-correlation functional quality.

The dependence of functional quality on the attenuation parameters--which control the limiting (r12-->0, infinity) values and the rate of attenuation--is investigated for a Coulomb-attenuated exchange-correlation functional. For the attenuation and functional form considered, satisfaction of an exact long-range condition is detrimental for properties such as atomisation energies and bond lengths, but does improve classical reaction barriers and small molecule electronic excitation energies. The functionals considered can provide high quality valence, Rydberg, intramolecular and asymptotic intermolecular charge transfer (CT) excitations, but none are able to provide a simultaneously optimal description of all classes; CT excitations are not necessarily improved compared to those from conventional functionals. The study highlights the need for further development of Coulomb-attenuated functionals.

Journal Article↗

Assessment of a Coulomb-attenuated exchange-correlation energy functional.

The recently proposed CAM-B3LYP exchange-correlation energy functional, based on a partitioning of the r operator in the exchange interaction into long- and short-range components, is assessed for the determination of molecular thermochemistry, structures, and second order response properties. Rydberg and charge transfer excitation energies and static electronic polarisabilities are notably improved over the standard B3LYP functional; classical reaction barriers also improve. Ionisation potentials, bond lengths, NMR shielding constants and indirect spin-spin coupling constants are comparable with the two functionals. CAM-B3LYP atomisation energies and diatomic harmonic vibrational wavenumbers are less accurate than those of B3LYP. Future research directions are outlined.

Journal Article↗

Computation of the hardness and the problem of negative electron affinities in density functional theory.

The absolute hardness in density functional theory (DFT) is discussed, emphasizing the charge-transfer excitation interpretation. Direct evaluation from the computed ionization potential and electron affinity is intrinsically problematic when the affinity is negative; the calculated affinity exhibits a strong basis set dependence, becoming near zero as diffuse functions are added. An alternative Koopmans-based approximation using local functional eigenvalues uniformly and significantly underestimates the hardness. A simple correction to the Koopmans expression is highlighted on the basis of a consideration of the integer discontinuity. The resulting hardness expression does not require the explicit computation of the affinity and has a straightforward interpretation in terms of the electronegativity. The correction eliminates the underestimation and gives hardness values that do not degrade as the electron affinity becomes more negative. For systems with large negative affinities, the values are an improvement over those from the other approaches. The success can be traced to an implicit, unconventional approximation for the electron affinity, which outperforms the standard approach when the affinity is significantly negative and which does not break down as the basis set becomes more diffuse.

Journal Article↗

Semiempirical hybrid functional with improved performance in an extensive chemical assessment.

It is demonstrated that there is still scope for improvement in the quality of conventional, semiempirical hybrid exchange-correlation functionals in density-functional theory. A new functional, denoted B97-3, is determined from a fit to eight chemical properties (316 data points). For a series of 25 chemical assessments (850 data points) including 17 assessments and 10 chemical properties absent from the fitting data, B97-3 provides the lowest or joint-lowest mean absolute error on 15 occasions, compared to 6, 5, and 4 occasions for B3LYP, PBE0, and B97-2, respectively [A. D. Becke, J. Chem. Phys. 98, 5648 (1993); M. Ernzerhof and G. E. Scuseria, J. Chem. Phys. 110, 5029 (1999); C. Adamo and V. Barone, J. Chem. Phys. 110, 6158 (1999); P. J. Wilson, T. J. Bradley, and D. J. Tozer, J. Chem. Phys. 115, 9233 (2001)]. Mean absolute errors from B97-3 are, on average, 21%, 18%, and 12% smaller than from these three functionals. The most notable improvements are obtained for classical reaction barriers, where the error reductions are 60%, 54%, and 27%.

Journal Article↗

Exchange methods in Kohn-Sham theory.

Differences between exchange methods in exchange-only Kohn-Sham theory are highlighted by calculations of diatomic molecule total energies, uncoupled isotropic NMR shieldings, and HOMO-LUMO eigenvalue differences. Optimised effective potential (OEP) and Wu-Yang (WY) results are very similar. Localised Hartree-Fock (LHF) and Krieger-Li-Iafrate (KLI) results are close to one another, but are different to OEP and WY. Becke 1988 exchange (B88X) is different again. Shieldings reduce from OEP/WY to LHF/KLI to B88X, which is consistent with an observed reduction in HOMO-LUMO gaps. LHF, KLI, and B88X shieldings and HOMO-LUMO gaps are closer to near-exact, correlated values, than are the OEP values. These variations arise entirely due to differences in the one-electron exchange potentials, which is clearly evident in potential difference plots, relative to OEP, for the N2 molecule. Density difference plots are also presented, which exhibit a spatial correlation with the potential differences. HOMO and LUMO probability density difference plots show a contraction of the LUMO relative to OEP, which is consistent with the NMR and HOMO-LUMO findings. Plots are also presented for near-exact, correlated Kohn-Sham calculations. The features are qualitatively similar to those observed in the LHF, KLI, and B88X plots, highlighting correlated character in these approximate exchange-only calculations.

Algorithms↗

Ground- and excited-state diatomic bond lengths, vibrational levels, and potential-energy curves from conventional and localized Hartree-Fock-based density-functional theory.

Ground- and excited-state diatomic bond lengths, vibrational levels, and potential-energy curves are determined using conventional and localized Hartree-Fock (LHF)-based density-functional theory. Exchange only and hybrid functionals (with various fractions of exchange) are considered, together with a standard generalized gradient approximation (GGA). Ground-state bond lengths and vibrational wave numbers are relatively insensitive to whether orbital exchange is treated using the conventional or LHF approach. Excited-state calculations are much more sensitive. For a standard fraction of orbital exchange, N2 and CO vertical excitation energies at experimental bond lengths are accurately described by both conventional and LHF-based approaches, providing an asymptotic correction is present. Excited-state bond lengths and vibrational levels are more accurate with the conventional approach. The best quality, however, is obtained with an asymptotically corrected GGA functional. For the ground and lowest four singlet excited states, the GGA mean absolute errors in bond lengths are 0.006 A (0.5%) and 0.011 A (0.8%) for N2 and CO, respectively. Mean absolute errors in fundamental vibrational wavenumbers are 49 cm(-1) (2.7%) and 68 cm(-1) (5.0%), respectively. The GGA potential-energy curves are compared with near-exact Rydberg-Klein-Rees curves. Agreement is very good for the ground and first excited state, but deteriorates for the higher states.

Journal Article↗

A semiempirical generalized gradient approximation exchange-correlation functional.

We describe our attempts to improve upon the quality of the KT1 and KT2 generalized gradient approximation (GGA) exchange-correlation functionals [T. W. Keal and D. J. Tozer, J. Chem. Phys. 119, 3015 (2003)], through the introduction of additional gradient-corrected exchange and correlation terms. A GGA functional, denoted KT3, is presented, which maintains the high quality main-group nuclear magnetic resonance shielding constants obtained with KT1 and KT2; results are 2-3 times more accurate than conventional GGA and hybrid functionals. For the extensive range of systems considered in this study, KT3 also provides atomization energies, ionization potentials, electron affinities, proton affinities, bond angles, and electronic polarizabilities that are comparable to, or that surpass, those of the best present-day GGAs. Furthermore, it provides equilibrium molecular bond lengths and diatomic harmonic vibrational wave numbers that are as accurate as those from the best hybrid functionals. Further improvements are required in the description of classical chemical reaction barriers.

Journal Article↗

The observation of a large gauche preference when 2-fluoroethylamine and 2-fluoroethanol become protonated.

The energies of the gauche and anti conformers of 2-fluoroethylamine, 2-fluoroethanol and their protonated analogues are calculated using density functional theory. Unlike the non protonated systems, the protonated systems show a strong gauche effect where the C-F and the C-(+)NH(3) or C-F and C-(+)OH(2) bonds are gauche rather than anti to each other. Single crystal X-ray diffraction studies of 2-fluoroethylammonium compounds identify the same conformational preference.

Journal Article↗