Stability, mobility and power currents in a two-dimensional model for waveguide arrays with nonlinear coupling

A two-dimensional nonlinear Schrödinger lattice with nonlinear coupling, modelling a square array of weakly coupled linear optical waveguides embedded in a nonlinear Kerr material, is studied. We find that despite a vanishing energy difference (Peierls-Nabarro barrier) of fundamental stationary modes the mobility of localized excitations is very poor. This is attributed to a large separation in parameter space of the bifurcation points of the involved stationary modes. At these points the stability of the fundamental modes is changed and an asymmetric intermediate solution appears that connects the points. The control of the power flow across the array when excited with plane waves is also addressed and shown to exhibit great flexibility that may lead to applications for power-coupling devices. In certain parameter regimes, the direction of a stable propagating plane-wave current is shown to be continuously tunable by amplitude variation (with fixed phase gradient). More exotic effects of the nonlinear coupling terms like compact discrete breathers and vortices, and stationary complex modes with nontrivial phase relations are also briefly discussed. Regimes of dynamical linear stability are found for all these types of solutions.     

Large scale dissipation and filament instability in two-dimensional turbulence

Coherent vortices in two-dimensional turbulence induce far-field effects that stabilize vorticity filaments and inhibit the generation of new vortices. We show that the large-scale energy sink often included in numerical simulations of statistically stationary two-dimensional turbulence reduces the stabilizing role of the vortices, leading to filament instability and to continuous formation of new coherent vortices. This counterintuitive effect sheds new light on the mechanisms responsible for vortex formation in forced-dissipated two-dimensional turbulence, and it has significant impact on the temporal evolution of the vortex population in freely decaying turbulence. The time dependence of vortex statistics in the presence of a large-scale energy sink can be approximately described by a modified version of the scaling theory developed for small-scale dissipation.     

Modeling of nonlinear and non-stationary multi-vortex behavior of CDWs at nanoscales in restricted geometries of internal junctions

We report on studies of stationary states and their transient dynamic for an incommensurate charge density wave (ICDW) in a restricted geometry of two spatial dimensions. The model takes into account multiple fields in mutual nonlinear interactions: the amplitude and the phase of the complex order parameter, and distributions of the electric and chemical potentials, of the density and the current of normal carriers. We observed spontaneous formation of vortices (the ICDW dislocations), and followed events of their creation and the subsequent evolution. The vortices appear when the voltage across, or the current through, the sample exceed a threshold. The number of vortices remnant in the reconstructed stationary state increases stepwise – in agreement with experiments, while a much greater number of vortices appears during the intermediate transient states. The vortex core concentrates the electric dipole leading to sharp drops of the electric and chemical potentials across the core. That can lead to enhanced inter-layer tunneling making the core to be a self-tuned microscopic tunneling junction. The results are applied to experiments on nano-fabricated mesa-junctions. They also appeal to modern efforts of the field-effect transformations in correlated electronic systems.     

A finite element method with mesh adaptivity for computing vortex states in fast-rotating Bose–Einstein condensates

Numerical computations of stationary states of fast-rotating Bose–Einstein condensates require high spatial resolution due to the presence of a large number of quantized vortices. In this paper we propose a low-order finite element method with mesh adaptivity by metric control, as an alternative approach to the commonly used high-order (finite difference or spectral) approximation methods. The mesh adaptivity is used with two different numerical algorithms to compute stationary vortex states: an imaginary time propagation method and a Sobolev gradient descent method. We first address the basic issue of the choice of the variable used to compute new metrics for the mesh adaptivity and show that refinement using simultaneously the real and imaginary part of the solution is successful. Mesh refinement using only the modulus of the solution as adaptivity variable fails for complicated test cases. Then we suggest an optimized algorithm for adapting the mesh during the evolution of the solution towards the equilibrium state. Considerable computational time saving is obtained compared to uniform mesh computations. The new method is applied to compute difficult cases relevant for physical experiments (large nonlinear interaction constant and high rotation rates).     

Topological defects in incommensurate magnetic and crystal structures and quasi-periodic solutions of the elliptic sine-Gordon equation

Vortex-type singular solutions with a topological charge of the elliptic sine-Gordon equation have been studied. One- and two-dimensional vortex lattices on a homogeneous and periodic background are constructed in the explicit form using the Bäcklund transformation. The interaction of vortices is investigated and finite energy configurations are found. On the basis of the obtained results new topological defects in incommensurate magnetic and crystal structures are predicted and described. The interaction of vortex magnetic structures with nonlinear spin waves is considered.     

Energy inconsistency of Galerkin approximationsfor plane Couette flow

