Structure and conformational stability of CH2CHCH2X (X=F, Cl and Br) molecules—post Hartree–Fock and density functional theory methods

The molecular structure and conformational stability of CH₂CHCH₂X (X=F, Cl and Br) molecules were studied using ab initio and density functional theory (DFT) methods. The molecular geometries of 3-fluoropropene were optimized employing BLYP and B3LYP levels of theory of DFT method implementing 6-311+G(d,p) basis set. The MP2/6-31G^*, BLYP and B3LYP levels of theory of ab initio and DFT methods were used to optimize the 3-chloropropene and 3-bromopropene molecules. The structural and physical parameters of the molecules are discussed with the available experimental values. The rotational potential energy surface of the above molecules were obtained at MP2/6-31G^* and B3LYP/6-311+G(d,p) levels of theory. The Fourier decomposition of the rotational potentials were analyzed. The HF/6-31G^* and MP2/6-31G^* levels of theory have predicted the cis conformer as the minimum energy structure for 3-fluoropropene, which is in agreement with the experimental values, whereas the BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory reverses the order of conformation. The ΔE values calculated for 3-chloropropene at MP2/6-31G^*, BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory show that the gauche form is more stable than the cis form, which is in agreement with the experimental value. The same levels of theory have also predicted that the gauche form is stable than cis for 3-bromopropene molecule. The maximum hardness principle has been able to predict the stable conformer of 3-fluoropropene at HF/6-31G^* level of theory, but the same level of theory reverses the conformational stability of 3-chloropropene and 3-bromopropene molecules and MP2/6-31G^* level of theory predicted the stable conformer correctly.     

Molecular structure and vibrational spectra of indole and 5-aminoindole by density functional theory and ab initio Hartree–Fock calculations

The molecular geometry and vibrational frequencies of indole and 5-aminoindole in the ground state have been calculated by using the Hartree–Fock and density functional method with 6-311++G(d,p) basis set. The optimized geometrical parameters obtained by DFT calculations are in good agreement with the experimental values. Comparison of the observed fundamental vibrational frequencies of indole and 5-aminoindole with the calculated results by density functional and Hartree–Fock methods indicates that B3LYP is superior to the scaled Hartree–Fock approach for molecular vibrational problems. The theoretical spectrograms for FT-IR spectrum of 5-aminoindole have been constructed.     

Quantum chemical studies on structure, conformation and isomerization of nitroso, nitro substituted benzene and 1,3-cyclopentadiene

The molecular structure, conformational stability and isomerization of nitroso, nitro substituted benzene and 1,3-cyclopentadiene in gas phase have been investigated using ab initio and density functional theory methods. The molecular geometries and energetics of possible conformers were obtained by employing MP2, B3LYP and B3PW91 levels of theory implementing 6-31G* basis set. The relative stabilities of the conformations were evaluated from the energy differences of the structure. Chemical hardness (η) and chemical potential (μ) were calculated at HF/6-31G* level of theory for all the positional and geometrical isomers to study the maximum hardness principle. Each optimized structure has been tested against the imaginary frequencies at MP2/6-31G* level of theory in order to be sure they are located at energy minimum.     

Structure and stability of XeF6 isomers: density functional theory study and the maximum hardness principle

The geometrical configurations of the four possible isomers with C_3v, O_h, C_s and C_2v symmetry on the potential energy surface of the XeF₆ molecule are optimized by using DFT-LDA/NL. Their relative energies, vibration frequencies, electronic chemical potential and hardness have been calculated. It is found that the C_2v configuration has one imaginary frequency. The relative energies of the four isomers increase in order of C_3v, O_h, C_s and C_2v, and the hardness values in same order. The isomer stability obeys the maximum hardness principle (MHP), while their hardness values are very close to each other. It is quite evident that the very close hardness is the main reason for the structure fluxionality of XeF₆.     

Molecular structure and vibrational spectra of lepidine and 2-chlorolepidine by density functional theory and ab initio Hartree–Fock Calculations

The molecular geometry and vibrational frequencies of lepidine and 2-chlorolepidine in the ground state have been calculated by using the Hartree–Fock and density functional methods (B3LYP) with 6-31G (d) as the basis set. The optimized geometric bond lengths obtained by using B3LYP and bond angles obtained by HF that correspond to the experimental values of 2-cl-lepidine molecule were given. Comparison of the observed fundamental vibrational frequencies of lepidine and 2-chlorolepidine, and calculated results by density functional B3LYP and Hartree–Fock methods indicates that B3LYP is superior to the scaled Hartree–Fock approach for molecular vibrational problems.     

