Energy-density relations in momentum space. III. Variational aspect

A variational formula for the momentum density is derived by using the Hellmann-Feynman theorem and by introducing a reference system whose Hamiltonian differs only in the kinetic energy part from that of the original system. As simple applications of the present results, the reduced mass correction and the relativistic correction for the hydrogen-like atom are discussed.     

Interatomic interactions in momentum space. Momentum density and kinetic energy in chemical bonding

On the basis of the virial theorem for a uniform scaling process of a polyatomic system, the total energy and its gradient are quantitatively related with the behavior of the electron density in momentum space through the kinetic energy of the system. For attractive and repulsive interactions, the behavior of the momentum density distribution and its effect on the stabilization energy and the interatomic force are examined. Some guiding principles are deduced for their interrelation. The results are used to clarify the role of kinetic energy in chemical bonding. Possible energy partitioning in this approach is also mentioned.     

Analytic geometrical derivatives of second-order molecular properties from perturbation theory

Assuming the Born-Oppenheimer approximation for molecular wavefunctions satisfies the Hellmann-Feynman theorem, Rayleigh-Schrödinger perturbation theory is employed to develop an analytic formula for derivatives of expectation values and second-order properties with respect to nuclear coordinates.     

Analytic derivative couplings between configuration-interaction-singles states with built-in electron-translation factors for translational invariance

We present a method for analytically calculating the derivative couplings between a pair of configuration-interaction-singles (CIS) excited states obtained in an atom-centered basis. Our theory is exact and has been derived using two completely independent approaches: one inspired by the Hellmann-Feynman theorem and the other following from direct differentiation. (The former is new, while the latter is in the spirit of existing approaches in the literature.) Our expression for the derivative couplings incorporates all Pulay effects associated with the use of an atom-centered basis, and the computational cost is minimal, roughly comparable to that of a single CIS energy gradient. We have validated our method against CIS finite-difference results and have applied it to the lowest lying excited states of naphthalene; we find that naphthalene derivative couplings include Pulay contributions sufficient to have a qualitative effect. Going beyond standard problems in analytic gradient theory, we have also constructed a correction, based on perturbative electron-translation factors, for including electronic momentum and eliminating spurious components of the derivative couplings that break translational symmetry. This correction is general and can be applied to any level of electronic structure theory.     

Force in scf theories

A new force approach is reported in which we improve the wavefunction so that the Hellmann-Feynman theorem is satisfied. A sufficient condition for the Hellmann-Feynman theorem to be satisfied is that the basis set includes AO derivatives ?x_r/?x_rfor any basis x_r. Here we test a procedure in which only the first derivative AOs are added to the “parent” AOs. The results are very encouraging.     

Atoms-in-molecules in momentum space: a Hirshfeld partitioning of electron momentum densities

A direct application of the Hirshfeld atomic partitioning (HAP) scheme is implemented for molecular electron momentum densities (EMDs). The momentum density contributions of individual atoms in diverse molecular systems are analyzed along with their topographical features and the kinetic energies of the atomic partitions. The proposed p-space HAP-based charge scheme does seem to possess the desirable attributes expected of any atoms in molecules partitioning. In addition to this, the main strength of the p-space HAP is the exact knowledge of the kinetic energy functional and the inherent ease in computing the kinetic energy. The charges derived from HAP in momentum space are found to match chemical intuition and the generally known chemical characteristics such as electronegativity, etc.     

On the local representation of the electronic momentum operator in atomic systems

The local quantum theory is applied to the study of the momentum operator in atomic systems. Consequently, a quantum-based local momentum expression in terms of the single-electron density is determined. The limiting values of this function correctly obey two fundamental theorems: Kato's cusp condition and the Hoffmann-Ostenhof and Hoffmann-Ostenhof exponential decay. The local momentum also depicts the electron shell structure in atoms as given by its local maxima and inflection points. The integration of the electron density in a shell gives electron populations that are in agreement with the ones expected from the Periodic Table of the elements. The shell structure obtained is in agreement with the higher level of theory computations, which include the Kohn-Sham kinetic energy density. The average of the local kinetic energy associated with the local momentum is the Weizsacker kinetic energy. In conclusion, the local representation of the momentum operator provides relevant information about the electronic properties of the atom at any distance from the nucleus.     