For classical solutions of the incompressible Navier-Stokes equations (NSE) the energybalance between kinetic energy, work done by external forces, and viscous dissipation holds rigorously true. It is shown in this paper that standard Galerkin approximations violate energy balance in the case of plane Couette flow, whereas Poiseuille flow turns out to be energy consistent at any cutoff. The main reason for this discrepancy is seen in the different boundary conditions between the stationary linear shear flow and its disturbances. In our analysis, essentially, we introduce an auxiliary external force field which enforces the finite dimensional Galerkin approximation to fulfil the NSE. It is exemplarily demonstrated how the energy discrepancy decreases when the number of disturbed modes is increased which couple to the stationary shear flow.     

A note on Chern–Simons solitons—A type III vortex from the wall vortex

We study some properties of topological Chern–Simons vortices in 2+1 dimensions. As has already been understood in the past, in the large magnetic flux limit, they are well described by a Chern–Simons domain wall, which has been compactified on a circle with the symmetric phase inside and the asymmetric phase on the outside.Our goal is two-fold. First we want to explore how the tension depends on the magnetic flux discretized by the integer n. The BPS case is already known, but not much has been explored about the non-BPS potentials. A generic renormalizable potential has two dimensionless parameters that can be varied. Variation of only one of them leads to a type I and type II vortex, very similar to the Abrikosov–Nielsen–Olesen (ANO) case. Variation of both the parameters leads to a much richer structure. In particular we have found a new type of vortex, which is type I-like for small flux and then turns type II-like for larger flux. We could tentatively denote it a type III vortex. This results in a stable vortex with number of fluxes which can be greater than one.Our second objective is to study the Maxwell–Chern–Simons theory and understand how the large n limit of the CS vortex is smoothly connected with the large n limit of the ANO vortex.     

Statistical mechanics of geophysical turbulence: application to jovian flows and Jupiter’s great red spot

We propose a parameterization of 2D geophysical turbulence in the form of a relaxation equation similar to a generalized Fokker–Planck equation [P.H. Chavanis, Phys. Rev. E 68 (2003) 036108]. This equation conserves circulation and energy and increases a generalized entropy functional determined by a prior vorticity distribution fixed by small-scale forcing [R. Ellis, K. Haven, B. Turkington, Nonlinearity 15 (2002) 239]. We discuss applications of this formalism to jovian atmosphere and Jupiter’s great red spot. We show that, in the limit of small Rossby radius where the interaction becomes short-range, our relaxation equation becomes similar to the Cahn–Hilliard equation describing phase ordering kinetics. This strengthens the analogy between the jet structure of the great red spot and a “domain wall”. Our relaxation equation can also serve as a numerical algorithm to construct arbitrary nonlinearly dynamically stable stationary solutions of the 2D Euler equation. These solutions can represent jets and vortices that emerge in 2D turbulent flows as a result of violent relaxation. Due to incomplete relaxation, the statistical prediction may fail and the system can settle on a stationary solution of the 2D Euler equation which is not the most mixed state. In that case, it can be useful to construct more general nonlinearly dynamically stable stationary solutions of the 2D Euler equation in an attempt to reproduce observed phenomena.     

Non-BCS superconductivity for underdoped cuprates by spin-vortex attraction

Within a gauge approach to the t-J model, we propose a new, non-BCS mechanism of superconductivity for underdoped cuprates. The gluing force of the superconducting mechanism is an attraction between spin vortices on two different Néel sublattices, centered around the empty sites described in terms of fermionic holons. The spin fluctuations are described by bosonic spinons with a gap generated by the spin vortices. Due to the no-double occupation constraint, there is a gauge attraction between holon and spinon binding them into a physical hole. Through gauge interaction the spin vortex attraction induces the formation of spin-singlet (RVB) spinon pairs with a lowering of the spinon gap. Lowering the temperature, the approach exhibits two crossover temperatures: at the higher crossover a finite density of incoherent holon pairs are formed leading to a reduction of the hole spectral weight, while at the lower crossover a finite density of incoherent spinon RVB pairs are also formed, giving rise to a gas of incoherent preformed hole pairs, and magnetic vortices appear in the plasma phase. Finally, at a even lower temperature the hole pairs become coherent, the magnetic vortices becoming dilute and superconductivity appears. The superconducting mechanism is not of BCS-type since it involves a gain in kinetic energy (for spinons) coming from the spin interactions.     

Towards a Landau–Ginzburg-Type Theory for Granular Fluids

