A self-consistent molecular field theory for aggregates of neutral molecules. I

A self-consistent theory is presented for aggregates of neutral molecules. According to the LCAO Hartree-Fock formalism a set of effective Hartree-Fock equations for molecules in the aggregate is derived. The molecular orbitals of each molecule are to be determined from the effective H-F equation for the molecule in which the interactions between the molecule and the surrounding ones are included as an intermolecular interaction field (molecular field). A self-consistent treatment leads to the molecular orbitals which are self-consistent with the molecular field. By this method, then-molecule problem becomesn times of one-molecule problem.     

Grid-based Thomas-Fermi-Amaldi equation with the molecular cusp condition

First, the Thomas-Fermi-Amaldi (TFA) equation was formulated with a newly derived condition to remove the singularities at the nuclei, which coincided with the molecular cusp condition. Next, the collocation method was applied to the TFA equation using the grid-based density functional theory. In this paper, the electron densities and the radial probabilities for specific atoms (He, Be, Ne, Mg, Ar, Ca) were found to agree with those from the Thomas-Fermi-Dirac (TFD) method. Total energies for specific atoms (He, Ne, Ar, Kr, Xe, Rn) and molecules (H2,CH4) were also found to be close to those from the Hartree-Fock method using the Pople basis set 6-311G relative to the TFD method. In addition, the computational expense to determine the electron density and its corresponding energy for a large scale structure, such as a carbon nanotube, is shown to be much more efficient compared to the conventional Hartree-Fock method using the 6-31G Pople basis set.     

Infrared and Raman spectra and theoretical study of methyl trifluoromethyl sulfone, CF3SO2CH3

The infrared and Raman spectra were obtained for liquid CF3SO2CH3, as well as the infrared spectrum of the gaseous substance. The molecular geometry was optimized by means of the Hartree-Fock (HF), second order electron correlation (MP2) and density functional theory (DFT) procedures of quantum chemistry, resulting in a structure with Cs symmetry. The wavenumbers corresponding to the normal modes of vibration were calculated using the DFT (B3LYP/6-31G**) approximation and their agreement with the measured values improved after scaling of the associated force field. An assignment of bands is proposed on the basis of such calculations and the comparison with related molecules.     

Matrix isolation study of the interaction of excited argon atoms with methyl cyanide,vibrational and electronic spectra of ketenimine

When methyl cyanide mixed with an excess of argon is codeposited at 14 K with a beam of argon atoms that has been excited by a low power microwave discharge, infrared and ultraviolet absorptions of several previously unidentified products appear. The most prominent set of absorptions is assigned to ketenimine, previously tentatively identified as the product of the reaction of NH with C₂H₂ in an argon matrix. Using a molecular geometry resulting from a recent ab initio calculation and a valence force field with four interaction constants, it has been possible to obtain a satisfactory least-squares force constant fit to the infrared data for seven isotopic species of ketenimine. Two electronic band systems are also reported for ketenimine, which is photodissociated by 2537-A radiation. The mechanism by which ketenimine is formed may involve an initial electron transfer from excited argon to methyl cyanide. Spectrospic data are also considered for the other products, one of them tentatively identified as CH₂CH, which differ in their behavior on mercury-are photolysis of the sample.     

On the calculation of dipole moment and polarizability derivatives by the analytical energy gradient method: Application to the formaldehyde molecule

Following the suggestion of Komornicki and McIver we have implemented an efficient computational scheme for the evaluation of dipole moment and polarizability derivatives at the Hartree-Fock SCF level. The derivatives are obtained by utilizing the analytical gradients of the molecular energy, calculated in the presence of an external electric field, with respect to the atomic cartesian coordinates, which are differentiated numerically with respect to the field. The implementation of the method within the framework of the MOLECULE program is discussed, concentrating on such aspects as numerical accuracy, utilization of molecular symmetry and computational efficiency. As an application, the dipole moment and polarizability derivatives of the formaldehyde molecule have been calculated, yielding infrared intensities and Raman scattering activities in the double harmonic approximation. The theoretical results are compared with the available experimental data; the agreement is satisfactory given the inherent restrictions of the SCF model.     

New approach to the correlation problem: Electron interaction potential as a variable in solving the many-particle Schrödinger equation

The many-electron wave function is represented as the product of the wave function of the independent particles and the function that depends only on the value of the interelectron interaction potential. The function defines the electron correlation effects; a standard linear differential equation was derived to define the function. The equation depends on the functions of independent particles; a generalization of the Hartree-Fock equations including electron correlation was obtained for these functions. The total energy calculation of two-electron ions shows that even solving an ordinary differential equation for the function of independent particles represented by the functions of noninteracting electrons leads to higher accuracy than the one achieved in the Hartree-Fock theory.     

Calculation of the electric hypershielding at the nuclei of molecules in a strong magnetic field

