Scattering Center Analysis

• Scattering Center Analysis is a set of analytical and computational procedures designed to identify effective points or frames that simplify scattering dynamics across various physical systems.
• It involves methods ranging from geometric shadow analysis in Rutherford scattering to optimal multipole expansion origins and sparse inverse techniques in SAR imaging.
• The approach highlights the coordinate-dependent nature of scattering centers, offering practical insights for modeling relativistic shocks, non-Hermitian transport, and multi-center scattering problems.

Scattering center analysis designates a family of analytical and computational procedures that identify the effective point, frame, origin, support, or central structure that organizes a scattering phenomenon. In the literature, the expression “scattering center” appears in several technically distinct senses: the focus of the projectile shadow in repulsive Rutherford scattering, the local “Weibel frame” in relativistic collisionless shocks, the optimal origin of a multipole expansion, the sparse support of attributed scattering centers in SAR, and the central non-Hermitian cluster in waveguide transport (Žugec et al., 2020, Pelletier et al., 2019, Kildishev et al., 2023, Yang et al., 2024, Jin et al., 2011). Taken together, these uses show that scattering center analysis is less a single formalism than a recurring strategy for compressing scattering dynamics into a geometrically, kinematically, or inversion-theoretically privileged representation.

1. Conceptual scope

Across the cited works, a scattering center is the entity with respect to which scattering becomes simplest. In classical Coulomb scattering, simplicity appears as a paraboloidal shadow whose focus coincides with the relevant physical center. In relativistic shock microphysics, it appears as a local frame in which the transverse current-filamentation instability becomes purely magnetic. In multipole electrodynamics, it appears as the expansion origin that minimizes residual quadrupole content. In SAR, it appears as a sparse coefficient support in a physics-informed dictionary. In non-Hermitian transport, it appears as the finite central graph attached to asymptotic leads (Žugec et al., 2020, Pelletier et al., 2019, Kildishev et al., 2023, Yang et al., 2024, Jin et al., 2011).

A central methodological distinction is therefore required. Some analyses seek a geometric center in real space, some a scattering-center frame in velocity space, some an optimal origin for truncating a field expansion, and some a latent sparse representation in an inverse problem. This distinction matters because the corresponding invariants are different: focal structure and scaling in Rutherford scattering, dielectric-tensor properties in Weibel turbulence, quadrupole norms in long-wave multipole theory, and regularized data-fidelity objectives in imaging.

The literature also makes clear that the preferred center is often not obvious. In the long-wave approximation, the optimal electric and magnetic scattering centers are “not co-local with the centers of mass” (Kildishev et al., 2023). In relativistic shocks, the scattering center is not the background-plasma rest frame but the Weibel frame (Pelletier et al., 2019). In repulsive Rutherford scattering, the scattering center depends on the coordinate frame: it is the target in the fixed-target frame and the center of mass in the CM frame (Žugec et al., 2020). A plausible implication is that scattering center analysis is fundamentally a problem of representation choice under dynamical constraints.

2. Geometric scattering centers in repulsive Rutherford scattering

For repulsive Rutherford scattering, the relevant geometric object is the “shadow” of the scattering process: the portion of space entirely shielded from admitting any particle trajectory. The analysis introduces the natural length scale

$$b_0 \equiv \frac{Z_p Z_t e^2}{4\pi \epsilon_0 \mu v_0^2},$$

described as the only combination of charges, masses, and initial relative speed $v_0$ that carries dimensions of length (Žugec et al., 2020).

In the fixed-target frame, one asks which trajectory just grazes a given polar ray. Minimization with respect to impact parameter yields

$\tilde b_0(\theta)=\frac{2b_0}{\tan(\theta/2)}.$

In cylindrical coordinates $(\rho,z)$, the shadow surface is

$$\rho(\theta)=\frac{4b_0}{\tan(\theta/2)}, \qquad z(\theta)=2b_0\left[\frac{1}{\tan^2(\theta/2)}-1\right],$$

and eliminating $\theta$ gives the single-valued paraboloid

$z=\frac{\rho^2}{8b_0}-2b_0.$

When expressed in scaled coordinates $\bar\rho=\rho/b_0$ and $\bar z=z/b_0$, the shadow becomes

$\bar z=\frac{\bar\rho^2}{8}-2,$

