Molecular structure of thiourea

The molecular structure of thiourea has been investigated under C(s), C(2), and C(2v) symmetry constraints. At the coupled-cluster level in conjunction with a triple-ζ basis set, only the C(2) conformer has been found to be a real minimum on the potential energy surface. Its equilibrium structure has therefore been accurately evaluated using both theoretical and experimental data. With respect to the former, high-level quantum-chemical calculations at the coupled-cluster level in conjunction with correlation-consistent basis sets ranging in size from triple- to quintuple-zeta have been carried out. Extrapolation to the complete basis-set limit as well as core-correlation effects and inclusion of full treatment of triple excitations in the cluster operator have been considered. On the basis of the vibrational ground-state rotational constants available for five isotopic species and the corresponding computed vibrational corrections, the semiexperimental equilibrium geometry of thiourea has also been determined for the first time.     

Rotational spectra of rare isotopic species of bromofluoromethane: determination of the equilibrium structure from ab initio calculations and microwave spectroscopy

Guided by theoretical predictions, the rotational spectra of the mono- and bideuterated species of bromofluoromethane, CDH(79)BrF, CDH(81)BrF, CD(2) (79)BrF, and CD(2) (81)BrF, have been recorded for the first time. Assignment of a few hundred rotational transitions led to the accurate determination of the ground-state rotational constants, all of the quartic and most of the sextic centrifugal distortion constants, as well as the full bromine quadrupole-coupling tensor for both (79)Br and (81)Br, in good agreement with corresponding theoretical predictions based on high-level coupled-cluster calculations. The rotational spectra of the (13)C containing species (13)CH(2) (79)BrF and (13)CH(2) (81)BrF have been observed in natural abundance and have been assigned, thus allowing the determination of the rotational and centrifugal distortion constants as well as the bromine quadrupole-coupling tensor. Furthermore, empirical equilibrium structures have been obtained within a least-squares fit procedure using the available experimental ground-state rotational constants for various isotopic species. Vibrational effects have been accounted for in the analysis using vibration-rotation interaction constants derived from anharmonic force fields computed at the second-order Moller-Plesset perturbation theory as well as coupled-cluster (CC) levels. The empirical equilibrium geometries obtained in this way agree well with the corresponding theoretical predictions obtained from CC calculations [at the CCSD(T) level] after extrapolation to the complete basis set limit and inclusion of core-valence correlation corrections and relativistic effects.     

Equilibrium structure and spectroscopic constants of maleic anhydride

The quadratic, cubic, and semi-diagonal quartic force fields of maleic anhydride have been calculated at the MP2 level of theory employing the cc-pVTZ basis set. The spectroscopic constants derived from the force field are in excellent agreement with the corresponding experimental values. The semi-experimental equilibrium structure has been derived from experimental ground state rotational constants and rovibrational corrections calculated from the cubic force field. This semi-experimental equilibrium structure is in excellent agreement with the ab initio structures computed at the CCSD(T) level of theory and it is closer to the ab initio structure than the purely experimental (or empirical) structures r ₀, r _m⁽¹⁾, and r _m⁽²⁾ obtained by microwave spectroscopy as well as the equilibrium structure derived from gas-phase electron diffraction data.     

Potential energy surface and rovibrational states of the ground Ar-HI complex

A potential energy surface for the ground electronic state of the Ar-HI van der Waals complex is calculated at the coupled-cluster with single and double excitations and a noniterative perturbation treatment of triple excitations [CCSD(T)] level of theory. Calculations are performed using for the iodine atom a correlation consistent triple-zeta valence basis set in conjunction with large-core Stuttgart-Dresden-Bonn relativistic pseudopotential, whereas specific augmented correlation consistent basis sets are employed for the H and Ar atoms supplemented with an additional set of bond functions. In agreement with previous studies, the equilibrium structure is found to be linear Ar-I-H, with a well depth of 205.38 cm(-1). Another two secondary minima are also predicted at a linear and bent Ar-H-I configurations with well depths of 153.57 and 151.57 cm(-1), respectively. The parametrized CCSD(T) potential is used to calculate rovibrational bound states of Ar-HI/Ar-DI complexes, and the vibrationally averaged structures of the different isomers are determined. Spectroscopic constants are also computed from the CCSD(T) surface and their comparison with available experimental data demonstrates the quality of the present surface in the corresponding configuration regions.     

