
Citation: Nádia C. M. Okuyama, Fernando Cezar dos Santos, Kleber Paiva Trugilo, Karen Brajão de Oliveira. Involvement of CXCL12 Pathway in HPV-related Diseases[J]. AIMS Medical Science, 2016, 3(4): 417-440. doi: 10.3934/medsci.2016.4.417
[1] | Nickolas Giardetti, Amy Shapiro, Stephen Windle, J. Douglas Wright . Metastability of solitary waves in diatomic FPUT lattices. Mathematics in Engineering, 2019, 1(3): 419-433. doi: 10.3934/mine.2019.3.419 |
[2] | Federico Cluni, Vittorio Gusella, Dimitri Mugnai, Edoardo Proietti Lippi, Patrizia Pucci . A mixed operator approach to peridynamics. Mathematics in Engineering, 2023, 5(5): 1-22. doi: 10.3934/mine.2023082 |
[3] | Panayotis G. Kevrekidis . Instabilities via negative Krein signature in a weakly non-Hamiltonian DNLS model. Mathematics in Engineering, 2019, 1(2): 378-390. doi: 10.3934/mine.2019.2.378 |
[4] | Jonathan A. D. Wattis . Asymptotic approximations to travelling waves in the diatomic Fermi-Pasta-Ulam lattice. Mathematics in Engineering, 2019, 1(2): 327-342. doi: 10.3934/mine.2019.2.327 |
[5] | Lorenzo Pistone, Sergio Chibbaro, Miguel D. Bustamante, Yuri V. Lvov, Miguel Onorato . Universal route to thermalization in weakly-nonlinear one-dimensional chains. Mathematics in Engineering, 2019, 1(4): 672-698. doi: 10.3934/mine.2019.4.672 |
[6] | Simone Paleari, Tiziano Penati . Hamiltonian lattice dynamics. Mathematics in Engineering, 2019, 1(4): 881-887. doi: 10.3934/mine.2019.4.881 |
[7] | Tiziano Penati, Veronica Danesi, Simone Paleari . Low dimensional completely resonant tori in Hamiltonian Lattices and a Theorem of Poincaré. Mathematics in Engineering, 2021, 3(4): 1-20. doi: 10.3934/mine.2021029 |
[8] | Carlo Danieli, Bertin Many Manda, Thudiyangal Mithun, Charalampos Skokos . Computational efficiency of numerical integration methods for the tangent dynamics of many-body Hamiltonian systems in one and two spatial dimensions. Mathematics in Engineering, 2019, 1(3): 447-488. doi: 10.3934/mine.2019.3.447 |
[9] | Andrea Kubin, Lorenzo Lamberti . Variational analysis in one and two dimensions of a frustrated spin system: chirality and magnetic anisotropy transitions. Mathematics in Engineering, 2023, 5(6): 1-37. doi: 10.3934/mine.2023094 |
[10] | Raimondo Penta, Ariel Ramírez-Torres, José Merodio, Reinaldo Rodríguez-Ramos . Effective governing equations for heterogenous porous media subject to inhomogeneous body forces. Mathematics in Engineering, 2021, 3(4): 1-17. doi: 10.3934/mine.2021033 |
Peridynamics is a modern branch of solid mechanics and materials science, which models the physical interactions in a material continuum not by partial differential equations but in terms of nonlocal integro-differential equations, see [40,42] for an introduction and overview. The simplest dynamical model for a spatially one-dimensional and infinitely extended medium is
∂2tu(t,y)=+∞∫−∞f(u(t,y+ξ)−u(t,y),ξ)dξ | (1.1) |
with time t>0, material space coordinate y∈R, bond variable ξ∈R, and scalar displacement field u. The elastic force function f is usually supposed to satisfy Newton's third law of motion via
f(r,ξ)=−f(−r,−ξ) | (1.2) |
so that (1.1) can equivalently be written as
∂2tu(t,y)=∞∫0f(u(t,y+ξ)−u(t,y),ξ)−f(u(t,y)−u(t,y−ξ),ξ)dξ. | (1.3) |
The symmetry condition (1.2) further ensures that (1.3) admits both a Lagrangian and Hamiltonian structure and can hence be regarded as a nonlocal wave equation without dissipation, see also §3. For energetic considerations it is also useful to write
f(r,ξ)=∂rΨ(r,ξ), |
where the micro-potential Ψ(r,ξ) quantifies the contribution to the potential energy coming from the elastic deformation of the bond ξ. The corresponding energy density is given by
e(t,y)=12(∂tu(t,y))2+∞∫0Ψ(u(t,y+ξ)−u(t,y),ξ)dξ |
and the total energy is conserved according to ddt∫+∞−∞e(t,y)dy=0.
Peridynamics and lattice models In this paper we focus on two particular settings. In the simplified continuous-coupling case we have
Ψ(r,ξ)=α(ξ)Φ(β(ξ)r) | (1.4) |
with smooth reference potential Φ:R→R and scaling coefficients provided by two sufficiently regular functions α,β:R+→R+. This is a typical choice in peridynamics and [41] proposes for instance
a(ξ)=χ[0,H](ξ)ξ,β(ξ)=χ[0,H](ξ)ξ−1,Φ(r)={c2r2+c3r3for r≤0, c2r2for r>0, | (1.5) |
where H is the horizon and χJ is shorthand for the indicator function of the interval J. The discrete-coupling case corresponds to
Ψ(r,ξ)=M∑m=1Φm(r)δξm(ξ), | (1.6) |
where
0<ξ1<...<ξM<∞ |
represent a finite number of active bonds and δξm abbreviates a Dirac distribution centered at ξm. The integrals in the force term can hence be replaced by sums and the wave equation (1.3) is actually a nonlocal lattice differential equation. For instance, for M=1 and ξ1=1 we recover a variant of the classical Fermi-Pasta-Ulam-Tsingou chain (FPUT) while M=2 with ξ1=1 and ξ2=2 describe an atomic chain with spring-like bonds between the nearest and the next-to-nearest neighbors.
Of course, both the discrete-coupling case and the continuous-coupling case are closely related and (1.3) can be viewed as a generalized lattice equation with a continuum of active bonds. Moreover, (1.3)+(1.6) with
Φm(x)=α(εm)Φ(β(εm)r),M=ε−2 | (1.7) |
is a discretized version of (1.3)+(1.4), in which the integrals have been approximated by Riemann sums with grid size ε>0 on the interval (0,ε−1).
Traveling waves and eigenvalue problem for the wave profile A traveling wave is a special solution to (1.3) that satisfies
u(t,y)=U(x),x=y−σt | (1.8) |
with wave speed σ and profile function U depending on x, the spatial variable in the comoving frame. Combining this ansatz with (1.3) we obtain in the discrete-coupling case the advance-delay-differential equation
σ2U′′(x)=M∑m=1Φ′m(U(x+ξm)−U(x))−Φ′m(U(x)−U(x−ξm)) | (1.9) |
and a similar formula with integrals over infinitely many shift terms can be derived in the continuous-coupling case. As explained below in greater detail, the existence of solution to (1.9) with M=1 has been established by several authors using different methods but very little is known about the uniqueness and dynamical stability of lattice waves. For the peridynamical analogue with continuous coupling we are only aware of the rigorous existence result in [36] although there have been some attempts in the engineering community to construct traveling waves numerically or approximately, see [4,41]. For our purposes it is more convenient to reformulate the traveling wave equation as an integral equation for the wave profile
W(x):=U′(x), |
which provides the velocity component in a traveling wave via ∂tu(t,y)=σW(x). To this end we write
f(U(x+ξ)−U(x),ξ)−f(U(x)−U(x−ξ),ξ)=∂xF(x,ξ) | (1.10) |
with
F(⋅,ξ):=Aξ∂rΨ(AξW,ξ), | (1.11) |
where the operator Aξ defined by
(AξW)(x):=x+ξ/2∫x−ξ/2W(˜x)d˜x=(χξ∗W)(x) | (1.12) |
is the x-convolution with the indicator function χξ of the interval [−ξ/2+ξ/2]. Thanks to (1.10) and (1.11) we infer from (1.3) and (1.8) that any traveling waves must satisfy the nonlinear eigenvalue problem
σ2W=∞∫0Aξ∂rΨ(AξWξ)dξ+η, | (1.13) |
where η is a constant of integration. The eigenvalue problem (1.13) is very useful for analytical investigations because the convolution operator (1.12) exhibits some nice invariance properties. Moreover, it also gives rise to an approximation scheme that can easily be implemented and shows — at least in practice — good convergence properties. Both issues will be discussed below in greater detail and many of the key arguments have already been exploited in the context of FPUT chains, see [6,15].
To simplify the presentation, we restrict our considerations in this paper to the special case
η=0 | (1.14) |
but emphasize that this condition can always be guaranteed by means of elementary transformations applied to the wave profile and the micro-potentials, see for instance [15]. For solitary waves one can alternatively relate η to waves that are homoclinic with respect to a non-vanishing asymptotic state limx→±∞W(x)≠0.
Variational setting for solitary and periodic waves In this paper we construct solitary waves for systems with convex and super-quadratic interaction potential but for completeness we also discuss the existence of periodic waves. Other types of traveling waves or potentials are also very important but necessitate a more sophisticated analysis and are beyond the scope of this paper. Examples are heteroclinic waves or phase transition waves with oscillatory tails in systems with double-well potential as studied in [4,16,21,23,39,46] for atomic chains.
For periodic waves we fix a length parameter 0<L<∞ and look for 2L-periodic wave profiles that are square-integrable on the periodicity cell. This reads
W∈L2(IL):={W:R→RwithW(⋅)=W(⋅+2L)and‖W‖2<∞} |
with
‖W‖22:=∫ILW(x)2dx,IL:=(−L,+L). |
For solitary waves we formally set L=∞ and show in §2.5 that the corresponding waves can in fact be regarded as the limit of periodic waves as L→∞. We also mention that any Aξ from (1.12) is a well-defined, bounded, and symmetric operator on L2(IL), see Lemma 2 below.
Both for periodic and solitary waves we consider the potential energy functional
P(W):=∞∫0∫ILΨ((AξW)(x),ξ)dxdξ |
as well as the functional
K(W):=12‖W‖22, |
where the latter quantifies after multiplication with σ2 the integrated kinetic energy of a traveling wave. Using these energies we can reformulate the traveling wave equation (1.13) with (1.14) as
σ2∂K(W)=∂P(W), | (1.15) |
where ∂ denotes the Gâteaux derivative with respect to W and the inner product in L2(IL). In particular, in the discrete-coupling case we have
P(W)=M∑m=1∫ILΦm(AξmW)dx,∂P(W)=M∑m=1AξmΦ′m(AξmW) | (1.16) |
while the continuous-coupling case corresponds to
P(W)=∞∫0∫ILα(ξ)Φ(β(ξ)AξW)dxdξ,∂P(W)=∞∫0α(ξ)β(ξ)AξΦ′(β(ξ)AξW)dξ, | (1.17) |
where we omitted the x-dependence of W and AξW to ease the notation.
Existence result The existence of periodic and solitary traveling waves in FPUT chains has been studied intensively over the last two decades and any of the proposed methods can also be applied to peridynamical media although the discussion of the technical details might be more involved.
1. The Mountain Pass Theorem allows to construct nontrivial critical points of the Lagrangian action functional
W↦σ2K(W)−P(W) |
provided that the potential energy is super-quadratic and that the prescribed value of the wave speed σ is sufficiently large. For details in the context of FPUT chains we refer to [34] and the references therein.
2. Another variational setting minimizes K subject to prescribed value of P, where σ−2 plays the role of a Lagrange multiplier. This approach was first described in [13] for FPUT chains with super-quadratic potentials and has recently been generalized to peridynamical media in [36].
3. It is also possible to construct traveling waves by maximizing P under the constraint of prescribed K. This idea was introduced in [6], has later been refined in [15] and provides also the base for our approach as it allows to impose additional shape constraints as discussed in §2. Moreover, a homogeneous constraint is more easily imposed than a non-homogeneous one and this simplifies the numerical computation of traveling waves.
4. The concepts of spatial dynamics and center-manifold reduction have been exploited in [29,30] for lattice waves. This non-variational approach is restricted to small amplitude waves but provides — at least in principle — a complete picture on all solutions to (1.13). Moreover, explicit or approximate solutions are available for some special potentials or asymptotic regimes. We refer to [43,47] and the more detailed discussion in §3.
Our key findings on the existence of peridynamical waves can informally be summarized as follows. The precise statements can be found in Assumption 1, Corollary 8, Proposition 9, and Proposition 10.
Main result (existence of solitary waves). Suppose that the coefficients and interaction potentials in (1.4) or (1.6) are sufficiently smooth and satisfy natural integrability and super-quadraticity assumptions. Then there exists a family of solitary waves which is parameterized by K=K(W) and comes with unimodal, even, and nonnegative profile functions W. Moreover, each solitary wave can be approximated by periodic ones.
Notice that almost nothing is known about the uniqueness or dynamical stability of solitary traveling waves in peridynamical media and it remains a very challenging task solve the underlying mathematical problems. In §3 we thus discuss the asymptotic regimes of near-sonic and high-speed waves, in which one might hope for first rigorous results in this context.
Improvement dynamics and numerical simulations A particular ingredient to our analysis is the improvement operator
T(W):=μ(W)∂P(W),μ(W):=‖W‖2/‖∂P(W)‖2, | (1.18) |
which conserves the norm constraint K but increases the variational objective function P, see Proposition 4. Moreover, the operator T respects — under natural convexity and normalization assumptions on the involved micro-potentials — the unimodality, the evenness, and the positivity (or negativity) of W. This enables us in §2 to restrict the constrained optimization problem to a certain shape cone in L2(IL) without changing the Euler-Lagrange equation for maximizers, see Corollary 5.
By iterating (1.18) in the discrete dynamics
W0⇝W1=T(W0)⇝...⇝Wn+1=T(Wn)⇝... | (1.19) |
we can also construct sequences (Wn)n∈N of profiles functions such that K(Wn) remains conserved while P(Wn) is bounded and strictly increasing in n. Using weak compactness for periodic waves as well as a variant of concentration compactness for solitary waves one can also show that there exist strongly convergent subsequences and that any accumulation point must be a traveling wave, see the proofs of Lemma 6 and Theorem 7. However, due to the lack of uniqueness results we are not able to prove the uniqueness of accumulation points or the convergence of the whole sequence. Moreover, the set of accumulation points might depend on the choice of the initial datum W0.
The improvement mapping is nonetheless very useful for numerical purposes. Fixing a periodicity length L it can easily be implemented:
1. divide the x-domain IL and the ξ-domain [0,∞) into a large but finite number of subintervals of length 0<ε≪1, and
2. replace all integrals with respect to x and ξ by Riemann sums as in (1.7).
The resulting scheme has been used to compute the numerical data presented in this paper. It shows very good and robust convergence properties in practice, and this indicates for a wide class of peridynamical media that there exists in fact a unique stable solution to the constrained optimization problem. A first example is shown in Figure 1 and corresponds to the peridynamical medium (1.5), for which approximate solitary waves are constructed in [41] by means of formal asymptotic expansions and ODE arguments. Notice that W0 was always chosen to be non-positive in order to pick up the super-quadratic branch in the potential and that there are no solitary waves extending into the harmonic branch, cf. the discussion at the beginning of §2.3.
We finally mention that variants of (1.19) have already been used in [5,6,15] for the numerical computations of periodic waves in FPUT chains. Moreover, the improvement dynamics shares some similarities with the Petviashvili iteration for traveling waves in Hamiltonian PDEs, see for instance [37]. Both approximation schemes combine linear pseudo-differential operators with pointwise nonlinearities and a dynamical normalization rule, but the details are rather different. In particular, (1.19) comes with a Lyapunov function but does not allow to prescribe the wave speed.
We now study the existence problem for a class of potentials in a certain cone of L2(IL) and begin by clarifying the precise setting.
A smooth function Φ:[0,∞)→[0,∞) is called weakly super-quadratic if it does not vanish identically and satisfies
Φ(0)=Φ′(0)=0,Φ′′(r)r≥Φ′(r)≥0for allr≥0. | (2.1) |
Such a function is increasing, convex, and grows at least quadratically since elementary arguments — including differentiation with respect to λ — show that
Φ′(r)r≥2Φ(r),Φ(λr)≥λ2Φ(r)forλ>1, | (2.2) |
Examples are the harmonic potential Φ(r)=cr2 and any analytic function with non-negative Taylor-coefficients. In what follows we rely on the following standing assumption.
Assumption 1. (admissible potentials) In the discrete-coupling case, each potential Φm is supposed to be super-quadratic in the sense of (2.1). In the continuous-coupling case, we assume that Φ is super-quadratic and that α, β are nonnegative functions on (0,∞) such that
∞∫0(ξ+ξ2)α(ξ)β2(ξ)ϕ(ξ1/2β(ξ)√2K)dξ<∞ |
holds for all K>0 with ϕ(r):=r−1|Φ′(r)|. Moreover, in the proofs we suppose that each potential is three times continuously differentiable on the interval [0,∞).
Both for finite and infinite L, we denote by C the positive cone of all L2(IL)-functions that are even, nonnegative and unimodal, where the latter means increasing and decreasing for −L<x<0 and 0>x>L, respectively. This reads
C:=¯{W∈C∞c(IL):W(x)=W(−x)≥0 and −W′(x)=W′(−x)≥0 for all x≥0}, |
where the closure has to be taken with respect to the norm in L2(IL), and we readily verify that C is a convex cone and closed under both strong and weak convergence. Moreover, the decay estimate
0≤W(x)≤‖W‖2√2|x| | (2.3) |
holds for any W∈C and any x∈IL.
We next collect some important properties of the convolution operators Aξ and recall that sinc(y):=y−1sin(y).
Lemma 2. (properties of Aξ) For given ξ∈(0,∞), the following statements are satisfied for both finite and infinite L:
1. Aξ maps L2(IL) to L2(IL)∩L∞(IL) with ‖AξW‖2≤ξ‖W‖2 and ‖AξW‖∞≤ξ1/2‖W‖2.
2. Aξ is symmetric.
3. Aξ diagonalizes in Fourier space and has symbol function ˆAξ(k)=ξsinc(kξ/2).
4. Aξ respects the evenness, the non-negativity, and the unimodality of functions.
Moreover, Aξ is compact for L<∞ as it maps weakly converging into strongly converging sequences.
Proof. All assertions follow from elementary computations and standard arguments, see for instance [15,Lemma 2.5].
We are now able to show that our assumptions and definitions allow for a consistent setting of the variational traveling wave problem inside the cone C.
Lemma 3. (properties of the potential energy functional) The functional P is well-defined, strongly continuous and Gâteaux differentiable on the cone C, where the derivative maps C continously into itself. Moreover, P is convex, super-quadratic in the sense of
〈∂P(W),W〉≥2P(W)≥0,P(λW)≥λ2P(W)for λ≥1, | (2.4) |
and non-degenerate as ∂P(W)≠0 holds for any W∈C with P(W)>0.
Proof. We present the proof for the continuous-coupling case (1.4) only but emphasize that the arguments for (1.6) are similar.
Integrability: For given W∈C with K(W)=K we estimate
0≤Φ′(β(ξ)AξW)≤β(ξ)(sup0≤r≤β(ξ)‖AξW‖∞ϕ(r))AξW≤β(ξ)ϕ(ξ1/2β(ξ)√2K)AξW |
with ϕ as in Assumption 1, where we used that the super-quadraticity of Φ implies the monotonicity of ϕ. The estimates from Lemma 2 imply via
‖AξΦ′(β(ξ)AξW)‖2≤ξ2β(ξ)ϕ(ξ1/2β(ξ)√2K)‖W‖2 |
the estimate
‖∂P(W)‖2≤∞∫0α(ξ)β(ξ)‖AξΦ′(β(ξ)AξW)‖2dξ<∞, |
and similarly we show P(W)<∞ using (2.2). Moreover, since C is invariant under Aξ, multiplication with non-negative constants, and the composition with Φ′, we easily show that ∂P(W)∈C.
Continuity: Let (Wn)n∈N⊂C be a sequence that converges strongly in L2(IL) to some limit W∞, where the closedness properties of C imply W∞∈C. Extracting a subsequence we can assume that Wn converges pointwise to W∞ and that there exists a dominating function ¯W∈L2(Il) such that |Wn(x)|≤¯W(x) for all n∈N and almost all x∈IL. Using the similar arguments as above we readily demonstrate that the functions
Vn(ξ,⋅):=α(ξ)β(ξ)AξΦ′(β(ξ)AξWn) |
converge pointwise in (ξ,x) to V∞ and are moreover dominated by a function ¯V∈L1((0,∞),L2(IL)). Lebesgue's theorem ensures the strong convergence ∂P(Wn)→∂P(W∞) in L2(IL) along the choosen subsequence and due to the uniqueness of the accumulation point we finally establish the convergence of the entire sequence by standard arguments. The proof of P(Wn)→P(W∞) is analogous.
Further properties: The convexity and super-quadraticity of P is a direct consequence of the corresponding properties of Φ, and the convexity inequality
P(˜W)−P(W)≥〈∂P(W),˜W−W〉 | (2.5) |
ensures that any critical point of P must be a global minimizer. In particular, we find the implication
P(W)>0=P(0)=minP⟹∂P(W)≠0 | (2.6) |
and the proof is complete.
We emphasize that all arguments presented below hold analogously with W∈−C for potentials that are super-quadratic on (−∞,0], see Figure 1 for an application. If Φ is super-quadratic on R, our results imply under certain conditions the existence of two family of waves, which are usually called expansive (W≥0) or compressive (W≤0).
In this paper, we prove the existence of traveling waves with unimodal profile functions for 0<L≤∞ by solving the constrained optimization problem
maximize P(W) subject to W∈CK | (2.7) |
with
CK:={W∈C:P(W)>0 and K(W)=K}. | (2.8) |
We further set
P(K):=sup{P(W):W∈CK} | (2.9) |
and recall that σ2, the square of wave speed, can be regarded as the Lagrange multiplier to the norm constraint K(W)=K. The second important observation, which we next infer from the properties of the improvement dynamics, is that the shape constraint W∈C does not contribute to the Euler-Lagrange equation for solutions to (2.7). The analogous observation for FPUT chains has been reported in [15] and a similar result has been derived in [44] using different techniques.
Proposition 4. (properties of the improvement dynamics) The set CK in (2.8) is invariant under the action of the improvement operator T in (1.18) and we have
P(T(W))≥P(W) |
for any W∈CK, where equality holds if and only if W=T(W).
Proof. The operator T is well defined on C due to (2.6) and the claimed invariance property is a direct consequence of Lemma 3. The convexity inequality (2.5) implies
P(T(W))−P(W)≥〈∂P(W),T(W)−W〉=μ−1(W)〈T(W),T(W)−W〉=12μ−1(W)‖T(W)−W‖22, | (2.10) |
where we used μ(W)>0 and that ‖W‖2=‖T(W)‖2 holds by construction.
Corollary 5. Any solution W∈CK to the optimization problem (2.7) satisfies the traveling wave equation (1.15) for some Lagrange multiplier σ2≥K−1P(K)>0.
Proof. Any maximizer W belongs to CK, so Proposition 1 ensures in combination with (1.18) the validity of (1.15) with σ2=μ−1(W). Testing this equation with W we finally obtain
σ22K=〈σ2W,W〉=〈∂P(W),W〉≥2P(W)=2P(K) |
thanks to (2.4).
It remains to show the existence of maximizers and here we have to distinguish between the periodic case (L<0) and the solitary one (L=∞). The latter is more involved and will be discussed in the subsequent section.
Lemma 6. (existence of periodic maximizers) For any K>0 and any 0<L<∞, the constrained optimization problem (2.7) admits at least one solution.
Proof. Since the discrete coupling-case is rather simple, we discuss the continuous-coupling case only.
Weak continuity of P_: Suppose that (Wn)n∈N⊂CK converges weakly to W∞ and notice that the properties of convolution operators imply the pointwise convergence Vn(x,ξ)→V∞(x,ξ) for any x∈IL and any ξ with
Vn(x,ξ):=α(ξ)Φ(β(ξ)(AξWn)(x)). |
By (2.2) and Lemma 2 we also have
|(AξWn)(x)|≤ξ1/2‖Wn‖2 |
as well as
0≤Vn(x,ξ)≤12α(ξ)ϕ(ξ1/2√2K)(β(ξ)(AξWn)(x))2≤ξα(ξ)β2(ξ)ϕ(ξ1/2√2K)K, |
and since the right hand side is an integrable majorant on IL×[0,∞) according to Assumption 1, we conclude that P(Wn)→P(W∞).
Direct Method: Now let (Wn)n∈N⊂CK be a maximizing sequence for (2.7). By weak compactness we can assume that Wn converges weakly to W∞ and find
W∞∈C,‖W∞‖2≤√2K,P(W∞)=P(K)>0. |
Since (2.4) implies
P(K)≥P(√2K‖W∞‖2W∞)≥2K‖W∞‖22P(W∞) |
we finally conclude that ‖W∞‖2=√2K, i.e., W∞ belongs to CK. The assertion now follows from Corollary 5.
The arguments in the proofs of Lemma 4 and Lemma 6 can also be used to show that any discrete orbit of the improvement dynamics (1.19) contains a strongly convergent subsequence and that any accumulation point must be a traveling wave. In particular, the estimate (2.10) can be viewed as a discrete analogue to LaSalle's invariance principle. It remains open under which conditions the optimization problem (2.7) has a unique solution and whether there exists further local maxima or unstable saddle points, see also the discussion in §1.
We further emphasize that Lemma 6 does not automatically imply that the maximizer has a non-constant profile function. For purely harmonic potentials one can show by means of Fourier transform — see the discussion in §3.1 — that the periodic maximizer is always constant and in numerical simulations with super-quadratic functions that grow rather weakly near the origin we observe that the maximizer is constant for small K but non-constant for large K. A typical example with discrete coupling is presented in Figure 2 and relies on
M=2,L=5,Φ1(r)=Φ2(r)=12r2+16r6, | (2.11) |
i.e., on potentials with positive second but vanishing third derivative at r=0. Similar energetic localization thresholds can be found in [15] for FPUT chains and in [48] for coherent structures in other Hamiltonian systems. Below — see the comments to Proposition 2 and Proposition 10 — we discuss sufficient conditions to guarantee that the maximizer is non-constant.
Another, more degenerate and less smooth example with localization threshold is presented in Figure 3 and concerns the chain
\begin{align} M = 1\, , \qquad L = 4\, , \qquad \Phi{\left({r}\right)} = \left\{\begin{array}{lcl} \tfrac12r^2&&\text{for $0\leq r\leq1$}\\ r-\tfrac12+\tfrac52{\left({r-1}\right)}^2&&\text{for $r \gt 1$} \end{array}\right. \end{align} | (2.12) |
with piecewise linear stress strain relation \Phi^\prime , for which explicit approximation formulas have been derived in [47].
We finally mention that there is no localization threshold in Hertzian chains with monomial interaction potential due to an extra scaling symmetry. We refer to Figure 4 for a numerical example with
\begin{align} M = 1\, , \qquad L = 4\, , \qquad \Phi{\left({r}\right)} = r^p\, , \end{align} | (2.13) |
and to [5,31,32] for a variant of the corresponding improvement dynamics, the super-exponential decay of the wave profiles, and asymptotic results for p\gtrapprox 1 .
Solitary traveling waves do not exist in the harmonic case in which the involved potentials are purely quadratic. This follows by Fourier analysis arguments as described in §3, which reveal that the linear variant of the operator \partial \mathcal{P} has no proper eigenfunctions but only continuous spectrum. In the variational setting, the non-existence of harmonic solitary waves becomes manifests in the lack of strong compactness for maximizing sequences. The usual strategy for both FPUT chains and peridynamical media is to require the potentials to grow super-quadratically in a proper sense and to exploit concentration compactness arguments to exclude the crucial vanishing and splitting scenarios for weakly convergent sequences. We follow a similar approach but the compactness conditions simplify since we work with unimodal functions.
In what follows we call the optimization problem (2.7) genuinely super-quadratic if
\begin{align} P{\left({K}\right)} \gt Q{\left({K}\right)}\, . \end{align} | (2.14) |
Here, P{\left({K}\right)} is defined in (2.9) and
\begin{align*} Q{\left({K}\right)}: = \sup\big\{ \mathcal{Q}{\left({W}\right)}\;:\, W\in \mathcal{C}_K\big\} \end{align*} |
with
\begin{align*} \mathcal{Q}{\left({W}\right)} = \sum\limits_{m = 1}^{M}\tfrac12\Phi^{\prime\prime}_m{\left({0}\right)}\int\limits_{I_L}{\big( \mathcal{A}_{\xi_m} W\big)}^2\, \mathrm{d}x\qquad \text{resp.}\qquad \mathcal{Q}{\left({W}\right)} = \tfrac12 \Phi^{\prime\prime}{\left({0}\right)}\int\limits_0^\infty {\alpha}(\xi)\beta^2(\xi)\int\limits_{I_L}{\big( \mathcal{A}_\xi W\big)}^2\, \mathrm{d}x{\rm d}\xi \end{align*} |
represents the harmonic analogue to (2.7). Notice that (2.4) combined with the 2 -homogeneity of \mathcal{Q} ensures that if (2.14) is satisfied for some K = K_* , then it also holds for all K\geq K_* .
We first show that (2.14) guarantees the existence of solitary waves and identify afterwards sufficient conditions for the nonlinear potential functions. Similar ideas have been used in [15] in the context of FPUT chains.
Theorem 7. (strong compactness of maximizing sequences) Let L = \infty and suppose that K satisfies (2.14). Then, any sequence {\left({W_n}\right)}_{n\in {\mathbb{N}}}\subset \mathcal{C}_K with \mathcal{P}{\left({W_n}\right)}\to P{\left({K}\right)} admits a strongly convergent subsequence.
Proof. Preliminaries: To elucidate the key ideas, we start with the discrete-coupling case (1.6) and discuss the necessary modification for (1.4) afterwards. Passing to a (not relabeled) subsequence we can assume that
\begin{align*} W_n\xrightarrow{n\to\infty}W_\infty\quad\text{weakly}\;{\rm in}\; {\mathsf{L}}^2({\mathbb{R}}) \end{align*} |
for some W_\infty\in \mathcal{C} with \left\|{W_\infty}\right\|_2^2\leq 2K . Our goal is to show
\begin{align} \left\|{W_\infty}\right\|_2^2\geq 2K \end{align} | (2.15) |
because this implies in combination with the weak convergence the desired strong convergence in the Hilbert space {\mathsf{L}}^2({\mathbb{R}}) .
Truncation in x : For given X > \xi_M we write
\begin{align*} W_n = \widetilde{W}_n+{\overline{W}}_n\qquad\text{with}\qquad \widetilde{W}_n: = \chi_{[{-X}, {+X}]}W_n\, , \quad {\overline{W}}_n: = \chi_{ {\mathbb{R}}\setminus[{-X}, {+X}]}W_n \end{align*} |
and observe that
\begin{align} \big\|{\widetilde{W}_n}\big\|_2^2+\big\|{{\overline{W}}_n}\big\|_2^2 = 2K\, . \end{align} | (2.16) |
Moreover, by (1.16) we have
\begin{align*} \mathcal{P}{\left({W_n}\right)} = \sum\limits_{m = 1}^M \mathcal{P}_m{\left({W_n}\right)} \end{align*} |
with
\begin{align*} \mathcal{P}_m{\left({W_n}\right)}&: = \int\limits_{x\in {\mathbb{R}}}\!\!\! \Phi_m{\left({ \mathcal{A}_{\xi_m}W_n}\right)}\, \mathrm{d}x \\& = \int\limits_{\left|{x}\right|\leq X-\xi_m/2}\!\!\! \Phi_m{\left({ \mathcal{A}_{\xi_m}W_n}\right)}\, \mathrm{d}x \;\;+\!\!\!\!\!\int\limits_{\left|{\left|{x}\right|- X}\right|\leq \xi_m/2} \!\!\!\Phi_m{\left({ \mathcal{A}_{\xi_m}W_n}\right)}\, \mathrm{d}x\;\;+\!\!\!\!\! \int\limits_{\left|{x}\right|\geq X+\xi_m/2} \!\!\!\Phi_m{\left({ \mathcal{A}_{\xi_m}W_n}\right)}\, \mathrm{d}x \, . \end{align*} |
The unimodality and the evenness of \mathcal{A}_{\xi_m}W_n combined with Lemma 2 and (2.3) imply
\begin{align*} 0\leq \int\limits_{\left|{\left|{x}\right|- X}\right|\leq \xi_m/2} \Phi_m{\left({ \mathcal{A}_{\xi_m}W_n}\right)}\, \mathrm{d}x\leq 2\xi_m \Phi_m{\left({\frac{\xi_m\sqrt{2K}}{\sqrt{2{\left({X-\xi_m/2}\right)}}}}\right)}\leq o_X{\left({1}\right)}\, , \end{align*} |
where o_X{\left({1}\right)} denotes a quantity that does not depend on n or m and becomes arbitrarily small as X\to\infty . In the same way we derive
\begin{align*} 0\leq \int\limits_{\left|{x}\right|\leq X-\xi_m/2} \Phi_m{\big( \mathcal{A}_{\xi_m}W\big)}\, \mathrm{d}x\leq \int\limits_{ {\mathbb{R}}} \Phi_m{\big( \mathcal{A}_{\xi_m}\widetilde{W}_n\big)}\, \mathrm{d}x+o_X{\left({1}\right)} = \mathcal{P}_m{\big(\widetilde{W}_n\big)} +o_X{\left({1}\right)} \end{align*} |
as well as
\begin{align*} 0\leq \int\limits_{\left|{x}\right|\geq X+\xi_m/2} \Phi_m{\big( \mathcal{A}_{\xi_m}W\big)}\, \mathrm{d}x\leq \mathcal{P}_m{\big({\overline{W}}_n\big)} +o_X{\left({1}\right)}\leq \mathcal{Q}_m{\big({\overline{W}}_n\big)} +o_X{\left({1}\right)} \end{align*} |
due to the smoothness of \Phi_m and since the convolution kernel corresponding to \mathcal{A}_{\xi_m} is supported in [{-\xi_m/2}, {+\xi_m/2}] . In summary, we have shown that
\begin{align} \mathcal{P}{\big(W_n\big)}\leq \mathcal{P}{\big(\widetilde{W}_n\big)}+\mathcal{Q}{\big({\overline{W}}_n\big)}+o_X{\left({1}\right)} \end{align} | (2.17) |
holds for any given X > \xi_M and all n\in {\mathbb{N}} .
Further estimates: By construction, \widetilde{W}_n converges weakly as n\to\infty to \widetilde{W}_\infty and \mathcal{A}_{\xi_m}\widetilde{W}_n converges (for every m ) pointwise to \mathcal{A}_{\xi_m}\widetilde{W}_\infty . Moreover, \mathcal{A}_{\xi_m}\widetilde{W}_n is pointwise bounded by \sqrt{\xi_m}\sqrt{2K} and compactly supported in [{-X-\xi_M/2}, {+X+\xi_M/2}] . We thus conclude that A_{\xi_m}\widetilde{W}_n converges strongly in {\mathsf{L}}^2{\left({ {\mathbb{R}}}\right)} to \mathcal{A}_{\xi_m}\widetilde{W}_\infty and this implies
\begin{align*} \mathcal{P}{\big(\widetilde{W}_n\big)}\leq \mathcal{P}{\big(\widetilde{W}_\infty\big)}+o_n{\left({1}\right)}\, , \end{align*} |
where o_n{\left({1}\right)} is allowed to depend on X but vanishes in the limit n\to\infty . On the other hand, the super-quadraticity of \mathcal{P} , see (2.4), implies
\begin{align*} \mathcal{P}{\big(\widetilde{W}_\infty\big)}\leq\frac{\big\|{\widetilde{W}_\infty}\big\|_2^2}{2K} \mathcal{P}{\left({\frac{\sqrt{2K}}{\big\|{\widetilde{W}_\infty}\big\|_2}\widetilde{W}_\infty}\right)}\leq \frac{\big\|{\widetilde{W}_\infty}\big\|_2^2}{2K} P{\left({K}\right)} \end{align*} |
and in the same way we demonstrate
\begin{align*} \mathcal{Q}{\big({\overline{W}}_n\big)}\leq \frac{\big\|{{\overline{W}}_n}\big\|_2^2}{2K}Q{\left({K}\right)}\, . \end{align*} |
We thus deduce
\begin{align*} \mathcal{P}{\left({W_n}\right)}\leq \frac{\big\|{\widetilde{W}_\infty}\big\|_2^2}{2K} P{\left({K}\right)}+ \frac{\big\|{{\overline{W}}_n}\big\|_2^2}{2K}Q{\left({K}\right)}+o_X{\left({1}\right)}+o_n{\left({1}\right)} \end{align*} |
from (2.17), and recalling that \mathcal{P}{\left({W_n}\right)}\geq P{\left({K}\right)}-o_n{\left({1}\right)} holds by construction, we arrive at
\begin{align} P{\left({K}\right)}+\frac{\left\|{{\overline{W}}_n}\right\|_2^2}{2K}{\big(P{\left({K}\right)}-Q{\left({K}\right)}\big)}\leq \frac{\big\|{\widetilde{W}_\infty}\big\|_2^2+\big\|{{\overline{W}}_n}\big\|_2^2}{2K}P{\left({K}\right)}+o_X{\left({1}\right)}+o_n{\left({1}\right)} \end{align} | (2.18) |
after writing Q{\left({K}\right)} = P{\left({K}\right)}-{\big(P{\left({K}\right)}-Q{\left({K}\right)}\big)} and rearranging terms.
Justification of (2.15): Since the weak convergence \widetilde{W}_n\to \widetilde{W}_\infty combined with (2.16) ensures
\begin{align*} \frac{\big\|{\widetilde{W}_\infty}\big\|_2^2+\big\|{{\overline{W}}_n}\big\|_2^2}{2K}\leq 1+o_n{\left({{1}}\right)}\, \end{align*} |
we find the uniform tightness estimate
\begin{align*} \limsup\limits_{n\to\infty}\left\|{{\overline{W}}_n}\right\|_2^2 = \frac{2K}{P{\left({K}\right)}-Q{\left({K}\right)}}o_X{\left({1}\right)} = o_X{\left({1}\right)} \end{align*} |
because letting n\to\infty for fixed X would otherwise produce a contradiction in (2.18). Combining this with \|{\widetilde{W}_\infty}\|_2\leq \|{W_\infty}\|_2 and letting n\to\infty we thus obtain
\begin{align*} P{\left({K}\right)}\leq \frac{\big\|{W_\infty}\big\|_2^2}{2K}P{\big(K\big)}+o_X{\left({1}\right)}\, , \end{align*} |
and this yields the desired result in the limit X\to\infty .
Concluding remarks: All arguments can be generalized to the continuous-coupling case provided that the truncation in x is accompanied by an appropriate cut-off in \xi -space. More precisely, choosing two parameters 0 < \Xi_1 < \Xi_2 < \infty we write
\begin{align*} \mathcal{P}{\left({W_n}\right)} = \mathcal{P}_{1}{\left({W_n}\right)}+\mathcal{P}_{2}{\left({W_n}\right)}+\mathcal{P}_{3}{\left({W_n}\right)}\, , \end{align*} |
where the three contributions on the right hand side stem from splitting the \xi -integrating in (1.17) into three integral corresponding to \xi\in({0}, {\Xi_1}) , \xi\in({\Xi_1}, {\Xi_2}) , and \xi\in({\Xi_2}, {\infty}) , respectively. Using the identity 0\leq \Phi{\left({r}\right)}\leq \tfrac12 \phi{\left({r}\right)}r^2 , which follows from (2.2), as well as Assumption 1 and Lemma 2 we then estimate
\begin{align*} 0\leq \mathcal{P}_1{\left({W_n}\right)}&\leq K\int\limits_{0}^{\Xi_1} {\alpha}(\xi)\beta^2(\xi)\phi{\left({\xi^{1/2}\beta(\xi) \sqrt{2K}}\right)}\, \mathrm{d}\xi = o_{\Xi_1}{\left({1}\right)}\, , \end{align*} |
where o_{\Xi_1} is independent of n and vanishes for fixed K in the limit \Xi_1\to0 . Moreover, in the same way we verify
\begin{align*} 0\leq \mathcal{P}_3{\left({W_n}\right)}&\leq K\int\limits_{\Xi_2}^\infty {\alpha}(\xi)\beta^2(\xi)\phi{\left({\xi^{1/2}\beta(\xi) \sqrt{2K}}\right)}\, \mathrm{d}\xi = o_{\Xi_2}{\left({1}\right)} \end{align*} |
with o_{\Xi_2}{\left({1}\right)}\to0 as \Xi_2\to\infty . To deal with the remaining term \mathcal{P}_2{\left({W_n}\right)} , we split the x -integration as above and repeat all asymptotic arguments by using a combined error term o_{X, \Xi_1, \Xi_2}{\left({1}\right)} instead of o_X{\left({1}\right)} .
Corollary 8. (existence of solitary waves) The constrained optimization problem (2.7) admits for L = \infty and any K > 0 with (2.14) at least one solution W\in \mathcal{C}_K , which provides a solitary traveling wave.
Proof. The existence of a maximizer is granted by Theorem 7 and Corollary 5 guarantees that this maximizer is a traveling wave.
Notice that any solitary wave is expected to decay exponentially, where the heuristic decay rate {\lambda} depends only on {\sigma} and will be identified below in (3.5) using the imaginary variant of the underlying dispersion relation. The exponential decay can also be deduced rigorously from the traveling wave equation and the properties of the convolution operator \mathcal{A}_\xi , see [23] for the details in an FPUT setting.
The natural strategy to verify (2.14) for sufficiently large K is to fix a test profile W and to show that the function
\begin{align*} \lambda\mapsto \lambda^{-2} \mathcal{P}{\left({ {\lambda} W}\right)} \end{align*} |
is unbounded as {\lambda}\to\infty . A special application of this idea is the following result, which guarantees the existence of solitary waves for all K > 0 under mild assumptions. For simplicity we restrict our considerations to the case of continuous but finite-range interactions and mention that similar results hold if \xi\mapsto\xi^2 {\alpha}(\xi)\beta^2(\xi) decays sufficiently fast or if \Phi_m^{\prime\prime\prime}{\left({0}\right)} > 0 holds for some m in the discrete-coupling case.
Proposition 9. (criterion for genuine super-quadraticity) Condition (2.14) holds for L = \infty in the continuous-coupling case provided that \Phi^{\prime\prime\prime}{\left({0}\right)} > 0 and {{\rm{supp}}}\, { {\alpha}}\subset[{0}, {\Xi}] for some 0 < \Xi < \infty .
Proof. Using Plancharel's theorem and the properties of the \mathrm{sinc} function we estimate
\begin{align} \begin{split} Q{\left({K}\right)}& = \sup\limits_{W\in \mathcal{C}_K}\tfrac{1}{2}\Phi^{\prime\prime}{\left({0}\right)} \int\limits_0^{\Xi}\xi^2 {\alpha}(\xi) {\beta}^2{\left({\xi}\right)}\int\limits_{ {\mathbb{R}}}\mathrm{sinc}^2{\big(k\xi/2\big)}\big|{\widehat{W}{\left({k}\right)}}\big|^2\, \mathrm{d}k{\rm d}\xi \\&\leq\Phi^{\prime\prime}{\left({0}\right)} K \int\limits_0^{\Xi}\xi^2 {\alpha}(\xi) {\beta}^2{\left({\xi}\right)}{\rm d}\xi\, , \end{split} \end{align} | (2.19) |
and define for any large \Lambda a piecewise constant function W_\Lambda\in \mathcal{C} by
\begin{align*} W_\Lambda{\left({x}\right)} = \sqrt{\frac{K}{\Lambda}}\chi_{I_\Lambda}{\left({x}\right)}\, . \end{align*} |
Since the unimodal function \mathcal{A}_\xi W_\Lambda attains the value \xi K^{1/2}\Lambda^{-1/2} for \left|{x}\right|\leq \Lambda-\xi/2 and vanishes for \left|{x}\right|\geq \Lambda+\xi/2 , we arrive at the lower bound
\begin{align*} \mathcal{Q}{\left({W_\Lambda}\right)}\geq \tfrac{1}{2}\Phi^{\prime\prime}{\left({0}\right)} \int\limits_0^{\Xi} {\alpha}(\xi) {\beta}^2{\left({\xi}\right)}\int\limits_{-\Lambda+\xi/2}^{+\Lambda-\xi/2} \xi^2 K \Lambda^{-1}\, \mathrm{d}x\, \mathrm{d}\xi \end{align*} |
and conclude with \Lambda\to\infty that the upper bound in (2.19) is actually sharp. Consequently, we have
\begin{align*} \mathcal{Q}{\left({W_\Lambda}\right)}\geq Q{\left({K}\right)}-C\Lambda^{-1} \end{align*} |
for some positive constant C independent of \Lambda . Moreover, using the smoothness of \Phi and thanks to Lemma 2 we estimate
\begin{align*} \mathcal{P}{\left({W_\Lambda}\right)}- \mathcal{Q}{\left({W_\Lambda}\right)}&\geq {\left({1-o_\Lambda{\left({1}\right)}}\right)} \tfrac16 \Phi^{\prime\prime\prime}{\left({0}\right)} \int\limits_{0}^{\Xi} {\alpha}(\xi) {\beta}^3(\xi)\int\limits_{-\Lambda+\xi/2}^{+\Lambda-\xi/2}{\big({\left({ \mathcal{A}_\xi W_\Lambda}\right)}{\left({x}\right)}\big)}^3\, \mathrm{d}x{\rm d}\xi \\&\geq {\left({1-o_\Lambda{\left({1}\right)}}\right)} \tfrac16 \Phi^{\prime\prime\prime}{\left({0}\right)} \int\limits_{0}^{\Xi} \xi^3 {\alpha}(\xi) {\beta}^3(\xi) \frac{K^{3/2}{\left({2 \Lambda -\xi}\right)}}{\Lambda^{3/2}} {\rm d}\xi\ \, , \end{align*} |
where the small quantity o_\Lambda{\left({1}\right)} is chosen such that
\begin{align*} \sup\limits_{0\leq\xi\leq\Xi}\;\sup\limits_{0\leq r\leq \xi^{1/2} {\beta}{\left({\xi}\right)}\sqrt{2K}}\left|{\frac{\Phi^{\prime\prime\prime}{\left({r}\right)}}{\Phi^{\prime\prime\prime}{\left({0}\right)}}}\right|\geq 1-o_\Lambda{\left({1}\right)}\, . \end{align*} |
In summary, we obtain
\begin{align*} P{\left({K}\right)}\geq \mathcal{P}{\left({W_\Lambda}\right)}- \mathcal{Q}{\left({W_\Lambda}\right)}+\mathcal{Q}{\left({W_\Lambda}\right)}\geq Q{\left({K}\right)}-C\Lambda^{-1}+c \Lambda^{-1/2} \end{align*} |
for every large \Lambda and some constant c > 0 , so the claim follows by choosing \Lambda sufficiently large.
Similar ideas reveal for large L and K that the constant function cannot be the maximizer for problem (2.7) provided that \Phi grows sufficiently fast. In fact, we have
\begin{align} \mathcal{P}{\left({K^{1/2}\Lambda^{-1/2}}\right)} = 2L \int\limits_0^{\Xi}\Phi{\left({\beta{\left({\xi}\right)}\xi K^{1/2}\Lambda^{-1/2}}\right)}{\rm d}\xi \end{align} | (2.20) |
and verify for L > \Lambda+\Xi the estimate
\begin{align} \mathcal{P}{\left({W_\Lambda}\right)}&\geq {\left({2M-\Xi}\right)} \int\limits_0^{\Xi} {\alpha}(\xi)\Phi{\left({\beta{\left({\xi}\right)}\xi K^{1/2}\Lambda^{-1/2}}\right)}{\rm d}\xi\, . \end{align} | (2.21) |
The right hands side in (2.21), however, is typically larger than that from (2.20) provided that r^{-2}\Phi{\left({r}\right)}\to\infty as r\to\infty .
We finally show that the solitary waves provided by Corollary 8 can be approximated by periodic ones. To this end we write the L -dependence of P and \mathcal{P} explicitly.
Proposition 10. (convergence of maximizers) Suppose that (2.14) is satisfied. Then we have
\begin{align} P_L{\left({K}\right)}\quad\xrightarrow{L\to\infty}\quad P_\infty{\left({K}\right)}\, . \end{align} | (2.22) |
Moreover, any family {\left({W_L}\right)}_{L\geq 1} with W_L being a maximizer of problem (2.7) on I_L converges along subsequences and in a strong {\mathsf{L}}^2 -sense to a maximizer W_\infty of problem (2.7) on I_\infty = {\mathbb{R}} .
Proof. Convergence of maxima: For finite L , we denote by
\begin{align*} \widetilde{W}_L{\left({x}\right)}: = W_L{\left({x}\right)}\quad \text{for}\quad \left|{x}\right| \lt L\, , \qquad \widetilde{W}_L{\left({x}\right)}: = 0\quad\text{else} \end{align*} |
the trivial continuation of W_L to a function from {\mathsf{L}}^2({\mathbb{R}}) , which satisfies \widetilde{W}_L\in \mathcal{C} as well as \mathcal{K}_\infty{\big(\widetilde{W}_L\big)} = \mathcal{K}_L{\left({W_L}\right)} . Moreover, given a solitary maximizer W_\infty we define for any 0 < L < \infty the 2L -periodic function {\overline{W}}_L to be the scaled periodic continuation of W_\infty restricted to I_L . This reads
\begin{align*} {\overline{W}}_L{\left({x}\right)}: = c_L\, W_\infty{\left({ x}\right)}\quad \text{for}\quad \left|{x}\right| \lt L\, , \qquad {\overline{W}}_L{\left({x+2L}\right)} = {\overline{W}}_L{\left({x}\right)}\quad \text{for}\quad x\in {\mathbb{R}} \end{align*} |
with
\begin{align*} c_L^{-2}: = {\left({2K}\right)}^{-1} \int\limits_{I_L}\big|{W_\infty{\left({x}\right)}}\big|^2\, \mathrm{d}x \end{align*} |
and implies {\overline{W}}_L\in \mathcal{C} as well as \mathcal{K}_L{\big({\overline{W}}_L\big)} = K . Since the kernel function in the convolution operator \mathcal{A}_\xi has compact support and since the decay estimate
\begin{align*} \big|{W_L{\left({x}\right)}}\big|\leq \sqrt{\frac{K}{\left|{x}\right|}} \end{align*} |
holds by (2.3) both for finite and infinite L , we readily show
\begin{align*} \big|{ \mathcal{P}_\infty{\big(\widetilde{W}_L\big)}- \mathcal{P}_L{\left({W_L}\right)}}\big|\leq o_L{\left({1}\right)}\, , \end{align*} |
with o_L{\left({1}\right)}\to0 as L\to\infty , and deduce
\begin{align} P_\infty{\left({K}\right)}\geq \limsup\limits_{L\to\infty} \mathcal{P}_\infty{\big(\widetilde{W}_L\big)}\geq \limsup\limits_{L\to\infty} \mathcal{P}_L{\big(W_L\big)} = \limsup\limits_{L\to\infty} P_L{\left({K}\right)}\, . \end{align} | (2.23) |
Similarly, we find
\begin{align*} \big|{ \mathcal{P}_L{\big({\overline{W}}_L\big)}- \mathcal{P}_\infty{\left({W_\infty}\right)}}\big|\leq o_L{\left({1}\right)}\, \end{align*} |
thanks to \left\|{{\overline{W}}_L-\chi_{I_L} W_\infty}\right\|_{ {\mathsf{L}}^2{\left({I_L}\right)}}\to0 , and this implies
\begin{align} \liminf\limits_{L\to\infty} P_L{\left({K}\right)} \geq \liminf\limits_{L\to\infty} \mathcal{P}_L{\big({\overline{W}}_L\big)} = \mathcal{P}_\infty{\left({W_\infty}\right)} = P_\infty{\left({K}\right)}\, . \end{align} | (2.24) |
The claim (2.22) now follows from (2.23) and (2.24).
Convergence of maximizers: By construction and (2.22) we have \mathcal{P}_\infty{\big(\widetilde{W}_L\big)}\to P{\left({K}\right)} , so Theorem 7 ensures the existence of strongly {\mathsf{L}}^2 -convergent subsequences. In particular, any accumulation point W_\infty belongs to \mathcal{C} and satisfies \mathcal{K}_\infty{\left({W_\infty}\right)} = K as well as \mathcal{P}_\infty{\big(W_\infty\big)} = P{\left({K}\right)} .
Proposition 10 implies that the periodic solutions to the optimization problem (2.7) come with non-constant profile functions provided that L is sufficiently large and that condition (2.14) is satisfied. We further mention that the techniques from the proof of Proposition 10 can be used to characterize the continuous-coupling case as the scaling limit of the discrete coupling case.
In this section we briefly discuss some asymptotic properties of traveling waves in peridynamical media. We focus again on the continuous-coupling case (1.4) but the essential arguments can easily be modified to cover the lattice case (1.6) as well. To ease the notation we suppose
\begin{align*} \Phi^{\prime\prime}{\left({0}\right)} = 1 \end{align*} |
and assume that the smooth and nonnegative coupling coefficients {\alpha} and \beta sufficiently regular for \xi\to 0 and \xi\to\infty so that all integrals appearing below are well-defined.
In the harmonic case \Phi{\left({r}\right)} = \tfrac12r^2 we can solve the linear equation (1.13) by Fourier transform. More precisely, with
\begin{align*} W_{n, L, K}: = \sqrt{2KL^{-1}}\cos{\big( n L^{-1}\pi x\big)}\, , \qquad {\sigma}^2_{n, L}: = \Theta^2{\big(n L^{-1}\pi\big)} \end{align*} |
and
\begin{align} \Theta^2{\left({k}\right)}: = \int\limits_0^\infty \xi^2\alpha{\left({\xi}\right)}\beta^2(\xi) \mathrm{sinc}^2{\left({k\xi/2}\right)}{\rm d}\xi\, . \end{align} | (3.1) |
we obtain the family of all even and periodic solutions, where the parameters n\in {\mathbb{N}} , K > 0 , and 0 < L < \infty can be chosen arbitrarily. Of course, the waves for n > 0 are neither unimodal nor positive and not relevant for the constrained optimization problem from §2 since we have
\begin{align*} \mathcal{P}_L{\big(W_{n, K, L}\big)} = \mathcal{Q}_L{\big(W_{n, K, L}\big)} = \Theta^2{\big(nL^{-1}\pi\big)}K \end{align*} |
but
\begin{align} P_L{\left({K}\right)} = Q_L{\left({K}\right)} = K\int\limits_0^\infty\xi^2 a(\xi)\beta^2{\left({\xi}\right)}{\rm d}\xi = \Theta^2{\left({0}\right)}K \end{align} | (3.2) |
for the harmonic standard potential. The dispersion relation for the peridynamical wave equation (1.1) is given by
\begin{align} {\omega}^2 = {\Omega}^2{\left({k}\right)}: = k^2\Theta^2{\left({k}\right)} \end{align} | (3.3) |
and follows from inserting the ansatz u({t}, {x}) = \exp{\left({ \mathtt{i} kx- \mathtt{i} {\omega} t}\right)} in the linear variant of the dynamical model (1.3). Elementary asymptotic arguments reveal
\begin{align*} \Theta^2{\left({k}\right)}\quad\xrightarrow{k\to0}\quad c_0: = \int\limits_0^\infty \xi^2 {\alpha}(\xi)\beta^2(\xi){\rm d}\xi\, , \qquad \Omega^2{\left({k}\right)}\quad\xrightarrow{k\to\pm\infty}\quad c_\infty: = 2\int\limits_0^\infty {\alpha}(\xi)\beta^2(\xi){\rm d}\xi \end{align*} |
as well as
\begin{align*} \Theta^2{\left({k}\right)}\quad\xrightarrow{k\to0}\quad\int\limits_0^\infty \xi^2 {\alpha}(\xi)\beta^2(\xi){\rm d}\xi+O{\big(k^2\big)}, \end{align*} |
provided that the coefficient functions {\alpha} and \beta from (1.4) are sufficiently non-singular such that {\alpha}\beta^2 is integrable in [{0}, {\infty}) . We refer to Figure 5 for an illustration and mention that the choice in (1.5) is less regular due to {\alpha}(\xi)\beta^2{\left({\xi}\right)} = \xi^{-1} for \xi\approx0 .
The simplest linear PDE model that comprises the same asymptotic properties in its dispersion relation is
\begin{align*} c_0\partial_t^2 \partial_x^2 u({t}, {x})-c_\infty\partial_t^2u({t}, {x})+c_0c_\infty \partial_x^2u({t}, {x}) = 0\, . \end{align*} |
This equation is well-posed as it can be written as the Banach-valued ODE
\begin{align*} \partial_t^2u({t}, {x}) = \mathcal{D} u({t}, {x})\, , \qquad \mathcal{D} = {\big(c_0^{-1}-c_{\infty}^{-2}\partial_x^2\big)}^{-1}\partial_x^2\, , \end{align*} |
where the bounded pseudo-differential operator \mathcal{D} has the symbol function
\begin{align*} \widehat{ \mathcal{D}}{\left({k}\right)} = \frac{c_0c_\infty k^2}{c_\infty+c_0k^2}\, . \end{align*} |
A nonlinear analogue would be the PDE
\begin{align} \partial_t^2u({t}, {x}) = {\big(c_0^{-1}-c_{\infty}^{-2}\partial_x^2\big)}^{-1}\partial_x\Psi^\prime{\big(\partial_xu({t}, {x})\big)}\, , \end{align} | (3.4) |
which can be viewed as a variant of the regularized Boussinesq equation, see for instance [35]. The corresponding traveling wave ODE
\begin{align*} -c_0W^{\prime\prime} = {\sigma}^{-2}c_0c_\infty\Psi^\prime{\left({W}\right)}-c_\infty W \end{align*} |
is of Hamiltonian type and admits homoclinic solutions provided that \Psi grows super-linearly and that {\sigma} is sufficiently large. It is desirable to explore the similarities and differences between the peridynamical model (1.1) and the PDE substitute (3.4) in greater detail.
We finally observe that the complexified version of \Theta^2 predicts the spatial decay of unimodal solitary waves with non-harmonic \Phi . In fact, the exponential ansatz
\begin{align} W{\left({x}\right)}\sim \exp{\left({- {\lambda}\left|{x}\right|}\right)}\quad\text{for}\quad \left|{x}\right|\to\infty \end{align} | (3.5) |
is compatible with (1.13) if and only if the rate parameter {\lambda} satisfies the transcendental equation
\begin{align*} {\sigma}^2 = \int\limits_0^\infty \xi^2\alpha{\left({\xi}\right)}\beta^2(\xi) \, \mathrm{sinch}^2{\left({\xi {\lambda}/2}\right)}{\rm d}\xi = \Theta^2{\left({ \mathtt{i} {\lambda}}\right)}\, . \end{align*} |
This equation admits a unique positive and real solution for {\sigma}^2 > \Theta^2{\left({0}\right)} , and all solitary waves from §2 meet this condition due to the lower bound for the wave speed in Corollary 5 and since Proposition 10 combined with (3.2) implies K^{-1}P_\infty{\left({K}\right)} > K^{-1}Q_\infty{\left({K}\right)} = \Theta^2{\left({0}\right)} .
It is well-known that the Korteweg-deVries (KdV) equation is naturally related to many nonlinear and spatially extended Hamiltonian systems as it governs the effective dynamics in the large wave-length regime, in which traveling waves propagate with near sonic speed, have small amplitudes, and carry small energies. For solitary waves in FPUT chains, this was first observed in [49] and has later been proven rigorously in [9]. We also refer to [8] for periodic KdV waves, to [3] for two-dimensional lattices, and to [14,28] for the more complicated case of polyatomic chains. Related results on initial value problems can be found in [25,26,38].
The KdV equation also governs the near sonic limit of nonlocal lattices and the corresponding existence problem for solitary lattice wave has been investigated rigorously in [22]. In what follows we discuss how the underlying ideas and asymptotic techniques can be applied in the peridynamical setting (1.4) provided that
\begin{align} \eta: = \tfrac12\Phi^{\prime\prime\prime}{\left({0}\right)} \gt 0\, . \end{align} |
The starting point for the construction of KdV waves is the ansatz
\begin{align} W{\left({x}\right)} = {\varepsilon}^2 {\overline{W}}{\big( {\varepsilon} x\big)}\, , \qquad {\sigma}^2 = \Theta^2{\left({0}\right)}+{\varepsilon}^2 \end{align} | (3.6) |
with sound speed \Theta{\left({0}\right)} from (3.1), natural scaling parameter {\varepsilon} > 0 , and rescaled space variable {\overline{x}}: = {\varepsilon}{x} . Inserting this ansatz into (1.13), using the integral identity
\begin{align*} {\big( \mathcal{A}_\xi W\big)}{\left({x}\right)} = {\varepsilon} \int\limits_{{\overline{x}}- {\varepsilon}\xi/2}^{{\overline{x}}+{\varepsilon}\xi/2}{\overline{W}}{\big({\overline{y}}\big)}\, \mathrm{d}\overline{y} = {\big( {\varepsilon} \mathcal{A}_{ {\varepsilon}\xi}{\overline{W}}\big)}{\left({{\overline{x}}}\right)}\, , \end{align*} |
and dividing by {\varepsilon}^4 we obtain a transformed integral equation for {\overline{W}} . This equations can be written as
\begin{align} \mathcal{B}_ {\varepsilon} V = \mathcal{N}_ {\varepsilon}{\left({V}\right)}+{\varepsilon}^2 \mathcal{G}_ {\varepsilon}{\left({V}\right)}\, , \end{align} | (3.7) |
where the operators
\begin{align*} \mathcal{B}_ {\varepsilon} {\overline{W}}: = {\overline{W}}+{\varepsilon}^{-2}\Theta^2{\left({0}\right)}{\overline{W}}- {\varepsilon}^{-4}\int\limits_0^{\infty} {\alpha}(\xi)\beta^2(\xi) \mathcal{A}_{ {\varepsilon}\xi}^2{\overline{W}}{\rm d}\xi \end{align*} |
and
\begin{align*} \mathcal{N}_ {\varepsilon}{\big({\overline{W}}\big)}: = {\varepsilon}^{-5}\eta\int\limits_0^{\infty} {\alpha}(\xi)\beta^3(\xi) \mathcal{A}_{ {\varepsilon}\xi}{\big( {\varepsilon} \mathcal{A}_{ {\varepsilon}\xi}{\overline{W}}\big)}^2{\rm d}\xi \end{align*} |
collect all terms that are linear and quadratic, respectively, with respect to {\overline{W}} , while
\begin{align*} \mathcal{G}_ {\varepsilon}{\big({\overline{W}}\big)}: = {\varepsilon}^{-7}\int\limits_0^{\infty} {\alpha}(\xi)\beta(\xi) \mathcal{A}_{ {\varepsilon}\xi}G{\big( {\varepsilon}\beta(\xi) \mathcal{A}_{ {\varepsilon}\xi}{\overline{W}}\big)}{\rm d}\xi \end{align*} |
with G{\left({r}\right)}: = \Phi^\prime{\left({r}\right)}-r-\eta r^2 = O{\big(r^3\big)} stems from the cubic and the higher order terms and does not contribute to the KdV limit.
The key asymptotic observations for the limit {\varepsilon}\to0 is the formal asymptotic expansion
\begin{align} \mathcal{A}_\delta U = \delta U+\tfrac{1}{24}\delta^3U^{\prime\prime}+{...}\, , \end{align} | (3.8) |
which provides in combination with the above formulas the formal limit equation
\begin{align*} {\overline{W}}-c_1{\overline{W}}^{\prime\prime} = c_2 {\overline{W}}^2 \end{align*} |
with positive coefficients
\begin{align} c_1 : = \tfrac12\int\limits_0^\infty\xi^4 {\alpha}(\xi)\beta^2(\xi){\rm d}\xi \gt 0\, , \qquad c_2: = \eta\int\limits_0^\infty\xi^3 {\alpha}(\xi)\beta^3(\xi){\rm d}\xi \gt 0\, . \end{align} | (3.9) |
This planar Hamiltonian ODE is precisely the traveling wave equation of the KdV equation and admits the homoclinic and even solution
\begin{align} {\overline{W}}_{\text{KdV}}{\left({{\overline{x}}}\right)}: = \frac{3}{2c_2}\mathrm{sech}^2{\left({\frac{{\overline{x}}}{2\sqrt{c_1}}}\right)}\, , \end{align} | (3.10) |
which is moreover unique, unimodal, positive, and exponentially decaying. However, the asymptotic expansion in (3.8) is not regular but singular because the error terms involve higher derivatives. The rigorous justification of the limit problem (3.9) is therefore not trivial and necessitates the use of refined asymptotic techniques. Following [9], we use the corrector ansatz
\begin{align} {\overline{W}} = {\overline{W}}_{\text{KdV}}+{\varepsilon}^2{\overline{V}}\, , \end{align} | (3.11) |
which transforms (3.7) into
\begin{align*} \mathcal{B}_ {\varepsilon} {\overline{V}}- \mathcal{M}_ {\varepsilon} {\overline{V}}: = \mathcal{F}_ {\varepsilon}{\big({\overline{V}}\big)}\, , \end{align*} |
where the linear operator
\begin{align*} \mathcal{M}_ {\varepsilon} {\overline{W}} = {\varepsilon}^{-5}2\eta\int\limits_0^{\infty} {\alpha}(\xi)\beta^3(\xi) \mathcal{A}_{ {\varepsilon}\xi}{\Big({\big( {\varepsilon} \mathcal{A}_{ {\varepsilon}\xi}{\overline{W}}_{\text{KdV}}\big)}{\big( {\varepsilon} \mathcal{A}_{ {\varepsilon}\xi}{\overline{V}}\big)}\Big)}{\rm d}\xi \end{align*} |
represents the additional linear term that stems from the linearization of \mathcal{N}_ {\varepsilon} around {\overline{W}}_{\text{KdV}} . Moreover, the nonlinear operator
\begin{align*} \mathcal{F}_ {\varepsilon}{\big({\overline{V}}\big)} = {\varepsilon}^{-2}{\Big( \mathcal{B}_ {\varepsilon} {\overline{W}}_{\text{KdV}}- \mathcal{Q}_ {\varepsilon}{\big({\overline{W}}_{\text{KdV}}\big)}\Big)}+{\varepsilon}^2 \mathcal{Q}_ {\varepsilon}{\big({\overline{V}}\big)}+\mathcal{G}_ {\varepsilon}{\big({\overline{W}}_{\text{KdV}}+{\varepsilon}^2{\overline{V}}\big)} \end{align*} |
is composed of small nonlinear terms in {\overline{V}} as well as residual terms depending on {\varepsilon} and {\overline{W}}_{\text{KdV}} only.
The rigorous justification of the KdV limit can be achieved along the lines of [22] and hinges on the following crucial ingredients:
1. The pseudo-differential operator \mathcal{B}_ {\varepsilon} is uniformly invertible and its inverse is almost compact in the sense that it can be written as the sum of a compact operator and a small bounded one. These statements follow from an asymptotic analysis of the corresponding symbol function
\begin{align} \widehat{ \mathcal{B}}_ {\varepsilon} {\left({k}\right)}: = 1+\int\limits_0^{\infty}\xi^2 {\alpha}(\xi)\beta^2(\xi){\big(1- {\varepsilon}^{-2}\mathrm{sinc}^2{\big( {\varepsilon} \xi k/2\big)}\big)}{\rm d}\xi\, , \end{align} | (3.12) |
see Figure 6 for an illustration.
2. The symmetric operator \mathcal{L}_ {\varepsilon}: = \mathcal{B}_ {\varepsilon}- \mathcal{M}_ {\varepsilon} is uniformly invertible on {\mathsf{L}}^2_{\text{even}}({\mathbb{R}}) as it satisfies
\begin{align*} \big\|{ \mathcal{L}_ {\varepsilon} {\overline{V}}}\big\|_2\geq c \big\|{{\overline{V}}}\big\|_2 \end{align*} |
on that space, where the constant c > 0 can be chosen independently of {\varepsilon} . This can be proven by contradiction using the almost compactness of \mathcal{B}_ {\varepsilon}^{-1} and the fact, that the limiting operator
\begin{align*} \mathcal{L}_0{\overline{V}} = {\overline{V}}-c_1{\overline{V}}^{\prime\prime}-2c_2 {\overline{W}}_{\text{KdV}}{\overline{V}} \end{align*} |
has a one-dimensional nullspace in {\mathsf{L}}^2({\mathbb{R}}) according to the Sturm-Liouville theory, which is spanned by the odd derivative of the KdV profile (3.10).
3. For sufficiently small {\varepsilon} > 0 , the operator \mathcal{L}_ {\varepsilon}^{-1} \mathcal{F}_ {\varepsilon} maps a certain ball in {\mathsf{L}}^2_ {{\rm{\, even}}}({\mathbb{R}}) contractively into itself so that there exists a locally unique fixed point.
In particular, for any sufficiently small {\varepsilon} > 0 there exists via (3.6) and (3.11) a near-sonic wave solution to (1.13). The nonlinear stability for such waves has been shown in the paper series [10,11,12] for FPUT chains, see also [27]. It remains a challenging task to generalize the proofs to the peridynamical case.
Another well-known asymptotic regime concerns atomic chains with super-polynomially growing interaction potential and waves which propagate much faster than the sound speed and carry a huge amount of energy. It has been observed in [7,45] for FPUT chains with singular Lennard-Jones-types potential that both the distance and the velocity profiles of those high-speed waves converge to simple limit functions which are naturally related to the traveling waves in the hard-sphere model. In other words, the particles interact asymptotically by elastic collisions and traveling waves come with strongly localized profile functions. The asymptotic analysis of the high-energy limit of FPUT chains and the underlying advance-delay-differential equation has later been refined by the authors. In [17,18] they derive accurate and almost explicit approximation formulas for the wave profiles by combining a nonlinear shape ODE with local scaling arguments and natural matching conditions. Moreover, the local uniqueness, smooth parameter dependence, and nonlinear orbital stability of high-speed waves are established in [19,20] using similar two-scale techniques, the non-asymptotic part of the Friesecke-Pego theory from [10,11], and the enhancement by Mizumachi in [33].
Preliminary computations as well as numerical simulations indicate that the high-energy waves in peridynamical media can likewise be approximated by simple profile functions. For two examples with discrete coupling we refer to Figure 7, which relies on
\begin{align} M = 5\, , \qquad L = 10\, , \qquad \Phi_m{\left({r}\right)} = \cosh{\left({\beta_m r}\right)}\, , \end{align} | (3.13) |
and illustrates that the wave profile W converges as K\to\infty to the scaled indicator function of an interval provided that the interaction potentials grow faster than a polynomial. This implies that the corresponding distance profiles are basically piecewise affine but the details depend on the choice of the parameters.
A rigorous asymptotic analysis of this numerical observation lies beyond the scope of this paper but first steps are already made in [2] for the nonlocal advance-delay-differential equations that govern solitary waves in two-dimensional FPUT lattices. More precisely, assuming exponentially growing potentials and using the multi-scale techniques from [17] one can derive a reduced Hamiltonian ODE system that can be expected to describe the local rescaling of the different distance profiles. This ODE, however, has more than two degrees of freedom and a qualitative or quantitative analysis of the relevant solutions is hence not trivial. It remains a challenging task to characterize the high-energy limit for peridynamical media or chains with more than nearest-neighbor interactions, and to justify the numerical observations in Figure 7.
Some aspects of our work build upon analytical or numerical investigations in the theses of Benedikt Hewer and Fanzhi Chen, see [1,24]. We also acknowledge the financial support by the Deutsche Forschungsgemeinschaft (DFG individual grant HE 6853/2-1).
The authors declare no conflict of interest.
[1] | BRASIL. Estimativa 2016: Incidência de Câncer no Brasil. INCA 2015. |
[2] |
Torre LA, Bray F, Siegel RL, et al. (2015) Global Cancer Statistics, 2012. CA a cancer J Clin 65: 87-108. doi: 10.3322/caac.21262
![]() |
[3] | Walboomers JMM, Jacobs MV, Manos MM, et al. (1999) Human Papillomavirus Is a Necessary Cause Of Invasive Cervical Cancer Worldwide. J Patholgy 189: 12-19. |
[4] |
Woodman CBJ, Collins SI, Young LS (2007) The natural history of cervical HPV infection: unresolved issues. Nat Rev Cancer 7: 11-22. doi: 10.1038/nrc2050
![]() |
[5] |
Smith JS, Lindsay L, Hoots B, et al. (2007) Human papillomavirus type distribution in invasive cervical cancer and high-grade cervical lesions: A meta-analysis update. Int J Cancer 121: 621-632. doi: 10.1002/ijc.22527
![]() |
[6] |
Lee LJ, Howitt B, Catalano P, et al. (2016) Gynecologic Oncology Prognostic importance of human papillomavirus (HPV) and p16 positivity in squamous cell carcinoma of the vulva treated with radiotherapy. Gynecol Oncol 142: 293-298. doi: 10.1016/j.ygyno.2016.05.019
![]() |
[7] |
Levovitz C, Chen D, Ivansson E, et al. (2014) TGFβ Receptor 1: An immune susceptibility gene in HPV-associated cancer. Cancer Res 74: 6833-6844. doi: 10.1158/0008-5472.CAN-14-0602-T
![]() |
[8] |
Mai S, Welzel G, Ottstadt M, et al. (2015) Prognostic relevance of HPV infection and p16 overexpression in squamous cell anal cancer. Radiat Oncol Biol 93: 819-827. doi: 10.1016/j.ijrobp.2015.08.004
![]() |
[9] |
Egawa N, Egawa K, Griffin H, et al. (2015) Human Papillomaviruses; epithelial tropisms, and the development of neoplasia. Viruses 7: 3863-3890. doi: 10.3390/v7072802
![]() |
[10] | Arbyn M, Sanjosé S de, Saraiya M, et al. (2011) EUROGIN 2011 roadmap on prevention and treatment of HPV- related disease. Int J Cancer 131: 1969-1982. |
[11] |
Doorbar J (2006) Molecular biology of human papillomavirus infection and cervical cancer. Clin Sci 110: 525-541. doi: 10.1042/CS20050369
![]() |
[12] |
Moody CA, Laimins LA (2010) Human papillomavirus oncoproteins : pathways to transformation. Nat Rev Cancer 10: 550-560. doi: 10.1038/nrc2886
![]() |
[13] |
Doorbar J (2013) The E4 protein —structure , function and patterns of expression. Virology 445: 80-98. doi: 10.1016/j.virol.2013.07.008
![]() |
[14] |
DiMaio D, Petti LM (2013) The E5 proteins. Virology 445: 99-114. doi: 10.1016/j.virol.2013.05.006
![]() |
[15] | Beaudenon S, Huibregtse JM (2008) HPV E6 , E6AP and cervical cancer. BMC Biochem 7: 1-7. |
[16] | Zhou W, Slingerland JM (2014) Links between oestrogen receptor activation and proteolysis : relevance to hormone-regulated cancer therapy. Nat Rev Cancer 14: 26-38. |
[17] |
zur Hausen H (2009) Papillomaviruses in the causation of human cancers—a brief historical account. Virology 384: 260-265. doi: 10.1016/j.virol.2008.11.046
![]() |
[18] |
Schiffman M, Castle PE, Jeronimo J, et al. (2007) Human papillomavirus and cervical cancer. Lancet 370: 890-907. doi: 10.1016/S0140-6736(07)61416-0
![]() |
[19] |
Schiffman M, Wentzensen N (2013) Human Papillomavirus Infection and the Multistage Carcinogenesis of Cervical Cancer. Cancer Epidemiol Biomarkers Prev 22: 553-560. doi: 10.1158/1055-9965.EPI-12-1406
![]() |
[20] | Castellsagué X, Muñoz N (2003) Chapter 3: Cofactors in human papillomavirus carcinogenesis—role of parity, oral contraceptives, and tobacco smoking. J Nation Cancer Institute Monographs 31: 20-28. |
[21] | Patel S, Chiplunkar S (2009) Host immune responses to cervical cancer. Curr Opin Obstet Gynecol 21: 54-59. |
[22] | Mbeunkui F, Johann DJJr (2010) Cancer and the tumor microenvironment: a review of an essential relationship. Cancer Chemoter Pharmacol 63: 571-582. |
[23] | Vandercappellen J, Van Damme J, Struyf S (2008) The role of CXC chemokines and their receptors in cancer. Cancer Lett 267: 226-244. |
[24] | Baggiolini M (2001) Chemokines in pathology and medicine. J Intern Med 250: 91-104. |
[25] |
Griffith JW, Sokol CL, Luster AD (2014) Chemokines and chemokine receptors: positioning cells for host defense and immunity. Annu Rev Immunol 32: 659-702. doi: 10.1146/annurev-immunol-032713-120145
![]() |
[26] |
Pozzobon T, Goldoni G, Viola A, et al. (2016) CXCR4 signaling in health and disease. Immunol Lett 177: 6-15. doi: 10.1016/j.imlet.2016.06.006
![]() |
[27] | Hanahan D, Weinberg RA (2011) Hallmarks of cancer: the next generation. Cell 144: 646-674. |
[28] | Wani N, Nasser MW, Ahirwar DK, et al. (2014) C-X-C motif chemokine 12/C-X-C chemokine receptor type 7 signaling regulates breast cancer growth and metastasis by modulating the tumor microenvironment. Breast Cancer Res16: R54. |
[29] |
Bergers G, Benjamin LE (2003) Tumorigenesis and the angiogenic switch. Nat Rev Cancer 3: 401-410. doi: 10.1038/nrc1093
![]() |
[30] |
Tachibana K, Hirota S, Iizasa H, et al. (1998). The chemokine receptor CXCR4 is essential for vascularization of the gastrointestinal tract. Nature 393: 591-594. doi: 10.1038/31261
![]() |
[31] |
Funakoshi A, Jimi A, Yasunami Y, et al. (1989) CXCL8/IL8 stimulates vascular endothelial growth factor (VEGF) expression and the autocrine activation of VEGFR2 in endothelial cells by activating NFkappaB through the CBM (Carma3/Bcl10/Malt1) complex. Biochem Biophys Res Commun 159: 913-918. doi: 10.1016/0006-291X(89)92195-5
![]() |
[32] |
Martin D, Galisteo R, Gutkind JS (2009) CXCL8/IL8 stimulates vascular endothelial growth factor (VEGF) expression and the autocrine activation of VEGFR2 in endothelial cells by activating NF-κB through the CBM (Carma3/Bcl10/Malt1) complex. J Biol Chem 284: 6038-6042. doi: 10.1074/jbc.C800207200
![]() |
[33] |
Nguyen DX, Bos PD, Massagué J (2009) Metastasis: from dissemination to organ-specific colonization. Nat Rev Cancer 9: 274-284. doi: 10.1038/nrc2622
![]() |
[34] | Zlotnik A, Burkhardt AM, Homey B (2011) Homeostatic chemokine receptors and organ-specific metastasis. Nat Rev Immunol 11: 597-606. |
[35] | Muller A, Homey B, Soto H, et al. (2001) Involvement of chemokine receptors in breast cancer metastasis. Nature 410: 50-56. |
[36] | Terasaki M, Sugita Y, Arakawa F, et al. (2011) CXCL12/CXCR4 signaling in malignant brain tumors: a potential pharmacological therapeutic target. Brain Tumor Pathol 28: 89-97. |
[37] |
Balkwill F (2004) Cancer and the chemokine network. Nat Rev Cancer 4: 540-550. doi: 10.1038/nrc1388
![]() |
[38] |
Zou W, Machelon V, Coulomb-L'Hermin A, et al. (2001) Stromal-derived factor-1 in human tumors recruits and alters the function of plasmacytoid precursor dendritic cells. Nat Med 7: 1339-1346. doi: 10.1038/nm1201-1339
![]() |
[39] |
Maksym RB, Tarnowski M, Grymula K, et al. (2009) The role of stromal derived factor-1-CXCR7 axis in development and cancer. Eur J Pharmacol 625: 31-40. doi: 10.1016/j.ejphar.2009.04.071
![]() |
[40] | Wang J, Knaut H (2014) Chemokine signaling in development and disease. Development 12: 4199-4205. |
[41] |
Wang Z, Ma Q, Liu Q, et al. (2008) Blockade of SDF-1/CXCR4 signalling inhibits pancreatic cancer progression in vitro via inactivation of canonical Wnt pathway. Br J Cancer 99: 1695-1703. doi: 10.1038/sj.bjc.6604745
![]() |
[42] |
Ganju RK, Brubaker SA, Meyer J, et al. (1998) The a-Chemokine, stromal cell-derived factor-1a, binds to the transmembrane G-protein-coupled CXCR-4 receptor and activates multiple signal transduction pathways. J Biol Chem 273: 23169-23175. doi: 10.1074/jbc.273.36.23169
![]() |
[43] |
Bhandari D, Robia SL, Marchese A (2009) The E3 ubiquitin ligase atrophin interacting protein 4 binds directly to the chemokine receptor CXCR4 via a novel WW domain-mediated interaction. Mol Biol Cell 20: 1324-1339. doi: 10.1091/mbc.E08-03-0308
![]() |
[44] |
Mines MA, Goodwin JS, Limbird LE, et al. (2009) Deubiquitination of CXCR4 by USP14 is critical for both CXCL12-induced CXCR4 degradation and chemotaxis but not ERK activation. J Biol Chem 284: 5742-5752. doi: 10.1074/jbc.M808507200
![]() |
[45] |
Busillo JM, Benovic JL (2007) Regulation of CXCR4 signaling. Biochim Biophys Acta 1768: 952-963. doi: 10.1016/j.bbamem.2006.11.002
![]() |
[46] | Cojoc M, Peitzsch C, Trautmann F, et al. (2013) Emerging targets in cancer management: role of the CXCL12/CXCR4 axis. Onco Targets Ther 6: 1347-1361. |
[47] |
Mellado M, Rodríguez-Frade JM, Mañes S, et al. (2001) Chemokine signaling and functional responses: the role of receptor dimerization and TK pathway activation. Annu Rev Immunol 19: 397-421. doi: 10.1146/annurev.immunol.19.1.397
![]() |
[48] | Teicher BA, Fricker SP (2010) CXCL12 (SDF-1)/CXCR4 Pathway in Cancer. Clin Cancer Res 16: 2927-2931. |
[49] |
Graham GJ, Locati M, Mantovani A, et al. (2012) The biochemistry and biology of the atypical chemokine receptors. Immunol Lett 145: 30-38. doi: 10.1016/j.imlet.2012.04.004
![]() |
[50] |
Freitas C, Desnoyer A, Meuris F, et al. (2014) The relevance of the chemokine receptor ACKR3/CXCR7 on CXCL12-mediated effects in cancers with a focus on virus-related cancers. Cytokine Growth Factor Rev 25: 307-316. doi: 10.1016/j.cytogfr.2014.04.006
![]() |
[51] | Veldkamp CT, Peterson FC, Pelzek AJ, et al. (2005) The monomer—dimer equilibrium of stromal cell-derived factor-1 ( CXCL 12 ) is altered by pH , phosphate , sulfate , and heparin. Protein Sci 1: 1071-1081. |
[52] | Ray P, Lewin SA, Mihalko LA, et al. (2012) Secreted CXCL12 (SDF-1) forms dimers under physiological conditions. Biochem J 442: 433-442. |
[53] | Wu B, Chien EY, Mol CD, et al. (2010) Structures of the CXCR4 Chemokine. Science 330: 1066-1071. |
[54] | Kramp BK, Sarabi A, Koenen RR, et al. (2011) Heterophilic chemokine receptor interactions in chemokine signaling and biology. Exp Cell Res 317: 655-663. |
[55] |
Thelen M, Thelen S (2008) CXCR7, CXCR4 and CXCL12: An eccentric trio? J Neuroimmunol 198: 9-13. doi: 10.1016/j.jneuroim.2008.04.020
![]() |
[56] |
Laguri C, Arenzana-Seisdedos F, Lortat-Jacob H (2008) Relationships between glycosaminoglycan and receptor binding sites in chemokines-the CXCL12 example. Carbohydr Res 343: 2018-2023. doi: 10.1016/j.carres.2008.01.047
![]() |
[57] |
Rueda P, Balabanian K, Lagane B, et al. (2008) The CXCL12γ Chemokine Displays Unprecedented Structural and Functional Properties that Make It a Paradigm of Chemoattractant Proteins. Gold JA, editor. PLoS One 3: e2543. doi: 10.1371/journal.pone.0002543
![]() |
[58] | Sadir R, Imberty A, Baleux F, et al. (2004) Heparan Sulfate/Heparin oligosaccharides protect stromal cell-derived factor-1 (SDF-1)/CXCL12 against proteolysis induced by CD26/dipeptidyl peptidase IV. J Biol Chem 279: 43854-43860. |
[59] | Brule S, Friand V, Sutton A, et al. (2009) Glycosaminoglycans and syndecan-4 are involved in SDF-1/CXCL12-mediated invasion of human epitheloid carcinoma HeLa cells. Biochim Biophys Acta-Gen Subj 1790: 1643-1650. |
[60] | Shirozu M, Nakano T, Inazawa J, et al. (1995) Struture and chromosomal localization of the human stromal cell-derived factor 1 (SDF1) gene. Genomics 28: 495-500. |
[61] | Winkler C, Modi W, Smith MW, et al. (1998) Genetic restriction of AIDS pathogenesis by an SDF-1 chemokine gene variant. ALIVE Study, Hemophilia Growth and Development Study (HGDS), Multicenter AIDS Cohort Study (MACS), Multicenter Hemophilia Cohort Study (MHCS), San Francisco City Cohort (SFCC). Science 279: 389-393. |
[62] |
Feng L, Nian S, Hao Y, et al. (2014) A Single Nucleotide Polymorphism in the Stromal Cell-Derived Factor 1 Gene Is Associated with Coronary Heart Disease in Chinese Patients. Int J Mol Sci 15: 11054-11063. doi: 10.3390/ijms150611054
![]() |
[63] | Reiche EMV, Ehara Watanabe MA, et al. (2006) The effect of stromal cell-derived factor 1 (SDF1/CXCL12) genetic polymorphism on HIV-1 disease progression. Int J Mol Med 18 785-793. |
[64] |
de Oliveira KB, Oda JMM, Voltarelli JC, et al. (2009) CXCL12 rs1801157 polymorphism in patients with breast cancer, hodgkin’s lymphoma, and non-hodgkin’s lymphoma. J Clin Lab Anal 23: 387-393. doi: 10.1002/jcla.20346
![]() |
[65] |
Cai C, Wang L-H, Dong Q, et al. (2013) Association of CXCL12 and CXCR4 gene polymorphisms with the susceptibility and prognosis of renal cell carcinoma. Tissue Antigens 82: 165-170. doi: 10.1111/tan.12170
![]() |
[66] | Schimanski C, Jordan M, Schlaegel F (2011) SNP rs1801157 significantly correlates with distant metastasis in CXCL12 expressing esophagogastri c cancer. Int J Oncol 39: 515-920. |
[67] | Razmkhah M, Doroudchi M, Ghayumi SMA, et al. ( 2005) Stromal cell-derived factor-1 (SDF-1) gene and susceptibility of Iranian patients with lung cancer. Lung Cancer 49: 311-315. |
[68] | Hassan S, Baccarelli A, Salvucci O, et al. (2008) Plasma stromal cell-derived factor-1: host derived marker predictive of distant metastasis in breast cancer. Clin Cancer Res 14: 446-454. |
[69] | Phillips RJ, Burdick MD, Lutz M, et al. ( 2003) The stromal derived factor–1/CXCL12–CXC chemokine receptor 4 biological axis in non–small cell lung cancer metastases. Am J Respir Crit Care Med 167: 1676-1686. |
[70] | Dehghani M, Kianpour S, Zangeneh A, et al. (2014) CXCL12 modulates prostate cancer cell adhesion by altering the levels or activities of β1-containing integrins. Int J Cell Biol 2014: 981750. |
[71] |
Beider K, Bitner H, Leiba M, et al. (2014) Multiple myeloma cells recruit tumor-supportive macrophages through the CXCR4/CXCL12 axis and promote their polarization toward the M2 phenotype. Oncotarget 5: 11283-11296. doi: 10.18632/oncotarget.2207
![]() |
[72] | Sei S, O’Neill DP, Stewart SK, et al. (2001) Increased level of stromal cell-derived factor-1 mRNA in peripheral blood mononuclear cells from children with AIDS-related lymphoma. Cancer Res 61: 5028-5037. |
[73] |
de Oliveira KB, Guembarovski RL, Oda JMM, et al. (2011) CXCL12 rs1801157 polymorphism and expression in peripheral blood from breast cancer patients. Cytokine 55: 260-265. doi: 10.1016/j.cyto.2011.04.017
![]() |
[74] |
Leiding JW, Holland SM (2012) Warts and all: human papillomavirus in primary immunodeficiencies. J Allergy Clin Immunol 130: 1030-1048. doi: 10.1016/j.jaci.2012.07.049
![]() |
[75] | Chow KYC, Brotin É, Ben Khalifa Y, et al. (2010) A Pivotal Role for CXCL12 Signaling in HPV-Mediated Transformation of Keratinocytes: Clues to Understanding HPV-Pathogenesis in WHIM Syndrome. Cell Host Microbe 8: 523-533. |
[76] |
Bachelerie F (2010) CXCL12/CXCR4-axis dysfunctions: Markers of the rare immunodeficiency disorder WHIM syndrome. Dis Markers 29: 189-198. doi: 10.1155/2010/475104
![]() |
[77] | Maciejewski-Duval A, Meuris F, Bignon A, et al. (2016) Altered chemotactic response to CXCL12 in patients carrying GATA2 mutations. J Leukoc Biol 99: 1065-1076. |
[78] |
Bignon A, Biajoux V, Bouchet-Delbos L, et al. (2011) CXCR4, a therapeutic target in rare immunodeficiencies? Med Sci (Paris) 27: 391-397. doi: 10.1051/medsci/2011274015
![]() |
[79] |
Maley SN, Schwartz SM, Johnson LG, et al. (2009) Genetic variation in CXCL12 and risk of cervical carcinoma: a population-based case-control study. Int J Immunogenet 36: 367-375. doi: 10.1111/j.1744-313X.2009.00877.x
![]() |
[80] |
Tee Y, Yang S-F, Wang P, et al. (2012) G801A Polymorphism of Human Stromal Cell–Derived Factor 1 Gene Raises No Susceptibility to Neoplastic Lesions of Uterine Cervix. Int J Gynecol Cancer 22: 1297-1302. doi: 10.1097/IGC.0b013e318265d334
![]() |
[81] | Roszak A, Misztal M, Sowińska A, et al. (2015) Stromal cell-derived factor-1 G801A polymorphism and the risk factors for cervical cancer. Mol Med Rep 11: 4633-4638. |
[82] |
Jaafar F, Righi E, Lindstrom V, et al. (2009) Correlation of CXCL12 expression and FoxP3+ cell infiltration with human papillomavirus infection and clinicopathological progression of cervical cancer. Am J Pathol 175: 1525-1535. doi: 10.2353/ajpath.2009.090295
![]() |
[83] | Huang Y, Zhang J, Cui Z-M, et al. (2013) Expression of the CXCL12/CXCR4 and CXCL16/CXCR6 axes in cervical intraepithelial neoplasia and cervical cancer. Chin J Cancer 32: 289-296. |
[84] |
Kurban S, Tursun M, Kurban G, et al. (2014) Role of CXCR7 and effects on CXCL12 in SiHa cells and upregulation in cervical squamous cell carcinomas in Uighur women. Asian Pac J Cancer Prev 15: 9211-9216 doi: 10.7314/APJCP.2014.15.21.9211
![]() |
[85] | Zanotta N, Tornesello ML, Annunziata C, et al. (2016) Candidate soluble immune mediators in young women with high-risk human papillomavirus infection: high expression of chemokines promoting angiogenesis and cell proliferation. PLoS One 11: e0151851. |
[86] |
Baker R, Dauner JG, Rodriguez AC, et al. (2011) Increased plasma levels of adipokines and inflammatory markers in older women with persistent HPV infection. Cytokine 53: 282-285. doi: 10.1016/j.cyto.2010.11.014
![]() |
[87] | Shen X, Wang S-H, Liang M, et al. (2008) The role and mechanism of CXCR4 and its ligand SDF-1 in the development of cervical cancer metastasis. Ai Zheng 27: 1044-1049. |
[88] | Yadav SS, Prasad SB, Das M, et al. (2014) Epigenetic silencing of CXCR4 promotes loss of cell adhesion in cervical cancer. Biomed Res Int 2014: 581403. |
[89] |
Donzella GA, Schols D, Lin SW, et al. (1998) AMD3100, a small molecule inhibitor of HIV-1 entry via the CXCR4 co-receptor. Nat Med 4: 72-77. doi: 10.1038/nm0198-072
![]() |
[90] | Burger JA, Stewart DJ (2009) CXCR4 chemokine receptor antagonists: perspectives in SCLC. Expert Opin Investig Drugs 18: 481-490. |
[91] | Sayyed SG, Hägele H, Kulkarni OP, et al. (2009) Podocytes produce homeostatic chemokine stromal cell-derived factor-1/CXCL12, which contributes to glomerulosclerosis, podocyte loss and albuminuria in a mouse model of type 2 diabetes. Diabetologia 52: 2445-2454. |
[92] |
Miao Z, Luker KE, Summers BC, et al. (2007) CXCR7 (RDC1) promotes breast and lung tumor growth in vivo and is expressed on tumor-associated vasculature. Proc Natl Acad Sci USA 104: 15735-15740. doi: 10.1073/pnas.0610444104
![]() |
[93] |
Wang J, Wang J, Sun Y, et al. (2005) Diverse signaling pathways through the SDF-1/CXCR4 chemokine axis in prostate cancer cell lines leads to altered patterns of cytokine secretion and angiogenesis. Cell Signal 17: 1578-1592. doi: 10.1016/j.cellsig.2005.03.022
![]() |
[94] |
Sowińska A, Jagodzinski PP (2007) RNA interference-mediated knockdown of DNMT1 and DNMT3B induces CXCL12 expression in MCF-7 breast cancer and AsPC1 pancreatic carcinoma cell lines. Cancer Lett 255: 153-159. doi: 10.1016/j.canlet.2007.04.004
![]() |
[95] |
Uchida D, Onoue T, Kuribayashi N, et al. (2011) Blockade of CXCR4 in oral squamous cell carcinoma inhibits lymph node metastases. Eur J Cancer 47: 452-459. doi: 10.1016/j.ejca.2010.09.028
![]() |
[96] |
Gupta S, Gupta S (2015) Role of human papillomavirus in oral squamous cell carcinoma and oral potentially malignant disorders: A review of the literature. Indian J Dent 6: 91-98. doi: 10.4103/0975-962X.155877
![]() |
[97] | Meuris F, Gaudin F, Aknin M-L, et al. (2016) Symptomatic Improvement in Human Papillomavirus-Induced Epithelial Neoplasia by Specific Targeting of the CXCR4 Chemokine Receptor. J Invest Dermatol 136: 473-480. |
[98] | von Knebel Doeberitz M (2001) New molecular tools for efficient screening of cervical cancer. Dis Markers 17: 123-128. |
[99] | Sun Y, Schneider A, Jung Y, et al. (2005) Skeletal localization and neutralization of the SDF-1(CXCL12)/CXCR4 axis blocks prostate cancer metastasis and growth in osseous sites in vivo. J Bone Miner Res 20: 318-329. |
[100] | Huang EH, Singh B, Cristofanilli M, et al. (2009) A CXCR4 Antagonist CTCE-9908 Inhibits Primary Tumor Growth and Metastasis of Breast Cancer. J Surg Res 155: 231-236. |
[101] |
Matsusue R, Kubo H, Hisamori S, et al. (2009) Hepatic stellate Cells cromote liver metastasis of colon cancer cells by the action of SDF-1/CXCR4 axis. Ann Surg Oncol 16: 2645-2653. doi: 10.1245/s10434-009-0599-x
![]() |
[102] | Rubin JB, Kung AL, Klein RS, et al. (2003) A small-molecule antagonist of CXCR4 inhibits intracranial growth of primary brain tumors. Proc Natl Acad Sci USA 100: 13513-13518. |
[103] |
Porvasnik S, Sakamoto N, Kusmartsev S, et al. (2009) Effects of CXCR4 antagonist CTCE-9908 on prostate tumor growth. Prostate 69: 1460-1469. doi: 10.1002/pros.21008
![]() |
[104] |
Epstein RJ (2004) The CXCL12-CXCR4 chemotactic pathway as a target of adjuvant breast cancer therapies. Nat Rev Cancer 4: 901-909. doi: 10.1038/nrc1473
![]() |
[105] |
Batchelor TT, Sorensen AG, di Tomaso E, et al. (2007) AZD2171, a pan-VEGF receptor tyrosine kinase inhibitor, normalizes tumor vasculature and alleviates edema in glioblastoma patients. Cancer Cell 11: 83-95. doi: 10.1016/j.ccr.2006.11.021
![]() |
[106] |
Kamoun WS, Ley CD, Farrar CT, et al. (2009) Edema control by cediranib, a vascular endothelial growth factor receptor-targeted kinase inhibitor, prolongs survival despite persistent brain tumor growth in mice. J Clin Oncol 27: 2542-2552. doi: 10.1200/JCO.2008.19.9356
![]() |
[107] | Shaked Y, Henke E, Roodhart JML, et al. (2008) Rapid chemotherapy-induced acute endothelial progenitor cell mobilization: implications for antiangiogenic drugs as chemosensitizing agents. Cancer Cell 14: 263-273. |
[108] |
Shaked Y, Tang T, Woloszynek J, et al. (2009) Contribution of granulocyte colony-stimulating factor to the acute mobilization of endothelial precursor cells by vascular disrupting agents. Cancer Res 69: 7524-7528. doi: 10.1158/0008-5472.CAN-09-0381
![]() |
[109] |
Kozin S V., Kamoun WS, Huang Y, et al. (2010) Recruitment of myeloid but not endothelial precursor cells facilitates tumor regrowth after local irradiation. Cancer Res 70: 5679-5685. doi: 10.1158/0008-5472.CAN-09-4446
![]() |
[110] |
Kioi M, Vogel H, Schultz G, et al. (2010) Inhibition of vasculogenesis, but not angiogenesis, prevents the recurrence of glioblastoma after irradiation in mice. J Clin Invest 120: 694-705. doi: 10.1172/JCI40283
![]() |
[111] | Hattermann K, Held-Feindt J, Lucius R, et al. (2010) The Chemokine Receptor CXCR7 Is Highly Expressed in Human Glioma Cells and Mediates Antiapoptotic Effects. Cancer Res 70: 3299-3308. |
[112] |
Zagzag D, Esencay M, Mendez O, et al. (2008) Hypoxia and vascular endothelial growth factor-induced stromal cell-derived factor-1α/CXCR4 expression in glioblastomas. Am J Pathol 173: 545-560. doi: 10.2353/ajpath.2008.071197
![]() |
[113] |
Batchelor TT, Duda DG, di Tomaso E, et al. (2010) Phase II study of cediranib, an oral pan-vascular endothelial growth factor receptor tyrosine kinase inhibitor, in patients with recurrent glioblastoma. J Clin Oncol 28: 2817-2823. doi: 10.1200/JCO.2009.26.3988
![]() |
1. | Michael Herrmann, Karsten Matthies, Nonlinear and nonlocal eigenvalue problems: variational existence, decay properties, approximation, and universal scaling limits, 2020, 33, 0951-7715, 4046, 10.1088/1361-6544/ab8350 | |
2. | Michael Herrmann, Karsten Matthies, A Uniqueness Result for a Simple Superlinear Eigenvalue Problem, 2021, 31, 0938-8974, 10.1007/s00332-021-09683-8 | |
3. | Anna Vainchtein, Solitary waves in FPU-type lattices, 2022, 434, 01672789, 133252, 10.1016/j.physd.2022.133252 | |
4. | Michael Herrmann, Katia Kleine, Korteweg–de Vries waves in peridynamical media, 2024, 152, 0022-2526, 376, 10.1111/sapm.12644 |