Dirac Pulse Boundary Condition

• Dirac Pulse Boundary Condition is a term used to describe localized boundary data in PDEs and Dirac systems, encompassing Dirac delta pulses, self-adjoint extension conditions, and singular reductions.
• It covers various models such as the Gurtin–Pipkin equation, 1D Dirac operators, and radial Schrödinger reductions, each demonstrating distinct mathematical treatments and physical interpretations.
• The concept differentiates between literal impulsive inputs and abstract self-adjoint junction conditions, thereby clarifying common misconceptions in boundary treatments.

Searching arXiv for papers on "Dirac pulse boundary condition" and closely related Dirac boundary-condition literature. arXiv search query: "Dirac pulse boundary condition" In the literature surveyed here, the expression “Dirac pulse boundary condition” is best understood as a non-standard umbrella label rather than a single canonical boundary prescription. It can denote, depending on context, a literal Dirac delta boundary input such as $u(t)=\delta(t)$, a point-supported or zero-range coupling encoded by self-adjoint extension data for a Dirac operator, or a distributional $\delta$-term generated at a singular point unless an auxiliary boundary condition is imposed. The common structure is localization at a boundary, interface, or distinguished point, but the mathematical realizations differ substantially across PDEs, Dirac operators, and radial reductions (Ivanov, 2013, Angelone, 2023, Etxebarria, 2013).

1. Terminological scope and principal meanings

The sources considered here suggest that the phrase does not identify a unique standard object. In one line of work, the boundary datum is literally a Dirac delta pulse in time; in another, the relevant object is a self-adjoint matching condition at a point defect; in a third, the issue is a distributional $\delta$-source that appears when a singular coordinate reduction is treated incorrectly. This suggests that the phrase is best used only with an explicit specification of the operator, domain, and type of boundary localization (Ivanov, 2013, Angelone, 2023, Etxebarria, 2013).

Setting	Boundary object	Representative formulation
Gurtin–Pipkin equation	Dirac delta boundary input	$\theta(0,t)=\delta(t)$
1D Dirac operator with junction	Point-junction self-adjoint extension	$\Psi_{\mathrm D,-}=U\Psi_{\mathrm D,+}$, $U\in U(2)$
Radial Schrödinger reduction	Distributional point term at origin	$\Delta(1/r)=-4\pi\delta(\mathbf r)$
Half-line Dirac evolution	Non-homogeneous boundary trace	$w_1(0,t)=h_1(t)$
Dirac equation on AdS	Boundary data at conformal infinity	$U(N/2)$-family, MIT-bag as a special case

A central distinction runs through the literature. A literal pulse potential is a singular coefficient in the differential equation or a prescribed distributional trace, whereas an abstract self-adjoint point interaction is a domain condition for the operator. The one-dimensional ring analysis makes this distinction explicit: the $U(2)$ family gives the complete set of self-adjoint point/junction conditions for the Dirac operator, but the paper does not identify every such condition with a literal $\delta$0-potential or short pulse (Angelone, 2023).

2. Dirac delta pulses as boundary data in evolutionary PDEs

The clearest literal boundary pulse in the sources appears in the one-dimensional Gurtin–Pipkin equation on the half-line,

$\delta$1

with

$\delta$2

Here the boundary value itself is a Dirac delta distribution concentrated at $\delta$3 (Ivanov, 2013).

Under the high-frequency assumption

$\delta$4

the solution has the singular decomposition

$\delta$5

with $\delta$6 continuous near the front $\delta$7. In the smooth-kernel formulation summarized in the same paper, the leading singularity is

$\delta$8

The result is therefore not instantaneous diffusive spreading but a moving Dirac mass concentrated on the characteristic line $\delta$9, with exponentially decreasing amplitude (Ivanov, 2013).

This boundary-pulse mechanism is sharply different from a Dirac-operator boundary condition. The boundary datum is inhomogeneous and distribution-valued, and the analysis is carried out by Laplace transform rather than by self-adjoint extension theory. Its conceptual importance lies in showing that a boundary $\delta$0-pulse can propagate as an interior singular front, rather than being merely a formal impulsive input.

3. Non-homogeneous boundary forcing for Dirac systems

For genuine Dirac evolution on the half-line, the quarter-plane analysis studies

$\delta$1

on

$\delta$2

with initial data

$\delta$3

and boundary condition

$\delta$4

Only one boundary trace is prescribed, which reflects the characteristic structure of the first-order system in this representation (Babaei et al., 17 Jun 2026).

