Preprint
Article

This version is not peer-reviewed.

Real Differentiation: A Fiberwise, Stack-Enriched Calculus of Change

Submitted:

24 December 2025

Posted:

25 December 2025

You are already at the latest version

Abstract
We develop Real Differentiation: a fibre- and stack-sensitive language for infinitesimal change whose basic objects are labelled (“energy-tagged”) fields and operators. The formalism is designed to enforce a no branch mixing principle: differentiation, integration, and PDE operators act fibrewise on each labelled branch, while still permitting monodromy-aware transport when labels are organised as a covering space or stack of analytic continuations. The core construction is a jet-theoretic endofunctor D, defined via the first infinitesimal neighbourhood of the diagonal. In the constant-label case \( \mathbb{E}\mathcal{x} \) , we show that the labelled topos (Sh(\( \mathcal{B} \))/\( \mathcal{X} \) )/\( \mathbb{E}\mathcal{x} \) decomposes canonically as a product of copies of Sh(\( \mathcal{B} \))/\( \mathcal{X} \) , so that labelled objects are literally families indexed by energy values. This makes “no branch mixing” a theorem: all constructions in the labelled topos act componentwise. In the stacky/local-system case, labels are organised by an étale map (or stack) \( \mathcal{E} \)\( \mathcal{X} \) ; labelled objects are sheaves on the total space and carry monodromy by descent. We further define directional derivatives by contraction of the universal derivation with vector fields; we prove Leibniz and chain rules, including a smooth/analytic chain rule under explicit functional-calculus hypotheses. Higher jets are related to differential operators, and an infinite-jet comonad governs the differential-operator calculus. Finally, we introduce controlled energy mixing via correspondences of label objects, providing a principled way to model coupled branch dynamics.
Keywords: 
;  ;  ;  ;  ;  ;  ;  

1. Motivation and Outline

Conventional calculus differentiates numbers. In spectral, multivalued, or analytic-continuation settings—where “a function” may in fact live on multiple branches and is only single-valued after choosing a sheet—naïve differentiation on a chosen principal branch is not stable under monodromy, and can break fundamental identities (e.g. pathwise versions of the fundamental theorem).
Real Differentiation replaces unlabelled scalars by labelled fields: objects living in a slice topos over a label object. In the simplest case, labels form a constant discrete object E ̲ X ; then labelled objects are equivalent to families indexed by energy values, and “no branch mixing” is enforced categorically. In the monodromy-aware case, labels form a covering space or stack E X ; then labelled objects are sheaves on the total space, and monodromy is encoded by descent along E X rather than being ignored.
After setting up the labelled foundations (Section 2), we define the Real Differential functor D via first jets (Section 3), prove fibrewise base-change properties, and define directional derivatives by universal differentials (Section 4). We then extend the chain rule beyond polynomials under explicit hypotheses, develop higher jets and a jet comonad governing differential operators (Sections 5 and 6), construct the Real de Rham complex and fibrewise integration (Section 7), and introduce controlled energy mixing via label correspondences (Section 9). Computational demonstrations are provided as in-document listings (Appendix A.2 and Appendix A.3).

Positioning and Scope

Jets, Kähler differentials, and the jet/differential-operator adjunction are classical (EGA, Hartshorne, Illusie). The contribution here is the systematic labelled packaging: (i) a decomposition theorem showing constant-labelled objects are families; (ii) a cover/stack label model capturing monodromy; (iii) a jet-based differentiation functor and identities that commute with passage to label fibres; and (iv) a controlled mixing extension via correspondences.

2. Energy–Fibred Foundations

2.1. Energy FIELD, Base Topos, Ringed Structure

Fix a commutative real-closed energy field  ( E , + , · ) . Let B be a Grothendieck site, and work in the ringed topos ( Sh ( B ) , O B ) . Unless stated otherwise, all internal algebra is O B -linear.
Remark 1
(Size convention). To form coproducts indexed by underlying sets of energy values, we assume the underlying set of E is small in a fixed Grothendieck universe, or we work universe-relatively. This is standard in topos-theoretic constructions and does not affect the geometric content.
Definition 1
(Energy-affine space). For n 0 set
A E n : = Spec E [ t 1 , , t n ] ,
viewed as an object of B .

2.2. Geometry over X and Relative Differentials

Fix a morphism π : X A E n in B , assumed energy-smooth when needed (Definition 9). We write ( Sh ( B ) / X , O X ) for the induced ringed topos over X . Relative Kähler differentials Ω X / E 1 are an O X -module, and we set
T X / E : = Hom ̲ O X ( Ω X / E 1 , O X ) .
When discussing vector fields, derivations, jets, and de Rham theory, we work in ( Sh ( B ) / X , O X ) .

2.3. Label Objects and Labelled Module Objects

