First-Order Green Operator

• A first-order Green operator is defined as the inverse or generalized inverse of a differential operator, ensuring unique solutions that honor boundary data and causal support.
• It is applied in perturbation theory, boundary-value problems, and hyperbolic systems, with roles in both ODE/PDE settings and electromagnetic formulations.
• Its implementations range from solving classical differential equations on bounded domains to enabling causal propagation on globally hyperbolic manifolds and quantization in curved spacetimes.

Searching arXiv for the cited papers to ground the article in current arXiv records. A first-order Green operator is an operator-valued solution construct associated with a first-order differential equation, but the literature uses the term in several distinct senses. For first-order linear PDEs and ODEs it can denote the inverse or generalized inverse that maps a source to the unique solution satisfying boundary data; on globally hyperbolic Lorentzian manifolds it denotes the advanced and retarded Green operators of prenormally hyperbolic or Green-hyperbolic first-order operators; in perturbation theory it denotes the linear term $G_1$ in the expansion of a resolvent; and in boundary $\psi$do theory it can denote a first-order boundary operator entering Green’s formula. Across these settings, the construction organizes inversion, causal propagation, or boundary correction at first order (Mkrtchian et al., 2020, Muehlhoff, 2010, Martin, 2011, Korporal et al., 2013, Grubb, 2016, Baer, 2013, Agarwal et al., 29 Mar 2026).

1. Terminological scope and core operator identities

In the Lorentzian setting, a first-order operator on a vector bundle $\mathcal E\to M$ is written locally as

$P \;=\; A^\mu(x)\,\nabla_\mu \;+\; B(x),$

or equivalently $P=\nabla_X+\Phi$. When $P$ is prenormally hyperbolic on a globally hyperbolic spacetime, its advanced and retarded Green operators are the unique linear maps

$G_\pm:\Gamma_0(\mathcal E)\to\Gamma(\mathcal E)$

satisfying

$P\,G_\pm=\Id,\qquad G_\pm\,P=\Id,$

on compact support together with

$\supp(G_\pm f)\subseteq J^\pm(\supp f).$

These identities characterize causal inversion for first-order hyperbolic dynamics (Muehlhoff, 2010).

In bounded-domain boundary-value theory, the same expression refers instead to the integral operator

$u(x)=\int_\Omega G(x,y)\,f(y)\,dy,$

where $\psi$0 is obtained from a fundamental solution $\psi$1 by a boundary correction that enforces homogeneous linear, possibly nonlocal, boundary conditions. In ordinary differential boundary problems, the Green operator is a generalized inverse $\psi$2 satisfying abstract projector identities of the form $\psi$3 and $\psi$4 (Mkrtchian et al., 2020, Korporal et al., 2013).

A different usage appears in perturbation theory. If $\psi$5 and $\psi$6, then the “first-order Green operator” is

$\psi$7

the linear sensitivity of the resolvent to the perturbation. In fractional-order boundary theory, a first-order Green operator is instead a boundary $\psi$8do $\psi$9 of order $\mathcal E\to M$0 appearing in the boundary pairing of Green’s formula (Martin, 2011, Grubb, 2016).

2. Boundary-integral construction for first-order linear PDEs

For a bounded domain $\mathcal E\to M$1 with smooth boundary, consider the general first-order linear differential operator

$\mathcal E\to M$2

together with homogeneous linear boundary conditions $\mathcal E\to M$3. Mkrtchian and Henkel construct the Green function from a fundamental solution $\mathcal E\to M$4 on all of $\mathcal E\to M$5, defined distributionally by

$\mathcal E\to M$6

Under fairly general conditions, such an $\mathcal E\to M$7 exists and is unique up to addition of a homogeneous solution, and one chooses the standard “causal” or “principal-part” fundamental solution (Mkrtchian et al., 2020).

The boundary operator is allowed to be nonlocal. With $\mathcal E\to M$8 boundary conditions,

$\mathcal E\to M$9

the adjoint boundary operators $P \;=\; A^\mu(x)\,\nabla_\mu \;+\; B(x),$0 have kernels $P \;=\; A^\mu(x)\,\nabla_\mu \;+\; B(x),$1. The key auxiliary object is the $P \;=\; A^\mu(x)\,\nabla_\mu \;+\; B(x),$2 boundary-to-boundary operator $P \;=\; A^\mu(x)\,\nabla_\mu \;+\; B(x),$3, whose kernel is

