An image system and surface singularity solutions for potential flow past an elliptical cylinder

The Milne-Thomson circle theorem is extended to give a simplegeneral expression for the image system in an elliptical cylinderintroduced into an otherwise specified unbounded potential flowwhich contains no singularities in the region to be occupiedby the ellipse. This image system is used to obtain an expressionfor the corresponding source-sink surface singularity distributionon the ellipse, thus providing new benchmark test cases forsource-sink solutions as obtained numerically by a panel method.Several typical examples are given to illustrate the generaltheoretical approach.     

Coarse master equation from Bayesian analysis of replica molecular dynamics simulations

We use Bayesian inference to derive the rate coefficients of a coarse master equation from molecular dynamics simulations. Results from multiple short simulation trajectories are used to estimate propagators. A likelihood function constructed as a product of the propagators provides a posterior distribution of the free coefficients in the rate matrix determining the Markovian master equation. Extensions to non-Markovian dynamics are discussed, using the trajectory "paths" as observations. The Markovian approach is illustrated for the filling and emptying transitions of short carbon nanotubes dissolved in water. We show that accurate thermodynamic and kinetic properties, such as free energy surfaces and kinetic rate coefficients, can be computed from coarse master equations obtained through Bayesian inference.     

Interlaced solitons and vortices in coupled DNLS lattices

In the present work, we propose a new set of coherent structures that arise in nonlinear dynamical lattices with more than one component, namely interlaced solitons. In the anti-continuum limit of uncoupled sites, these are waveforms whose one component has support where the other component does not. We illustrate systematically how one can combine dynamically stable unary patterns to create stable ones for the binary case of two-components. For the one-dimensional setting, we provide a detailed theoretical analysis of the existence and stability of these waveforms, while in higher dimensions, where such analytical computations are far more involved, we resort to corresponding numerical computations. Lastly, we perform direct numerical simulations to showcase how these structures break up, when they are exponentially or oscillatorily unstable, to structures with a smaller number of participating sites.     

Discrete solitons in nonlinear Schrödinger lattices with a power-law nonlinearity

We study the discrete nonlinear Schrödinger lattice model with the onsite nonlinearity of the general form, |u|^2σu. We systematically verify the conditions for the existence and stability of discrete solitons in the one-dimensional version of the model predicted by means of the variational approximation (VA), and demonstrate the following: monostability of fundamental solitons (FSs) in the case of the weak nonlinearity, 2σ+1<3.68; bistability, in a finite range of values of the soliton’s power, for 3.68<2σ+1<5; and the presence of a threshold (minimum norm of the FS), for 2σ+1≥5. We also perform systematic numerical simulations to study higher-order solitons in the same general model, i.e., bound states of the FSs. While all in-phase bound states are unstable, stability regions are identified for antisymmetric double solitons and their triple counterparts. These numerical findings are supplemented by an analytical treatment of the stability problem, which allows quantitively accurate predictions for the stability features of such multipulses. When these waveforms are found to be unstable, we show, by means of direct simulations, that they self-trap into a persistent lattice breather, or relax into a stable FS, or sometimes decay completely.