Equilibrium structures of heterocyclic molecules with large principal axis rotations upon isotopic substitution

Equilibrium structures, r(e), of the heterocyclic molecules oxirane, furazan, furan, ethylene ozonide, and 1,3,4-oxadiazole have been determined using three different, somewhat complementary techniques: a completely experimental technique (r(m)), a semiexperimental technique (r(e)(SE), whereby equilibrium rotational constants are derived from experimental effective ground-state rotational constants and corrections based principally on an ab initio cubic force field), and an ab initio technique (r(e)(BO), whereby geometry optimizations are usually performed at the coupled cluster level of theory including single and double excitations augmented by a perturbational estimate of the effects of connected triple excitations [CCSD(T)] using quadruple-ζ Gaussian basis sets). All these molecules are asymmetric tops with the moment of inertia I(c) much larger than the other two moments of inertia, I(a) and I(b). Molecules of this shape experience a large rotation of the principal axis system upon certain isotopic substitutions. For such isotopologues it is difficult to obtain a good structural fit to the semiexperimental moments of inertia I(a) and I(b), which may significantly reduce the accuracy of the r(e)(SE) structural parameters. The origin of this difficulty is explained. For the heavy-atom skeleton of these molecules it was possible to determine a rather accurate empirical mass-dependent structure without a priori knowledge of the equilibrium structure.     

CVRQD ab initio ground-state adiabatic potential energy surfaces for the water molecule

The high accuracy ab initio adiabatic potential energy surfaces (PESs) of the ground electronic state of the water molecule, determined originally by Polyansky et al. [Science 299, 539 (2003)] and called CVRQD, are extended and carefully characterized and analyzed. The CVRQD potential energy surfaces are obtained from extrapolation to the complete basis set of nearly full configuration interaction valence-only electronic structure computations, augmented by core, relativistic, quantum electrodynamics, and diagonal Born-Oppenheimer corrections. We also report ab initio calculations of several quantities characterizing the CVRQD PESs, including equilibrium and vibrationally averaged (0 K) structures, harmonic and anharmonic force fields, harmonic vibrational frequencies, vibrational fundamentals, and zero-point energies. They can be considered as the best ab initio estimates of these quantities available today. Results of first-principles computations on the rovibrational energy levels of several isotopologues of the water molecule are also presented, based on the CVRQD PESs and the use of variational nuclear motion calculations employing an exact kinetic energy operator given in orthogonal internal coordinates. The variational nuclear motion calculations also include a simplified treatment of nonadiabatic effects. This sophisticated procedure to compute rovibrational energy levels reproduces all the known rovibrational levels of the water isotopologues considered, H(2) (16)O, H(2) (17)O, H(2) (18)O, and D(2) (16)O, to better than 1 cm(-1) on average. Finally, prospects for further improvement of the ground-state adiabatic ab initio PESs of water are discussed.     

Femtosecond degenerate four-wave mixing of carbon disulfide: high-accuracy rotational constants

Femtosecond degenerate four-wave mixing (fs-DFWM) rotational coherence spectroscopy (RCS) has been used to determine the rotational and centrifugal distortion constants of the 00 (0)0 ground and 01 (1)0 vibrationally excited states of gas-phase CS(2). RCS transients were recorded over the 0-3300 ps optical delay range, allowing the observation of 87 recurrences. The fits yield rotational constants B(00 (0)0)=3.271 549 2(18) GHz for (12)C(32)S(2) and B(00 (0)0)=3.175 06(21) GHz for the (12)C(32)S(34)S isotopomer. The rotational constants of the degenerate 01 (1)0 bending level of (12)C(32)S(2) are B(01 (1)0)=3.276 72(40) and 3.279 03(40) GHz for the e and f substrates, respectively. These fs-DFWM rotational constants are ten times more accurate than those obtained by CO(2) laser/microwave heterodyne measurements and are comparable to those obtained by high-resolution Fourier transform infrared spectroscopy. Ab initio calculations were performed at two levels, second-order Moller-Plesset theory and coupled-cluster singles, doubles, and iterative triples [CCSD(T)]. The equilibrium and vibrationally averaged C=S distances were calculated using large Dunning basis sets. An extrapolation procedure combining the ab initio rotational constants with the experiment yields an equilibrium C=S bond length of 155.448 pm to an accuracy of +/-20 fm. The theoretical C=S bond length obtained by a complete basis set extrapolation at the CCSD(T) level is r(e)(C=S)=155.579 pm, or 0.13 pm longer than that in the experiment.     

