Skip to content
Shadow Theory

Section 2 4 October 2026

Physical model and statistical assumptions

Reading position 3 of 15

2 Physical model and statistical assumptions

2.1 Hamiltonians and a common smooth domain

Let q=(q1,…,qN)∈(R3)Nq=(q_1,\ldots,q_N)\in(\R^3)^N, with distinguishable particles of fixed masses mi>0m_i>0, and set ℏ=1\hbar=1. Initially the wavefunction is scalar. The Hamiltonians are

Hv=−12∑i=1Nmi−1Δi+∑iUi(qi)+∑{i,j}∈EWij(qi,qj)+∑ivi(qi). H_v=-\frac12\sum_{i=1}^N m_i^{-1}\Delta_i+ \sum_i U_i(q_i)+\sum_{\{i,j\}\in E}W_{ij}(q_i,q_j)+\sum_i v_i(q_i). (2.1)

All potentials are real. The graph ({1,…,N},E)(\{1,\ldots,N\},E) is connected, and every edge satisfies

DqiDqjWij≢0. D_{q_i}D_{q_j}W_{ij}\not\equiv0. (2.2)

This condition excludes an interaction that is only a sum of one-body terms. It requires neither full rank nor nonvanishing of the mixed derivative everywhere.

For definiteness we use the following common-domain class throughout. The uncontrolled potential is a positive definite confining quadratic form plus a real smooth bounded perturbation with bounded derivatives of every positive order. Each Ui,WijU_i,W_{ij} is correspondingly quadratic plus such a bounded smooth function. The initial control class is vi∈Cc∞(R3;R)v_i\in\Cc(\R^3;\R). Finite smooth time-dependent schedules and, subsequently, finite sums of translates of Schwartz profiles are also allowed. On every finite time interval their coefficients and all required spatial derivatives are bounded. There are no singular configuration sets or vector potentials in this model.

The Hamiltonians are self-adjoint on the domain of the confining quadratic oscillator, with common invariant core S(R3N)\Sch(\R^{3N}). These statements also hold for the finite-dimensional matrix-valued bounded perturbations used below. Here is a direct justification. A bounded symmetric perturbation preserves the oscillator's self-adjoint domain, and its interaction-picture Dyson series converges in operator norm on finite intervals. Write H(t)=A+B(t)H(t)=A+B(t), with AA the positive confining oscillator. Leibniz' rule and oscillator estimates give

∥[Am,B(t)]ϕ∥2≤Cm,T∥(A+1)mϕ∥2(0≤t≤T, ϕ∈S). \norm{[A^m,B(t)]\phi}_2\le C_{m,T}\norm{(A+1)^m\phi}_2 \quad(0\le t\le T,\ \phi\in\Sch).

Indeed, the commutator is a finite sum of position–derivative monomials of degree at most 2m2m, multiplied by bounded derivatives of BB. The graph seminorms ∥(A+1)mϕ∥2\norm{(A+1)^m\phi}_2 are equivalent to the Schwartz topology: oscillator raising and lowering operators bound the monomials in one direction, and expansion of AmA^m, followed by Sobolev embedding, gives the reverse seminorm comparisons. Spectral projections of AA commute with AmA^m. For the projected equation the symmetric BAmB A^m term cancels in the energy derivative, so the displayed bound and Gronwall give estimates independent of the projection. Strong unitary convergence and weak compactness in each graph domain identify the limit; lower semicontinuity preserves the uniform graph bounds. Applying a higher graph bound and interpolation gives continuity in every lower graph norm. The equation then gives differentiability in those norms, proving preservation of S\Sch and smooth time dependence there. Thus the tangent vectors and differentiations used below are defined on the stated core.

Virtual coordinate conjugations in the proof are evaluated on this common smooth core. We do not require all conjugated unbounded operators to have the original oscillator's full operator domain.

For every fixed admitted control, the lower bound is H(t)≥(inf⁡σ(A)−∥B(t)∥)IH(t)\ge (\inf\sigma(A)-\norm{B(t)})I. It is uniform on a specified finite smooth schedule, but is not asserted to be uniform over all unrestricted control amplitudes. The Schwartz assertion means S=⋂m≥0D((A+1)m)\Sch=\bigcap_{m\ge0}D((A+1)^m), with the equivalent family of graph seminorms just described; it does not identify any single finite graph domain with S\Sch.

2.2 Assignments, regularity, and consistency

Write Ss=S(Rd;Cs)∖{0}\mathscr S_s=\Sch(\R^d;\C^s)\setminus\{0\}, where d=3Nd=3N and initially s=1s=1. Let

