Papers

We develop an intrinsic theory of weighted, inhomogeneous Besov spaces on general homogeneous Lie groups without recourse to group-specific arguments. Starting from a definition in terms of localised test functions, we establish an equivalent, wavelet-like, multiscale characterisation. This provides a unified mechanism for deriving Besov embeddings, a Taylor-remainder characterisation, Young-type product estimates, Schauder estimates for convolution semigroups, and a weighted Kolmogorov criterion for random distributions. We apply this framework to parabolic Anderson-type equations associated with positive Rockland operators and with singular initial data. For space-time noise, our results cover the Young regime. For purely spatial noise, we also establish well-posedness in the first singular regime using a variant of the Cole-Hopf transform for Rockland operators. Both applications rely on a mild sewing lemma that accommodates time-dependent Banach spaces and increments with a singularity at the initial time.

For each nonnegative integer $m$, we construct smooth symmetric $3\times 3$ coefficient matrices $A_m$ satisfying the fixed ellipticity bound \[ I\leq A_m\leq 2^{81}I \] for which the smooth solutions of uniformly elliptic equations in nondivergence form \[ \text{tr}(A_m(x)D^2 u_m)=A_m(x):D^2u_m=0\qquad\text{in }B_2\subset {\mathbb R}^3 \] have common Dirichlet data, satisfy $\|u_m\|_{L^\infty(B_2)}\leq1$, but \[ \lim_{m\to \infty}\|Du_m\|_{L^1(B_1)}=\infty. \] Thus, there is no interior $W^{1,1}$ estimate depending only on ellipticity in dimension three, and consequently no such $W^{1,p}$ estimate for any $p\geq1$. This resolves in the negative an open question raised by Nadirashvili, Tkachev, and Vlăduţ. The construction also gives a uniformly convergent limit $u\notin \text{BV}_{\rm loc}(B_1)$ for a measurable uniformly elliptic coefficient matrix obtained as an $L^1$ limit of the $A_m$.

We study the Vlasov–Poisson–Boltzmann system on the three-dimensional torus under uniform shear flow. For Maxwell molecules with Grad's angular cutoff and sufficiently small shear rate, we consider perturbations around the spatially homogeneous self-similar profile of the sheared Boltzmann equation. In self-similar variables, we prove the global stability of this profile, including global existence and uniqueness, the exponential decay of nonzero-spatial derivatives in weighted $L^\infty$ spaces, and the uniform boundedness of the renormalized perturbation. The analysis combines a Caflisch's decomposition, Guo's $L^\infty$–$L^2$ estimates, macro–micro analysis, and a spectral study of the zero-frequency mode. We further derive a closed zero-frequency system for a renormalized total energy and suitable second-order moments, which yields a precise large-time description of the shear-induced energy growth.

We study a viscous one-velocity Baer–Nunziato system for two barotropic compressible fluids in a bounded three-dimensional domain. In contrast with models in which the volume fraction is merely transported or determined by an algebraic equilibrium constraint, it satisfies a differential closure driven by the pressure gap between the phases. For arbitrary finite-energy initial data and adiabatic exponents $\gamma^\pm>1$, we construct global-in-time dissipative measure-valued solutions and prove dissipative measure-valued–strong uniqueness relative to any sufficiently regular solution with the same initial data. The principal difficulties are the singular and non-continuous behavior of the pressure-relaxation source at vanishing volume fractions, the associated concentration defects, and the lack of direct coercivity of the standard thermodynamic relative energy with respect to the volume fraction. These are resolved by an endpoint cutoff argument and an augmented relative energy coupled to the renormalized volume-fraction equation.

The existence of global classical solutions for reaction-diffusion systems arising from chemical reaction networks remains a major open problem in the deterministic setting, especially for reactions with high polynomial growth. We prove that a suitably chosen, physically motivated transport noise yields unique strong solutions for complex balanced chemical reaction networks that are global in time with arbitrarily high probability. These solutions possess paths in $C^{\theta}_t C^{\infty}_x$ for all $\theta<1/2$ and are, in particular, classical in space. Furthermore, we show that a suitable transport noise can enhance the dissipation of spatial fluctuations at an arbitrarily prescribed exponential rate. Our proofs rely on a combination of scaling-limit arguments, maximal $L^p(L^q)$-regularity, and entropy-entropy dissipation estimates.

