Shawn T Brown | PubMed Health

PubMed2006

Advances in methods and algorithms in a modern quantum chemistry program package.

Advances in theory and algorithms for electronic structure calculations must be incorporated into program packages to enable them to become routinely used by the broader chemical community. This work reviews advances made over the past five years or so that constitute the major improvements contained in a new release of the Q-Chem quantum chemistry package, together with illustrative timings and applications. Specific developments discussed include fast methods for density functional theory calculations, linear scaling evaluation of energies, NMR chemical shifts and electric properties, fast auxiliary basis function methods for correlated energies and gradients, equation-of-motion coupled cluster methods for ground and excited states, geminal wavefunctions, embedding methods and techniques for exploring potential energy surfaces.

Algorithms↗

PubMed2006

Efficient computation of the exchange-correlation contribution in the density functional theory through multiresolution.

A new algorithm is presented to improve the efficiency of the computation of exchange-correlation contributions, a major time-consuming step in a density functional theory (DFT) calculation. The new method, called multiresolution exchange correlation (mrXC), takes advantage of the variation in resolution among the Gaussian basis functions and shifts the calculation associated with low-resolution (smooth) basis function pairs to an even-spaced cubic grid. The cubic grid is much less dense in the vicinity of the nuclei than the atom-centered grid and the computation on the former is shown to be much more efficient than on the latter. MrXC does not alter the formalism of the current standard algorithm based on the atom-centered grid (ACG), but instead employs two fast and accurate transformations between the ACG and the cubic grid. Preliminary results with local density approximation have shown that mrXC yields three to five times improvement in efficiency with negligible error. The extension to DFT functionals with generalized gradient approximation is also briefly discussed.

Algorithms↗

PubMed2003

Cyclopentadiene annulated polycyclic aromatic hydrocarbons: investigations of electron affinities.

The adiabatic electron affinities of cyclopentadiene and 10 associated benzannelated derivatives have been predicted with both density functional and Hartree-Fock theory. These systems can also be regarded as benzenoid polycyclic aromatic hydrocarbons (PAHs) augmented with five-membered rings. Like the PAHs, the electron affinities of the present systems generally increase with the number of rings. To unequivocally bind an electron, cyclopentadiene must have at least two conventionally fused benzene rings. 1H-Benz[f]indene, a naphthalene-annulated cyclopentadiene, is predicted to have a zero-point energy corrected adiabatic electron affinity of 0.13 eV. Since the experimental E(A) of naphthalene is negative (-0.19 eV), the five-membered ring appendage contributes to the stability of the naphthalene-derived 1H-benz[f]indene radical anion significantly. The key to binding the electron is a contiguous sequence of fused benzenes, since fluorene, the isomer of 1H-benz[f]indene, with separated six-membered rings, has an electron affinity of -0.07 eV. Each additional benzene ring in the sequence fused to cyclopentadiene increases the electron affinity by 0.15-0.65 eV: the most reliable predictions are cyclopentadiene (-0.63 eV), indene (-0.49 eV), fluorene (-0.07 eV), 1H-benz[f]indene (0.13 eV), 1,2-benzofluorene (0.25 eV), 2,3-benzofluorene (0.26 eV), 12H-dibenzo[b,h]fluorene (0.65 eV), 13H-indeno[1,2-b]anthracene (0.82 eV), and 1H-cyclopenta[b]naphthacene (1.10 eV). In contrast, if the six-membered ring-fusion is across the C(2)-C(3) cyclopentadiene single bond, only a single benzene is needed to bind an electron. The theoretical electron affinity of the resulting molecule, isoindene, is 0.49 eV, and this increases to 1.22 eV for 2H-benz[f]indene. The degree of aromaticity is responsible for this behavior. While the radical anions are stabilized by conjugation, which increases with the size of the system, the regular indenes, like PAHs in general, suffer from the loss of aromatic stabilization in forming their radical anions. While indene is 21 kcal mol(-1) more stable than isoindene, the corresponding radical anion isomers have almost the same energy. Nucleus-independent chemical shift calculations show that the highly aromatic molecules lose almost all aromaticity when an extra electron is present. The radical anions of cyclopentadiene and all of its annulated derivatives have remarkably low C-H bond dissociation energies (only 18-34 kcal mol(-1) for the mono-, bi-, and tricyclics considered). Hydrogen atom loss leads to the restoration of aromaticity in the highly stabilized cyclopentadienyl anion congeners.

Journal Article↗