---
title: 'MP²A: MP2 Predictor & Multi-Precision Arithmetic'
url: https://www.emergentmind.com/topics/mp-2a
type: topic
---

# MP²A: MP2 Predictor & Multi-Precision Arithmetic

Searching arXiv for the cited works and the term to ground the article.
arXiv search query: 2001.06364 MP2 accuracy predictor weak interactions
“MP$^2$A” is not presented in the cited arXiv sources as the formal name of a single standardized method or software system. In the materials associated with this label, two distinct contexts appear. One is **MAP**, the **MP2 accuracy predictor** for weak interactions, which estimates the reliability of second-order Møller–Plesset perturbation theory from the structure of the Hartree–Fock adiabatic connection [2001.06364]. The other is **MP**, Brent’s portable Fortran package for **multiple-precision floating-point arithmetic** and special functions [1004.3173]. This suggests that “MP$^2$A” is best understood as an ambiguous designation rather than a canonical technical term.

## 1. Nomenclature and scope

The most direct chemically relevant interpretation is **MAP**, introduced as the **MP2 accuracy predictor**. MAP is described as a simple, physically motivated diagnostic for how reliable **second-order Møller–Plesset perturbation theory (MP2)** is for **weak interactions**, especially **noncovalent binding energies** [2001.06364]. Its purpose is not to replace MP2, but to indicate when the MP2 approximation is likely to be trustworthy and when it should be treated with caution.

A separate and unrelated interpretation arises from Brent’s **MP** package. That work does **not** describe a separate package called “MP$^2$A”; rather, it describes **MP** itself, a package of ANSI Standard Fortran subroutines for **multiple-precision floating-point arithmetic** and the evaluation of elementary and special functions [1004.3173]. In this reading, the notation “MP$^2$A” can only be treated as an interpretive label associated with MP and its extended arithmetic capabilities, not as the formal name of a distinct package.

## 2. Adiabatic-connection basis of the MP2 accuracy predictor

MAP is built from the **Hartree–Fock (HF) adiabatic connection (AC)**. In that framework one defines a coupling-constant dependent Hamiltonian
\[
\hat{H}_{\lambda}^{\rm HF}=\hat{T}+\hat{V}_{\rm ext}+\lambda \hat{V}_{ee}+ \Big(1-\lambda \Big) \Big(\hat{J}+\hat{K} \Big),
\]
where \(\hat{T}\) is kinetic energy, \(\hat{V}_{\rm ext}\) is the external potential, \(\hat{V}_{ee}\) is the electron-electron repulsion, and \(\hat{J}\), \(\hat{K}\) are the Hartree and exchange operators computed from the physical HF reference and held fixed along the path [2001.06364]. The corresponding correlation energy is obtained from
\[
E_{\rm c}= \int_0^1 W_{c,\lambda}\, d \lambda,
\]
with HF-AC integrand
\[
W_{c, \lambda}= \langle \Psi_\lambda | \hat{V}_{ee} - \hat{J}-\hat{K} | \Psi_\lambda \rangle - \langle \Psi_0 | \hat{V}_{ee} - \hat{J}-\hat{K} | \Psi_0\rangle .
\]

Expanding the integrand around \(\lambda=0\) gives the perturbation series, and truncation at second order yields
\[
W_{\lambda,c}^{\rm HF} \approx 2 E_c^{\rm MP2}\,\lambda.
\]
In this sense, **MP2 is exactly the linear approximation to the exact HF AC curve near \(\lambda=0\)** [2001.06364]. The central physical interpretation follows directly: if the exact HF AC curve is close to a straight line between \(\lambda=0\) and \(\lambda=1\), then MP2 should be good; if the exact AC curve has strong curvature, MP2 will be less reliable.

