Quadratic Interaction Functional for General Systems of Conservation Laws

For the Glimm scheme approximation \({u_\varepsilon}\) to the solution of the system of conservation laws in one space dimension

$$u_t + f(u)_x = 0, \qquad u(0, x) = u_0(x) \in \mathbb{R}^n,$$

with initial data u ₀ with small total variation, we prove a quadratic (w.r.t. Tot. Var. (u ₀)) interaction estimate, which has been used in the literature for stability and convergence results. No assumptions on the structure of the flux f are made (apart from smoothness), and this estimate is the natural extension of the Glimm type interaction estimate for genuinely nonlinear systems.

More precisely, we obtain the following results:

• a new analysis of the interaction estimates of simple waves;
• a Lagrangian representation of the derivative of the solution, i.e., a map \({\mathtt{x}(t, w)}\) which follows the trajectory of each wave w from its creation to its cancellation;
• the introduction of the characteristic interval and partition for couples of waves, representing the common history of the two waves;
• a new functional \({\mathfrak{Q}}\) controlling the variation in speed of the waves w.r.t. time.

This last functional is the natural extension of the Glimm functional for genuinely nonlinear systems.

The main result is that the distribution \({D_{t} \hat \sigma_k(t,w)}\) is a measure with total mass \({\leq \mathcal{O}(1) {\rm Tot. Var.} (u_0)^2}\) , where \({\hat{\sigma}_k(t, w)}\) is the speed given to the wave w by the Riemann problem at the grid point \({(i\varepsilon, \mathtt{x}(i\varepsilon, w)), t \in [i\varepsilon, (i + 1)\varepsilon)}\).     

Spectral and Scattering Theory of Space-cutoff Charged $${P(\varphi)_{2}}$$ Models

We consider in this paper space-cutoff charged \({P(\varphi)_{2}}\) models arising from the quantization of the non-linear charged Klein–Gordon equation:

$(\partial_{t}+\i V(x))^{2}\phi(t, x)+ (-\Delta_{x}+ m^{2})\phi(t,x)+ g(x)\partial_{\overline{z}}P(\phi(t,x), \overline{\phi}(t,x))=0,$

where V(x) is an electrostatic potential, g(x) ≥ 0 a space-cutoff, and \({P(\lambda, \overline{\lambda})}\) a real bounded below polynomial. We discuss various ways to quantize this equation, starting from different CCR representations. After describing the construction of the interacting Hamiltonian H we study its spectral and scattering theory. We describe the essential spectrum of H, prove the existence of asymptotic fields and of wave operators, and finally prove the asymptotic completeness of wave operators. These results are similar to the case when V = 0.

     

Domain and Range of the Modified Wave Operator for Schrödinger Equations with a Critical Nonlinearity

We study the final problem for the nonlinear Schrödinger equation

$i{\partial }_{t}u+\frac{1}{2}\Delta u=\lambda|u|^{\frac{2}{n}}u,\quad (t,x)\in {\mathbf{R}}\times \mathbf{R}^{n},$

where\(\lambda \in{\bf R},n=1,2,3\). If the final data\(u_{+}\in {\bf H}^{0,\alpha }=\left\{ \phi \in {\bf L}^{2}:\left( 1+\left\vert x\right\vert \right) ^{\alpha }\phi \in {\bf L}^{2}\right\} \) with\(\frac{ n}{2} < \alpha < \min \left( n,2,1+\frac{2}{n}\right) \) and the norm\(\Vert \widehat{u_{+}}\Vert _{{\bf L}^{\infty }}\) is sufficiently small, then we prove the existence of the wave operator in L ². We also construct the modified scattering operator from H ^0,α to H ^0,δ with\(\frac{n}{2} < \delta < \alpha\).

     

Normal Forms for Semilinear Quantum Harmonic Oscillators

We consider the semilinear harmonic oscillator