Rotational spectra of rare isotopic species of fluoroiodomethane: determination of the equilibrium structure from rotational spectroscopy and quantum-chemical calculations

Supported by accurate quantum-chemical calculations, the rotational spectra of the mono- and bi-deuterated species of fluoroiodomethane, CHDFI and CD(2)FI, as well as of the (13)C-containing species, (13)CH(2)FI, were recorded for the first time. Three different spectrometers were employed, a Fourier-transform microwave spectrometer, a millimeter/submillimter-wave spectrometer, and a THz spectrometer, thus allowing to record a huge portion of the rotational spectrum, from 5 GHz up to 1.05 THz, and to accurately determine the ground-state rotational and centrifugal-distortion constants. Sub-Doppler measurements allowed to resolve the hyperfine structure of the rotational spectrum and to determine the complete iodine quadrupole-coupling tensor as well as the diagonal elements of the iodine spin-rotation tensor. The present investigation of rare isotopic species of CH(2)FI together with the results previously obtained for the main isotopologue [C. Puzzarini, G. Cazzoli, J. C. Lo?pez, J. L. Alonso, A. Baldacci, A. Baldan, S. Stopkowicz, L. Cheng, and J. Gauss, J. Chem. Phys. 134, 174312 (2011); G. Cazzoli, A. Baldacci, A. Baldan, and C. Puzzarini, Mol. Phys. 109, 2245 (2011)] enabled us to derive a semi-experimental equilibrium structure for fluoroiodomethane by means of a least-squares fit procedure using the available experimental ground-state rotational constants together with computed vibrational corrections. Problems related to the missing isotopic substitution of fluorine and iodine were overcome thanks to the availability of an accurate theoretical equilibrium geometry (computed at the coupled-cluster singles and doubles level augmented by a perturbative treatment of triple excitations).     

Rotational spectra of 1-chloro-2-fluoroethylene. II. Equilibrium structures of the cis and trans isomer

Equilibrium structures for the cis and trans isomer of 1-chloro-2-fluoroethylene are reported. The structures are obtained within a least-squares fit procedure using the available experimental ground-state rotational constants for various isotopic species of both forms. Vibrational effects were eliminated before the analysis using vibration-rotation interaction constants derived from computed quadratic and cubic force fields with the required quantum chemical calculations carried out using second-order Moller-Plesset perturbation as well as coupled-cluster (CC) theory. The semiexperimental or empirical equilibrium geometries obtained in this way agree well with the corresponding theoretical predictions obtained from CC calculations [at the CCSD(T) level] after extrapolation to the complete basis-set limit and inclusion of core-valence correlation corrections. The present results allow a detailed analysis of the geometrical differences between the two forms of 1-chloro-2-fluoroethylene. They are also compared to the structural data available for other halogenated ethylenes.     

Equilibrium structure and torsional barrier of BH3NH3

Born-Oppenheimer equilibrium structures, r(e)(BO), of the electronic ground state of the borazane (BH3NH3) molecule of C3v point-group symmetry are computed ab initio using the CCSD(T) method with basis sets up to quintuple-zeta quality. Inclusion of the counterpoise correction and extrapolation of the structural parameters to the complete basis set limit yield a best estimate of r(e)(BO) of BH3NH3. The anharmonic force field of BH3NH3, computed at the CCSD(T) level of theory with a basis set of triple-zeta quality, allows the determination of semi-experimental equilibrium rotational constants, which in turn result in a semi-experimental equilibrium structure, r(e)(SE). The r(e)(BO) and r(e)(SE) structures are in excellent agreement, indicating the validity of the methods used for their determination. The empirical mass-dependent structure, r(m)(1), of BH3NH3 is also determined. Although it is inferior in quality to the previous two structures, it is much more accurate than the standard empirical r0 and r(s) structures reported earlier for BH3NH3. The semi-experimental r(e)(SE) as well as the empirical r(m)(1) structures determined are based on experimental ground-state rotational constants available from the literature for nine isotopologues of borazane. The effective barrier to the internal rotation of BH3NH3, a molecule isoelectronic with CH3CH3, has been computed ab initio, employing the focal-point analysis (FPA) approach, to be 699 +/- 11 cm(-1). This compares favorably with an empirical redetermination of the effective barrier based on the above r(e)(SE) structure, V3 = 718(17) cm(-1).     