$P \;=\; A^\mu(x)\,\nabla_\mu \;+\; B(x),$4

The construction assumes that $P \;=\; A^\mu(x)\,\nabla_\mu \;+\; B(x),$5 is invertible as an integral operator on $P \;=\; A^\mu(x)\,\nabla_\mu \;+\; B(x),$6 (Mkrtchian et al., 2020).

The boundary density $P \;=\; A^\mu(x)\,\nabla_\mu \;+\; B(x),$7 is then defined by

$P \;=\; A^\mu(x)\,\nabla_\mu \;+\; B(x),$8

and the Green function satisfying

$P \;=\; A^\mu(x)\,\nabla_\mu \;+\; B(x),$9

is

$P=\nabla_X+\Phi$0

One checks directly that $P=\nabla_X+\Phi$1, that the boundary correction is $P=\nabla_X+\Phi$2-homogeneous for $P=\nabla_X+\Phi$3, and that $P=\nabla_X+\Phi$4 (Mkrtchian et al., 2020).

This construction isolates boundary nonlocality in the kernels $P=\nabla_X+\Phi$5 and their adjoints. The operator $P=\nabla_X+\Phi$6 collects multiple boundary reflections, its inverse sums them to enforce $P=\nabla_X+\Phi$7, and the density $P=\nabla_X+\Phi$8 implements the correction of the fundamental solution required by the boundary conditions. The resulting Green operator is

$P=\nabla_X+\Phi$9

3. Ordinary first-order boundary problems and generalized inverses

For the first-order ordinary differential operator

$P$0

on an interval $P$1, with boundary condition $P$2, the Green operator appears in the generalized boundary-problem framework of Korporal and Regensburger. With $P$3, $P$4, $P$5, boundary functional $P$6, and exceptional space $P$7, the triple $P$8 is a regular boundary problem admitting a unique Green operator

$P$9

that maps each $G_\pm:\Gamma_0(\mathcal E)\to\Gamma(\mathcal E)$0 to the unique $G_\pm:\Gamma_0(\mathcal E)\to\Gamma(\mathcal E)$1 satisfying $G_\pm:\Gamma_0(\mathcal E)\to\Gamma(\mathcal E)$2 and $G_\pm:\Gamma_0(\mathcal E)\to\Gamma(\mathcal E)$3 (Korporal et al., 2013).

In the notation of the general theory,

$G_\pm:\Gamma_0(\mathcal E)\to\Gamma(\mathcal E)$4

Here $G_\pm:\Gamma_0(\mathcal E)\to\Gamma(\mathcal E)$5, because $G_\pm:\Gamma_0(\mathcal E)\to\Gamma(\mathcal E)$6 and $G_\pm:\Gamma_0(\mathcal E)\to\Gamma(\mathcal E)$7, while $G_\pm:\Gamma_0(\mathcal E)\to\Gamma(\mathcal E)$8 is the projector onto $G_\pm:\Gamma_0(\mathcal E)\to\Gamma(\mathcal E)$9 along $P\,G_\pm=\Id,\qquad G_\pm\,P=\Id,$0. Since

$P\,G_\pm=\Id,\qquad G_\pm\,P=\Id,$1

the projector is

$P\,G_\pm=\Id,\qquad G_\pm\,P=\Id,$2

and therefore

$P\,G_\pm=\Id,\qquad G_\pm\,P=\Id,$3

No compatibility condition on $P\,G_\pm=\Id,\qquad G_\pm\,P=\Id,$4 is needed; every $P\,G_\pm=\Id,\qquad G_\pm\,P=\Id,$5 yields a unique solution (Korporal et al., 2013).

Variation of constants gives the integral representation

$P\,G_\pm=\Id,\qquad G_\pm\,P=\Id,$6

so the Green kernel is

$P\,G_\pm=\Id,\qquad G_\pm\,P=\Id,$7

This is the elementary first-order prototype of a causal kernel on an interval (Korporal et al., 2013).

Within the same framework, the first-order Green operator is also an atomic factor in higher-order factorizations. If a higher-order operator factors into first-order pieces, the overall Green operator satisfies a reverse-order law, and in the second-order case $P\,G_\pm=\Id,\qquad G_\pm\,P=\Id,$8 the composite regular problem has Green operator $P\,G_\pm=\Id,\qquad G_\pm\,P=\Id,$9. This places the single-step first-order Green operator inside the multiplicative structure of generalized inverses and boundary problems (Korporal et al., 2013).

