Frequency-dependent polarizability of open-shell atomic systems

The time-dependent variational principle is used to calculate the frequency-dependent dipole polarizabilities of 2p open-shell atomic systems. A general theory is built up to include frequency-dependent perturbations using the Roothaan–Hartree–Fock (HFR) formalism. The excitation energies corresponding to the low-lying excited states of the system are obtained from the poles of the dynamic polarizability values. A fairly accurate knowledge of transition oscillator strengths is also obtained from a knowledge of the polarizability values near the poles. The excitation energies and oscillator strengths are compared with available data.     

Perturbation treatment of the Hartree–Fock equations of 1s2p3s 4Pθ state of the lithium isoelectronic sequence

The first order Hartree–Fock equations of the 1s2p3s ⁴P⁰ state of the three-electron atomic systems have been solved exactly. These solutions are used to evaluate Hartree–Fock energy up to third order with high accuracy. The third order Hartree–Fock energies for Li to Ne⁷⁺ are compared with those derived from experiment and other theoretical calculations.     

Dispersion interaction between helium pair from hydrodynamic analogy to quantum mechanics

Frequency-dependent multipole polarizabilities of the He sequence have been calculated from a hydrodynamic model of quantum mechanics and using an independent-particle model Hamiltonian. Our present scheme is parallel to the uncoupled Hartree–Fock approximation so our values of polarizabilities, multipole transition energies and dispersion force coefficients between the He? He pair are comparable with earlier works using the uncoupled Hartree–Fock approximation.     

Electron correlation and relativistic effects in xenon tetrafluoride

Ab initio all-electron fully relativistic Dirac–Fock self-consistent field and Dirac–Fock–Breit calculations are reported for the XeF₄ molecule at various internuclear distances assuming the experimental D_4h geometry with our recently developed relativistic universal Gaussian basis set. The nonrelativistic limit Hartree–Fock calculations were also performed for XeF₄ at various internuclear distances. The calculated relativistic correction to the total energy of molecule at the Dirac–Fock level is ～ ?5856 eV, whereas the magnetic part of the Breit correction to the electron-electron interaction is calculated as ～ 177 eV. The electron correlation effects were included in the nonrelativistic Hartree–Fock calculations using the second-order Møller-Plesset (MP 2) theory, and the calculated correlation energy for XeF₄ is ?71 eV. The basis-set superposition error (BSSE ) was estimated by using the counterpoise method for Xe and F. The inclusion of both the relativistic and electron correlation effects in the calculated total energies of F, Xe, and XeF₄ predicts the Xe—F bond length and dissociation energy of XeF₄ as 1.952 Å and 5.59 eV, respectively, which are in excellent agreement with the experimental values of 1.953 Å and 5.69 eV, respectively, for XeF₄. The contribution of the electron correlation and relativistic effects to the dissociation energy of XeF₄ is 8.11 and 0.05 eV, respectively. The Breit interaction, however, contributes only 0.02 eV to the dissociation energy of XeF₄. Electron correlation is most significant for the prediction of an accurate value of dissociation energy, whereas relativistic effects are very important for the prediction of spin-orbital splitting as well as the energies of the orbitals, especially the inner orbitals of XeF₄. © 1995 John Wiley & Sons, Inc.     

Ground-state energies of closed-shell atoms

Ground-state energies of He, Be, and Ne isoelectronic series have been calculated. The Dirac–Hartree–Fock energy values have been corrected by adding Breit, vacuum polarization, self-energy, nuclear mass, and electron correlation corrections. The resulting energies are compared with the experimental values. As a result, an estimate of the correlation-relativistic cross-term energy is obtained. The effect, for large-Z atoms, proves to be quite substantial.     

Kβ/Kα X-ray intensity ratios for elements in the range 16⩽Z⩽92 excited by 5.9, 59.5 and 123.6 keV photons

The K shell intensity ratios K_β/K_α for 59 elements in the atomic region 16⩽Z⩽92 have been measured at excitation energies of 5.9, 59.5 and 123.6 keV. K X-rays emitted by samples have been counted by a Si(Li) detector with resolution 160 eV at 5.9 keV. The measured values were compared with the theoretical values calculated using Scofield's tables based on the Hartree–Slater and Hartree–Fock theories and available experimental values. Reasonable agreement is typically obtained between present and theoretical values.     

Excited state properties from the equations of motion method. Application of the MCTDHF-MCRPA to the dipole moments and oscillator strengths of the A 1Π, a′3Π, a′3Σ+ and d3Δ low-lying valence states of CO