$i\psi_t=(-\Delta +{|x|}^{2} +M)\psi +\partial_2 g(\psi,\bar \psi), \quad x\in \mathbb{R}^{d},\, t\in \mathbb{R},$

where M is a Hermite multiplier and g a smooth function globally of order 3 at least. We prove that such a Hamiltonian equation admits, in a neighborhood of the origin, a Birkhoff normal form at any order and that, under generic conditions on M related to the non resonance of the linear part, this normal form is integrable when d = 1 and gives rise to simple (in particular bounded) dynamics when d ≥ 2. As a consequence we prove the almost global existence for solutions of the above equation with small Cauchy data. Furthermore we control the high Sobolev norms of these solutions.

     

Krein formula with compensated singularities for the ND-mapping and the generalized Kirchhoff condition at the Neumann Schrödinger junction

The Neumann Schrödinger operator \(\mathcal{L}\) is considered on a thin 2D star-shaped junction, composed of a vertex domain Ω_int and a few semi-infinite straight leads ω _m , m = 1, 2, ..., M, of width δ, δ ? diam Ω_int, attached to Ω_int at Γ ? ?Ω_int. The potential of the Schrödinger operator l ^ω on the leads vanishes, hence there are only a finite number of eigenvalues of the Neumann Schrödinger operator L _int on Ω_int embedded into the open spectral branches of l ^ω with oscillating solutions χ _±(x, p) = \(e^{ \pm iK_ + x} e_m \) of l ^ω χ _± = p ² χ _±. The exponent of the open channels in the wires is

$K_ + (\lambda ) = p\sum\limits_{m = 1}^M {e^m } \rangle \langle e^m = \sqrt \lambda P_ + $

, with constant e _m , on a relatively small essential spectral interval Δ ? [0, π ² δ ^?2). The scattering matrix of the junction is represented on Δ in terms of the ND mapping

$\mathcal{N} = \frac{{\partial P_ + \Psi }}{{\partial x}}(0,\lambda )\left| {_\Gamma \to P_ + \Psi _ + (0,\lambda )} \right|_\Gamma $

as

$S(\lambda ) = (ip\mathcal{N} + I_ + )^{ - 1} (ip\mathcal{N} - I_ + ), I_ + = \sum\limits_{m = 1}^M {e^m } \rangle \langle e^m = P_ + $

. We derive an approximate formula for \(\mathcal{N}\) in terms of the Neumann-to-Dirichlet mapping \(\mathcal{N}_{\operatorname{int} } \) of L _int and the exponent K _? of the closed channels of l ^ω . If there is only one simple eigenvalue λ ₀ ∈ Δ, L _intφ0 = λ _0φ0 then, for a thin junction, \(\mathcal{N} \approx |\vec \phi _0 |^2 P_0 (\lambda _0 - \lambda )^{ - 1} \) with

$\vec \phi _0 = P_ + \phi _0 = (\delta ^{ - 1} \int_{\Gamma _1 } {\phi _0 (\gamma )} d\gamma ,\delta ^{ - 1} \int_{\Gamma _2 } {\phi _0 (\gamma )} d\gamma , \ldots \delta ^{ - 1} \int_{\Gamma _M } {\phi _0 (\gamma )} d\gamma )$

and \(P_0 = \vec \phi _0 \rangle |\vec \phi _0 |^{ - 2} \langle \vec \phi _0 \),

$S(\lambda ) \approx \frac{{ip|\vec \phi _0 |^2 P_0 (\lambda _0 - \lambda )^{ - 1} - I_ + }}{{ip|\vec \phi _0 |^2 P_0 (\lambda _0 - \lambda )^{ - 1} + I_ + }} = :S_{appr} (\lambda )$