The third-rank electric hypershielding at the nuclei of 14 small molecules has been evaluated at the Hartree-Fock level of accuracy, by a pointwise procedure for the geometrical derivatives of magnetic susceptibilities and by a straightforward use of its definition within the Rayleigh-Schrodinger perturbation theory. The connection between these two quantities is provided by the Hellmann-Feynman theorem. The magnetically induced hypershielding at the nuclei accounts for distortion of molecular geometry caused by strong magnetic fields and for related changes of magnetic susceptibility. In homonuclear diatomics H(2), N(2), and F(2), a field along the bond direction squeezes the electron cloud toward the center, determining shorter but stronger bond. It is shown that constraints for rotational and translational invariances and hypervirial theorems provide a natural criterion for Hartree-Fock quality of computed nuclear electric hypershielding.     

Analytic first and second derivatives of the energy in the fragment molecular orbital method combined with molecular mechanics

Analytic first and second derivatives of the energy are developed for the fragment molecular orbital method interfaced with molecular mechanics in the electrostatic embedding scheme at the level of Hartree-Fock and density functional theory. The importance of the orbital response terms is demonstrated. The role of electrostatic embedding upon molecular vibrations is analyzed, comparing force field and quantum mechanical treatments for an ionic liquid and a solvated protein. The method is applied for 100 protein conformations sampled in molecular dynamics (MD) to take into account the complexity of a flexible protein structure in solution, and a good agreement with experimental data is obtained: Frequencies from an experimental infrared (IR) spectrum are reproduced within 17 cm⁻¹ .     

Some remarks on the SCRF theory of solvent effects and the calculation of proton potentials

Some aspects concerning the self-consistent reaction field theory of solvent effects are discussed. In particular, the variational solution to the non-linear Schrödinger equation is considered; a necessary and sufficient constraint to be added to the standard variational procedure is discussed. The exact solution of the non-linear equation is presented within the molecular orbital approach; correlation defaults to the Hartree-Fock like solutions are stated. Some thermodynamical correspondences are established with the magnitudes calculated with the self-consistent reaction field theory. Finally, we have commented upon the proton potentials calculated within this theory. An INDO calculation of a water trimer has been used as an example to discuss different types of proton translocation potentials.     

Local approach to coupled cluster evaluation of polarizabilities for long conjugated molecules

Electron correlation effect in conjugated polymer is a long-standing problem, especially for the nonlinear optical properties. We have implemented a spin-adapted Coupled Cluster with singles and doubles excitation with local molecular orbital approach. As a first application, we evaluate the static polarizability of conjugated polyene chains with finite field approach. It is found that the local molecular orbital approach can tremendously reduce the computational costs at sufficiently high accuracy. It is also found that the electron correlation can largely reduce the molecular polarizability with respect to the Hartree-Fock mean field results.     

Wave functions derived from experiment. V. Investigation of electron densities,electrostatic potentials,and electron localization functions for noncentrosymmetric crystals

The constrained Hartree-Fock method using experimental X-ray diffraction data is extended and applied to the case of noncentrosymmetric molecular crystals. A new way to estimate the errors in derived properties as a derivative with respect to added Gaussian noise is also described. Three molecular crystals are examined: ammonia [NH(3)], urea [CO(NH(2))(2)], and alloxan [(CO)(4)(NH)(2)]. The energetic and electrical properties of these molecules in the crystalline state are presented. In all cases, an enhancement of the dipole moment is observed upon application of the experimental constraint. It is found that the phases of the structure factors are robustly determined by the constrained Hartree-Fock model, even in the presence of simulated noise. Plots of the electron density, electrostatic potential, and the electron localization function for the molecules in the crystal are displayed. In general, relative to the Hartree-Fock model, there is a depletion of charge around hydrogen atoms and lone pair regions, and a build-up of charge within the molecular framework near nuclei, directed along the bonds. The electron localization function plots reveal an increase in the pair density between vicinal hydrogen atoms.     

The effect of electron correlation on one-electron properties in the 2 3S and 2 1S excited states of the helium atom

The density difference ΔD(r₁), defined as the difference between the exact and the Hartree-Fock radial density functions is reported for the 2 ³S and the 2 ¹S excited states of the helium atom. The density differences in both cases resemble the inverse of the helium ground-state density difference, but are more long-ranged due to the diffuseness of the excited states. The effect of electron correlation on one-electron distributions is also demonstrated by comparison of (r₁ⁿ) expectation values and several indices of the relative sizes of atoms. Electron correlation leads to a contraction of the Hartree-Fock electron clouds.     

The electron density distribution in CN−, LiCN and LiNC. The use of minimal and extended basis set SCF calculations