A central organising idea is to separate two layers: (i) a geometric base X , and (ii) a label object  E over X whose fibres represent branch/energy/sheet data.
Definition 2
(Label object and labelled topos). A label object over X   is an object E of Sh ( B ) / X . The associated labelled topos is the slice
( Sh ( B ) / X ) / E .
An  E -labelled object is a morphism F E in Sh ( B ) / X .
Definition 3
(Labelled structure sheaf and labelled modules). Let p : ( Sh ( B ) / X ) / E Sh ( B ) / X be the projection geometric morphism. Define the pulled-back structure sheaf
O X E : = p * O X ,
a ring object in the labelled topos. An  O X E -modulemeans a module object in ( Sh ( B ) / X ) / E over the ring O X E .
Remark 2
(Why modules must be internal). This formalism avoids a subtle but important pitfall: a label map is not naturally an O X -linear morphism to a ring of scalars. Instead, labels are tracked by working in the slice topos, where all algebra (including module structure) is internal and automatically respects the labelling.

2.4. Constant Energy Labels and the Slice-as-Family Decomposition

We now define the constant energy label object and prove that it yields a canonical product decomposition.
Definition 4
(Constant energy label object E ̲ X ). Let 1 X denote the terminal object of Sh ( B ) / X . The constant energy label object on X  is the constant object with fibre E ,
E ̲ X : = Δ X ( E ) E E 1 X in Sh ( B ) / X .
For each E E , write ι E : 1 X E ̲ X for the coproduct inclusion.
Definition 5
(Energy fibre (constant-label case)). Let F E ̲ X be an E ̲ X -labelled object in Sh ( B ) / X . For E E , its energy fibre is the pullback
F E : = F × E ̲ X 1 X along ι E : 1 X E ̲ X .
Theorem 1
(Slice-as-family decomposition for constant labels). There is an equivalence of toposes
( Sh ( B ) / X ) / E ̲ X E E ( Sh ( B ) / X ) ,
under which a labelled object F E ̲ X corresponds to its family of fibres ( F E ) E E , and morphisms correspond to families of morphisms F E G E .
Proof. 
In any Grothendieck topos, coproducts are disjoint and stable under pullback (extensivity). Since E ̲ X E E 1 X , slicing over E ̲ X decomposes as a product of slices over each summand:
( Sh ( B ) / X ) / E E 1 X E E ( Sh ( B ) / X ) / 1 X .
But ( Sh ( B ) / X ) / 1 X Sh ( B ) / X because 1 X is terminal. The equivalence is realised by pulling back along the inclusions ι E , yielding the family of fibres F E , and conversely by taking the coproduct E F E E 1 X = E ̲ X .    □
Corollary 1
(Canonical coproduct decomposition). For any E ̲ X -labelled object F E ̲ X there is a canonical isomorphism in Sh ( B ) / X :
F E E F E .
Proposition 1
(Modules in the constant-labelled topos are families of modules). The category of O X E ̲ X -modules in the labelled topos ( Sh ( B ) / X ) / E ̲ X is canonically equivalent to the product category
O X E ̲ X - Mod ( Sh ( B ) / X ) / E ̲ X E E O X - Mod Sh ( B ) / X .
In particular, any construction in the labelled topos defined by topos limits, colimits, and O X -module operations acts componentwise on fibres.
Proof. 
Apply Theorem 1 and observe that pulling back O X E ̲ X = p * O X to the Eth factor recovers O X . Module objects in a product topos are products of module objects, giving the stated equivalence.    □

2.5. Reformulated Tag Models: Model A and Model B

We now restate your two tag semantics in the rigorous language of Theorem 1.
Definition 6
(Model A: single-energy (concentrated) objects). An E ̲ X -labelled object F E ̲ X is concentrated at energy E 0 E  if F E = for all E E 0 . We write F @ E 0 for such an object.
Definition 7
(Model B: general labelled objects). A general labelled object is an arbitrary F E ̲ X . Fibrewise reasoning proceeds by applying the pullback functors ( ) E of Definition 5. By Theorem 1, this is equivalent to reasoning componentwise in the family ( F E ) E E .
Remark 3
(No branch mixing as a theorem). In Model B, any morphism F G in the labelled topos is equivalent to a family of morphisms F E G E . Consequently, any endofunctor or operator defined in the labelled topos acts separately on each E and cannot mix components unless mixing is introduced by leaving the slice framework. This is the precise form of “no branch mixing.”

2.6. Stacky/Local-System Labels: Covers and Monodromy