We establish the existence of an energy minimizer for a variational model of compressible magnetoelastic solids. The analysis is carried out in a new admissible class of deformations consisting of mappings of finite distortion, which extends previously available existence frameworks. A key ingredient is a compactness result under the critical integrability assumption on the outer distortion coefficient, which significantly weakens the regularity requirements imposed in earlier works. To obtain this result, we prove a diameter estimate for finite-distortion mappings satisfying the Ciarlet-Nečas condition and derive an open mapping theorem under the optimal integrability assumption, that is, the outer distortion is in $L^{n-1}$. This provides a partial positive result in the direction of the Iwaniec-Šverák conjecture and implies that admissible deformations are homeomorphisms. These topological and compactness properties allow us to apply the direct method of the calculus of variations and establish the existence of minimizers for compressible magnetoelastic solids within the admissible class of deformations consisting of mappings of finite distortions.

In this paper, we consider the Laplace operator on the half-space with Dirichlet and Neumann boundary conditions. These operators are studied on Sobolev spaces with power weights measuring the distance to the boundary. We prove optimal estimates for the resolvent operators and the corresponding heat semigroups. In addition, it is proved that the Dirichlet and Neumann Laplacians admit a bounded $H^\infty$-functional calculus on Sobolev spaces with certain compatibility conditions at the boundary. We show that these compatibility conditions cannot be omitted in general. The results in this paper are a direct extension of those obtained by Lindemulder, Lorist, the author, and Veraar in [J. Funct. Anal., 289(8):110985, 2025].

We are concerned with the global classical solution problem of 3D compressible isentropic Euler equations of Chaplygin gases \[ \begin{cases} \partial_t\rho + \mathrm{div}(\rho v) = 0,\\ \partial_t(\rho v) + \mathrm{div}(\rho v \otimes v) + \nabla p = 0,\\ \rho(0,x) = \bar\rho + \varepsilon\rho_0(x),\ v(0,x) = \varepsilon v_0(x). \end{cases} \] where $\bar\rho>0$ is a constant, $\varepsilon>0$ is small, the state equation is $p=p(\rho)=P_0-\frac{D}{\rho}$ with $P_0$ and $D$ being some positive constants. For the 3D compressible Euler equations of Chaplygin gases, which are a prototype of multidimensional nonlinear symmetric hyperbolic systems with totally linearly degenerate eigenvalues, there is a basic conjecture imposed by A. Majda: it typically has a global classical solution $(\rho, v)$ with $(\rho-\bar\rho, v)\in C([0,\infty), H^s(\Bbb R^3)) \cap C^1([0,\infty), H^{s-1}(\Bbb R^3))$ when $(\rho_0, v_0)\in H^s(\Bbb R^3)$ with $s>\frac52$ unless $(\rho, v)$ itself blows up in finite time. In this paper, under the assumptions that for any fixed constant $\mu$ with $0<\mu<1/2$, integer $N\geq 15$, $\mathrm{rot}\,v_0(x) \equiv 0$ and \[ \|(\rho_0, v_0)\|_{H^{N}(\mathbb{R}^3)}+\sum_{|a|\leq 13} \|\langle x\rangle^{1+\mu} \nabla^a(\rho_0,v_0)\|_{L^2(\mathbb{R}^3)} \leq 1, \] we show that the classical solution $(\rho, v)$ exists globally. Our main ingredients include: establishing a series of new decay estimates of energy bounds, weighted pointwise space-time $L^\infty$-$L^2$ estimates and weighted Strichartz-type estimates for the 3D potential flow equation of Chaplygin gases.

In this paper, we establish two improved Liouville-type theorems for the three-dimensional stationary tropical climate model without imposing any temperature assumptions. Owing to the special structure inherent to the equations, we fully exploit the $L^2$-estimate for the oscillation of the temperature, refined interpolation techniques and an iteration argument to establish a Liouville-type theorem under more flexible and relaxed power-law growth conditions. As an interesting consequence, we obtain the triviality of solutions under the integrability conditions imposed only on $(u,\nabla v)$ or $(\nabla u,\nabla v)$. Furthermore, by developing a systematic framework to handle the energy function associated with the non-trivial solutions and with the aid of delicate ODE analysis and a contradiction argument, we refine our Liouville-type theorem by introducing logarithmic corrections to the growth rates.

