Contracted well-tempered Gaussian basis sets in SCF calculations on the ground and excited electronic states of neutral and ionized diatomic molecules containing first-row atoms

Summary Calculations were done on ground and excited states of C₂, C ₂ ⁺ , C ₂ ^– , N₂, N ₂ ⁺ , O₂, O ₂ ⁺ , O ₂ ^– , CO, CO⁺, CO²⁺, and CO^– using contracted well-tempered basis sets. The (14s 10p) basis sets were augmented with threed, one or twof, and oneg functions. Total energies, orbital energies, and spectroscopic constants were compared with the best available computational data.     

Quantum Monte Carlo study of ground and first excited state of C,N, O,F, and Ne atoms using Slater-Jastrow-Backflow wave function

All-electron fixed-node diffusion quantum Monte Carlo energies of the two lowest-lying states of C, N, O, F, and Ne atoms are reported. The Slater-Jastrow form is used as the trial wave function. We will use single- and multideterminant wave functions as the Slater part. The single-determinant wave function has been computed by the Hartree-Fock method and the multideterminant wave functions have been computed by the complete active space self-consistent field, configuration interaction with single and double excitation, configuration interaction with single, double, triple, and quadruple excitation and second-order configuration interaction. For the ground- and first excited states, the multideterminant wave functions have computed more than 99% of the correlation energy. Significant improvements have been achieved using the backflow transformations and up to 99.8% of the correlation energy has been recovered. A very good agreement with the experimental data has been obtained for the excitation energies.     

Prediction of spectroscopic constants for diatomic molecules in the ground and excited states using time-dependent density functional theory

Spectroscopic constants of the ground and next seven low-lying excited states of diatomic molecules CO, N2, P2, and ScF were computed using the density functional theory SAOP/ATZP model, in conjunction with time-dependent density functional theory (TD-DFT) and a recently developed Slater type basis set, ATZP. Spectroscopic constants, including the equilibrium distances r(e), harmonic vibrational frequency omega(e), vibrational anharmonicity omega(e)x(e), rotational constant B(e), centrifugal distortion constant D(e), the vibration-rotation interaction constant alpha(e), and the vibrational zero-point energy E(n)0 were generated in an effort to establish a reliable database for electron spectroscopy. By comparison with experimental values and a similar model with an established larger Slater-type basis set, et-QZ3P-xD, it was found that this model provides reliably accurate results at reduced computational costs, for both the ground and excited states of the molecules. The over all errors of all eight lowest lying electronic states of the molecules under study using the effective basis set are r(e)(+/-4%), omega(e)(+/-5% mostly without exceeding +/-20%), omega(e)x(e)(+/-5% mostly without exceeding 20%, much more accurate than a previous study on this constant of +/-30%), B(e)(+/-8%), D(e)(+/-10%), alpha(e)(+/-10%), and E(n)0(+/-10%). The accuracy obtained using the ATZP basis set is very competitive to the larger et-QZ3P-xD basis set in particular in the ground electronic states. The overall errors in r(e), omega(e)x(e), and alpha(e) in the ground states were given by +/-0.7, +/-10.1, and +/-8.4%, respectively, using the efficient ATZP basis set, which is competitive to the errors of +/-0.5, +/-9.2, and +/-9.1%, respectively for those constants using the larger et-QZ3P-xD basis set. The latter basis set, however, needs approximately four times of the CPU time on the National Supercomputing Facilities (Australia). Due to the efficiency of the model (TD-DFT, SAOP and ATZP), it will be readily applied to study larger molecular systems.     

Geometrical shapes of molecules by iterative extended Hückel method

The iterative extended Hückel method proposed by Harriset al. has been employed to calculate the shapes of 21 molecules consisting of H, C, N, and O. It has given fairly reasonable results, especially when it has been combined with a new formula; the formula regards the orbital energies and the wave functions of the Hückel method as those of the Hartree-Fock method.     

A chelating resin containing S, N and O atoms: Synthesis and adsorption properties for Hg(II)

