General Boundary Conditions (GBC) in Quantum Systems

• General Boundary Conditions (GBC) are self-adjoint boundary prescriptions derived from conservation laws and operator symmetry in quantum and continuum models.
• They encompass a one-parameter Robin family that interpolates between Dirichlet and Neumann limits, thereby modifying quantization, spectral properties, and transition frequencies in systems like the gravitational quantum bouncer.
• GBC principles extend to diverse fields such as obstacle scattering, transport networks, and Casimir physics, providing a unified framework to ensure well-defined operator domains and physical observables.

General Boundary Conditions (GBC) denotes a family of admissible boundary prescriptions determined by the operator structure and conservation laws of a model rather than a single fixed endpoint rule. In the gravitational quantum bouncer relevant to qBounce, GBC is the one-parameter Robin family required by self-adjointness of the linear-potential Hamiltonian on the half-line; the corresponding real parameter $\lambda$ interpolates continuously between the Dirichlet and Neumann limits and modifies the Airy quantization of ultracold neutrons above a mirror (Sung et al., 17 Oct 2025). Across the broader literature, the same phrase is used in related but nonidentical senses—mixed Dirichlet/Neumann data in transport networks, smooth Robin perturbations of Neumann in obstacle scattering, or unitary-matrix parameterizations of self-adjoint extensions in continuum and Casimir models—but the common core is the systematic classification of boundary laws compatible with current conservation, self-adjointness, or an analogous variational principle (Kharitonov, 2022).

1. Foundational principle: conservation laws and self-adjointness

In one-dimensional continuum models with polynomial Hamiltonian

$H(p)=\sum_{n=0}^N h_n p^n,\qquad p=-i\partial_x,$

the boundary problem is controlled by norm conservation, equivalently probability-current conservation at the boundary. The formalism developed for continuum models derives the current as a Hermitian quadratic form in the boundary data, diagonalizes it into positive- and negative-current sectors, and rewrites it in the standardized form

$$j[\psi(0)] = \tilde\psi_+^\dagger[\psi(0)]\,\tilde\psi_+[\psi(0)]-\tilde\psi_-^\dagger[\psi(0)]\,\tilde\psi_-[\psi(0)].$$

When the positive and negative sectors have equal dimension, the most general admissible boundary conditions are

$$\tilde\psi_+[\psi(0)] = U\,\tilde\psi_-[\psi(0)],\qquad U\in U(N/2).$$

This is the universal “standardized” form of GBC in that formalism, and the paper emphasizes that admissibility is tied to self-adjointness rather than mere symmetricity (Kharitonov, 2022).

That distinction is essential. A boundary condition can nullify the quadratic current for a single state and still fail to define a self-adjoint operator. The continuum-model analysis therefore separates current-nullifying relations from admissible relations, and shows that only the unitary family associated with balanced positive and negative current sectors describes a well-defined boundary. This supplies a general mathematical template for later uses of GBC in more specialized settings, including the qBounce half-line Hamiltonian (Kharitonov, 2022).

2. GBC in the qBounce gravitational quantum bouncer

For qBounce, the system is the gravitational quantum bouncer: ultracold neutrons moving above a horizontal mirror in Earth’s gravitational field. The Hamiltonian is

$\hat H=\frac{\hat p^2}{2m}+F_0\hat x,$

with the neutron specialization $F_0=mg_c$, and the potential

$$U(x)= \begin{cases} mg_c x,& x\ge 0,\ \infty,& x<0. \end{cases}$$

The wavefunction lives in $L^2((0,\infty))$, vanishes at $+\infty$, and must satisfy a boundary condition at the mirror $x=0$. The central point of the qBounce analysis is that on a domain with boundary the Hamiltonian may be symmetric without being self-adjoint, whereas self-adjointness is what guarantees a real spectrum, unitary time evolution, and conservation of probability current (Sung et al., 17 Oct 2025).

Using the boundary form obtained by integration by parts and the von Neumann deficiency equations,

$H(p)=\sum_{n=0}^N h_n p^n,\qquad p=-i\partial_x,$0

