Inhomogeneous Biharmonic Schrödinger Equation

• The inhomogeneous biharmonic Schrödinger equation is a fourth-order dispersive PDE featuring a spatially weighted nonlinearity that modifies translation invariance and concentrates effects near the origin.
• It incorporates standard and mixed-dispersion models, with clearly defined criticality regimes (mass, energy, intercritical) and conserved quantities such as mass and fourth-order energy.
• Advanced analytical frameworks using weighted Sobolev, Lorentz, and Besov spaces underpin its well-posedness, scattering theory, concentration-compactness, and boundary control approaches.

The inhomogeneous biharmonic Schrödinger equation is a class of fourth-order dispersive equations in which the biharmonic propagator is coupled to a spatially weighted nonlinearity. The canonical whole-space model is

$$i\partial_t u+\Delta^2u=\lambda |x|^{-b}|u|^\sigma u,\qquad u(0,x)=u_0(x),$$

with $u:\mathbb R\times\mathbb R^d\to\mathbb C$, $u_0$ posed in $L^2$, $H^s$, or $\dot H^s$, and $b>0$ measuring the strength of the inhomogeneity near the origin. Closely related formulations include mixed-dispersion equations of the form $$i\partial_t u-\Delta^2u+\mu\Delta u+|x|^{-b}|u|^\alpha u=0$$, energy-critical focusing models, and equations with a spatially growing factor $|x|^b$. The weight breaks translation invariance and introduces singular or unbounded coefficients, so the subject couples fourth-order Strichartz theory with weighted Sobolev inequalities, variational thresholds, virial identities, and concentration-compactness. In the broader fourth-order Schrödinger literature, the same inhomogeneous paradigm also includes initial-boundary value problems with non-homogeneous boundary data, internal control, and equations with external potentials on waveguide manifolds (Guzmán et al., 2019, Dinh et al., 2022, Alqaied et al., 10 Nov 2025, Li et al., 2020).

1. Model equations, parameters, and criticality

The standard decaying-weight equation is parameterized by the spatial dimension $d$, Sobolev index $u:\mathbb R\times\mathbb R^d\to\mathbb C$0, inhomogeneity exponent $u:\mathbb R\times\mathbb R^d\to\mathbb C$1, power $u:\mathbb R\times\mathbb R^d\to\mathbb C$2 or $u:\mathbb R\times\mathbb R^d\to\mathbb C$3, and a sign parameter $u:\mathbb R\times\mathbb R^d\to\mathbb C$4 or $u:\mathbb R\times\mathbb R^d\to\mathbb C$5 distinguishing focusing and defocusing conventions. In whole-space Cauchy problems, the principal linear flow is $u:\mathbb R\times\mathbb R^d\to\mathbb C$6, while mixed-dispersion variants replace it by $u:\mathbb R\times\mathbb R^d\to\mathbb C$7. The mixed-dispersion model

$u:\mathbb R\times\mathbb R^d\to\mathbb C$8

is studied for $u:\mathbb R\times\mathbb R^d\to\mathbb C$9, $u_0$0, and $u_0$1, with the weight remaining singular at the origin but decaying at infinity (Guzmán et al., 2019).

For the decaying-weight equation, the scaling-critical Sobolev index is

$u_0$2

This yields the standard trichotomy: mass-critical when $u_0$3, equivalently $u_0$4; energy-critical when $u_0$5, equivalently $u_0$6 for $u_0$7; and intercritical, or mass-supercritical and energy-subcritical, when $u_0$8. In the mixed-dispersion setting the same regime is encoded by

$u_0$9

and the range $L^2$0 is exactly the mass-supercritical, energy-subcritical range. When $L^2$1, the equation is scale-invariant; when $L^2$2, the repulsive harmonic perturbation breaks the scaling symmetry while preserving the same basic criticality heuristics (Dinh et al., 2022).

Whole-space models typically conserve mass and a fourth-order energy. For the equation $L^2$3,

$L^2$4

and

$L^2$5