The constant-label model treats labels as a trivial discrete family. To encode analytic continuation and monodromy, labels must instead vary over X .
Definition 8
(Étale label object (cover of branches)). A label object E over X is calledétale if it is represented by an étale (or, in the present terminology, energy-smooth and unramified) morphism p : X ˜ X in B such that fibres are discrete sets of branches/sheets.
Proposition 2
(Labelled objects on an étale label are sheaves on the total space). If E is represented by a morphism p : X ˜ X , then there is an equivalence of toposes
( Sh ( B ) / X ) / E Sh ( B ) / X ˜ .
Under this equivalence, O X E -modules correspond to O X ˜ -modules.
Proof. 
This is the standard identification of a slice topos over a representable object with the topos over the representing object, using the fact that slicing is associative: ( E / X ) / X ˜ E / X ˜ for any topos E . Here E = Sh ( B ) .    □
Remark 4
(Monodromy as descent data). If p : X ˜ X is a covering with deck group Γ, then objects on X with branch data correspond to Γ-equivariant objects on X ˜ . In particular, a multivalued analytic function on X becomes a single-valued function on X ˜ , and monodromy is encoded by the Γ-action rather than ignored. Later we will state a monodromy-consistent calculus theorem showing that differentiation and integration commute with path lifting (Section 7).
Remark 5
(Stacks of labels (optional)). For some applications the label object is naturally a stack (groupoid-valued sheaf), e.g. when branches are identified up to symmetry. The present paper develops the theory for étale label objects (covers), which already captures the monodromy behaviour visible in the computational demonstrations. Extending from étale covers to stacks is formal: replace E by a groupoid object and work in the corresponding 2-topos of stacks.

2.7. Energy-Smooth Morphisms

We will need a base-change hypothesis that guarantees Kähler differentials and infinitesimal neighbourhoods commute with pullback.
Definition 9
(Energy-flat and energy-smooth). A morphism u : X X in B isenergy-flatif pullback along u is exact on O X -modules. It isenergy-smoothif, furthermore, the canonical map
u * Ω X / E 1 Ω X / E 1
is an isomorphism and the formation of the first infinitesimal neighbourhood of the diagonal commutes with pullback along u.

3. The Real Differential Functor (Via First Jets)

3.1. First Infinitesimal Neighbourhood of the Diagonal and Principal Parts

Let Y be an object of B equipped with a structure sheaf O Y , and assume Y is energy-smooth over E (in the sense appropriate to the geometric model: schemes, manifolds, etc.). Write Y × E Y for the fibre product over Spec E , and let
Δ Y : Y Y × E Y
be the diagonal. Let I O Y × E Y be the ideal sheaf of Δ Y .
Definition 10
(First infinitesimal neighbourhood). The first infinitesimal neighbourhood of the diagonal is the closed subobject
( Y × E Y ) [ 1 ] : = Spec O Y × E Y / I 2 ,
with the induced projections
p 1 , p 2 : ( Y × E Y ) [ 1 ] Y .
Definition 11
(Sheaf of principal parts of order 1). Define the sheaf of principal parts of order 1 on Y to be
P Y / E 1 : = p 1 * O ( Y × E Y ) [ 1 ] .
It is naturally an O Y -algebra (via p 1 ), and carries a canonical augmentation ϵ : P Y / E 1 O Y induced by restriction along the diagonal.
Remark 6.
The construction of P Y / E 1 and the exact sequences below are standard (principal parts / jets / differentials); see [1], IV or [2], II §8.

3.2. The First Jet Functor

Definition 12
(First jet functor). Let F be an O Y -module. Define
J Y / E 1 ( F ) : = p 1 * p 2 * F P Y / E 1 O Y F .
Proposition 3
(Fundamental exact sequence for first jets). There is a natural short exact sequence of O Y -modules
0 Ω Y / E 1 O Y F ι J Y / E 1 ( F ) ϵ F 0 ,
where ϵ is induced by restriction along the diagonal.
Reference. 
See ([2] II, §8) or ([1] IV).    □

3.3. Real Differentiation on a Labelled Base

We now specialise to the context of Section 2. There are two natural choices of “base” on which calculus is performed:
  • the geometric base X , yielding ordinary (unlabelled) calculus on X ;
  • a label base E X (constant family E ̲ X or an étale cover X ˜ X ), yielding labelled calculus.
Since the labelled topos ( Sh ( B ) / X ) / E is equivalent to Sh ( B ) / E when E is representable (Proposition 2), it is natural to perform jets on the label base itself.
Definition 13
(Real Differential functor on a labelled base). Assume the label object E is representable by an energy-smooth morphism p : E X (examples: E = E ̲ X or E = X ˜ étale over X ). Write O E for the induced structure sheaf.
Define theReal Differential functor on E  to be the endofunctor
D E : O E - Mod O E - Mod , D E ( F ) : = J E / E 1 ( F ) .
The augmentation ϵ : J E / E 1 ( F ) F is called thetag-preserving counit.
Remark 7
(Tag preservation is automatic). Working on E ensures that all constructions are intrinsically label-aware: a section of O E is already a label-tagged scalar. The map ϵ does not change the basepoint of E , hence does not change the label. In the constant-label case E = E ̲ X , this recovers “no branch mixing” in the strongest possible form: D acts componentwise on energy indices (Section 2.4).