which is independent of all physical parameters. The paper interprets this as a manifestation of self-similarity and as the revelation of a natural length scale for Rutherford scattering (Žugec et al., 2020).

The focal structure identifies the scattering center. Since a paraboloid of the form $v_0$0 has focal distance $v_0$1, the fixed-target shadow has $v_0$2. Because the vertex lies at $v_0$3 and the target sits at $v_0$4, the target is exactly at the focus. The conclusion is explicit: in the fixed-target frame, the scattering center is the target itself (Žugec et al., 2020).

In the center-of-mass frame, the projectile and target cast distinct shadows. With $v_0$5 and $v_0$6, the two bodies obey

$v_0$7

The projectile shadow is

$v_0$8

while the target shadow is

$v_0$9

Their focal distances are $\tilde b_0(\theta)=\frac{2b_0}{\tan(\theta/2)}.$0 and $\tilde b_0(\theta)=\frac{2b_0}{\tan(\theta/2)}.$1, but because their vertices lie at $\tilde b_0(\theta)=\frac{2b_0}{\tan(\theta/2)}.$2 and $\tilde b_0(\theta)=\frac{2b_0}{\tan(\theta/2)}.$3, both foci coincide at $\tilde b_0(\theta)=\frac{2b_0}{\tan(\theta/2)}.$4, the origin of the CM frame. The common focus is therefore the CM itself (Žugec et al., 2020).

This analysis resolves a common ambiguity: the scattering center is not an invariant point independent of coordinates. The paper states that the notion is coordinate-dependent but always sits at the focus of the shadow. In the target-rest frame it is the target; in the CM frame it is the CM. The same scattering dynamics therefore admits different effective centers, each selected by the geometry of forbidden trajectories (Žugec et al., 2020).

3. Scattering-center frames in relativistic collisionless shocks

In unmagnetized, relativistic collisionless pair shocks, the term “scattering center” refers not to a point in configuration space but to a local frame, the “Weibel frame” $\tilde b_0(\theta)=\frac{2b_0}{\tan(\theta/2)}.$5. This is defined as the frame in which the linear transverse current-filamentation instability is purely magnetic: $\tilde b_0(\theta)=\frac{2b_0}{\tan(\theta/2)}.$6 so that $\tilde b_0(\theta)=\frac{2b_0}{\tan(\theta/2)}.$7 and only $\tilde b_0(\theta)=\frac{2b_0}{\tan(\theta/2)}.$8 remains. In this frame, particles see almost magnetostatic filaments and scatter as off magnetic fluctuations (Pelletier et al., 2019).

The frame velocity is first derived kinetically. In the local plasma rest frame $\tilde b_0(\theta)=\frac{2b_0}{\tan(\theta/2)}.$9, the ratio of fields of the fastest-growing mode yields

$(\rho,z)$0

In the far precursor, where $(\rho,z)$1, the explicit scaling becomes

$(\rho,z)$2

and the corresponding four-velocity is sub-relativistic relative to the background plasma (Pelletier et al., 2019).

A quasistatic nonlinear model gives the same scaling. In $(\rho,z)$3, vanishing electrostatic contribution implies

$(\rho,z)$4

To first order in $(\rho,z)$5, the result again is $(\rho,z)$6, with the leading nonlinear correction multiplying $(\rho,z)$7 by $(\rho,z)$8 (Pelletier et al., 2019).

Transformation to the shock frame $(\rho,z)$9 gives

$$\rho(\theta)=\frac{4b_0}{\tan(\theta/2)}, \qquad z(\theta)=2b_0\left[\frac{1}{\tan^2(\theta/2)}-1\right],$$0

for $$\rho(\theta)=\frac{4b_0}{\tan(\theta/2)}, \qquad z(\theta)=2b_0\left[\frac{1}{\tan^2(\theta/2)}-1\right],$$1 and $$\rho(\theta)=\frac{4b_0}{\tan(\theta/2)}, \qquad z(\theta)=2b_0\left[\frac{1}{\tan^2(\theta/2)}-1\right],$$2. Thus $$\rho(\theta)=\frac{4b_0}{\tan(\theta/2)}, \qquad z(\theta)=2b_0\left[\frac{1}{\tan^2(\theta/2)}-1\right],$$3 moves almost at the upstream speed but slightly more slowly toward the shock than the plasma itself (Pelletier et al., 2019).