Anharmonic force field of cis- and trans-formic acid from high-level ab initio calculations, and analysis of resonance polyads

The quadratic, cubic, and semidiagonal quartic force fields of cis- and trans-formic acid have been calculated using three different levels of theory. They all give satisfactory results, including the one at the lowest level of theory which is the MP2 method employing a basis set of triple-zeta quality. The results are used to theoretically analyze resonance polyads, including the one involving the 4(1), 5(1), 6(1), 8(1), 7(1)9(1), and 9(2) vibrational states. A semiexperimental equilibrium structure is derived from experimental ground state rotational constants and rovibrational interaction parameters calculated from the ab initio force field. The ab initio structure calculated at the CCSD(T) level of theory using a basis set of quintuple-zeta quality is in excellent agreement with the semiexperimental structure.     

Rotational spectra and equilibrium structures of H2SiS and Si2S

The rotational spectra of two small silicon sulfides, silanethione H(2)SiS and the disilicon sulfide ring Si(2)S, have been detected in the centimeter band by Fourier transform microwave spectroscopy of a molecular beam; lines of H(2)SiS were also observed in the millimeter band up to 377 GHz in a glow discharge. Precise rotational and centrifugal distortionconstants have been determined for the normal and a number of the more abundant rare isotopic species of both closed-shell molecules. Theoretical equilibrium (r(e)) structures of H(2)SiS and Si(2)S were derived from coupled-cluster calculations that included triple and quadruple excitations, core correlation, and extrapolation to the basis-set limit. The r(e) structures agree to within 5×10(-4) A? and 0.1(°) with empirical equilibrium (r(e)(emp)) structures derived from the experimental rotational constants, combined with theoretical vibrational and electronic corrections. Both H(2)SiS and Si(2)S are good candidates for radioastronomical detection in the circumstellar shells of evolved carbon-rich stars such as IRC+10216, because they are fairly polar and are similar in composition to the abundant astronomical molecule SiS.     

The case of the weak N-X bond: Ab initio, semi-experimental, and experimental equilibrium structures of XNO (X = H, F, Cl, OH) and FNO2

The equilibrium structures of FNO, ClNO, HONO, and FNO2 have been determined using three different, somewhat complementary methods: a completely experimental, a semi-experimental (where the equilibrium rotational constants are derived from the experimental effective ground-state rotational constants and an ab initio cubic force field), and an ab initio, where geometry optimizations are usually performed at the coupled cluster level of nonrelativistic electronic structure theory using small to very large Gaussian basis sets. For the sake of comparison, the equilibrium structures of HNO and N2O have also been redetermined, confirming and extending earlier results. The semi-experimental method gives structural parameters in good agreement with the reliable experimental results for each compound investigated. Because of inadequate treatment of electron correlation, the single-reference CCSD(T) method gives N-X (X[double bond]F, Cl, OH) bonds that are too strong and associate bond lengths that are significantly too short. The discrepancy increases with increase in the size of the basis set. A much more elaborate treatment of electron correlation at the CCSDTQ level solves this problem and results in increased bond lengths, correctly representing the weakness of the N-X bond in these XNO and XNO2 species. The equilibrium structures determined are accurate to better than 0.001 A and 0.1 degrees .     

Spectroscopic investigation of fluoroiodomethane, CH2FI: Fourier-transform microwave and millimeter-/submillimeter-wave spectroscopy and quantum-chemical calculations

Guided by theoretical predictions, the rotational spectrum of fluoroiodomethane, CH(2)FI, has been recorded and assigned. Accurate values are reported for the ground-state rotational constants, all quartic, sextic, and two octic centrifugal-distortion constants. The hyperfine structure of the rotational spectrum was thoroughly investigated using a Fourier-transform microwave spectrometer and the Lamb-dip technique in the millimeter-/submillimeter-wave region, thus allowing the accurate determination of the complete iodine quadrupole-coupling tensor and of the diagonal elements of the iodine spin-rotation tensor. Relativistic effects turned out to be essential for the accurate theoretical prediction of the dipole moment and quadrupole-coupling constants and were accounted for by direct perturbation theory and a spin-free four-component treatment based on the Dirac-Coulomb Hamiltonian. The relativistic corrections to the dipole moment amount to up to 34% and to the iodine quadrupole-coupling tensor to about 15-16% of the total values.     