the paper introduces the dimensionless variables

$H(p)=\sum_{n=0}^N h_n p^n,\qquad p=-i\partial_x,$1

and reduces the deficiency problem to Airy equations with normalizable solutions

$H(p)=\sum_{n=0}^N h_n p^n,\qquad p=-i\partial_x,$2

The deficiency indices are therefore $H(p)=\sum_{n=0}^N h_n p^n,\qquad p=-i\partial_x,$3, so the Hamiltonian has a $H(p)=\sum_{n=0}^N h_n p^n,\qquad p=-i\partial_x,$4 family of self-adjoint extensions. In the convention adopted in the paper, the general self-adjoint boundary condition is

$H(p)=\sum_{n=0}^N h_n p^n,\qquad p=-i\partial_x,$5

or, dimensionally,

$H(p)=\sum_{n=0}^N h_n p^n,\qquad p=-i\partial_x,$6

This is a Robin-type condition with real, dimensionless $H(p)=\sum_{n=0}^N h_n p^n,\qquad p=-i\partial_x,$7. The special cases are

$H(p)=\sum_{n=0}^N h_n p^n,\qquad p=-i\partial_x,$8

and finite real $H(p)=\sum_{n=0}^N h_n p^n,\qquad p=-i\partial_x,$9 gives the general Robin family. The paper notes a sign/convention subtlety: $$j[\psi(0)] = \tilde\psi_+^\dagger[\psi(0)]\,\tilde\psi_+[\psi(0)]-\tilde\psi_-^\dagger[\psi(0)]\,\tilde\psi_-[\psi(0)].$$0 is introduced earlier through $$j[\psi(0)] = \tilde\psi_+^\dagger[\psi(0)]\,\tilde\psi_+[\psi(0)]-\tilde\psi_-^\dagger[\psi(0)]\,\tilde\psi_-[\psi(0)].$$1, but the final adopted convention is the equivalent Robin form $$j[\psi(0)] = \tilde\psi_+^\dagger[\psi(0)]\,\tilde\psi_+[\psi(0)]-\tilde\psi_-^\dagger[\psi(0)]\,\tilde\psi_-[\psi(0)].$$2 (Sung et al., 17 Oct 2025).

Applied to the boundary current, the same condition gives

$$j[\psi(0)] = \tilde\psi_+^\dagger[\psi(0)]\,\tilde\psi_+[\psi(0)]-\tilde\psi_-^\dagger[\psi(0)]\,\tilde\psi_-[\psi(0)].$$3

so no probability current leaks through the mirror. In this setting, GBC is therefore not an optional phenomenological modification of a hard wall; it is the exact self-adjoint-extension structure of the half-line Hamiltonian (Sung et al., 17 Oct 2025).

3. Spectral structure, matrix elements, and sum rules

With

$$j[\psi(0)] = \tilde\psi_+^\dagger[\psi(0)]\,\tilde\psi_+[\psi(0)]-\tilde\psi_-^\dagger[\psi(0)]\,\tilde\psi_-[\psi(0)].$$4

the stationary Schrödinger equation becomes the Airy equation

$$j[\psi(0)] = \tilde\psi_+^\dagger[\psi(0)]\,\tilde\psi_+[\psi(0)]-\tilde\psi_-^\dagger[\psi(0)]\,\tilde\psi_-[\psi(0)].$$5

whose normalizable solution is

$$j[\psi(0)] = \tilde\psi_+^\dagger[\psi(0)]\,\tilde\psi_+[\psi(0)]-\tilde\psi_-^\dagger[\psi(0)]\,\tilde\psi_-[\psi(0)].$$6

Imposing the general boundary condition yields the exact quantization condition

$$j[\psi(0)] = \tilde\psi_+^\dagger[\psi(0)]\,\tilde\psi_+[\psi(0)]-\tilde\psi_-^\dagger[\psi(0)]\,\tilde\psi_-[\psi(0)].$$7