The isomers of phenalene and their singlet and triplet states: A Hartree–Fock and density functional computational investigation

Each of the isomers of phenalene, 1H-, 2H-, 3aH-, and 9bH-phenalene, as well as the cation, neutral radical, and anion in the phenalenyl system, have been examined at the Hartree–Fock 6-31G(d) and the density functional B3LYP/6-31G(d) levels of theory. The structures and properties of the phenalenes were determined as both singlets and triplets. While the data indicate that both the 9bH- and the 1H-isomer will exist as ground state singlets, the 2H- and 3aH-phenalenes are predicted to exist as ground state triplets; only the synthesis of the 1H-isomer has been reported in the literature. Structurally, 1H- and 2H-phenalene are planar systems, 3aH- and 9bH-phenalene are non-planar systems puckered at the saturated carbon, and the cation, neutral radical, and anion of phenalenyl are planar D_3h systems.     

A density functional study of the glycine molecule: Comparison with post-Hartree–Fock calculations and experiment

The potential energy surface of un-ionized glycine has been explored with density functional theory. The performance of several nonlocal functionals has been evaluated and the results are presented in the context of available experimental information and post-Hartree–Fock quantum chemical results. The zero-point and thermal vibrational energies along with vibrational entropies play a very important role in determining the relative stability of glycine conformers; the realization of this has led to some revision and reinterpretation of the experimental results. Uncertainties in the vibrational contributions to the energy differences of several tenths of a kilocalorie/mole remain. The uncertainty in the vibrational free energy is even larger, about 1 kcal/mol. In the final analysis, we suggest that the best estimate of the electronic energy difference between the two lowest glycine conformers should be revised downward from 1.4 to 1.0 kcal/mol. Thirteen stationary points on the potential energy surface have been localized. For the majority of these, there is close agreement among various nonlocal density functionals and the post-Hartree–Fock methods. However, the second conformer (IIn), which has a strong hydrogen bond between the hydroxyl hydrogen and the nitrogen of the amine group, presents a distinct challenge. The relative energy of this conformer is extremely sensitive to the basis set, the level of correlation, or the functional used. The widely used BP86, PP86, and BP91 nonlocal functionals overestimate the strength of the hydrogen bond and predict that this conformer is the lowest energy structure. This contradicts both experiment and high-level post-Hartree–Fock studies. The adiabatic connection method (ACM) and the BLYP functional yield the correct order. The ACM method, in particular, gives energies which are in reasonable agreement with MP2, although these are somewhat low as compared with experiment. Based on this study, ACM should perform well for this type of bioorganic application, with typical errors of a few tenths of a kilocalorie/mole and only rarely exceeding 0.5 kcal/mol. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1609–1631, 1997     

Molecular structure, vibrational spectra and NBO analysis of phenylisothiocyanate by density functional method

The molecular geometry, vibrational frequencies and NBO analysis of phenylisothiocyanate (PITC) in the ground state have been calculated by using density functional theory calculation (B3LYP) with 6-311++G(d,p) basis set. The optimized geometrical parameters obtained by DFT calculations are in good agreement with experimental values. Comparison of the observed fundamental vibrational frequencies of the PITC and calculated result by density functional theory (B3LYP) indicates B3LYP is superior for molecular vibrational problems. The entropy of the title compound was also performed at HF/B3LYP/6-311++G(d,p) levels of theory. Natural bond orbital (NBO) analysis of title molecule is also carried out. A detailed interpretation of the IR and Raman spectra of PITC is reported on the basis of the calculated potential energy distribution (PED). The theoretical spectrogram for IR spectrum of the title molecule has been constructed.     

Density functional theory investigations on the conformational stability and vibrational spectra of the model compound of vitamin K

Becke 3-Lee–Yang–Parr density functional theory (DFT/B3LYP) using 6-31G(d) and 6-311G(d) basis sets and the scaled quantum mechanics (SQM) force field method are used to study molecular conformations and vibrational spectra of a model compound of vitamin K (VK). In this molecule, there are six conformers on the torsional potential energy map of the dihedral angles C₈C₁₄C₁₅C₁₆ () and C₇C₈C₁₄C₁₅ (β). It is shown that the VK_1 conformer ( = 237.7° and β = 274.2°) is the most stable form. For this lowest energy conformer, the harmonic force fields calculated by B3LYP/6-311G(d) and B3LYP/6-31G(d) methods are scaled with one scale factor of 0.9623 and a set of different scale factors transferred from the previous studies for the other similar molecules, respectively. The vibrational frequencies, IR intensities, and Raman activities are obtained for the lowest energy conformer. On the basis of B3LYP calculations, the normal mode analysis is performed to assign the vibrational fundamental frequencies according to the potential energy distributions. The computational frequencies are in good agreement with the observed results.     

