A Revisit on the Derivation of the Particular Solution for the Differential Operator ∆2 ± λ2

In this paper, we applied the polyharmonic splines as the basis functions to derive particular solutions for the differential operator ∆² ± λ². Similar to the derivation of fundamental solutions, it is non-trivial to derive particular solutions for higher order differential operators. In this paper, we provide a simple algebraic factorization approach to derive particular solutions for these types of differential operators in 2D and 3D. The main focus of this paper is its simplicity in the sense that minimal mathematical background is required for numerically solving higher order partial differential equations such as thin plate vibration. Three numerical examples in both 2D and 3D are given to validate particular solutions we derived.     

Monte Carlo simulation of spin-wave excitation for the magnetic films with structure and interaction fluctuations

Monte Carlo simulation studies are performed to examine influence of structure and interaction fluctuations on magnetic properties of a ferromagnetic system modelled with a Heisenberg Hamiltonian. It is found that the spontaneous magnetization at low temperature for the multilayered films decreases with temperature in a Bloch law of spin-wave excitations. Both Bloch coefficient B and exponent b vary evidently because of a strong surface and size effect in the finite magnetic films with free boundaries. For the disordered bulk FCC magnet with periodic boundary, the Bloch T^3/2 law is followed at low temperature and B is greatly influenced by the structure and interaction fluctuations. At the same time, Bloch coefficient B of the amorphous magnet with the coordination and interaction fluctuations has been derived. The simulated results are in good agreement with the theoretical predictions of spin-wave excitation, and explain the experimental facts well.     

A New Treatment for Some Periodic Schr?dinger Operators Ⅱ: The Wave Function

Following the approach of our previous paper we continue to study the asymptotic solution of periodic Schrödinger operators. Using the eigenvalues obtained earlier the corresponding asymptotic wave functions are derived. This gives further evidence in favor of the monodromy relations for the Floquet exponent proposed in the previous paper. In particular, the large energy asymptotic wave functions are related to the instanton partition function of N=2 supersymmetric gauge theory with surface operator. A relevant number theoretic dessert is appended.     

周期性磁谐振材料本构参数的理论分析

将薄的磁谐振介质板等效为面磁流, 利用周期性边界条件, 给出了面磁流的指数形式. 通过计算无穷个面磁流在不同空间位置上产生的总电场和总磁场, 推导出了周期性磁谐振人工材料的色散关系和布洛赫阻抗, 进而获取了布洛赫本构参数的理论计算公式. 由于考虑了磁谐振人工材料中的电反谐振对布洛赫介电常数和磁导率的影响, 所以基于仿真实验的布洛赫本构参数的提取值和布洛赫本构参数理论预测值之间的误差很小, 这说明本文推导的布洛赫本构参数的理论计算公式在描述周期性磁谐振材料的电磁特性方面是十分有效的. 这些理论公式将在解释磁谐振现象、设计和优化周期性磁谐振材料等方面提供重要的理论依据. 关键词： 周期性结构 磁谐振 布洛赫本构参数 面磁流     

Spectral properties of reduced Bloch Hamiltonians

Bloch Hamiltonians are defined, and the existence of bands is proven for a large class of periodic operators. The results are strong enough to include most of the reasonable physical models of a single electron in crystals. A notable exception is the Dirac Bloch Hamiltonian for a Coulombic crystal with high charge. Properties of the Bloch waves are briefly described and it is shown that “simple” Bloch Hamiltonians do not have Bloch waves with a finite number of Fourier coefficients. The asymptotic distribution of the bands is determined, and it is shown that for a large class of Hamiltonians, it is determined by the kinetic energy alone.     

Hofstadter spectrum in electric and magnetic fields

The problem of Bloch electrons in two dimensions subjected to magnetic and intense electric fields is investigated. Magnetic translations, electric evolution, and energy translation operators are used to specify the solutions of the Schrödinger equation. For rational values of the magnetic flux quanta per unit cell and commensurate orientations of the electric field relative to the original lattice, an extended superlattice can be defined and a complete set of mutually commuting space-time symmetry operators is obtained. Dynamics of the system is governed by a finite difference equation that exactly includes the effects of: an arbitrary periodic potential, an electric field orientated in a commensurable direction of the lattice, and coupling between Landau levels. A weak periodic potential broadens each Landau level in a series of minibands, separated by the corresponding minigaps. The addition of the electric field induces a series of avoided and exact crossing of the quasienergies, for sufficiently strong electric field the spectrum evolves into equally spaced discreet levels, in this “magnetic Stark ladder” the energy separation is an integer multiple of hE/aB, with a the lattice parameter.     

Finite-size scaling studies of one-dimensional reaction-diffusion systems. Part I. Analytical results