By examining the exact operator O_λ⁺ which is the solution of the equations of motion-Green's function method, we rederive expressions for non-reference (usually excited) state properties. Hence, additional useful information such as state expectation values, oscillator strengths, and frequency dependent and independent polarizabilities may be easily obtained from an equation of motion-Green's function calculation. With the multiconfigurational random phase approximation (MCRPA), which is equivalent to the multiconfigurational time dependent Hartree-Fock (MCTDHF), excitation energies, oscillator strengths, and excitation operators from the ground states are obtained for the low-lying valence (under 10 eV above the ground state) states of CO at the experimental ground state equilibrium geometry. We apply these techniques to obtain the excited state dipole moments for and oscillator strengths between the A ¹Π, a ³Π, a′ ³Σ⁺, and d ³Δ states of CO and compare our results to other calculations and experiments.     

Ab initio study of organic mixed valency

A series of six radical cations of the type (D L D)⁺ was investigated at the ab initio unrestricted Hartree–Fock level. One localized and one delocalized conformation were systematically searched by full geometry optimization. At both nuclear arrangements, mostly found as being minima in the symmetry‐restrained Hartree–Fock framework, excitation energies were calculated through the expansion of the wave function on single electronic excitations of the Hartree–Fock fundamental determinant and at the unrestricted Hartree–Fock or at the multiconfigurational self consistent field levels. Few calculations were also performed by taking into account some part of the electronic correlation. Except for N,N,N′,N′‐tetramethyl p‐phenylenediamine, all the studied compounds are localized stable cations, at the symmetry‐restrained Hartree–Fock level. However, the reoptimization of their wave function changes this observation since only three of them seem to conserve a localized stable conformation. Most of the studied systems are characterized by one or two excited electronic states very close to the fundamental one and should thus present an unresolved broadened first absorption band in the near‐infrared region. These features are in agreement with the available experimental data. Strong Hartree–Fock instabilities are found for the delocalized structure and put in relation with the existence of the large nonadiabatic coupling in this conformational region. The solvent influence is discussed in the Onsager dipolar reaction field framework. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 76: 552–573, 2000     

New algorithm for high-order time-dependent hartree–fock theory for nonlinear optical properties

A new formalism and algorithm is developed for solving the general-order time-dependent Hartree–Fock (TDHF) problem. It is shown that for any order a generalization of the TDHF equations can be derived where all lower-order solutions constitute a constant term. This makes it very easy to obtain high-order solutions. As the space required for the mapping of density matrices to Fock matrices in a problem of a giver order is largely reduced, we can perform the most time-consuming steps within the core memory of the machine and easily manipulate vector products via optimum routines. The second hyperpolarizability γ is obtained from the secondorder TDHF solution via a 2n rule. The formalism also allows for expressing all terms in the equation diagrammatically, which provides additional physical insight and a more systematic evaluation of terms. To illustrate the method, TDHF results are presented for trans-butadiene and carbon monoxide for several optical processes, including correlation corrections to their static hyperpolarizabilities obtained via coupled cluster (CCSD) and many-body perturbation theory. The hybrid TDHF/CCSD method provides excellent agreement with the DC–SHG experiments (χ||⁽²⁾ = 11.4 x 10^?32 esu/mol compared to 12.9 ± 11.4 × 10^?32 esu/mol and χ||⁽³⁾ = 149 compared to 144 ± 4 × 10^?39 esu/mol).     

Extended average energy perturbation calculation of the 2s2 and 2p21S resonance state energies of two-electron atoms

An extended average energy (EAE) calculation of the 2s² and 2p²¹S resonance state energies of two-electron atoms is carried out. We take the bare-nucleus hamiltonian as the initial approximation and treat the zeroth-order degeneracy by van Vleck perturbation theory, For both H^? and He the lower state energies agree quite well with experiment while the upper state is about 1–2 eV too high. Our remits for He are comparable in accuracy to those obtained by the 1/Z Hartree—Fock perturbation method. A brief concluding discussion of ways to improve the simple EAE technique is presented.     

Self-consistent-field method including correlation effect for atomic systems. Three- and four-electron atomic systems

The correlation effect for three- and foru-electron atomic systems has been taken into account by modifying the potentials of the electron interactions appearing in the Hartree–Fock equations. The correlation energies obtained for Li, Be⁺, B²⁺, Li^?, Be and B⁺ differ by less than 25 percent from the exact values.     