Extending the molecular size in accurate quantum-chemical calculations: the equilibrium structure and spectroscopic properties of uracil

The equilibrium structure of uracil has been investigated using both theoretical and experimental data. With respect to the former, quantum-chemical calculations at the coupled-cluster level in conjunction with a triple-zeta basis set have been carried out. Extrapolation to the basis set limit, performed employing the second-order M?ller-Plesset perturbation theory, and inclusion of core-correlation and diffuse-function corrections have also been considered. Based on the available rotational constants for various isotopic species together with corresponding computed vibrational corrections, the semi-experimental equilibrium structure of uracil has been determined for the first time. Theoretical and semi-experimental structures have been found in remarkably good agreement, thus pointing out the limitations of previous experimental determinations. Molecular and spectroscopic properties of uracil have then been studied by means of the composite computational approach introduced for the molecular structure evaluation. Among the results achieved, we mention the revision of the dipole moment. On the whole, it has been proved that the computational procedure presented is able to provide parameters with the proper accuracy to support experimental investigations of large molecules of biological interest.     

Accurate ab initio structure, dissociation energy, and vibrational spectroscopy of the F(-)-CH4 anion complex

Accurate equilibrium structure, dissociation energy, global potential energy surface (PES), dipole moment surface (DMS), and the infrared vibrational spectrum in the 0-3000 cm(-1) range of the F(-)-CH4 anion complex have been obtained. The equilibrium electronic structure calculations employed second-order M?ller-Plesset perturbation theory (MP2) and coupled-cluster (CC) method up to single, double, triple, and perturbative quadruple excitations using the aug-cc-p(C)VXZ [X = 2(D), 3(T), 4(Q), and 5] correlation-consistent basis sets. The best equilibrium geometry has been obtained at the all-electron CCSD(T)/aug-cc-pCVQZ level of theory. The dissociation energy has been determined based on basis set extrapolation techniques within the focal-point analysis (FPA) approach considering (a) electron correlation beyond the all-electron CCSD(T) level, (b) relativistic effects, (c) diagonal Born-Oppenheimer corrections (DBOC), and (d) variationally computed zero-point vibrational energies. The final D(e) and D0 values are 2398 +/- 12 and 2280 +/- 20 cm(-1), respectively. The global PES and DMS have been computed at the frozen-core CCSD(T)/aug-cc-pVTZ and MP2/aug-cc-pVTZ levels of theory, respectively. Variational vibrational calculations have been performed for CH4 and F(-)-CH4 employing the vibrational configuration interaction (VCI) method as implemented in Multimode.     

Two-electron reduced density matrices from the anti-Hermitian contracted Schrodinger equation: enhanced energies and properties with larger basis sets

Two-electron reduced density matrices (2-RDMs) have recently been directly determined from the solution of the anti-Hermitian contracted Schrodinger equation (ACSE) to obtain 95%-100% of the ground-state correlation energy of atoms and molecules, which significantly improves upon the accuracy of the contracted Schrodinger equation (CSE) [D. A. Mazziotti, Phys. Rev. Lett. 97, 143002 (2006)]. Two subsets of the CSE, the ACSE and the contraction of the CSE onto the one-particle space, known as the 1,3-CSE, have two important properties: (i) dependence upon only the 3-RDM and (ii) inclusion of all second-order terms when the 3-RDM is reconstructed as only a first-order functional of the 2-RDM. The error in the 1,3-CSE has an important role as a stopping criterion in solving the ACSE for the 2-RDM. Using a computationally more efficient implementation of the ACSE, the author treats a variety of molecules, including H2O, NH3, HCN, and HO3-, in larger basis sets such as correlation-consistent polarized double- and triple-zeta. The ground-state energy of neon is also calculated in a polarized quadruple-zeta basis set with extrapolation to the complete basis-set limit, and the equilibrium bond length and harmonic frequency of N2 are computed with comparison to experimental values. The author observes that increasing the basis set enhances the ability of the ACSE to capture correlation effects in ground-state energies and properties. In the triple-zeta basis set, for example, the ACSE yields energies and properties that are closer in accuracy to coupled cluster with single, double, and triple excitations than to coupled cluster with single and double excitations. In all basis sets, the computed 2-RDMs very closely satisfy known N-representability conditions.     

