It is known that weak interactions of two solitary waves in generalized nonlinear Schrödinger (NLS) equations exhibit fractal dependence on initial conditions, and the dynamics of these interactions is governed by a universal two-degree-of-freedom ODE system [Y. Zhu J. Yang, Universal fractal structures in the weak interaction of solitary waves in generalized nonlinear Schrödinger equations, Phys. Rev. E 75 (2007) 036605]. In this paper, this ODE system is analyzed comprehensively. Using asymptotic methods along separatrix orbits, a simple second-order map is derived. This map does not have any free parameters after variable rescalings, and thus is universal for all weak interactions of solitary waves in generalized NLS equations. Comparison between this map’s predictions and direct simulations of the ODE system shows that the map can capture the fractal-scattering phenomenon of the ODE system very well both qualitatively and quantitatively.

Solitons of the power-energy saturation (PES) equation are studied using adiabatic perturbation theory. In the anomalous regime individual soliton pulses are found to be well approximated by solutions of the classical nonlinear Schrödinger (NLS) equation with the key parameters of the soliton changing slowly as they evolve. Evolution equations are found for the pulse amplitude(s), velocity(ies), position(s), and phase(s) using integral relations derived from the PES equation. The results from the integral relations are shown to agree with multi-scale perturbation theory. It is shown that the single soliton case exhibits mode-locking behavior for a wide range of parameters, while the higher states form effective bound states. Using the fact that there is weak overlap between tails of interacting solitons, evolution equations are derived for the relative amplitudes, velocities, positions, and phase differences. Comparisons of interacting soliton behavior between the PES equation and the classical NLS equation are also exhibited.

We consider two-stage hybrid protocols that combine quantum resources and classical resources to generate classical correlations shared by two separated players. Our motivation is twofold. First, in the near future, the scale of quantum information processing is quite limited, and when quantum resource available is not sufficient for certain tasks, a possible way to strengthen the capability of quantum schemes is introducing extra classical resources. We analyze the mathematical structures of these hybrid protocols, and characterize the relation between the amount of quantum resources and classical resources needed. Second, a fundamental open problem in communication complexity theory is to describe the advantage of sharing prior quantum entanglement over sharing prior randomness, which is still widely open. It turns out that our quantum and classical hybrid protocols provide new insight into this important problem.

Many problems in nature, being characterized by a parameter, are of interest both with a fixed parameter value and with the parameter approaching an asymptotic limit. Numerical schemes that are convergent in both regimes offer robust discretizations, which can be highly desirable in practice. The asymptotically compatible schemes studied in an earlier published version of this paper meet such objectives for a class of parametrized problems. An extended version of the abstract mathematical framework is established rigorously here, together with applications to the numerical solution of both nonlocal models and their local limits. In particular, the framework can be applied to nonlocal models of diffusion and a general state-based peridynamic system parametrized by the horizon radius. Recent findings have exposed the risks associated with some discretizations of nonlocal models when the horizon radius is proportional to the discretization parameter. Thus, it is desirable to develop asymptotically compatible schemes for such models so as to offer robust numerical discretizations to problems involving nonlocal interactions on multiple scales. This work provides new insight in this regard through a careful analysis of related conforming finite element discretizations, and the finding is valid under minimal regularity assumptions on exact solutions. It reveals that for the nonlocal models under consideration and their local limit, as long as the finite element space contains continuous piecewise linear functions, the Galerkin finite element approximation is always asymptotically compatible. For piecewise constant finite elements, whenever applicable, it is shown that a correct local limit solution can also be obtained as long as the discretization (mesh) parameter decreases faster than the modeling (horizon) parameter does. These results can be used to guide future computational studies of nonlocal problems. Some other applications, such as the fractional PDE limit of nonlocal models, and open questions are also presented

Matrix concentration inequalities provide information about the probability that a random matrix is close to its expectation with respect to the l2 operator norm. This paper uses semigroup methods to derive sharp nonlinear matrix inequalities. The main result is that the classical Bakry–Émery curvature criterion implies subgaussian concentration for “matrix Lipschitz” functions. This argument circumvents the need to develop a matrix version of the log-Sobolev inequality, a technical obstacle that has blocked previous attempts to derive matrix concentration inequalities in this setting. The approach unifies and extends much of the previous work on matrix concentration. When applied to a product measure, the theory reproduces the matrix Efron–Stein inequalities due to Paulin et al. It also handles matrix-valued functions on a Riemannian manifold with uniformly positive Ricci curvature.