BAC-MP4 predictions of thermochemistry for gas-phase indium compounds in the In-H-C-O-Cl system

The BAC-MP4 methodology has been used to calculate the heats of formation and associated thermodynamic parameters for a family of 51 different one-, two-, and three-coordinate indium compounds. The ligands explored were H, CH3, C2H5, OH, OCH3, and Cl. The BAC-MP4 methodology was calibrated by both experimental data and high-level coupled-cluster calculations. The results agree well with the small amount of published experimental data on InHn and InCln species. A linear variation in compound heats of formation with ligand substitution is identified. For example, the heat of formation gradient, Delta(DeltaHf degrees)/Deltan, in the series InHn-3Cln (n=0-3), is -89.6 kcal/mol, with R=1.000. Additionally, trends in bond strengths, bond angles, and bond dissociation energies (BDEs) depending on ligand identity are identified. The BDEs decrease monotonically with increasing atomic number in group III and show similar variations among ligand systems for different metals. These observations and results can be combined with the data from previous studies to guide the selection of CVD precursors and experimental parameters for thin-film deposition of indium-containing materials for semiconductor and optical thin-film applications.     

Daubechies wavelets as a basis set for density functional pseudopotential calculations

Daubechies wavelets are a powerful systematic basis set for electronic structure calculations because they are orthogonal and localized both in real and Fourier space. We describe in detail how this basis set can be used to obtain a highly efficient and accurate method for density functional electronic structure calculations. An implementation of this method is available in the ABINIT free software package. This code shows high systematic convergence properties, very good performances, and an excellent efficiency for parallel calculations.     

Universal systematic sequence of even-tempered Gaussian primitive functions in electronic correlation studies

The possibility of developing a universal systematic sequence of eventempered Gaussian primitive functions for atomic and molecular electronic structure studies is examined. The radial beryllium-like ions are used to demonstrate this approach both within the Hartree-Fock model and by including correlation effects. Correlation energies are computed using the diagrammatic many-body perturbation theory. The Hartree extrapolation procedure is used to obtain empirical upper bounds to the basis set limit and the procedure of Schmidt and Ruedenberg is employed to obtain empirical lower bounds for the basis set limit. The convergence properties of the calculations with respect to the size of the basis set are examined.     

Opto-electronic properties and molecular design of new materials based on pyrrole studied by DFT

In this work, we investigate oligopyrroles and derivatives, which serve as models for corresponding polymers. In order to discuss these materials, we carried out DFT calculations and used DFT methods to calculate ground state electronic structures. We are particularly interested in exploring the potential of several substituent groups as electron donors with numerous ties to electronic materials by exploring and comparing the energies of HOMO, LUMO, Gap energies, and structural properties. Results are discussed in comparison with the properties of the doped oligomers. The theoretical ground-state geometry and electronic structure of the studied molecules were obtained by the DFT method at B3LYP level with 6-31G(d) basis set. The opto-electronic properties of these materials were determined by ZINDO/s and TD//B3LYP/6-31G(d) calculations performed on the B3LYP/6-31(d) optimized geometries. The results of this study demonstrate how electronic properties can be tuned by the backbone ring or side group and suggest these compounds as good candidates for opto-electronic applications.     

Carbene proton attachment energies: theoretical study

The geometries and electronic energies of six singlet carbenes, with methyl and phenyl substituents, and the corresponding carbenium ions were obtained using several density functional theory (DFT) variants and the second-order M?ller-Plesset method for electron correlation and compared with G3 results, with the aim to determine a relatively low-cost computational protocol that is sufficiently accurate for the specific molecules and ions of interest. Some additional calculations were performed at the CCSD(T) level. Results for diphenylcarbene, methylphenylcarbene, and their cations, which were not previously investigated by ab initio methods, are reported as are calculations on methylene, methylcarbene, dimethylcarbene, and phenylcarbene. The MPW3LYP/6-311+G(d,p) hybrid DFT level was found to give results that were in close agreement with those obtained using G3 theory, with a mean absolute deviation (MAD) of 1.76 kcal/mol for the calculated proton attachment energies (PAEs). Equilibrium geometries obtained with this method were compared with those obtained at the MP2/6-311G(d,p) level of theory, and bond lengths and bond angles had MADs of 0.005 A and 1.0 degrees, respectively. Harmonic vibrational frequencies of all the carbene molecules and the corresponding ions were computed to verify that the stationary points were true minima, to obtain zero-point corrected energies, to assist in infrared studies of the molecules. The recommended combination of method and basis set is expected to be a useful framework that uses modest amounts of computer resources to obtain usable thermochemical data on moderate-sized hydrocarbons and hydrocarbon cations, including coal-mimetic species.     

