We consider a broad class of systems of nonlinear integro-differential equations posed on the real line that arise as Euler–Lagrange equations to energies involving nonlinear nonlocal interactions. Although these equations are not readily cast as dynamical systems, we develop a calculus that yields a natural Hamiltonian formalism. In particular, we formulate Noether’s theorem in this context, identify a degenerate symplectic structure, and derive Hamiltonian differential equations on finite-dimensional center manifolds when those exist. Our formalism yields new natural conserved quantities. For Euler–Lagrange equations arising as traveling-wave equations in gradient flows, we identify Lyapunov functions. We provide several applications to pattern-forming systems including neural field and phase separation problems.

We consider the unique recovery of a non-compactly supported and non-periodic perturbation of a Schrödinger operator in an unbounded cylindrical domain, also called waveguide, from boundary measurements. More precisely, we prove recovery of a general class of electric potentials from the partial Dirichlet-to-Neumann map, where the Dirichlet data is supported on slightly more than half of the boundary and the Neumann data is taken on the other half of the boundary. We apply this result in different contexts including recovery of some general class of non-compactly supported coefficients from measurements on a bounded subset and recovery of an electric potential, supported on an unbounded cylinder, of a Schrödinger operator in a slab.

We consider a class of optimal control problems for measure-valued nonlinear transport equations describing traffic flow problems on networks. The objective is to minimise/maximise macroscopic quantities, such as traffic volume or average speed, controlling few agents, e.g. smart traffic lights and automated cars. The measure theoretic approach allows to study in a same setting local and non-local drivers interactions and to consider the control variables as additional measures interacting with the drivers distribution. We also propose a gradient descent adjoint-based optimisation method, obtained by deriving first-order optimality conditions for the control problem, and we provide some numerical experiments in the case of smart traffic lights for a 2–1 junction.

Estimates for sample paths of fast–slow stochastic ordinary differential equations have become a key mathematical tool relevant for theory and applications. In particular, there have been breakthroughs by Berglund and Gentz to prove sharp exponential error estimates. In this paper, we take the first steps in order to generalise this theory to fast–slow stochastic partial differential equations. In a simplified setting with a natural decomposition into low- and high-frequency modes, we demonstrate that for a short-time period the probability for the corresponding sample path to leave a neighbourhood around the stable slow manifold of the system is exponentially small as well.

Critical transitions (or tipping points) are drastic sudden changes observed in many dynamical systems. Large classes of critical transitions are associated with systems, which drift slowly towards a bifurcation point. In the context of stochastic ordinary differential equations, there are results on growth of variance and autocorrelation before a transition, which can be used as possible warning signs in applications. A similar theory has recently been developed in the simplest setting for stochastic partial differential equations (SPDEs) for self-adjoint operators in the drift term. This setting leads to real discrete spectrum and growth of the covariance operator via a certain scaling law. In this paper, we develop this theory substantially further. We cover the cases of complex eigenvalues, degenerate eigenvalues as well as continuous spectrum. This provides a fairly comprehensive theory for most practical applications of warning signs for SPDE bifurcations.

This work is devoted to the study of an integro-differential system of equations modelling the genetic adaptation of a pathogen by taking into account both mutation and selection processes. First, we study the asymptotic behaviour of the system and prove that it eventually converges to a stationary state. Next, we more closely investigate the behaviour of the system in the presence of multiple EAs. Under suitable assumptions and based on a small mutation variance asymptotic, we describe the existence of a long transient regime during which the pathogen population remains far from its asymptotic behaviour and highly concentrated around some phenotypic value that is different from the one described by its asymptotic behaviour. In that setting, the time needed for the system to reach its large time configuration is very long and multiple evolutionary attractors may act as a barrier of evolution that can be very long to bypass.

This paper is concerned with the global dynamics and spreading speeds of a partially degenerate non-local dispersal system with monostable nonlinearity in periodic habitats. We first obtain the existence of the principal eigenvalue for a periodic eigenvalue problem with partially degenerate non-local dispersal. Then we study the coexistence and extinction dynamics. Finally, the existence and characterization of spreading speeds are considered. In particular, we show that the spreading speed is linearly determinate. Overall, we extend the existing results on global dynamics and spreading speeds for the degenerate reaction–diffusion system to the degenerate non-local dispersal case. The extension is non-trivial and meaningful.

