No Arabic abstract
The Hartree-Fock problem was recently recast as a semidefinite optimization over the space of rank-constrained two-body reduced-density matrices (RDMs) [Phys. Rev. A 89, 010502(R) (2014)]. This formulation of the problem transfers the non-convexity of the Hartree-Fock energy functional to the rank constraint on the two-body RDM. We consider an equivalent optimization over the space of positive semidefinite one-electron RDMs (1-RDMs) that retains the non-convexity of the Hartree-Fock energy expression. The optimized 1-RDM satisfies ensemble $N$-representability conditions, and ensemble spin-state conditions may be imposed as well. The spin-state conditions place additional linear and nonlinear constraints on the 1-RDM. We apply this RDM-based approach to several molecular systems and explore its spatial (point group) and spin ($S^2$ and $S_3$) symmetry breaking properties. When imposing $S^2$ and $S_3$ symmetry but relaxing point group symmetry, the procedure often locates spatial-symmetry-broken solutions that are difficult to identify using standard algorithms. For example, the RDM-based approach yields a smooth, spatial-symmetry-broken potential energy curve for the well-known Be--H$_2$ insertion pathway. We also demonstrate numerically that, upon relaxation of $S^2$ and $S_3$ symmetry constraints, the RDM-based approach is equivalent to real-valued generalized Hartree-Fock theory.
This work presents an algorithm to evaluate Coulomb and exchange matrices in Fock operator using range separation techniques at various aspects. This algorithm is particularly favorable for the scenario of (1) all-electron calculations or (2) computing exchange matrix for a large number of $mathbf{k}$-point samples. An all electron Hartree-Fock calculation with 110k basis functions is demonstrated in this work.
Restoration of pseudo-spin symmetry (PSS) along the $N=32$ and $34$ isotonic chains and the physics behind are studied by applying the relativistic Hartree-Fock theory with effective Lagrangian PKA1. Taking the proton pseudo-spin partners $(pi2s_{1/2},pi1d_{3/2})$ as candidates, systematic restoration of PSS along both isotonic chains is found from sulphur (S) to nickel (Ni), while distinct violation from silicon (Si) to sulphur is discovered near the drip lines. The effects of the tensor-force components introduced naturally by the Fock terms are investigated, which can only partly interpret the systematics from calcium to nickel, but fail for the overall trends. Further analysis following the Schr{o}dinger-like equation of the lower component of Dirac spinor shows that the contributions from the Hartree terms dominate the overall systematics of the PSS restoration, and such effects can be self-consistently interpreted by the evolution of the proton central density profiles along both isotonic chains. Specifically the distinct PSS violation is found to tightly relate with the dramatic changes from the bubble-like density profiles in silicon to the central-bumped ones in sulphur.
The relation between the correlation energy and the entanglement is analytically constructed for the Moshinskys model of two coupled harmonic oscillators. It turns out that the two quantities are far to be proportional, even at very small couplings. A comparison is made also with the 2-point Ising model.
We apply Projected Hartree-Fock theory (PHF) for approximating ground states of Heisenberg spin clusters. Spin-rotational, point-group and complex-conjugation symmetry are variationally restored from a broken-symmetry mean-field reference, where the latter corresponds to a product of local spin states. A fermionic formulation of the Heisenberg model furnishes a conceptual connection to PHF applications in quantum chemistry and detailed equations for a self-consistent field optimization of the reference state are provided. Different PHF variants are benchmarked for ground-state energies and spin-pair correlation functions of antiferromagnetic spin rings and three different polyhedra, with various values of the local spin-quantum number s. The low computational cost and the compact wave-function representation make PHF a promising complement to existing approaches for ground states of molecular spin clusters, particularly for large s and moderately large N. The present work may also motivate future explorations of more accurate post-PHF methods for Heisenberg spin clusters.
Four-component Dirac Hartree--Fock is an accurate mean-field method for treating molecular systems where relativistic effects are important. However, the computational cost and complexity of the two-electron interaction makes this method less common, even though we can consider the Dirac Hartree--Fock Hamiltonian the ground truth of electronic structure, barring explicit quantum-electrodynamical effects. Being able to calculate these effects is then vital to the design of lower scaling methods for accurate predictions in computational spectroscopy and properties of heavy element complexes that must include relativistic effects for even qualitative accuracy. In this work, we present a Pauli quaternion formalism of maximal component- and spin-separation for computing the Dirac-Coulomb-Gaunt Hartree--Fock ground state, with a minimal floating-point-operation count algorithm. This approach also allows one to explicitly separate different spin physics from the two-body interactions, such as spin-free, spin-orbit, and the spin-spin contributions. Additionally, we use this formalism to examine relativistic trends in the periodic table, and analyze the basis set dependence of atomic gold and gold dimer systems.