The paper compares these predictions with dedicated large-scale 2D3V PIC simulations. Over the well-defined precursor region $$\rho(\theta)=\frac{4b_0}{\tan(\theta/2)}, \qquad z(\theta)=2b_0\left[\frac{1}{\tan^2(\theta/2)}-1\right],$$4, the kinetic-linear estimate and quasistatic-nonlinear estimate bracket and track the PIC data. Far upstream, $$\rho(\theta)=\frac{4b_0}{\tan(\theta/2)}, \qquad z(\theta)=2b_0\left[\frac{1}{\tan^2(\theta/2)}-1\right],$$5 as predicted; near the shock, where $$\rho(\theta)=\frac{4b_0}{\tan(\theta/2)}, \qquad z(\theta)=2b_0\left[\frac{1}{\tan^2(\theta/2)}-1\right],$$6, one finds $$\rho(\theta)=\frac{4b_0}{\tan(\theta/2)}, \qquad z(\theta)=2b_0\left[\frac{1}{\tan^2(\theta/2)}-1\right],$$7–$$\rho(\theta)=\frac{4b_0}{\tan(\theta/2)}, \qquad z(\theta)=2b_0\left[\frac{1}{\tan^2(\theta/2)}-1\right],$$8. Residual discrepancies of factor $$\rho(\theta)=\frac{4b_0}{\tan(\theta/2)}, \qquad z(\theta)=2b_0\left[\frac{1}{\tan^2(\theta/2)}-1\right],$$9 are attributed to the marginally-kinetic nature of the beam and the onset of damping and oblique modes at large $\theta$0 (Pelletier et al., 2019).

The physical significance of the scattering-center frame is operational. In $\theta$1, pitch-angle scattering off $\theta$2 dominates; the effective scattering frequency is evaluated naturally in this frame; the noninertial electric field $\theta$3 leads to plasma heating and deceleration; and the spatial dependence $\theta$4 couples turbulence growth to the deceleration profile and the shock-transition width (Pelletier et al., 2019).

4. Optimal scattering centers in multipole electrodynamics

For multipole decompositions of scattered radiation, the scattering center is the chosen expansion origin. Because multipoles re-mix under translation, the set $\theta$5 is not unique. The problem addressed in “The art of finding the optimal scattering center(s)” is therefore to determine the origin $\theta$6 for which the truncated multipolar spectrum is maximally compact (Kildishev et al., 2023).

The paper defines the optimal electric and magnetic scattering centers separately by minimizing residual poloidal quadrupole norms: $\theta$7 with

$\theta$8

The analysis is carried out within the long-wave approximation, where dipoles are $\theta$9, the electric and magnetic quadrupoles are $z=\frac{\rho^2}{8b_0}-2b_0.$0, and toroidal terms are retained explicitly (Kildishev et al., 2023).

For the electric problem, the shifted poloidal quadrupole is

$z=\frac{\rho^2}{8b_0}-2b_0.$1

Its norm can be written as

$z=\frac{\rho^2}{8b_0}-2b_0.$2

and minimization gives the linear system $z=\frac{\rho^2}{8b_0}-2b_0.$3. The paper further derives a compact closed-form solution and an axial simplification for axisymmetric cases (Kildishev et al., 2023).

For the magnetic problem, the translation law is more involved: $z=\frac{\rho^2}{8b_0}-2b_0.$4 The corresponding norm takes the form $z=\frac{\rho^2}{8b_0}-2b_0.$5, and the stationarity condition becomes a system of three coupled cubic equations. In the axisymmetric case, this reduces to a single real cubic whose unique real root yields the optimal axial displacement (Kildishev et al., 2023).

Two conclusions are emphasized. First, the optimal electric and magnetic scattering centers are not generally the center of mass. Second, they are not generally co-located with one another. For a dielectric cone of height $z=\frac{\rho^2}{8b_0}-2b_0.$6, base radius $z=\frac{\rho^2}{8b_0}-2b_0.$7, and $z=\frac{\rho^2}{8b_0}-2b_0.$8, the optimal electric center lies approximately $z=\frac{\rho^2}{8b_0}-2b_0.$9 above the CM and the magnetic center approximately $\bar\rho=\rho/b_0$0 below over $\bar\rho=\rho/b_0$1. Re-expansion about these points reduces residual quadrupole contributions by more than $\bar\rho=\rho/b_0$2, and dipoles plus toroidal terms reproduce the far field to within $\bar\rho=\rho/b_0$3 (Kildishev et al., 2023).