The theory is rigorous for smooth data,

$\delta$5

and it makes boundary regularity depend on explicit compatibility conditions. Continuity up to the corner requires

$\delta$6

while $\delta$7-regularity for the homogeneous problem requires

$\delta$8

For the forced system

$\delta$9

the corrected second condition is

$\theta(0,t)=\delta(t)$0

The paper is explicit that a literal boundary condition such as $\theta(0,t)=\delta(t)$1 is not covered by this framework, but smooth pulse approximations are naturally accommodated through the exact integral formulas obtained by the Fokas unified transform (Babaei et al., 17 Jun 2026).

A complementary numerical perspective is provided by the staggered-grid leap-frog scheme with discrete transparent boundary conditions for the $\theta(0,t)=\delta(t)$2D Dirac equation. In the exterior semi-infinite leads, the discrete boundary laws are

$\theta(0,t)=\delta(t)$3

These DTBCs are derived from the fully discrete exterior problem by Z-transform and enforce outgoing behavior without spurious reflection. In the massless zero-potential case with $\theta(0,t)=\delta(t)$4, the kernel simplifies to $\theta(0,t)=\delta(t)$5, so the transparent boundary condition becomes local: $\theta(0,t)=\delta(t)$6 Within the language of pulse propagation, these are exact outflow conditions for localized Dirac wave packets rather than point-supported pulse sources (Hammer et al., 2013).

4. Point-supported Dirac boundary conditions and zero-range couplings

For the one-dimensional free Dirac operator on a ring cut at one point, the admissible boundary data are classified by self-adjoint extension theory. The bulk operator is

$\theta(0,t)=\delta(t)$7

and the deficiency indices are

$\theta(0,t)=\delta(t)$8

Hence the self-adjoint extensions are parameterized by $\theta(0,t)=\delta(t)$9, with domain

$\Psi_{\mathrm D,-}=U\Psi_{\mathrm D,+}$0

Equivalently, in components,

$\Psi_{\mathrm D,-}=U\Psi_{\mathrm D,+}$1

This is the complete four-parameter family of self-adjoint point/junction conditions for the model (Angelone, 2023).

The same paper emphasizes a distinction that is crucial for the present topic. These $\Psi_{\mathrm D,-}=U\Psi_{\mathrm D,+}$2 relations are the correct abstract framework for a sharply localized defect, but they are more general than any one specific $\Psi_{\mathrm D,-}=U\Psi_{\mathrm D,+}$3 potential or pulse model. The anti-diagonal family

$\Psi_{\mathrm D,-}=U\Psi_{\mathrm D,+}$4

yields the pure phase-jump gluing

$\Psi_{\mathrm D,-}=U\Psi_{\mathrm D,+}$5

which is point-like and transmitting, whereas the diagonal family $\Psi_{\mathrm D,-}=U\Psi_{\mathrm D,+}$6 gives local chiral endpoint conditions and is closer to bag-wall confinement (Angelone, 2023).

The general local theory of Dirac-type operators sharpens this point. A local smooth boundary condition is a smooth subbundle

$\Psi_{\mathrm D,-}=U\Psi_{\mathrm D,+}$7

and for $\Psi_{\mathrm D,-}=U\Psi_{\mathrm D,+}$8 the operator is self-adjoint iff $\Psi_{\mathrm D,-}=U\Psi_{\mathrm D,+}$9 is the graph of a unitary map

$U\in U(2)$0

Regularity is governed by the Shapiro–Lopatinski condition

$U\in U(2)$1

For transmission conditions on an interface $U\in U(2)$2, defined by endomorphisms $U\in U(2)$3, self-adjointness and regularity are likewise reduced to explicit fiberwise algebraic conditions (Große et al., 2024).

The codimension of the singular set can obstruct such constructions. For the free $U\in U(2)$4D Dirac operator on truncated Fock space, there is no self-adjoint interior-boundary-condition Hamiltonian that couples a point source at the origin to lower sectors; every self-adjoint extension is block diagonal. By contrast, after adding a sufficiently strong Coulomb singularity,

$U\in U(2)$5

the paper proves the existence of self-adjoint IBC Hamiltonians with nontrivial particle creation, using the asymptotic coefficients $U\in U(2)$6 and $U\in U(2)$7 in the singular expansion near the origin (Henheik et al., 2020).

5. Canonical Dirac boundary-condition families in geometry and materials