Statistical thermodynamics of 1-butanol, 2-methyl-1-propanol, and butanal

The purpose of the present investigation is to calculate partition functions and thermodynamic quantities, viz., entropy, enthalpy, heat capacity, and Gibbs free energies, for 1-butanol, 2-methyl-1-propanol, and butanal in the vapor phase. We employed the multi-structural (MS) anharmonicity method and electronic structure calculations including both explicitly correlated coupled cluster theory and density functional theory. The calculations are performed using all structures for each molecule and employing both the local harmonic approximation (MS-LH) and the inclusion of torsional anharmonicity (MS-T). The results obtained from the MS-T calculations are in excellent agreement with experimental data taken from the Thermodynamics Research Center data series and the CRC Handbook of Chemistry and Physics, where available. They are also compared with Benson's empirical group additivity values, where available; in most cases, the present results are more accurate than the group additivity values. In other cases, where experimental data (but not group additivity values) are available, we also obtain good agreement with experiment. This validates the accuracy of the electronic structure calculations when combined with the MS-T method for estimating the thermodynamic properties of systems with multiple torsions, and it increases our confidence in the predictions made with this method for molecules and temperatures where experimental or empirical data are not available.     

乙基溴及其阳离子的低能激发态从头算研究

应用ANO-S基组以及ECP基组在CASSCF理论水平下计算了乙基溴及其阳离子的低能激发态几何构型, 并应用CASPT2方法对动态相关能进行单点能校正. 根据乙基溴基态能量和相应阳离子电子态的能量差对光电子谱(Photoelectron spectrum)的谱线进行了理论指认. 在乙基溴的基态预测几何下, 进行了谐振频率计算, 对各个振动频率进行了理论指认. 计算结果与实验值符合得较好.     

Ab initio calculations of the lowest electronic states in the CuNO system

The lowest singlet and triplet electronic levels of the A' and A" symmetry species of the neutral copper-nitrosyl (CuNO) system are calculated by ab initio methods at the multi-reference configuration interaction (MRCI) level of theory with single and double excitations, and at the coupled cluster level of theory with both perturbational (CCSD(T)) and full inclusion of triple excitations (CCSDT). Experimental data are difficult to obtain, hence the importance of carrying out calculations as accurate as possible to address the structure and dynamics of this system. This paper aims at validating a theoretical protocol to develop global potential energy surfaces for transition metal nitrosyl complexes. For the MRCI calculations, the comparison of level energies at linear structures and their values from C(2v) and C(s) symmetry restricted calculations has allowed to obtain clear settings regarding atomic basis sizes, active orbital spaces and roots obtained at the multi-configurational self-consistent field (MCSCF) level of theory. It is shown that a complete active space involving 18 valence electrons, 11 molecular orbitals and the prior determination of 12 roots in the MCSCF calculation is needed for overall qualitatively correct results from the MRCI calculations. Atomic basis sets of the valence triple-zeta type are sufficient. The present calculations yield a bound singlet A' ground state for CuNO. The CCSD(T) calculations give a quantitatively more reliable account of electronic correlation close to equilibrium, while the MRCI energies allow to ensure the qualitative assessment needed for global potential energy surfaces. Relativistic coupled cluster calculations using the Douglas-Kroll-Hess Hamiltonian yield a dissociation energy of CuNO into Cu and NO to be (59 ± 5) kJ mol(-1) ((4940 ± 400) hc?cm(-1)). Favorable comparison is made with some of previous theoretical results and a few known experimental data.     

What are the most efficient basis set strategies for correlated wave function calculations of reaction energies and barrier heights?