3.4. Constant Energy Labels: E = E ̲ X

Let E = E ̲ X E E 1 X as in Definition 4. Then E is a disjoint union of copies of X , and O E is the corresponding coproduct of copies of O X .
Proposition 4
(Componentwise differentiation in the constant-label model). Under the equivalence
( Sh ( B ) / X ) / E ̲ X E E ( Sh ( B ) / X )
of Theorem 1, the functor D E ̲ X corresponds to the product of jet functors on X :
D E ̲ X E E J X / E 1 .
Equivalently, for O E ̲ X -modules F ,
D E ̲ X ( F ) E J X / E 1 ( F E )
naturally in F and in E E .
Proof. 
Since E ̲ X is a disjoint coproduct of copies of 1 X , the labelled base E ̲ X is a disjoint union of copies of X . Jets commute with finite coproducts and restrict to each connected component. Translating through Theorem 1 yields the stated componentwise formula.    □

3.5. Étale Label Covers: E = X ˜ X

Let p : X ˜ X be an étale label cover encoding branches (e.g. a Riemann surface of analytic continuations). Then labelled objects over X are equivalently objects over X ˜ (Proposition 2), and Real differentiation is simply differentiation on X ˜ .
Proposition 5
(Label-cover Real differentiation). For an étale label cover p : X ˜ X , the labelled Real Differential functor is
D X ˜ = J X ˜ / E 1 on O X ˜ - Mod .
In particular, deck transformations act by automorphisms commuting with D .
Proof. 
This is immediate from Definition 13 applied to E = X ˜ . Deck transformations are automorphisms of X ˜ over X and hence preserve the diagonal and its infinitesimal neighbourhood, commuting with jet formation.    □

3.6. Base Change for First Jets

Proposition 6
(Base change for first jets). Let u : Y Y be an energy-smooth morphism. Then for any O Y -module F there is a canonical isomorphism of O Y -modules
u * J Y / E 1 ( F ) J Y / E 1 ( u * F ) ,
functorial in F .
Proof sketch 
Form the cartesian square for diagonals and their first infinitesimal neighbourhoods:
Preprints 191283 i017
Energy-smoothness ensures formation of ( Y × E Y ) [ 1 ] commutes with base change, and that Ω 1 behaves well. Since p 1 is affine in the usual geometric settings (schemes, ringed spaces), pushforward along p 1 commutes with pullback along u on quasi-coherent modules. Apply this to p 2 * F to obtain the desired isomorphism. See [1], IV or [2], II, §8 for details in the scheme case.    □
Corollary 2
(Base change for Real differentiation). Under the hypotheses of Proposition 6, the Real Differential functor commutes with base change:
u * D Y ( F ) D Y ( u * F ) .
In particular, in the constant-label model this yields fibrewise identities ( D E ̲ X F ) E D X ( F E ) .

4. Directional Derivatives and Universal Identities

4.1. Directional Derivatives on Scalar Functions

Let Y be an energy-smooth base (either X , E ̲ X , or an étale label cover X ˜ ). The universal derivation
d : O Y Ω Y / E 1
is an E -linear derivation. If v Γ ( T Y / E ) is a vector field (i.e. an O Y -linear map Ω Y / E 1 O Y ), define the directional derivative of a scalar f O Y by contraction.
Definition 14
(Directional derivative). Let v Γ ( T Y / E ) . Define
v ( f ) : = ι v ( d f ) O Y .
Proposition 7
(Leibniz rule). For all f , g O Y ,
v ( f g ) = v ( f ) g + f v ( g ) .
Proof. 
This is immediate from d ( f g ) = f d g + g d f and O Y -linearity of ι v .    □

4.2. Polynomial Chain Rule (Always Valid)

Proposition 8
(Polynomial chain rule). Let F ( t ) E [ t ] and a O Y . Then
v F ( a ) = F ( a ) v ( a ) .
Proof. 
Expand F ( t ) = k c k t k and use repeated Leibniz to obtain d ( a k ) = k a k 1 d a , hence d ( F ( a ) ) = F ( a ) d a . Contract with v.    □

4.3. Upgraded Chain Rule: Smooth/Analytic Functional Calculus

