Relativistic and correlation effects on molecular properties. I. The dihalogens F 2,Cl 2,Br 2,I 2, and At 2
|
|
- Katherine Shelton
- 6 years ago
- Views:
Transcription
1 Relativistic and correlation effects on molecular properties. I. The dihalogens F 2,Cl 2,Br 2,I 2, and At 2 L. Visscher Laboratory of Chemical Physics and Material Science Center, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands K. G. Dyall Thermosciences Institute, NASA Ames Research Center, Moffett Field, California Received 17 November 1995; accepted 6 March 1996 A benchmark study of a number of relativistic correlation methods is presented. Bond lengths, harmonic frequencies, and dissociation energies of the molecules F 2,Cl 2,Br 2,I 2, and At 2 are calculated at various levels of theory, using both the Schrödinger and the Dirac Coulomb Gaunt Hamiltonian American Institute of Physics. S I. INTRODUCTION The additivity of relativistic and electron correlation effects has long been questioned. Until recently few methods were available to quantitatively assess the interplay of both deviations of the nonrelativistic mean-field theory in molecular calculations. The so-called scalar relativistic methods 1,2 usually assume that corrections of higher order than c 2 can be neglected for chemical accuracy. Furthermore it is assumed that the effect of the spin orbit coupling on the form of the orbitals may be neglected, allowing a partition of the Hamiltonian into a spin-independent and a spin-dependent part. The latter part is then only used in the final stage of the calculation to couple the correlated many-electron functions. Since all these methods can be derived from the manyelectron Dirac Coulomb or Dirac Coulomb Breit equation, a straightforward way to provide a benchmark for the different approximations involved is to solve this equation directly. We here use relativistic analogues of a number of standard ab initio methods, which allows us to also study the influence of the choice of method to the calculated relativistic effects. The starting point is the Dirac Hartree Fock method, from which we get a set of 4-component molecular spinors. These spinors provide a natural basis for the socalled no-pair treatment of electron correlation. This means we do not consider states that have one or more positron electron pairs. In our treatment we do not include additional terms derived from the underlying theory, quantum electrodynamics QED. From experience with atomic calculations these corrections are known to have only very small differential effects on the valence shell properties. We, however, do consider the effect of the two-electron Gaunt interaction that describes the magnetic interaction between the electrons. Estimates of the importance of this term on the properties calculated will be given on basis of perturbation theory. The methods used are Hartree Fock HF, second-order Mo ller Plesset perturbation theory MP2, configuration interaction with single and double substitutions CISD, coupled cluster with single and double substitutions CCSD and the latter method, perturbatively corrected for the effect of triple excitations CCSD T. This gives a range of correlation treatments, from no electron correlation in the HF method to a fairly high level of correlation in the CCSD T method. II. COMPUTATIONAL DETAILS Nonrelativistic calculations were performed using GAUSSIAN92/DFT, 3 MOLPRO94 Ref. 4, and the nonrelativistic option of MOLFDIR. 5 All 4-component relativistic calculations were also performed with the latter code, using a newly-developed relativistic CCSD T code. 6 Two different models of the nuclei were used. The nonrelativistic calculations were done using the conventional point-nucleus model. Since this unphysical model gives rise to a weak singularity at the origin for the relativistic wave function, we used a Gaussian charge distribution to represent the spatial extent of the nucleus in the relativistic calculations. The exponents used in this work were E 09, E 09, E 09, E 09, and E 09 for F, Cl, Br, I, and At, respectively. The speed of light in atomic units was taken to be All molecular calculations were performed using D 4h symmetry, except for MOLPRO94, where D 2h symmetry was used. The atomic calculations were carried out in O h.to prevent spurious discrepancies between the nonrelativistic and relativistic dissociation energies, we calculated both the nonrelativistic and the relativistic atomic asymptotes in a basis of spinors optimized for the spherical average of the four 2 P 3/2 states. In the MP2, CI, and CC calculations the halogens were treated as 7 valence electron atoms, i.e., correlating only the valence s and p-electrons. The virtual space used was truncated slightly by removing the virtual spinors that consisted of corelike functions. Spectroscopic constants were in general obtained by fitting the potential energy curve to a 4th order polynomial in the internuclear distance. For the nonrelativistic calculations on F 2,Cl 2, and Br 2 we used the frequency option of GAUSSIAN J. Chem. Phys. 104 (22), 8 June /96/104(22)/9040/7/$ American Institute of Physics
2 L. Visscher and K. G. Dyall: Relativistic effects on dihalogens. I 9041 TABLE I. Summary of the basis set sizes. L, large component basis; S, small component basis. pvdz pvtz Primitive Contracted Primitive Contracted F L: 9s,5p,1d L: 3s,2p,1d L: 10s,6p,2d,1f L: 4s,3p,2d,1f S: 5s,10p,5d,1f S: 2s,4p,2d,1f S: 6s,12p,7d,2f,1g S: 3s,6p,4d,2f,1g Cl L: 12s,8p,1d L: 4s,4p,1d L: 15s,9p,2d,1f L: 5s,5p,2d,1f S: 8s,13p,8d,1f S: 3s,5p,4d,1f S: 9s,17p,10d,2f,1g S: 4s,7p,6d,2f,1g Br L: 15s,12p,6d L: 5s,5p,2d L: 17s,14p,9d,1f L: 6s,7p,3d,1f S: 12s,21p,12d,6f S: 4s,7p,6d,3f S: 14s,26p,15d,9f,1g S: 5s,9p,8d,4f,1g I L: 16s,12p,7d L: 6s,7p,3d L: 21s,17p,12d,1f L: 7s,8p,4d,1f S: 12s,23p,12d,7f S: 5s,9p,8d,4f S: 17s,33p,18d,12f,1g S: 6s,11p,10d,5f,1g At L: 23s,20p,15d,10f L: 7s,9p,5d,1f L: 25s,22p,17d,11f L: 7s,10p,6d,2f S: 20s,23p,20d,15f,10g S: 6s,9p,11d,6f,2g S: 22s,25p,22d,17f,11g S: 7s,11p,11d,7f,3g III. BASIS SETS The basis sets used for F and Cl were the correlation consistent polarized valence cc-pvdz and cc-pvtz basis sets of Dunning and co-workers. 7,8 The fluorine basis sets were augmented by a tight p function, with exponents and for the cc-pvdz and cc-pvtz basis sets, respectively. This addition reduced the error made in the 2p spin orbital splitting as compared to the numerical Dirac Hartree Fock value by a factor of 10 to less than 1%. For the Cl basis no additional functions were required to obtain this accuracy in the 3p spin orbital splitting. Both basis sets for both atoms were recontracted on basis of numerical Dirac Hartree Fock calculations using a modification 9 of GRASP. 10 For Cl one extra contracted large component p-function relative to the nonrelativistic contraction scheme was used to describe the 2 p spin orbit splitting. The small component primitive basis set was generated using unrestricted kinetic balance and contracted according to the atomic balance scheme. For Br we used the basis sets of Schäfer et al. 11 and F,gri 12 to supply a double and triple zeta description of the valence basis. The former basis set was extended with a tight s and a tight p function, exponents E 06 and E 04, while the latter, in addition to a tight s and a tight p function with exponents E 07 and 5.728E 05 was also extended with a correlating d and f function, exponents 0.24 and These exponents were obtained by optimizing the average of the nonrelativistic CISD energies of Br and Br. The outermost primitive function in the original d set was also used for correlation to produce a (2d,1f ) correlating set. The basis sets were contracted using a similar scheme as for the F and Cl basis sets. The iodine double zeta basis was that of Stromberg et al. 13 extended with a tight s and p function and a correlating d function, with exponents 4.35E 06, 1.85E 04, and For the iodine triple zeta description we used F,gri s basis set, 14 again extended in a similar manner to the Br basis with a tight s and p, exponents E 08 and E 05, and a correlating d and f, exponents 0.22 and 0.39, and recontracted as was done for Br. Since the use of extended nonrelativistic primitive basis sets for very heavy atoms is questionable and the number of small component primitive functions generated by the kinetic balance procedure becomes a burden, we followed a somewhat different procedure to obtain the At basis sets. We optimized two sets of exponents at the Dirac Hartree Fock level using a version of GRASP developed by Dyall and F,gri 15 under the constraints that the d exponents from subset of the s exponents and that the f exponents from a subset of the p exponents. This ensures a reasonably compact small component basis set without relaxing the kinetic balance concept. The resulting basis sets were contracted according to the cc-pvxz scheme described above. After completion of the calculations a flaw was discovered in the open-shell exponent optimization routine. This means that the exponents were not completely optimized, though the description is still better than would be obtained with a pure nonrelativistic basis set. A fully optimized set of exponents is to be published shortly, the basis sets used in this work are available upon request. For most of the nonrelativistic calculations we used the primitive basis sets as described above, but without the additional tight s and p-functions. For F 2 we also used the extended cc-pvdz basis, because comparing the nonrelativistic results in the original basis with the relativistic results in the extended basis gave rise to a spurious difference of 3 cm 1 in the frequency. The sizes of the primitive and contracted basis sets are summarized in Table I. IV. RESULTS AND DISCUSSION Calculated spectroscopic constants for the molecules F 2, Cl 2,Br 2,I 2, and At 2 are given in Tables II VI. Table VII gives the relativistic effect on these properties, calculated as the difference x rel x nonrel. The relativistic decrease of the dissociation energy of the molecules is clearly seen in all molecules. This is driven by the spin orbit coupling in the valence p-shell. For the lighter elements the atomic spin orbit coupling is almost completely quenched in the molecule. This gives a decrease of the atomic asymptotes relative to the molecule of two times 1/3 of the valence p-splitting. For all halogen molecules we find a relativistic destabilization of the bond. This destabilization, which is also found
3 9042 L. Visscher and K. G. Dyall: Relativistic effects on dihalogens. I TABLE II. Properties of F 2, calculated at various levels of theory. TABLE IV. Properties of Br 2, calculated at various levels of theory. NR-HF DC-HF DC G-HF NR-MP DC-MP NR-CISD DC-CISD NR-CCSD DC-CCSD NR-CCSD T DC-CCSD T Experiment a NR-HF DC-HF DC G-HF NR-MP DC-MP NR-CISD DC-CISD NR-CCSD DC-CCSD NR-CCSD T DC-CCSD T Experiment a in other p-block elements like thallium and bismuth, 16 can be explained by considering the hybridization of the orbitals involved. The six possible linear combinations of the atomic p-orbitals may be labeled 1/2g, 1/2u, 3/2u, 1/2u *, 1/2g *, and 3/2g *. The orbitals are degenerate in the nonrelativistic limit, but since in the relativistic case the 1/2 and 1/2 are allowed to mix, six nondegenerate orbitals will be found. The bonding and antibonding orbitals are partly and partly with the amount of mixing determined by the spin orbit splitting. If the energy difference between the p 1/2 and p 3/2 is small, the 5 occupied linear combinations will be a pure -bond 1/2g, two pure -bonds 1/2u and 3/2u and 2 pure -antibonds 1/2g * and 3/2g *. The LUMO is the pure -antibonding 1/2u * orbital and the total bond order is one -bond. This is the case for F 2. On the other hand, if the spin orbit splitting is large the hybridization energy will make the formation of pure - or -bonds less favorable. In the hyper-relativistic limit no hybridization will occur and we combine the atomic p 1/2 s with each other and the atomic p 3/2 s also separately. The highest (1/2u) combination is the TABLE III. Properties of Cl 2, calculated at various levels of theory. Molecular NR-MP2, NR-CCSD, and NR-CCSD T values were taken from Ref. 18. NR-HF DC-HF DC G-HF NR-MP DC-MP NR-CISD DC-CISD NR-CCSD DC-CCSD NR-CCSD T DC-CCSD T Experiment a p 3/2,1/2 p 3/2,1/2 which is one-third -bonding and two-thirds -antibonding. This makes the 3/2g * a more favorable LUMO since it is pure -antibonding. The bond order is now one -bond, which means that the effect of relativity is to change the bond from pure into pure. In astatine the hyper-relativistic limit is not yet reached; the lowest determinant that can be formed is still the one with the (1/2u) LUMO. The amount of mixing is, however, appreciable as can be seen by looking at the amount of *-character in the occupied (1/2u) orbital. This percentage grows from 0.00% in F, 0.01% in Cl, 0.19% in Br, 1.10% in I to 10.20% in At. A. Fluorine For the lighter elements comparison may be made with a wealth of nonrelativistic calculations. We initially compared our relativistic results for F 2 to the nonrelativistic calculations done by Peterson et al., 17 since they employed the same basis set. This gave rise to a surprisingly large relativistic effect of minus 1 to 2 cm 1 in the frequency for the ccpvdz calculations. Since errors in the fitting procedure are of the order of 0.05 cm 1, the difference between fits of TABLE V. Properties of I 2, calculated at various levels of theory. NR-HF DC-HF DC G-HF NR-MP DC-MP NR-CISD DC-CISD NR-CCSD DC-CCSD NR-CCSD T DC-CCSD T Experiment a
4 L. Visscher and K. G. Dyall: Relativistic effects on dihalogens. I 9043 TABLE VI. Properties of At 2, calculated at various levels of theory. the data to 4th and 6th order polynomials, the difference should originate from another source. We found this deviation to depend to a large extent on the extra tight p function that was added in the relativistic calculation. Redoing the nonrelativistic calculations with the extended basis gave frequencies that were 3 cm 1 smaller than reported by Peterson et al. The relativistic effect changed sign and was found to be maximally 1.7 cm 1 on the CCSD T level. This relativistic frequency shift is significantly larger than at the CCSD T /cc-pvtz level, where it is 0.3 cm 1. Adding the tight p-function to the triple zeta basis and redoing the calculation only alters the result by 0.1 cm 1. The relativistic effect on the frequency is thus exaggerated due to the basis set incompleteness and no conclusion can be drawn on basis of the double zeta results. On the triple zeta level the relativistic frequency shift is more reliable and found to be smaller than 1 cm 1. The relativistic effect on the dissociation energy is 0.8 to 0.9 kcal mol 1 which is in good agreement with the value of 0.77 kcal mol 1 that is obtained by correcting the nonrelativistic result only for the spin orbit effect on the asymptote this approximation will hereafter be referred to as ASO-only. For these light elements the two-electron relativistic effects are of the same order as the one-electron effects. Adding the Gaunt interaction as a perturbative correction to the two-electron integrals at the HF level shifts the frequency by 0.3 cm 1 and increases the bond length by Å, resulting in total relativistic shift of 0.8 cm 1 and Å. Its effect on the D e is 0.05 kcal mol 1. B. Chlorine NR-HF DC-HF DC G-HF NR-MP DC-MP NR-CISD DC-CISD NR-CCSD DC-CCSD NR-CCSD T DC-CCSD T The nonrelativistic Cl 2 results were taken from Woon and Dunning, 18 except for the dissociation energies. For this number we recalculated the atomic Cl energy to avoid spurious effects from the use of double-group symmetry in the relativistic calculation. This affects the D e on the CCSD level by 0.2 kcal mol 1 ; we find 39.5 kcal mol 1, while Woon and Dunning report 39.3 kcal mol 1. Again the relativistic effect is slightly dependent on the basis set used. The cc-pvdz basis shows no significant effect on the bond length and frequency, while the cc-pvtz basis gives a small TABLE VII. Relativistic effects on the properties of the halogens at various levels of theory. Molecule but consistent shift for both of these spectroscopic constants. The shift in D e is slightly larger than the ASO-only value of 1.68 kcal mol 1, and quite consistent for all methods. Twoelectron relativistic effects are of the same order of magnitude as in F 2, we find a frequency shift of 0.2 cm 1, a bond expansion of Å and an increase of the D e by 0.07 kcal mol 1 upon adding the perturbative Gaunt correction. C. Bromine r e Å e cm 1 D e kcal mol 1 F 2 HF MP CISD CCSD CCSD T Cl 2 HF MP CISD CCSD CCSD T Br 2 HF MP CISD CCSD CCSD T I 2 HF MP CISD CCSD CCSD T At 2 HF MP CISD CCSD CCSD T In bromine a slight bond contraction is found. This is not an indication of a relativistically enhanced bond as the frequency decreases by 3 cm 1. The lowering of the dissociation energy, with values ranging from 6.6 to 7.7 kcal mol 1, is still close to the ASO-only value of 7.0 kcal mol 1. Adding the Gaunt interaction again gives only minor changes; a bond length increase of Å, a frequency shift of 0.1 cm 1, and an increase of the D e by 0.12 kcal mol 1. The results may be compared with the results of Lee and Lee, 19 who did RECP-MP2 calculations with and without spin orbit coupling. Using Christiansen and Ermler s RECP Ref. 20 they find a bond length of 2.31 Å which is increased upon inclusion of the spin orbit effect by 0.01 Å to 2.32 Å. The effect of the spin orbit operator on the frequency is 6 cm 1, twice as large as the total relativistic effect calculated by us. Schwerdtfeger et al. 21 used different pseudopotentials and did HF and CI calculations. At the nonrelativistic HF level they find an r e of Å that decreases to Å upon addition of scalar relativistic effects and to Å upon addition of both scalar and spin orbit effects. The total
5 9044 L. Visscher and K. G. Dyall: Relativistic effects on dihalogens. I contraction of Å agrees well with our HF/DF results. At the CI level the only semiempirical pseudopotential that they used includes core-polarization effects. This gives an r e of Å, which is in between our double and triplet zeta CISD results. D. Iodine As the relativistic effects increase on going to iodine the dependence on the quality of the basis set becomes relatively less important. For iodine both the DZP and the TZP basis give the same shifts in frequency and a similar shift in the bond length. The method used to calculate the relativistic effect is more important the relativistic effect on the r e value ranges from Å at the HF level to Å at the CCSD T level of theory. Relativity leads to significant weakening of the bond, causing longer bond lengths and smaller dissociation energies than those calculated at the nonrelativistic level of theory. The calculated relativistic effect on the D e ranges from 14.7 to 11.9 kcal mol 1, whereas 14.5 kcal mol 1 is the atomic effect of 2/3 of the SO splitting. This indicates that the quenching of the spinorbit coupling in the molecule is still fairly complete. Twoelectron relativistic effects may be neglected since the absolute size of the two-electron Gaunt correction to the properties remains almost the same as for the previous molecules. We find a bond length increase of Å, a frequency decrease of 0.04 cm 1, and a D e increase by 0.12 kcal mol 1. We can again compare with several RECP-MP2 results. Lee and Lee 19 use Christiansens and Ermlers RECP Ref. 22 and find a spin orbit induced bond length change of 0.02 Å. The RECP-MP2 value is 2.73 Å. The RECP-MP2 spin orbit induced change in frequency is 11 cm 1, slightly larger than the all electron DC-MP2 total relativistic value of 9 cm 1. Schwerdtfeger s results are also available for this molecule. They find a relativistic expansion of Å at the HF level, slightly larger than our calculated value. Their CISD r e value is Å, which is slightly shorter than our CISD values. This can be due to the increasing importance of corevalence correlation that is neglected in our approach, since we do not correlate the d-electrons. When the d-electrons are included in the valence space a contraction of 0.02 Å both at the MP2/pVTZ level and CCSD/pVTZ of theory is found. 23 An extensive treatment of the ground and excited states of I 2 is given by Teichteil and Pelissier. 24 They give results for the ground state using a relativistic averaged pseudopotential and an effective one-electron spin orbit operator with a valence basis of triple zeta quality plus 2s, 2p, and 2d polarization functions. At the MCSCF/MRSDCI Q level of theory they find a bond distance of 2.75 Å without the SOoperator and 2.77 with the SO-operator. The e values are 199 and 185 cm 1, respectively. Again these changes are larger than the total relativistic effect found by us which may either point to an overestimation of the spin orbit effect by the effective one-electron spin orbit operator or to a cancellation of scalar relativistic and spin orbit effects. Teichteil and Pelissier also studied the importance of core-valence polarization effects by comparing 7-electron pseudopotential results with 17-electron pseudopotential results and found a small effect Å at the MCSCF/ MRSDCI Q level of theory, but a 0.02 Å contraction when using second-order perturbation theory after the MCSCF step. The latter result is consistent with the finding of Schwerdtfeger et al. who report a contraction of 0.02 Å due to core-valence polarization and correlation effects. Recent work by Dolg, 25 who compared Douglas Kroll MP2 results with and without correlating the d-electrons also gives a contraction of 0.02 Å. This makes it likely that the core-valence correlation is at least partly responsible for the discrepancy between our CCSD T results and experiment. Teichteil and Pelissier estimated the effect of adding f and g functions to the valence basis and found the f function to decrease the bond length by 0.06 Å and the g function by an additional 0.02 Å. The effect of adding an f function agrees well with the difference of 0.05 Å between our double and triple zeta basis set results, which also differ by an f function, in addition to the extra s, p, and d function. E. Astatine Finally we consider the heaviest halogen treated here, astatine. The bond length expansion is 0.1 Å and the molecule is unbound at the Dirac Fock level of theory. The different methods again give somewhat different relativistic shifts for the r e values but show a rather consistent relativistic shift for the e and D e. As we saw for iodine, the relativistic effect is largest at the CCSD T level of theory. We see that the spin orbit effect is no longer completely quenched as the expected relativistic effect on the D e,45 kcal mol 1, is significantly larger than the calculated values. Again two-electron relativistic effects are very small, giving a bond length increase of Å, a frequency decrease of 0.15 cm 1, and a D e increase by 0.14 kcal mol 1. No literature values for the spectroscopic constants could be found. Though the diatomic form is found to be the ground state of neutral At, 21 all isotopes are radioactive and experimental information on the spectroscopic constants seems hard to obtain. Theoretical data is also very scarce and limited to semi-empirical extrapolation of values for other elements. Dirac Fock 26 and pseudopotential calculations 27 are performed for HAt but not for At 2. On the basis of our CCSD and CCSD T results we may predict the properties of At 2. The calculated bond length is consistently too long in our CCSD T /pvtz approach. Using the correction for I 2 and assuming a conservative error margin of 0.1 Å, we predict a bond length of Å. The calculated frequency is consistently too low at this level of theory for Cl 2,Br 2, and I 2, so we again correct with the error found in I 2 and obtain an e value of cm 1. For the D e the error is 7.5, 6.7, and 6.3 kcal mol 1 for Cl 2,Br 2, and I 2, respectively. Correction with the I 2 error of 6.3 kcal mol 1 and assuming an error bar of twice the range of the errors for the D e in these molecules gives an estimated value of 20 3 kcal mol 1.
6 L. Visscher and K. G. Dyall: Relativistic effects on dihalogens. I 9045 Summarizing the trends, we see that the D e is decreased by relativity. Up to bromine, this effect can be modelled by considering the spin orbit coupling as completely quenched in the molecule and correcting only the atomic asymptote. For iodine and especially astatine this simple model fails and scalar relativistic and molecular spin orbit effects should be taken into account. The influence of spin orbit coupling and other relativistic effects on the bond length and harmonic frequency is somewhat less predictable. In chlorine, iodine, and astatine we find an expansion of the bond length relative to the nonrelativistic result, while for bromine we find a small contraction. The effect on the harmonic frequency is more consistent. In all cases we find a decrease, that is more pronounced when we allow for a higher level of electron correlation. The relativistic correction to the bond length and frequency is fairly independent of the basis set employed. The remaining errors, when comparing our highest level CCSD T /pvtz results with experiment, are probably mostly due to the one-particle basis set incompleteness and the neglect of d-correlation for the heavier elements. Higher order relativistic effects are likely to be small, as can be seen from the magnitude of the Gaunt interaction corrections. The comparison with other theoretical work is not always straightforward due to method and basis set differences and the use of a semiempirical core-valence polarization potential in some schemes. 21 The question whether the relativistic pseudopotential schemes do reproduce the effects of relativity in a reliable way is therefor not always unambiguously answered. Comparison with the 2-component RECP- MP2 scheme of Lee and Lee 19 with our DC-MP2/pVDZ results is well possible, since both the basis set and the method are similar. The results for Br 2 and I 2 are in good agreement for the r e and the e and in reasonable agreement for the D e. The larger discrepancy in the D e is understandable since Lee and Lee used a different method to treat the open-shell atoms than we did. Schwerdtfeger et al. 21 also treated both molecules. Their HF results using a 2-component pseudopotential should be comparable to our pvdz or pvtz results since they employed a (9s,6p,2d) valence basis. For both Br 2 and I 2 the PP-HF results are in between our DC-HF/ pvdz and DC-HF/pVTZ results, with the D e s close to the pvtz result and the r e close to the pvdz result. Teichteil and Pelissier 24 employed a contracted [5s,5p,2d] basis and MCSCF followed by a selected CI method. Because this method is quite different from all the methods that we used a straightforward comparison is not possible. Their r e and e for I 2 are 2.77 Å and 185 cm 1, which is closest to our DC-CCSD T /pvdz result. V. CONCLUSIONS The effect of relativity on spectroscopic properties of the dihalogens has been considered by comparing nonrelativistic with relativistic all-electron methods. The halogen halogen bond is weakened by the increasing relativistic hybridization energy of the p-orbitals which leads to an increase in bond length and a decrease in the frequency and the dissociation energy. The additivity of relativistic and correlation effects that is often implicitly used in calculations on heavy systems has been confirmed for the lighter elements, but for the heaviest halogens the nonadditivity is apparent. The similarity of the relativistic corrections, obtained with different methods to treat electron correlation, implies that the relativistic-correlation cross term is well described at low order in the correlation treatment. Consequently one may estimate the relativistic effect at a certain level of theory by a correction calculated on a lower, more economical, level of theory. For the dihalogens treated here it has been found that such a scheme would lead to a small underestimation of the relativistic effects, as the higher level methods tend to show somewhat larger relativistic corrections. The effect of two-electron relativistic corrections, i.e, the Gaunt interaction, has been found to be unimportant at the level of accuracy reached in the calculated spectroscopic constants. It may be of interest for highly accurate calculations of the ligther elements since the effect on the frequency of F 2 and Cl 2 is of the order of a few tenths of a wave number. One has to be careful when comparing results from different methods, especially when the basis set limit is not reached as is usually the case in compounds like the ones described here. When comparison was possible we found the relativistic pseudopotential calculations in good agreement with our all-electron calculations. ACKNOWLEDGMENTS L.V. was supported by a National Research Council postdoctoral fellowship and a Cray Research Grant. K.G.D. was supported by a prime contract NAS from NASA to Eloret and by contract BPNL A-A3 from Pacific Northwest Laboratory to Eloret. Pacific Northwest Laboratory is a multiprogram laboratory operated by the Batelle Memorial Institute for the U.S. Department of Energy under Contract No. DE-AC06-76RLO-1830; the work was performed under the auspices of the Office of Scientific Computing and under the Office of Health and Environmental Research, which funds the Environmental and Molecular Sciences Laboratory Project, D-384. L.V. wishes to thank Michael Dolg and Bert de Jong for providing him with their preliminary I 2 results. 1 L. L. Foldy and S. A. Wouthuysen, Phys. Rev. 78, M. Douglas and N. M. Kroll, Ann. Phys. NY 82, GAUSSIAN 92/DFT, Revision G.2, M. J. Frisch et al., Gaussian, Inc., Pittsburgh, Pennsylvania, MOLPRO94 is a package of ab initio programs written by H.-J. Werner and P. J. Knowles, with contributions from J. Almlöf, R. D. Amos, M. J. O. Deegan, S. T. Elbert, C. Hampel, W. Meyer, K. Peterson, R. Pitzer, A. J. Stone, and P. R. Taylor. 5 L. Visscher, O. Visser, P. J. C. Aerts, H. Merenga, and W. C. Nieuwpoort, Comp. Phys. Comm. 81, L. Visscher, K. G. Dyall, and T. J. Lee, Int. J. Quantum Chem. Symp. 29, A more extensive description, including the open-shell formalism used and the details of the implementation is to be published elsewhere. 7 T. H. Dunning, Jr., J. Chem. Phys. 90, D. Woon and T. Dunning, Jr., J. Chem. Phys. 98, K. G. Dyall, L. Visscher, and T. Saue unpublished ; K. G. Dyall, P. R. Taylor, K. F,gri, Jr., and H. Partridge, J. Chem. Phys. 95,
7 9046 L. Visscher and K. G. Dyall: Relativistic effects on dihalogens. I 10 K. G. Dyall, I. P. Grant, C. T. Johnson, E. P. Plummer, and F. Parpia, Comp. Phys. Comm. 50, A. Schäfer, C. Huber, and R. Ahlrichs, J. Chem. Phys. 100, K. F,gri, Jr. unpublished. The extended Br basis set used in this work is available upon request. 13 A. Stromberg, O. Gropen, and U. Wahlgren, J. Comput. Chem. 4, K. F,gri, Jr. unpublished. The extended I basis set used in this work is available upon request. 15 K. G. Dyall and K. F,gri, Jr., Theor. Chim. Acta in press. 16 W. C. Ermler, R. B. Ross, and P. A. Christiansen, Adv. Quantum Chem. 19, K. A. Peterson, R. A. Kendall, and T. H. Dunning, J. Chem. Phys. 99, D. E. Woon and T. H. Dunning, J. Chem. Phys. 101, S. Y. Lee and Y. S. Lee, Chem. Phys. Lett. 187, M. M. Hurley, L. F. Pacios, P. A. Christiansen, R. B. Ross, and W. C. Ermler, J. Chem. Phys. 84, P. Schwerdtfeger, L. v. Szentpály, K. Vogel, H. Silberbach, H. Stoll, and H. Preuss, J. Chem. Phys. 84, L. A. LaJohn, P. A. Christiansen, R. B. Ross, T. Atashroo, and W. C. Ermler, J. Chem. Phys. 87, W. A. de Jong private communication. 24 C. Teichteil and M. Pelissier, Chem. Phys. 180, M. Dolg private communication. 26 O. Matsuoka, J. Chem. Phys. 97, M. Dolg, W. Kuchle, H. Stoll, H. Preuss, and P. Schwerdtfeger, Mol. Phys. 74, K. P. Huber and G. Herzberg, Constants of Diatomic Molecules Van Nostrand Reinhold, New York, 1979.
Relativistic and correlation effects on molecular properties. II. The hydrogen halides HF, HCl, HBr, HI, and HAt
Relativistic and correlation effects on molecular properties. II. The hydrogen halides HF, HCl, HBr, HI, and HAt L. Visscher Laboratory of Chemical Physics and Materials Science Center, University of Groningen,
More informationRelativistic and correlated calculations on the ground, excited, and ionized states of iodine
Relativistic and correlated calculations on the ground, excited, and ionized states of iodine W. A. de Jong, L. Visscher, a) and W. C. Nieuwpoort Laboratory for Chemical Physics and Materials Science Centre,
More informationAn ab initio two-component relativistic method including spin orbit coupling using the regular approximation
JOURNAL OF CHEMICAL PHYSICS VOLUME 113, NUMBER 10 8 SEPTEMBER 2000 An ab initio two-component relativistic method including spin orbit coupling using the regular approximation S. Faas, J. H. van Lenthe,
More informationDetlev Figgen, Erich Goll, and Hermann Stoll b) Institut für Theoretische Chemie, Universität Stuttgart, D Stuttgart, Germany
JOURNAL OF CHEMICAL PHYSICS VOLUME 119, NUMBER 21 1 DECEMBER 2003 Systematically convergent basis sets with relativistic pseudopotentials. II. Small-core pseudopotentials and correlation consistent basis
More informationApproximating the basis set dependence of coupled cluster calculations: Evaluation of perturbation theory approximations for stable molecules
JOURNAL OF CHEMICAL PHYSICS VOLUME 113, NUMBER 18 8 NOVEMBER 2000 Approximating the basis set dependence of coupled cluster calculations: Evaluation of perturbation theory approximations for stable molecules
More informationQUANTUM CHEMISTRY FOR TRANSITION METALS
QUANTUM CHEMISTRY FOR TRANSITION METALS Outline I Introduction II Correlation Static correlation effects MC methods DFT III Relativity Generalities From 4 to 1 components Effective core potential Outline
More informationA fully relativistic Dirac Hartree Fock and second-order Mo ller Plesset study of the lanthanide and actinide contraction
JOURNAL OF CHEMICAL PHYSICS VOLUME 109, NUMBER 24 22 DECEMBER 1998 A fully relativistic Dirac Hartree Fock and second-order Mo ller Plesset study of the lanthanide and actinide contraction J. K. Laerdahl
More informationAccurate multireference configuration interaction calculations on the lowest 1 and 3 electronic states of C 2,CN, BN, and BO
Accurate multireference configuration interaction calculations on the lowest 1 and 3 electronic states of C 2,CN, BN, and BO Kirk A. Peterson a) Department of Chemistry, Washington State University and
More informationBenchmark calculations with correlated molecular wave functions
Theor Chem Acc (1997) 97:251±259 Benchmark calculations with correlated molecular wave functions XII. Core correlation e ects on the homonuclear diatomic molecules B 2 -F 2 Kirk A. Peterson 1, Angela K.
More informationApplications of Newly Developed spdsmcps for First-Row Transition Metal Atoms
1st WSEAS Int. Conf. on COMPUTATIONAL CHEMISTRY, Cairo, Egypt, December 29-31, 2007 14 Applications of Newly Developed spdsmcps for First-Row Transition Metal Atoms E. MIYOSHI, 1 Y. OSANAI, 2 M. S. MON,
More informationLec20 Fri 3mar17
564-17 Lec20 Fri 3mar17 [PDF]GAUSSIAN 09W TUTORIAL www.molcalx.com.cn/wp-content/uploads/2015/01/gaussian09w_tutorial.pdf by A Tomberg - Cited by 8 - Related articles GAUSSIAN 09W TUTORIAL. AN INTRODUCTION
More informationTheoretical determination of the heat of formation of methylene
Theoretical determination of the heat of formation of methylene Nikos L. Doltsinis and Peter J. Knowles School of Chemistry, University of Birmingham, Edgbaston, Birmingham B5 2TT, United Kingdom The heat
More informationFormulation and implementation of the relativistic Fock-space coupled cluster method for molecules
JOURNAL OF CHEMICAL PHYSICS VOLUME 115, NUMBER 21 1 DECEMBER 2001 Formulation and implementation of the relativistic Fock-space coupled cluster method for molecules Lucas Visscher a) Department of Theoretical
More informationSystematically convergent basis sets for transition metals. I. All-electron correlation consistent basis sets for the 3d elements Sc Zn
THE JOURNAL OF CHEMICAL PHYSICS 123, 064107 2005 Systematically convergent basis sets for transition metals. I. All-electron correlation consistent basis sets for the 3d elements Sc Zn Nikolai B. Balabanov
More informationPerformance of Hartree-Fock and Correlated Methods
Chemistry 460 Fall 2017 Dr. Jean M. Standard December 4, 2017 Performance of Hartree-Fock and Correlated Methods Hartree-Fock Methods Hartree-Fock methods generally yield optimized geomtries and molecular
More informationSolution of the Electronic Schrödinger Equation. Using Basis Sets to Solve the Electronic Schrödinger Equation with Electron Correlation
Solution of the Electronic Schrödinger Equation Using Basis Sets to Solve the Electronic Schrödinger Equation with Electron Correlation Errors in HF Predictions: Binding Energies D e (kcal/mol) HF Expt
More informationRelativistic double-zeta, triple-zeta, and quadruple-zeta basis sets for the lanthanides La Lu
Theor Chem Acc (2010) 127:369 381 DOI 10.1007/s00214-009-0725-7 REGULAR ARTICLE Relativistic double-zeta, triple-zeta, and quadruple-zeta basis sets for the lanthanides La Lu André S. P. Gomes Kenneth
More informationDensity functional theory predictions of anharmonicity and spectroscopic constants for diatomic molecules
JOURNAL OF CHEMICAL PHYSICS VOLUME 115, NUMBER 6 8 AUGUST 2001 Density functional theory predictions of anharmonicity and spectroscopic constants for diatomic molecules Mutasem Omar Sinnokrot and C. David
More informationJack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah
1. Born-Oppenheimer approx.- energy surfaces 2. Mean-field (Hartree-Fock) theory- orbitals 3. Pros and cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usually does HF-how? 6. Basis sets and notations
More informationTheoretical study of spin-orbit coupling constants for O 2
JOURNAL OF CHEMICAL PHYSICS VOLUME 115, NUMBER 16 22 OCTOBER 2001 Theoretical study of spin-orbit coupling constants for O 2 A 2 3Õ2,1Õ2u, v Ä0 17 and a 4 5Õ2,3Õ2,1Õ2,À1Õ2u, v Ä0 25 D. G. Fedorov, M. S.
More informationHighly accurate quantum-chemical calculations
1 Highly accurate quantum-chemical calculations T. Helgaker Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Oslo, Norway A. C. Hennum and T. Ruden, University
More informationAb initio treatment of electron correlations in polymers: Lithium hydride
JOURNAL OF CHEMICAL PHYSICS VOLUME 112, NUMBER 10 8 MARCH 2000 Ab initio treatment of electron correlations in polymers: Lithium hydride chain and beryllium hydride polymer Ayjamal Abdurahman a) Max-Planck-Institut
More informationAN INTRODUCTION TO QUANTUM CHEMISTRY. Mark S. Gordon Iowa State University
AN INTRODUCTION TO QUANTUM CHEMISTRY Mark S. Gordon Iowa State University 1 OUTLINE Theoretical Background in Quantum Chemistry Overview of GAMESS Program Applications 2 QUANTUM CHEMISTRY In principle,
More informationA coupled cluster study of the spectroscopic properties and electric dipole moment functions of nitrous sulfide
A coupled cluster study of the spectroscopic properties and electric dipole moment functions of nitrous sulfide Youngshang Pak and R. Claude Woods Department of Chemistry, University of Wisconsin Madison,
More informationExercise 1: Structure and dipole moment of a small molecule
Introduction to computational chemistry Exercise 1: Structure and dipole moment of a small molecule Vesa Hänninen 1 Introduction In this exercise the equilibrium structure and the dipole moment of a small
More informationThe Rigorous Calculation of Molecular Properties to Chemical Accuracy. T. Helgaker, Department of Chemistry, University of Oslo, Norway
1 The Rigorous Calculation of Molecular Properties to Chemical Accuracy T. Helgaker, Department of Chemistry, University of Oslo, Norway A. C. Hennum and T. Ruden, University of Oslo, Norway S. Coriani,
More informationMO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij )
MO Calculation for a Diatomic Molecule Introduction The properties of any molecular system can in principle be found by looking at the solutions to the corresponding time independent Schrodinger equation
More informationI. INTRODUCTION JOURNAL OF CHEMICAL PHYSICS VOLUME 118, NUMBER 3 15 JANUARY ; Electronic mail:
JOURNAL OF CHEMICAL PHYSICS VOLUME 118, NUMBER 3 15 JANUARY 2003 Quasirelativistic theory for the magnetic shielding constant. I. Formulation of Douglas Kroll Hess transformation for the magnetic field
More informationJack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah
1. Born-Oppenheimer approx.- energy surfaces 2. Mean-field (Hartree-Fock) theory- orbitals 3. Pros and cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usually does HF-how? 6. Basis sets and notations
More informationThe Accurate Calculation of Molecular Energies and Properties: A Tour of High-Accuracy Quantum-Chemical Methods
1 The Accurate Calculation of Molecular Energies and Properties: A Tour of High-Accuracy Quantum-Chemical Methods T. Helgaker Centre for Theoretical and Computational Chemistry Department of Chemistry,
More informationDFT calculations of NMR indirect spin spin coupling constants
DFT calculations of NMR indirect spin spin coupling constants Dalton program system Program capabilities Density functional theory Kohn Sham theory LDA, GGA and hybrid theories Indirect NMR spin spin coupling
More informationStuart Carter Department of Chemistry, University of Reading, Reading RG6 2AD, United Kingdom
JOURNAL OF CHEMICAL PHYSICS VOLUME 117, NUMBER 4 22 JULY 2002 The ab initio potential energy surface and vibrational-rotational energy levels of X 2 MgOH Jacek Koput a) Department of Chemistry, Adam Mickiewicz
More informationPseudopotential. Meaning and role
Pseudopotential. Meaning and role Jean-Pierre Flament jean-pierre.flament@univ-lille.fr Laboratoire de Physique des Lasers, Atomes et Molécules (PhLAM) Université de Lille-Sciences et technologies MSSC2018
More informationOVERVIEW OF QUANTUM CHEMISTRY METHODS
OVERVIEW OF QUANTUM CHEMISTRY METHODS Outline I Generalities Correlation, basis sets Spin II Wavefunction methods Hartree-Fock Configuration interaction Coupled cluster Perturbative methods III Density
More informationChemistry 334 Part 2: Computational Quantum Chemistry
Chemistry 334 Part 2: Computational Quantum Chemistry 1. Definition Louis Scudiero, Ben Shepler and Kirk Peterson Washington State University January 2006 Computational chemistry is an area of theoretical
More informationConvergence properties of the coupled-cluster method: the accurate calculation of molecular properties for light systems
1 Convergence properties of the coupled-cluster method: the accurate calculation of molecular properties for light systems T. Helgaker Centre for Theoretical and Computational Chemistry, Department of
More informationBeyond the Hartree-Fock Approximation: Configuration Interaction
Beyond the Hartree-Fock Approximation: Configuration Interaction The Hartree-Fock (HF) method uses a single determinant (single electronic configuration) description of the electronic wavefunction. For
More informationarxiv:cond-mat/ v1 10 May 1996
Cohesive energies of cubic III-V semiconductors Beate Paulus, Peter Fulde Max-Planck-Institut für Physik komplexer Systeme, Bayreuther Str. 40, 01187 Dresden, Germany arxiv:cond-mat/9605064v1 10 May 1996
More informationFull configuration interaction potential energy curves for breaking bonds to hydrogen: An assessment of single-reference correlation methods
JOURNAL OF CHEMICAL PHYSICS VOLUME 118, NUMBER 4 22 JANUARY 2003 Full configuration interaction potential energy curves for breaking bonds to hydrogen: An assessment of single-reference correlation methods
More informationCorrelation effects in MgO and CaO: Cohesive energies and lattice constants
PHYSICAL REVIEW B VOLUME 54, NUMBER 19 15 NOVEMBER 1996-I Correlation effects in MgO and CaO: Cohesive energies and lattice constants Klaus Doll and Michael Dolg Max-Planck-Institut für Physik Komplexer
More informationThird-order Douglas Kroll ab initio model potential for actinide elements
JOURNAL OF CHEMICAL PHYSICS VOLUME 117, NUMBER 8 22 AUGUST 2002 Third-order Douglas Kroll ab initio model potential for actinide elements Jozef Paulovič, Takahito Nakajima, and Kimihiko Hirao a) Department
More informationOslo node. Highly accurate calculations benchmarking and extrapolations
Oslo node Highly accurate calculations benchmarking and extrapolations Torgeir Ruden, with A. Halkier, P. Jørgensen, J. Olsen, W. Klopper, J. Gauss, P. Taylor Explicitly correlated methods Pål Dahle, collaboration
More informationEnergy-adjusted pseudopotentials for the rare earth elements
Theor Chim Acta (1989) 75:173-194 O Springer-Ver!ag 1989 Energy-adjusted pseudopotentials for the rare earth elements M. Dolg, H. Stoll, A. Savin, and H. Preuss Institut f/it Theoretische Chemie, Universit/it
More informationT. Helgaker, Department of Chemistry, University of Oslo, Norway. T. Ruden, University of Oslo, Norway. W. Klopper, University of Karlsruhe, Germany
1 The a priori calculation of molecular properties to chemical accuarcy T. Helgaker, Department of Chemistry, University of Oslo, Norway T. Ruden, University of Oslo, Norway W. Klopper, University of Karlsruhe,
More informationBasis sets for electron correlation
Basis sets for electron correlation Trygve Helgaker Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway The 12th Sostrup Summer School Quantum Chemistry
More informationIntroduction to computational chemistry Exercise I: Structure and electronic energy of a small molecule. Vesa Hänninen
Introduction to computational chemistry Exercise I: Structure and electronic energy of a small molecule Vesa Hänninen 1 Introduction In this exercise the equilibrium structure and the electronic energy
More informationThis is a very succinct primer intended as supplementary material for an undergraduate course in physical chemistry.
1 Computational Chemistry (Quantum Chemistry) Primer This is a very succinct primer intended as supplementary material for an undergraduate course in physical chemistry. TABLE OF CONTENTS Methods...1 Basis
More informationPerformance of coupled cluster theory in thermochemical calculations of small halogenated compounds
JOURNAL OF CHEMICAL PHYSICS VOLUME 118, NUMBER 8 22 FEBRUARY 2003 Performance of coupled cluster theory in thermochemical calculations of small halogenated compounds David Feller, a) Kirk A. Peterson,
More informationAb initio study of the BiSe and BiTe electronic spectra: What happens with X 2 X 1 emission in the heavier Bi chalcogenides?
JOURNAL OF CHEMICAL PHYSICS VOLUME 120, NUMBER 16 22 APRIL 2004 Ab initio study of the BiSe and BiTe electronic spectra: What happens with X 2 X 1 emission in the heavier Bi chalcogenides? Rainer M. Lingott,
More informationChemistry Publications
Chemistry Publications Chemistry 2007 Accurate Ab Initio Potential Energy Curve of F2. I. Nonrelativistic Full Valence Configuration Interaction Energies Using the Correlation Energy Extrapolation by Intrinsic
More informationEffect of the self-interaction error for three-electron bonds: On the development of new exchange-correlation functionals
RESEARCH PAPER Effect of the self-interaction error for three-electron bonds: On the development of new exchange-correlation functionals Jürgen Gräfenstein, Elfi Kraka and Dieter Cremer* PCCP www.rsc.org/pccp
More informationShould bromoform absorb at wavelengths longer than 300 nm?
JOURNAL OF CHEMICAL PHYSICS VOLUME 117, NUMBER 13 1 OCTOBER 2002 Should bromoform absorb at wavelengths longer than 300 nm? Kirk A. Peterson a) Department of Chemistry, Washington State University, Pullman,
More informationLec20 Wed 1mar17 update 3mar 10am
564-17 Lec20 Wed 1mar17 update 3mar 10am Figure 15.2 Shows that increasing the diversity of the basis set lowers The HF-SCF energy considerably, but comes nowhere near the exact experimental energy, regardless
More informationTheoretical study of the low-lying excited singlet states of furan
JOURNAL OF CHEMICAL PHYSICS VOLUME 119, NUMBER 2 8 JULY 2003 Theoretical study of the low-lying excited singlet states of furan E. V. Gromov, A. B. Trofimov, and N. M. Vitkovskaya Laboratory of Quantum
More informationAb initio calculations on the ground and low-lying excited states of InI
MOLECULAR PHYSICS, 1OCTOBER 23, VOL. 11, NO. 19, 2963 2968 Ab initio calculations on the ground and low-lying excited states of InI WENLI ZOU, MEIRONG LIN*, XINZHENG YANG and BAOZHENG ZHANG Institute of
More informationAlgebraic Studies for Electronic Structures and Photodissociation of Bromine Molecule
International Journal of Modern Physics and Applications Vol. 1, No. 1, 2015, pp. 1-5 http://www.publicscienceframework.org/journal/ijmpa Algebraic Studies for Electronic Structures and Photodissociation
More informationAccurate theoretical near-equilibrium potential energy and dipole moment surfaces of HgClO and HgBrO
JOURNAL OF CHEMICAL PHYSICS VOLUME 120, NUMBER 14 8 APRIL 2004 Accurate theoretical near-equilibrium potential energy and dipole moment surfaces of HgClO and HgBrO Nikolai B. Balabanov a) and Kirk A. Peterson
More informationDivergence in Møller Plesset theory: A simple explanation based on a two-state model
JOURNAL OF CHEMICAL PHYSICS VOLUME 112, NUMBER 22 8 JUNE 2000 Divergence in Møller Plesset theory: A simple explanation based on a two-state model Jeppe Olsen and Poul Jørgensen a) Department of Chemistry,
More informationComputational Material Science Part II. Ito Chao ( ) Institute of Chemistry Academia Sinica
Computational Material Science Part II Ito Chao ( ) Institute of Chemistry Academia Sinica Ab Initio Implementations of Hartree-Fock Molecular Orbital Theory Fundamental assumption of HF theory: each electron
More informationPotential energy curves for neutral and multiply charged carbon monoxide
PRAMANA c Indian Academy of Sciences Vol. 74, No. 1 journal of January 2010 physics pp. 49 55 Potential energy curves for neutral and multiply charged carbon monoxide PRADEEP KUMAR 1 and N SATHYAMURTHY
More informationSystematic ab initio calculations on the energetics and stability of covalent O 4
JOURNAL OF CHEMICAL PHYSICS VOLUME 120, NUMBER 21 1 JUNE 2004 Systematic calculations on the energetics and stability of covalent O 4 Ramón Hernández-Lamoneda a) Centro de Investigación en Química, Universidad
More informationChemistry 4560/5560 Molecular Modeling Fall 2014
Final Exam Name:. User s guide: 1. Read questions carefully and make sure you understand them before answering (if not, ask). 2. Answer only the question that is asked, not a different question. 3. Unless
More informationCorrelated ab initio calculations for ground-state properties of II-VI semiconductors
PHYSICAL REVIEW B VOLUME 56, NUMBER 12 15 SEPTEMBER 1997-II Correlated ab initio calculations for ground-state properties of II-VI semiconductors Martin Albrecht and Beate Paulus Max-Planck-Institut für
More informationLecture 10: Relativistic effects. Einstein s great idea as applied to chemistry: why gold is yellow and lead batteries function
Lecture 10: Relativistic effects Einstein s great idea as applied to chemistry: why gold is yellow and lead batteries function Course exam Examination on Monday May 7, 2.00 pm 4.00 pm Content The whole
More informationRadiative Transition Probabilities and Lifetimes for the Band Systems A 2 Π X 2 Σ + of the Isovalent Molecules BeF, MgF and CaF
950 Brazilian Journal of Physics, vol. 35, no. 4A, December, 2005 Radiative Transition Probabilities and Lifetimes for the Band Systems of the Isovalent Molecules BeF, MgF and CaF Marina Pelegrini a, Ciro
More informationComputational Chemistry. An Introduction to Molecular Dynamic Simulations
Computational Chemistry An Introduction to Molecular Dynamic Simulations Computational chemistry simulates chemical structures and reactions numerically, based in full or in part on the fundamental laws
More informationLecture 4: methods and terminology, part II
So theory guys have got it made in rooms free of pollution. Instead of problems with the reflux, they have only solutions... In other words, experimentalists will likely die of cancer From working hard,
More informationAb initio calculations for potential energy surfaces. D. Talbi GRAAL- Montpellier
Ab initio calculations for potential energy surfaces D. Talbi GRAAL- Montpellier A theoretical study of a reaction is a two step process I-Electronic calculations : techniques of quantum chemistry potential
More informationThe calculation of the universal density functional by Lieb maximization
The calculation of the universal density functional by Lieb maximization Trygve Helgaker, Andy Teale, and Sonia Coriani Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry,
More informationOrbital entanglement in bond-formation processes
Orbital entanglement in bond-formation processes Katharina Boguslawski a, Pawel Tecmer a, Gergely Barcza b, Örs Legezab,1, and Markus Reiher a,2 arxiv:10.20v2 [physics.chem-ph] 11 May 201 a ETH Zurich,
More informationThe ab initio model potential method: Third-series transition metal elements
JOURNAL OF CHEMCAL PHYSCS VOLUME 110, NUMBER 2 8 JANUARY 1999 The ab initio model potential method: Third-series transition metal elements Marcos Casarrubios and Luis Seijo a) Departamento de Química,
More informationComputational Methods. Chem 561
Computational Methods Chem 561 Lecture Outline 1. Ab initio methods a) HF SCF b) Post-HF methods 2. Density Functional Theory 3. Semiempirical methods 4. Molecular Mechanics Computational Chemistry " Computational
More informationElectron Correlation Methods
Electron Correlation Methods HF method: electron-electron interaction is replaced by an average interaction E HF c = E 0 E HF E 0 exact ground state energy E HF HF energy for a given basis set HF E c
More informationBasis Set for Molecular Orbital Theory
Basis Set for Molecular Orbital Theory! Different Types of Basis Functions! Different Types of Atom Center Basis Functions! Classifications of Gaussian Basis Sets! Pseudopotentials! Molecular Properties
More informationG3-RAD and G3X-RAD: Modified Gaussian-3 (G3) and Gaussian-3X (G3X) procedures for radical thermochemistry
G3-RAD and G3X-RAD: Modified Gaussian-3 (G3) and Gaussian-3X (G3X) procedures for radical thermochemistry David J. Henry, Michael B. Sullivan, and Leo Radom Citation: The Journal of Chemical Physics 118,
More informationDensity Functional Theory. Martin Lüders Daresbury Laboratory
Density Functional Theory Martin Lüders Daresbury Laboratory Ab initio Calculations Hamiltonian: (without external fields, non-relativistic) impossible to solve exactly!! Electrons Nuclei Electron-Nuclei
More information4 Diatomic molecules
s manual for Burrows et.al. Chemistry 3 Third edition 4 Diatomic molecules Answers to worked examples WE 4.1 The Lewis model (on p. 174 in Chemistry 3 ) Use the Lewis model to describe the bonding in (a)
More informationAn ab initio study of some noble gas monohalides
JOURNAL OF CHEMICAL PHYSICS VOLUME 114, NUMBER 5 1 FEBRUARY 2001 An ab initio study of some noble gas monohalides Gerald J. Hoffman a) and Mitchell Colletto Department of Chemistry, The College of New
More information( R)Ψ el ( r;r) = E el ( R)Ψ el ( r;r)
Born Oppenheimer Approximation: Ĥ el ( R)Ψ el ( r;r) = E el ( R)Ψ el ( r;r) For a molecule with N electrons and M nuclei: Ĥ el What is E el (R)? s* potential surface Reaction Barrier Unstable intermediate
More informationCOUPLED-CLUSTER CALCULATIONS OF GROUND AND EXCITED STATES OF NUCLEI
COUPLED-CLUSTER CALCULATIONS OF GROUND AND EXCITED STATES OF NUCLEI Marta Włoch, a Jeffrey R. Gour, a and Piotr Piecuch a,b a Department of Chemistry,Michigan State University, East Lansing, MI 48824 b
More informationChemistry 3502/4502. Final Exam Part I. May 14, 2005
Advocacy chit Chemistry 350/450 Final Exam Part I May 4, 005. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle
More informationAtom-molecule molecule collisions in spin-polarized polarized alkalis: potential energy surfaces and quantum dynamics
Atom-molecule molecule collisions in spin-polarized polarized alkalis: potential energy surfaces and quantum dynamics Pavel Soldán, Marko T. Cvitaš and Jeremy M. Hutson University of Durham with Jean-Michel
More informationInvestigation of Spectroscopic Properties and Spin-Orbit Splitting in the X 2 Π and A 2 Π Electronic States of the SO + Cation
Int. J. Mol. Sci. 2012, 13, 8189-8209; doi:10.3390/ijms13078189 Article OPEN ACCESS International Journal of Molecular Sciences ISSN 1422-0067 www.mdpi.com/journal/ijms Investigation of Spectroscopic Properties
More informationGaussian Basis Sets for Solid-State Calculations
Gaussian Basis Sets for Solid-State Calculations K. Doll Molpro Quantum Chemistry Software Institute of Theoretical Chemistry, D-70569 Stuttgart, Germany MW-MSSC 2017, Minneapolis, July 10, 2017 Introduction
More informationSystematic study of the static electrical properties of the CO molecule: Influence of the basis set size and correlation energy
Systematic study of the static electrical properties of the CO molecule: Influence of the basis set size and correlation energy Josep M. Luis, Josep Martí, Miquel Duran, and José L. Andrés a) Institute
More informationCOMPUTATIONAL STUDIES OF BONDING AND PHOSPHORESCENT PROPERTIES OF GROUP 12 OLIGOMERS AND EXTENDED EXCIMERS. John J. Determan, B.A.
COMPUTATIONAL STUDIES OF BONDING AND PHOSPHORESCENT PROPERTIES OF GROUP 12 OLIGOMERS AND EXTENDED EXCIMERS John J. Determan, B.A. Thesis Prepared for the Degree of MASTER OF SCIENCE UNIVERSITY OF NORTH
More information154 J. Chem. Phys. 108 (1), 1 January /98/108(1)/154/23/$ American Institute of Physics
An examination of intrinsic errors in electronic structure methods using the Environmental Molecular Sciences Laboratory computational results database and the Gaussian-2 set David Feller and Kirk A. Peterson
More informationHomologation of Boronic Esters with Organolithium Compounds: A Computational Assessment of Mechanism
Homologation of Boronic Esters with Organolithium Compounds: A Computational Assessment of Mechanism Stéphanie Essafi,*,1 Simone Tomasi, 2 Varinder K. Aggarwal, 1 Jeremy Harvey*,1 1 School of Chemistry,
More informationAb initio study of spectroscopic and radiative characteristics of ion-pair states of the Cl 2 molecule
JOURNAL OF CHEMICAL PHYSICS VOLUME 115, NUMBER 20 22 NOVEMBER 2001 Ab initio study of spectroscopic and radiative characteristics of ion-pair states of the Cl 2 molecule D. B. Kokh, a) A. B. Alekseyev,
More informationBasis Sets and Basis Set Notation
Chemistry 46 Fall 215 Dr. Jean M. Standard November 29, 217 Basis Sets and Basis Set Notation Using the LCAO-MO approximation, molecular orbitals can be represented as linear combinations of atomic orbitals,
More informationCHEM1901/ J-5 June 2013
CHEM1901/3 2013-J-5 June 2013 Oxygen exists in the troposphere as a diatomic molecule. 4 (a) Using arrows to indicate relative electron spin, fill the left-most valence orbital energy diagram for O 2,
More informationLowest vibrational states of 4 He 3 He + : Non-Born-Oppenheimer calculations
Lowest vibrational states of He 3 He + : Non-Born-Oppenheimer calculations Monika Stanke,, * Dariusz Kȩdziera, Sergiy Bubin, Marcin Molski, 3 and Ludwik Adamowicz, Department of Chemistry, University of
More informationQuantum Chemistry Methods
1 Quantum Chemistry Methods T. Helgaker, Department of Chemistry, University of Oslo, Norway The electronic Schrödinger equation Hartree Fock theory self-consistent field theory basis functions and basis
More informationVol. 9 COMPUTATIONAL CHEMISTRY 319
Vol. 9 COMPUTATIONAL CHEMISTRY 319 COMPUTATIONAL QUANTUM CHEMISTRY FOR FREE-RADICAL POLYMERIZATION Introduction Chemistry is traditionally thought of as an experimental science, but recent rapid and continuing
More informationElectric properties of molecules
Electric properties of molecules For a molecule in a uniform electric fielde the Hamiltonian has the form: Ĥ(E) = Ĥ + E ˆµ x where we assume that the field is directed along the x axis and ˆµ x is the
More informationAn Accurate Calculation of Potential Energy Curves and Transition Dipole Moment for Low-Lying Electronic States of CO
Commun. Theor. Phys. 59 (2013) 193 198 Vol. 59, No. 2, February 15, 2013 An Accurate Calculation of Potential Energy Curves and Transition Dipole Moment for Low-Lying Electronic States of CO LU Peng-Fei
More informationIntroduction to Electronic Structure Theory
CSC/PRACE Spring School in Computational Chemistry 2017 Introduction to Electronic Structure Theory Mikael Johansson http://www.iki.fi/~mpjohans Objective: To get familiarised with the, subjectively chosen,
More informationCorrelation consistent basis sets for lanthanides. The atoms La Lu
Correlation consistent basis sets for lanthanides. The atoms La Lu Qing Lu and Kirk A. Peterson a) Department of Chemistry, Washington State University, Pullman, WA 99164-4630, USA Abstract Using the 3rd-order
More informationBreit interaction in heavy atoms
LETTER TO THE EDITOR Breit interaction in heavy atoms M G Kozlov, SGPorsev, and I I Tupitsyn Petersburg Nuclear Physics Institute, 188350, Gatchina, Russia E-mail: mgk@mf1309.spb.edu St. Petersburg State
More informationHighly Accurate Ab Initio Computation of Thermochemical Data
Chapter 1 Highly Accurate Ab Initio Computation of Thermochemical Data Trygve Helgaker Department of Chemistry, University of Oslo, P. O. Box 1033 Blindern, N-0315 Oslo, Norway Wim Klopper and Asger Halkier
More information