As electronic structure methods are being used to obtain quantitatively accurate reaction energies and barrier heights for increasingly larger systems, the choice of an efficient basis set is becoming more critical. The optimum strategy for achieving basis set convergence can depend on the way that electron correlation is treated and can take advantage of flexibility in the order in which basis functions are added. Here we study several approaches for estimating accurate reaction energies and barrier heights from post-Hartree-Fock electronic structure calculations. First and second, we evaluate methods of estimating the basis set limit of second order Mo?ller-Plesset perturbation theory and of coupled cluster theory with single and double excitations and a quasiperturbative treatment of connected triple excitations by using explicitly correlated basis functions (in the F12a implementation) along with valence, polarization, and diffuse one-electron basis functions. Third, we test the scheme of adding a higher-order correction to MP2 results (sometimes called MP2∕CBS + ΔCCSD(T)). Finally, we evaluate the basis set requirements of these methods in light of comparisons to Weizmann-3.2, Weizmann-4, and CCSDT(2)(Q)∕CBS+CV+R results.     

An ab initio potential surface describing abstraction and exchange for H+CH4

We report an ab initio-based global potential energy surface for H+CH4 that describes the abstraction and exchange reactions. The PES, which is invariant with respect to any permutation of five H atoms, is a fit to 20,728 electronic energies calculated using the partially spin-restricted coupled-cluster method (RCCSD(T)) with a moderately large basis (aug-cc-pVTZ). A first set of quasiclassical trajectory calculations using this PES are reported for the H+CD4-->HD+CD3 reaction at collision energies of 0.65 and 1.52 eV and are compared to experiment and recent direct dynamics calculations done with density functional theory.     

Energy-consistent relativistic pseudopotentials for the 4d elements: atomic and molecular applications

Recently reported energy-consistent relativistic pseudopotentials have been used with series of matching correlation consistent basis sets in benchmark calculations of various atomic and molecular properties. The basis set convergence of the 4d metal electron affinities and 5s2-->5s0 excitation energies are reported at the CCSD(T) level of theory, and the effects of valence and 4s4p correlation are investigated. In addition the impact of correlating the low-lying 3d electrons was also studied in all-electron Douglas-Kroll-Hess (DKH) calculations, which also included the ionization potentials and 5s2-->5s1 excitation energies. For all four atomic properties, higher order coupled cluster calculations through CCSDTQ are reported. The final calculated values are generally all within 1 kcal/mol of experiment. A notable exception is the ionization potential of Tc, the currently accepted experimental value of which is suggested to be too high by about 3 kcal/mol. Molecular calculations are also reported for the low-lying electronic states of ZrO and RuF, as well as the ground electronic state of Pd2. The effects of spin-orbit coupling are investigated for these cases in pseudopotential calculations. Wherever possible, the pseudopotential results have been calibrated against DKH calculations with correlation consistent basis sets of triple-zeta quality. In all cases the calculated data for these species are in very good agreement with experiment. In particular, the correct electronic ground state for the RuF molecule (4Phi92) was obtained, which was made possible by utilizing systematic sequences of correlation consistent basis sets.     

Density functional theory study on molecular structure and vibrational IR spectra of isocytosine

Density functional theory with the combined Becke3-LYP exchange-correlation energy functional [DFT(B3-LYP) method] using the 6-31G(d, p) basis set is applied to predict molecular parameters (geometries, rotational constants, dipole moments) and vibrational IR spectra (harmonic wavenumbers, absolute intensities) of six tautomers of the isocytosine molecule. The results are compared with the corresponding data calculated at the conventional ab initio Hartree-Fock (HF) level using the same basis set and with available experimental data. Calculations show that (a) three amino tautomers are slightly nonplanar species with, evidently, a distorted amino group, (b) the DFT (B3-LYP)/6-31G(d, p) method predicts better molecular parameters, than do the HF calculations, and (c) the DFT(B3-LYP)-calculated vibrational IR spectra of isocytosine agree well with the available recorded IR spectra, and they show marked improvement over the IR spectra predicted at the HF/6-31G(d, p) level. Tautomeric stabilities of isocytosine are discussed on the basis of computed electronic energies by the DFT(B3-LYP) and ab initio approaches [including the MP2 and MP4(SDQ) calculations of electronic energies] and predicted zero-point vibrational energies by DFT(B3-LYP) and HF methods. This relative energies at 0 K of the tautomeric forms of isocytosine predicted by both conventional ab initio and DFT(B3-LYP) methods correlate well with the experimental data, showing the predominance of the aminohydroxy tautomer of isocytosine for an isolated molecule. © 1997 John Wiley & Sons, Inc.     

