Patterns in passive optical systems with boundary effects

Pattern formation is closely related to system boundary conditions in nonlinear dynamic systems.apart from the system control parameters,To avoid complexity,system boundary conditions are usually considered to be infinite or periodic,and the initial conditions spatially homogeneous,But it is not always the case in real situations,or sometimes periodic boundary conditions are not exact.To show the important and interesting boundary effects in real pattern formation,we suggest a simple universal boundary condition in a typical optical pattern formation system.Numerical simulations of the passive optical system show that pattern characteristics such as distribution symmetry,peak number,structure strength,evolution course and stability are all greatly influenced by the system boundary conditions.     

Symmetry boundary conditions

A simple approach to energy conserving boundary conditions using exact symmetries is described which is especially useful for numerical simulations using the finite difference method. Each field in the simulation is normally either symmetric (even) or antisymmetric (odd) with respect to the simulation boundary. Another possible boundary condition is an antisymmetric perturbation about a nonzero value. One of the most powerful aspects of this approach is that it can be easily implemented in curvilinear coordinates by making the scale factors of the coordinate transformation symmetric about the boundaries. The method is demonstrated for magnetohydrodynamics (MHD), reduced MHD, and a hybrid code with particle ions and fluid electrons. These boundary conditions yield exact energy conservation in the limit of infinite time and space resolution. Also discussed is the interpretation that the particle charge reverses sign at a conducting boundary with boundary normal perpendicular to the background magnetic field.     

The discrete Frenkel-Kontorova model and its extensions: I. Exact results for the ground-states

We present a rigorous study of the classical ground-states under boundary conditions of a class of one-dimensional models generalizing the discrete Frenkel-Kontorova model. The extremalization equations of the energy of these models turn out to define area preserving twist maps which exhibits periodic, quasi-periodic and chaotic orbits. For all boundary conditions, we select among all the extremum solutions of the energy of the model, those which correspond to the ground-states of the infinite system. We prove that these ground-states are either periodic (commensurate) or quasi-periodic (incommensurate) but are never chaotic. We also prove the existence of elementary discommensurations which are minimum energy configuration of the model for certain special boundary conditions. The topological structure of the whole set of ground-states is described in details. In addition to physical applications, consequences for twist map homeomorphisms are mentioned. Part II (S. Aubry, P.Y. LeDaeron and G. Andre) will be mostly devoted to exact results on the transition by breaking of analyticity which occurs on the incommensurate ground states when the model parameters vary and on its connection with the stochasticity threshold in the corresponding twist map.     

Symmetries and Boundary Conditions with a Twist

Interest in finite-size systems has risen in the last decades, due to the focus on nanotechnological applications and because they are convenient for numerical treatment that can subsequently be extrapolated to infinite lattices. Independently of the envisioned application, special attention must be given to boundary condition, which may or may not preserve the symmetry of the infinite lattice. Here, we present a detailed study of the compatibility between boundary conditions and conservation laws. The conflict between open boundary conditions and momentum conservation is well understood, but we examine other symmetries, as well: we discuss gauge invariance, inversion, spin, and particle-hole symmetry and their compatibility with open, periodic, and twisted boundary conditions. In the interest of clarity, we develop the reasoning in the framework of the one-dimensional half-filled Hubbard model, whose Hamiltonian displays a variety of symmetries. Our discussion includes analytical and numerical results. Our analytical survey shows that, as a rule, boundary conditions break one or more symmetries of the infinite-lattice Hamiltonian. The exception is twisted boundary condition with the special torsion Θ = πL/2, where L is the lattice size. Our numerical results for the ground-state energy at half-filling and the energy gap for L = 2–7 show how the breaking of symmetry affects the convergence to the L → ∞ limit. We compare the computed energies and gaps with the exact results for the infinite lattice drawn from the Bethe-Ansatz solution. The deviations are boundary-condition dependent. The special torsion yields more rapid convergence than open or periodic boundary conditions. For sizes as small as L = 7, the numerical results for twisted condition are very close to the L → ∞ limit. We also discuss the ground-state electronic density and magnetization at half filling under the three boundary conditions.     

Derivation of the amplitude equation for a hydromagnetic convective problem

The amplitude equation for a convective system under a vertical magnetic field is derived. The coefficients in these equations have been numerically calculated for infinite Prandtl number fluids and for boundary conditions both free and rigid top and bottom. The results confirm that, for realistic parameter values only stationary convection can be present and the pattern is made by convective rolls.     