The paper also discusses an earlier exact indicator based on comparing the endpoint \(W_{c,1}\) with the initial slope \(W'_{c,0}\). For a perfectly linear AC curve, the indicator equals \(1\), and MP2 is exact; convexity or concavity shifts the value away from \(1\), thereby quantifying nonlinearity. The helium isoelectronic series illustrates this idea: as nuclear charge \(Z\) increases, the AC curve becomes more linear and the indicator tends to \(1\), consistent with the increasing accuracy of MP2 in the large-\(Z\) limit [2001.06364].

## 3. Interpolated construction of MAP

The exact indicator is not practical for routine use because it requires the fully interacting wavefunction \(\Psi\), through \(W_{c,1}\). MAP is obtained by replacing the exact endpoint with an interpolated estimate along the HF adiabatic connection using the **Seidl–Perdew–Levy (SPL)** form
\[
W_{c,\lambda}^{\rm SPL}(\mathbf{W})= W_{c,\infty} \left( 1- \left( 1+\frac{4 E_c^{\rm PT2}\lambda}{W_{c,\infty} \right)^{-1/2} \right),
\]
with input set
\[
\mathbf{W}=\{W_0,W'_0,W_\infty\}, \qquad W_0=E_x,\quad W_0'=2E_c^{\rm PT2},\quad W_{c,\infty}=W_\infty-W_0.
\]
For the HF case, \(E_c^{\rm PT2}=E_c^{\rm MP2}\) [2001.06364].

For a complex \(M\) composed of fragments \(F_i\), the interaction AC curve is defined as
\[
W_{\lambda,c}^{\rm int}(M)=W_{\lambda,c}(M)-\sum_{i=1}^N W_{\lambda,c}(F_i).
\]
The interpolated interaction curve is then
\[
W_{c,\lambda}^{\rm SPL,int}(M)= W_{c,\lambda}^{\rm SPL}(M) - W_{c,\lambda}^{\rm SPL}\Big(\sum_i^N F_i\Big).
\]
The construction works cleanly for systems whose fragments have **non-degenerate ground states**, because under that condition the interaction AC vanishes at infinite separation [2001.06364].

The practical predictor is defined through
\[
X^{\rm SPL}=\frac{W_{c,\lambda=1}^{\rm SPL,int}(M)}{2E_c^{\rm MP2}(M)-2 \sum_{i}^N E_c^{\rm MP2}(F_i)},
\qquad
MAP = 1-X^{\rm SPL}.
\]
In operational terms, MAP is a **dimensionless measure of how far the HF AC interaction curve deviates from the MP2 straight-line picture**, with the practical endpoint supplied by interpolation rather than by the exact wavefunction [2001.06364].

The interpolation requires a large-\(\lambda\) input. For that purpose the paper uses the **point-charge-plus-continuum (PC)** model,
\[
W_{\infty}^{\rm HF} \approx W_{\infty}^{\rm PC}[\rho^{\rm HF}] = \int \left[ A\rho^{\rm HF}(\mathbf{r})^{4/3} + B\frac{|\nabla \rho^{\rm HF}(\mathbf{r})|^2}{\rho^{\rm HF}(\mathbf{r})^{4/3}} \right]\, d\mathbf{r},
\]
with
\[
A=-1.451, \qquad B=5.317\times10^{-3}.
\]
This is a semilocal approximation evaluated on the HF density, and the authors note that the “bare” PC model works best in this context [2001.06364].

## 4. Computational characteristics, tested systems, and observed behavior

A notable property of MAP is that it comes at **negligible extra cost** beyond the MP2 calculation itself [2001.06364]. The required ingredients are the usual **MP2 correlation energy**, the HF density \(\rho^{\rm HF}\), its gradient \(\nabla \rho^{\rm HF}\), and the PC-model estimate \(W_{\infty}^{\rm PC}[\rho^{\rm HF}]\). This makes the predictor suitable as a screening tool for MP2 reliability in weakly bound complexes.

The method is illustrated on the **helium isoelectronic series**—H\(^-\), He, Be\(^{2+}\), Ne\(^{8+}\)—and on noncovalent dimers including the **benzene dimer** and **acetic acid dimer** [2001.06364]. Its main practical assessment is on the **S22** and **S66** noncovalent benchmark sets. The reported trend is monotonic in the intended sense: as **MAP decreases toward 0**, the AC curve becomes more linear and **MP2 becomes more accurate**.

The benchmark observations are specific. When **MAP < 0.20**, the relative MP2 errors are **always below 25%** in the reported tests; around **MAP \(\sim 0.25\)**, MP2 errors can become very large, **up to about 80%** [2001.06364]. The large-MAP cases are often **stacking complexes**, which are known failure cases for MP2, whereas **hydrogen-bonded systems** tend to have smaller MAP values and better MP2 performance.

For the **S66** subsets, mean absolute errors and mean absolute relative errors are reported separately for **hydrogen-bonded**, **mixed/others**, and **dispersion-dominated** complexes. The averaged MAP increases in the order
\[
\text{H-bonded} \;\to\; \text{others} \;\to\; \text{dispersion},
\]
and MP2 accuracy decreases correspondingly [2001.06364]. This aligns with the established observation that MP2 is less reliable for dispersion-driven and stacked systems than for many hydrogen-bonded complexes.

## 5. Applicability, limits, and interpretive status

MAP is intended for **interaction energies / weak interactions**, not for arbitrary strongly correlated or multireference problems [2001.06364]. The predictor is expected to be useful when the system is a **weakly interacting complex**, the fragments have **non-degenerate ground states**, the interaction can be reasonably described by the **HF AC plus interpolation**, and the MP2 error is mainly governed by the **curvature** of the AC curve. Under those conditions, a small MAP functions as a compact warning signal of near-linearity and probable MP2 reliability, whereas a large MAP signals pronounced curvature and probable unreliability.

The stated limitations are explicit. The current predictor is **restricted to systems dissociating into fragments with non-degenerate ground states**. The SPL interpolation is **not exact** for the HF AC curve itself, but is used as a practical surrogate. For fragments with degeneracy or near-degeneracy, a different treatment is needed, and future work is mentioned on small covalently bonded diatomics to extend the approach [2001.06364].

The interpretive status of MAP is therefore diagnostic rather than foundational. It does not alter the formal definition of MP2; instead, it uses AC geometry to anticipate when the MP2 approximation is likely to succeed or fail. A plausible implication is that MAP is most informative in exactly those regimes where weak interaction energies are inexpensive enough to compute at MP2 level, yet sufficiently delicate that qualitative misranking or substantial overbinding remains a practical concern.

## 6. Multiple-precision arithmetic: the separate MP context

In a separate line of work relevant to the designation “MP$^2$A,” **MP** is a package of ANSI Standard Fortran subroutines for **multiple-precision floating-point arithmetic** and the evaluation of elementary and special functions [1004.3173]. It is intended to let a Fortran program work with real numbers at arbitrarily high precision, subject to storage limitations and the chosen number of digits. Brent emphasizes machine independence or near machine independence: the package should run on any machine with an **ANSI Fortran (1966)** compiler, enough memory, and at least a **16-bit integer wordlength**.

The key runtime parameters are the base \(B\) and the number of digits \(T\), both user-selectable and dynamically changeable. A multiple-precision number is stored in an integer array in unpacked form as
\[
\text{word 1} = \text{sign},\qquad \text{word 2} = \text{exponent},\qquad \text{words } 3 \text{ to } T+2 = \text{fraction digits in base } B.
\]
Word 1 contains the sign \(0\), \(-1\), or \(+1\); word 2 contains the exponent to base \(B\); and the remaining words store the normalized fraction with one base-\(B\) digit per word [1004.3173]. A packed format is also supported, with two base-\(B\) digits per word, to save storage.

The package supports **truncated arithmetic**, **rounded-to-nearest arithmetic**, and **directed rounding**. The rounding mode is controlled by `RNDRL` in `COMMON /MPCOM/`, where `0` means truncate/chop, `1` means round to nearest with ties to even, `2` means round down, and `3` means round up [1004.3173]. For rounded arithmetic, the basic operations produce correctly rounded \(T\)-digit base-\(B\) results; for truncated arithmetic, the guide states a relative error on the order of \(B^{1-T}\).

Global control is organized through the labeled common block
\[
\texttt{COMMON /MPCOM/},
\]
which stores parameters such as `BASE`, `NUMDIG`, `MAXEXP`, `LUN`, workspace pointers, and controls for precision, formatting, rounding, and underflow [1004.3173]. Exponents range from \(1-M\) to \(+M\), so representable numbers lie roughly between \(B^{-M}\) and almost \(B^M\). Underflow is handled by `MPUNFL`, which usually sets the result to zero, while overflow calls `MPOVFL` and terminates execution with an error. The package also records the minimum and maximum exponents encountered in a run through `MNEXPN` and `MXEXPN`.

The technical restrictions are explicit: \(B \ge 2\), \(T \ge 2\), \(M > T\), \(8B^2 - 1\) must fit in the machine’s largest integer, and \(B^{T-1} \ge 10^7\) is recommended for practical decimal accuracy [1004.3173]. MP avoids dependence on standard library functions such as `SIN`, `EXP`, or `ALOG`, because one goal is to use MP to **test** those routines rather than depend on them. It also avoids reliance on `REAL` and `DOUBLE PRECISION` except in a few conversion routines.

Usability is extended by the **Augment precompiler**, which lets Fortran code declare variables as `MULTIPLE` or `MULTIPAK` and then write expressions such as
\[
X = Y + Z*\exp(X+1)/Y
\]
instead of manually invoking routines like `MPADDI`, `MPEXP`, `MPMUL`, `MPDIV`, and `MPADD` [1004.3173]. The package also provides a broad function library, including arithmetic, powers and roots, elementary functions, special functions, constants, input/output, conversions, and comparisons. Representative routines include `MPADD`, `MPMUL`, `MPDIV`, `MPEXP`, `MPLN`, `MPPI`, `MPBESJ`, `MPERF`, `MPGAM`, and `MPZETA`. Algorithmic details are given for many routines; for example, `MPEXP` and `MPLN` have time \(O(\sqrt{T}\,m(T))\), whereas `MPMUL` uses the straightforward \(O(T^2)\) multiplication algorithm with guard digits [1004.3173].

Taken together, these two contexts show that “MP$^2$A” does not denote a single unified formal object in the cited literature. In chemical electronic-structure theory it aligns most naturally with **MAP**, a diagnostic for MP2 reliability derived from HF adiabatic-connection curvature [2001.06364]. In numerical computation it may evoke the **MP** ecosystem for multiple-precision arithmetic, but that work describes **MP**, not a distinct system formally named “MP$^2$A” [1004.3173].

Source: https://www.emergentmind.com/topics/mp-2a