Structural and spectroscopic characterization of 2,3-difluorobenzoic acid and 2,4-difluorobenzoic acid with experimental techniques and quantum chemical calculations

In this study, the molecular conformation, vibrational and electronic transition analysis of 2,3-difluorobenzoic acid and 2,4-difluorobenzoic acid (C7H4F2O2) were presented using experimental techniques (FT-IR, FT-Raman and UV) and quantum chemical calculations. FT-IR and FT-Raman spectra in solid state were recorded in the region 4000-400 cm(-1) and 4000-5 cm(-1), respectively. The UV absorption spectra of the compounds that dissolved in ethanol were recorded in the range of 200-800 nm. The structural properties of the molecules in the ground state were calculated using density functional theory (DFT) and second order M?ller-Plesset perturbation theory (MP2) employing 6-311++G(d,p) basis set. Optimized structure of compounds was interpreted and compared with the earlier reported experimental values. The scaled vibrational wavenumbers were compared with experimental results. The complete assignments were performed on the basis of the experimental data and total energy distribution (TED) of the vibrational modes, calculated with scaled quantum mechanics (SQM) method. A study on the electronic properties, such as absorption wavelength, excitation energy, dipole moment and frontier molecular orbital energy, were performed by time dependent DFT (TD-DFT) approach. Based on the UV spectra and TD-DFT calculations, the electronic structure and the assignments of the absorption bands of steady compounds were discussed. The calculated HOMO and LUMO energies show that charge transfer occurs within the molecules.     

Extrapolating to the one-electron basis set limit in polarizability calculations

We report calculations of polarizabilities using total energies extrapolated to the complete basis set limit. A dual-level scheme has been employed, with the complete basis set limit of the correlation energy determined by the recently reported uniform singlet- and triplet-pair extrapolation method. The finite field approach has been employed, with tensors and averaged polarizabilities for the ground electronic states of H 2, N 2, CO, and H 2O reported and compared with available experimental data in the literature. Exploratory results are also presented for C 6H 4NO 2NH 2.     

An ab initio study of the vibronic, spin-orbit, and magnetic hyperfine structure in the X2Pi electronic state of NCO

In the present study we give the results of the ab initio calculations on the vibronic, spin-orbit, and magnetic hyperfine structure in the X (2)Pi electronic state of the NCO radical. The calculations of the potential surfaces and the electronic mean values of the hyperfine coupling constants are carried out by means of the density functional theory approach (B3LYP functional combined with an atomic orbital basis set suitable for calculations of the hyperfine structure). The vibronic levels, spin-orbit splitting, and the vibronic mean values of the components of the hyperfine tensor in the vibronic species are calculated using a variational method. The results of the calculations are in good agreement with the available experimental data.     

Vibrational energies of PH3 calculated variationally at the complete basis set limit

The potential energy surface for the electronic ground state of PH(3) was calculated at the CCSD(T) level using aug-cc-pV(Q+d)Z and aug-cc-pVQZ basis sets for P and H, respectively, with scalar relativistic corrections included. A parametrized function was fitted through these ab initio points, and one parameter of this function was empirically adjusted. This analytical PES was employed in variational calculations of vibrational energies with the newly developed program TROVE. The convergence of the calculated vibrational energies with increasing vibrational basis set size was improved by means of an extrapolation scheme analogous to the complete basis set limit schemes used in ab initio electronic structure calculations. The resulting theoretical energy values are in excellent agreement with the available experimentally derived values.     

BAC-MP4 predictions of thermochemistry for gas-phase antimony compounds in the Sb-H-C-O-Cl system

Calibrated by both experimental data and high-level coupled-cluster calculations, the BAC-MP4 methodology was applied to 51 SbL(n) (L = H, CH(3), C(2)H(5), Cl, and OH, n = 1-5) molecules, providing calculated heats of formation and associated thermodynamic parameters. These data identify a linear variation in heats of formation with ligand substitution, trends in bond dissociation energies (BDEs) with ligand identity [BDE(Sb-C(2)H(5)) < BDE(Sb-CH(3)) < BDE(Sb-H) < BDE(Sb-Cl) < BDE(Sb-OH)], and a monotonic decrease in BDE upon successive ligand elimination. The linear variation in BDE is consistent with the behavior of other group V elements, in contrast to the characteristic high-low-high trend of adjacent group III (In) and group IV (Sn) elements. Additionally, these data complement those of previous studies of metal-organic species and provide a foundation of thermochemical data that can aid in the selection of CVD precursors and deposition conditions for the growth of antimony-containing materials.     