A comparative post-Hartree–Fock and density functional theory study of monochalcogenide diatomic molecules

HOMO and LUMO (FMOs) play important role in the optical properties of meridianal isomer of tris(8-hydroxyquinolino)aluminum (mer-Alq3) and its derivatives. The frontier molecular orbitals (FMOs) also play a vital role in the process of charge transport. It is urgent to find the reason of FMO distribution pattern among the ligands. The structures of mer-Alq3 and its “CH”/N substituted derivatives have been optimized at the B3LYP/6-31G^* level. Energy decomposition analysis has been performed at the B3LYP/DZP level. The results of energy-partitioning analysis of ground states are discussed. It has been explained that HOMOs are on A-ligands due to weaker electrostatic interaction energy between L_a-AlL_bL_c fragments while LUMOs are on B-ligands due to weaker orbital interaction energy between L_b-AlL_aL_c fragments.     

Post Hartree–Fock and density functional theory studies on Di-Protonated Allopurinol

Post Hartree–Fock and density functional theory (DFT) methods have been employed to study the molecular properties of Di-Protonated Allopurinol²⁺ tautomers in gaseous and aqueous phase environments. The tautomers in gaseous phase have been optimized at MP2/6-311G(2d,2p) and B3LYP/6-311G(2d,2p) levels of theory. The self-consistent reaction field theory (SCRF) has been employed to optimize the tautomers in aqueous phase (ε = 78.5) at B3LYP/6-311G(2d,2p) level of theory and the solvent effect has been studied. The structure, energetics and relative stabilities of the tautomers have been analyzed both in gaseous and aqueous phases. The principle of maximum hardness (MHP) has been tested at B3LYP/6-311G(2d,2p) level of theory. The condensed Fukui functions have been calculated using the atomic charges obtained through Natural population analysis to identify the relative change in the most reactive site of the optimized structures. NMR studies have been carried out, on the basis of Cheeseman coworker’s method, to analyze the molecular environment as well as the delocalization activities of electron clouds.     

Analysis of vibrational spectra of 2 and 3-methylpiperidine based on density functional theory calculations

The experimental and theoretical vibrational spectra of 2 and 3-methylpiperidine (abbreviated as 2-MP and 3-MP) were studied. The FT-Infrared spectra of 2-MP and 3-MP molecules were recorded in the liquid phase. The structural and spectroscopic analysis of the title molecules were made by using density functional harmonic calculations. For the title molecules, only one form was found most stable structure by using B3LYP level with the 6-311G (d,p) basis set. Selected experimental bands were assigned and characterized based on the scaled theoretical wave numbers by their total energy distribution (TED).     

Structure and conformational equilibrium of thiacalix[4]arene by density functional theory