A novel chelating resin containing S, N and O atoms (PSME-EDA) was synthesized by using poly(2-hydroxyethylmercaptomethylstyrene) (PSME) and diethanolamine (EDA) as materials. Its structure was characterized by elemental analysis, Fourier transform-infrared spectra (FTIR), scanning electron microscopy (SEM) and X-ray diffraction (XRD). The adsorption of the resin for Hg²⁺ was investigated. The saturated adsorption capacity of PSME-EDA for Hg²⁺ could reach to about 1.1 mmol/g at 25 °C when the initial Hg²⁺ concentration was 0.02 mol/l. Some factors affecting the adsorption such as temperature, reaction time and ion concentration were also studied. The results showed that adsorption was controlled by liquid film diffusion. The increasing of temperature was beneficial to adsorption. The Langmuir model was better than the Freundlich model to describe the isothermal process. The values of ΔG, ΔH, and ΔS calculated at 25 °C were −7.99 kJ mol⁻¹, 22.5 kJ mol⁻¹ and 34.4 J mol⁻¹ K⁻¹, respectively. The adsorption mechanism of PSME-EDA resin for Hg(II) was confirmed by X-ray photoelectron spectroscopy (XPS).     

Comparison of the capacity of nitrogen, oxygen, sulfur, and selenium atoms to exist in onium states. quantum-chemical investigation of isomeric model systems with two heteroatoms

HF, B3LYP, and MP2 calculations with the 6-31+G(d) basis set with correction to the energy of zero-point vibrations were carried out to determine the energy characteristics of model molecules containing two heteroatoms in the sp ³- and sp ²-hybrid states; different combinations of N, O, S, and Se atoms were studied. The stability of the onium states of the nitrogen atom was found to be greater than for its chalcogen analogs and the relative stability of onium states of the chalcogen analogs was found to depend on the hybridization of these atoms. Analysis of these results permitted us to construct a stability series of onium derivatives and to interpret the positional selectivity in electrophilic substitution reactions of five-membered heterocyclic compounds with one heteroatom. To J. Stradins, an outstanding and tireless scientist, with our deep respect and sincere affection. Translated from Khimiya Geterotsiklicheskikh Soedinenii, No. 12, pp. 1801–1808, December, 2008.     

Proton affinities of molecules containing nitrogen and oxygen: Comparing density functional results to experiment

Summary Proton affinities were calculated using density functional theory for 11 small molecules whose primary protonation site is on nitrogen, and eight small molecules that protonate on oxygen. Calculations were performed using both the local spin density approximation and nonlocal gradient corrections to the exchange correlation functional. The results were not sensitive to whether the nonlocal gradient correction was implemented on the final local spin density optimized geometry or whether the correction was included in the self-consistent calculation of the energy at each optimization step. Although negligible basis set dependence was found using the analytic Gaussian basis sets, numerical basis sets required augmentation by a double set of polarization functions to achieve reasonable agreement with experiment. All calculations systematically underestimated oxygen proton affinities.     

Extension of a simplified method for molecular correlation energy calculations to molecules containing third row atoms II. Application to HCl,HCl+, ClO and NCl molecules

A simplified method for molecular correlation energy calculations developed in I is applied to the determination of the potential curves of some diatomic chlorinated molecules HCl(¹⁺), HCl⁺(²II), ClO(²II) and NCl(X ^3–, a ¹, b ¹⁺). Dissociation energies, vibrational frequencies and equilibrium internuclear distances are derived from these curves. The ionisation potential (I _p) of HCl, the proton affinity (PA) of Cl, and the term energies of the excited states of the NCl radical are also calculated. It is shown that the results are very sensitive to correlation effects and that the correlated results converge to the corresponding experimental values within 10% for PA, D _e , T _e and _e, 2% for R _e and 0.3% for I _p. This agreement allows us to predict the following lower limits for the dissociation energies of the NCl radical: 2.14, 3.28 and 2.47 eV respectively for the X^3–, a ¹ and b ¹⁺ states. Results on HF and HF⁺ are also discussed and compared with those obtained for HCl and HCl⁺.Boursier IRSIA     

The structure of CX2YNO (X, Y=F, Cl) molecules in the ground and lowest excited singlet electronic states