Atomic Xα calculations based on EHFS(α) = Eexp

The total energies and one-electron energies for first- and second-row atoms were calculated by using the Hartree–Fock and the Hartree–Fock-Slater Hamiltonian with X_α orbitals, u_i(α_exp); α was parametrized from E^HFS (α_exp) = E^exp. The E^HF (α_exp) total energies are always higher than the Hartree–Fock energies for the atoms. The relation of the calculated ionization potential to the experimental ionization potential depends on the α used to define u_i(α), α_exp, or α_HF.     

Ab initio calculations of relativistic and electron correlation effects in polyatomics using the universal Gaussian basis set: XeF2

Ab initio accurate all-electron relativistic molecular orbital Dirac–Fock self-consistent field calculations are reported for the linear symmetric XeF₂ molecule at various internuclear distances with our recently developed relativistic universal Gaussian basis set. The nonrelativistic limit Hartree–Fock calculations were also performed for XeF₂ at various internuclear distances. The relativistic correction to the electronic energy of XeF₂ was calculated as ～ ?215 hartrees (?5850 eV) by using the Dirac–Fock method. The dominant magnetic part of the Breit interaction correction to the nonrelativistic interelectron Coulomb repulsion was included in our calculations by both the Dirac–Fock–Breit self-consistent field and perturbation methods. The calculated Breit correction is ～6.5 hartrees (177 eV) for XeF₂. The relativistic Dirac–Fock as well as the nonrelativistic HF wave functions predict XeF₂ to be unbound, due to neglect of electron correlation effects. These effects were incorporated for XeF₂ by using various ab initio post Hartree–Fock methods. The calculated dissociation energy obtained using the MP 2(full) method with our extensive basis set of 313 primitive Gaussians that included d and f polarization functions on Xe and F is 2.77 eV, whereas the experimental dissociation energy is 2.78 eV. The calculated correlation energy is ～ ?2 hartrees (?54 eV) at the predicted internuclear distance of 1.986 Å, which is in excellent agreement with the experimental Xe—F distance of 1.979 Å in XeF₂. In summary, electron correlation effects must be included in accurate ab initio calculations since it has been shown here that their inclusion is crucial for obtaining theoretical dissociation energy (D_e) close to experimental value for XeF₂. Furthermore, relativistic effects have been shown to make an extremely significant contribution to the total energy and orbital binding energies of XeF₂. © 1995 John Wiley & Sons, Inc.     

Ab initio studies of the electronic structure and density of states of metallic beryllium

Hartree–Fock–Roothaan studies are reported for low-lying electronic states of metallic beryllium as modeled by a moiety of 135 beryllium atoms. The system corresponds to 16 coordination shells of a central Be with internuclear separations derived from the lattice constants of the bulk metal. The calculations become tractable by use of the full D_3h symmetry of the system at both the integrals and self-consistent-field stages and by employing ab initio effective potentials for the 1s electrons of each beryllium atom. Ionization potentials, binding energies, orbital energies, electric field gradients, nuclear-electrostatic potentials, diamagnetic shielding constants, second moments, and Mulliken populations are calculated for selected electronic states. The calculated ionization potential for the lowest state agrees to within 10% of the experimental bulk work function. A density-of-states analysis for that state is reported and compared with band structure calculations.     

One-Parameter double-zeta atomic functions for the hartree–fock energies of helium isoelectronic sequence

The double-zeta atomic functions are characterized by the nuclear charge z of the two-electron atomic system. The Hartree–Fock total energies and the corresponding orbital energies are calculated using various atomic wave functions for the helium isoelectronic sequence. The expectation values rⁿ of various wave functions are also examined. It is found that the accuracy of our one-parameter double-zeta functions corresponds to the accuracy of the usual five-parameter double-zeta functions.     

Ab initio calculation of the effective valence shell hamiltonian of carbon: Simultaneous treatment of neutral and ion states

The n = 2 effective valence shell hamiltonian, $\text{H}$^v, of carbon is evaluated through second order using ³P Hartree—Fock orbitals (5s4p) with added d functions to provide results within a few percent of the spd convergence limits. The calculated $\text{H}$^v is employed to evaluate the n = 2 valence states of C, C^?, C⁺, C²⁺ and C³⁺ with an average deviation of the 21 excitation energies, ionization potentials and electron affinity from experimental values of 0.32 eV. Three-electron parts of $\text{H}$^v contribute substantially to a number of these excitation energies.     