For the 2-D incompressible Euler equation, we prove the existence of $m$-fold symmetric uniformly rotating vortex patches with 90-degree corners for all $m\geq 12$. Properties of these patches and the induced stationary flows in the co-rotating frame are characterized. Our construction is based on a novel fixed-point method.

We prove uniqueness in the Lorentzian Calderón problem in a semiglobal setting, comparing a potentially large perturbation of the Minkowski metric against a small one. Our proof uses a reduced amount of data, solving a formally determined version of the problem. In contrast to traditional approaches, it employs neither control-theoretic arguments nor a reduction to a geometric inverse problem via microlocal methods or high-frequency solutions. Instead, it relies on distorted plane waves and weighted $L^2$-estimates.

This paper derives the second-order positive Burgers' equation from the standard Burgers' hierarchy to explore its complete integrability and exact analytical solutions. We construct this higher-order nonlinear evolution equation by systematically applying the recursion operator to the classical Burgers' equation. Expanding on this structural framework, we then derive the explicit third-order equation and use Complete Bell Polynomials to generalize the $n$-th order hierarchy. Through the Cole-Hopf transformation, the second-order nonlinear equation rigorously maps to the linear third-order dispersion equation. We establish complete integrability by explicitly formulating the scalar Lax pair (zero-curvature representation), which allows us to directly derive the associated differential and algebraic Darboux transformations. To systematically classify explicit stationary, time-dependent traveling wave, and self-similar profiles, we apply several techniques to the linear domain: separation of variables, the Hirota perturbation method, traveling wave reduction via the Complete Discrimination System for Polynomial Method (CDSPM), and Lie similarity reduction. Crucially, our analysis demonstrates the exact finite truncation of the Hirota perturbation series. We conclude by outlining how this equation impacts the theoretical understanding of fluid dynamics, nonlinear transport phenomena, and higher-order wave propagation.

We audit the exact-arithmetic certificate in arXiv:2310.19781v2, its 2025 Communications in Mathematical Physics version, and their shared MATLAB archive. CMP corrects none of the audited items; we therefore treat both versions together and use the openly accessible arXiv source for exact provenance. We found 11 proof-affecting defects (alongside typographical slips) and give reproducible exact counterexamples to several implemented bounds and reconstruction steps. To the best of our knowledge, this is the first independent, version-pinned, source-level audit of a published computer-assisted proof (CAP) in analysis to identify multiple proof-affecting defects in its released computational certificate. The findings do not refute the intended theorem or its analytic reduction, but they show that the published certificate does not prove the claimed conclusion. Future work will give a corrected, structurally simpler proof with substantially less machine assistance.

We study the multiscale behavior of a coupled visco-poroelastic system arising in the modelling of fibre-reinforced hydrogels (FIHs) used in tissue engineering scaffolds. The composite material consists of a periodic fibre scaffold, which is governed by quasi-static linear elasticity, and a hydrogel phase saturating the interstitial space, which is modelled as a Biot linear poroelastic medium enhanced with Kelvin–Voigt structural damping. The two phases are coupled through continuity of displacement and traction across their shared interface. Both the fibre scaffold and the hydrogel phase are connected, so that mechanical forces can be transmitted through the composite and interstitial fluid can flow directly through the hydrogel network. Starting from a microscopic ($\varepsilon$-scale) model, we derive uniform a-priori estimates and establish well-posedness via a Rothe time-discretisation argument for both the case of standard Biot fluid content $\eta=0$ and the case of viscous fluid content $\eta=\alpha\delta>0$. We then perform a rigorous two-scale homogenization in the limit $\varepsilon \to 0$ using periodic unfolding. In addition to the usual effective elasticity, storage, coupling, and permeability coefficients, the homogenized constitutive laws contain nonlocal-in-time memory terms generated by the microscopic viscoelastic relaxation. All effective coefficients and memory kernels are explicitly characterized in terms of the microscale geometry and material parameters.

In the present paper, we consider the Cauchy problem of a system of three nonlinear Schrödinger equations with quadratic nonlinearity. We first prove the existence of ground states in the form of solitary wave solutions with two frequency parameters. We then show that the conditions on the frequency parameters for the existence of ground states depend on the resonance structure of the system. Next, we give the two results for global solutions. The first is an improvement of global well-posedness for initial data below the ground state threshold. The second is an improvement of global well-posedness for oscillating initial data. We also prove the orbital stability of the ground state sets with small speed parameter.