We consider two single-species reaction-diffusion models on one-dimensional lattices of lengthL: the coagulation-decoagulation model and the annihilation model. For the coagulation model the system of differential equations describing the time evolution of the empty interval probabilities is derived for periodic as well as for open boundary conditions. This system of differential equations grows quadratically withL in the latter case. The equations are solved analytically and exact expressions for the concentration are derived. We investigate the finite-size behavior of the concentration and calculate the corresponding scaling functions and the leading corrections for both types of boundary conditions. We show that the scaling functions are independent of the initial conditions but do depend on the boundary conditions. A similarity transformation between the two models is derived and used to connect the corresponding scaling functions.     

Damped low-frequency surface-plasmon resonance in metal-dielectric periodic structures of planar semi-infinite extent

In the low-frequency limit with respect to the bulk plasma frequency of metal, damped surface-plasmon resonance is examined for a periodic semi-infinite structure with metal-dielectric unit cells in slab geometry. In comparison to the author’s earlier results in [1], the additional material damping is found to alter the resonance characteristics in many nontrivial ways. In particular, the damped Bloch waves propagating in the direction normal to the slab planes are induced, thereby altering wave stability with respect to the ratio of dielectric constants.     

Ruijsenaars' Commuting Difference Operators as Commuting Transfer Matrices

For Belavin's elliptic quantum R-matrix, we construct an L-operator as a set of difference operators acting on functions on the type A weight space. According to the fundamental relation RLL=LLR, taking the trace of the L-operator gives a set of commuting difference operators. We show that for the above mentioned L-operator this approach gives Macdonald type operators with elliptic theta function coefficient, actually equivalent to Ruijsenaars' operators. The relationship between the difference L-operator and Krichever's Lax matrix is given, and an explicit formula for elliptic commuting differential operators is derived. We also study the invariant subspace for the system which is spanned by symmetric theta functions on the weight space. Received: 27 December 1995 / Accepted: 11 November 1996     

Solution of Time-Dependent Diffusion Equations with Variable Coefficients Using Multiwavelets

A new numerical algorithm is developed for the solution of time-dependent differential equations of diffusion type. It allows for an accurate and efficient treatment of multidimensional problems with variable coefficients, nonlinearities, and general boundary conditions. For space discretization we use the multiwavelet bases introduced by Alpert (1993,SIAM J. Math. Anal.24, 246–262), and then applied to the representation of differential operators and functions of operators presented by Alpert, Beylkin, and Vozovoi (Representation of operators in the multiwavelet basis, in preparation). An important advantage of multiwavelet basis functions is the fact that they are supported only on non-overlapping subdomains. Thus multiwavelet bases are attractive for solving problems in finite (non periodic) domains. Boundary conditions are imposed with a penalty technique of Hesthaven and Gottlieb (1996,SIAM J. Sci. Comput., 579–612) which can be used to impose rather general boundary conditions. The penalty approach was extended to a procedure for ensuring the continuity of the solution and its first derivative across interior boundaries between neighboring subdomains while time stepping the solution of a time dependent problem. This penalty procedure on the interfaces allows for a simplification and sparsification of the representation of differential operators by discarding the elements responsible for interactions between neighboring subdomains. Consequently the matrices representing the differential operators (on the finest scale) have block-diagonal structure. For a fixed order of multiwavelets (i.e., a fixed number of vanishing moments) the computational complexity of the present algorithm is proportional to the number of subdomains. The time discretization method of Beylkin, Keiser, and Vozovoi (1998, PAM Report 347) is used in view of its favorable stability properties. Numerical results are presented for evolution equations with variable coefficients in one and two dimensions.     

THE NONLINEAR SCHRÖDINGER EQUATION FOR THE DELTA-COMB POTENTIAL: QUASI-CLASSICAL CHAOS AND BIFURCATIONS OF PERIODIC STATIONARY SOLUTIONS

The nonlinear Schrödinger equation is studied for a periodic sequence of delta-potentials (a delta-comb) or narrow Gaussian potentials. For the delta-comb the time-independent nonlinear Schrödinger equation can be solved analytically in terms of Jacobi elliptic functions and thus provides useful insight into the features of nonlinear stationary states of periodic potentials. Phenomena well-known from classical chaos are found, such as a bifurcation of periodic stationary states and a transition to spatial chaos. The relation to new features of nonlinear Bloch bands, such as looped and period doubled bands, are analyzed in detail. An analytic expression for the critical nonlinearity for the emergence of looped bands is derived. The results for the delta-comb are generalized to a more realistic potential consisting of a periodic sequence of narrow Gaussian peaks and the dynamical stability of periodic solutions in a Gaussian comb is discussed.     

The hilbert problem for matrices and a new class of integrable equations

A series of new integrable nonlinear differential equations is derived as compatibility conditions between generalized Lax pairs of operators which are meromorphic functions of the spectral parameter on the Riemann surface S of genus 1. On employing the Hilbert problem for the surface S, a general method of integration of these equations is proposed. The method is applied to obtain soliton solutions for asymmetric chiral SU(2) theory.     

