Source author record

Piotr Gniewek

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

4works
3topics
1close collaborators

Published work

4 publications

preprint2016arXiv

Convergence properties of the multipole expansion of the exchange contribution to the interaction energy

The conventional surface integral formula Jsurf[Φ]J_{\rm surf}[Φ] and an alternative volume integral formula Jvar[Φ]J_{\rm var}[Φ] are used to compute the asymptotic exchange splitting of the interaction energy of the hydrogen atom and a proton employing the primitive function ΦΦ in the form of its truncated multipole expansion. Closed-form formulas are obtained for the asymptotics of Jsurf[ΦN]J_{\rm surf}[Φ_N] and Jvar[ΦN]J_{\rm var}[Φ_N], where ΦNΦ_N is the multipole expansion of ΦΦ truncated after the 1/RN1/R^N term, RR being the internuclear separation. It is shown that the obtained sequences of approximations converge to the exact results with the rate corresponding to the convergence radius equal to 2 and 4 when the surface and the volume integral formulas are used, respectively. When the multipole expansion of a truncated, KKth order polarization function is used to approximate the primitive function the convergence radius becomes equal to unity in the case of Jvar[Φ]J_{\textrm{var}}[Φ]. At low order the observed convergence of Jvar[ΦN]J_{\rm var}[Φ_N] is, however, geometric and switches to harmonic only at certain value of N=NcN=N_c dependent on KK. An equation for NcN_c is derived which very well reproduces the observed KK-dependent convergence pattern. The results shed new light on the convergence properties of the conventional SAPT expansion used in applications to many-electron diatomics.

preprint2016arXiv

Determination of the exchange interaction energy from the polarization expansion of the wave function

The exchange contribution to the energy of the hydrogen atom interacting with a proton is calculated from the polarization expansion of the wave function using the conventional surface-integral formula and two formulas involving volume integrals: the formula of the symmetry-adapted perturbation theory (SAPT) and the variational formula recommended by us. At large internuclear distances RR, all three formulas yield the correct expression −(2/e)Re−R-(2/e)Re^{-R}, but approximate it with very different convergence rates. In the case of the SAPT formula, the convergence is geometric with the error falling as 3−K3^{-K}, where KK is the order of the applied polarization expansion. The error of the surface-integral formula decreases exponentially as aK/(K+1)!a^K/(K+1)!, where a=ln⁡2−12a=\ln2 -\tfrac{1}{2}. The variational formula performs best, its error decays as K1/2[aK/(K+1)!]2K^{1/2} [a^{ K}/(K+1)!]^2. These convergence rates are much faster than those resulting from approximating the wave function through the multipole expansion. This shows the efficiency of the partial resummation of the multipole series effected by the polarization expansion. Our results demonstrate also the benefits of incorporating the variational principle into the perturbation theory of molecular interactions.

preprint2015arXiv

Asymptotics of the exchange splitting energy for a diatomic molecular ion from a volume integral formula of symmetry-adapted perturbation theory

The exchange splitting energy JJ of the lowest \emph{gerade} and \emph{ungerade} states of the H2+_2^+ molecular ion was calculated using a volume integral expression of symmetry-adapted perturbation theory and standard basis set techniques of quantum chemistry. The performance of the proposed expression was compared to the well known surface integral formula. Both formulas involve the primitive function which we calculated employing either the Hirschfelder-Silbey perturbation theory or the conventional Rayleigh-Schrödinger perturbation theory (the polarization expansion). Our calculations show that very accurate values of JJ can be obtained using the proposed volume integral formula. When the Hirschfelder-Silbey primitive function is used in both formulas the volume formula gives much more accurate results than the surface integral expression. We also show that using the volume integral formula with the primitive function approximated by Rayleigh-Schrödinger perturbation theory, one correctly obtains only the first four terms in the asymptotic expansion of the exchange splitting energy.

preprint2015arXiv

Exchange splitting of the interaction energy and the multipole expansion of the wave function

The exchange splitting JJ of the interaction energy of the hydrogen atom with a proton is calculated using the conventional surface-integral formula Jsurf[φ]J_{\textrm{surf}}[φ], the volume-integral formula of the symmetry-adapted perturbation theory JSAPT[φ]J_{\textrm{SAPT}}[φ], and a variational volume-integral formula Jvar[φ]J_{\textrm{var}}[φ]. The calculations are based on the multipole expansion of the wave function φφ, which is divergent for any internuclear distance RR. Nevertheless, the resulting approximations to the leading coefficient j0j_0 in the large-RR asymptotic series J(R)=2e−R−1R(j0+j1R−1+j2R−2+⋯ )J(R) = 2 e^{-R-1} R ( j_0 + j_1 R^{-1} + j_2 R^{-2} +\cdots ) converge, with the rate corresponding to the convergence radii equal to 4, 2, and 1 when the Jvar[φ]J_{\textrm{var}}[φ], Jsurf[φ]J_{\textrm{surf}}[φ], and JSAPT[φ]J_{\textrm{SAPT}}[φ] formulas are used, respectively. Additionally, we observe that also the higher jkj_k coefficients are predicted correctly when the multipole expansion is used in the Jvar[φ]J_{\textrm{var}}[φ] and Jsurf[φ]J_{\textrm{surf}}[φ] formulas. The SAPT formula JSAPT[φ]J_{\textrm{SAPT}}[φ] predicts correctly only the first two coefficients, j0j_0 and j1j_1, gives a wrong value of j2j_2, and diverges for higher jnj_n. Since the variational volume-integral formula can be easily generalized to many-electron systems and evaluated with standard basis-set techniques of quantum chemistry, it provides an alternative for the determination of the exchange splitting and the exchange contribution of the interaction potential in general.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.