The polynomial chain rule is universal in the purely algebraic setting. To state a chain rule for genuinely smooth or analytic functions, one must assume that O Y is equipped with a corresponding functional calculus.
Definition 15
( C 1 -functional calculus (one variable)). Assume E = R . A  C 1 -functional calculuson a sheaf of commutative R -algebras O Y is an assignment that to each C 1 function φ : R R associates a sheaf morphism
φ O : O Y O Y ,
such that:
(a)
φ φ O respects composition and identities: ( φ ψ ) O = φ O ψ O and ( id ) O = id ;
(b)
for polynomials F ( t ) R [ t ] , one has F O ( a ) = F ( a ) under the algebra structure.
A C -functional calculus is defined similarly with φ C ( R ) .
Remark 8.
A sheaf of C -rings in the sense of Dubuc–Moerdijk–Reyes (or in C -algebraic geometry) canonically carries such a calculus, and the chain rule below is built into that theory. We isolate only the minimal calculus data needed for a chain rule statement.
Theorem 2
(Smooth chain rule under functional calculus). Assume E = R and that O Y carries a C 1 -functional calculus in the sense of Definition 15. Let φ C 1 ( R ) and a O Y . Then for any vector field v Γ ( T Y / R ) ,
v φ O ( a ) = ( φ ) O ( a ) v ( a ) .
Proof. 
The identity is local on Y. In the standard geometric models motivating this paper (smooth manifolds, analytic spaces), φ O ( a ) is literally the composition φ a , and the formula is the classical chain rule. Abstractly, one can reduce to the polynomial case by approximating φ by polynomials on compacta and using continuity of derivations, provided the functional calculus is continuous in the appropriate topology. Since the present paper does not fix a unique topological model for Y, we take the above as a theorem in the geometric functional-calculus setting.    □
Remark 9
(Multi-variable version). If O Y carries an n-variable calculus φ : R n R (as in a C -ring), then for a 1 , , a n O Y ,
d φ ( a 1 , , a n ) = i = 1 n φ x i ( a 1 , , a n ) d a i ,
and hence the corresponding directional derivative identity holds after contraction with v.
Remark 10
(Fibrewise/labelwise validity). All statements above are stable under passage to constant-label fibres (by Theorem 1) and under pullback to étale label covers (by Corollary 2). Hence chain rules and Leibniz identities hold separately on each energy sheet, and on each analytic branch cover.

5. Higher Jets and Differential Operators

5.1. Higher Infinitesimal Neighbourhoods and Principal Parts

Let Y be an energy-smooth base over E (in the sense of Definition 9), with structure sheaf O Y .
Let Δ Y : Y Y × E Y be the diagonal, with ideal sheaf I O Y × E Y .
Definition 16
(jth infinitesimal neighbourhood). For j 0 , the jth infinitesimal neighbourhood of the diagonalis
( Y × E Y ) [ j ] : = Spec O Y × E Y / I j + 1 ,
with induced projections p 1 [ j ] , p 2 [ j ] : ( Y × E Y ) [ j ] Y .
Definition 17
(Principal parts of order j). Define the principal parts of order j to be the sheaf of O Y -algebras
P Y / E j : = ( p 1 [ j ] ) * O ( Y × E Y ) [ j ] .
There is a canonical augmentation ϵ [ j ] : P Y / E j O Y induced by restricting along Δ Y .
Definition 18
(Jets of order j). Let F be an O Y -module. Define thejth jet module
J Y / E j ( F ) : = ( p 1 [ j ] ) * ( p 2 [ j ] ) * F P Y / E j O Y F .
Proposition 9
(Exact sequences and associated graded in the smooth case). Assume Y is smooth over E (in the sense that Kähler differentials are locally free and the diagonal is a regular immersion). Then J Y / E j ( F ) admits a natural filtration whose graded pieces satisfy
gr J Y / E j ( F ) i = 0 j Sym i Ω Y / E 1 O Y F .
Moreover there are short exact sequences
0 Sym j Ω Y / E 1 F J Y / E j ( F ) J Y / E j 1 ( F ) 0 .
Reference. 
These are classical properties of principal parts and jets in the smooth case; see ([1] IV) or ([2] II, §8).    □

5.2. Differential Operators and the Jet Adjunction

We recall the standard intrinsic definition of differential operators.
Definition 19
(Differential operators of order j ). Let F , G be O Y -modules. Define inductively:
  • Diff Y / E 0 ( F , G ) : = Hom ̲ O Y ( F , G ) ;
  • for j 1 , an E -linear map D : F G is of order j if for all local sections a O Y , the commutator
    [ D , a ] : F G , [ D , a ] ( s ) : = D ( a s ) a D ( s ) ,
    is a differential operator of order j 1 .
Write Diff Y / E j ( F , G ) for the sheaf of such operators.
Theorem 3
(Jet–differential operator adjunction). For all O Y -modules F , G and all j 0 there is a natural isomorphism
Hom ̲ O Y J Y / E j ( F ) , G Diff Y / E j ( F , G ) ,
functorial in F and G .
Reference. 
This is standard; see ([1] IV). Concretely, J j ( F ) is the universal recipient for order j Taylor expansions, and P j is the universal algebra of order-j principal parts.    □