Geometry,basis set,and correlation energy dependence of computed protonation energies of imino bases

The nitrogen protonation energies of the imino bases HN?CHR, where R is H, CH₃, NH₂, OH, and F, have been evaluated to determine the dependence of absolute and relative protonation energies on geometry, basis set, and correlation effects. Reliable absolute protonation energies require a basis set larger than a split-valence plus polarization basis, the inclusion of correlation, and optimized geometries of at least Hartree–Fock 4-31G quality. Consistent relative protonation energies can be obtained at the Hartree–Fock level with smaller basis sets. Extending the split-valence basis set by the addition of polarization functions on all atoms decreases the computed absolute Hartree–Fock nitrogen protonation energies of the imino bases HN?CHR except when R is F, but increases the oxygen protonation energies of the carbonyl bases O?CHR.     

Comprehensive understanding of size‐, shape‐, and composition‐dependent polarizabilities of SimCn (m,n = 1–4) clusters

We performed a comprehensive study of the size‐, shape‐, and composition‐dependent polarizabilities of Si_mC_n (m, n = 1–4) clusters on the basis of the density‐functional‐based coupled perturbed Hartree–Fock calculations. We found better correlations between the polarizabilities and both the binding energies (E_b) and change in charge distribution (Δq) than the energy gaps. The α values exhibit overall decreasing and increasing trends with increases in the E_b and Δq values, respectively. For isomers with the same E_b values and different polarizabilities, Δq can well explain the difference in polarizabilities. The π‐electron delocalization effect is the best factor for understanding the shape‐dependence. For a given m/n value, the linear clusters have an obviously larger polarizability than both the prolate and compact clusters, irrespective of the cluster size. We fit a quantitative expression [α = A ? (A ? B) × exp(?k(m/n))] to describe the composition‐dependent polarizabilities. © 2012 Wiley Periodicals, Inc.     

An improved generator coordinate Hartree–Fock method applied to the choice of contracted Gaussian basis sets for first‐row diatomic molecules

Accurate Gaussian basis sets (18s for Li and Be and 20s11p for the atoms from B to Ne) for the first‐row atoms, generated with an improved generator coordinate Hartree–Fock method, were contracted and enriched with polarization functions. These basis sets were tested for B₂, C₂, BeO, CN⁻, LiF, N₂, CO, BF, NO⁺, O₂, and F₂. At the Hartree–Fock (HP), second‐order Møller–Plesset (MP2), fourth‐order Møller–Plesset (MP4), and density functional theory (DFT) levels, the dipole moments, bond lengths, and harmonic vibrational frequencies were studied, and at the MP2, MP4, and DFT levels, the dissociation energies were evaluated and compared with the corresponding experimental values and with values obtained using other contracted Gaussian basis sets and numerical HF calculations. For all diatomic molecules studied, the differences between our total energies, obtained with the largest contracted basis set [6s5p3d1f], and those calculated with the numerical HF methods were always less than 3.2 mhartree. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 15–23, 2000     

Comparison of ground and excited state polarizabilities of thiophene,cyclopentadiene, and fulvene oligomers and their cyano substituted derivatives—Ab initio study

Ab Initio study of the ground and excited state polarizabilities of thiophene, fulvene, and cyclopentadiene based conducting oligomers and their cyano derivatives have been performed using the restricted Hartree–Fock (RHF) and the configuration interaction singles (CIS) approaches, respectively, with 3‐21G* basis set. For comparison purposes, for some small oligomers (monomers and dimers), higher basis sets (6‐31G*, 6‐31+G*, aug‐cc‐pVTZ) were also employed in the computations of polarizabilities. The trends in polarizability as a function of oligomer length were investigated. For all systems, the RHF polarizability increases as n^1.2–1.3 as n gets larger and the CIS polarizability increases as n^1.4–1.6 for n less than seven or eight rings and then increases approximately linearly with n for larger n. For the thiophene based systems the dependence of the polarizability on bond length alternation (BLA) along the backbone of the oligomers was also investigated using the RHF, density functional (DFT), and CIS theories (with 3‐21G* basis set). For thiophene dimer, we also performed RHF/aug‐cc‐pVTZ calculations of polarizabilities versus BLA. We found that the polarizability is largest when BLA is near zero (for both ground and excited states), which correlates with the lowest excitation energy. Comparison with experimental results has been made where possible. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 1983–1995, 2007