In the mixed-dispersion model the energy acquires the additional quadratic term $L^2$6. For the growing-weight equation $L^2$7, the critical index becomes

$L^2$8

with mass-critical exponent $L^2$9 and energy-critical exponent $H^s$0 for $H^s$1 (Alqaied et al., 10 Nov 2025).

2. Cauchy theory in $H^s$2, $H^s$3, and $H^s$4

The early whole-space theory established local and global well-posedness for the decaying-weight equation in $H^s$5 and $H^s$6, together with a small-data intercritical theory and a long-time perturbation principle. In the mass-subcritical regime, local $H^s$7 well-posedness extends globally by mass conservation. In $H^s$8, global continuation was obtained in the mass-subcritical regime and in the mass-critical regime under smallness of the $H^s$9-mass, while sufficiently small intercritical data yield global solutions and scattering. A later refinement enlarged the intercritical parameter range, especially in dimensions $\dot H^s$0, and also developed a $\dot H^s$1-critical local/global theory and stability result for $\dot H^s$2 under the restriction $\dot H^s$3 (Guzmán et al., 2019, Guzmán et al., 2021).

Fractional Sobolev local theory was then extended well beyond $\dot H^s$4 and $\dot H^s$5. For

$\dot H^s$6

and $\dot H^s$7, the Cauchy problem is locally well-posed in $\dot H^s$8, with maximal lifespan, uniqueness in the natural Strichartz class, and the standard blow-up alternative. This extended earlier work of Guzmán–Pastor and Liu–Zhang by enlarging the admissible ranges of $\dot H^s$9 and $b>0$0. Standard continuous dependence in the full $b>0$1 topology, rather than only in $b>0$2, was subsequently proved on compact subintervals of the lifespan, and the flow map was shown to be locally Lipschitz under additional regularity assumptions on the nonlinearity (An et al., 2022, An et al., 2023).

The critical theory was developed in Sobolev-Lorentz and Lorentz-type settings. For

$b>0$3

local well-posedness in $b>0$4 was obtained not only in the subcritical case $b>0$5 but also in the critical case $b>0$6, together with small-data global well-posedness and scattering when $b>0$7. A further refinement treated the $b>0$8-critical equation

$b>0$9

for $$i\partial_t u-\Delta^2u+\mu\Delta u+|x|^{-b}|u|^\alpha u=0$$0, $$i\partial_t u-\Delta^2u+\mu\Delta u+|x|^{-b}|u|^\alpha u=0$$1, and $$i\partial_t u-\Delta^2u+\mu\Delta u+|x|^{-b}|u|^\alpha u=0$$2, obtaining local well-posedness and small-data scattering in both $$i\partial_t u-\Delta^2u+\mu\Delta u+|x|^{-b}|u|^\alpha u=0$$3 and $$i\partial_t u-\Delta^2u+\mu\Delta u+|x|^{-b}|u|^\alpha u=0$$4 under less restrictive regularity assumptions on the nonlinear map $$i\partial_t u-\Delta^2u+\mu\Delta u+|x|^{-b}|u|^\alpha u=0$$5 than in earlier critical results (An et al., 2022, Jang et al., 2024).

3. Variational thresholds, scattering, and blow-up

In the intercritical focusing mixed-dispersion problem, a threshold theory is organized around the action

$$i\partial_t u-\Delta^2u+\mu\Delta u+|x|^{-b}|u|^\alpha u=0$$6

the Pohozaev functional

$$i\partial_t u-\Delta^2u+\mu\Delta u+|x|^{-b}|u|^\alpha u=0$$7

and the constrained minimization level

$$i\partial_t u-\Delta^2u+\mu\Delta u+|x|^{-b}|u|^\alpha u=0$$8

The invariant set

$$i\partial_t u-\Delta^2u+\mu\Delta u+|x|^{-b}|u|^\alpha u=0$$9

