Properties of atoms in electronically excited molecules within the formalism of TDDFT.

The topological analysis of the electron density for electronic excited states under the formalism of the quantum theory of atoms in molecules using time-dependent density functional theory (TDDFT) is presented. Relaxed electron densities for electronic excited states are computed by solving a Z-vector equation which is obtained by means of the Sternheimer interchange method. This is in contrast to previous work in which the electron density for excited states is obtained using DFT instead of TDDFT, that is, through the imposition of molecular occupancies in accordance with the electron configuration of the excited state under consideration. Once the electron density of the excited state is computed, its topological characterization and the properties of the atoms in molecules are obtained in the same manner that for the ground state. The analysis of the low-lying π→π* singlet and triplet vertical excitations of CO and C6H6 are used as representative examples of the application of this methodology. Altogether, it is shown how this procedure provides insights on the changes of the electron density following photoexcitation and it is our hope that it will be useful in the study of different photophysical and photochemical processes.

Subshell fitting of relativistic atomic core electron densities for use in QTAIM analyses of ECP-based wave functions.

Scalar-relativistic, all-electron density functional theory (DFT) calculations were done for free, neutral atoms of all elements of the periodic table using the universal Gaussian basis set. Each core, closed-subshell contribution to a total atomic electron density distribution was separately fitted to a spherical electron density function: a linear combination of s-type Gaussian functions. The resulting core subshell electron densities are useful for systematically and compactly approximating total core electron densities of atoms in molecules, for any atomic core defined in terms of closed subshells. When used to augment the electron density from a wave function based on a calculation using effective core potentials (ECPs) in the Hamiltonian, the atomic core electron densities are sufficient to restore the otherwise-absent electron density maxima at the nuclear positions and eliminate spurious critical points in the neighborhood of the atom, thus enabling quantum theory of atoms in molecules (QTAIM) analyses to be done in the neighborhoods of atoms for which ECPs were used. Comparison of results from QTAIM analyses with all-electron, relativistic and nonrelativistic molecular wave functions validates the use of the atomic core electron densities for augmenting electron densities from ECP-based wave functions. For an atom in a molecule for which a small-core or medium-core ECPs is used, simply representing the core using a simplistic, tightly localized electron density function is actually sufficient to obtain a correct electron density topology and perform QTAIM analyses to obtain at least semiquantitatively meaningful results, but this is often not true when a large-core ECP is used. Comparison of QTAIM results from augmenting ECP-based molecular wave functions with the realistic atomic core electron densities presented here versus augmenting with the limiting case of tight core densities may be useful for diagnosing the reliability of large-core ECP models in particular cases. For molecules containing atoms of any elements of the periodic table, the production of extended wave function files that include the appropriate atomic core densities for ECP-based calculations, and the use of these wave functions for QTAIM analyses, has been automated.

Atomic partitioning of the dissociation energy of the P-O(H) bond in hydrogen phosphate anion (HPO4(2-)): disentangling the effect of Mg2+.

This paper has three goals: (1) to provide a first step in understanding the atomic basis of the role of magnesium in facilitating the dissociation of the P-O bond in phosphorylated biochemical fuel molecules (such as ATP or GTP), (2) to compare second-order Møller-Plesset perturbation theory (MP2) results with those obtained at the more economical density functional theory (DFT) level for a future study of larger more realistic models of ATP/GTP, and (3) to examine the calculation of atomic total energies from atomic kinetic energies within a Kohn-Sham implemention of DFT, as compared to ab initio methods. A newly described method based on the quantum theory of atoms in molecules (QTAIM), which is termed the "atomic partitioning of the bond dissociation energy" (APBDE), is applied to a simple model of phosphorylated biological molecules (HPO42-). The APBDE approach is applied in the presence and in the absence of magnesium. It is found that the P-O(H) bond in the magnesium complex is shorter, exhibits a higher stretching frequency, and has a higher electron density at the bond critical point than in the magnesium-free hydrogen phosphate anion. Though these data would seem to suggest a stronger P-O(H) bond in the magnesium complex compared to the magnesium-free case, the homolytic breaking of the P-O(H) bond in the complex is found to be easier, i.e., has a lower BDE. This effect is the result of the balance of several atomic contributions to the BDE induced by the magnesium cation, which stabilizes the dissociation product more than it stabilizes the intact model molecule.

NMR chemical shifts. Substituted acetylenes.

The MPW1PW91/6-311+G(2d,p) and MP2/6-311+G(2d,p) GIAO nuclear shieldings for a series of monosubstituted acetylenes have been calculated using the MP2/6-311G(2d,p) geometries. Axially symmetric substituents such as fluorine may lead to large changes in the isotropic shielding but have little effect on the tensor component (zz) about the C[triple bond]C bond axis. On the other hand, substituents such as vinyl and aldehyde groups lead to essentially no difference in the isotropic shielding but are calculated to give a large zz paramagnetic shift to the terminal carbon of the acetylene group, without having much effect on the inner carbon. The tensor components of the chemical shifts for trimethylsilylacetylene, methoxyacetylene, and propiolaldehyde have been measured and are in reasonable agreement with the calculations. The downfield shift at the terminal carbon of propiolaldehyde along with a small upfield shift at the adjacent carbon has been found to result from the coupling of the in-plane pi MO of the acetylene with the pi* orbital that has a node near the central carbon. The tensor components for acetonitrile also have been measured, and the shielding of cyano and acetylenic carbons are compared.