This directly corrects a widespread simplification that uses the center of mass or geometric center as the default multipole origin. The cited analysis shows that such choices can leave large residual quadrupole content and degrade truncation efficiency. The optimal scattering center is instead defined variationally by the compactness of the translated multipolar spectrum (Kildishev et al., 2023).

5. Attributed scattering centers and scattering-structure inversion

In SAR, scattering center analysis is cast as a sparse inverse problem. The continuous attributed scattering center model writes the complex echo as

$\bar\rho=\rho/b_0$4

where each scatterer is described by amplitude $\bar\rho=\rho/b_0$5, frequency-dependency exponent $\bar\rho=\rho/b_0$6, and location $\bar\rho=\rho/b_0$7. After sampling and discretization on an $\bar\rho=\rho/b_0$8 grid with $\bar\rho=\rho/b_0$9, the measurement is written as

$\bar z=z/b_0$0

with a physics-informed dictionary $\bar z=z/b_0$1 and sparse coefficient vector $\bar z=z/b_0$2. Canonical extraction is then

$\bar z=z/b_0$3

The deep-unfolding approach treats each ISTA iteration as one layer with learnable $\bar z=z/b_0$4: $\bar z=z/b_0$5 initialized by $\bar z=z/b_0$6 and trained with

$\bar z=z/b_0$7

The paper uses AdamW with weight-decay $\bar z=z/b_0$8, OneCycleLR with peak LR $\bar z=z/b_0$9, $\bar z=\frac{\bar\rho^2}{8}-2,$0 epochs, batch size $\bar z=\frac{\bar\rho^2}{8}-2,$1, and an RTX 3090. On D-15, $\bar z=\frac{\bar\rho^2}{8}-2,$2 gives residual $\bar z=\frac{\bar\rho^2}{8}-2,$3 and time $\bar z=\frac{\bar\rho^2}{8}-2,$4; the reported baselines are AMP $\bar z=\frac{\bar\rho^2}{8}-2,$5 and $\bar z=\frac{\bar\rho^2}{8}-2,$6, OMP $\bar z=\frac{\bar\rho^2}{8}-2,$7 and $\bar z=\frac{\bar\rho^2}{8}-2,$8, and ISTA $\bar z=\frac{\bar\rho^2}{8}-2,$9 and $v_0$00. The paper states that the method cuts the ISTA residual by approximately $v_0$01 and speeds inference by nearly two orders of magnitude, while preserving interpretability through the physically meaningful dictionary (Yang et al., 2024).

A related inverse problem appears in VLBI observations of the Galactic Center, where the scattering structure is a stochastic phase screen rather than a sparse target signature. The unknowns are

$v_0$02

including the wavelet coefficients $v_0$03 of the intrinsic image $v_0$04 and the phase-screen realizations $v_0$05. The data-fidelity objective is

$v_0$06

where $v_0$07 applies the scattering forward model before comparison to visibilities and closure quantities. The screen prior is

$v_0$08

and the full multiobjective problem is

$v_0$09

with scalarization

$v_0$10

The scattering screen is parameterized by a power-law spectrum $v_0$11 with $v_0$12, an ensemble-average kernel $v_0$13, and a real-space forward model

$v_0$14

The paper reports that at $v_0$15, for moving screens with $v_0$16 and $v_0$17, the recovered speeds are $v_0$18 and $v_0$19, with errors below $v_0$20, and intrinsic ring $v_0$21. At $v_0$22, MOEA/D yields a cluster of approximately $v_0$23 individuals, ring-morphology clusters containing about $v_0$24–$v_0$25 of solutions depending on the prior, and MO-PSO refinement gives a recovered ring with $v_0$26 and screen with $v_0$27 (Mus et al., 22 Apr 2025).

These two inverse formulations show complementary versions of scattering center analysis: one extracts attributed pointlike support from SAR echoes, while the other reconstructs both intrinsic source and scattering screen in a highly nonconvex, degenerate interferometric inverse problem.