controls the global dynamics. If $|x|^b$0, then the corresponding solution is global and uniformly bounded in $|x|^b$1. The same framework yields spacetime estimates for $|x|^b$2, with exponent $|x|^b$3 in general and exponent $|x|^b$4 for radial data. Scattering in $|x|^b$5 then follows for all $|x|^b$6, for $|x|^b$7 with $|x|^b$8, and in the remaining low-dimensional cases under the stated radial hypotheses. When $|x|^b$9, this extends the scattering results of Saanouni and of Campos–Guzmán to dimensions three and four; when $d$0, the energy scattering result is new. The same paper also identifies $d$1 through a sharp Gagliardo–Nirenberg inequality and the ground state $d$2 solving $d$3 (Dinh et al., 2022).

At the energy-critical level, the non-radial focusing equation

$d$4

admits a sharp threshold formulation in terms of the ground state $d$5 solving

$d$6

For $d$7, if

$d$8

then the solution exists globally and scatters in both time directions in $d$9. The proof establishes energy trapping, a Palais–Smale compactness statement, and a concentration-compactness/rigidity argument in the Kenig–Merle style, adapted to the lack of translation invariance and to the fact that the kinetic energy is not itself conserved in the weighted setting (Guzmán et al., 4 Aug 2025).

Complementary blow-up results identify the opposite side of the ground-state barrier. For the focusing energy-critical equation, finite-time blow-up of the $u:\mathbb R\times\mathbb R^d\to\mathbb C$00-norm is proved for radial data when either $u:\mathbb R\times\mathbb R^d\to\mathbb C$01 or

$u:\mathbb R\times\mathbb R^d\to\mathbb C$02

For non-radial data, the conclusion depends on the size of $u:\mathbb R\times\mathbb R^d\to\mathbb C$03: when $u:\mathbb R\times\mathbb R^d\to\mathbb C$04, negative energy yields blow-up in the norm-divergence sense; when $u:\mathbb R\times\mathbb R^d\to\mathbb C$05, negative energy and the same positive-energy/large-kinetic condition imply finite-time blow-up. The mechanism is a localized virial argument with error terms estimated either by Strauss decay in the radial case or by the decay of the weight $u:\mathbb R\times\mathbb R^d\to\mathbb C$06 in the non-radial case (Scarpelli et al., 7 Jul 2025).

4. Analytical frameworks and proof technology

The earliest whole-space proofs relied on standard biharmonic Strichartz estimates, weighted Gagliardo–Nirenberg or Caffarelli–Kohn–Nirenberg inequalities, Hardy–Littlewood estimates, and explicit splitting of space into the unit ball and its complement to treat $u:\mathbb R\times\mathbb R^d\to\mathbb C$07. This framework supports local theory, small-data global theory, scattering criteria based on finiteness of critical Strichartz norms, and long-time perturbation results, but it generates dimension-dependent restrictions in low regularity or near-critical regimes (Guzmán et al., 2019, Guzmán et al., 2021).

A different route uses bilinear Strichartz-type estimates in Besov spaces. For coefficients $u:\mathbb R\times\mathbb R^d\to\mathbb C$08 behaving like $u:\mathbb R\times\mathbb R^d\to\mathbb C$09, new Besov-space bilinear estimates for the Duhamel operator were used to prove local well-posedness in the whole $u:\mathbb R\times\mathbb R^d\to\mathbb C$10-subcritical case with $u:\mathbb R\times\mathbb R^d\to\mathbb C$11. The key gain is that the product structure is handled through time-difference Besov norms rather than only pointwise Sobolev control, which is particularly effective when $u:\mathbb R\times\mathbb R^d\to\mathbb C$12 has limited differentiability and $u:\mathbb R\times\mathbb R^d\to\mathbb C$13 is singular (Liu et al., 2021).

The Lorentz-space program was developed to treat the weight more systematically. Sobolev-Lorentz spaces

$u:\mathbb R\times\mathbb R^d\to\mathbb C$14