4. Prenormally hyperbolic first-order operators on globally hyperbolic spacetimes

Mühlhoff studies first-order linear differential operators on vector bundles over globally hyperbolic Lorentzian manifolds. A first-order operator $\supp(G_\pm f)\subseteq J^\pm(\supp f).$0 is called prenormally hyperbolic if there exists another first-order operator $\supp(G_\pm f)\subseteq J^\pm(\supp f).$1 such that the composition $\supp(G_\pm f)\subseteq J^\pm(\supp f).$2 is normally hyperbolic of second order, with principal symbol

$\supp(G_\pm f)\subseteq J^\pm(\supp f).$3

This condition is the mechanism by which a first-order system inherits well-posed hyperbolic propagation from a second-order normally hyperbolic operator (Muehlhoff, 2010).

The main theorem states that if $\supp(G_\pm f)\subseteq J^\pm(\supp f).$4 is prenormally hyperbolic on a globally hyperbolic $\supp(G_\pm f)\subseteq J^\pm(\supp f).$5, then there exist unique advanced and retarded Green operators

$\supp(G_\pm f)\subseteq J^\pm(\supp f).$6

satisfying

$\supp(G_\pm f)\subseteq J^\pm(\supp f).$7

for compactly supported data, together with

$\supp(G_\pm f)\subseteq J^\pm(\supp f).$8

These are the defining identities of first-order Green operators in the causal sense (Muehlhoff, 2010).

The proof reduces the problem to second order. One chooses $\supp(G_\pm f)\subseteq J^\pm(\supp f).$9 so that $u(x)=\int_\Omega G(x,y)\,f(y)\,dy,$0 is normally hyperbolic, uses the known Green operators $u(x)=\int_\Omega G(x,y)\,f(y)\,dy,$1 of $u(x)=\int_\Omega G(x,y)\,f(y)\,dy,$2, and sets

$u(x)=\int_\Omega G(x,y)\,f(y)\,dy,$3

The identities $u(x)=\int_\Omega G(x,y)\,f(y)\,dy,$4 and $u(x)=\int_\Omega G(x,y)\,f(y)\,dy,$5 are verified using the adjoint problem and Lorentzian integration by parts, while the support property follows from $u(x)=\int_\Omega G(x,y)\,f(y)\,dy,$6 and the causal support of $u(x)=\int_\Omega G(x,y)\,f(y)\,dy,$7. Uniqueness follows because any difference of two candidates would be supported in a compact set both in the future and the past of itself, forcing it to vanish (Muehlhoff, 2010).

This causal formulation is also the entry point to algebraic quantization. The advanced and retarded Green operators furnish the basic building blocks for constructing global solutions and for canonical CAR/CCR constructions on curved spacetimes (Muehlhoff, 2010).

5. Green-hyperbolic systems, Dirac-type operators, and structural closure

A more systematic notion is that of a Green-hyperbolic operator. For a globally hyperbolic Lorentzian manifold $u(x)=\int_\Omega G(x,y)\,f(y)\,dy,$8 and a finite-dimensional vector bundle $u(x)=\int_\Omega G(x,y)\,f(y)\,dy,$9, a first-order operator

$\psi$00

is called Green-hyperbolic if both $\psi$01 and its formal dual $\psi$02 admit advanced and retarded Green’s operators. On compactly supported smooth sections,

$\psi$03

The cited systematic study proves that such $\psi$04 are unique (Baer, 2013).

A major class of first-order examples is given by symmetric hyperbolic systems. Their construction proceeds through a Cauchy temporal function $\psi$05, spacelike slices $\psi$06, an energy estimate, and the resulting uniqueness, finite speed of propagation, local existence, global existence by patching local solutions, and stability. The advanced Green operator is then defined by solving

$\psi$07

for $\psi$08 chosen so that $\psi$09. By finite speed and uniqueness this yields a well-defined $\psi$10; the retarded operator $\psi$11 is analogous. No closed-form kernel exists in general on curved backgrounds, but in special cases such as Minkowski space one recovers integral-kernel representations in terms of light-cone integrals (Baer, 2013).

The Dirac operator is a canonical first-order example. On a spin spacetime, with spinor bundle $\psi$12,

$\psi$13

satisfies

$\psi$14