Its roots are denoted $$j[\psi(0)] = \tilde\psi_+^\dagger[\psi(0)]\,\tilde\psi_+[\psi(0)]-\tilde\psi_-^\dagger[\psi(0)]\,\tilde\psi_-[\psi(0)].$$8, and the levels are

$$j[\psi(0)] = \tilde\psi_+^\dagger[\psi(0)]\,\tilde\psi_+[\psi(0)]-\tilde\psi_-^\dagger[\psi(0)]\,\tilde\psi_-[\psi(0)].$$9

In the Dirichlet limit, $$\tilde\psi_+[\psi(0)] = U\,\tilde\psi_-[\psi(0)],\qquad U\in U(N/2).$$0, the zeros of $$\tilde\psi_+[\psi(0)] = U\,\tilde\psi_-[\psi(0)],\qquad U\in U(N/2).$$1; in the Neumann limit, $$\tilde\psi_+[\psi(0)] = U\,\tilde\psi_-[\psi(0)],\qquad U\in U(N/2).$$2, the zeros of $$\tilde\psi_+[\psi(0)] = U\,\tilde\psi_-[\psi(0)],\qquad U\in U(N/2).$$3. The paper also derives perturbative asymptotics for small $$\tilde\psi_+[\psi(0)] = U\,\tilde\psi_-[\psi(0)],\qquad U\in U(N/2).$$4 and large $$\tilde\psi_+[\psi(0)] = U\,\tilde\psi_-[\psi(0)],\qquad U\in U(N/2).$$5, showing that in the Dirichlet regime the transition gaps are uniformly rescaled to the displayed order, whereas in the Neumann regime the transition shifts are state-dependent already at leading orders in $$\tilde\psi_+[\psi(0)] = U\,\tilde\psi_-[\psi(0)],\qquad U\in U(N/2).$$6 (Sung et al., 17 Oct 2025).

A substantial part of the work is the derivation of closed formulas beyond the spectrum itself. The normalization is

$$\tilde\psi_+[\psi(0)] = U\,\tilde\psi_-[\psi(0)],\qquad U\in U(N/2).$$7

and the paper introduces boundary amplitudes $$\tilde\psi_+[\psi(0)] = U\,\tilde\psi_-[\psi(0)],\qquad U\in U(N/2).$$8, in particular $$\tilde\psi_+[\psi(0)] = U\,\tilde\psi_-[\psi(0)],\qquad U\in U(N/2).$$9 and $\hat H=\frac{\hat p^2}{2m}+F_0\hat x,$0, together with identities such as

$\hat H=\frac{\hat p^2}{2m}+F_0\hat x,$1

It then derives a generalized recursion relation for matrix elements $\hat H=\frac{\hat p^2}{2m}+F_0\hat x,$2, explicit diagonal moments $\hat H=\frac{\hat p^2}{2m}+F_0\hat x,$3 and $\hat H=\frac{\hat p^2}{2m}+F_0\hat x,$4, and explicit off-diagonal formulas for $\hat H=\frac{\hat p^2}{2m}+F_0\hat x,$5 and $\hat H=\frac{\hat p^2}{2m}+F_0\hat x,$6 (Sung et al., 17 Oct 2025).

The momentum sector is more delicate. The paper emphasizes that $\hat H=\frac{\hat p^2}{2m}+F_0\hat x,$7 is not self-adjoint on the half-line with these boundary conditions, even though

$\hat H=\frac{\hat p^2}{2m}+F_0\hat x,$8

Using a modified Ehrenfest theorem, it derives

$\hat H=\frac{\hat p^2}{2m}+F_0\hat x,$9

which yields closed expressions for $F_0=mg_c$0 and $F_0=mg_c$1. These formulas feed into generalized sum rules. The Thomas–Reiche–Kuhn dipole sum rule,

$F_0=mg_c$2

and the Bethe sum rule remain unchanged, while the diagonal moments entering closure and monopole identities acquire explicit $F_0=mg_c$3-dependence. The paper’s conclusion is therefore not that GBC destroys standard operator identities, but that it changes the detailed moments and matrix elements on which observable spectroscopy depends (Sung et al., 17 Oct 2025).