We consider a parabolic relaxation formulation of the compressible Navier-Stokes-Korteweg system and prove the global-in-time existence of finite energy weak solutions to the associated initial-boundary-value-problem. Our proof is based on a three-level approximation scheme, a weak compactness property of the effective viscous flux, and parabolic regularity estimates. Our result holds for a broad variety of non-monotone pressure functions and generalizes the corresponding results known for the compressible Navier-Stokes equations to the relaxation system.

We establish exact anisotropic $L^p$-Hardy and Caffarelli-Kohn-Nirenberg identities associated with the Minkowski functional $\|\cdot\|_K$ and the anisotropic radial derivative $\mathcal R_K$, where $1<p<\infty$ and $K\subset\RN$ is a smooth and strictly convex body containing the origin in its interior, not necessarily origin-symmetric. These identities contain explicit nonnegative $R_p$-remainders and yield sharp anisotropic radial $L^p$-Hardy and $L^p$-Caffarelli-Kohn-Nirenberg inequalities. In the noncritical parameter regions, we characterize the full family of extremal functions in the natural weighted completion space, while in the critical case the sharp constant is not attained by any nonzero function in that space. As applications, we establish sharp anisotropic $L^p$-Heisenberg uncertainty principles and obtain explicit anisotropic Gaussian-type optimizers. Finally, letting $K^*$ denote the polar body of $K$, we derive norm-based anisotropic gradient inequalities associated with the origin-symmetric body $K^*\cap(-K^*)$. When $K$ is origin-symmetric, the corresponding gradient constants are sharp, including those in the critical case.

We study transparent boundary conditions (TBCs) for the time-dependent heat equation in a branching, quasi-one-dimensional domain modeled as a metric star graph. By combining the classical concept of TBCs with the theory of partial differential equations (PDEs) on networks, we derive an exact vertex condition that ensures unobstructed thermal flow across junctions. Specifically,we derive a sum rule for the diffusion coefficients that eliminates thermal backflow at the vertex.We demonstrate the validity of this analytical model numerically using a Crank-Nicolson finite-difference method, which confirms smooth, unobstructed heat propagation through the network.These results provide a practical mathematical framework for tunable control and optimization of thermal diffusion in low-dimensional structures. In particular, combining the well-known concept of the TBCs and theory for PDE on metric graphs, we propose a mathematical model providing a control tool for thermal diffusion in networks.

We generalize the classical vector-valued tangential and normal trace theory on Lipschitz domains to the setting of non-Lipschitz $H^1$-extension domains, which includes domains with fractal boundaries such as the Koch snowflake. We define and study these operators based on the surjectivity of the trace operator for elements of $H^1(\Omega)$ and the Hilbert structure of the associated trace space. The normal trace operator is defined on $\boldsymbol{H}({\rm div},\Omega)$ in $\mathbb{R}^n$. A generalized Stokes formula allows us to introduce the tangential trace operator on $\boldsymbol{H}({\bf curl},\Omega)$ and $\boldsymbol{H}^1(\Omega)$ in two and three dimensions. Following the approach of Buffa, Costabel and Sheen (2002) for Lipschitz domains, we define two abstract tangential boundary spaces as images of these trace operators, establish their Hilbert structure, and construct an abstract rotation operator linking them, in place of the geometric rotation based on the normal vector. We also extend Costabel's (1991) coercive bilinear form approach for Maxwell's equations to non-Lipschitz domains, yielding a theory that supports the treatment of the Hodge-Dirac operator and Green's formulas and enables the solution of boundary-value problems for the $\operatorname{\textbf{curl}}\operatorname{\textbf{curl}}+\textbf{I}$ operator within the $H^1$-extension domains framework.

Given $s\in(0,1)$ and $n>2s$, we construct nontrivial solutions of the problem \begin{equation*} \begin{cases} (-\Delta)^s u=u^{2^*_s-1} &{\mbox{ in }}\mathbb{R}^n_+,\\ {\mathcal{N}}_s u=0&{\mbox{ in }}\mathbb{R}^n_-,\end{cases} \end{equation*}where ${\mathcal{N}}_s$ represents the nonlocal Neumann condition of exterior type.