Derivation of the amplitude equation for a hydromagnetic convective problem

The amplitude equation for a convective system under a vertical magnetic field is derived. The coefficients in these equations have been numerically calculated for infinite Prandtl number fluids and for boundary conditions both free and rigid top and bottom. The results confirm that, for realistic parameter values only stationary convection can be present and the pattern is made by convective rolls.     

Derivation of the amplitude equation for a hydromagnetic convective problem

The amplitude equation for a convective system under a vertical magnetic field is derived. The coefficients in these equations have been numerically calculated for infinite Prandtl number fluids and for boundary conditions both free and rigid top and bottom. The results confirm that, for realistic parameter values only stationary convection can be present and the pattern is made by convective rolls.     

Influence of boundary conditions on square bond percolation nearp c

The influence of boundary conditions on square bond percolation for system sizes ranging from 10×10 to 240×240 is studied for the quantitiesP _∞, χ, the effective percolation threshold and the finite-size scaling relations forP _∞ and χ. The Monte Carlo simulations suggest that free edges approximate the infinite system as well as the more complicated periodic boundary conditions.     

General Boundary Conditions for a Majorana Single-Particle in a Box in (1 + 1) Dimensions

We consider the problem of a Majorana single-particle in a box in (1 + 1) dimensions. We show that the most general set of boundary conditions for the equation that models this particle is composed of two families of boundary conditions, each one with a real parameter. Within this set, we only have four confining boundary conditions—but infinite not confining boundary conditions. Our results are also valid when we include a Lorentz scalar potential in this equation. No other Lorentz potential can be added. We also show that the four confining boundary conditions for the Majorana particle are precisely the four boundary conditions that mathematically can arise from the general linear boundary condition used in the MIT bag model. Certainly, the four boundary conditions for the Majorana particle are also subject to the Majorana condition.     

超对称WKB近似与一维无限深势阱

证明了一维无限深势阱为非形状不变势.采用超对称WKB近似,得到了一维无限深势阱的精确能谱公式.计算结果表明,把势场具有形状不变性作为超对称WKB近似能够给出精确能谱的充要条件,是一个过强的提法,需要加以弱化.     

Droplet States in the XXZ Heisenberg Chain

We consider the ground states of the ferromagnetic XXZ chain with spin up boundary conditions. The ground state of this model, restricted to a sector with a fixed number of down spins, describes a droplet of down spins in an environment of up spins. We find the exact energy and the states that describe these droplets in the limit of an infinite number of down spins. We prove that there is a gap in the spectrum above the droplet states. As the XXZ Hamiltonian has a gap above the fully magnetized ground states as well, this means that the droplet states (for sufficiently large droplets) form an isolated band. The width of this band tends to zero in the limit of infinitely large droplets. We also prove the analogous results for finite chains with periodic boundary conditions and for the infinite chain. Received: 5 September 2000 / Accepted: 8 December 2000     

Scattering by cavities of arbitrary shape in an infinite plate and associated vibration problems

This paper presents a solution for the displacement of a uniform elastic thin plate with an arbitrary cavity, modelled using the biharmonic plate equation. The problem is formulated as a system of boundary integral equations by factorizing the biharmonic equation, with the unknown boundary values expanded in terms of a Fourier series. At the edge of the cavity we consider free-edge, simply-supported and clamped boundary conditions. Methods to suppress ill-conditioning which occurs at certain frequencies are discussed, and the combined boundary integral equation method is implemented to control this problem. A connection is made between the problem of an infinite plate with an arbitrary cavity and the vibration problem of an arbitrarily shaped finite plate, using the jump discontinuity present in single-layer distributions at the boundary. The first few frequencies and modes of displacement are computed for circular and elliptic cavities, which provide a check on our numerics, and results for the displacement of an infinite plate are given for four specific cavity geometries and various boundary conditions.     

Topological phase boundary in a generalized Kitaev model

We study the effects of the next-nearest-neighbor hopping and nearest-neighbor interactions on topological phases in a one-dimensional generalized Kitaev model. In the noninteracting case, we define a topological number and calculate exactly the phase diagram of the system. With addition of the next-nearest-neighbor hopping, the change of phase boundary between the topological and trivial regions can be described by an effective shift of the chemical potential. In the interacting case, we obtain the entanglement spectrum, the degeneracies of which correspond to the topological edge modes, by using the infinite time-evolving block decimation method. The results show that the interactions change the phase boundary as adding an effective chemical potential which can be explained by the change of the average number of particles.     