Anharmonic rovibrational analysis for disilacyclopropenylidene (Si2CH2)

The global minimum on the Si(2)CH(2) electronic singlet potential energy surface has been theoretically predicted to be a peculiar hydrogen bridged (Si···H···Si) disilacyclopropenylidene structure (Si(2)CH(2)). An accurate quartic force field for Si(2)CH(2) has been determined employing ab initio coupled-cluster theory with single and double excitations and a perturbative treatment for triple excitations [CCSD(T)], in combination with the correlation consistent core-valence quadruple zeta (cc-pCVQZ) basis set. The vibration-rotation coupling constants, equilibrium and zero-point vibration corrected rotational constants, centrifugal distortion constants, and harmonic and fundamental vibrational frequencies for six isotopologues of Si(2)CH(2) are predicted using vibrational second-order perturbation theory (VPT2). The anharmonic corrections for the vibrational motions involving the H bridged bonds are found to be more than 5% with respect to the corresponding harmonic vibrational frequencies. In this light, an experimental detection and characterization of disilacyclopropenylidene (Si(2)CH(2)) is highly desired.     

Lowest electronic states of neutral and ionic LiN

We have investigated the potential energy curves (PECs) of the LiN heteronuclear diatomic molecule, including its ionic species LiN⁺ and LiN⁻, using explicitly correlated multi-reference configuration interaction (MRCI-F12) calculations in conjunction with the correlation consistent quintuple-𝜁 basis set. The effect of core–valence correlation, scalar relativistic effects, and the size of the basis sets has been investigated. A comprehensive set of spectroscopic constants determined based on the above-mentioned calculations are also reported for the lowest electronic states and all systems, including dissociation energies, harmonic and anharmonic vibrational frequencies, and rotational constants. Additional parameters, such as the dipole moments, equilibrium spin-orbit constants, excitation energies, and rovibrational energy levels, are also documented. We found that the three triplet states of LiN, namely, X ³∑⁻, A ³Π, and 2 ³∑⁻, exhibit substantial potential wells in the PEC diagrams, while the quintet states are repulsive in nature. The ground state of the anion also shows a deep potential well in the vicinity of its equilibrium geometry. In contrast, the ground and excited states of the cation are very loosely bound. Charge transfer properties of each of these states are also analyzed to obtain an in-depth understanding of the interatomic interactions. We found that the core–valence correlation has a substantial effect on the calculated spectroscopic constants.     

Spectroscopic constants and molecular properties of the X2Π and a4Σ− electronic states of the SiF radical

The potential energy curves (PECs) of the X²Π and a⁴Σ^? electronic states of the SiF radical have been studied by an ab initio quantum chemical method. The calculations have been made using the complete active space self‐consistent field (CASSCF) method, which is followed by the valence internally contracted multireference configuration interaction (MRCI) approach in combination with several correlation‐consistent basis sets. The effects on the PECs by the core‐valence correlation and relativistic corrections are included. The way to consider the relativistic correction is to use the third‐order Douglas–Kroll Hamiltonian approximation. The relativistic corrections are made at the level of cc‐pV5Z basis set. The core‐valence correlation corrections are performed using the cc‐pCV5Z basis set. To obtain more reliable results, the PECs determined by the MRCI calculations are also corrected for size‐extensivity errors by means of the Davidson modification (MRCI+Q). These PECs are extrapolated to the complete basis set limit by the total‐energy extrapolation scheme. Using these PECs, the spectroscopic parameters are determined and compared with those reported in the literature. With these PECs obtained by the MRCI+Q/CV+DK+56 calculations, the vibrational levels, inertial rotation, and centrifugal distortion constants of the first 20 vibrational state of each electronic state are calculated when the rotational quantum number J equals zero. Comparison with the Rydberg‐Klein‐Rees (RKR) data shows that the present results are reliable and accurate. The molecular constants of the X²Π and a⁴Σ^? electronic states determined by the MRCI+Q/CV+DK+56 calculations should be good prediction for future laboratory experiment. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011