Contents · Previous · Next

only b=τb0b=\tau b_0, retaining the associated κ\kappa and all other coefficients. Use ρ=(1+R2)−8\rho=(1+R^2)^{-8}, Q=π5/2520Q=\pi^5/2520, dμr=210R9(1+R2)−8dRd\mu_r=210R^9(1+R^2)^{-8}dR, and FB(z)=21z5−35z6+15z7F_B(z)=21z^5-35z^6+15z^7, z=R2/(1+R2)z=R^2/(1+R^2). Let SS be the preceding smooth step, δ=10−12\delta=10^{-12}, hs=10−4h_s=10^{-4}, β=δ+(1−δ)S(R−1000)\beta=\delta+(1-\delta)S(R-1000), χ=e−1/(1−R2)\chi=e^{-1/(1-R^2)} for R<1R<1 and zero otherwise. For each pinned rational radius R0∈(2.47414,2.47883)R_0\in(2.47414,2.47883) and weight ww, set

B=S(2(R−R0+hs)/hs)S(2(R0+hs−R)/hs),I=∫B dμr,j1=∫χ dμr,r=β+1−w−∫β dμrj1χ+wIB. B=S(2(R-R_0+h_s)/h_s)S(2(R_0+h_s-R)/h_s),\quad I=\int B\,d\mu_r,\quad j_1=\int\chi\,d\mu_r,\quad r=\beta+\frac{1-w-\int\beta\,d\mu_r}{j_1}\chi+\frac wI B.

Exact integrals define this density: its core coefficient is positive, ∫r dμr=1\int r\,d\mu_r=1, r≥δr\ge\delta, and r=1r=1 for R≥1001R\ge1001. The complete energy, including wall and the appropriate pair-charge or affine-charge core compensation, is

Eτ/Q−D=Ewall+∫(H−D)r dμr+12Cd(∫tr dμr)2+τ∫b0(rlog⁡r−r+1) dμr. E_\tau/Q-D=E_{\rm wall}+\int(H-D)r\,d\mu_r+ \tfrac12 C_d(\int tr\,d\mu_r)^2+\tau\int b_0(r\log r-r+1)\,d\mu_r.

Here D=(d+2d)D=\binom{d+2}d, HH includes the actual flattened mass, saturated drift and occupation/phase-edge potential, and tt is the carrier feature, not squared radius. The previous core construction gives

Contents · Previous · Next