Density functional theory calculations for the structure and conformational equilibrium of thiacalix[4]arene are reported. The conformational equilibrium of thiacalix[4]arene, a heterocalixarene in which the phenol groups are bridged by sulphur atoms is compared to the conformational equilibrium of calix[4]arene. Thiacalix[4]arene conformational energies relative to the cone conformer (ΔE's) are reduced in comparison with calix[4]arene. This conformational change is in qualitative agreement with recent NMR spectroscopy measurements of the conformational equilibrium for a tetraethylether of thiacalix[4]arene in a CDCl₃ solution which indicates an enhanced chemical exchange of thiacalixarene conformers in comparison with similar methylene bridged structures. Density functional theory results for the structure of thiacalix[4]arene are in good agreement with recent X-ray diffraction measurements. The electrostatic potentials in the cone conformers of thiacalix[4]arene and calix[4]arene suggest that their complexation or recognition abilities can be significantly different. Dipole moments of the four thiacalix[4]arene conformers are in the order: cone>1,2-alternate>partial-cone>1,3-alternate.     

Ab initio and density functional studies on the structure and vibrational spectra of 2-hydroxy-1,4-naphthoquinone-1-oxime derivatives

Structure and vibrational frequencies of lawsoneoxime and its C₃-substituted (R=CH₃, NH₂, Cl, NO₂) derivatives in keto and nitrosophenol forms have been obtained employing the Hartree–Fock and density functional methods. Charge distributions in different conformers have been studied using the molecular electrostatic potential topography as a tool. For all these derivatives except for nitrolawsoneoxime the amphi conformer in the keto form is predicted to be of lowest energy, which can partly be attributed to hydrogen bonding through the oximino nitrogen. In the nitro derivative, however, the preference to form a six membered ring owing to O–H–O hydrogen-bonded interactions makes the anti conformer (keto) the stablest. Further one of the nitrosophenol conformers of nitrolawsoneoxime turns out to be very close in energy (0.21 kJ mol^–1 higher) to this anti conformer. The consequences of hydrogen bonding on charge distribution and vibrational spectra are discussed.     

Comparative studies on the structures,infrared spectrum,and thermodynamic properties of phthalocyanine using ab initio Hartree–Fock and density functional theory methods

The molecular structure, vibrational spectrum, standard thermodynamic functions, and enthalpy of formation of free base phthalocyanine (Pc) have been studied using the density functional theory B3LYP procedure, as well as the ab initio Hartree–Fock method. Various basis sets 3‐21G, 6‐31G*, and LANL2DZ have been employed. The results obtained at various levels are discussed and compared with each other and with the available experimental data. It is shown that calculations performed at the Hartree–Fock level cannot produce a reliable geometry and related properties such as the dipole moment of Pc and similar porphyrin‐based systems. Electron correlation must be included in the calculations. The basis set has comparatively less effect on the calculated results. The results derived at the B3LYP level using the smaller 3‐21G and LANL2DZ basis sets are very close to those produced using the medium 6‐31G* basis set. The geometry of Pc obtained at the B3LYP level has D_2h symmetry and the diameter of the central macrocycle is about 4 Å. The enthalpy of formation of Pc in the gas phase has been predicted to be 1518.50 kJ/mol at the B3LYP/6‐311G(2d,2p)//B3LYP/6‐31G* level via an isodesmic reaction. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

The exact Fermi potential yielding the Hartree–Fock electron density from orbital‐free density functional theory

The exact expression for the Fermi potential yielding the Hartree–Fock electron density within an orbital‐free density functional formalism is derived. The Fermi potential, which is defined as that part of the potential that depends on the particles’ nature, is in this context given as the sum of the Pauli potential and the exchange potential. The exact exchange potential for an orbital‐free density functional formalism is shown to be the Slater potential.     

Gas-phase properties and Fukui indices of sulfine (CH2SO). Potential energy surface and maximum hardness principle for its protonation process. A density functional study

Gas-phase thermochemical properties of sulfine (CH₂SO) and the potential energy surface of its protonation process were studied by the density functional method employing different exchange-correlation potentials. All calculations showed that the most stable protonated isomer is planar with the proton bonded to the oxygen atom in a trans arrangement of the skeleton. Three transition states were located that allow interconversion between the different isomers. Hardnesses and Fukui indices were calculated to follow the reactivity trend along the protonation path and to explain the preference for a particular protonation site on neutral sulfine. Proton affinity, gas-phase basicity and heat of formation values, obtained for the first time fully quantum mechanically, agree well with those derived by a recent mass spectrometry experimental study. Good agreement between density functional theory and previous high-level theoretical and experimental data was also found for the heat of formation of sulfine and its most stable protonated form. Received: 12 October 1998 / Accepted: 24 November 1998 / Published online: 16 March 1999     

Synthesis, crystal structure and density functional theoretical studies on phenyl-thiocarbamic acid-O-pyridin-4-ylmethyl ester

The title compound of phenyl-thiocarbamic acid-O-pyridin-4-ylmethyl ester has been synthesized and characterized by elemental analysis, IR, electronic spectroscopy and X-ray single crystal diffraction. Density functional theory (DFT) method calculations of the structure, atomic charge distributions and the thermodynamic properties at different temperatures have been performed. Calculated results show that DFT method at B3LYP/6-311G^** level can well reproduce the structure of the title compound. The predicted vibrational frequencies are compared with the experimental ones and they support each other on the whole. The atoms of sulfur, oxygen and the atom of nitrogen from pyridine ring all have bigger negative charges, which make the title compound become a multidentate organic ligand. The correlation equations of the thermodynamic properties of and with temperature are also obtained.     

Studies of chemical hardness and Fukui function using the exact solution of the density functional theory

Neal's procedure has been applied to determine the electron density ρ(x) for the H₂ molecule. The chemical hardness has been calculated employing the ab initio and density functional theory methods and the values are found to be reasonably good. The principle of maximum hardness (PMH) was tested. Fukui functions and the distribution of electron density around the internuclear distance were studied employing the electron density of the H₂ molecule. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 4–10, 2001