Simple explicit formulas for Gaussian path integrals with time-dependent frequencies

Quadratic fluctuations require an evaluation of ratios of functional determinants of second-order differential operators. We relate these ratios to the Green functions of the operators for Dirichlet, periodic and antiperiodic boundary conditions on a line segment. This permits us to take advantage of Wronski's construction method for Green functions without knowledge of eigenvalues. Our final formula expresses the ratios of functional determinants in terms of an ordinary 2 × 2 determinant of a constant matrix constructed from two linearly independent solutions of the homogeneous differential equations associated with the second-order differential operators. For ratios of determinants encountered in semiclassical fluctuations around a classical solution, the result can further be expressed in terms of this classical solution. In the presence of a zero mode, our method allows for a simple universal regularization of the functional determinants. For Dirichlet's boundary condition, our result is equivalent to Gelfand-Yaglom's. Explicit formulas are given for a harmonic oscillator with an arbitrary time-dependent frequency.     

Collective excitations of a periodic Bose condensate in the Wannier representation

We study the dispersion relation of the excitations of a dilute Bose-Einstein condensate confined in a periodic optical potential and its Bloch oscillations in an accelerated frame. The problem is reduced to one-dimensionality through a renormalization of the s-wave scattering length and the solution of the Bogolubov-de Gennes equations is formulated in terms of the appropriate Wannier functions. Some exact properties of a periodic one-dimensional condensate are easily demonstrated: (i) the lowest band at positive energy refers to phase modulations of the condensate and has a linear dispersion relation near the Brillouin zone centre; (ii) the higher bands arise from the superposition of localized excitations with definite phase relationships; and (iii) the wavenumber-dependent current under a constant force in the semiclassical transport regime vanishes at the zone boundaries. Early results by Slater [Phys. Rev. 87, 807 (1952)] on a soluble problem in electron energy bands are used to specify the conditions under which the Wannier functions may be approximated by on-site tight-binding orbitals of harmonic-oscillator form. In this approximation the connections between the low-lying excitations in a lattice and those in a harmonic well are easily visualized. Analytic results are obtained in the tight-binding scheme and are illustrated with simple numerical calculations for the dispersion relation and semiclassical transport in the lowest energy band, at values of the system parameters which are relevant to experiment. Received 3 December 1999 and Received in final form 22 March 2000     

Nonlinear Bloch waves

The nonlinear Bloch theorem for the temporal and spatial Schrödinger solitons in dispersive and nonlinear periodic structures is proved. It is shown that bright and dark solitary nonlinear Bloch waves exist only under certain conditions and that the parameter functions describing dispersion and nonlinearity periodic inhomogeneities cannot be chosen independently.     

Third-order constants of motion in quantum mechanics

We investigate the general form of a third-order linear differential operator that is required to commute with the Schrödinger Hamiltonian in two dimensions, and find that the third-order part must be a polynomial of third degree in the generators of the Euclidean group. Partial differential equations that the potentialV must satisfy are derived, and solved for the special cases where the Schrödinger equation separates in polar or Cartesian coordinates. The functionsV thus obtained are nonsingular, but are periodic through elliptic functions. After separation of variables, the Schrödinger equation gives Lame's equation.     

Symplectic Symmetries of Hamiltonian Systems

Abstract

We prove a generalization to the case of s × s matrix linear differential operators of the classical theorem of E. Cotton giving necessary and sufficient conditions for equivalence of eigenvalue problems for scalar linear differential operators. The conditions for equivalence to a matrix Schrödinger operator are derived and formulated geometrically in terms of vanishing conditions on the curvature of a gl(s, R)-valued connection. These conditions are illustrated on a class of matrix differential operators of physical interest, arising by symmetry reduction from Dirac’s equation for a spinor field minimally coupled with a cylindrically symmetric magnetic field.     

Floquet–Bloch solutions in a sawtooth photonic crystal

Band structure of a sawtooth photonic crystal for optical wave propagation along the axis of periodicity is investigated. Floquet–Bloch solutions are found and illustrated for the bandgaps, allowed bands, and bandedges of the crystal. Special attention is given to the cases where Floquet–Bloch solutions become periodic functions.     

Form factors and correlation functions of an interacting spinless fermion model

Introducing the fermionic R-operator and solutions of the inverse scattering problem for local fermion operators, we derive a multiple integral representation for zero-temperature correlation functions of a one-dimensional interacting spinless fermion model. Correlation functions particularly considered are the one-particle Green's function and the density–density correlation function both for any interaction strength and for arbitrary particle densities. In particular for the free fermion model, our formulae reproduce the known exact results. Form factors of local fermion operators are also calculated for a finite system.