5.3. Labelwise Behaviour of Higher Jets

Proposition 10
(Constant-label case: componentwise jets). In the constant-label model E = E ̲ X , formation of J j is componentwise: for any O E ̲ X -module F ,
J E ̲ X / E j ( F ) E J X / E j ( F E )
naturally in E E .
Proof. 
Same argument as Proposition 4, using that E ̲ X is a disjoint union of copies of X and jets restrict to components.    □
Proposition 11
(Étale-label case: jets commute with pullback). Let p : X ˜ X be an étale label cover. Then for any O X -module F on X ,
p * J X / E j ( F ) J X ˜ / E j ( p * F ) ,
and the identification is compatible with the jet adjunction (Theorem 3).
Proof. 
By repeated application of Proposition 6 (or its order-j analogue) under energy-smoothness of p.    □

6. Infinite Jets, the Jet Comonad, and the coKleisli Calculus

6.1. Infinite Jets

Finite jets control differential operators of bounded order. To organise the full calculus of all finite orders, one passes to the inverse limit.
Definition 20
(Infinite principal parts and infinite jets). Assume the inverse limits exist in the chosen category of modules (e.g. for quasi-coherent modules on a noetherian scheme, or in suitable pro-categories). Define
P Y / E : = lim j P Y / E j , J Y / E ( F ) : = lim j J Y / E j ( F ) .
We write J ( F ) when Y / E is understood.
Remark 11
(Pro-object caveat). If inverse limits do not exist strictly in O Y -modules, the correct replacement is the pro-object { J j ( F ) } j 0 . All statements below may be interpreted in that pro-setting. To keep notation light, we state results as if the limit exists.

6.2. Coalgebra and Comonad Structure

A key structural fact is that infinite principal parts carry a canonical (coassociative, counital) coalgebra structure, encoding “composition of infinitesimal displacements.”
Proposition 12
(Coalgebra structure on P ). There are canonical maps of O Y -modules
ϵ : P Y / E O Y , Δ : P Y / E P Y / E O Y P Y / E ,
such that ( P Y / E , Δ , ϵ ) is a counital, coassociative coalgebra object in O Y -modules.
Construction sketch 
The counit ϵ is induced by restriction along the diagonal. The comultiplication Δ is induced by the inclusion of the formal groupoid of infinitesimal neighbourhoods: informally, the ( ) -neighbourhood of the diagonal in Y × Y composes via the middle factor in Y × Y × Y . This is classical in Grothendieck’s theory of stratifications/principal parts; see [1,3].    □
Theorem 4
(Jet comonad). Define the endofunctor J : O Y - Mod O Y - Mod by J ( F ) = P O Y F . Then with counit
ε F : J ( F ) F
induced from ϵ : P O Y , and comultiplication
δ F : J ( F ) J ( J ( F ) )
induced from Δ, the triple ( J , ε , δ ) is a comonad on O Y -modules.
Remark 12
(Relation to connections and stratifications). Coalgebras over the jet comonad (i.e. J -comodules) are precisely stratified modules. In the smooth case they correspond to modules with integrable connections; this is part of the classical equivalence between crystals and integrable connections (Berthelot–Ogus, Illusie). We will not use this equivalence in proofs, but it motivates why J governs differential operator calculus.

6.3. CoKleisli Category and Differential Operators

Definition 21
(CoKleisli morphisms). The coKleisli category of the comonad J has the same objects as O Y -modules, and morphisms
Hom coKl ( J ) ( F , G ) : = Hom O Y J ( F ) , G .
Composition is defined by the comonad structure: for f : J ( F ) G and g : J ( G ) H ,
g f : = g J ( f ) δ F .
Theorem 5
(CoKleisli morphisms are differential operators). Assume Y is energy-smooth over E . Then for O Y -modules F , G there is a natural identification
Hom ̲ O Y J ( F ) , G j 0 Diff Y / E j ( F , G ) ,
compatible with composition: coKleisli composition corresponds to composition of differential operators.
Proof sketch 
By Theorem 3, morphisms out of J j ( F ) correspond to order j differential operators. Passing to the inverse limit over j, morphisms out of J ( F ) correspond to compatible families, i.e. operators of some finite order. Compatibility of composition follows from the coalgebra/comonad map Δ encoding composition of infinitesimal expansions.    □

6.4. Labelwise Jet Comonad

All constructions in this section apply verbatim to the labelled bases Y = E ̲ X and Y = X ˜ (an étale label cover). In the constant-label case the comonad decomposes componentwise over E E ; in the étale-label case it is the usual jet comonad on X ˜ and is automatically monodromy-aware.

7. Real de Rham Theory and Monodromy-Consistent Calculus

7.1. De Rham Complex on a Labelled Base

