Time-Periodic Cahn-Hilliard-Gurtin System

• Time-Periodic Cahn-Hilliard-Gurtin System is a class of mixed-order, time-periodic PDEs that extend classical phase-field models with higher-order spatial and temporal derivatives.
• It introduces a novel mixed-order complementing condition that refines traditional criteria like the Lopatinskiĭ–Shapiro conditions for boundary operator admissibility.
• The framework employs Newton-polygon Sobolev spaces and factorization techniques to achieve explicit function space estimates and well-posedness in half-space domains.

The time-periodic Cahn-Hilliard-Gurtin (CHG) system describes a class of mixed-order, time-periodic partial differential equations (PDEs) posed on the half-space, where the unknowns and data exhibit time-periodicity. The system generalizes classical CH and Gurtin-type phase-field equations by incorporating higher-order mixed derivatives and allowing for highly general boundary conditions, including those expressing Neumann and complex trace-coupling. Recent mathematical analysis has revealed that classical boundary regularity and solvability criteria, such as the Lopatinskiĭ–Shapiro conditions, are insufficient for ensuring well-posedness in this mixed-order, time-periodic context. Instead, a novel "mixed-order complementing condition" is required, yielding improved criteria for maximal regularity and explicit function space estimates in strong Sobolev norms (Neuttiens et al., 29 Dec 2025).

1. System Definition and Time-Periodicity

The time-periodic CHG system is posed for unknowns $u=(u_1,u_2)$ in the domain $(t,x)\in\T\times\R^n_+$, with

$\T = \R/T\Z,\quad \R^n_+ = \{x\in\R^n\,:\,x_n>0\}.$

The evolution equations read

$\begin{cases} \partial_tu_1-\Delta_xu_2=f_1,\ \Delta_xu_1-\partial_tu_1+u_2=f_2, \end{cases} \qquad \text{in }\T\times\R^n_+,$

subject to boundary operators

$B_1u = g_1,\quad B_2u = g_2 \qquad \text{on }\T\times\R^{n-1},$

which may involve derivatives up to third order in $x_n$. Time-periodicity is imposed:

$$u_i(t+T,x) = u_i(t,x),\quad f_i(t+T,x) = f_i(t,x),\quad g_i(t+T,x') = g_i(t,x').$$

This structure models processes that attain periodic steady states under cyclic forcing or thermal conditions.

2. Boundary Operators and Complementing Conditions

General boundary operators $B_i$ are defined by