Momentum density interpretation of the H(1s)–H+ long-range force

The recently proposed method of using momentum densities for interatomic interactions is applied to the long-range force between the ground-state hydrogen atom and the proton, and the results are compared with those from using the position density based on the electrostatic Hellmann–Feynaman theorem. A new physical interpretation of the long-range force is obtained, which is complementary to that in position space. It is found that some perturbative changes in the position density do not accompany changes in the momentum density.     

Electron momentum density and band structure calculations of α- and β-GeTe

We have measured isotropic experimental Compton profile of α-GeTe by employing high energy (662 keV) γ-radiation from a ¹³⁷Cs isotope. To compare our experiment, we have also computed energy bands, density of states, electron momentum densities and Compton profiles of α- and β-phases of GeTe using the linear combination of atomic orbitals method. The electron momentum density is found to play a major role in understanding the topology of bands in the vicinity of the Fermi level. It is seen that the density functional theory (DFT) with generalised gradient approximation is relatively in better agreement with the experiment than the local density approximation and hybrid Hartree–Fock/DFT.     

The outer valance orbital electron densities of cyclopentane by binary (e,2e) spectroscopy

The binding energy spectra and electron distributions in momentum space of the valence orbitals of cyclopentane (C(5)H(10)) are studied by Electron Momentum Spectroscopy (EMS) in a noncoplanar symmetric geometry. The impact energy was 1200 eV plus binding energy and energy resolution of the EMS spectrometer was 1.2 eV. The experimental momentum profiles of the outer valence orbitals are compared with the theoretical momentum distributions calculated using Hartree-Fock and density functional theory (DFT) methods. The shapes of the experimental momentum distributions are generally quite well described by both the Hartree-Fock and DFT calculations when the large and diffuse basis sets are used.     

Electron Momentum Spectroscopy Study on Inner Orbitals of Methyl iodide

The binding energy spectrum and electron momentum profiles of the inner orbitals of methyl iodide have been measured using an electron momentum spectrometer at the impact energy of 1200 eV plus binding energy. Two peaks in the binding energy spectrum, arising from the spin-orbit splitting, are observed and the corresponding electron momentum profiles are obtained. Relativistic density functional calculations are performed to elucidate the experimental electron momentum profiles of two spin-orbit splitting components, showing agreement with each other except for the intensity in low momentum region. The measured high intensity in the low momentum region can be further explained by the distorted wave calculation.     

Vibronic coupling in cyclopentadienyl radical: a method for calculation of vibronic coupling constant and vibronic coupling density analysis

A method of calculation of vibronic or electron-phonon coupling constant is presented for a Jahn-Teller molecule, cyclopentadienyl radical. It is pointed out that symmetry breaking at degenerate point and violation of Hellmann-Feynman theorem occur in the calculations based on a single Slater determinant. In order to overcome these difficulties, the electronic wave functions are calculated using generalized restricted Hartree-Fock and complete active space self-consistent-field method and the couplings are computed as matrix elements of the electronic operator of the vibronic coupling. Our result agrees well with the experimental and theoretical values. A concept of vibronic coupling density is proposed in order to explain the order of magnitude of the coupling constant from view of the electronic and vibrational structures. It illustrates the local properties of the coupling and enables us to control the interaction. It could open a way to the engineering of vibronic interactions.     

Rototranslational sum rules for electromagnetic hypershielding at the nuclei and related atomic Cartesian derivatives of the optical rotatory power