For $N\geq 1$, $p\geq 1$, and $\delta>0$, consider the nonlocal threshold functional \[ I_{\delta,p}(u)=\iint_{\{|u(x)-u(y)|>\delta\}} \frac{\delta^p}{|x-y|^{N+p}}\,\dd x\,\dd y. \] Nguyen and Squassina \cite{NguyenSquassina} asked whether $I_{\delta,p}$ decreases under Schwarz rearrangement. We give a negative answer in every dimension. For each $N\geq1$ and $\delta>0$, one and the same explicit counterexample works for every $p\geq1$: it is nonnegative, bounded, compactly supported, and takes only five values. Its four nontrivial superlevel sets are nested balls whose centers alternate between two points. Once the active threshold interactions are isolated, the energy difference reduces to comparing the interaction of a unit ball with a concentric annulus and with an eccentric shell. The sign is strict because the potential generated by the unit ball decreases with the radius. We also compute the energy gap in closed form in dimension one. This settles Open Problem~2.2 of Nguyen and Squassina.

In this paper, we prove the local null controllability of the Boussinesq system in dimensions two and three with a reduced number of localized controls. More precisely, the controls act on the temperature equation and on only $N-2$ components of the velocity equation. Our proof is based on a spectral approach in the spirit of the Lebeau–Robbiano method. The main difficulty comes from the coupled structure of the system. The velocity and temperature equations are governed by the Stokes operator and the Dirichlet Laplacian, respectively. Their natural spectral decompositions are different, and there is no common spectral localization naturally adapted to the coupled system. Therefore, the usual componentwise spectral argument cannot be applied directly. We show that the cascade structure of the linearized system makes it possible to overcome this difficulty. The main ingredient is a mixed observability estimate in which only the Stokes component is spectrally localized, while the heat component is treated without any frequency restriction. Combining this estimate with the dissipation of the high Stokes frequencies allows us to carry out a Lebeau–Robbiano iteration for the linearized system. The resulting linear controllability estimate is transferred to the nonlinear Boussinesq system through a time-iteration argument. As an important consequence, the controls have a small-time cost bounded by $C\exp(C/T)$, which recovers the expected parabolic order while preserving the reduced number of controls.

This paper studies the exterior Dirichlet problem for Hessian quotient equations with nonconstant right-hand sides. We prove the existence of classical admissible solutions with prescribed asymptotic Hessians and establish convergence of the Hessian at infinity. The main difficulties are obtaining second-order estimates on expanding annuli that are uniform in the outer radius and deriving Hessian convergence under an integral tail condition with no prescribed decay rate. These are resolved through a radius-independent boundary-to-interior estimate and a blow-down argument. We also allow nonradial perturbations of the source. Under stronger pointwise assumptions, we obtain higher-order asymptotic expansions and solutions for every sufficiently large prescribed asymptotic constant.

We study positive solutions of the critical fractional Hartree equation with a non-removable isolated singularity at the origin. This equation is doubly nonlocal, involving both the fractional Laplacian and a Riesz convolution potential. We first prove that every positive singular solution is radially symmetric about the origin, by combining the Caffarelli–Silvestre extension with the method of moving spheres. We then establish the existence of Delaunay-type periodic singular solutions. After the Emden–Fowler transformation, the Hartree convolution survives as a genuinely nonlocal integral term, so that the resulting periodic equation cannot be reduced to an ordinary differential equation. We construct nonconstant periodic solutions for all sufficiently large periods by minimizing a Rayleigh-type quotient in a periodic fractional Sobolev space.

We investigate the long-time behavior of small solutions to the Dirac-Hartree equation in two spatial dimensions. This model describes the mean-field dynamics of relativistic fermions interacting through the three-dimensional Coulomb potential $|x|^{-1}$, which gives rise to long-range effects in the scattering dynamics. We prove global well-posedness and long-range scattering (modified scattering) for small initial data in weighted Sobolev spaces. In this setting, long-range scattering means that, unlike linear scattering, an additional logarithmic phase correction is required to describe the precise asymptotics of nonlinear solutions. Our approach relies on the space-time resonance method, combined with special null structures inherent in the equation. Compared to the three-dimensional case \cite{CKLY2022,cloos}, the novelty lies in overcoming the weaker time decay inherent to the two dimensional problem.