4. Experimental relevance for qBounce: inferred $F_0=mg_c$4, spectroscopy, and fifth-force systematics

For qBounce the characteristic scales used are

$F_0=mg_c$5

with $F_0=mg_c$6 and

$F_0=mg_c$7

The measured $F_0=mg_c$8 transition quoted in the paper is

$F_0=mg_c$9

which in the quoted qBounce analysis corresponds to

$$U(x)= \begin{cases} mg_c x,& x\ge 0,\ \infty,& x<0. \end{cases}$$0

slightly above $$U(x)= \begin{cases} mg_c x,& x\ge 0,\ \infty,& x<0. \end{cases}$$1. To test whether boundary physics can account for this, the paper fits $$U(x)= \begin{cases} mg_c x,& x\ge 0,\ \infty,& x<0. \end{cases}$$2 through

$$U(x)= \begin{cases} mg_c x,& x\ge 0,\ \infty,& x<0. \end{cases}$$3

and finds

$$U(x)= \begin{cases} mg_c x,& x\ge 0,\ \infty,& x<0. \end{cases}$$4

For this fitted positive $$U(x)= \begin{cases} mg_c x,& x\ge 0,\ \infty,& x<0. \end{cases}$$5, transition frequencies shift upward relative to the Dirichlet case, which can mimic an increase in the inferred $$U(x)= \begin{cases} mg_c x,& x\ge 0,\ \infty,& x<0. \end{cases}$$6 (Sung et al., 17 Oct 2025).

The paper interprets positive $$U(x)= \begin{cases} mg_c x,& x\ge 0,\ \infty,& x<0. \end{cases}$$7 as making the mirror “more repulsive than the usual Dirichlet mirror.” Its tabulated values show the generalized energies $$U(x)= \begin{cases} mg_c x,& x\ge 0,\ \infty,& x<0. \end{cases}$$8 lower than the Dirichlet values $$U(x)= \begin{cases} mg_c x,& x\ge 0,\ \infty,& x<0. \end{cases}$$9 by about $L^2((0,\infty))$0-$L^2((0,\infty))$1 peV in absolute level position for the first few states, while the spectroscopic gaps shift in the direction needed to account for the observed transition change. The plotted $L^2((0,\infty))$2, $L^2((0,\infty))$3, and $L^2((0,\infty))$4 make explicit that the boundary parameter modifies directly measurable transition frequencies (Sung et al., 17 Oct 2025).

The same mechanism creates a systematic issue for short-range new-physics searches. The paper argues that generalized boundary effects can mimic or mask short-range fifth-force signals, because exotic near-mirror interactions are also expected to perturb neutron wavefunctions and transition frequencies. An unmodeled nonzero $L^2((0,\infty))$5 could therefore be misidentified as new physics, or conversely could hide it. The paper is explicit that the fitted $L^2((0,\infty))$6 is treated phenomenologically: it may “soak up” true boundary physics together with unmodeled systematics in the simplified Hamiltonian, and it is not claimed to be a uniquely determined microscopic mirror parameter (Sung et al., 17 Oct 2025).

To estimate under-mirror penetration, the paper extends the wavefunction for $L^2((0,\infty))$7 as

$L^2((0,\infty))$8

This gives the estimate

$L^2((0,\infty))$9

and for $+\infty$0,

$+\infty$1

The paper explicitly calls this only an estimate. Its broader interpretive claim is narrower: the mathematically rigorous part is the existence of the one-parameter Robin family; the physical meaning of $+\infty$2 may include imperfect reflection, finite penetration, surface interaction details, mirror composition or coating dependence, and omitted systematic shifts (Sung et al., 17 Oct 2025).

5. Field-dependent meanings of “general boundary conditions”

The term GBC is used across several research areas, but its precise content is discipline-specific. In convex-obstacle scattering, “general boundary conditions” does not mean arbitrary local elliptic boundary operators in the full Boutet de Monvel sense; it means Neumann and general smooth Robin conditions

$+\infty$3