On globally hyperbolic Lorentzian manifolds with timelike boundary, a rigorous class of boundary conditions is given by families $U\in U(2)$8 that are local in time and non-local in the spatial directions. The prototype is the slicewise APS condition

$U\in U(2)$9

where $\Delta(1/r)=-4\pi\delta(\mathbf r)$0 is the induced boundary Dirac operator on $\Delta(1/r)=-4\pi\delta(\mathbf r)$1. The resulting Cauchy problem is well posed for admissible families, but the paper is explicit that it does not treat boundary data of the form $\Delta(1/r)=-4\pi\delta(\mathbf r)$2 or inhomogeneous distributional traces (Baer et al., 2022).

On $\Delta(1/r)=-4\pi\delta(\mathbf r)$3, the conformal boundary is timelike and the mass window determines whether boundary data are needed. Writing

$\Delta(1/r)=-4\pi\delta(\mathbf r)$4

the classification is: $\Delta(1/r)=-4\pi\delta(\mathbf r)$5 In the nontrivial window $\Delta(1/r)=-4\pi\delta(\mathbf r)$6, the generalized boundary conditions in the Poincaré patch take the form

$\Delta(1/r)=-4\pi\delta(\mathbf r)$7

with MIT bag as a distinguished special case. The paper also shows that generalized choices can support normalizable bound states, while the proper MIT-bag example does not exhibit that feature in the worked $\Delta(1/r)=-4\pi\delta(\mathbf r)$8 construction (Dappiaggi et al., 13 Nov 2025).

In graphene and other Dirac-like condensed-matter models, local boundary conditions are likewise organized by algebraic projector data. For graphene nanoribbons and nanodots, one admissible family is

$\Delta(1/r)=-4\pi\delta(\mathbf r)$9

with MIT bag corresponding to

$w_1(0,t)=h_1(t)$0

equivalently

$w_1(0,t)=h_1(t)$1

A different but related line of work shows that a symmetry-allowed Wilson mass term is equivalent to Berry–Mondragon discontinuous boundary conditions, so that the regularized $w_1(0,t)=h_1(t)$2 model can be solved with the simple hard-wall condition $w_1(0,t)=h_1(t)$3 while reproducing the same low-energy boundary physics (Beneventano et al., 2010, Araújo et al., 2019).

6. Distributional point terms, origin conditions, and persistent misconceptions

A particularly instructive analogue to pulse-like boundary localization occurs in the radial Schrödinger equation. For a central potential one writes

$w_1(0,t)=h_1(t)$4

and for $w_1(0,t)=h_1(t)$5 obtains the reduced radial equation

$w_1(0,t)=h_1(t)$6

The subtlety is that the substitution $w_1(0,t)=h_1(t)$7 is singular at the origin, and the distributional identity

$w_1(0,t)=h_1(t)$8

produces an extra point-supported term unless

$w_1(0,t)=h_1(t)$9

For $U(N/2)$0, the exact equation contains

$U(N/2)$1

with the understood meaning that the distribution acts through $U(N/2)$2. The conclusion is that the usual reduced radial equation is exactly equivalent to the three-dimensional Schrödinger equation only if

$U(N/2)$3

Otherwise one has implicitly introduced an additional contact interaction at the origin (Etxebarria, 2013).

This example corrects a recurring misconception: point-supported singular terms are not always optional regularity artifacts, and normalizability alone does not settle the issue. The same caution applies in Dirac problems. A point-supported interaction may be a literal boundary pulse, an interface jump, or an abstract self-adjoint extension; these possibilities are mathematically distinct, and the literature does not identify them automatically with one another. The safest usage therefore distinguishes among three cases: prescribed inhomogeneous boundary data, self-adjoint boundary or junction conditions, and spurious or induced $U(N/2)$4-terms generated by singular reductions (Etxebarria, 2013, Angelone, 2023).

In that sense, the phrase “Dirac pulse boundary condition” is most precise when it is anchored to a concrete model. In boundary-forced evolutionary PDEs it means a Dirac delta input such as $U(N/2)$5; in one-dimensional Dirac theory it usually refers more safely to zero-range point/junction conditions parameterized by $U(N/2)$6; in geometric Dirac problems it refers to a self-adjoint extension at timelike or spatial boundary; and in singular radial reductions it names the appearance of a point-supported $U(N/2)$7-source that forces an auxiliary boundary condition. The unifying theme is localization, but the governing notions are different: distributional forcing, current-conserving self-adjointness, and equivalence of singular reductions.