support Hölder, Hardy, embedding, product, and chain-rule estimates that are well matched to the fact that $u:\mathbb R\times\mathbb R^d\to\mathbb C$15. A central advance was the extension of the Lorentz-space fractional chain rule from $u:\mathbb R\times\mathbb R^d\to\mathbb C$16 to all $u:\mathbb R\times\mathbb R^d\to\mathbb C$17, which made critical $u:\mathbb R\times\mathbb R^d\to\mathbb C$18-analysis possible in a unified framework (An et al., 2022).

In the intercritical scattering theory, Lorentz estimates also reshape the nonlinear decay mechanism. Dispersive and Strichartz estimates for $u:\mathbb R\times\mathbb R^d\to\mathbb C$19 were proved in Lorentz spaces, including derivative-gain estimates crucial in dimensions three and four. The weighted nonlinearity could then be estimated without splitting space into near and far regions, and localized virial/Morawetz functionals supplied the spacetime bounds needed for the scattering criterion (Dinh et al., 2022).

At the critical level, Lorentz-type refinements were pushed further to Besov-Lorentz and Triebel–Lizorkin–Lorentz spaces. Regular Strichartz estimates with gains in the spatial regularity index were established for the biharmonic propagator and then inserted into $u:\mathbb R\times\mathbb R^d\to\mathbb C$20-critical fixed-point arguments. In the non-radial energy-critical theory, this functional-analytic machinery is complemented by linear profile decomposition, nonlinear profile analysis, and rigidity via localized virial identities, producing a full concentration-compactness/rigidity scheme for the weighted fourth-order problem (Jang et al., 2024, Guzmán et al., 4 Aug 2025).

5. Boundary-driven and controlled fourth-order flows

The inhomogeneous biharmonic Schrödinger equation also appears as an initial-boundary value problem with prescribed boundary traces. On the half-line $u:\mathbb R\times\mathbb R^d\to\mathbb C$21, the problem with inhomogeneous Dirichlet–Neumann data

$u:\mathbb R\times\mathbb R^d\to\mathbb C$22

was analyzed by combining the Fokas unified transform with Fourier analysis. This yielded an explicit representation formula for the linear problem, fractional Sobolev space-time estimates, local well-posedness for $u:\mathbb R\times\mathbb R^d\to\mathbb C$23, low-regularity local well-posedness for $u:\mathbb R\times\mathbb R^d\to\mathbb C$24 through boundary Strichartz estimates, and global $u:\mathbb R\times\mathbb R^d\to\mathbb C$25 well-posedness for the defocusing model up to cubic nonlinearities (Özsarı et al., 2018).

On a bounded interval $u:\mathbb R\times\mathbb R^d\to\mathbb C$26, local well-posedness has been established for non-homogeneous Navier and Dirichlet boundary conditions. In the Navier case, the boundary traces are $u:\mathbb R\times\mathbb R^d\to\mathbb C$27 and $u:\mathbb R\times\mathbb R^d\to\mathbb C$28; in the Dirichlet case, they are $u:\mathbb R\times\mathbb R^d\to\mathbb C$29 and $u:\mathbb R\times\mathbb R^d\to\mathbb C$30. The admissible boundary data lie in the optimal trace spaces $u:\mathbb R\times\mathbb R^d\to\mathbb C$31, with compatibility conditions at $u:\mathbb R\times\mathbb R^d\to\mathbb C$32 when the Sobolev regularity is high enough. These results use boundary integral operators and a fixed-point scheme in $u:\mathbb R\times\mathbb R^d\to\mathbb C$33, with an additional $u:\mathbb R\times\mathbb R^d\to\mathbb C$34-based component at low regularity in the Navier problem (Li et al., 2020).

A low-regularity quarter-plane theory was developed for the cubic equation with inhomogeneous Dirichlet and Neumann data by introducing a Duhamel boundary forcing operator adapted to the fourth-order flow. The solution is constructed on the whole line and then corrected by boundary forcing terms whose coefficients solve an explicit linear system determined by boundary traces. Bourgain-type spaces $u:\mathbb R\times\mathbb R^d\to\mathbb C$35 and the associated mixed trace space $u:\mathbb R\times\mathbb R^d\to\mathbb C$36 permit local well-posedness for $u:\mathbb R\times\mathbb R^d\to\mathbb C$37, analytic dependence on the data, and a formulation that also suggests extensions to star graphs (Capistrano-Filho et al., 2018).

