CNDO calculation of the diamagnetic susceptibilities of cyclic conjugated molecules

We present the results obtained within the CNDO /2 method for the diamagnetic susceptibilities of a series of conjugated molecules including benzene, fulvene, pyridine, pyrrole, and furan. The calculations have been carried out, according to Bley's method, in the framework of coupled Hartree–Fock perturbation theory, using gauge invariant atomic orbitals and the London approximation. The two usual CNDO parametrizations, due respectively to Pople, Santry, and Segal and to Del Bene and Jaffe, did not enable one to obtain the principal susceptibilities of the molecules under consideration. We investigated the effect of the CNDO parametrizations on the different contributions to the susceptibility, and gave a first improvement to Bley's method, which permits the calculation of the mean susceptibility of conjugated molecules with less than 9% error.     

Instability of stable negative ions in the Xα method or other local density functional schemes

Stable negative ions such as O^? and F^? are calculated to be unstable if exchange—correlation potentials of the Xα type are used, whereas the Hartree—Fock method correctly predicts their stability. The reason for this difference in stability is the difference in the one-electron energies between the two schemes.     

Generalization of thermodynamics in of fractional-order derivatives and calculation of heat-transfer properties of noble gases

Journal of Thermal Analysis and Calorimetry - This paper performs an entropy generation analysis for nanofluid turbulent flow around a rotating cylinder. A finite volume method based on SST...     

A density functional theory calculation of the electronic properties of several high-spin and low-spin iron(II) pyrazolylborate complexes

Density functional theory has been used to study the electronic spin-state properties of low-spin Fe[HB(pz)3]2, high-spin Fe[HB(3-Mepz)3]2, high-spin Fe[HB(3,5-Me 2pz)3]2, and high-spin Fe[HB(3,4,5-Me 3pz)3]2 complexes that exhibit very different iron(II) electronic spin-sate crossover behaviors with changing temperature and pressure. Excellent agreement is obtained between the experimentally observed M?ssbauer-effect quadrupole splittings and isomer shifts of these complexes and those calculated with the B3LYP functional and various different basis sets for both the high-spin and low-spin states of iron(II). The calculations for Fe[HB(pz)3]2 that use the LANL2DZ, 6-31++G(d,p), and 6-311++G(d,p) basis sets for iron all lead to very similar electric field gradients and thus quadrupole splittings. The initial calculations, which were based upon the known X-ray structures, were followed by structural optimization, an optimization that led to small increases in the Fe-N bond distances. Optimization led to at most trivial changes in the intraligand bond distances and angles. The importance of the 3-methyl-H...H-3-methyl nonbonded intramolecular interligand interactions in controlling the minimum Fe-N bond distances and determining the iron(II) spin state both in Fe[HB(3-Mepz)3]2 and in the related methyl-substituted complexes has been identified.     

An investigation into the use of expansion methods in the calculation of chemical shifts and diamagnetic susceptibilities

The chemical shift and the diamagnetic susceptibility of the hydrogenic atom with magnetic dipole and origin of the external magnetic field vector potential noncoincident with the hydrogenic nucleus have been calculated from perturbation theory using a set of expansion functions whose radial parts are single exponent associated Laguerre functions. In contrast to hydrogenic expansion functions these functions give rapid convergence to the exact values of the second-order energy summations when centred at the hydrogenic nucleus. The rate of convergence is fairly insensitive to the choice of expansion function exponent.     

Calculating the thermodynamic and structural properties of liquid noble gases

The structural and energy parameters in the liquid branch of the saturation line of noble gases, including ⁴He, are determined.     

Calculation of diamagnetic susceptibilities from the Thomas–Fermi–Amaldi theory for ions

The importance of the Fermi–Amaldi correction for the Thomas–Fermi theory for ions is analyzed. An analytical trial electron density with a proper limiting behavior has been chosen and the Thomas–Fermi–Amaldi energy expression is minimized in a variational way for several closed-shell ions. The obtained optimum electron densities have been used for calculating the diamagnetic susceptibilities, and a comparison with other theoretical results and available experimental data is performed.     

Some applications of local density functional theory to the calculation of reaction energetics

Summary We report the results of a local density functional investigation of the energetics of some isomerization reactions, involving the conversions of several unsaturated systems to highly strained molecules related to triprismane and tetrahedrane. The program DMol was used at the DNP level to compute the activation barriers and total energy changes associated with these processes. We also show, for more than 70 first- and second-row atoms and molecules, that the errors (non-local corrections) in their energies correlate very well with the number of electrons, within isonuclear series. This should provide a useful empirical means for improving dissociation energies obtained within the local approximation.     

Suitability of density functional methods for calculation of electrostatic properties

A systematic analysis was performed on the suitability of the molecular electrostatic potential (MEP) and MEP-derived properties determined by means of density functional (DFT) methods. Attention was paid to the electrostatic potential (ESP) derived charges, the ESP and exact quantum mechanical dipole moments, the depth of MEP minima, and the MEP distribution in layers around the molecule for a large series of molecules. The electrostatic properties were determined at either local or nonlocal DFT levels using different functionals. The results were compared with the values estimated from quantum mechanical calculations performed at Hartree–Fock, Møller–Plesset up to fourth order, and CIPSI levels. The suitability of the MEP-derived properties estimated from DFT methods is discussed for application in different areas of chemical interest. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 980–991, 1997     

Density functional calculation of the electronic absorption spectrum of Cu+ and Ag+ aqua ions

