Basis set dependence of the molecular electrostatic potential topography. A case study of substituted benzenes

The basis set dependence of the topographical structure of the molecular electrostatic potential (MESP), as well as the effect of substituents on the MESP distribution, has been investigated with substituted benzenes as test cases. The molecules are studied at HF-SCF 3�21G and 6�31G** levels, with a further MESP topographical investigation at the 3�21G, double-zeta, 6�31G*, 6�31G**, double-zeta polarized and triple-zeta polarized levels. The MESP critical points for a 3�21G optimized/6�31G** basis are similar to the corresponding 6�31G** optimized/6�31G** ones. More generally, the qualitative features of the MESP topography computed at the polarized level are independent of the level at which optimization is carried out. For a proper representation of oxygen lone pairs, however, optimization using a polarized basis set is required. The nature of the substituent drastically changes the MESP distribution over the phenyl ring. The values and positions of MESP minima indicate the most active site for electrophilic attack. This point is strengthened by a study of disubstituted benzenes.

Intramolecular hydrogen bond: Can it be part of the basis set of valence internal coordinates in normal mode analysis?

It has been shown earlier1] that the relaxed force constants (RFCs) could be used as a measure of bond strength only when the bonds form a part of the complete valence internal coordinates (VIC) basis. However, if the bond is not a part of the complete VIC basis, its RFC is not necessarily a measure of bond strength. Sometimes, it is possible to have a complete VIC basis that does not contain the intramolecular hydrogen bond (IMHB) as part of the basis. This means the RFC of IMHB is not necessarily a measure of bond strength. However, we know that IMHB is a weak bond and hence its RFC has to be a measure of bond strength. We resolve this problem of IMHB not being part of the complete basis by postulating `equivalent' basis sets where IMHB is part of the basis at least in one of the equivalent sets of VIC. As long as a given IMHB appears in one of the equivalent complete VIC basis sets, its RFC could be used as a measure of bond strength parameter.

Ground and Excited State Intramolecular Proton Transfer in Salicylic Acid: an Ab Initio Electronic Structure Investigation

Energetics of the ground and excited state intramolecular proton transfer in salicylic acid have been studied by ab initio molecular orbital calculations using the 6-31G** basis set at the restricted Hartree-Fock (RHF) and configuration interaction-single excitation (CIS) levels and also using the semiempirical method AM1 at the RHF level as well as with single and pair doubles excitation configuration interaction spanning eight frontier orbitals (PECI = 8). The ab initio potential energy profile for intramolecular proton transfer in the ground state reveals a single minimum corresponding to the primary form, in the first excited singlet state, however, there are two minima corresponding to the primary and tautomeric forms, separated by a barrier of similar to 6 kcal/mol, thus accounting for dual emission in salicylic acid. Electron density changes with electronic excitation and tautomerism indicate no zwitterion formation. Changes in spectral characteristics with change in pH, due to protonation and deprotonation of salicylic acid, are also accounted for, qualitatively. Although the AM1 calculations suggest a substantial barrier for proton transfer in the ground as well as the first excited state of SA, it predicts the transition wavelength in near quantitative accord with the experimental results for salicylic acid and its protonated and deprotonated forms.

Theoretical studies of lithium bonding in lithium chloride/aliphatic amine complexes

In the systematic study of amine … LiCl [amines = NH3, CH3NH2, (CH3)2NH] complexes the possibility of an ion-pair structure and the effect of methylation on the stabilization energy is investigated. ΔEis evaluated by the SCF/4-31G method and augmented by the approximate dispersion energy calculated perturbationally. The interaction energy decreases with the increasing number of methyl groups in the amine. The dispersion energy plays a negligible role in the stabilization of complexes. None of the systems studied are ion pairs; their Li bonds are of a so-called molecular type. Due to the divergence of the multipole expansion, the attempt to correct the 4-31G stabilization energies via the electrostatic energy fails. The relative order of the ΔE in the series of complexes is verified instead in the extended basis set calculation. The lithium bonds are compared with their H-bonded analogues.

Investigation of inter-ion interactions in N,N,N',N'-tetramethylethylenediammonium dithiocyanate via experimental and theoretical charge density studies

The crystal structure of the N,N,N',N'-tetramethylethylenediammonium dithiocyanate salt has been examined by experimental charge density studies from high-resolution X-ray diffraction data. The corresponding results are compared with multipole refinements, using theoretical structure factors obtained from a periodic density functional theory calculation at the B3LYP level with a 6-31G** basis set. The salt crystallizes in space group P (1) over bar and contains only a single ion pair with an inversion center in the cation. The salt has thus one unique classical N+-H center dot center dot center dot(NCS)(-) hydrogen bond but also has six other weaker interactions: four C-H center dot center dot center dot S, one C-H center dot center dot center dot N, and one C-H center dot center dot center dot C-pi. The nature of all these interactions has been examined topologically using Bader's quantum theory of "atoms in molecules" and all eight of the Koch-Popelier criteria. The experimental and theoretical approaches agree well and both show that the inter-ion interactions, even in this simplest of systems, play an integrated and complex role in the packing of the ions in the crystal. Electrostatic potential maps are derived from experimental charge densities. This is the first time such a system has been examined in detail by these methods.

Enhanced dielectric response of ZrO2 upon Ti doping and introduction of O vacancies

We determine the electronic properties and dielectric response of zirconia (ZrO2) with oxygen vacancies (O vacancies) and Ti doping using first-principles density functional theory calculations based on pseudopotentials and a plane wave basis. We find significantly enhanced static dielectric response in zirconia with Ti doping and introduction of oxygen vacancies. Softening of phonon modes are responsible for the enhanced dielectric response of doped samples compared to pure zirconia.

InMnO3: A biferroic

InMnO3 which has a hexagonal structure similar to that of YMnO3 is found to show a canted antiferromagnetic behavior below 50 K (T-N) and a ferroelectric (FE) transition at 500 K accompanied by hystersis. We have determined the structure, polarization, and energetics of the FE and paraelectric (PE) phases of InMnO3 using first-principles density functional theory calculations based on pseudopotentials and a plane-wave basis, and find the polarization of the PE phase to be a half-integer quantum. The difference in polarization of the FE and PE phases calculated along a simple path is different from the absolute value of polarization of the FE phase. A weak piezoelectric response is exhibited by InMnO3 to uniaxial strain.

Topological analysis off charge density distribution in concomitant polymorphs of 3-acetylcoumarin, a case off packing polymorphism

Detailed investigation of the charge density distribution in concomitant polymorphs of 3-acetylcoumarin in terms of experimental and theoretical densities shows significant differences in the intermolecular features when analyzed based on the topological properties via the quantum theory of atoms in molecules. The two forms, triclinic and monoclinic (Form A and Form B), pack in the crystal lattice via weak C-H---O and C-H---pi interactions. Form A results in a head-to-head molecular stack, while Form B generates a head-to-tail stack. Form A crystallizes in PI (Z' = 2) and Form B crystallizes in P2(1)/n (Z = 1). The electron density maps of the polymorphs demonstrate the differences in the nature of the charge density distribution in general. The charges derived from experimental and theoretical analysis show significant differences with respect to the polymorphic forms. The molecular dipole moments differ significantly for the two forms. The lattice energies evaluated at the HF and DFT (B3LYP) methods with 6-31G** basis set for the two forms clearly suggest that Form A is the thermodynamically stable form as compared to Form B. Mapping of electrostatic potential over the molecular surface shows dominant variations in the electronegative region, which bring out the differences between the two forms.

Bond length distortions associated with torsion potentials; ab-initio studies on methanedio

The 4-31G basis set is used to study the bond length variations as functions of dihedral angels in methanediol. This study is compared with O---C---O bond angle optimization studies by Gorenstein and Kar and the possible reason for bond length shorteing in the trans---trans configuration is analysed.

Effects of Zr and Ti doping on the dielectric response of CeO2: A comparative first-principles study

Zr doping in ceria (CeO2) results in enhanced static dielectric response compared to pure ceria. On the other hand, Ti doping in ceria keeps its dielectric constant unchanged. We use first-principles density functional theory calculations based on pseudopotentials and a plane wave basis to determine electronic properties and dielectric response of Zr/Ti-doped and oxygen-vacancy-introduced ceria. Softening of phonon modes is responsible for the enhancement in dielectric response of Zr-doped ceria compared to that of pure ceria. The ceria-zirconia mixed oxides should have potential use as high-k materials in the semiconductor industry. (c) 2010 Elsevier Ltd. All rights reserved.

An examination of the BEBO model with the results of ab initio calculations of a reaction series

The BEBO (bond energy-bond order) model of Johnston and Parr is examined with the results of ab initio MO calculations on a series of metathetic reactions which involve hydrogen transfer. Energies are calculated at the 6-31G**/PMP2 = full//6-31G** level while the bond orders are estimated using the 6-31G** basis set with the geometry optimisation at the single configuration unrestricted Hartree-Fock frame. Our analysis reveals that the bond-order exponent in the BEBO theory is greater than unity for the reaction series and the entropy term becomes implicitly present in the BEBO model.

Ab initio studies of the structure and vibrational spectra of methyl isocyanide-borane and acetonitrile-borane

The SCF/DZP and MP2/DZP methods of ab initio quantum chemistry have been utilized to study the structure, vibrational spectra, binding energy, and barrier to internal rotation of methyl isocyanide-borane and acetonitrile-borane adducts. The eclipsed conformation of the complexes was predicted to be a minimum, and the staggered form is a transition state with a barrier height of about 10 cal/mol. The vibrational analyses of CH3NC-BH3 and CH3CN-BH3 and several of their isotopomers have been carried out by the GF matrix method. Computations have also been carried out for free CH3NC and CH3CN in order to investigate the changes in CH3NC and CH3CN as a result of their complex formation with BH3. To obtain an acceptable set of force constants, a recently proposed procedure ''RECOVES'' has been utilized. The increase in the N=C/C=N stretching force constant of CH3NC/CH3CN on adduct formation is interpreted with the help of Parr and Borkman's model. The binding energies for the two adducts have been determined taking basis set superposition error (BSSE) into consideration. The effect of the BSSE on structure, dipole moment, and vibrational frequencies of CH3CN and CH3NC is also evaluated. The predicted infrared band intensities for the two complexes are in good agreement with the experimentally observed features, and they have been utilized in the assignment of vibrational frequencies.

Electron donor-acceptor complexes of I-2 with diethyl ether and diethyl sulphide. An ab initio MO study

The charge-transfer complexes of I-2 with the n-donors diethyl ether and diethyl sulfide were studied at the Hartree-Fock and MP2 levels. The structures were fully optimized using the 3-21G((*)) basis set as well as with effective core potentials. The calculations consistently yield a C-2v structure for the ether-I-2 complex, but an unsymmetrical form for the sulfide-I-2 complex. A natural bond orbital analysis and the BSSE-corrected complexation energies reveal stronger interactions in the sulfide complex. The computed orbital energies of the monomers and complexes reproduce the trends in experimentally observed vertical ionization potentials.

SCF and electron correlation studies on structure, force constants and vibrational spectra of borane monoammoniate complex

Effects of basis set and electron correlation on the equilibrium geometry, force constants and vibrational spectra of BH3NH3 have been studied. A series of basis sets ranging from double zeta to triple zeta including polarization and diffuse functions have been utilized. All the SCF based calculations overestimate the dative B-N bond distance and considerable improvement occurs when the treatment for electron correlation is introduced. Detailed vibrational analysis for BH3NH3 has been carried out. The mean absolute percentage deviation of the ab initio predicted vibration frequencies of (BH3NH3)-B-11 from the experiment is about 10% for the SCF based calculations and the MP2 method shows better agreement, the overall deviation being 5-6%. The ground state effective force constants of BH3NH3 were obtained using RECOVES procedure. The RECOVES sets of force constants are found to be highly satisfactory for the prediction of the vibrational frequencies of different isotopomers of BH3NH3. The mean absolute percentage deviation of the calculated frequencies of different isotopomers from the experiment is much less than 1%. The RECOVES-MP2/augDZP set of force constants was found to be the best set among the different sets for this molecule. Theoretical infrared intensities are in fair agreement with the observed spectral features.

Isomeric identification of methylated naphthalenes using gas phase infrared spectroscopy

We report gas phase mid-infrared spectra of 1- and 2- methyl naphthalenes at 0.2 cm(-1) resolution. Assignment of observed bands have been made using scaled quantum mechanical (SQM) calculations where the force fields rather the frequencies are scaled to find a close fit between observed and calculated bands. The structure of the molecules has been optimized using B3LYP level of theory in conjunction with standard 6-311G** basis set to obtain the harmonic frequencies. Using the force constants in Cartesian coordinates from the Gaussian output, scaled force field calculations are carried out using a modified version of the UMAT program in the QCPE package. Potential energy distributions of the normal modes obtained from such calculations helped us assign the observed bands and identify the unique features of the spectra of 1- and 2-MNs which are important for their isomeric identification.