Let Y denote an energy-smooth base (again: X , E ̲ X , or an étale label cover X ˜ ). Define
Ω Y / E p : = O Y p Ω Y / E 1 , p 0 .
The universal derivation induces the de Rham differential
d : Ω Y / E p Ω Y / E p + 1 , d 2 = 0 .
Remark 13
(Fibrewise validity). In the constant-label model, Ω E ̲ X / E decomposes into the product of Ω X / E over energies. On an étale label cover, Ω X ˜ / E is the usual de Rham complex on the cover, encoding branch data intrinsically.

7.2. Integration Along Labelled Paths (Geometric Models)

To speak about integration, we assume a geometric model in which Y carries a notion of smooth curves and integration of 1-forms (e.g. Y is a smooth manifold or complex analytic curve space). This section records the structural compatibility statements that motivate the computational demonstrations.
Definition 22
(Pathwise integral). Let γ : [ 0 , 1 ] Y be a C 1 path and ω Γ ( Y , Ω Y / E 1 ) . Define
γ ω : = 0 1 ω γ ( t ) ( γ ˙ ( t ) ) d t ,
whenever the right-hand side is defined in the ambient model.
Theorem 6
(Fundamental theorem on a labelled base). Let f Γ ( Y , O Y ) and let γ : [ 0 , 1 ] Y be a C 1 path. Then
γ d f = f ( γ ( 1 ) ) f ( γ ( 0 ) ) .
Proof. 
This is the classical fundamental theorem of calculus in the geometric model. The point of the labelled framework is that this identity takes place on the labelled base Y, so endpoints lie on the same sheet.    □

7.3. Monodromy-Consistent Calculus on an Étale Label Cover

Let p : X ˜ X be an étale label cover capturing analytic branches.
Definition 23
(Lifted path). Let γ : [ 0 , 1 ] X be a path and let x ˜ 0 X ˜ satisfy p ( x ˜ 0 ) = γ ( 0 ) . Aliftof γ starting at x ˜ 0 is a path γ ˜ : [ 0 , 1 ] X ˜ such that p γ ˜ = γ and γ ˜ ( 0 ) = x ˜ 0 .
Theorem 7
(Monodromy-consistent fundamental theorem). Let f Γ ( X ˜ , O X ˜ ) be a single-valued branch function on the cover, and let γ : [ 0 , 1 ] X be a path with a fixed lift γ ˜ : [ 0 , 1 ] X ˜ . Then
γ ˜ d f = f ( γ ˜ ( 1 ) ) f ( γ ˜ ( 0 ) ) .
If γ is a loop based at x 0 and γ ˜ ( 0 ) = x ˜ 0 , then γ ˜ ( 1 ) = g · x ˜ 0 for a unique deck transformation g; thus analytic continuation of f along γ is given by f ( g · x ˜ 0 ) .
Remark 14
(Connection to the code demonstrations). The scripts in Appendix A.2 implement this theorem in the basic case f ( z ) = z along loops around 0, contrasting lifted-path integration with principal-branch evaluation, which exhibits discontinuities at branch cuts. Appendix A.3 performs the analogous construction for complete elliptic integrals via a Gauss–Manin ODE system, numerically capturing logarithmic monodromy around m = 1 .

8. Energy Laplacian and Branchwise PDE

8.1. Gradient, Divergence, Laplacian on a Labelled Base

Assume Y carries a Riemannian metric (or, in algebraic settings, a suitable metric structure) on T Y / E inducing an isomorphism : Ω Y / E 1 T Y / E . Define for f O Y :
f : = ( d f ) , div ( V ) : = tr ( V ) , Δ E f : = div ( f ) ,
where on vector fields is the Levi–Civita connection when available.
Proposition 13
(Labelwise Laplacian is fibrewise). In the constant-label model Y = E ̲ X , Δ E acts componentwise on energy fibres. In the étale-label model Y = X ˜ , Δ E acts on each branch on the cover and is compatible with deck transformations.

8.2. Heat Equation Without Branch Mixing

Given u ( t , ) Γ ( Y , O Y ) , consider the heat equation
t u = Δ E u .
In the constant-label model this is a family of independent heat equations indexed by E E . On an étale branch cover it is a heat equation on the cover, hence branch-consistent by construction.

9. Controlled Energy/Label Mixing via Correspondences

The labelled slice framework enforces no mixing. To model controlled mixing between labels (e.g. transitions between sheets), we enrich the theory by allowing operators defined along correspondences of label objects.

9.1. Label Correspondences