This paper is concerned with the stability issue in determining absorption and diffusion coefficients in photoacoustic imaging. Assuming that the medium is layered and the acoustic wave speed is known, we derive global Hölder stability estimates of the photoacoustic inversion. These results show that the reconstruction is stable in the region close to the optical illumination source, and deteriorate exponentially far away. Several experimental pointed out that the resolution depth of the photoacoustic modality is about tens of millimeters. Our stability estimates confirm these observations and give a rigorous quantification of this depth resolution.

The scope of the paper is twofold. We show that for a large class of measurable vector fields in the sense of Weaver (i.e. derivations over the algebra of Lipschitz functions), called in the paper laminated, the notion of integral curves may be naturally defined and characterized (when appropriate) by an ordinary differential equation. We further show that for such vector fields the notion of a flow of the given positive Borel measure similar to the classical one generated by a smooth vector field (in a space with smooth structure) may be defined in a reasonable way, so that the measure ‘flows along’ the appropriately understood integral curves of the given vector field and the classical continuity equation is satisfied in the weak sense.

We construct travelling waves in the Burgers equation with the fractional Laplacian $(D^{2})^{\unicode[STIX]{x1D6FC}}$ , $\unicode[STIX]{x1D6FC}\in (1/2,1)$ . This is done by first constructing odd solutions $u_{\unicode[STIX]{x1D700}}$ of $uu^{\prime }=K_{\unicode[STIX]{x1D700}_{1}}\ast u-k_{\unicode[STIX]{x1D700}_{1}}u+\unicode[STIX]{x1D700}_{2}u^{\prime \prime }$ , $u(-\infty )=u_{c}>0$ , with $K_{\unicode[STIX]{x1D700}_{1}}\ast u-k_{\unicode[STIX]{x1D700}_{1}}u$ nonsingular, and then passing to the limit $\unicode[STIX]{x1D700}_{1},\unicode[STIX]{x1D700}_{2}\rightarrow 0$ , to give $K_{\unicode[STIX]{x1D700}_{1}}\ast u_{\unicode[STIX]{x1D700}}-k_{\unicode[STIX]{x1D700}_{1}}u_{\unicode[STIX]{x1D700}}\rightarrow (D^{2})^{\unicode[STIX]{x1D6FC}}u_{0}$ pointwise. The proof relies on operator splitting.

This article deals with numerical inversion for the initial distribution in the multi-term time-fractional diffusion equation using final observations. The inversion problem is of instability, but it is uniquely solvable based on the solution's expression for the forward problem and estimation to the multivariate Mittag-Leffler function. From view point of optimality, solving the inversion problem is transformed to minimizing a cost functional, and existence of a minimum is proved by the weakly lower semi-continuity of the functional. Furthermore, the homotopy regularization algorithm is introduced based on the minimization problem to perform numerical inversions, and the inversion solutions with noisy data give good approximations to the exact initial distribution demonstrating the efficiency of the inversion algorithm.

We consider the free linear Schrödinger equation on a torus 𝕋d , perturbed by a Hamiltonian nonlinearity, driven by a random force and subject to a linear damping:

Here u = u(t, x), x ∈ 𝕋d , 0 < ν ≪ 1, q∗ ∈ ℕ, f is a positive continuous function, ρ is a positive parameter and are standard independent complex Wiener processes. We are interested in limiting, as ν → 0, behaviour of distributions of solutions for this equation and of its stationary measure. Writing the equation in the slow time τ = νt, we prove that the limiting behaviour of them both is described by the effective equation

where the nonlinearity F(u) is made out of the resonant terms of the monomial |u| 2q∗ u.

Pinning effect on current-induced magnetic transverse domain wall dynamics in nanostrip is studied for its potential application to new magnetic memory devices. In this study, we carry out a series of calculations by solving generalized Landau-Lifshitz equation involving a current spin transfer torque in one and two dimensional models. The critical current for the transverse wall depinning in nanostrip depends on the size of artificial rectangular defects on the edges of nanostrip. We show that there is intrinsic pinning potential for a defect such that the transverse wall oscillates damply around the pinning site with an intrinsic frequency if the applied current is below critical value. The amplification of the transverse wall oscillation for both displacement and maximum value of m3 is significant by applying AC current and current pulses with appropriate frequency. We show that for given pinning potential, the oscillation amplitude as a function of the frequency of the AC current behaves like a Gaussian distribution in our numerical study, which is helpful to reduce strength of current to drive the transverse wall motion.