直接模拟中不同边界条件的实施及对沉降规律的影响

在牛顿流体中, 对颗粒在4种不同边界的垂直通道中的沉降运动进行了直接数值模拟. 计算结果表明:通过计算区域随颗粒运动而移动构建的无限长通道能准确模拟颗粒自由下落到稳定沉降的发展过程; 周期性边界条件由于流场变化, 对颗粒沉降产生了影响, 不能模拟颗粒的自由沉降过程; 底部封闭边界适合模拟封闭容器内颗粒与固壁的相互作用过程, 若颗粒达到稳定沉降, 也能模拟无限长通道内的沉降过程; 流化边界适合模拟流化床内气固两相流动. 计算结果有助于更好地理解和使用不同边界条件. 关键词： 直接数值模拟 边界条件 沉降 任意拉格朗日-欧拉方法     

A new integral equation method for direct electromagnetic scattering in homogeneous media and its numerical confirmation

In this paper, we derive a new integral equation method for direct electromagnetic scattering in homogeneous media and present a numerical confirmation of the new method via a computer simulation. The new integral equation method is based on a paper written by DeSanto [1], originally for scattering from an infinite rough surface separating homogeneous dielectric half-spaces. Here, it is applied to a bounded scatterer, which can be an ohmic conductor or a dielectric, with some simplification of the continuity conditions for the fields. The new integral equation method is developed by choosing the electric field and its normal derivative as boundary unknowns, which are not the usual boundary unknowns. The new integral equation method may provide significant computational advantages over the standard Stratton-Chu method [2] because it leads to a 50% sparse, rather than 100% dense, impedance (collocation) matrix. Our theoretical development of the new integral equation method is exact.     

Soliton and boundary condition induced fractional fermion number

We show that for a fermion in a bounded background potential in a finite box, eigenvalues of the total charge are independent of whether the potential is solitonic and depend only on the boundary condition: half-odd integral or integral for charge conjugation (C) invariant boundary conditions and an arbitrary fraction forC non-invariant boundary conditions. Fractional fermion numbers for infinite space Jackiw-Rebbi and Goldstone-Wilczek Hamiltonians are reproduced in finite space by using boundary conditions different from the periodic ones of Rajaraman and Bell.     

Magnetic and dynamic properties of the Hubbard model in infinite dimensions

An essentially exact solution of the infinite dimensional Hubbard model is made possible by using a self-consistent mapping of the Hubbard model in this limit to an effective single impurity Anderson model. Solving the latter with quantum Monte Carlo procedures enables us to obtain exact results for the one and two-particle properties of the infinite dimensional Hubbard model. In particular, we find antiferromagnetism and a pseudogap in the single-particle density of states for sufficiently large values of the intrasite Coulomb interaction at half filling. Both the antiferromagnetic phase and the insulating phase above the Néel temperature are found to be quickly suppressed on doping. The latter is replaced by a heavy electron metal with a quasiparticle mass strongly dependent on doping as soon asn<1. At half filling the antiferromagnetic phase boundary agrees surprisingly well in shape and order of magnitude with results for the three dimensional Hubbard model.     

Infinite Number of Integrals of Motion in Classically Integrable System With Boundary (Ⅱ)

In Affine Toda field theory, links among three generating functions for integrals of motion derived from P. (Ⅰ) are studied, and some classically integrable boundary conditions are obtained. An infinite number of integrals of motion are calculated in ZMS model with quasi-periodic condition. We find the classically integrable boundary conditions and K± matrices of ZMS model with independent boundary conditions on each end. It is identified that an infinite number of integrals of motion does exist and one of them is the Hamiltonian, so this system is completely integrable.     

A Monte Carlo study of dipolar hard spheres The pair distribution function and the dielectric constant

This paper is a systematic investigation of the effects of boundary conditions upon Monte Carlo calculations for dipolar fluids. Results are reported for the minimum image (MI), spherical cut-off (SC) and uniform reaction field methods. All three approximations are shown to give different pair distribution functions, g(12), and none yields the infinite system result. It is concluded that theories giving g(12) for an infinite system should not be evaluated by direct comparison with Monte Carlo results. Two alternative methods by which meaningful comparisons can be made are described in the text. The dependence of the thermodynamic properties upon boundary conditions is important only at large values of the dipole moment. For small to moderate dipoles both MI and SC are found to give reasonable estimates of the dielectric constant.