In this paper, we investigate the semilinear parabolic equation $\partial_t f-\Delta f=a(x,t)f^p$ on a complete Riemannian manifold with Ricci curvature bounded from below. By means of Nash-Moser iteration, we establish Li-Yau type gradient estimates for positive solutions to this equation, where the coefficient function $a$ can be either strictly sign-definite or sign-changing. As an application, we derive Liouville theorems for ancient and eternal solutions on manifolds with nonnegative Ricci curvature, generalizing a number of classical results. Our proof features the incorporation and tuning of two parameters in the auxiliary quantities to accommodate the variable coefficient $a$ and to extend the admissible range of $p$.

We study the simultaneous limits of homogenisation $(\varepsilon\to0)$ and dimension reduction $(h\to0)$ for thin heterogeneous nematic liquid crystal films in the Landau–de Gennes $Q$-tensor framework. The elastic energy density is governed by a tensor $\mathbf{A}(x'/\varepsilon,x_3/h)$ that is periodic in the in-plane fast variable and measurable in the normalised thickness variable. Surface anchoring on the top and bottom faces is modelled by a weak anchoring energy of strength $h^\gamma$, $\gamma\geq0$, in a general set-valued framework covering the principal classical anchoring geometries. The simultaneous limit reveals two principal features. First, the anchoring scaling yields a hierarchy of effective behaviours depending on $\gamma$: a hard constraint $\mathcal{Q}\in H^1(\omega;\mathscr{A})$ for $0\leq\gamma<1$, a finite surface density contribution for $\gamma=1$, and vanishing anchoring for $\gamma>1$. For $0\leq \gamma<1$ with uniaxial anchoring, the homogenised energy reduces on the constrained class to anisotropic Oseen–Frank and Ericksen energies. Second, the effective elastic response depends on the scale ratio $\rho=\lim h/\varepsilon\in[0,\infty]$: each regime has a distinct corrector structure and yields a two-dimensional LdG energy with homogenised tensor $\mathbf{A}^{\mathrm{hom}}_\rho$, characterised by a regime-dependent cell problem. To the best of our knowledge, this is the first rigorous treatment of the simultaneous homogenisation and dimension reduction limit for the Landau–de Gennes energy, and the first in which the anchoring strength enters as a scaling parameter, producing a hierarchy of qualitatively distinct effective models.

We find an explicit rank-one convex non-quasiconvex integrand in $\mathbb{R}^{2\times 4}$: to falsify the quasiconvexity inequality, we exhibit a map $\mathbb{T}^4\to \mathbb{R}^2$ with $12$ non-zero Fourier modes. In fact, this map is obtained from a scalar potential, so we also find a rank one convex integrand in $\mathbb{R}^{4\times 4}_\text{sym}$ which is not quasiconvex. These examples are obtained by transpositions and restrictions of Grabovsky's example of a rank-one convex, non-quasiconvex integrand in $\mathbb{R}^{8 \times 2}.$ We also modify Šverák's example to construct a rank-one convex non-quasiconvex integrand in $\mathbb{R}^{3\times 3}_\text{sym}$.

We develop a general principle for passing from the analysis of partial differential operators on the whole space to boundary value problems on the half-space. The reduction separates the boundary problem into restricted invertibility of the interior operator and a complementing problem formulated entirely on the trace space. For a broad class of operators admitting suitable factorizations, the restricted invertibility is governed by the whole-space theory of the factors, while the complementing problem determines both well-posedness and the admissible combinations of interior and boundary data. In this way, boundary value theory is reduced to whole-space operator theory together with an explicit analysis on the trace space. We realize this principle for time-periodic partial differential equations and systems of mixed order, developing the required Newton-polygon function spaces, interpolation theory, and $L^p$-trace theory. For mixed-order differential operators under general boundary conditions, well-posedness and the admissible data space are determined by an explicit complemented boundary matrix. This yields new $L^p$-well-posedness results for Cahn–Hilliard equations with dynamic boundary conditions and for Cahn–Hilliard–Gurtin systems. Most notably, we obtain the first maximal $L^p$-regularity result for a genuinely mixed-order interior system coupled to dynamic boundary conditions carrying a distinct mixed-order structure.

The goal of this paper is to investigate the traveling wave solutions with one-sided far field conditions for the one-dimensional hyperbolic Keller-Segel equations with quorum sensitivity. The traveling wave solutions are piece-wise smooth in the moving coordinate, and satisfy the entropy inequality at the discontinuity.