. The related boundary condition for the components P ₊Ψ(0) and P ₊Ψ′(0) of the scattering Ansatz in the open channel \(P_ + \Psi (0) = (\bar \Psi _1 ,\bar \Psi _2 , \ldots ,\bar \Psi _M ), P_ + \Psi '(0) = (\bar \Psi '_1 , \bar \Psi '_2 , \ldots , \bar \Psi '_M )\) includes the weighted continuity (1) of the scattering Ansatz Ψ at the vertex and the weighted balance of the currents (2), where

$\frac{{\bar \Psi _m }}{{\bar \phi _0^m }} = \frac{{\delta \sum\nolimits_{t = 1}^M { \bar \Psi _t \bar \phi _0^t } }}{{|\vec \phi _0 |^2 }} = \frac{{\bar \Psi _r }}{{\bar \phi _0^r }} = :\bar \Psi (0)/\bar \phi (0), 1 \leqslant m,r \leqslant M$

(1)

,

$\sum\limits_{m = 1}^M {\bar \Psi '_m } \bar \phi _0^m + \delta ^{ - 1} (\lambda - \lambda _0 )\bar \Psi /\bar \phi (0) = 0$

(1)

. Conditions (1) and (2) constitute the generalized Kirchhoff boundary condition at the vertex for the Schrödinger operator on a thin junction and remain valid for the corresponding 1D model. We compare this with the previous result by Kuchment and Zeng obtained by the variational technique for the Neumann Laplacian on a shrinking quantum network.

     

Super-Exponential Decay of Diffraction Managed Solitons

This is the second part of a series of papers where we develop rigorous decay estimates for breather solutions of an averaged version of the non-linear Schrödinger equation. In this part we study the diffraction managed discrete non-linear Schrödinger equation, an equation which describes coupled waveguide arrays with periodic diffraction management geometries. We show that, for vanishing average diffraction, any solution \({f \in l^2(\mathbb{Z})}\) of the non-linear and non-local diffraction management equation decays super-exponentially. More precisely, we have the bound

$\limsup_{|x|\to\infty}\big((|x|+1)\ln(|x|+1)\big)^{-1}\, \ln|f(x)|\le -\frac{1}{4}$

for any diffraction management soliton.

     

A Centre-Stable Manifold for the Focussing Cubic NLS in $${\mathbb{R}}^{1+3}$$

Consider the focussing cubic nonlinear Schrödinger equation in \({\mathbb{R}}^3\) :

$i\psi_t+\Delta\psi = -|\psi|^2 \psi. \quad (0.1) $

It admits special solutions of the form e ^itα ?, where \(\phi \in {\mathcal{S}}({\mathbb{R}}^3)\) is a positive (? > 0) solution of

$-\Delta \phi + \alpha\phi = \phi^3. \quad (0.2)$

The space of all such solutions, together with those obtained from them by rescaling and applying phase and Galilean coordinate changes, called standing waves, is the 8-dimensional manifold that consists of functions of the form \(e^{i(v \cdot + \Gamma)} \phi(\cdot - y, \alpha)\) . We prove that any solution starting sufficiently close to a standing wave in the \(\Sigma = W^{1, 2}({\mathbb{R}}^3) \cap |x|^{-1}L^2({\mathbb{R}}^3)\) norm and situated on a certain codimension-one local Lipschitz manifold exists globally in time and converges to a point on the manifold of standing waves. Furthermore, we show that \({\mathcal{N}}\) is invariant under the Hamiltonian flow, locally in time, and is a centre-stable manifold in the sense of Bates, Jones [BatJon]. The proof is based on the modulation method introduced by Soffer and Weinstein for the L ²-subcritical case and adapted by Schlag to the L ²-supercritical case. An important part of the proof is the Keel-Tao endpoint Strichartz estimate in \({\mathbb{R}}^3\) for the nonselfadjoint Schrödinger operator obtained by linearizing (0.1) around a standing wave solution. All results in this paper depend on the standard spectral assumption that the Hamiltonian

$\mathcal H = \left ( \begin{array}{cc}\Delta + 2\phi(\cdot, \alpha)^2 - \alpha &;\quad \phi(\cdot, \alpha)^2 \\ -\phi(\cdot, \alpha)^2 &;\quad -\Delta - 2 \phi(\cdot, \alpha)^2 + \alpha \end{array}\right ) \quad (0.3)$

has no embedded eigenvalues in the interior of its essential spectrum \((-\infty, -\alpha) \cup (\alpha, \infty)\) .

     

Noncommutative Lévy Processes for Generalized (Particularly Anyon) Statistics

Let \({T=\mathbb R^d}\) . Let a function \({QT^2\to\mathbb C}\) satisfy \({Q(s,t)=\overline{Q(t,s)}}\) and \({|Q(s,t)|=1}\). A generalized statistics is described by creation operators \({\partial_t^\dagger}\) and annihilation operators ? _t , \({t\in T}\), which satisfy the Q-commutation relations: \({\partial_s\partial^\dagger_t = Q(s, t)\partial^\dagger_t\partial_s+\delta(s, t)}\) , \({\partial_s\partial_t = Q(t, s)\partial_t\partial_s}\), \({\partial^\dagger_s\partial^\dagger_t = Q(t, s)\partial^\dagger_t\partial^\dagger_s}\). From the point of view of physics, the most important case of a generalized statistics is the anyon statistics, for which Q(s, t) is equal to q if s < t, and to \({\bar q}\) if s > t. Here \({q\in\mathbb C}\) , |q| = 1. We start the paper with a detailed discussion of a Q-Fock space and operators \({(\partial_t^\dagger,\partial_t)_{t\in T}}\) in it, which satisfy the Q-commutation relations. Next, we consider a noncommutative stochastic process (white noise) \({\omega(t)=\partial_t^\dagger+\partial_t+\lambda\partial_t^\dagger\partial_t}\) , \({t\in T}\) . Here \({\lambda\in\mathbb R}\) is a fixed parameter. The case λ = 0 corresponds to a Q-analog of Brownian motion, while λ ≠ 0 corresponds to a (centered) Q-Poisson process. We study Q-Hermite (Q-Charlier respectively) polynomials of infinitely many noncommutatative variables \({(\omega(t))_{t\in T}}\) . The main aim of the paper is to explain the notion of independence for a generalized statistics, and to derive corresponding Lévy processes. To this end, we recursively define Q-cumulants of a field \({(\xi(t))_{t\in T}}\). This allows us to define a Q-Lévy process as a field \({(\xi(t))_{t\in T}}\) whose values at different points of T are Q-independent and which possesses a stationarity of increments (in a certain sense). We present an explicit construction of a Q-Lévy process, and derive a Nualart–Schoutens-type chaotic decomposition for such a process.     

A global attracting set for the Kuramoto-Sivashinsky equation

New bounds are given for the L²-norm of the solution of the Kuramoto-Sivashinsky equation $$\partial _t U(x,t) = - (\partial _x^2 + \partial _x^4 )U(x,t) - U(x,t)\partial _x U(x,t)$$ , for initial data which are periodic with periodL. There is no requirement on the antisymmetry of the initial data. The result is $$\mathop {\lim \sup }\limits_{t \to \infty } \left\| {U( \cdot ,t)} \right\|_2 \leqslant const. L^{8/5} $$ .     

The Scaling Limit of the KPZ Equation in Space Dimension 3 and Higher

We study in the present article the Kardar–Parisi–Zhang (KPZ) equation

$$\begin{aligned} \partial _t h(t,x)=\nu \Delta h(t,x)+\lambda |\nabla h(t,x)|^2 +\sqrt{D}\, \eta (t,x), \qquad (t,x)\in \mathbb {R}_+\times \mathbb {R}^d \end{aligned}$$

in \(d\ge 3\) dimensions in the perturbative regime, i.e. for \(\lambda >0\) small enough and a smooth, bounded, integrable initial condition \(h_0=h(t=0,\cdot )\). The forcing term \(\eta \) in the right-hand side is a regularized space-time white noise. The exponential of h—its so-called Cole-Hopf transform—is known to satisfy a linear PDE with multiplicative noise. We prove a large-scale diffusive limit for the solution, in particular a time-integrated heat-kernel behavior for the covariance in a parabolic scaling. The proof is based on a rigorous implementation of K. Wilson’s renormalization group scheme. A double cluster/momentum-decoupling expansion allows for perturbative estimates of the bare resolvent of the Cole-Hopf linear PDE in the small-field region where the noise is not too large, following the broad lines of Iagolnitzer and Magnen (Commun Math Phys 162(1):85–121, 1994). Standard large deviation estimates for \(\eta \) make it possible to extend the above estimates to the large-field region. Finally, we show, by resumming all the by-products of the expansion, that the solution h may be written in the large-scale limit (after a suitable Galilei transformation) as a small perturbation of the solution of the underlying linear Edwards–Wilkinson model (\(\lambda =0\)) with renormalized coefficients \(\nu _{eff}=\nu +O(\lambda ^2),D_{eff}=D+O(\lambda ^2)\).

     

Parabolic Anderson Model in a Dynamic Random Environment: Random Conductances

The parabolic Anderson model is defined as the partial differential equation ? u(x, t)/? t = κ Δ u(x, t) + ξ(x, t)u(x, t), x ∈ ? ^d , t ≥ 0, where κ ∈ [0, ∞) is the diffusion constant, Δ is the discrete Laplacian, and ξ is a dynamic random environment that drives the equation. The initial condition u(x, 0) = u ₀(x), x ∈ ? ^d , is typically taken to be non-negative and bounded. The solution of the parabolic Anderson equation describes the evolution of a field of particles performing independent simple random walks with binary branching: particles jump at rate 2d κ, split into two at rate ξ ∨ 0, and die at rate (?ξ) ∨ 0. In earlier work we looked at the Lyapunov exponents

$$ \lambda _{p}(\kappa ) = \lim\limits _{t\to \infty } \frac {1}{t} \log \mathbb {E} ([u(0,t)]^{p})^{1/p}, \quad p \in \mathbb{N} , \qquad \lambda _{0}(\kappa ) = \lim\limits _{t\to \infty } \frac {1}{t}\log u(0,t). $$

For the former we derived quantitative results on the κ-dependence for four choices of ξ : space-time white noise, independent simple random walks, the exclusion process and the voter model. For the latter we obtained qualitative results under certain space-time mixing conditions on ξ. In the present paper we investigate what happens when κΔ is replaced by Δ^??, where ?? = {??(x, y) : x, y ∈ ? _d , x ～ y} is a collection of random conductances between neighbouring sites replacing the constant conductances κ in the homogeneous model. We show that the associated annealed Lyapunov exponents λ _p (??), p ∈ ?, are given by the formula

$$ \lambda _{p}(\mathcal{K} ) = \text{sup} \{\lambda _{p}(\kappa ) : \, \kappa \in \text{Supp} (\mathcal{K} )\}, $$

where, for a fixed realisation of ??, Supp(??) is the set of values taken by the ??-field. We also show that for the associated quenched Lyapunov exponent λ ₀(??) this formula only provides a lower bound, and we conjecture that an upper bound holds when Supp(??) is replaced by its convex hull. Our proof is valid for three classes of reversible ξ, and for all ?? satisfying a certain clustering property, namely, there are arbitrarily large balls where ?? is almost constant and close to any value in Supp(??). What our result says is that the annealed Lyapunov exponents are controlled by those pockets of ?? where the conductances are close to the value that maximises the growth in the homogeneous setting. In contrast our conjecture says that the quenched Lyapunov exponent is controlled by a mixture of pockets of ?? where the conductances are nearly constant. Our proof is based on variational representations and confinement arguments.

     

Random Walk with Barycentric Self-interaction

We study the asymptotic behaviour of a d-dimensional self-interacting random walk (X _n ) _n∈? (?:={1,2,3,…}) which is repelled or attracted by the centre of mass \(G_{n} = n^{-1} \sum_{i=1}^{n} X_{i}\) of its previous trajectory. The walk’s trajectory (X ₁,…,X _n ) models a random polymer chain in either poor or good solvent. In addition to some natural regularity conditions, we assume that the walk has one-step mean drift

$\mathbb{E}[X_{n+1} - X_n \mid X_n - G_n = \mathbf{x}] \approx\rho\|\mathbf{x}\|^{-\beta}\hat{ \mathbf{x}}$

for ρ∈? and β≥0. When β<1 and ρ>0, we show that X _n is transient with a limiting (random) direction and satisfies a super-diffusive law of large numbers: n ^?1/(1+β) X _n converges almost surely to some random vector. When β∈(0,1) there is sub-ballistic rate of escape. When β≥0 and ρ∈? we give almost-sure bounds on the norms ‖X _n ‖, which in the context of the polymer model reveal extended and collapsed phases.

Analysis of the random walk, and in particular of X _n ?G _n , leads to the study of real-valued time-inhomogeneous non-Markov processes (Z _n ) _n∈? on [0,∞) with mean drifts of the form

$ \mathbb{E}[ Z_{n+1} - Z_n \mid Z_n = x ] \approx\rho x^{-\beta} - \frac {x}{n},$

(0.1)

where β≥0 and ρ∈?. The study of such processes is a time-dependent variation on a classical problem of Lamperti; moreover, they arise naturally in the context of the distance of simple random walk on ? ^d from its centre of mass, for which we also give an apparently new result. We give a recurrence classification and asymptotic theory for processes Z _n satisfying (0.1), which enables us to deduce the complete recurrence classification (for any β≥0) of X _n ?G _n for our self-interacting walk.

     

Periodic and quasi-periodic solutions of nonlinear wave equations via KAM theory

In this paper the nonlinear wave equation $$u_u - u_{xx} + v(x)u(x,t) + \varepsilon u^3 (x,t) = 0$$ is studied. It is shown that for a large class of potentials,v(x), one can use KAM methods to construct periodic and quasi-periodic solutions (in time) for this equation.     

Front solutions for the Ginzburg-Landau equation

We prove the existence of front solutions for the Ginzburg-Landau equation $$\partial _t u(x,t) = \partial _x^2 u(x,t) + (1 - |u(x,t)|^2 )u(x,t)$$ , interpolating between two stationary solutions of the form \(u(x) = \sqrt {1 - q^2 } e^{iqx}\) with different values ofq atx=±∞. Such fronts are shown to exist when at least one of theq is in the Eckhaus-unstable domain.     

Entanglement Involved in Time Evolution of Two-Mode Squeezed State in Single-Mode Diffusion Channel

We derive the evolution law of an initial two-mode squeezed vacuum state \( \text {sech}^{2}\lambda e^{a^{\dag }b^{\dagger }\tanh \lambda }\left \vert 00\right \rangle \left \langle 00\right \vert e^{ab\tanh \lambda }\) (a pure state) passing through an a-mode diffusion channel described by the master equation

$$\frac{d\rho \left( t\right) }{dt}=-\kappa \left[ a^{\dagger}a\rho \left( t\right) -a^{\dagger}\rho \left( t\right) a-a\rho \left( t\right) a^{\dagger}+\rho \left( t\right) aa^{\dagger}\right] , $$

since the two-mode squeezed state is simultaneously an entangled state, the final state which emerges from this channel is a two-mode mixed state. Performing partial trace over the b-mode of ρ(t) yields a new chaotic field, \(\rho _{a}\left (t\right ) =\frac {\text {sech}^{2}\lambda }{1+\kappa t \text {sech}^{2}\lambda }:\exp \left [ \frac {- \text {sech}^{2}\lambda }{1+\kappa t\text {sech}^{2}\lambda }a^{\dagger }a \right ] :,\) which exhibits higher temperature and more photon numbers, showing the diffusion effect. Besides, measuring a-mode of ρ(t) to find n photons will result in the collapse of the two-mode system to a new Laguerre polynomial-weighted chaotic state in b-mode, which also exhibits entanglement.

     

On the Exact Evaluation of the Face-Centred Cubic Lattice Green Function

The mathematical properties of the lattice Green function

Open image in new window

are investigated, where w=w ₁+iw ₂ lies in a complex plane which is cut from w=?1 to w=3, and {? ₁,? ₂,? ₃} is a set of integers with ? ₁+? ₂+? ₃ equal to an even integer. In particular, it is proved that G(2n,0,0;w), where n=0,1,2,…, is a solution of a fourth-order linear differential equation of the Fuchsian type with four regular singular points at w=?1,0,3 and ∞. It is also shown that G(2n,0,0;w) satisfies a five-term recurrence relation with respect to the integer variable n. The limiting function

$G^{-}(2n,0,0;w_1)\equiv\lim_{\epsilon\rightarrow0+}G(2n,0,0;w_1-\mathrm{i}\epsilon) =G_{\mathrm{R}}(2n,0,0;w_1)+\mathrm{i}G_{\mathrm {I}}(2n,0,0;w_1) ,\nonumber $

where w ₁∈(?1,3), is evaluated exactly in terms of ₂ F ₁ hypergeometric functions and the special cases G ^?(2n,0,0;0), G ^?(2n,0,0;1) and G(2n,0,0;3) are analysed using singular value theory. More generally, it is demonstrated that G(? ₁,? ₂,? ₃;w) can be written in the form

Open image in new window

where Open image in new window are rational functions of the variable ξ, K(k _?) and E(k _?) are complete elliptic integrals of the first and second kind, respectively, with

$k_{-}^2\equiv k_{-}^2(w)={1\over2}- {2\over w} \biggl(1+{1\over w} \biggr)^{-{3\over2}}- {1\over2} \biggl(1-{1\over w} \biggr ) \biggl(1+{1\over w} \biggr)^{-{3\over2}} \biggl(1-{3\over w} \biggr)^{1\over2}\nonumber $

and the parameter ξ is defined as

$\xi\equiv\xi(w)= \biggl(1+\sqrt{1-{3\over w}} \,\biggr)^{-1} \biggl(-1+\sqrt{1+{1\over w}} \,\biggr) .\nonumber $

This result is valid for all values of w which lie in the cut plane. The asymptotic behaviour of G ^?(2n,0,0;w ₁) and G(2n,0,0;w ₁) as n→∞ is also determined. In the final section of the paper a new ₂ F ₁ product form for the anisotropic face-centred cubic lattice Green function is given.

     

Meixner Class of Non-Commutative Generalized Stochastic Processes with Freely Independent Values I. A Characterization

Let T be an underlying space with a non-atomic measure σ on it (e.g. \({T=\mathbb R^d}\) and σ is the Lebesgue measure). We introduce and study a class of non-commutative generalized stochastic processes, indexed by points of T, with freely independent values. Such a process (field), ω = ω(t), \({t\in T}\) , is given a rigorous meaning through smearing out with test functions on T, with \({\int_T \sigma(dt)f(t)\omega(t)}\) being a (bounded) linear operator in a full Fock space. We define a set CP of all continuous polynomials of ω, and then define a non-commutative L ²-space L ²(τ) by taking the closure of CP in the norm \({\|P\|_{L^2(\tau)}:=\|P\Omega\|}\) , where Ω is the vacuum in the Fock space. Through procedure of orthogonalization of polynomials, we construct a unitary isomorphism between L ²(τ) and a (Fock-space-type) Hilbert space \({\mathbb F=\mathbb R\oplus\bigoplus_{n=1}^\infty L^2(T^n,\gamma_n)}\) , with explicitly given measures γ _n . We identify the Meixner class as those processes for which the procedure of orthogonalization leaves the set CP invariant. (Note that, in the general case, the projection of a continuous monomial of order n onto the n ^th chaos need not remain a continuous polynomial.) Each element of the Meixner class is characterized by two continuous functions λ and η ≥ 0 on T, such that, in the \({\mathbb F}\) space, ω has representation \({\omega(t)=\partial_t^\dagger+\lambda(t)\partial_t^\dagger\partial_t+\partial_t+\eta(t)\partial_t^\dagger\partial^2_t}\) , where \({\partial_t^\dagger}\) and ? _t are the usual creation and annihilation operators at point t.     

A Lax pair structure for the half-wave maps equation

We consider the half-wave maps equation

$$\begin{aligned} \partial _t \vec {S} = \vec {S} \wedge |\nabla | \vec {S}, \end{aligned}$$

where \(\vec {S}= \vec {S}(t,x)\) takes values on the two-dimensional unit sphere \(\mathbb {S}^2\) and \(x \in \mathbb {R}\) (real line case) or \(x \in \mathbb {T}\) (periodic case). This an energy-critical Hamiltonian evolution equation recently introduced in Lenzmann and Schikorra (2017, arXiv:1702.05995v2), Zhou and Stone (Phys Lett A 379:2817–2825, 2015) which formally arises as an effective evolution equation in the classical and continuum limit of Haldane–Shastry quantum spin chains. We prove that the half-wave maps equation admits a Lax pair and we discuss some analytic consequences of this finding. As a variant of our arguments, we also obtain a Lax pair for the half-wave maps equation with target \(\mathbb {H}^2\) (hyperbolic plane).

     

Higher recurrences for Apostol-Bernoulli-Euler numbers

We present explicit formulas for sums of products of Apostol-Bernoulli and Apostol-Euler numbers of the form

$\sum\limits_{_{m_1 , \cdots ,m_N \geqslant n}^{m_1 + \cdots + m_N = n} } {\left( {_{m_1 , \cdots m_N }^n } \right)B_{m_1 } (q) \cdots B_{m_N } (q),} \sum\limits_{_{m_1 , \cdots ,m_N \geqslant n}^{m_1 + \cdots + m_N = n} } {\left( {_{m_1 , \cdots m_N }^n } \right)E_{m_1 } (q) \cdots E_{m_N } (q),}$

where N and n are positive integers, B _m (q) _n stand for the Apostol-Bernoulli numbers, E _m (q) for the Apostol-Euler numbers, and \(\left( {\begin{array}{*{20}c} n \\ {m_1 , \cdots ,m_N } \\ \end{array} } \right) = \frac{{n!}}{{m_1 ! \cdots m_N !}}.\) Our formulas involve Stirling numbers of the first kind. We also derive results for Apostol-Bernoulli and Apostol-Euler polynomials. As an application, for q = 1 we recover results of Dilcher, and our paper can be regarded as a q-extension of that of Dilcher.

     

Thermodynamic properties of water–aliphatic alcohol systems in a wide range of parameters

The results of thermal and thermodynamic (phase diagram) property calculations of water–aliphatic alcohol (methanol, ethanol, n-propanol) systems in liquid and vapor phases, as well as supercritical fluid water–methanol systems have been presented. The calculations are based on the polynomial equation of state, represented by expansion of the compressibility factor into a power series of reduced density (ω = ρ/ρ_cr and reduced temperature (τ = T/T _cr)

$$Z = \frac{p}{{RT{\rho _m}}} = 1 + \sum\limits_{i = 1}^m {\sum\limits_{j = 0}^{{n_i}} {\frac{{{a_{ij}}{\omega ^i}}}{{{\tau ^j}}}} } $$

, which describes experimental p,ρ,T,x-dependencies with an average relative error of 1.2%.