Control and stabilization theory on the torus $u:\mathbb R\times\mathbb R^d\to\mathbb C$38 concerns the fourth-order nonlinear Schrödinger equation

$u:\mathbb R\times\mathbb R^d\to\mathbb C$39

with internal forcing or damping supported in an arbitrary nonempty open subset. Using Bourgain spaces, observability for the linearized adjoint system, propagation of compactness and regularity, and unique continuation, one obtains local exact controllability near the origin, global exponential stabilization under localized damping, and global exact controllability on bounded $u:\mathbb R\times\mathbb R^d\to\mathbb C$40-sets (Filho et al., 2018).

6. Growing weights, geometric settings, and present limitations

A distinct branch of the theory replaces the decaying weight $u:\mathbb R\times\mathbb R^d\to\mathbb C$41 by a growing factor $u:\mathbb R\times\mathbb R^d\to\mathbb C$42. For

$u:\mathbb R\times\mathbb R^d\to\mathbb C$43

the main difficulty is that the source term is unbounded and translation invariance is lost in the opposite direction. The available results are radial: local well-posedness in the radial energy space $u:\mathbb R\times\mathbb R^d\to\mathbb C$44 for $u:\mathbb R\times\mathbb R^d\to\mathbb C$45 and in $u:\mathbb R\times\mathbb R^d\to\mathbb C$46 for $u:\mathbb R\times\mathbb R^d\to\mathbb C$47, together with small-data global well-posedness and scattering in $u:\mathbb R\times\mathbb R^d\to\mathbb C$48 in the energy-subcritical range. The decisive estimates are Strauss-type radial inequalities, which replace the Hardy-type tools used in decaying-weight problems; the same results also emphasize that radial symmetry is not merely cosmetic but structurally built into the argument (Alqaied et al., 10 Nov 2025).

Geometric variants arise on product manifolds. On the waveguide manifold $u:\mathbb R\times\mathbb R^d\to\mathbb C$49, normalized standing waves have been studied for

$u:\mathbb R\times\mathbb R^d\to\mathbb C$50

under the mass constraint $u:\mathbb R\times\mathbb R^d\to\mathbb C$51, in the mixed regime

$u:\mathbb R\times\mathbb R^d\to\mathbb C$52

The resulting variational picture depends on the sign of $u:\mathbb R\times\mathbb R^d\to\mathbb C$53: for $u:\mathbb R\times\mathbb R^d\to\mathbb C$54, mountain-pass normalized solutions are produced; for $u:\mathbb R\times\mathbb R^d\to\mathbb C$55, both local-minimum and mountain-pass branches occur for small masses. Orbital stability is then proved for the minimizer-type branches by the standard Cazenave–Lions strategy on the mass shell (Wang et al., 2024).

Several current limitations are explicit. In the mixed-dispersion scattering theory, the case $u:\mathbb R\times\mathbb R^d\to\mathbb C$56 is not covered because the linear dispersive and Strichartz theory is only local in time. In the growing-weight model, low-dimensional cases and further scattering problems are left open, and the radial assumption remains essential. In the earlier $u:\mathbb R\times\mathbb R^d\to\mathbb C$57-critical decaying-weight theory, the restrictions $u:\mathbb R\times\mathbb R^d\to\mathbb C$58 and $u:\mathbb R\times\mathbb R^d\to\mathbb C$59 are technical inputs of the contraction argument. A plausible implication is that future progress will continue to depend on replacing symmetry-dependent decay arguments and weight-splitting estimates by non-radial critical techniques of the kind already developed for the energy-critical whole-space equation (Dinh et al., 2022, Guzmán et al., 2021).