$$B_i u = \sum_{j=0}^3\,op\bigl[b_j^{(i)}(\tau,\xi')\bigr]\,Tr_j\,u,\quad Tr_j\,u = \partial_n^j u\big|_{x_n=0},\quad i=1,2,$$

with $b_j^{(i)}$ smooth in $(\tau,\xi')$, the dual variables to $t$ and $x'$. For mixed-order systems, the well-posedness of boundary problems hinges on a new complementing condition: the invertibility of an extended boundary matrix $\CB^{(m)}(\tau,\xi')$ (see Def. 3.11), with rows involving both the boundary symbols and factors from the interior symbol $D^+$. Each row and column of $\CB^{(m)}$ is assigned specific upper and lower orders according to the Agmon–Douglis–Nirenberg mixed-order formalism, adapted here to Newton polygons. Explicitly, for nonzero $(\tau, \xi')$, invertibility of $\CB^{(m)}(\tau,\xi')$ is required. This condition strictly refines the classical Lopatinskiĭ–Shapiro condition, which examines independence modulo $D^+$ only; the mixed-order criterion demands full invertibility and order alignment.

3. Function Spaces and Maximal Regularity Theorems

Analysis of the time-periodic CHG system is carried out in Newton-polygon Sobolev spaces $H^\mu_\perp(\T\times\R^n_+)$, which encode the mixed-order character of the symbol and use the purely oscillatory, mean-free subspace with respect to time. For instance,

$H^{\mu}_\perp(\T\times\R^n_+) = \left\{u:\exists\,U\in H^\mu_\perp(\T\times\R^n),\,U|_{\R^n_+}=u\right\}.$

The principal result (Thm 1.4/Thm 1.1) asserts that the trace-coupled operator

$$S(u_1, u_2) = \left(\partial_tu_1-\Delta u_2,\, \Delta u_1-\partial_tu_1+u_2,\, Tr_1u_1,\, Tr_1u_2\right)$$

defines an isomorphism

$S : \E \xrightarrow{\sim} \F \times \G,$

where \begin{align*} \E_1 &= \overline{H}^{(0,3)}_\perp \cap \overline{H}^{(1,1)}_\perp, \ \E_2 &= \overline{H}^{(0,2)}_\perp, \ \E &= \E_1 \times \E_2, \ \F_1 &= L^2_\perp, \ \F_2 &= \overline{H}^{(0,1)}_\perp, \ \F &= \F_1 \times \F_2, \ \G_1 &= H^{(0,\frac{3}{2})}_\perp \cap H^{(\frac{3}{4},0)}_\perp, \ \G_2 &= H^{(0,\frac{1}{2})}_\perp, \ \G &= \G_1 \times \G_2. \end{align*} The a priori estimate follows:

$\|(u_1, u_2)\|_{\E}\lesssim \|(f_1, f_2)\|_{\F} + \|(g_1, g_2)\|_{\G}.$

4. Symbol Analysis, Factorization, and Solution Strategy

The interior symbol for the CHG system is

$$D(\tau, \xi) = i\tau + |\xi|^2\left(i\tau + |\xi|^2\right),$$

with Newton polygon $\{(0,0), (4,0), (2,1), (0,1)\}$, a structure necessary for $\NN$-ellipticity. The factorization $D = D^- D^+$, based on Paley–Wiener theory, allows construction of right-inverses that respect support in $x_n$, crucial for boundary problems on half-spaces. On the full space $\R^n$, solutions are characterized by the isomorphism

$op[D]:H^{\mu_D+\nu}_\perp\to H^\nu_\perp,$

while on the half-space, the trace theorem (Thm 3.15) provides control over the boundary data within Newton-polygon spaces on $\Gamma$. The boundary-value problem is resolved by solving the boundary system via invertibility of the extended $\CB$ matrix, then extending solutions inward using Poisson-type kernel techniques and the factor $D^+$.

5. Boundary Condition Admissibility and Counterexamples

The admissibility of boundary conditions is dictated by the mixed-order complementing condition. For Neumann-type conditions, e.g.

$B_1u = \partial_nu_1,\qquad B_2u = \partial_nu_2,$

the extended $\CB$ satisfies the mixed-order complementing condition with $m=0$, guaranteeing well-posedness. In contrast, classical Dirichlet boundary conditions,

$B_1u = Tr_0u,\qquad B_2u = Tr_1u,$

are generally non-admissible: even when the Lopatinskiĭ–Shapiro condition is met, Theorem 1.8 demonstrates existence of $f\in L^2_\perp$ for which the problem admits no solution in $H^{(1,2)}_\perp \cap H^{(0,4)}_\perp$. This failure motivates the introduction and necessity of the mixed-order criterion.

BC type	Extended $\CB$ satisfies complementing?	Classical Lopatinskiĭ–Shapiro sufficient?
Neumann	Yes	Yes
Dirichlet	No	May fail (counterexample exists)

6. Interior and Boundary Symbols; Explicit Matrix Construction

The detailed analysis of interior and boundary symbols is essential for mixed-order ellipticity. The interior symbol (as determinant) is

$$D(\tau, \xi', \xi_n) = i\tau + (|\xi'|^2+\xi_n^2)\left(i\tau+|\xi'|^2+\xi_n^2\right),$$

and upon factorization,

$$D^+(\tau, \xi', \xi_n) = (\xi_n-\rho_1^+)(\xi_n-\rho_2^+) = (-1)\xi_n^2 + i(\rho_1^++\rho_2^+)\xi_n + \rho_1^+\rho_2^+,$$

where $\rho_j^+(\tau,\xi')$ are the two roots in $\Im\xi_n>0$. The extended boundary matrix for $m = ord\nu$ is

$\CB^{(m)}(\tau,\xi') = \begin{pmatrix} b_0^{(1)} & b_1^{(1)} & b_2^{(1)} & b_3^{(1)} & 0 & \cdots \ b_0^{(2)} & b_1^{(2)} & b_2^{(2)} & b_3^{(2)} & 0 & \cdots \ \rho_1^+\rho_2^+ & \rho_1^+ + \rho_2^+ & -1 & 0 & 0 & \cdots \ 0 & \rho_1^+\rho_2^+ & \rho_1^+ + \rho_2^+ & -1 & \cdots \ \vdots & \vdots & \vdots & \vdots & \ddots \end{pmatrix}.$

Verification of complementing involves showing $\det\CB^{(m)}(\tau,\xi')\neq 0$ for all nonzero $(\tau, \xi')$, and correct order alignment in all entries, confirming $\NN$-elliptic character. Classical independence modulo $D^+$ yields only partial guarantees.

7. Implications and Scope

The advanced framework for the time-periodic Cahn-Hilliard-Gurtin system on the half-space yields a comprehensive $L^2$-theory for mixed-order time-periodic systems. It specifies sharp function space estimates, rigorously addresses well-posedness limitations of the classical theory, and provides concrete criteria for boundary data admissibility. A plausible implication is that similar mixed-order extensions may be required for other PDE systems exhibiting high-order temporal or spatial structure, periodicity, and complex boundary constraints (Neuttiens et al., 29 Dec 2025). The methodology anchors in deep symbol analysis, factorization, and precise Sobolev class geometry, setting a new standard for elliptic and parabolic boundary problem theory in complex domains.