We construct countably infinitely many non-radial singular solutions of the problem

of the form

where v(σ) depends only on σ ∈𝕊N−1. To this end we construct countably infinitely many solutions of

using ordinary differential equation techniques.

The fractional derivatives include nonlocal information and thus their calculation requires huge storage and computational cost for long time simulations. We present an efficient and high-order accurate numerical formula to speed up the evaluation of the Caputo fractional derivative based on the L2-1σ formula proposed in [A. Alikhanov, J. Comput. Phys., 280 (2015), pp. 424-438], and employing the sum-of-exponentials approximation to the kernel function appeared in the Caputo fractional derivative. Both theoretically and numerically, we prove that while applied to solving time fractional diffusion equations, our scheme not only has unconditional stability and high accuracy but also reduces the storage and computational cost.

Due to the increasing application of fractional calculus in engineering and biomedical processes, we analyze a new method for the numerical simulation of a large class of coupled systems of fractional-order partial differential equations. In this paper, we study shifted Jacobi polynomials in the case of two variables and develop some new operational matrices of fractional-order integrations as well as fractional-order differentiations. By the use of these operational matrices, we present a new and easy method for solving a generalized class of coupled systems of fractional-order partial differential equations subject to some initial conditions. We convert the system under consideration to a system of easily solvable algebraic equation without discretizing the system, and obtain a highly accurate solution. Also, the proposed method is compared with some other well-known differential transform methods. The proposed method is computer oriented. We use MatLab to perform the necessary calculation. The next two parts will appear soon.

This article is intended to fill in the blank of the numerical schemes with second-order convergence accuracy in time for nonlinear Stokes’ first problem for a heated generalized second grade fluid with fractional derivative. A linearized difference scheme is proposed. The time fractional-order derivative is discretized by second-order shifted and weighted Gr¨unwald-Letnikov difference operator. The convergence accuracy in space is improved by performing the average operator. The presented numerical method is unconditionally stable with the global convergence order of in maximum norm, where τ and h are the step sizes in time and space, respectively. Finally, numerical examples are carried out to verify the theoretical results, showing that our scheme is efficient indeed.

We consider a possibly anisotropic integrodifferential semilinear equation, driven by a non-decreasing nonlinearity. We prove that if the solution grows less than the order of the operator at infinity, then it must be affine (possibly constant).

This paper is concerned with the invisibility cloaking in acoustic wave scattering from a new perspective. We are especially interested in achieving the invisibility cloaking by completely regular and isotropic mediums. It is shown that an interior transmission eigenvalue problem arises in our study, which is the one considered theoretically in Cakoni et al. (Transmission eigenvalues for inhomogeneous media containing obstacles, Inverse Problems and Imaging, 6 (2012), 373–398). Based on such an observation, we propose a cloaking scheme that takes a three-layer structure including a cloaked region, a lossy layer and a cloaking shell. The target medium in the cloaked region can be arbitrary but regular, whereas the mediums in the lossy layer and the cloaking shell are both regular and isotropic. We establish that if a certain non-transparency condition is satisfied, then there exists an infinite set of incident waves such that the cloaking device is nearly invisible under the corresponding wave interrogation. The set of waves is generated from the Herglotz approximation of the associated interior transmission eigenfunctions. We provide both theoretical and numerical justifications.

We consider the unique solvability of an inverse-source problem with integral transmitting condition for a time-fractional mixed type equation in rectangular domain where the unknown source term depends only on the space variable. The solution is based on a series expansion using a bi-orthogonal basis in space, corresponding to a non-self-adjoint boundary value problem. Under certain regularity conditions on the given data, we prove the uniqueness and existence of the solution for the given problem. The influence of the transmitting condition on the solvability of the problem is also demonstrated. Two different transmitting conditions are considered — viz. a full integral form and a special case. In order to simplify the bulky expressions appearing in the proof of our main result, we establish a new property of the recently introduced Mittag-Leffler type function in two variables.