Nonhomogeneous Cauchy–Riemann equation is a complex PDE prescribing the antiholomorphic derivative, forming a cornerstone in modern complex analysis.
It is solved in scalar and vector-valued settings using techniques like Hörmander-type L² methods, functional-analytic splitting, and weighted smooth space approaches.
Extensions to nonlinear cases and geometric formulations on Riemann surfaces enable applications in Hardy theory, boundary value problems, and elliptic regularity.
The nonhomogeneous Cauchy–Riemann equation is the first-order complex PDE obtained by prescribing the antiholomorphic derivative of an unknown function. In the planar scalar setting, with z=x+iy, it is written
∂u=f,∂:=21(∂1+i∂2),
and in a vector-valued setting it becomes
∂Eu:=21(∂1E+i∂2E)u=f,
where u takes values in a locally convex Hausdorff space E (Kruse, 2019). Across recent literature, the equation appears in several analytically distinct regimes: weighted smooth Fréchet and locally convex spaces, Lp-Dolbeault theory on complex manifolds, Hardy classes on the disk, compact-support problems in several variables, nonlinear equations of the form ∂zˉu=E(z,u), and geometric formulations on open Riemann surfaces with varying complex structures (Kruse, 2018).
1. Basic formulations and geometric meaning
In R2≅C, the operator
∂=21(∂x∂+i∂y∂)
annihilates holomorphic functions, so the homogeneous equation ∂u=0 is the Cauchy–Riemann system, while the nonhomogeneous equation ∂u=f,∂:=21(∂1+i∂2),0 prescribes the antiholomorphic component of the differential (Blair, 2023). For scalar-valued ∂u=f,∂:=21(∂1+i∂2),1, this is the standard planar model; for ∂u=f,∂:=21(∂1+i∂2),2-valued smooth maps ∂u=f,∂:=21(∂1+i∂2),3, the same operator is defined by iterated partial derivatives in the locally convex topology of ∂u=f,∂:=21(∂1+i∂2),4 (Kruse, 2019).
On a Riemann surface ∂u=f,∂:=21(∂1+i∂2),5, the differential decomposes as
∂u=f,∂:=21(∂1+i∂2),6
and the nonhomogeneous equation takes the intrinsic form
∂u=f,∂:=21(∂1+i∂2),7
where ∂u=f,∂:=21(∂1+i∂2),8 is a ∂u=f,∂:=21(∂1+i∂2),9-form. In a ∂Eu:=21(∂1E+i∂2E)u=f,0-holomorphic coordinate ∂Eu:=21(∂1E+i∂2E)u=f,1, this is exactly ∂Eu:=21(∂1E+i∂2E)u=f,2 when ∂Eu:=21(∂1E+i∂2E)u=f,3 (Forstneric, 19 Aug 2025). This coordinate realization is one reason the planar equation remains the local model even in geometric settings.
A related real-variable formulation arises in generalized Cauchy–Riemann systems
∂Eu:=21(∂1E+i∂2E)u=f,4
which are linked to Beltrami equations ∂Eu:=21(∂1E+i∂2E)u=f,5 by taking ∂Eu:=21(∂1E+i∂2E)u=f,6. In the special case ∂Eu:=21(∂1E+i∂2E)u=f,7, the classical Cauchy–Riemann equations are recovered. The same framework is tied to the divergence-form ∂Eu:=21(∂1E+i∂2E)u=f,8-harmonic equation ∂Eu:=21(∂1E+i∂2E)u=f,9 through the relation u0, where u1 is the u2-rotation matrix (Gutlyanskii et al., 2024).
The equation also admits nonlinear variants. A basic family studied in the literature is
u3
with separable special case
u4
In these forms, the equation remains elliptic but loses linear superposition, and regularity as well as explicit solvability depend strongly on the structure of the right-hand side (Coffman et al., 2014).
2. Solvability in weighted smooth and vector-valued spaces
A major modern line of work studies u5 on spaces u6 of weighted u7-smooth u8-valued functions, where growth is controlled on an exhaustion u9 by weights E0. For E1, the seminorms are
E2
and E3 is the corresponding projective limit space (Kruse, 2019). In the strip model
E4
typical weights are
E5
which control growth in the real direction uniformly on horizontal strips (Kruse, 2019).
Within this setting, the scalar operator
E6
is shown to be surjective under a combination of weight hypotheses, density conditions, and subharmonicity of E7. One of the central statements is that if Condition 3.3 with E8, Condition 4.2 with E9, and subharmonicity of Lp0 hold, then Lp1 is surjective on Lp2; in the strip case with Lp3 and
Lp4
this yields surjectivity on Lp5 (Kruse, 2019). An earlier paper derives the corresponding scalar surjectivity in weighted smooth spaces through a Hörmander-type Lp6 argument, a weakened form of weak reducibility for the holomorphic subspace, and a Mittag-Leffler patching procedure (Kruse, 2018).
The vector-valued problem is handled by functional-analytic splitting theory. A key hypothesis is that the kernel Lp7, identified with a weighted holomorphic space Lp8, has Vogt’s property Lp9. Under suitable strip and weight assumptions, including weights of the form ∂zˉu=E(z,u)0, this kernel indeed has property ∂zˉu=E(z,u)1 (Kruse, 2019). One then combines scalar surjectivity, nuclearity of ∂zˉu=E(z,u)2, and the tensor product representation
∂zˉu=E(z,u)3
to deduce surjectivity of
∂zˉu=E(z,u)4
for two large classes of targets: ∂zˉu=E(z,u)5 with ∂zˉu=E(z,u)6 Fréchet satisfying ∂zˉu=E(z,u)7, and ultrabornological (PLS)-spaces with ∂zˉu=E(z,u)8 (Kruse, 2019).
A closely related extension treats domains of the form ∂zˉu=E(z,u)9, where R2≅C0 is compact and the exhaustion is given by strips with holes R2≅C1. In that case, weighted duality for
R2≅C2
is used to prove that R2≅C3 has property R2≅C4, and hence that
R2≅C5
is surjective under the same structural assumptions on R2≅C6 (Kruse, 2019).
3. R2≅C7, compact-support, and exact-support theories
In complex manifolds, the nonhomogeneous equation is formulated distributionally on spaces of R2≅C8 forms. For a complex manifold R2≅C9, holomorphic vector bundle ∂=21(∂x∂+i∂y∂)0, and ∂=21(∂x∂+i∂y∂)1, the equation ∂=21(∂x∂+i∂y∂)2 means that the current ∂=21(∂x∂+i∂y∂)3 is represented by an ∂=21(∂x∂+i∂y∂)4 form (Laurent-Thiébaut, 2013). The ∂=21(∂x∂+i∂y∂)5-Dolbeault complex is a fine resolution of ∂=21(∂x∂+i∂y∂)6, yielding Dolbeault isomorphisms
∂=21(∂x∂+i∂y∂)7
so local solvability of ∂=21(∂x∂+i∂y∂)8 for ∂=21(∂x∂+i∂y∂)9-closed data is built into the cohomological framework (Laurent-Thiébaut, 2013).
On strictly ∂u=00-convex domains with smooth boundary, global ∂u=01 solvability is obtained by homotopy formulas. For ∂u=02, there are continuous linear operators ∂u=03 and compact operators ∂u=04 such that
∂u=05
In completely strictly ∂u=06-convex domains this leads to vanishing of ∂u=07 for ∂u=08, hence global solvability of ∂u=09 with ∂u=f,∂:=21(∂1+i∂2),00 for every ∂u=f,∂:=21(∂1+i∂2),01-closed ∂u=f,∂:=21(∂1+i∂2),02 in that range (Laurent-Thiébaut, 2013).
A sharper variant is solvability with exact support: given ∂u=f,∂:=21(∂1+i∂2),03 supported in ∂u=f,∂:=21(∂1+i∂2),04, find ∂u=f,∂:=21(∂1+i∂2),05 also supported in ∂u=f,∂:=21(∂1+i∂2),06 such that ∂u=f,∂:=21(∂1+i∂2),07. This is characterized by ∂u=f,∂:=21(∂1+i∂2),08-Serre duality. For ∂u=f,∂:=21(∂1+i∂2),09, the complexes ∂u=f,∂:=21(∂1+i∂2),10 and ∂u=f,∂:=21(∂1+i∂2),11 are dual under the integral pairing, and vanishing or Hausdorff properties of the dual cohomology groups control solvability with support constraints (Laurent-Thiébaut, 2013). In top degree ∂u=f,∂:=21(∂1+i∂2),12, an additional orthogonality condition against holomorphic dual forms is necessary and sufficient.
Compact-support solvability is also studied by explicit convolution methods in ∂u=f,∂:=21(∂1+i∂2),13. For one complex variable, compactly supported solvability of
∂u=f,∂:=21(∂1+i∂2),14
is characterized by moment conditions
∂u=f,∂:=21(∂1+i∂2),15
In higher dimensions, Amar and Mongodi construct compactly supported solutions in ∂u=f,∂:=21(∂1+i∂2),16 and on domains ∂u=f,∂:=21(∂1+i∂2),17, where ∂u=f,∂:=21(∂1+i∂2),18 is a polydisc and ∂u=f,∂:=21(∂1+i∂2),19 is the zero locus of a holomorphic function, by repeated convolution with the one-variable Cauchy kernel ∂u=f,∂:=21(∂1+i∂2),20 and by enforcing higher-dimensional structure conditions (Amar et al., 2011). Their results give ∂u=f,∂:=21(∂1+i∂2),21-bounded solution operators and compact-support control for ∂u=f,∂:=21(∂1+i∂2),22-forms under explicit integrability and moment-type conditions.
4. Boundary classes, Hardy theory, and integral representations
On the unit disk ∂u=f,∂:=21(∂1+i∂2),23, the equation ∂u=f,∂:=21(∂1+i∂2),24 can be studied inside Hardy-type spaces. For ∂u=f,∂:=21(∂1+i∂2),25,
∂u=f,∂:=21(∂1+i∂2),26
and Blair focuses on the regime ∂u=f,∂:=21(∂1+i∂2),27 with ∂u=f,∂:=21(∂1+i∂2),28, ∂u=f,∂:=21(∂1+i∂2),29 (Blair, 2023). In this setting Vekua’s representation of the second kind yields
∂u=f,∂:=21(∂1+i∂2),30
where ∂u=f,∂:=21(∂1+i∂2),31 is holomorphic and
∂u=f,∂:=21(∂1+i∂2),32
This decomposes a nonhomogeneous solution into a holomorphic Hardy component and a singular-integral correction (Blair, 2023).
The mapping properties of ∂u=f,∂:=21(∂1+i∂2),33 govern boundary regularity. If ∂u=f,∂:=21(∂1+i∂2),34 with ∂u=f,∂:=21(∂1+i∂2),35, then ∂u=f,∂:=21(∂1+i∂2),36 has radial boundary values in ∂u=f,∂:=21(∂1+i∂2),37 for ∂u=f,∂:=21(∂1+i∂2),38; if ∂u=f,∂:=21(∂1+i∂2),39, then ∂u=f,∂:=21(∂1+i∂2),40 with ∂u=f,∂:=21(∂1+i∂2),41 (Blair, 2023). Blair uses this to prove that every ∂u=f,∂:=21(∂1+i∂2),42 has a boundary value in the sense of distributions and that
∂u=f,∂:=21(∂1+i∂2),43
where ∂u=f,∂:=21(∂1+i∂2),44 are ∂u=f,∂:=21(∂1+i∂2),45-atoms and ∂u=f,∂:=21(∂1+i∂2),46 (Blair, 2023). Thus the boundary traces of nonhomogeneous solutions are represented as an atomic Hardy component plus a regular error term.
This boundary decomposition is stable under the Hilbert transform. Hoepfner–Hounie had proved continuity of the Hilbert transform on the atomic Hardy space ∂u=f,∂:=21(∂1+i∂2),47, and Blair shows that the Hilbert transform remains continuous on the boundary-value space ∂u=f,∂:=21(∂1+i∂2),48 equipped with the quasi-norm built from the atomic part and the ∂u=f,∂:=21(∂1+i∂2),49-norm of ∂u=f,∂:=21(∂1+i∂2),50 (Blair, 2023). The same framework yields solvability of Schwarz-type boundary value problems for ∂u=f,∂:=21(∂1+i∂2),51 with boundary data prescribed as Hardy-space distributions.
Integral-operator methods also underpin the weighted smooth theory. In Kruse’s weighted setting, the scalar surjectivity proof starts from Hörmander-type weighted ∂u=f,∂:=21(∂1+i∂2),52 solutions on each exhausted subdomain and then upgrades them to smooth weighted solutions through equivalence of seminorm systems and projective limit arguments (Kruse, 2018). This suggests a common structural pattern across function classes: first obtain a local or coarse right inverse, then improve regularity and compatibility through analytic or functional-analytic constraints.
5. Nonlinear equations and elliptic regularity
The nonlinear inhomogeneous equation
∂u=f,∂:=21(∂1+i∂2),53
is studied under minimal assumptions on ∂u=f,∂:=21(∂1+i∂2),54. One results framework assumes that ∂u=f,∂:=21(∂1+i∂2),55 is continuous, the classical partial derivatives ∂u=f,∂:=21(∂1+i∂2),56 exist outside a countable set, and the equation holds almost everywhere. If the right-hand side reduces to a prescribed function ∂u=f,∂:=21(∂1+i∂2),57, ∂u=f,∂:=21(∂1+i∂2),58, then ∂u=f,∂:=21(∂1+i∂2),59 for every rectangle ∂u=f,∂:=21(∂1+i∂2),60, and if ∂u=f,∂:=21(∂1+i∂2),61, then ∂u=f,∂:=21(∂1+i∂2),62 with ∂u=f,∂:=21(∂1+i∂2),63 (Coffman et al., 2014).
When ∂u=f,∂:=21(∂1+i∂2),64 is continuous, the same paper derives
∂u=f,∂:=21(∂1+i∂2),65
If ∂u=f,∂:=21(∂1+i∂2),66, then ∂u=f,∂:=21(∂1+i∂2),67 for some ∂u=f,∂:=21(∂1+i∂2),68; if ∂u=f,∂:=21(∂1+i∂2),69, then ∂u=f,∂:=21(∂1+i∂2),70, and if ∂u=f,∂:=21(∂1+i∂2),71 is smooth, then ∂u=f,∂:=21(∂1+i∂2),72 is smooth (Coffman et al., 2014). These are elliptic regularity statements for the nonhomogeneous Cauchy–Riemann system under weak initial hypotheses.
The paper also exhibits a borderline counterexample. There exists a differentiable function ∂u=f,∂:=21(∂1+i∂2),73 such that ∂u=f,∂:=21(∂1+i∂2),74 exists everywhere and is continuous, while ∂u=f,∂:=21(∂1+i∂2),75 exists everywhere but is not locally bounded near ∂u=f,∂:=21(∂1+i∂2),76 (Coffman et al., 2014). This shows that continuity of the antiholomorphic derivative does not force comparable regularity of the holomorphic derivative. A plausible implication is that regularity theories for ∂u=f,∂:=21(∂1+i∂2),77 must distinguish sharply between control of the equation itself and control of the full first-order gradient.
For separable nonlinear equations
∂u=f,∂:=21(∂1+i∂2),78
the solution set can sometimes be described explicitly. When ∂u=f,∂:=21(∂1+i∂2),79 is holomorphic and nonvanishing, one introduces a primitive ∂u=f,∂:=21(∂1+i∂2),80 and a function ∂u=f,∂:=21(∂1+i∂2),81 with ∂u=f,∂:=21(∂1+i∂2),82, obtaining locally
∂u=f,∂:=21(∂1+i∂2),83
Near simple zeros of ∂u=f,∂:=21(∂1+i∂2),84, the solution takes forms such as
∂u=f,∂:=21(∂1+i∂2),85
and for higher multiplicity zeros an implicit formula of the form
∂u=f,∂:=21(∂1+i∂2),86
is derived (Coffman et al., 2014). These formulas make explicit how the nonhomogeneous term interacts with the zero structure of the nonlinear factor.
6. Families of complex structures, Beltrami reduction, and applications
On a fixed smooth open orientable surface ∂u=f,∂:=21(∂1+i∂2),87, the equation ∂u=f,∂:=21(∂1+i∂2),88 for a family of complex structures ∂u=f,∂:=21(∂1+i∂2),89 and ∂u=f,∂:=21(∂1+i∂2),90-forms ∂u=f,∂:=21(∂1+i∂2),91 can be transformed into a nonhomogeneous Beltrami equation with respect to a fixed background structure. Writing ∂u=f,∂:=21(∂1+i∂2),92 through a Beltrami coefficient ∂u=f,∂:=21(∂1+i∂2),93, one obtains an equation of the form
∂u=f,∂:=21(∂1+i∂2),94
and solves it by Cauchy and Beurling operators ∂u=f,∂:=21(∂1+i∂2),95 and ∂u=f,∂:=21(∂1+i∂2),96 together with the inverse
∂u=f,∂:=21(∂1+i∂2),97
for small ∂u=f,∂:=21(∂1+i∂2),98 (Forstneric, 19 Aug 2025). The resulting local theorem states that if the family ∂u=f,∂:=21(∂1+i∂2),99 is of class ∂Eu:=21(∂1E+i∂2E)u=f,00, then there exists
∂Eu:=21(∂1E+i∂2E)u=f,01
such that ∂Eu:=21(∂1E+i∂2E)u=f,02 for every parameter ∂Eu:=21(∂1E+i∂2E)u=f,03 (Forstneric, 19 Aug 2025).
The main regularity feature is a gain of one derivative in the space variable and no loss of regularity in the parameter. Global solvability on the whole open surface is obtained by exhaustion with Runge domains and approximation; the paper notes that, except perhaps for the condition ∂Eu:=21(∂1E+i∂2E)u=f,04 in the global corollary, the result is “optimal” in this sense (Forstneric, 19 Aug 2025). The same solvability feeds into a sheaf-theoretic resolution
∂Eu:=21(∂1E+i∂2E)u=f,05
from which vanishing of higher cohomology ∂Eu:=21(∂1E+i∂2E)u=f,06 for ∂Eu:=21(∂1E+i∂2E)u=f,07 is deduced and then an Oka–Grauert principle for line bundles on families of open Riemann surfaces (Forstneric, 19 Aug 2025).
The weighted smooth vector-valued theory also has a parameter-dependence consequence. If ∂Eu:=21(∂1E+i∂2E)u=f,08 depends smoothly, holomorphically, or distributionally on a parameter and the relevant target space ∂Eu:=21(∂1E+i∂2E)u=f,09 belongs to the admissible classes above, then there exists a family ∂Eu:=21(∂1E+i∂2E)u=f,10 with the same type of parameter dependence satisfying
∂Eu:=21(∂1E+i∂2E)u=f,11
(Kruse, 2019). This is an explicit preservation-of-regularity statement under the solution operator.
A broader survey of generalized Cauchy–Riemann systems emphasizes that the homogeneous Beltrami and ∂Eu:=21(∂1E+i∂2E)u=f,12-harmonic theories provide the operator-theoretic backbone for nonhomogeneous equations of the schematic types
∂Eu:=21(∂1E+i∂2E)u=f,13
The survey does not treat these nonhomogeneous equations explicitly, but it clarifies how existence, representation, and regularity results for the homogeneous operators organize the corresponding inhomogeneous theories (Gutlyanskii et al., 2024).
The contemporary picture is therefore not a single solvability theorem but a family of structurally different theories. In weighted smooth spaces the equation is treated by nuclear Fréchet and (PLS) techniques; in ∂Eu:=21(∂1E+i∂2E)u=f,14 theory by Dolbeault resolutions, homotopy formulas, and duality; in Hardy classes by singular-integral representations and atomic boundary decompositions; in nonlinear settings by elliptic regularity and explicit factorization; and on families of Riemann surfaces by Beltrami reduction and parametric integral operators (Kruse, 2019).