In this paper we show how, under certain restrictions, the hydrodynamic equations for the freely evolving granular fluid fit within the framework of the time dependent Landau–Ginzburg (LG) models for critical and unstable fluids. The granular fluid, which is usually modeled as a fluid of inelastic hard spheres (IHS), exhibits two instabilities: the spontaneous formation of vortices and of high density clusters. We suppress the clustering instability by imposing constraints on the system sizes, in order to illustrate how LG-equations can be derived for the order parameter, being the rate of deformation or shear rate tensor, which controls the formation of vortex patterns. From the shape of the energy functional we obtain the stationary patterns in the flow field. Quantitative predictions of this theory for the stationary states agree well with molecular dynamics simulations of a fluid of inelastic hard disks.     

Vortex motion in superconducting ladders

We study vortex motion in superconducting ladders. Since the single-vortex excitation energy is finite in this system, even in the absence of an external magnetic field there may exist excess vortices that are stabilized by external currents. The motion of the excess vortices induces localized voltage drops, given by a nonlinear power law, at currents less than the critical current of the junctions. The localized voltage drop moves with constant velocity and its amplitude is proportional to the velocity. Above the critical current of the junctions there is continuous creation and annihilation of vortex-antivortex pairs.     

Self-consistent turbulence in the two-dimensional nonlinear Schrödinger equation with a repulsive potential

The dynamics of dark solitons (vortices) with the same topological charge (vorticity) in the two-dimensional nonlinear Schr?dinger (NLS) equation in a defocusing medium is studied. The dynamics differ from those in incompressible media due to the possibility of energy and angular momentum radiation. The problem of the breakup of a multicharged dark soliton, which is a local decrease of the wave function intensity, into a number of chaotically moving vortices with single charge, is studied both analytically and numerically. After an initial period of intensive wave radiation, there emerges a nonuniform, steady turbulent self-organized motion of these vortices which is restricted in space by the size of the potential well of the initial multicharged dark soliton. Separate orbits of finite widths arise in this turbulent motion. That is, the statistical probability to observe a vortex in a given point has maxima near certain points (orbit positions). In spite of the fact that numerical calculations were performed in a finite region, the turbulent distributions of the vortices do not depend on the size of the container when its radius is larger than the size of the potential well of the primary multicharged dark soliton. The steady turbulent distribution of vortices on these orbits can be obtained as the extremal of the Lyapunov functional of the NLS equation, and obeys some simple rules. The first is the absence of Cherenkov resonance with linear (sound) waves. The second is the condition of a potential energy maximum in the region of vortex motion. These conditions give an approximately equidistant disposition of orbits of the same number of vortices on each orbit, which corresponds to a constant rotating velocity. The magnitude of this velocity is mainly determined by the sound velocity. An integral estimation of the self-consistent rotation of the vortex zone is given.     

Exact results for the 1D asymmetric exclusion process and KPZ fluctuations

We discuss several properties of the current of the one-dimensional asymmetric simple exclusion process (ASEP) through exact solutions. First we explain the stationary measure for the finite system with boundaries and its average current. Then we study the fluctuation properties of the current and the position of a tagged particle for the infinite system in the Kardar-Parisi-Zhang (KPZ) scaling region.     

Stability of giant vortices in quantum liquids

We show how giant vortices can be stabilized for strong external potentials in Bose-Einstein condensates. We illustrate the formation of these vortices thanks to the Ginzburg-Landau dissipative dynamics for two typical potentials in two spatial dimensions. The giant vortex stability is studied for the particular case of a rotating cylindrical hard wall. Due to axial symmetry the minimization of the perturbed energy is simplified into a one dimensional relaxation dynamics. Solving this 1D minimization problem, we observe that giant vortices are either never stable, or only stable in a finite frequency range. Finally we obtain the marginal curve for the minimum frequency needed to observe a giant vortex.     

相邻双侧边盖驱动方腔流动的三维线性整体稳定性

文章考察了相邻双侧边盖驱动方腔流动（即上壁面向右运动和左侧壁面向下运动）的三维线性整体稳定性.首先,采用Taylor-Hood有限元方法并经由Newton迭代过程计算得到双侧边盖驱动方腔流动的二维稳态基本流.其次,Taylor-Hood有限元在Chebyshev Gauss配置点上进行离散,同时Gauss配置点也可以用于线性稳定性方程的高阶有限差分格式离散.然后,离散得到的矩阵形式的广义特征值问题可以结合shift-and-invert算法采用隐式重启Arnoldi方法计算.最后,通过对线性稳定性方程特征值的计算,发现了一个最不稳定的驻定模态和两对对称行波模态.最不稳定的三维驻定模态的临界Reynolds数为Re_c=261.5,远远小于二维不稳定的临界Reynolds数Re_c2d=1 061.7.通过画出这3类三维不稳定模态的流向扰动速度和扰动涡量的空间等值面图像,可以发现不稳定扰动位于稳态基本流的两个主涡区域,因此可以认为主涡区域是三维扰动失稳的主要能量来源地.