nDot.io physics
08-ilang-propagation / 07-bond-inertia
07bond inertiaverified

Determines how the cost of a bond made by oscillator cells depends on separation and joint motion, with which factor it enters the inertia of the bound pair, and whether the pair falls like its parts.

Code of version 1

check_bound_pair.py9 KB · raw
"""Check (ii) of 07-bond-inertia: bound level E_b(P) and its curvature from the exact
diagonalization of the effective joint contract (7.2)-(7.4) in the sector P, against
the formulas (7.10)-(7.14) of item 2.

Ring of L cells; bodies (sigma1, sigma2) = (+,+) with rest costs M1, M2 of their levels,
tau = 1; field c0 = tau/r, m_h = c0/ell.  Grips: square (4 F0 F2 = F1^2) and linear (F2 = 0).
F1 is set (proportional to c0) so that the static bond is the same for every r; two
binding strengths: a 'deep' one (a_b < ell) and a 'shallow' one (a_b > ell).
Self-energies of the single bodies are not included (they cancel in D at this order).
"""
import json
import time
from pathlib import Path

import numpy as np
from scipy.linalg import eigh

OUT = Path(__file__).parent / "output"
OUT.mkdir(exist_ok=True)
AIRY1 = 1.0187929716  # first zero of Ai'


def setup(L, ell, r, tau, M1, M2, F0, F1, F2, sg1=1, sg2=1):
    c0 = tau / r; mh = c0 / ell; kap = c0**2
    n = np.arange(L); qs = 2 * np.pi * n / L
    om = np.sqrt(mh**2 + 4 * kap * np.sin(qs / 2)**2)
    w0 = np.sum(1 / (2 * L * om)); abar = F0 + F2 * w0
    Mb1, Mb2 = abar * M1, abar * M2
    I1, I2 = sg1 * Mb1 / tau**2, sg2 * Mb2 / tau**2          # signed inertias (7.1)
    return dict(L=L, n=n, om=om, c0=c0, mh=mh, kap=kap, tau=tau, r=r, M1=M1, M2=M2,
                F0=F0, F1=F1, F2=F2, Mb1=Mb1, Mb2=Mb2, I1=I1, I2=I2, I=I1 + I2,
                mred=I1 * I2 / (I1 + I2), mu12=sg1 * sg2 * M1 * M2, sg1=sg1, sg2=sg2,
                rest=sg1 * Mb1 + sg2 * Mb2)


def hamiltonian(P, s, static=False):
    L, n, om, tau = s["L"], s["n"], s["om"], s["tau"]
    ks = 2 * np.pi * n / L
    p1 = P * s["I1"] / s["I"] + ks; p2 = P * s["I2"] / s["I"] - ks
    e1 = s["sg1"] * np.sqrt(s["Mb1"]**2 + tau**2 * np.sin(p1)**2)
    e2 = s["sg2"] * np.sqrt(s["Mb2"]**2 + tau**2 * np.sin(p2)**2)
    th1 = np.arctan2(tau * np.sin(p1), s["Mb1"]); th2 = np.arctan2(tau * np.sin(p2), s["Mb2"])
    if static:
        C = np.ones((L, L)); D1 = np.zeros((L, L)); D2 = np.zeros((L, L))
    else:
        C = np.cos((th1[:, None] + th1[None, :]) / 2) * np.cos((th2[:, None] + th2[None, :]) / 2)
        D1 = e1[None, :] - e1[:, None]; D2 = e2[:, None] - e2[None, :]   # rows k', cols k
    m = (n[None, :] - n[:, None]) % L                                   # q = k - k'
    wq2 = om[m]**2
    U = -s["mu12"] * s["F1"]**2 * C / (2 * L) * (1 / (wq2 - D1**2) + 1 / (wq2 - D2**2))
    ratio = max(np.max(np.abs(D1) / np.sqrt(wq2)), np.max(np.abs(D2) / np.sqrt(wq2)))
    if s["F2"] != 0:
        for mi in range(L):
            kp = (n - mi) % L
            ob = om[(mi - n) % L]; Om = om + ob; w = Om / (om * ob)
            d2 = np.concatenate([D1[kp, n], D2[kp, n]])**2
            S = (1 / (Om[None, :]**2 - d2[:, None])) @ w
            U[kp, n] += -s["mu12"] * s["F2"]**2 * C[kp, n] / (2 * L**2) * (S[:L] + S[L:])
            ratio = max(ratio, np.sqrt(d2.max()) / Om.min())
    if static:
        return U, ratio                      # static mutual kernel only
    return U + np.diag(e1 + e2), ratio