so by the square-root property it is Green-hyperbolic. Concretely,

$\psi$15

where $\psi$16 are the Klein–Gordon Green’s operators (Baer, 2013).

The class is structurally stable. It is closed under composition, square roots, direct sums, and adjoint-factorizations. Thus, once one has seed examples such as wave operators and Dirac-type operators, one can build a large class of first-order Green-hyperbolic operators with unique advanced and retarded Green operators and controlled causal support (Baer, 2013).

6. Perturbative and boundary $\psi$17do meanings of “first-order Green operator”

In perturbation theory for elliptic boundary-value problems, the unperturbed Green operator $\psi$18 is the inverse of the Laplacian with Dirichlet boundary condition on a bounded $\psi$19 or real-analytic domain,

$\psi$20

As an operator on $\psi$21,

$\psi$22

is bounded $\psi$23, and $\psi$24 in $\psi$25 with zero boundary values. For perturbations

$\psi$26

the first-order Green operator is

$\psi$27

For the Helmholtz, Schrödinger, and Laplace–Beltrami perturbations, this gives the first-order correction kernels explicitly, and $\psi$28 encodes the linear sensitivity of the resolvent to the perturbation. The same paper also reviews Hadamard’s first-variation formula for domain deformations (Martin, 2011).

A distinct boundary-theoretic usage appears in Grubb’s analysis of fractional-order classical $\psi$29do’s $\psi$30 of order $\psi$31 with even symbol. The Green formula is

$\psi$32

where $\psi$33 is a classical $\psi$34do on $\psi$35 of order $\psi$36. In this setting, $\psi$37 is the first-order Green operator entering the boundary pairing (Grubb, 2016).

The operator $\psi$38 emerges after introducing the order-reducing operators $\psi$39, forming the conjugated order-zero operator $\psi$40, factoring $\psi$41 in Boutet de Monvel’s calculus, and extracting the singularity of the Poisson-trace composition $\psi$42 at the boundary. Its symbol is obtained from the jump of a bounded part of the symbol across $\psi$43, and as a first-order boundary $\psi$44do it acts continuously

$\psi$45

In the same framework, the Dirichlet-to-Neumann operator

$\psi$46

is a classical $\psi$47do of order $\psi$48, and the Neumann problem is Fredholm solvable exactly when $\psi$49 is elliptic (Grubb, 2016).

These two usages are mathematically adjacent but conceptually different. In one case, “first-order Green operator” denotes the first derivative of a resolvent with respect to an operator perturbation; in the other, it denotes an order-one boundary operator in a Green identity.

7. First-order Maxwell Green operators in macroscopic quantum electrodynamics

A recent electromagnetic usage rewrites Maxwell’s equations as a first-order operator equation for the dual field

$\psi$50

with

$\psi$51

$\psi$52

Maxwell’s equations become

$\psi$53

The retarded first-order Green operator is

$\psi$54

so that $\psi$55 (Agarwal et al., 29 Mar 2026).

Its position-space kernel $\psi$56 satisfies

$\psi$57

with the Sommerfeld radiation condition at infinity. Love’s equivalence yields the Stratton–Chu representation

$\psi$58

where

$\psi$59

This makes the first-order Green operator a propagator of the electromagnetic state between surfaces (Agarwal et al., 29 Mar 2026).

Two adjoint symmetries organize the theory. Under the energy inner product one obtains the generalized optical theorem

$\psi$60

and therefore

$\psi$61

which decomposes $\psi$62 into bulk-absorption and boundary-radiation channels. Under the reciprocal inner product for reciprocal media, $\psi$63, yielding Lorentz reciprocity at kernel level,

$\psi$64

These identities are derived with minimal vector calculus from operator symmetries (Agarwal et al., 29 Mar 2026).

The same formalism supports Heisenberg–Langevin quantization. After introducing dual polarization and magnetization operators coupled to local bosonic baths, one obtains

$\psi$65

with a bulk Langevin commutator fixed by the fluctuation-dissipation theorem. Inverting Maxwell’s equation gives

$\psi$66

Adding bulk and boundary contributions yields the exact commutator

$\psi$67

which is consistent with the fluctuation-dissipation theorem even when dielectrics extend to the boundary. The associated exact transfer law between two nonintersecting surfaces uses

$\psi$68

and yields pseudo-unitarity in the lossless limit (Agarwal et al., 29 Mar 2026).