Electron density maps are reported for the CN^?ion and the LiCN and LiNC molecules, calculated from molecular wave-functions near the Hartree-Fock limit. The electron density distribution derived from CNDO/ 2 wavefunctions does not resemble the ab initio results. The ultimate ability of a minimal basis set to represent the electron density near the Hartree-Fock limit, has been tested. The requirement of N-representability of the trial electron density has been satisfied. It is found that the molecular valence density cannot be reproduced to a satisfactory extent by a minimal set of Slater orbitals, even when the exponents of the basis orbitals are optimized.     

On the theory of solvent effects

A simple electrostatic analysis is given of the virtual charge (solvaton) model to represent the environment effect on the electronic wave function of a solute immersed in a polarizable surrounding. New features of this model are found. The classical aspects are discussed and secondly the quantal implications are considered. A correct Hartree-Fock-like operator is derived which represents an electron in a molecular orbital subjected to the average effect of the other electrons and to the reaction field produced by the virtual charges on the atomic centers.A general formalism based on the preceding model is presented in appendix. The final equations have a form similar to Newton's equation to represent a solvated electron. Unlike some other theories in this field, there is no cut-off involved in the evaluation of the molecular integrals.     

外电场作用下寡聚苯分子导线的性质

利用Hartree-Fock 方法在6-31G*水平上对聚苯分子进行了计算研究. 分别从几何构型、分子轨道空间分布和分子轨道能级三个方面讨论了外电场对寡聚苯分子导线的影响, 给出了分子导线的性质与外电场的关系. 进一步, 连接硫原子于聚苯分子的两端, 并共价结合在金电极上. 利用非平衡格林函数方法对其在0-2.0 V 偏压下电子输运特征进行了深入研究.     

A chemometric analysis of ab initio vibrational frequencies and infrared intensities of methyl fluoride

Factorial design and principal component analyses are applied to CH₃F infrared frequencies and intensities calculated from ab initio wave functions. In the factorial analysis, the quantitative effects of changing from a 6–31G to a 6–311G basis, of including polarization and diffuse orbitals, and of correcting for electron correlation using the second-order Møller-Plesset procedure are determined for all frequencies and intensities. The most significant main effect observed for the frequencies corresponds to the shift from Hartree-Fock to MP2 calculations, which tends to lower all frequency values by approximately 100 cm⁻¹. For the intensities, the main effects are larger for the CF stretching and the CH₃ asymmetric stretching modes. Interaction effects between two or more of the four factors are found to be of minor importance, except for the interaction between correlation and polarization. The principal component analysis indicates that wave functions with polarization and diffuse orbitals at the second-order Møller-Plesset level provide the best estimates for the harmonic frequencies, but not for the intensities. For the frequencies, the first principal component distinguishes between MP2 and Hartree-Fock calculations, while the second component separates the wave functions with polarization orbitals from those without these orbitals. For the intensities, the separation is similar but less well defined. This analysis also shows that wave function optimization to calculate accurate intensities is more difficult than an optimization for frequencies. © 1996 by John Wiley & Sons, Inc.     

A study of the effect of the synthesis conditions of titanium dioxide on its morphology and cell toxicity properties

In this study, titanium dioxide nanoparticles (NPs) were synthesized using the home microwave method, and the effect of the microwave irradiation time on the structure of NPs was investigated. In addition, the morphological effect of these NPs on the toxicity of HDMSCs cells was investigated. The crystalline structure and morphology of the NPs were analyzed using X-ray diffraction (XRD), Fourier transform infrared spectroscopy (FTIR), and field emission scanning electron microscopy (FE-SEM); the cytotoxicity was determined by the methyl thiazolyl tetrazolium (MTT) assay. X-ray diffraction analysis revealed that all thin films had a polycrystalline nature with an anatase phase of TiO₂. It was also found that the crystallite size increased with increasing microwave radiation time. The FTIR spectrum showed Ti-O-Ti properties by the peak in the range between 527 and 580 cm^?1. Further, the FE-SEM images showed that the grain size increased with increasing irradiation time. The MTT assay results showed that the accumulation of NPs leads to toxicity.     

Collective solvent coordinates for the infrared spectrum of HOD in D2O based on an ab initio electrostatic map

An ab initio MP2 vibrational Hamiltonian of HOD in an external electrostatic potential parametrized by the electric field and its gradient-tensor is constructed. By combining it with the fluctuating electric field induced by the D(2)O solvent obtained from molecular dynamics simulations, we calculate the infrared absorption of the O-H stretch. The resulting solvent shift and infrared line shape for three force fields (TIP4P, SPC/E, and SW) are in good agreement with the experiment. A collective coordinate response for the solvent effect is constructed by identifying the main electrostatic field and gradient components contributing to the line shape. This allows a realistic stochastic Liouville equation simulation of the line shapes which is not restricted to Gaussian frequency fluctuations.