Theoretical investigation of the dissociation dynamics of vibrationally excited vinyl bromide on an ab initio potential-energy surface obtained using modified novelty sampling and feed-forward neural networks

The reaction dynamics of vibrationally excited vinyl bromide have been investigated using classical trajectory methods on a neural network potential surface that is fitted to an ab initio database of 12 122 configuration energies obtained from electronic structure calculations conducted at the MP4(SDQ) level of theory using a 6-31G(d,p) basis set for the carbon and hydrogen atoms and Huzinaga's (43334334) basis set augmented with split outer s and p orbitals (4332143214) and a polarization f orbital with an exponent of 0.5 for the bromine atom. The sampling of the 12-dimensional configuration hyperspace of vinyl bromide prior to execution of the electronic structure calculations is accomplished by combining novelty-sampling methods, chemical intuition, and trajectory sampling on empirical and neural network surfaces. The final potential is obtained using a two-layer feed-forward neural network comprising 38 and 1 neurons, respectively, with hyperbolic tangent sigmoid and linear transfer functions in the hidden and output layers, respectively. The fitting is accomplished using the Levenberg-Marquardt algorithm with early stopping and Bayesian regularization methods to avoid overfitting. The interpolated potentials have a standard deviation from the ab initio results of 0.0578 eV, which is within the range generally regarded as "chemical accuracy" for the purposes of electronic structure calculations. It is shown that the potential surface may be easily and conveniently transferred from one research group to another. The files required for transfer of the vinyl bromide surface can be obtained from the Electronic Physics Auxiliary Publication Service. Total dissociation rate coefficients for vinyl bromide are obtained at five different excitation energies between 4.50 and 6.44 eV. Branching ratios into each of the six open reaction channels are computed at 24 vibrational energies in the range between 4.00 and 6.44 eV. The distribution of vibrational energies in HBr formed via three-center dissociation from vinyl bromide is determined and compared with previous theoretical and experimental results. It is concluded that the combination of ab initio electronic structure calculations, novelty sampling with chemical intuition and trajectories on empirical analytic surfaces, and feed-forward neural networks provides a viable framework in which to execute purely ab initio molecular-dynamics studies on complex systems with multiple open reaction channels.     

Using many-particle basis sets for electronic energy calculations of molecules

A method of calculating the electronic energies of molecules using a many-particle basis set is proposed. In this case, the many-particle Schrödinger equation may be solved without resorting to the one-electron approximation. The results of the electronic energy calculations of some two-electron systems are in good agreement with experimental results.     

Accurate ab initio binding energies of alkaline earth metal clusters

The effects of basis set superposition error (BSSE) and core-correlation on the electronic binding energies of alkaline earth metal clusters Y(n) (Y = Be, Mg, Ca; n = 2-4) at the Moller-Plesset second-order perturbation theory (MP2) and the single and double coupled cluster method with perturbative triples correction (CCSD(T)) levels are examined using the correlation consistent basis sets cc-pVXZ and cc-pCVXZ (X = D, T, Q, 5). It is found that, while BSSE has a negligible effect for valence-electron-only-correlated calculations for most basis sets, its magnitude becomes more pronounced for all-electron-correlated calculations, including core electrons. By utilizing the negligible effect of BSSE on the binding energies for valence-electron-only-correlated calculations, in combination with the negligible core-correlation effect at the CCSD(T) level, accurate binding energies of these clusters up to pentamers (octamers in the case of the Be clusters) are estimated via the basis set extrapolation of ab initio CCSD(T) correlation energies of the monomer and cluster with only the cc-pVDZ and cc-pVTZ sets, using the basis set and correlation-dependent extrapolation formula recently devised. A comparison between the CCSD(T) and density functional theory (DFT) binding energies is made to identify the most appropriate DFT method for the study of these clusters.