6. Non-Hermitian scattering centers in waveguide transport

In tight-binding transport, the scattering center is the finite central subsystem coupled to semi-infinite leads. “Hermitian scattering behavior for the non-Hermitian scattering center” considers a one-dimensional setup with two Hermitian leads attached to a central non-Hermitian cluster

$v_0$28

where

$v_0$29

with $v_0$30 and $v_0$31 Hermitian and the inter-cluster coupling purely anti-Hermitian, $v_0$32. For an incoming plane wave of energy $v_0$33, Bethe-ansatz matching yields reflection and transmission amplitudes $v_0$34 and $v_0$35, and a key property of the inverse truncated matrix,

$v_0$36

implies $v_0$37 and $v_0$38. The central result is

$v_0$39

so the Dirac flux is conserved exactly, even though $v_0$40 is non-Hermitian. The paper further shows that any parity-symmetric real Hermitian graph with additional $v_0$41-symmetric potentials can be transformed into this structure, and in a four-site example flux conservation holds if and only if the gain and loss strengths are balanced, $v_0$42 (Jin et al., 2011).

A different non-Hermitian scattering-center analysis appears in the flux-controlled triangular-ring model. The central three sites form a non-Hermitian ring threaded by Aharonov-Bohm flux $v_0$43, with onsite potential $v_0$44. The exact amplitudes satisfy

$v_0$45

$v_0$46

$v_0$47

with spectral singularities at $v_0$48. A closed-form solution exists at

$v_0$49

At this point,

$v_0$50

and the scattering matrix becomes

$v_0$51

The paper states that a $v_0$52-symmetric non-Hermitian scattering center always has symmetric transmission although the dynamics within the isolated center can be unidirectional, while the flux-controlled triangular ring realizes perfect unidirectionality at the spectral singularity (Li et al., 2014).

These results show that “scattering center analysis” in non-Hermitian transport can mean either structural criteria for Hermitian behavior despite non-Hermiticity or parameter tuning to produce asymmetric transmission and reflectionless absorption.

7. Multi-center, fixed-center, and internal-resonance generalizations

The term also appears in problems where the center is fixed, multiple, or spatially extended. In Euler’s two-center problem, a particle moves in

$v_0$53

with fixed centers at $v_0$54 and $v_0$55. Together with $v_0$56, the system possesses a Runge-Lenz-type integral $v_0$57, and the resulting Liouville-integrable scattering dynamics carries nontrivial topology identified as scattering monodromy. For small loops $v_0$58 around critical lines $v_0$59, the scattering-monodromy matrices are

$v_0$60

The paper emphasizes that $v_0$61 carries “pure” scattering monodromy, while $v_0$62 carry mixed scattering and Hamiltonian monodromy (Martynchuk et al., 2018).

In the three-body fixed-center approximation with attraction, the heavy pair acts as two static scatterers and the light-particle three-body amplitude is ambiguous up to one real parameter. In coordinate space, the divergent kernel is replaced by

$v_0$63

so that the renormalized multiple-scattering term depends on the undetermined constant $v_0$64. In momentum space, the same freedom appears as a homogeneous solution

$v_0$65

The cited analysis states that coordinate-space, momentum-space, and finite-cutoff treatments are equivalent and that the parameter must be fixed by one three-body datum such as the three-body scattering length (Kudryavtsev et al., 2016).

A different extension concerns a small acoustic scattering center compared with the wavelength. For a fluid sphere of radius $v_0$66, the internal and scattered fields are expanded in partial waves with coefficients

$v_0$67

$v_0$68

The internal stored energy is

$v_0$69

with each term proportional to $v_0$70, so resonances arise when the denominator of $v_0$71 vanishes. In the small-particle regime, the monopole asymptotics yield a resonance condition

$v_0$72

which produces a closed-form small-sphere resonance frequency $v_0$73. The paper interprets this as an internal resonance driven by impedance mismatch, even when the incident wavelength is much larger than the sphere (Rodrigues et al., 2024).

These generalizations suggest that scattering center analysis extends naturally from single effective centers to fixed-center approximations, multi-center integrable scattering, and finite-size resonant inclusions. The unifying theme is the extraction of a reduced structure—topological, renormalized, or resonant—that governs the observable scattering response.