The UV absorption of aqueous Cu+ and Ag+ has been studied using Time Dependent Density Functional Theory (TDDFT) response techniques. The TDDFT electronic spectrum was computed from finite temperature dynamical trajectories in solution generated using the Density Functional Theory (DFT) based Ab Initio Molecular Dynamics (AIMD) method. The absorption of the two ions is shown to arise from similar excitation mechanisms, namely transitions from d orbitals localized on the metal center to a rather delocalized state originating from hybridization of the metal s orbital to the conduction band edge of the solvent. The ions differ in the way the spectral profile builds up as a consequence of solvent thermal motion. The Cu+ absorption is widely modulated, both in transition energies and intensities by fluctuations in the coordination environment which is characterized by the formation of strong coordination bonds to two water molecules in an approximately linear geometry. Though, on average, absorption intensities are typical of symmetry forbidden transitions of metal ions in the solid state, occasionally very short (<100 fs) bursts in intensity are observed, associated with anomalous Cu-H interactions. Absorption by the Ag+ complex is in comparison relatively stable in time, and can be interpreted in terms of the energy splitting of the metal 4d manifold in an average crystal field corresponding to a fourfold coordination in a distorted tetrahedral arrangement. Whereas the spectral profile of the Ag+ aqua ion is in good agreement with experiment, the overall position of the band is underestimated by 2 eV in the BLYP approximation to DFT. The discrepancy with experiment is reduced to 1.3 eV when a hybrid functional (PBE0) is used. The remaining inaccuracy of TDDFT in this situation is related to the delocalized character of the target state in d-->s transitions.     

A density functional study of the electronic spectrum of permanganate

Multiplet splittings for six excited electronic configurations of the permanganate ion, MnO₄⁻, are calculated. Earlier density functional calculations on the same subject are improved upon by the numerical evaluation of some two-electron integrals to resolve certain multideterminantal states. Excellent agreement with the experimental spectrum is obtained, and a reassignment of bands in the 25,000–35,000 cm⁻¹ range is proposed. Fully symmetric (a₁) vibrational frequencies are calculated, and the origin and magnitude of the most significant Jahn-Teller distortions of the excited states are discussed. © 1996 John Wiley & Sons, Inc.     

A local density functional—LCAO method with effective potentials for obtaining the electronic structure of molecules and clusters

The introduction of effective core potentials into the Xα LCAO method can greatly expand the size of the molecular system which may be feasibly treated. Therefore, we have incorporated norm-conserving effective potentials into this framework, using a gaussian basis set, for calculations of the electronic structure of molecules and cluster models of solids and surfaces. This combined methodology makes possible the determination of equilibrium geometries and certain aspects of potential surfaces required for a large class of electronic structure problems.     

A density matrix-based method for the linear-scaling calculation of dynamic second- and third-order properties at the Hartree-Fock and Kohn-Sham density functional theory levels

A density matrix-based time-dependent self-consistent field (D-TDSCF) method for the calculation of dynamic polarizabilities and first hyperpolarizabilities using the Hartree-Fock and Kohn-Sham density functional theory approaches is presented. The D-TDSCF method allows us to reduce the asymptotic scaling behavior of the computational effort from cubic to linear for systems with a nonvanishing band gap. The linear scaling is achieved by combining a density matrix-based reformulation of the TDSCF equations with linear-scaling schemes for the formation of Fock- or Kohn-Sham-type matrices. In our reformulation only potentially linear-scaling matrices enter the formulation and efficient sparse algebra routines can be employed. Furthermore, the corresponding formulas for the first hyperpolarizabilities are given in terms of zeroth- and first-order one-particle reduced density matrices according to Wigner's (2n+1) rule. The scaling behavior of our method is illustrated for first exemplary calculations with systems of up to 1011 atoms and 8899 basis functions.     

Growth pattern and electronic properties of acetonitrile clusters: a density functional study

We report a systematic theoretical study on the growth pattern and electronic properties of acetonitrile clusters [(CH(3)CN)(n) (n = 1, 9, 12)] using density functional approach at the B3LYP6-31++G(d,p) level. Although we have considered a large number of configurations for each cluster, the stability of the lowest energy isomer was verified from the Hessian calculation. It is found that the lowest energy isomer of the dimer adopts an antiparallel configuration. For trimer and tetramer, cyclic ring structures were found to be favored over the dipole stabilized structure. In general, it is found that the intermolecular CH...N interactions play a significant role in the stabilization of the cyclic layered geometry of acetonitrile clusters. A critical comparison between trimer and tetramer clusters suggests that the three member cyclic ring is more stable than four member rings. The growth motif for larger clusters (n = 5-9, 12) follows a layered pattern consisting of three or four membered rings, which, in fact, is used as the building block. Based on the stability analysis, it is found that clusters with an even number of molecular entities are more stable than the odd clusters, except trimer and nonamer. The exceptional stability of these two clusters is attributed to the formation of trimembered cyclic rings, which have been found to form the building blocks for larger clusters.     

Calculation of the electronic and optical properties of indium tin oxide by density functional theory

Density functional theory (DFT) was used to calculate the bulk electronic and optical properties of indium tin oxide (ITO). The ITO model was constructed replacing indium atoms with tin atoms in the cubic unit cell of indium oxide. To allow more possibilities for tin atom substitution than afforded by the forty-atom primitive cell of indium oxide all eighty atoms of the unit cell were included in the stoichiometry (In_32−xSn_xO₄₈) using periodic boundary conditions. A number of properties of ITO were calculated including the optical band gap, charge carrier density and plasma frequency. The dependence of the electronic and optical properties of ITO on a variety of parameters such as the tin content, cubic lattice parameter and the distance between adjacent tin atoms was investigated. The electronic and optical properties agreed well with experimental data and allowed insight into the origin of the electronic and optical properties of ITO.