Definition 24
(Label correspondence). Let E be a label object over X (constant or étale). A label correspondence is a span in Sh ( B ) / X :
E s R t E .
Intuitively, R encodes allowed transitions from a source label (via s) to a target label (via t).
Definition 25
(Transfer functor induced by a correspondence). Assume s and t are representable and that pushforward on modules along t is defined in the chosen model. Define the transfer (mixing) endofunctor on O E -modules by
T R ( F ) : = t * ( s * F ) .
Remark 15
(Energy-preserving operators as the diagonal correspondence). If R = E and s = t = id E (the diagonal span), then T R id : no mixing. More generally, a mixing pattern is determined by the fibres of s and t.

9.2. Constant-Label Case: Correspondences Are Relations on Energies

When E = E ̲ X E E 1 X , a correspondence R E ̲ X × X E ̲ X amounts to a relation on the energy index set (possibly varying over X ). In this case, T R acts like a matrix of functors/operators
( T R ( F ) ) E E : ( E E ) R F E ,
with the precise sum/product depending on whether t * or t ! is used.
Example 1
(Coupled branch heat system). Let energies be discrete indices E { 1 , , N } and let κ E , E 0 be coupling constants. A controlled-mixing heat system has the form
t u E = Δ u E + E κ E , E u E .
The Laplacian term is no-mixing; the coupling term is encoded by a correspondence R allowing transitions E E with weights κ E , E .

9.3. Composition of Correspondences (Optional Categorical Structure)

Under standard hypotheses (suitable base change and projection formula), spans compose by pullback:
( E s R t E ) ( E s R t E ) : = E s π 1 R × E R t π 2 E ,
and the induced transfer functors compose accordingly. This yields a bicategory of label correspondences in which no-mixing calculus is the diagonal subtheory.

Appendix A. Computational Demonstrations and Reproducibility

Appendix A.1. Reproducibility Notes

The scripts below are included verbatim for transparency and reproducibility. They are not required for the theoretical development, but they implement:
  • monodromy-consistent differentiation along lifted paths (vs. principal-branch failure);
  • numerical branchwise heat flow with no label mixing;
  • a basic descent/gluing check enforcing label compatibility;
  • (elliptic case) numerical Gauss–Manin transport for K ( m ) and E ( m ) and detection of logarithmic monodromy around m = 1 .
  • Dependencies. The first script uses numpy and matplotlib; it optionally uses sympy. The second script uses numpy, matplotlib, and mpmath.
  • Outputs. The scripts save figures to PNG files. The appendix figures below are representative outputs (filenames as provided).

Appendix A.2. Listing: real_diff_demo.py

Preprints 191283 i001
Preprints 191283 i002
Preprints 191283 i003
Preprints 191283 i004
Preprints 191283 i005
Preprints 191283 i006
Preprints 191283 i007
Preprints 191283 i008
Preprints 191283 i009
Preprints 191283 i010
Preprints 191283 i011

Appendix A.3. Listing: rd_elliptic_monodromy.py

Preprints 191283 i012
Preprints 191283 i013
Preprints 191283 i014
Preprints 191283 i015
Preprints 191283 i016

Appendix B. Figures (Representative Outputs)

Remark A1.
For journal submission, we strongly recommend renaming files to stable descriptive names (e.g.sqrt_monodromy_turns1.png) and updating captions to match exact content. The figures below are included with the filenames provided.
Figure A1. Representative computational output (1) from the accompanying scripts.
Figure A1. Representative computational output (1) from the accompanying scripts.
Preprints 191283 g0a1
Figure A2. Representative computational output (2) from the accompanying scripts.
Figure A2. Representative computational output (2) from the accompanying scripts.
Preprints 191283 g0a2
Figure A3. Representative computational output (3) from the accompanying scripts.
Figure A3. Representative computational output (3) from the accompanying scripts.
Preprints 191283 g0a3
Figure A4. Representative computational output (4) from the accompanying scripts.
Figure A4. Representative computational output (4) from the accompanying scripts.
Preprints 191283 g0a4
Figure A5. Representative computational output (5) from the accompanying scripts.
Figure A5. Representative computational output (5) from the accompanying scripts.
Preprints 191283 g0a5
Figure A6. Representative computational output (6) from the accompanying scripts.
Figure A6. Representative computational output (6) from the accompanying scripts.
Preprints 191283 g0a6
Figure A7. Representative computational output (7) from the accompanying scripts.
Figure A7. Representative computational output (7) from the accompanying scripts.
Preprints 191283 g0a7

References

  1. A. Grothendieck and J. Dieudonné. Éléments de Géométrie Algébrique IV. Publ. Math. IHÉS, 1964–1967.
  2. R. Hartshorne. Algebraic Geometry. Graduate Texts in Mathematics, Vol. 52. Springer, 1977.
  3. L. Illusie. Complexe Cotangent et Déformations I. Lecture Notes in Mathematics, Vol. 239. Springer, 1971.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

Disclaimer

Terms of Use

Privacy Policy

Privacy Settings

© 2026 MDPI (Basel, Switzerland) unless otherwise stated