rΨ=Ψ†Ψ,JiΨ=mi−1Im⁡(Ψ†∇iΨ),viΨ=JiΨ/rΨ. r^\Psi=\Psi^\dagger\Psi,\qquad J_i^\Psi=m_i^{-1}\operatorname{Im}(\Psi^\dagger\nabla_i\Psi), \qquad v_i^\Psi=J_i^\Psi/r^\Psi.

The velocity is defined on ΩΨ={rΨ>0}\Omega_\Psi=\{r^\Psi>0\}. A probability assignment is a map P:Ψ↦μΨ=pΨ(q) dq\Pcal:\Psi\mapsto \mu^\Psi=p^\Psi(q)\,dq. The unsuperscripted vi(qi)v_i(q_i) in (2.1) is a control potential; viΨv_i^\Psi is a guidance velocity. Derivatives of P\Pcal below denote derivatives of its local density representatives, or the resulting distributions on a nonzero patch. We impose the following standing assumptions.

  1. (A1)

    The density is nonnegative, integrable, and normalized; P(cΨ)=P(Ψ)\Pcal(c\Psi)=\Pcal(\Psi) for c∈C∖{0}c\in\C\setminus\{0\}. The assignment has no control-potential or preparation-history argument.

  2. (A2)

    It has compatible C1C^1 density representatives on compact nonzero patches. More precisely, for each closed ball KK, the set

    UK={Ψ∈Ss:min⁡KrΨ>0} \mathscr U_K=\{\Psi\in\mathscr S_s:\min_K r^\Psi>0\}

    is open in the underlying real Schwartz space, and Ψ↦pΨ∣K\Psi\mapsto p^\Psi|_K is Bastiani C2C^2 into C1(K)C^1(K). This means that its first and second real directional derivatives exist and are jointly continuous in the base point and directions; the second derivative is symmetric [2]. Representatives and derivatives agree on overlaps.

  3. (A3)

    For every state and every admitted Hamiltonian, the same assignment satisfies the infinitesimal consistency identity

    DPΨ[−iHvΨ]=−div⁡(pΨvΨ)on ΩΨ. D\Pcal_\Psi[-iH_v\Psi]=-\diver(p^\Psi v^\Psi) \quad\hbox{on }\Omega_\Psi. (2.3)

    The identity is local, and may equivalently be read against compactly supported test functions there.

  4. (A4)

    For extension to all states we assume the following approximation continuity: whenever normalized nowhere-zero states in Ss\mathscr S_s converge in L2L^2 to a normalized state Ψ\Psi, their assigned measures converge weakly to μΨ\mu^\Psi.

The domain is the full indicated Schwartz state space; in particular it includes all the scalar phases, compact coordinate transports, and internal variations used in the proofs. One may instead use a smaller domain closed under those operations and containing the approximation classes, with the same local calculus. The full-space formulation avoids an implicit richness hypothesis.

If a density assignment is equivariant under each admitted Schrödinger history, differentiation gives (2.3). We use this necessary infinitesimal condition as the definition needed for the theorem. It is a statistical assumption, not a consequence of the equation for the wavefunction. Spatial locality of vi(qi)v_i(q_i), finite-jet locality of pΨ(q)p^\Psi(q) as a functional of Ψ\Psi, and the absence of a history argument in P\Pcal are three different notions. Only the first and third occur here.

Born equilibrium satisfies all four assumptions. The map in (1.1) is smooth on Ss\mathscr S_s, is projective, and obeys (2.3) by the Schrödinger continuity equation and conservation of ∥Ψ∥2\norm{\Psi}_2. For normalized states,

∥∣Ψ∣2−∣Φ∣2∥1≤2∥Ψ−Φ∥2, \norm{|\Psi|^2-|\Phi|^2}_1\le 2\norm{\Psi-\Phi}_2, (2.4)

which proves the stronger continuity in total variation.

Theorem 2.1 (Local-control characterization)

For the scalar model (2.1) in the common-domain class above, suppose the interaction graph is connected and (2.2) holds on every edge. Under A1–A3, with all real compactly supported smooth one-body controls, the only assigned density on a nowhere-zero state is (1.1). Under A4 the same conclusion holds, almost everywhere, for every nonzero Schwartz state. For N=1N=1, the interaction hypotheses are unnecessary.

The proof occupies the next two sections. A profile-control version is given in theorem 5.1, and the precise internal-state extension in theorem 6.1. All of these are characterizations of density-valued assignments. No assertion about arbitrary singular configuration measures is part of the theorem.