handled uniformly through the augmented operator $+\infty$4. In that setting, the Robin term is lower order after semiclassical scaling, so the leading cubic resonance-band structure remains Neumann-like and is governed by the zeros of $+\infty$5 (Jin, 2014).

In transport-network optimization, GBC means the simultaneous allowance of prescribed inflows or outflows and prescribed pressures. The network variables satisfy

$+\infty$6

with vertices partitioned into Neumann vertices $+\infty$7 and Dirichlet vertices $+\infty$8. The paper’s central conclusion is that under these general boundary conditions, minimally dissipative networks need not be trees, while minimizers of the modified objective called complementary dissipation do recover tree-like structure (Chang et al., 2017).

In Casimir and continuum spectral problems, GBC is again tied to self-adjointness. For the scalar field in a slab $+\infty$9, admissible self-adjoint extensions of $x=0$0 are classified by $x=0$1, with

$x=0$2

subject to positivity constraints

$x=0$3

This family includes Dirichlet, Neumann, Robin, periodic, and antiperiodic cases and controls the sign of the thermal Casimir force, while preserving thermodynamic stability through positive entropy for all admissible boundary conditions (Munoz-Castaneda et al., 2020). The related vacuum-energy analysis for homogeneous parallel plates likewise treats the complete four-parameter family of self-adjoint, non-negative boundary conditions parameterized by $x=0$4, and classifies the parameter space into attractive, repulsive, and Casimirless sectors (Asorey et al., 2013).

Other literatures use the term in still different ways. In half-space kinetic equations, GBC means prescribed incoming data plus a reflection operator,

$x=0$5

with the reflection satisfying a weighted contraction property; this covers specular, diffuse, Maxwell, bounce-back, and linearized Cercignani–Lampis type conditions (Li et al., 2015). In zero-number theory for one-dimensional parabolic equations, the “general boundary conditions” of interest are those for which the boundary trace may be zero at some times and nonzero at others, subject to the structural condition $x=0$6; the result is an extension of the zero-number diminishing property to boundary-touching nodal events (Lou, 2018).

6. Interpretation, limitations, and common misconceptions

A common misconception is to read “general” as “arbitrary.” The surveyed literature does not support that reading. In the continuum-model formalism, admissible GBC are precisely the current-conserving, self-adjoint boundary relations parameterized by a unitary matrix, and current-nullifying relations with unbalanced positive and negative current sectors are ruled out as non-self-adjoint (Kharitonov, 2022). In obstacle scattering, the title’s phrasing is explicitly narrower than arbitrary local elliptic BCs and is concretely the class of Neumann and smooth Robin conditions (Jin, 2014). In qBounce, the rigorous statement is narrower still: self-adjointness on the half-line forces a one-parameter Robin family, not an unrestricted boundary operator algebra (Sung et al., 17 Oct 2025).

A second misconception is that the qBounce parameter $x=0$7 is already a microscopic material constant. The paper does not make that claim. What is exact is the self-adjoint-extension result and the boundary condition

$x=0$8

What remains phenomenological is the physical interpretation of $x=0$9, which may encode short-distance mirror physics, imperfect reflection, finite penetration, surface interaction details, mirror composition or coating dependence, and omitted systematic effects in the effective Hamiltonian. A plausible implication is that future microscopic neutron–mirror scattering models would be needed to decide how much of an experimentally fitted $H(p)=\sum_{n=0}^N h_n p^n,\qquad p=-i\partial_x,$00 is genuine boundary physics and how much is effective bookkeeping (Sung et al., 17 Oct 2025).

The broader significance of GBC is therefore methodological as much as physical. Enforcing self-adjointness or its field-specific analogue changes not only boundary values but spectra, matrix elements, effective Hamiltonians, and experimental inference. In the qBounce paper this point is made explicitly for gravitational quantum states, and the same work states that the formalism is broadly applicable to confined particles, Stark-like linear potentials, bouncing Bose condensates, quarkonium models, and more generally to quantum problems in which domain issues determine observables (Sung et al., 17 Oct 2025).