Two molecular properties, the nuclear electromagnetic hypershielding (psi(gamma,alphabeta) ('I)) and the gradient of the electric dipole-magnetic dipole polarizability (nabla(Igamma)G(alphabeta) (')), have been calculated using the time-dependent Hartree-Fock method. Provided the Hellmann-Feynman theorem is satisfied, these quantities are equivalent and are related through the nabla(Igamma)G(alphabeta) (')=eZ(I)psi(gamma,alphabeta) ('I) relation, where Z(I) is the atomic number of atom I and e the magnitude of the electron charge. In such a case, the determination of the nuclear electromagnetic hypershielding presents the computational advantage over the evaluation of the gradient of G(alphabeta) (') of requiring only the knowledge of nine mixed second-order derivatives of the density matrix with respect to both electric and magnetic fields (D(alpha,beta)(-omega,omega)) instead of the 3N (N is the number of atoms) derivatives of the density matrix with respect to the Cartesian coordinates (D(Igamma)). It is shown here for the H(2)O(2) molecule that very large basis sets such as the aug-cc-pVQZ or the R12 basis are required to satisfy the Hellmann-Feynman theorem. These basis set requirements have been substantiated by considering the corresponding rototranslational sum rules. The origin dependence of the rototranslational sum rules for the gradient of G(alphabeta) (') has then been theoretically described and verified for the H(2)O(2) molecule.     

Floating functions satisfying the hellmann—feynman theorem: Single floating scheme

The electrostatic calculation for molecules using approximated variational wave functions leads to well known difficulties connected with the application of the Hellmann-Feynman (H? F) theorem. This is due to the basis set inadequacies in the underlying calculations. This defect can easily be remedied by floating functions, whose centers are optimized in space. We can keep almost everything of the traditional wave function with a nuclear-fixed basis set, but we apply single floating to ensure the H? F theorem. Then, one can obtain a wave function obeying the H? F theorem. This provides a great conceptual simplification and may lead to practical advantages. The single floating scheme, which retains one expansion center per nucleus, is successfully applied to a series of small molecules using SCF and CASSCF wave functions with sufficiently polarized basis sets.     

Multiresolution quantum chemistry in multiwavelet bases: Analytic derivatives for Hartree-Fock and density functional theory

An efficient and accurate analytic gradient method is presented for Hartree-Fock and density functional calculations using multiresolution analysis in multiwavelet bases. The derivative is efficiently computed as an inner product between compressed forms of the density and the differentiated nuclear potential through the Hellmann-Feynman theorem. A smoothed nuclear potential is directly differentiated, and the smoothing parameter required for a given accuracy is empirically determined from calculations on six homonuclear diatomic molecules. The derivatives of N2 molecule are shown using multiresolution calculation for various accuracies with comparison to correlation consistent Gaussian-type basis sets. The optimized geometries of several molecules are presented using Hartree-Fock and density functional theory. A highly precise Hartree-Fock optimization for the H2O molecule produced six digits for the geometric parameters.     

Observation of Interference Effects in (e, 2e) Electron Momentum Spectroscopy of SF6 (cited: 1)

Two-dimensional electron density map (2D map) of binding energy and relative azimuthal angle (i.e., momentum) for the outer-valence molecular orbitals of SF₆ has been measured by a highly sensitive electron momentum spectrometer with noncoplanar symmetric geometry at the impact energy of 1.2 keV plus binding energy. The experimental electron momentum profiles for the relevant molecular orbitals have been extracted from the 2D map and interpreted on the basis of the quantitative calculations using the density functional theory with B3LYP hybrid functional. For the outermost F2p nonbonding orbitals of SF₆, the interference patterns are clearly observed in the ratios of the electron momentum profiles of molecular orbitals to that of atomic F2p orbital.     

Atomic quantum similarity indices in position and momentum spaces

Quantum similarity for atoms is investigated using electron densities in position and momentum spaces. Contrary to the results in position space, the analysis in the momentum space shows how the momentum density carries fundamental information about periodicity and structure of the system and reveals the pattern of Mendeleev's table. A global analysis in the joint r-p space keeps this result.     

Determination of noncovalent interaction energies from electronic densities

Starting from the Hellmann-Feynman theorem, an expression is derived for the interaction energy in forming a noncovalently bound complex. It invokes only classical electrostatics, involving the charge distributions of the components as they are in the complex. We propose a method for obtaining these from the total charge distribution. Integration over the electronic densities is carried out by a numerical procedure slightly modified from that of Gavezzotti. We calculate the interaction energies for four molecular dimers at a variety of computational levels. The results are analyzed and compared to the best estimated values available in the literature.