def analyse(s):
    L, tau, r = s["L"], s["tau"], s["r"]
    h = 0.015 * abs(s["I"]) * tau            # P <= 3h keeps the pair velocity P/I below 0.045 tau
    H0, ratio = hamiltonian(0.0, s)
    up = s["mred"] < 0                       # bound level above the band (mixed branches)
    sel = [L - 2, L - 1] if up else [0, 1]
    E, V = eigh(H0, subset_by_index=sel)
    if up:
        E, V = E[::-1], V[:, ::-1]
    E0, psi = E[0], V[:, 0]
    Ps = h * np.arange(1, 4); EP = []
    for P in Ps:
        HP, rt = hamiltonian(P, s); ratio = max(ratio, rt)
        j = L - 1 if up else 0
        EP.append(eigh(HP, eigvals_only=True, subset_by_index=[j, j])[0])
    EP = np.array(EP)
    A = np.stack([Ps**2, Ps**4, Ps**6], axis=1)
    a, b, c = np.linalg.solve(A, EP - E0)
    Ipair = 1 / (2 * a)
    # relative label measured from the nearest band minimum (0 or pi: the two valleys of Model D)
    ks = 2 * np.pi * s["n"] / L; ks = (ks + np.pi / 2) % np.pi - np.pi / 2
    T = np.sum(np.abs(psi)**2 * ks**2) / (2 * s["mred"])
    w_pi = np.sum(np.abs(psi[np.abs(2 * np.pi * s["n"] / L - np.pi) < np.pi / 2])**2)
    dE = E0 - s["rest"]
    Ust, _ = hamiltonian(0.0, s, static=True)     # static kernel U~(k-k')/L
    Uexp = psi @ Ust @ psi
    Uq = L * Ust[0, :]                             # U~(q) at q = 2 pi j / L
    Ur = np.real(np.fft.ifft(Uq))                  # U(r), r = 0..L-1
    rr = np.where(s["n"] > L // 2, s["n"] - L, s["n"])
    rUp = rr * (np.roll(Ur, -1) - np.roll(Ur, 1)) / 2
    psir = np.fft.ifft(psi) * np.sqrt(L)           # psi(r) = L^-1/2 sum_k psi_k e^{ikr}
    rUp_exp = np.sum(np.abs(psir)**2 * rUp)
    I1, I2 = s["I1"], s["I2"]
    D_num = Ipair - I1 - I2
    D_formula = ((2 - r**2) * dE + (1 - r**2) * T) / tau**2 + 4 * I1 * I2 * T
    D_general = (3 * T + 2 * Uexp) / tau**2 - (Uexp + rUp_exp) / s["c0"]**2 + 4 * I1 * I2 * T
    lam = s["mu12"] * s["F1"]**2 / (2 * s["kap"])
    g = -Uq[0]
    rad = np.sqrt(np.sum(np.abs(psir)**2 * rr**2))
    EPf = E0 + Ps**2 / (2 * (I1 + I2 + D_formula))
    return dict(r=r, F1=s["F1"], F2=s["F2"], ell=s["c0"] / s["mh"], L=L,
                branches=[s["sg1"], s["sg2"]], m_red=s["mred"],
                E_b0=E0, gap_to_next=abs(E[1] - E[0]), weight_pi_valley=w_pi, dE=dE, T=T,
                U_exp=Uexp, virial_2T=2 * T, virial_rUp=rUp_exp, rms_size=rad,
                dE_deep_asym=Ur[0] + np.sign(s["mred"]) * AIRY1
                * (lam**2 / (2 * abs(s["mred"])))**(1 / 3),
                dE_shallow_asym=-s["mred"] * g**2 / 2,
                I1=I1, I2=I2, I_pair=Ipair, D_num=D_num, D_formula=D_formula,
                D_general=D_general, D_deep_asym=(2 - r**2) * dE / tau**2,
                D_shallow_asym=dE / tau**2 - 4 * I1 * I2 * dE,
                rel_dev_D=abs(D_num - D_formula) / abs(D_formula),
                # D at first order in dE/M: I^2 (1/I - 1/I_pair) removes the trivial D^2/I term
                D_num_lin=(I1 + I2)**2 * (1 / (I1 + I2) - 1 / Ipair),
                dev_over_dE=abs(D_num - D_formula) / (abs(dE) / tau**2),
                dev_lin_over_dE=abs((I1 + I2)**2 * (1 / (I1 + I2) - 1 / Ipair) - D_formula)
                / (abs(dE) / tau**2),
                expected_scale=abs(dE) / abs(s["rest"]) + (tau / min(s["Mb1"], s["Mb2"]) / rad)**2,
                T_over_dE=T / dE,
                P=list(Ps), E_b_P_num=list(EP), E_b_P_formula=list(EPf),
                max_EP_dev_over_shift=float(np.max(np.abs(EP - EPf) / (EP - E0))),
                quartic_coeff=b, max_transfer_over_quantum=ratio,
                D_over_I=D_num / (I1 + I2))


import sys
L = int(sys.argv[1]) if len(sys.argv) > 1 else 1024
ell = float(sys.argv[2]) if len(sys.argv) > 2 else 30.0
tau, M1, M2, F0 = 1.0, 0.5, 0.4, 1.0
results = []
t0 = time.time()
runs = [((1, 1), g) for g in ("square", "linear")] + [((1, -1), "linear")]
# (name, ell, target size a_b): deep a_b = ell/3; shallow a_b = 3 ell; very shallow a_b = 11 ell
regimes = [("deep", ell, ell / 3), ("shallow", ell, 3 * ell), ("v-shallow", 8.0, 88.0)]
for (sg1, sg2), grip in runs:
    mred0 = abs((M1 * M2 / (sg1 * M1 + sg2 * M2))) / tau**2
    for strength, ell_, ab in regimes:
        if sg1 != sg2 and strength == "v-shallow":
            continue
        for r in (0.5, 1.0, 1.4):
            c0 = tau / r; mh = c0 / ell_
            if strength == "deep":
                lam = 1 / (2 * mred0 * ab**3); F1 = c0 * np.sqrt(2 * lam / (M1 * M2))
            else:
                gg = 1 / (mred0 * ab); F1 = mh * np.sqrt(gg / (M1 * M2))
            F2 = F1**2 / (4 * F0) if grip == "square" else 0.0
            s = setup(L, ell_, r, tau, M1, M2, F0, F1, F2, sg1, sg2)
            res = analyse(s); res.update(strength=strength, grip=grip)
            results.append(res)
            print(f"{sg1:+d}{sg2:+d} {strength:9s} {grip:6s} r={r:.2f} size={res['rms_size']:.1f} "
                  f"dE={res['dE']:.3e} T/dE={res['T_over_dE']:.3f} D_num={res['D_num']:.4e} "
                  f"D_lin={res['D_num_lin']:.4e} D_form={res['D_formula']:.4e} "
                  f"D_gen={res['D_general']:.4e} dev/dE={res['dev_over_dE']:.3f} "
                  f"devlin/dE={res['dev_lin_over_dE']:.3f} scale={res['expected_scale']:.3f} "
                  f"D/I={res['D_over_I']:.3f} EPdev={res['max_EP_dev_over_shift']:.3f} "
                  f"ratio={res['max_transfer_over_quantum']:.2f} gap={res['gap_to_next']:.2e} "
                  f"wpi={res['weight_pi_valley']:.2f} t={time.time() - t0:.0f}s", flush=True)
summary = {"parameters": {"L": L, "tau": tau, "M_k1": M1, "M_k2": M2, "F0": F0,
                          "regimes": regimes, "r": [0.5, 1.0, 1.4], "P": "h, 2h, 3h with h = 0.015 |I| tau"},
           "max_dev_lin_over_dE": max(x["dev_lin_over_dE"] for x in results),
           "max_dev_over_dE": max(x["dev_over_dE"] for x in results),
           "results": results}
(OUT / "bound_pair_check.json").write_text(json.dumps(summary, indent=1, default=float))
print("max dev/dE:", summary["max_dev_over_dE"], " max devlin/dE:", summary["max_dev_lin_over_dE"])
check_mutual_term.py7 KB · raw
"""Check (i) of 07-bond-inertia: the mutual term (7.3)-(7.4) against second-order
perturbation theory in the total contract, restricted to at most two field quanta.

Ring of L = 16 cells, two Model D bodies (clock level fixed, so S_i -> M_i sigma_3),
grip F(q) = F0 + F1 q + F2 q^2.  The unperturbed contract H0 contains the motion terms,
the first-order grip abar*M_i*sigma_3 (abar = F0 + F2 w0) and the field; the
perturbation is V = sum_i sum_x |x><x| M_i sigma_3 (dF(q_x) - F2 w0).  The mutual
term is the part of the symmetric (des Cloizeaux) second-order effective contract with
one vertex on each body.  Everything is built numerically in the position basis of the
bodies and the occupation basis of the field modes (<= 2 quanta).
"""
import json
from itertools import combinations_with_replacement
from pathlib import Path

import numpy as np

OUT = Path(__file__).parent / "output"
OUT.mkdir(exist_ok=True)


def fock(L):
    states = [()] + [(j,) for j in range(L)] + list(combinations_with_replacement(range(L), 2))
    index = {s: i for i, s in enumerate(states)}
    B = np.zeros((L, len(states), len(states)))
    for i, s in enumerate(states):
        for j in set(s):
            t = list(s); t.remove(j)
            B[j, index[tuple(sorted(t))], i] = np.sqrt(s.count(j))
    return states, B


def body(L, tau, Mbar):
    """Position-basis contract (index 2x+s) and its plane-wave eigenvectors."""
    s1 = np.array([[0, 1], [1, 0]], complex); s3 = np.diag([1.0, -1.0]).astype(complex)
    h = np.zeros((2 * L, 2 * L), complex)
    for x in range(L):
        y = (x + 1) % L
        h[2*y:2*y+2, 2*x:2*x+2] += 0.5 * tau * (-1j * s1)
        h[2*x:2*x+2, 2*y:2*y+2] += 0.5 * tau * (1j * s1)
        h[2*x:2*x+2, 2*x:2*x+2] += Mbar * s3
    W = np.zeros((2 * L, 2 * L), complex); e = np.zeros(2 * L); lab = []
    xs = np.arange(L)
    for n in range(L):
        p = 2 * np.pi * n / L
        pw = np.exp(1j * p * xs) / np.sqrt(L)
        blk = np.array([[np.vdot(np.kron(pw, np.eye(2)[a]), h @ np.kron(pw, np.eye(2)[b]))
                         for b in range(2)] for a in range(2)])
        ev, U = np.linalg.eigh(blk)
        for col, band in ((1, +1), (0, -1)):
            u = U[:, col]
            ref = u[0] if band == +1 else u[1]          # phase convention of Step 1
            u = u * abs(ref) / ref if abs(ref) > 1e-12 else u
            v = np.kron(pw, u); k = len(lab)
            W[:, k] = v; e[k] = ev[col]; lab.append((n, band))
    resid = np.max(np.abs(h @ W - W * e))
    return W, e, lab, resid


def run(case):
    L, kap, mh, tau = case["L"], case["kappa_h"], case["m_h"], case["tau"]
    M1, M2, sg1, sg2 = case["M1"], case["M2"], case["sigma1"], case["sigma2"]
    F0, F1, F2 = case["F0"], case["F1"], case["F2"]
    qs = 2 * np.pi * np.arange(L) / L
    om = np.sqrt(mh**2 + 4 * kap * np.sin(qs / 2)**2)
    w0 = np.sum(1 / (2 * L * om))
    abar = F0 + F2 * w0
    states, B = fock(L)
    EF = np.array([sum(om[j] for j in s) for s in states])
    vac = np.zeros(len(states)); vac[0] = 1
    X = []
    for x in range(L):
        qx = sum((B[j] * np.exp(1j * qs[j] * x) + B[j].T * np.exp(-1j * qs[j] * x))
                 / np.sqrt(2 * L * om[j]) for j in range(L))
        chi = F1 * (qx @ vac) + F2 * (qx @ (qx @ vac) - w0 * vac)
        X.append(chi)
    X = np.array(X)                                   # (x, fock)
    assert abs(X[:, 0]).max() < 1e-12                 # no vacuum component
    W1, e1, lab1, r1 = body(L, tau, abar * M1)
    W2, e2, lab2, r2 = body(L, tau, abar * M2)
    s3 = np.diag([1.0, -1.0])

    def Y(psi1, psi2):
        B1 = np.zeros((2 * L, L), complex); B2 = np.zeros((2 * L, L), complex)
        for x in range(L):
            B1[2*x:2*x+2, x] = M1 * s3 @ psi1[2*x:2*x+2]
            B2[2*x:2*x+2, x] = M2 * s3 @ psi2[2*x:2*x+2]
        V1 = np.einsum('ax,b,xc->abc', B1, psi2, X)
        V2 = np.einsum('a,bx,xc->abc', psi1, B2, X)
        tr = lambda T: np.einsum('ap,bq,abc->pqc', W1.conj(), W2.conj(), T)
        return tr(V1), tr(V2)

    En = e1[:, None, None] + e2[None, :, None] + EF[None, None, :]
    idx1 = {l: i for i, l in enumerate(lab1)}; idx2 = {l: i for i, l in enumerate(lab2)}

    def eps(M, sg, p):
        return sg * np.sqrt((abar * M)**2 + tau**2 * np.sin(p)**2)

    def th(M, p):
        return np.arctan2(tau * np.sin(p), abar * M)

    rows = []
    for (n1, n2, m) in case["labels"]:
        i1, i2 = idx1[(n1 % L, sg1)], idx2[(n2 % L, sg2)]
        f1, f2 = idx1[((n1 - m) % L, sg1)], idx2[((n2 + m) % L, sg2)]
        Yi = Y(W1[:, i1], W2[:, i2]); Yf = Y(W1[:, f1], W2[:, f2])
        Ei, Ef = e1[i1] + e2[i2], e1[f1] + e2[f2]
        cross = (Yf[0].conj() * Yi[1] + Yf[1].conj() * Yi[0])[:, :, 1:]   # Q: >= 1 quantum
        num = 0.5 * np.sum(cross * (1 / (Ei - En[:, :, 1:]) + 1 / (Ef - En[:, :, 1:])))
        # closed form (7.3)-(7.4)
        p1, p2, q = qs[n1 % L], qs[n2 % L], qs[m % L]
        p1f, p2f = p1 - q, p2 + q
        D1 = eps(M1, sg1, p1) - eps(M1, sg1, p1f); D2 = eps(M2, sg2, p2f) - eps(M2, sg2, p2)
        c1 = np.cos((th(M1, p1) + th(M1, p1f)) / 2); c2 = np.cos((th(M2, p2) + th(M2, p2f)) / 2)
        mu12 = sg1 * M1 * sg2 * M2
        wq = om[m % L]
        U1 = -mu12 * F1**2 * c1 * c2 / (2 * L) * (1 / (wq**2 - D1**2) + 1 / (wq**2 - D2**2))
        ob = om[(m - np.arange(L)) % L]; Om = om + ob
        U2 = -mu12 * F2**2 * c1 * c2 / (2 * L**2) * np.sum(Om / (om * ob) * (
            1 / (Om**2 - D1**2) + 1 / (Om**2 - D2**2)))
        form = U1 + U2
        rows.append({"labels_n1_n2_m": [n1, n2, m], "numerical": [num.real, num.imag],
                     "formula": form, "rel_dev": abs(num - form) / abs(form)})
    return {"case": {k: v for k, v in case.items() if k != "labels"}, "abar": abar,
            "eigvec_residual": max(r1, r2), "rows": rows,
            "max_rel_dev": max(r["rel_dev"] for r in rows)}


labels = [(0, 0, 1), (1, -2, 3), (3, -5, 2), (-2, 4, -1), (5, 2, 7), (2, 2, 0), (-3, 1, 0), (6, -6, 8)]
base = dict(L=16, kappa_h=1.0, m_h=0.3, tau=0.8, M1=0.7, M2=0.5, F0=1.0, F1=0.3, F2=0.2,
            sigma1=1, sigma2=1, labels=labels)
cases = [base,
         dict(base, F2=base["F1"]**2 / (4 * base["F0"])),          # square grip
         dict(base, sigma2=-1),                                     # mixed branches
         dict(base, tau=1.3)]                                       # r > 1
res = [run(c) for c in cases]
summary = {"max_rel_dev_all": max(r["max_rel_dev"] for r in res), "cases": res}
(OUT / "mutual_term_check.json").write_text(json.dumps(summary, indent=1))
for r in res:
    print(r["case"]["tau"], r["case"]["F2"], r["case"]["sigma2"], "max rel dev", r["max_rel_dev"],
          "resid", r["eigvec_residual"])
print("overall", summary["max_rel_dev_all"])
output/bound_pair_check.json39 KB · raw
{
 "parameters": {
  "L": 1024,
  "tau": 1.0,
  "M_k1": 0.5,
  "M_k2": 0.4,
  "F0": 1.0,
  "regimes": [
   [
    "deep",
    30.0,
    10.0
   ],
   [
    "shallow",
    30.0,
    90.0
   ],
   [
    "v-shallow",
    8.0,
    88.0
   ]
  ],
  "r": [
   0.5,
   1.0,
   1.4
  ],
  "P": "h, 2h, 3h with h = 0.015 |I| tau"
 },
 "max_dev_lin_over_dE": 0.15616620972696707,
 "max_dev_over_dE": 0.1705434907971377,
 "results": [
  {
   "r": 0.5,
   "F1": 0.3,
   "F2": 0.0225,
   "ell": 30.0,
   "L": 1024,
   "branches": [
    1,
    1
   ],
   "m_red": 0.22440277537589084,
   "E_b0": 0.8618975453057469,
   "gap_to_next": 9.59675365588275e-05,
   "weight_pi_valley": 0.49999999986346394,
   "dE": -0.04693369496661104,
   "T": 0.005528448192450895,
   "U_exp": -0.05297891874622328,
   "virial_2T": 0.01105689638490179,
   "virial_rUp": 0.011840173045768837,
   "rms_size": 10.262832842569932,
   "dE_deep_asym": -0.0447356151085497,
   "dE_shallow_asym": -1.842032138095726,
   "I1": 0.5049062445957544,
   "I2": 0.40392499567660356,
   "I_pair": 0.8432940210181585,
   "D_num": -0.06553721925419942,
   "D_formula": -0.07347764910719537,
   "D_general": -0.07457782554994448,
   "D_deep_asym": -0.08213396619156932,
   "D_shallow_asym": -0.008646265587440027,
   "rel_dev_D": 0.10806592139892478,
   "D_num_lin": -0.07063049277508511,
   "dev_over_dE": 0.1691839915575543,
   "dev_lin_over_dE": 0.06066337487674359,
   "expected_scale": 0.10983391538639795,
   "T_over_dE": -0.1177927328411682,
   "P": [
    0.013632468604085369,
    0.027264937208170737,
    0.04089740581225611
   ],
   "E_b_P_num": [
    0.8620077257264205,
    0.8623381587525126,
    0.8628885200932579
   ],
   "E_b_P_formula": [
    0.8620087821533566,
    0.8623424926961859,
    0.8628986769342346
   ],
   "max_EP_dev_over_shift": 0.010249343479491474,
   "quartic_coeff": -0.2613189881694042,
   "max_transfer_over_quantum": 0.33171097239229597,
   "D_over_I": -0.07211153880951458,
   "strength": "deep",
   "grip": "square"
  },
  {
   "r": 1.0,
   "F1": 0.15,
   "F2": 0.005625,
   "ell": 30.0,
   "L": 1024,
   "branches": [
    1,
    1
   ],
   "m_red": 0.2233124987990565,
   "E_b0": 0.8573279457483991,
   "gap_to_next": 9.72445210023487e-05,
   "weight_pi_valley": 0.49999999981014853,
   "dE": -0.04708767438777994,
   "T": 0.005698109946961936,
   "U_exp": -0.05310587878660787,
   "virial_2T": 0.011396219893923872,
   "virial_rUp": 0.011728780633827604,
   "rms_size": 10.158754985069127,
   "dE_deep_asym": -0.0446518213276928,
   "dE_shallow_asym": -1.832261971615496,
   "I1": 0.5024531222978772,
   "I2": 0.4019624978383018,
   "I_pair": 0.8648231488057178,
   "D_num": -0.03959247133046123,
   "D_formula": -0.042484346587957536,
   "D_general": -0.043137001779727256,
   "D_deep_asym": -0.04708767438777994,
   "D_shallow_asym": -0.009046990273949722,
   "rel_dev_D": 0.068069194650531,
   "D_num_lin": -0.04140505438656711,
   "dev_over_dE": 0.06141469705385992,
   "dev_lin_over_dE": 0.022920906912967495,
   "expected_scale": 0.1120361188789229,
   "T_over_dE": -0.12101064707584482,
   "P": [
    0.013566234302042686,
    0.027132468604085372,
    0.04069870290612806
   ],
   "E_b_P_num": [
    0.857434342527708,
    0.8577534363535444,
    0.8582849380407014
   ],
   "E_b_P_formula": [
    0.8574347075743045,
    0.8577549930520206,
    0.8582888021815476
   ],
   "max_EP_dev_over_shift": 0.004037797250061421,
   "quartic_coeff": -0.2375890572293668,
   "max_transfer_over_quantum": 0.6646285852304364,
   "D_over_I": -0.04377685485407665,
   "strength": "deep",
   "grip": "square"
  },
  {
   "r": 1.4,
   "F1": 0.10714285714285714,
   "F2": 0.002869897959183673,
   "ell": 29.999999999999996,
   "L": 1024,
   "branches": [
    1,
    1
   ],
   "m_red": 0.2230009912056753,
   "E_b0": 0.8557794394073392,
   "gap_to_next": 0.00010032928036673994,
   "weight_pi_valley": 0.4999999996477008,
   "dE": -0.04737457497564568,
   "T": 0.005972712405184229,
   "U_exp": -0.053360834697098,
   "virial_2T": 0.011945424810368458,
   "virial_rUp": 0.01154867743047816,
   "rms_size": 9.978341296151653,
   "dE_deep_asym": -0.04462783583579728,
   "dE_shallow_asym": -1.8294720007751788,
   "I1": 0.5017522302127694,
   "I2": 0.40140178417021555,
   "I_pair": 0.9037144667049654,
   "D_num": 0.0005604523219804314,
   "D_formula": -0.002817068487178634,
   "D_general": -0.0020399855152443682,
   "D_deep_asym": -0.0018949829990258394,
   "D_shallow_asym": -0.009208813823113872,
   "rel_dev_D": 1.1989487740646798,
   "D_num_lin": 0.0005601047489174638,
   "dev_over_dE": 0.07129395484593543,
   "dev_lin_over_dE": 0.07128661814553977,
   "expected_scale": 0.11478854819909934,
   "T_over_dE": -0.12607421614346262,
   "P": [
    0.013547310215744773,
    0.027094620431489545,
    0.040641930647234314
   ],
   "E_b_P_num": [
    0.8558809731472368,
    0.8561854775143813,
    0.8566926623012874
   ],
   "E_b_P_formula": [
    0.8558813621457753,
    0.8561871303610833,
    0.8566967440532633
   ],
   "max_EP_dev_over_shift": 0.004469611967702649,
   "quartic_coeff": -0.23976131535738196,
   "max_transfer_over_quantum": 0.9309629959354799,
   "D_over_I": 0.0006205501089017693,
   "strength": "deep",
   "grip": "square"
  },
  {
   "r": 0.5,
   "F1": 0.03333333333333333,
   "F2": 0.0002777777777777778,
   "ell": 30.0,
   "L": 1024,
   "branches": [
    1,
    1
   ],
   "m_red": 0.22224914263152679,
   "E_b0": 0.8999571716008505,
   "gap_to_next": 1.639294534561131e-07,
   "weight_pi_valley": 0.5000000036555611,
   "dE": -0.00015185605683298853,
   "T": 8.978240627243822e-05,
   "U_exp": -0.00024168392029264154,
   "virial_2T": 0.00017956481254487644,
   "virial_rUp": 0.00018071627482328689,
   "rms_size": 93.36927724447472,
   "dE_deep_asym": 0.00039118015422424967,
   "dE_shallow_asym": -0.0002778145004853527,
   "I1": 0.5000605709209353,
   "I2": 0.4000484567367482,
   "I_pair": 0.899982587464329,
   "D_num": -0.000126440193354449,
   "D_formula": -0.00012656796643167894,
   "D_general": -0.00012693538207890742,
   "D_deep_asym": -0.0002657480994577299,
   "D_shallow_asym": -3.0341775787887272e-05,
   "rel_dev_D": 0.0010095214518510102,
   "D_num_lin": -0.00012645795716754112,
   "dev_over_dE": 0.0008414091600604606,
   "dev_lin_over_dE": 0.0007244311911694731,
   "expected_scale": 0.0008854570382630171,
   "T_over_dE": -0.591233620475086,
   "P": [
    0.013501635414865253,
    0.027003270829730506,
    0.04050490624459576
   ],
   "E_b_P_num": [
    0.9000584407983854,
    0.9003621608954486,
    0.9008680697025722
   ],
   "E_b_P_formula": [
    0.9000584481073607,
    0.9003622776268912,
    0.9008686601594422
   ],
   "max_EP_dev_over_shift": 0.0006482139648203618,
   "quartic_coeff": -0.21953558609447288,
   "max_transfer_over_quantum": 0.33290331283867036,
   "D_over_I": -0.00014047208667985377,
   "strength": "shallow",
   "grip": "square"
  },
  {
   "r": 1.0,
   "F1": 0.016666666666666666,
   "F2": 6.944444444444444e-05,
   "ell": 30.0,
   "L": 1024,
   "branches": [
    1,
    1
   ],
   "m_red": 0.2222356824268745,
   "E_b0": 0.8999026426643069,
   "gap_to_next": 1.639776399109394e-07,
   "weight_pi_valley": 0.5000000004376254,
   "dE": -0.0001518711645347448,
   "T": 8.981186810629223e-05,
   "U_exp": -0.00024170709594243092,
   "virial_2T": 0.00017962373621258445,
   "virial_rUp": 0.00018071404696043262,
   "rms_size": 93.3664117135673,
   "dE_deep_asym": 0.000391211986731761,
   "dE_shallow_asym": -0.00027779613903641044,
   "I1": 0.5000302854604676,
   "I2": 0.4000242283683741,
   "I_pair": 0.8999734398808037,
   "D_num": -8.107394803796852e-05,
   "D_formula": -8.00129658060076e-05,
   "D_general": -8.112733985524965e-05,
   "D_deep_asym": -0.0001518711645347448,
   "D_shallow_asym": -3.035951409911759e-05,
   "rel_dev_D": 0.013260128796291404,
   "D_num_lin": -8.108125157016852e-05,
   "dev_over_dE": 0.006986067666045922,
   "dev_lin_over_dE": 0.007034157981428532,
   "expected_scale": 0.000885614869205544,
   "T_over_dE": -0.591368798556524,
   "P": [
    0.013500817707432623,
    0.027001635414865246,
    0.04050245312229787
   ],
   "E_b_P_num": [
    0.900003900625608,
    0.900307587031375,
    0.9008134397442658
   ],
   "E_b_P_formula": [
    0.9000039077993716,
    0.9003077032045655,
    0.9008140288798889
   ],
   "max_EP_dev_over_shift": 0.0006468352129557022,
   "quartic_coeff": -0.21954508359787697,
   "max_transfer_over_quantum": 0.6658215450122636,
   "D_over_I": -9.007670845744562e-05,
   "strength": "shallow",
   "grip": "square"
  },
  {
   "r": 1.4,
   "F1": 0.011904761904761906,
   "F2": 3.543083900226758e-05,
   "ell": 29.999999999999996,
   "L": 1024,
   "branches": [
    1,
    1
   ],
   "m_red": 0.2222318366541167,
   "E_b0": 0.899887041428078,
   "gap_to_next": 1.6405983616074593e-07,
   "weight_pi_valley": 0.499999996042689,
   "dE": -0.000151897021094749,
   "T": 8.985260003391622e-05,
   "U_exp": -0.00024174582460171219,
   "virial_2T": 0.00017970520006783245,
   "virial_rUp": 0.00018071548793573475,
   "rms_size": 93.35791785532854,
   "dE_deep_asym": 0.00039122108210056547,
   "dE_shallow_asym": -0.0002777908929430961,
   "I1": 0.5000216324717626,
   "I2": 0.4000173059774101,
   "I_pair": 0.900016306558192,
   "D_num": -2.2631890980673308e-05,
   "D_formula": -2.044607676639802e-05,
   "D_general": -2.242608912640842e-05,
   "D_deep_asym": -6.0758808437899995e-06,
   "D_shallow_asym": -3.036888908582351e-05,
   "rel_dev_D": 0.10690629010390638,
   "D_num_lin": -2.2632460084260243e-05,
   "dev_over_dE": 0.014390105865946105,
   "dev_lin_over_dE": 0.014393852506813935,
   "expected_scale": 0.0008858017854117481,
   "T_over_dE": -0.5915362881137001,
   "P": [
    0.01350058407673759,
    0.02700116815347518,
    0.04050175223021277
   ],
   "E_b_P_num": [
    0.899988291061647,
    0.9002919524840576,
    0.9007977635574365
   ],
   "E_b_P_formula": [
    0.8999882981088894,
    0.9002920681513237,
    0.9007983515553808
   ],
   "max_EP_dev_over_shift": 0.0006456392409711133,
   "quartic_coeff": -0.21956066236469862,
   "max_transfer_over_quantum": 0.932156130798505,
   "D_over_I": -2.5145457617277724e-05,
   "strength": "shallow",
   "grip": "square"
  },
  {
   "r": 0.5,
   "F1": 0.1264124742118579,
   "F2": 0.003995028409090909,
   "ell": 8.0,
   "L": 1024,
   "branches": [
    1,
    1
   ],
   "m_red": 0.22251581990111746,
   "E_b0": 0.9009555077046186,
   "gap_to_next": 3.518550482195515e-06,
   "weight_pi_valley": 0.5000000000839894,
   "dE": -0.00023356289490694682,
   "T": 0.00018890663520083084,
   "U_exp": -0.00042271882719686145,
   "virial_2T": 0.0003778132704016617,
   "virial_rUp": 0.0003743160168396989,
   "rms_size": 70.43815714038287,
   "dE_deep_asym": 0.0040477360369665615,
   "dE_shallow_asym": -0.0002909778755615992,
   "I1": 0.5006605947775142,
   "I2": 0.4005284758220114,
   "I_pair": 0.9010783872925595,
   "D_num": -0.0001106833069661195,
   "D_formula": -0.00011553018737303592,
   "D_general": -0.00011509214388844184,
   "D_deep_asym": -0.00040873506608715693,
   "D_shallow_asym": -4.621852345483579e-05,
   "rel_dev_D": 0.04195336748884784,
   "D_num_lin": -0.00011069690266942545,
   "dev_over_dE": 0.020751928121320203,
   "dev_lin_over_dE": 0.020693718090521562,
   "expected_scale": 0.0015155409760032662,
   "T_over_dE": -0.8088041350750197,
   "P": [
    0.013517836058992884,
    0.027035672117985767,
    0.04055350817697865
   ],
   "E_b_P_num": [
    0.901056896619534,
    0.9013609756968225,
    0.9018674822317514
   ],
   "E_b_P_formula": [
    0.9010569044738735,
    0.9013610947816384,
    0.9018680786279131
   ],
   "max_EP_dev_over_shift": 0.0006539614253848226,
   "quartic_coeff": -0.2189143056325244,
   "max_transfer_over_quantum": 0.3178633341909536,
   "D_over_I": -0.00012281918476051454,
   "strength": "v-shallow",
   "grip": "square"
  },
  {
   "r": 1.0,
   "F1": 0.06320623710592894,
   "F2": 0.0009987571022727273,
   "ell": 8.0,
   "L": 1024,
   "branches": [
    1,
    1
   ],
   "m_red": 0.22236902106166984,
   "E_b0": 0.9003609799724661,
   "gap_to_next": 3.5195335787996385e-06,
   "weight_pi_valley": 0.4999999999162469,
   "dE": -0.00023355532729674167,
   "T": 0.00018905380796325254,
   "U_exp": -0.00042271741111813126,
   "virial_2T": 0.0003781076159265051,
   "virial_rUp": 0.00037419997124650686,
   "rms_size": 70.4627185503051,
   "dE_deep_asym": 0.004049702608680365,
   "dE_shallow_asym": -0.00029076268231335276,
   "I1": 0.5003302973887571,
   "I2": 0.40026423791100574,
   "I_pair": 0.9005158973759179,
   "D_num": -7.863792384499302e-05,
   "D_formula": -8.211239419273771e-05,
   "D_general": -7.831302537087651e-05,
   "D_deep_asym": -0.00023355532729674167,
   "D_shallow_asym": -4.6464127236072584e-05,
   "rel_dev_D": 0.04231359202106894,
   "D_num_lin": -7.864479093428843e-05,
   "dev_over_dE": 0.014876433725402616,
   "dev_lin_over_dE": 0.014847031316239474,
   "expected_scale": 0.001516486144173158,
   "T_over_dE": -0.809460482667783,
   "P": [
    0.013508918029496442,
    0.027017836058992883,
    0.040526754088489325
   ],
   "E_b_P_num": [
    0.9004622984035786,
    0.900766166116787,
    0.9012723206690287
   ],
   "E_b_P_formula": [
    0.9004623060961741,
    0.9007662844672977,
    0.9012729150858372
   ],
   "max_EP_dev_over_shift": 0.0006522443370339234,
   "quartic_coeff": -0.21927429856651987,
   "max_transfer_over_quantum": 0.6358827365937227,
   "D_over_I": -8.73177892633097e-05,
   "strength": "v-shallow",
   "grip": "square"
  },
  {
   "r": 1.4,
   "F1": 0.045147312218520676,
   "F2": 0.0005095699501391466,
   "ell": 8.0,
   "L": 1024,
   "branches": [
    1,
    1
   ],
   "m_red": 0.22232707853611333,
   "E_b0": 0.9001909624967261,
   "gap_to_next": 3.5258174653218788e-06,
   "weight_pi_valley": 0.4999999999381641,
   "dE": -0.0002337055745329497,
   "T": 0.00018937559549962567,
   "U_exp": -0.0004229805202176989,
   "virial_2T": 0.00037875119099925133,
   "virial_rUp": 0.00037424013016275587,
   "rms_size": 70.44606782779596,
   "dE_deep_asym": 0.004050264743786341,
   "dE_shallow_asym": -0.00029070120433157605,
   "I1": 0.500235926706255,
   "I2": 0.40018874136500404,
   "I_pair": 0.9003892650707421,
   "D_num": -3.5403000516964767e-05,
   "D_formula": -3.950531249676491e-05,
   "D_general": -3.0659607264638686e-05,
   "D_deep_asym": -9.348222981318049e-06,
   "D_shallow_asym": -4.656463364318203e-05,
   "rel_dev_D": 0.10384203340085164,
   "D_num_lin": -3.540439255093125e-05,
   "dev_over_dE": 0.01755333388173745,
   "dev_lin_over_dE": 0.01754737752417403,
   "expected_scale": 0.0015177708832092482,
   "T_over_dE": -0.8103169805773116,
   "P": [
    0.013506370021068886,
    0.027012740042137773,
    0.04051911006320666
   ],
   "E_b_P_num": [
    0.9002922569550054,
    0.900596052770992,
    0.901102087565316
   ],
   "E_b_P_formula": [
    0.9002922647164269,
    0.9005961713755288,
    0.9011026824740322
   ],
   "max_EP_dev_over_shift": 0.000652938588562903,
   "quartic_coeff": -0.21938646862240738,
   "max_transfer_over_quantum": 0.8902982620204775,
   "D_over_I": -3.931811485440445e-05,
   "strength": "v-shallow",
   "grip": "square"
  },
  {
   "r": 0.5,
   "F1": 0.3,
   "F2": 0.0,
   "ell": 30.0,
   "L": 1024,
   "branches": [
    1,
    1
   ],
   "m_red": 0.22222222222222224,
   "E_b0": 0.8531745289388937,
   "gap_to_next": 9.478055590994572e-05,
   "weight_pi_valley": 0.4999999998738922,
   "dE": -0.04682547106110635,
   "T": 0.005530243666606173,
   "U_exp": -0.05287676904662079,
   "virial_2T": 0.011060487333212346,
   "virial_rUp": 0.011853004832551818,
   "rms_size": 10.311469814420844,
   "dE_deep_asym": -0.044567785091673036,
   "dE_shallow_asym": -1.8225,
   "I1": 0.5,
   "I2": 0.4,
   "I_pair": 0.8346130826192849,
   "D_num": -0.06538691738071511,
   "D_formula": -0.07337269667369654,
   "D_general": -0.07448267110662088,
   "D_deep_asym": -0.08194457435693611,
   "D_shallow_asym": -0.009365094212221266,
   "rel_dev_D": 0.10883856877301144,
   "D_num_lin": -0.07050958925536949,
   "dev_over_dE": 0.1705434907971377,
   "dev_lin_over_dE": 0.06114423098041569,
   "expected_scale": 0.1108095580243452,
   "T_over_dE": -0.11810332157447621,
   "P": [
    0.0135,
    0.027,
    0.0405
   ],
   "E_b_P_num": [
    0.8532837023531206,
    0.8536111156759874,
    0.8541564485602384
   ],
   "E_b_P_formula": [
    0.8532847660428864,
    0.8536154773548645,
    0.8541666628748281
   ],
   "max_EP_dev_over_shift": 0.01040239380867793,
   "quartic_coeff": -0.2684236796899897,
   "max_transfer_over_quantum": 0.33291823226568923,
   "D_over_I": -0.0726521304230168,
   "strength": "deep",
   "grip": "linear"
  },
  {
   "r": 1.0,
   "F1": 0.15,
   "F2": 0.0,
   "ell": 30.0,
   "L": 1024,
   "branches": [
    1,
    1
   ],
   "m_red": 0.22222222222222224,
   "E_b0": 0.852965736411823,
   "gap_to_next": 9.664675904330622e-05,
   "weight_pi_valley": 0.4999999998131306,
   "dE": -0.04703426358817697,
   "T": 0.00569964159298269,
   "U_exp": -0.05305544818100998,
   "virial_2T": 0.01139928318596538,
   "virial_rUp": 0.011734786751458033,
   "rms_size": 10.1824460472382,
   "dE_deep_asym": -0.044567785091673036,
   "dE_shallow_asym": -1.8225,
   "I1": 0.5,
   "I2": 0.4,
   "I_pair": 0.8604321676767906,
   "D_num": -0.03956783232320937,
   "D_formula": -0.04247455031379082,
   "D_general": -0.04313159687913379,
   "D_deep_asym": -0.04703426358817697,
   "D_shallow_asym": -0.009406852717635389,
   "rel_dev_D": 0.06843434407444876,
   "D_num_lin": -0.04138739859882264,
   "dev_over_dE": 0.06180001064823963,
   "dev_lin_over_dE": 0.023114037130188093,
   "expected_scale": 0.11254064521782219,
   "T_over_dE": -0.12118062786924151,
   "P": [
    0.0135,
    0.027,
    0.0405
   ],
   "E_b_P_num": [
    0.8530716344989586,
    0.853389232864403,
    0.8539182441691642
   ],
   "E_b_P_formula": [
    0.8530720014795634,
    0.8533907966827845,
    0.8539221220214864
   ],
   "max_EP_dev_over_shift": 0.004071202877211847,
   "quartic_coeff": -0.2407394521618218,
   "max_transfer_over_quantum": 0.6658364645313785,
   "D_over_I": -0.043964258136899304,
   "strength": "deep",
   "grip": "linear"
  },
  {
   "r": 1.4,
   "F1": 0.10714285714285714,
   "F2": 0.0,
   "ell": 29.999999999999996,
   "L": 1024,
   "branches": [
    1,
    1
   ],
   "m_red": 0.22222222222222224,
   "E_b0": 0.8526626298632534,
   "gap_to_next": 9.990123518632821e-05,
   "weight_pi_valley": 0.49999999965817127,
   "dE": -0.04733737013674666,
   "T": 0.005974809943829078,
   "U_exp": -0.05332574480969842,
   "virial_2T": 0.011949619887658155,
   "virial_rUp": 0.011552333205192283,
   "rms_size": 9.994444505282106,
   "dE_deep_asym": -0.04456778509167302,
   "dE_shallow_asym": -1.8224999999999993,
   "I1": 0.5,
   "I2": 0.4,
   "I_pair": 0.9005502975324838,
   "D_num": 0.0005502975324838122,
   "D_formula": -0.0028494643964825294,
   "D_general": -0.002071325088014323,
   "D_deep_asym": -0.0018934948054698786,
   "D_shallow_asym": -0.009467474027349328,
   "rel_dev_D": 1.1931231473406432,
   "D_num_lin": 0.0005499612632325924,
   "dev_over_dE": 0.0718198311216956,
   "dev_lin_over_dE": 0.07181272744757411,
   "expected_scale": 0.11516657952569045,
   "T_over_dE": -0.12621761467883072,
   "P": [
    0.0135,
    0.027,
    0.0405
   ],
   "E_b_P_num": [
    0.8527638099509272,
    0.85306725376015,
    0.8535726722776062
   ],
   "E_b_P_formula": [
    0.8527642014461524,
    0.8530689161948497,
    0.8535767741093451
   ],
   "max_EP_dev_over_shift": 0.004507297323872314,
   "quartic_coeff": -0.2421386677645948,
   "max_transfer_over_quantum": 0.9321710503439297,
   "D_over_I": 0.0006114417027597913,
   "strength": "deep",
   "grip": "linear"
  },
  {
   "r": 0.5,
   "F1": 0.03333333333333333,
   "F2": 0.0,
   "ell": 30.0,
   "L": 1024,
   "branches": [
    1,
    1
   ],
   "m_red": 0.22222222222222224,
   "E_b0": 0.899848156379562,
   "gap_to_next": 1.6390056067905334e-07,
   "weight_pi_valley": 0.5000000015010844,
   "dE": -0.00015184362043807464,
   "T": 8.977757480107558e-05,
   "U_exp": -0.00024166665584447004,
   "virial_2T": 0.00017955514960215117,
   "virial_rUp": 0.00018070749213300137,
   "rms_size": 93.37875106329811,
   "dE_deep_asym": 0.0003912438212355896,
   "dE_shallow_asym": -0.00027777777777777783,
   "I1": 0.5,
   "I2": 0.4,
   "I_pair": 0.8998735564712375,
   "D_num": -0.00012644352876256537,
   "D_formula": -0.00012657109482496344,
   "D_general": -0.0001269387365169857,
   "D_deep_asym": -0.00026572633576663063,
   "D_shallow_asym": -3.0368724087614918e-05,
   "rel_dev_D": 0.0010078609383483584,
   "D_num_lin": -0.00012646129566525223,
   "dev_over_dE": 0.0008401147313929538,
   "dev_lin_over_dE": 0.000723106834481636,
   "expected_scale": 0.0008854918602034809,
   "T_over_dE": -0.5912502253440998,
   "P": [
    0.0135,
    0.027,
    0.0405
   ],
   "E_b_P_num": [
    0.899949413313131,
    0.9002530966335309,
    0.9007589441969326
   ],
   "E_b_P_formula": [
    0.899949420620813,
    0.9002532133445659,
    0.9007595345508208
   ],
   "max_EP_dev_over_shift": 0.0006481793859774486,
   "quartic_coeff": -0.21960373901835506,
   "max_transfer_over_quantum": 0.33291823226568923,
   "D_over_I": -0.00014049280973618374,
   "strength": "shallow",
   "grip": "linear"
  },
  {
   "r": 1.0,
   "F1": 0.016666666666666666,
   "F2": 0.0,
   "ell": 30.0,
   "L": 1024,
   "branches": [
    1,
    1
   ],
   "m_red": 0.22222222222222224,
   "E_b0": 0.8998481350533298,
   "gap_to_next": 1.6396319146849692e-07,
   "weight_pi_valley": 0.4999999964930027,
   "dE": -0.00015186494667018824,
   "T": 8.980945324474055e-05,
   "U_exp": -0.00024169846451107543,
   "virial_2T": 0.0001796189064894811,
   "virial_rUp": 0.0001807096559493376,
   "rms_size": 93.37114808737516,
   "dE_deep_asym": 0.0003912438212355896,
   "dE_shallow_asym": -0.00027777777777777783,
   "I1": 0.5,
   "I2": 0.4,
   "I_pair": 0.8999189212201704,
   "D_num": -8.10787798296575e-05,
   "D_formula": -8.00173840743958e-05,
   "D_general": -8.113219813039896e-05,
   "D_deep_asym": -0.00015186494667018824,
   "D_shallow_asym": -3.037298933403764e-05,
   "rel_dev_D": 0.013264564538561681,
   "D_num_lin": -8.108608467491952e-05,
   "dev_over_dE": 0.006989076666696366,
   "dev_lin_over_dE": 0.007037177597307315,
   "expected_scale": 0.0008856322913935422,
   "T_over_dE": -0.5913771098197116,
   "P": [
    0.0135,
    0.027,
    0.0405
   ],
   "E_b_P_num": [
    0.899949386883004,
    0.9002530549015064,
    0.9007588769943133
   ],
   "E_b_P_formula": [
    0.899949394056086,
    0.9002531710643543,
    0.900759466078135
   ],
   "max_EP_dev_over_shift": 0.0006468174958270325,
   "quartic_coeff": -0.21957915912942547,
   "max_transfer_over_quantum": 0.6658364645313785,
   "D_over_I": -9.008753314406388e-05,
   "strength": "shallow",
   "grip": "linear"
  },
  {
   "r": 1.4,
   "F1": 0.011904761904761906,
   "F2": 0.0,
   "ell": 29.999999999999996,
   "L": 1024,
   "branches": [
    1,
    1
   ],
   "m_red": 0.22222222222222224,
   "E_b0": 0.8998481074198512,
   "gap_to_next": 1.6404951419524139e-07,
   "weight_pi_valley": 0.4999999970709411,
   "dE": -0.0001518925801488047,
   "T": 8.985087606726593e-05,
   "U_exp": -0.00024173966019562631,
   "virial_2T": 0.00017970175213453187,
   "virial_rUp": 0.00018071235192073924,
   "rms_size": 93.36130025370238,
   "dE_deep_asym": 0.00039124382123558983,
   "dE_shallow_asym": -0.00027777777777777783,
   "I1": 0.5,
   "I2": 0.4,
   "I_pair": 0.8999773617590303,
   "D_num": -2.26382409697079e-05,
   "D_formula": -2.0451843376714736e-05,
   "D_general": -2.2432467116863396e-05,
   "D_deep_asym": -6.075703205952227e-06,
   "D_shallow_asym": -3.0378516029760927e-05,
   "rel_dev_D": 0.10690467126706382,
   "D_num_lin": -2.2638810417208167e-05,
   "dev_over_dE": 0.01439436732756278,
   "dev_lin_over_dE": 0.014398116342160516,
   "expected_scale": 0.0008858142403654313,
   "T_over_dE": -0.591542233196919,
   "P": [
    0.0135,
    0.027,
    0.0405
   ],
   "E_b_P_num": [
    0.8999493526740125,
    0.9002530009636385,
    0.900758790167155
   ],
   "E_b_P_formula": [
    0.8999493597207359,
    0.9002531166233899,
    0.9007593781278133
   ],
   "max_EP_dev_over_shift": 0.000645626218364461,
   "quartic_coeff": -0.21958501031893798,
   "max_transfer_over_quantum": 0.9321710503439297,
   "D_over_I": -2.5153601077453223e-05,
   "strength": "shallow",
   "grip": "linear"
  },
  {
   "r": 0.5,
   "F1": 0.1264124742118579,
   "F2": 0.0,
   "ell": 8.0,
   "L": 1024,
   "branches": [
    1,
    1
   ],
   "m_red": 0.22222222222222224,
   "E_b0": 0.8997667238551352,
   "gap_to_next": 3.5112944269943824e-06,
   "weight_pi_valley": 0.499999999949246,
   "dE": -0.0002332761448647913,
   "T": 0.00018871399967667653,
   "U_exp": -0.0004222396067343977,
   "virial_2T": 0.00037742799935335306,
   "virial_rUp": 0.0003739423499263032,
   "rms_size": 70.52952591788096,
   "dE_deep_asym": 0.00405167058270851,
   "dE_shallow_asym": -0.0002905475206611571,
   "I1": 0.5,
   "I2": 0.4,
   "I_pair": 0.8998891131310601,
   "D_num": -0.00011088686893989674,
   "D_formula": -0.00011572655401453617,
   "D_general": -0.00011529170049540097,
   "D_deep_asym": -0.0004082332535133848,
   "D_shallow_asym": -4.665522897295826e-05,
   "rel_dev_D": 0.041820005061513614,
   "D_num_lin": -0.00011090053273188572,
   "dev_over_dE": 0.020746592316349153,
   "dev_lin_over_dE": 0.020688018851853274,
   "expected_scale": 0.0015156251104270018,
   "T_over_dE": -0.8089725581930233,
   "P": [
    0.0135,
    0.027,
    0.0405
   ],
   "E_b_P_num": [
    0.8998679790363803,
    0.9001716570789454,
    0.9006774957763861
   ],
   "E_b_P_formula": [
    0.8998679868760469,
    0.9001717759387817,
    0.9006780910433398
   ],
   "max_EP_dev_over_shift": 0.0006535850961582819,
   "quartic_coeff": -0.21965606937048449,
   "max_transfer_over_quantum": 0.31801940920764343,
   "D_over_I": -0.00012320763215544083,
   "strength": "v-shallow",
   "grip": "linear"
  },
  {
   "r": 1.0,
   "F1": 0.06320623710592894,
   "F2": 0.0,
   "ell": 8.0,
   "L": 1024,
   "branches": [
    1,
    1
   ],
   "m_red": 0.22222222222222224,
   "E_b0": 0.8997665880582952,
   "gap_to_next": 3.515903727491221e-06,
   "weight_pi_valley": 0.5000000000364732,
   "dE": -0.0002334119417047864,
   "T": 0.00018895748570474895,
   "U_exp": -0.0004224777895483311,
   "virial_2T": 0.0003779149714094979,
   "virial_rUp": 0.0003740131344727559,
   "rms_size": 70.50840800831219,
   "dE_deep_asym": 0.00405167058270851,
   "dE_shallow_asym": -0.0002905475206611571,
   "I1": 0.5,
   "I2": 0.4,
   "I_pair": 0.899921224270779,
   "D_num": -7.87757292209923e-05,
   "D_formula": -8.224595314098725e-05,
   "D_general": -7.845247834304108e-05,
   "D_deep_asym": -0.0002334119417047864,
   "D_shallow_asym": -4.6682388340957276e-05,
   "rel_dev_D": 0.042193248268960315,
   "D_num_lin": -7.87826249528556e-05,
   "dev_over_dE": 0.014867379512158826,
   "dev_lin_over_dE": 0.014837836328494199,
   "expected_scale": 0.0015165287325784869,
   "T_over_dE": -0.8095450657950383,
   "P": [
    0.0135,
    0.027,
    0.0405
   ],
   "E_b_P_num": [
    0.8998678396266231,
    0.9001715068345995,
    0.9006773274882196
   ],
   "E_b_P_formula": [
    0.8998678473118106,
    0.9001716250723567,
    0.9006779213399334
   ],
   "max_EP_dev_over_shift": 0.0006520544672729157,
   "quartic_coeff": -0.21964560507345132,
   "max_transfer_over_quantum": 0.6360388184152869,
   "D_over_I": -8.752858802332477e-05,
   "strength": "v-shallow",
   "grip": "linear"
  },
  {
   "r": 1.4,
   "F1": 0.045147312218520676,
   "F2": 0.0,
   "ell": 8.0,
   "L": 1024,
   "branches": [
    1,
    1
   ],
   "m_red": 0.22222222222222224,
   "E_b0": 0.8997663968366161,
   "gap_to_next": 3.5232230739534742e-06,
   "weight_pi_valley": 0.4999999999138649,
   "dE": -0.00023360316338394416,
   "T": 0.00018930683479706397,
   "U_exp": -0.000422809387859243,
   "virial_2T": 0.00037861366959412795,
   "virial_rUp": 0.0003741066979832965,
   "rms_size": 70.47867938118075,
   "dE_deep_asym": 0.00405167058270851,
   "dE_shallow_asym": -0.0002905475206611571,
   "I1": 0.5,
   "I2": 0.4,
   "I_pair": 0.899964465580154,
   "D_num": -3.553441984605055e-05,
   "D_formula": -3.9633220102888e-05,
   "D_general": -3.079553133278771e-05,
   "D_deep_asym": -9.344126535357827e-06,
   "D_shallow_asym": -4.672063267678882e-05,
   "rel_dev_D": 0.10341830025915008,
   "D_num_lin": -3.553582289582247e-05,
   "dev_over_dE": 0.017545996370352108,
   "dev_lin_over_dE": 0.017539990245470916,
   "expected_scale": 0.0015178020093314184,
   "T_over_dE": -0.8103778735475603,
   "P": [
    0.0135,
    0.027,
    0.0405
   ],
   "E_b_P_num": [
    0.8998676435394727,
    0.9001712961485586,
    0.9006770924624055
   ],
   "E_b_P_formula": [
    0.8998676512955497,
    0.9001714146723505,
    0.9006776869670187
   ],
   "max_EP_dev_over_shift": 0.000652802754696084,
   "quartic_coeff": -0.2196518047679511,
   "max_transfer_over_quantum": 0.8904543457814018,
   "D_over_I": -3.948268871783394e-05,
   "strength": "v-shallow",
   "grip": "linear"
  },
  {
   "r": 0.5,
   "F1": 0.09999999999999998,
   "F2": 0.0,
   "ell": 30.0,
   "L": 1024,
   "branches": [
    1,
    -1
   ],
   "m_red": -2.0000000000000004,
   "E_b0": 0.10521653408384209,
   "gap_to_next": 1.0793825338584195e-05,
   "weight_pi_valley": 0.4999999992703591,
   "dE": 0.005216534083842109,
   "T": -0.000634429852143054,
   "U_exp": 0.005895870045078229,
   "virial_2T": -0.001268859704286108,
   "virial_rUp": -0.00130277572088047,
   "rms_size": 10.192265055722883,
   "dE_deep_asym": 0.0049519761212970005,
   "dE_shallow_asym": 0.20249999999999987,
   "I1": 0.5,
   "I2": -0.4,
   "I_pair": 0.1091059972008151,
   "D_num": 0.009105997200815119,
   "D_formula": 0.009160656139330843,
   "D_general": 0.0092477208343923,
   "D_deep_asym": 0.00912893464672369,
   "D_shallow_asym": 0.009389761350915796,
   "rel_dev_D": 0.0059667056250532905,
   "D_num_lin": 0.008346009783545685,
   "dev_over_dE": 0.01047801809347455,
   "dev_lin_over_dE": 0.15616620972696707,
   "expected_scale": 0.11232960354207369,
   "T_over_dE": -0.12161903707447463,
   "P": [
    0.0014999999999999996,
    0.002999999999999999,
    0.004499999999999999
   ],
   "E_b_P_num": [
    0.10522684426630122,
    0.10525776412194315,
    0.10530926161716546
   ],
   "E_b_P_formula": [
    0.10522683999477912,
    0.10525775772759022,
    0.10530928728227539
   ],
   "max_EP_dev_over_shift": 0.0004143013095355112,
   "quartic_coeff": -176.10978507687184,
   "max_transfer_over_quantum": 0.3329182322656888,
   "D_over_I": 0.0910599720081512,
   "strength": "deep",
   "grip": "linear"
  },
  {
   "r": 1.0,
   "F1": 0.04999999999999999,
   "F2": 0.0,
   "ell": 30.0,
   "L": 1024,
   "branches": [
    1,
    -1
   ],
   "m_red": -2.0000000000000004,
   "E_b0": 0.10524086728926331,
   "gap_to_next": 1.104125421383384e-05,
   "weight_pi_valley": 0.49999999893655184,
   "dE": 0.005240867289263332,
   "T": -0.0006561156081344888,
   "U_exp": 0.005917248025808297,
   "virial_2T": -0.0013122312162689776,
   "virial_rUp": -0.0012884518264522974,
   "rms_size": 10.05619918156441,
   "dE_deep_asym": 0.0049519761212970005,
   "dE_shallow_asym": 0.20249999999999987,
   "I1": 0.5,
   "I2": -0.4,
   "I_pair": 0.10587354142730789,
   "D_num": 0.005873541427307893,
   "D_formula": 0.005765759775770923,
   "D_general": 0.005762245514364719,
   "D_deep_asym": 0.005240867289263332,
   "D_shallow_asym": 0.009433561120673999,
   "rel_dev_D": 0.018693399608824216,
   "D_num_lin": 0.005547695248619441,
   "dev_over_dE": 0.02056561358799845,
   "dev_lin_over_dE": 0.041608481023402026,
   "expected_scale": 0.11421206096108616,
   "T_over_dE": -0.12519218135491347,
   "P": [
    0.0014999999999999996,
    0.002999999999999999,
    0.004499999999999999
   ],
   "E_b_P_num": [
    0.10525149215812432,
    0.10528335458682603,
    0.10533641809573116
   ],
   "E_b_P_formula": [
    0.10525150400196112,
    0.10528341414005453,
    0.10533659770354356
   ],
   "max_EP_dev_over_shift": 0.0018797100625309036,
   "quartic_coeff": -200.6071184796449,
   "max_transfer_over_quantum": 0.6658364645313776,
   "D_over_I": 0.05873541427307895,
   "strength": "deep",
   "grip": "linear"
  },
  {
   "r": 1.4,
   "F1": 0.035714285714285705,
   "F2": 0.0,
   "ell": 29.999999999999996,
   "L": 1024,
   "branches": [
    1,
    -1
   ],
   "m_red": -2.0000000000000004,
   "E_b0": 0.10527675015136677,
   "gap_to_next": 1.150308363484509e-05,
   "weight_pi_valley": 0.4999999980247235,
   "dE": 0.005276750151366791,
   "T": -0.000693275318023022,
   "U_exp": 0.005950386722903798,
   "virial_2T": -0.001386550636046044,
   "virial_rUp": -0.0012657221845162694,
   "rms_size": 9.854889327650914,
   "dE_deep_asym": 0.004951976121297002,
   "dE_shallow_asym": 0.20249999999999982,
   "I1": 0.5,
   "I2": -0.4,
   "I_pair": 0.10111115005388334,
   "D_num": 0.0011111500538833474,
   "D_formula": 0.0014312345657751916,
   "D_general": 0.0011936252509173917,
   "D_deep_asym": 0.00021107000605467298,
   "D_shallow_asym": 0.009498150272460225,
   "rel_dev_D": 0.22364224533557056,
   "D_num_lin": 0.001098939190476233,
   "dev_over_dE": 0.060659402607670464,
   "dev_lin_over_dE": 0.06297349045659989,
   "expected_scale": 0.11712164499857652,
   "T_over_dE": -0.1313830100224565,
   "P": [
    0.0014999999999999996,
    0.002999999999999999,
    0.004499999999999999
   ],
   "E_b_P_num": [
    0.10528787542623215,
    0.105321238124094,
    0.10537679891705432
   ],
   "E_b_P_formula": [
    0.10528784140944734,
    0.10532111518368906,
    0.1053765714740919
   ],
   "max_EP_dev_over_shift": 0.0030576129774092992,
   "quartic_coeff": -216.2214284694246,
   "max_transfer_over_quantum": 0.9321710503439284,
   "D_over_I": 0.011111500538833476,
   "strength": "deep",
   "grip": "linear"
  },
  {
   "r": 0.5,
   "F1": 0.01111111111111111,
   "F2": 0.0,
   "ell": 30.0,
   "L": 1024,
   "branches": [
    1,
    -1
   ],
   "m_red": -2.0000000000000004,
   "E_b0": 0.10001687934133133,
   "gap_to_next": 1.8237506838048034e-08,
   "weight_pi_valley": 0.5000000001671265,
   "dE": 1.6879341331349496e-05,
   "T": -9.992617514748371e-06,
   "U_exp": 2.686912866799518e-05,
   "virial_2T": -1.9985235029496742e-05,
   "virial_rUp": -2.0081647943445517e-05,
   "rms_size": 93.3407287909861,
   "dE_deep_asym": -4.347153569284329e-05,
   "dE_shallow_asym": 3.08641975308642e-05,
   "I1": 0.5,
   "I2": -0.4,
   "I_pair": 0.10003006959877929,
   "D_num": 3.0069598779325357e-05,
   "D_formula": 3.0038478205599034e-05,
   "D_general": 3.0057628622406527e-05,
   "D_deep_asym": 2.9538847329861617e-05,
   "D_shallow_asym": 3.0382814396429093e-05,
   "rel_dev_D": 0.001036023646514916,
   "D_num_lin": 3.006055968961972e-05,
   "dev_over_dE": 0.0018437078269472656,
   "dev_lin_over_dE": 0.0013081958346132447,
   "expected_scale": 0.0008861542155032919,
   "T_over_dE": -0.5920028109265959,
   "P": [
    0.0014999999999999996,
    0.002999999999999999,
    0.004499999999999999
   ],
   "E_b_P_num": [
    0.10002812481278413,
    0.10006184747466067,
    0.1001180061319825
   ],
   "E_b_P_formula": [
    0.10002812596301733,
    0.10006186582807532,
    0.10011809893650529
   ],
   "max_EP_dev_over_shift": 0.0009177046180465777,
   "quartic_coeff": -226.54967986939957,
   "max_transfer_over_quantum": 0.3329182322656888,
   "D_over_I": 0.0003006959877932536,
   "strength": "shallow",
   "grip": "linear"
  },
  {
   "r": 1.0,
   "F1": 0.005555555555555555,
   "F2": 0.0,
   "ell": 30.0,
   "L": 1024,
   "branches": [
    1,
    -1
   ],
   "m_red": -2.0000000000000004,
   "E_b0": 0.10001688173572928,
   "gap_to_next": 1.824500049651956e-08,
   "weight_pi_valley": 0.500000001772491,
   "dE": 1.688173572930718e-05,
   "T": -9.996232598830308e-06,
   "U_exp": 2.6872705858050452e-05,
   "virial_2T": -1.9992465197660617e-05,
   "virial_rUp": -2.0081879159477375e-05,
   "rms_size": 93.33305136630152,
   "dE_deep_asym": -4.347153569284329e-05,
   "dE_shallow_asym": 3.08641975308642e-05,
   "I1": 0.5,
   "I2": -0.4,
   "I_pair": 0.10002501743781948,
   "D_num": 2.501743781951893e-05,
   "D_formula": 2.4878721808371426e-05,
   "D_general": 2.496287330010115e-05,
   "D_deep_asym": 1.688173572930718e-05,
   "D_shallow_asym": 3.038712431275292e-05,
   "rel_dev_D": 0.005575688824207425,
   "D_num_lin": 2.501118066293272e-05,
   "dev_over_dE": 0.008216928245517261,
   "dev_lin_over_dE": 0.007846281726312063,
   "expected_scale": 0.0008862961821976129,
   "T_over_dE": -0.5921329867447552,
   "P": [
    0.0014999999999999996,
    0.002999999999999999,
    0.004499999999999999
   ],
   "E_b_P_num": [
    0.10002812777496882,
    0.10006185213698804,
    0.10011801361718106
   ],
   "E_b_P_formula": [
    0.10002812893756922,
    0.10006187054308906,
    0.10011810655228877
   ],
   "max_EP_dev_over_shift": 0.0009189496563554751,
   "quartic_coeff": -226.60271216982906,
   "max_transfer_over_quantum": 0.6658364645313776,
   "D_over_I": 0.00025017437819518934,
   "strength": "shallow",
   "grip": "linear"
  },
  {
   "r": 1.4,
   "F1": 0.003968253968253968,
   "F2": 0.0,
   "ell": 29.999999999999996,
   "L": 1024,
   "branches": [
    1,
    -1
   ],
   "m_red": -2.0000000000000004,
   "E_b0": 0.10001688484159843,
   "gap_to_next": 1.825566756097352e-08,
   "weight_pi_valley": 0.5000000004373143,
   "dE": 1.6884841598452383e-05,
   "T": -1.000094502687185e-05,
   "U_exp": 2.6877343240542837e-05,
   "virial_2T": -2.00018900537437e-05,
   "virial_rUp": -2.0082165442448026e-05,
   "rms_size": 93.32309706918805,
   "dE_deep_asym": -4.3471535692843306e-05,
   "dE_shallow_asym": 3.0864197530864185e-05,
   "I1": 0.5,
   "I2": -0.4,
   "I_pair": 0.10001849825007054,
   "D_num": 1.849825007055328e-05,
   "D_formula": 1.8277056911232554e-05,
   "D_general": 1.8434058937701775e-05,
   "D_deep_asym": 6.753936639380997e-07,
   "D_shallow_asym": 3.039271487721429e-05,
   "rel_dev_D": 0.012102230703499459,
   "D_num_lin": 1.849482885086572e-05,
   "dev_over_dE": 0.013100102718227428,
   "dev_lin_over_dE": 0.012897481943397581,
   "expected_scale": 0.00088648030864014,
   "T_over_dE": -0.5923031595267384,
   "P": [
    0.0014999999999999996,
    0.002999999999999999,
    0.004499999999999999
   ],
   "E_b_P_num": [
    0.10002813161372072,
    0.10006185817190671,
    0.10011802330410911
   ],
   "E_b_P_formula": [
    0.10002813278580526,
    0.10006187661842578,
    0.10011811633945995
   ],
   "max_EP_dev_over_shift": 0.0009198810079805646,
   "quartic_coeff": -226.6435680344599,
   "max_transfer_over_quantum": 0.9321710503439284,
   "D_over_I": 0.00018498250070553285,
   "strength": "shallow",
   "grip": "linear"
  }
 ]
}
output/mutual_term_check.json9 KB · raw
{
 "max_rel_dev_all": 3.897751476969538e-14,
 "cases": [
  {
   "case": {
    "L": 16,
    "kappa_h": 1.0,
    "m_h": 0.3,
    "tau": 0.8,
    "M1": 0.7,
    "M2": 0.5,
    "F0": 1.0,
    "F1": 0.3,
    "F2": 0.2,
    "sigma1": 1,
    "sigma2": 1
   },
   "abar": 1.1044079964957103,
   "eigvec_residual": 1.2401028427486198e-15,
   "rows": [
    {
     "labels_n1_n2_m": [
      0,
      0,
      1
     ],
     "numerical": [
      -0.009095531629472589,
      -5.011369276867195e-17
     ],
     "formula": -0.009095531629472589,
     "rel_dev": 5.509704634117998e-15
    },
    {
     "labels_n1_n2_m": [
      1,
      -2,
      3
     ],
     "numerical": [
      -0.002055022405216982,
      -3.7776884649305124e-19
     ],
     "formula": -0.0020550224052169764,
     "rel_dev": 2.7496018845798712e-15
    },
    {
     "labels_n1_n2_m": [
      3,
      -5,
      2
     ],
     "numerical": [
      -0.0019647772798085067,
      1.4353025475214594e-17
     ],
     "formula": -0.0019647772798085297,
     "rel_dev": 1.3792099839515786e-14
    },
    {
     "labels_n1_n2_m": [
      -2,
      4,
      -1
     ],
     "numerical": [
      -0.00497183203244238,
      -1.8622863212809886e-17
     ],
     "formula": -0.004971832032442402,
     "rel_dev": 5.7490609418608994e-15
    },
    {
     "labels_n1_n2_m": [
      5,
      2,
      7
     ],
     "numerical": [
      -0.0007791912219067409,
      -2.31840925424854e-19
     ],
     "formula": -0.0007791912219067381,
     "rel_dev": 3.4913156748031588e-15
    },
    {
     "labels_n1_n2_m": [
      2,
      2,
      0
     ],
     "numerical": [
      -0.013263873307308011,
      6.47154058509828e-35
     ],
     "formula": -0.013263873307308263,
     "rel_dev": 1.8963910329122623e-14
    },
    {
     "labels_n1_n2_m": [
      -3,
      1,
      0
     ],
     "numerical": [
      -0.014873993450272818,
      -4.942551449792622e-35
     ],
     "formula": -0.014873993450272793,
     "rel_dev": 1.6327914050029304e-15
    },
    {
     "labels_n1_n2_m": [
      6,
      -6,
      8
     ],
     "numerical": [
      -0.0007598802139332297,
      3.766671451786338e-19
     ],
     "formula": -0.0007598802139332235,
     "rel_dev": 8.290311254847195e-15
    }
   ],
   "max_rel_dev": 1.8963910329122623e-14
  },
  {
   "case": {
    "L": 16,
    "kappa_h": 1.0,
    "m_h": 0.3,
    "tau": 0.8,
    "M1": 0.7,
    "M2": 0.5,
    "F0": 1.0,
    "F1": 0.3,
    "F2": 0.0225,
    "sigma1": 1,
    "sigma2": 1
   },
   "abar": 1.0117458996057673,
   "eigvec_residual": 1.2401028427486198e-15,
   "rows": [
    {
     "labels_n1_n2_m": [
      0,
      0,
      1
     ],
     "numerical": [
      -0.00786464976668061,
      -4.3709330242262885e-17
     ],
     "formula": -0.007864649766680611,
     "rel_dev": 5.5620709887315565e-15
    },
    {
     "labels_n1_n2_m": [
      1,
      -2,
      3
     ],
     "numerical": [
      -0.0014901074050707327,
      -1.7207426720664638e-19
     ],
     "formula": -0.001490107405070737,
     "rel_dev": 2.912690106592727e-15
    },
    {
     "labels_n1_n2_m": [
      3,
      -5,
      2
     ],
     "numerical": [
      -0.0014263056881132026,
      1.013924362726558e-17
     ],
     "formula": -0.0014263056881131915,
     "rel_dev": 1.0519079363062898e-14
    },
    {
     "labels_n1_n2_m": [
      -2,
      4,
      -1
     ],
     "numerical": [
      -0.00401346513145179,
      -1.494365719580665e-17
     ],
     "formula": -0.004013465131451822,
     "rel_dev": 8.820568839304457e-15
    },
    {
     "labels_n1_n2_m": [
      5,
      2,
      7
     ],
     "numerical": [
      -0.0004995978675767499,
      3.733667855684837e-19
     ],
     "formula": -0.0004995978675767546,
     "rel_dev": 9.577859537479806e-15
    },
    {
     "labels_n1_n2_m": [
      2,
      2,
      0
     ],
     "numerical": [
      -0.011404426886681245,
      2.337933810742404e-36
     ],
     "formula": -0.011404426886681016,
     "rel_dev": 2.0078474885604593e-14
    },
    {
     "labels_n1_n2_m": [
      -3,
      1,
      0
     ],
     "numerical": [
      -0.012960383691982028,
      -3.6687004238147716e-36
     ],
     "formula": -0.012960383691982028,
     "rel_dev": 2.8307035586333926e-34
    },
    {
     "labels_n1_n2_m": [
      6,
      -6,
      8
     ],
     "numerical": [
      -0.0004848820280437704,
      -1.301175741193455e-19
     ],
     "formula": -0.00048488202804377825,
     "rel_dev": 1.6213309975128075e-14
    }
   ],
   "max_rel_dev": 2.0078474885604593e-14
  },
  {
   "case": {
    "L": 16,
    "kappa_h": 1.0,
    "m_h": 0.3,
    "tau": 0.8,
    "M1": 0.7,
    "M2": 0.5,
    "F0": 1.0,
    "F1": 0.3,
    "F2": 0.2,
    "sigma1": 1,
    "sigma2": -1
   },
   "abar": 1.1044079964957103,
   "eigvec_residual": 1.2401028427486198e-15,
   "rows": [
    {
     "labels_n1_n2_m": [
      0,
      0,
      1
     ],
     "numerical": [
      0.009095531629472586,
      5.010945651449388e-17
     ],
     "formula": 0.009095531629472589,
     "rel_dev": 5.522428226449429e-15
    },
    {
     "labels_n1_n2_m": [
      1,
      -2,
      3
     ],
     "numerical": [
      0.0020550224052169807,
      2.2527785353840697e-19
     ],
     "formula": 0.0020550224052169764,
     "rel_dev": 2.1131914798815703e-15
    },
    {
     "labels_n1_n2_m": [
      3,
      -5,
      2
     ],
     "numerical": [
      0.0019647772798085067,
      -1.4431741308232534e-17
     ],
     "formula": 0.0019647772798085297,
     "rel_dev": 1.3813361791940207e-14
    },
    {
     "labels_n1_n2_m": [
      -2,
      4,
      -1
     ],
     "numerical": [
      0.004971832032442384,
      1.8499102246367158e-17
     ],
     "formula": 0.004971832032442402,
     "rel_dev": 5.221673751335779e-15
    },
    {
     "labels_n1_n2_m": [
      5,
      2,
      7
     ],
     "numerical": [
      0.0007791912219067406,
      9.53273701951327e-19
     ],
     "formula": 0.0007791912219067381,
     "rel_dev": 3.4261963570479904e-15
    },
    {
     "labels_n1_n2_m": [
      2,
      2,
      0
     ],
     "numerical": [
      0.013263873307308013,
      -3.4407581456318445e-35
     ],
     "formula": 0.013263873307308263,
     "rel_dev": 1.8833124740645916e-14
    },
    {
     "labels_n1_n2_m": [
      -3,
      1,
      0
     ],
     "numerical": [
      0.01487399345027281,
      7.665628986769642e-36
     ],
     "formula": 0.014873993450272793,
     "rel_dev": 1.1662795750020932e-15
    },
    {
     "labels_n1_n2_m": [
      6,
      -6,
      8
     ],
     "numerical": [
      0.0007598802139332297,
      -4.082649082797704e-19
     ],
     "formula": 0.0007598802139332235,
     "rel_dev": 8.292901431948084e-15
    }
   ],
   "max_rel_dev": 1.8833124740645916e-14
  },
  {
   "case": {
    "L": 16,
    "kappa_h": 1.0,
    "m_h": 0.3,
    "tau": 1.3,
    "M1": 0.7,
    "M2": 0.5,
    "F0": 1.0,
    "F1": 0.3,
    "F2": 0.2,
    "sigma1": 1,
    "sigma2": 1
   },
   "abar": 1.1044079964957103,
   "eigvec_residual": 1.898408963776082e-15,
   "rows": [
    {
     "labels_n1_n2_m": [
      0,
      0,
      1
     ],
     "numerical": [
      -0.009452768437474534,
      -5.180710986354804e-17
     ],
     "formula": -0.009452768437474474,
     "rel_dev": 8.443487036947806e-15
    },
    {
     "labels_n1_n2_m": [
      1,
      -2,
      3
     ],
     "numerical": [
      -0.002154974744517884,
      9.392057234842114e-20
     ],
     "formula": -0.0021549747445178684,
     "rel_dev": 7.245000815492042e-15
    },
    {
     "labels_n1_n2_m": [
      3,
      -5,
      2
     ],
     "numerical": [
      -0.0014044877488271026,
      1.0317629824255665e-17
     ],
     "formula": -0.0014044877488270796,
     "rel_dev": 1.7938637053735863e-14
    },
    {
     "labels_n1_n2_m": [
      -2,
      4,
      -1
     ],
     "numerical": [
      -0.0035197995011773575,
      -1.2775932186492922e-17
     ],
     "formula": -0.0035197995011773432,
     "rel_dev": 5.45043427745903e-15
    },
    {
     "labels_n1_n2_m": [
      5,
      2,
      7
     ],
     "numerical": [
      -0.0007897085393569411,
      -1.8911538743628514e-17
     ],
     "formula": -0.0007897085393569168,
     "rel_dev": 3.897751476969538e-14
    },
    {
     "labels_n1_n2_m": [
      2,
      2,
      0
     ],
     "numerical": [
      -0.0077984672498935224,
      2.8061683132107806e-35
     ],
     "formula": -0.007798467249893521,
     "rel_dev": 2.22444157343931e-16
    },
    {
     "labels_n1_n2_m": [
      -3,
      1,
      0
     ],
     "numerical": [
      -0.00946145479080685,
      6.923458992192756e-35
     ],
     "formula": -0.009461454790806904,
     "rel_dev": 5.683737749033285e-15
    },
    {
     "labels_n1_n2_m": [
      6,
      -6,
      8
     ],
     "numerical": [
      -0.0007598802139332206,
      3.041983685234327e-19
     ],
     "formula": -0.0007598802139332235,
     "rel_dev": 3.7312349346042454e-15
    }
   ],
   "max_rel_dev": 3.897751476969538e-14
  }
 ]
}
Analyse with ChatGPTnDot Physics Research plugin
ChatGPT plugin

Analyse this research in ChatGPT

The free plugin nDot Physics Research gives ChatGPT read-only access to everything published here.

  1. Open ChatGPT and select the plugin in the message box.
  2. Write your question after it:nDot Physics Research Explain what 08-ilang-propagation/07-bond-inertia establishes and what it builds on.
  3. Name the subproject when you cite an equation number, e.g. (17.7) of 03-ilang-space.
Open in ChatGPT

Interested in this research?

Join the readers' list and get an e-mail when a new stage or a correction is published. No password, and you can unsubscribe at any time.