Ab initio quantum-chemical calculations of equilibrium geometric parameters, vibrational frequencies, and potentials of internal rotation for CCIF₂NO and CCl₂FNO molecules in the ground (S₀) and lowest excited singlet (S₁) electronic states were performed. The results of calculations were compared with experimental data. A new interpretation of experimental spectra of the CCIF₂NO molecule was suggested. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 8, pp. 1453–1458, August, 1999.     

Comparative analysis of hydrogen bonding with participation of the nitrogen, oxygen and sulfur atoms in the 2(2'-heteroaryl)pyrroles and their trifluoroacetyl derivatives based on the 1H, 13C, 15N spectroscopy and DFT calculations

The N-H...X (X = N,O,S) intramolecular hydrogen bond in the series of 2(2'-heteroaryl)pyrroles and their trifluoroacetyl derivatives is examined by the (1)H, (13)C, (15)N spectroscopy and density functional theory (DFT) calculations. The influence of the hydrogen bond on coupling and shielding constants is considered. It is shown that the N-H...N intramolecular hydrogen bond causes a larger increase in the absolute size of the (1)J(N,H) coupling constant and a larger deshielding of the bridge proton than the N-H...O hydrogen bond. The effect of the N-H...S interaction on the (1)J(N,H) coupling constant and the shielding of the bridge proton is small. The NMR parameter changes in the series of the 2(2'-heteroaryl)pyrroles due to N-H...X hydrogen bond and the series of the 1-vinyl-2-(2'-heteroaryl)-pyrroles due to C-H...X hydrogen bond have the same order. The proximity of the nitrogen, oxygen or sulfur lone pair to the F...H hydrogen bridge quenches the trans-hydrogen bond spin-spin couplings (1h)J(F,H-1) and (2h)J(F,N).     

Theoretical study of structures of the X<Subscript>2</Subscript>CO and XYCO molecules (X and Y = H,F, or Cl) in the ground and lowest excited triplet electronic states

Systematic calculations of the structures of the H₂CO, F₂CO, Cl₂CO, HClCO, HFCO, and FClCO molecules in the S₀ and T₁ states were performed using the B3LYP and MP2 methods with different AO basis sets and also at the CCSD(T)/cc-pV(T+d)Z level of theory. The saturation of the correlation consistent sequence of basis sets cc-pV(N+d)Z (N = D, T, Q, and 5) and aug-cc-pV(N+d)Z (N = D, T, and Q) was studied. Recommendations for choosing the calculation method are given. The relativistic corrections were estimated. The influence of the number and type of halogen atoms on the geometric parameters of the molecules in the S₀ and T₁ states and the heights of inversion barriers in the T₁ state was investigated. Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 12, pp. 2625–2635, December, 2005.     

CEPA calculations on open-shell molecules. II. Singlet-triplet energy splitting in π2 configurations of diatomic molecules

The CEPA-PNO method is used for calculating the energy difference ΔE _ST between the³∑⁻ and the¹Δ states of diatomic molecules in electronic π² configurations. An analysis of the contribution of electron correlation to ΔE _ST is performed in terms of physically understandable effects such as direct correlation, dynamic spin polarization, semiinternal and internal excitations. It is shown that these effects are of completely different importance for the molecules treated in this study: For C₂ the direct correlation between the two singly occupied π-orbitals is the dominant correlation contribution to ΔE _ST; for O₂, S₂, SO the internal excitation π _u ² → π _g ² is predominant, whereas for NH and PH there is a close competition between the direct correlation and the spin polarization of the underlying σ-orbitals. The basis set dependence of these effects is investigated, in particular for NH. Our final results reproduce experimental values of ΔE _ST within 0.05–0.10 eV.     

Molecular parameters of tetraatomic carbonyls X2CO and XYCO (X,Y = H,F, Cl) in the ground and lowest excited electronic states,part 1: A test of ab initio methods

Geometrical parameters of tetraatomic carbonyl molecules X₂CO and XYCO (X, Y = H, F, Cl) in the ground (S₀) and lowest excited singlet (S₁) and triplet (T₁) electronic states as well as values of barriers to inversion in S₁ and T₁ states and S₁ ← S₀ and T₁ ← S₀ adiabatic transition energies were systematically investigated by means of various quantum‐chemical techniques. The following methods were tested: HF, MP2, CIS, CISD, CCSD, EOM‐CCSD, CCSD(T), CR‐EOM‐CCSD(T), CASSCF, MR‐MP2, CASPT2, CASPT3, NEVPT2, MR‐CISD, and MR‐AQCC within cc‐pVTZ and cc‐pVQZ basis sets. The accuracy of quantum‐chemical methods was estimated in comparison with experimental data and rather accurate structures of excited electronic states were obtained. MP2 and CASPT2 methods appeared to be the most efficient and CCSD(T), CR‐EOM‐CCSD(T), and MR‐AQCC the most accurate. It was found that at equilibrium all the molecules under study are nonplanar in S₁ and T₁ electronic states with CO out‐of‐plane angle ranging from 34° (H₂CO, S₁) to 52° (F₂CO, T₁), and height of barrier to inversion varying from 300 (H₂CO, S₁) to 11,000 (F₂CO, T₁) cm^?1. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Temperature dependence of X- and Q-band EPR spectra of the dinuclear manganese(II) complex [(NO2Bpmp)Mn2(mu-OAc)2]+: determination of the exchange constant and of the spin parameters for the S=1, 2, and 3 spin states

A new dinuclear manganese(II) complex was synthesised with the biscompartimental ligand 2,6-bis[bis(2-pyridylmethyl)aminomethyl]-4-nitrophenol (NO(2)BpmpH) and characterised by X-ray crystallography. Magnetic susceptibility measurements revealed that the two high-spin Mn(II) ions are antiferromagnetically coupled with a singlet-to-triplet separation of 7.2 cm(-1). The powder EPR spectra were recorded for both X- and Q-bands between 1.8 K and 35 K. A detailed analysis of these spectra led to the determination of three out of five individual spin-state zero-field splitting parameters. From the proposed simulations, the exchange coupling constant J and the intermetallic distance have been computed.     

Synthesis of monohydridohexamethylbenzeneruthenium(II) complexes containing group V donor ligands. Isomerism arising from cyclometallation of a tertiary phosphine at aliphatic and aromatic carbon atoms

The η-hexamethylbenzenehydridoruthenium(II) complexes RuHCl(η-C₆Me₆)L (L = PPh₃ (11), AsPh₃ (12), P(C₆H₄-p-F)₃ (14), P(C₆H₄-p-Me)₃ (15), P(C₆H₄-p-OMe)₃ (16), P-t-BuPh₂ (17), P-i-PrPh₂ (18), P-i-Pr₃ (19), PCy₃ (20) and P-t-BuMe₂ (21)) have been made by heating [RuCl₂(η-C₆Me₆)]₂, the ligand and sodium carbonate in propan-2-ol. The triarylphosphine complexes 11, 14 and 15 react with methyllithium to give aryl ortho-metallated hydridoruthenium(II) complexes such as $\text{RuH(}\text{o}\text{-C}{}_6\text{H}{}_4\text{P}$Ph₂)(η-C₆Me₆) (22) and 19 similarly gives the isopropyl cyclometallated complex $\text{RuH(CH}{}_2\text{CHMeP}$-i-Pr₂(η-C₆Me₆) (29) as a mixture of diastereomers. Reaction of 17 with methyllithium gives initially the t-butyl cyclometallated complex $\text{RuH(CH}{}_2\text{CMe}{}_2\text{P}$Ph₂)(η-C₆Me₆) (25) which isomerizes by a first order process (k$\text{0}\text{?}$.2 h^?1 in C₆D₆ or THF-d₈ at 50°C) to the aryl ortho-metallated complex $\text{RuH(}\text{o}\text{-C}{}_6\text{H}{}_4\text{P}$-t-BuPh)(η-C₆Me₆) (26). The similarly generated isopropyl cyclometallated complex $\text{RuH(CH}{}_2\text{CHMeP}$Ph₂)(η-C₆Me₆) (27) has not been isolated in a pure state owing to rapid isomerization to $\text{RuH(}\text{o}\text{-C}{}_6\text{H}{}_4\text{P}$-i-PrPh)(η-C₆Me₆) (28); both 27 and 28 exist as a pair of diastereomers. The formation of the cyclometallated complexes and the isomerizations are thought to involve intermediate 16-electron ruthenium(O) complexes Ru(η-C₆Me₆)L.