A posteriori error analysis of a quadratic finite volume method for nonlinear elliptic problems

In this article, we construct and analyze a residual-based a posteriori error estimator for a quadratic finite volume method (FVM) for solving nonlinear elliptic partial differential equations with homogeneous Dirichlet boundary conditions. We shall prove that the a posteriori error estimator yields the global upper and local lower bounds for the norm error of the FVM. So that the a posteriori error estimator is equivalent to the true error in a certain sense. Numerical experiments are performed to illustrate the theoretical results.     

Potts models with magnetic field: Arithmetic,geometry, and computation

We give a sheaf theoretic interpretation of Potts models with external magnetic field, in terms of constructible sheaves and their Euler characteristics. We show that the polynomial countability question for the hypersurfaces defined by the vanishing of the partition function is affected by changes in the magnetic field: elementary examples suffice to see non-polynomially countable cases that become polynomially countable after a perturbation of the magnetic field. The same recursive formula for the Grothendieck classes, under edge-doubling operations, holds as in the case without magnetic field, but the closed formulae for specific examples like banana graphs differ in the presence of magnetic field. We give examples of computation of the Euler characteristic with compact support, for the set of real zeros, and find a similar exponential growth with the size of the graph. This can be viewed as a measure of topological and algorithmic complexity. We also consider the computational complexity question for evaluations of the polynomial, and show both tractable and NP-hard examples, using dynamic programming.     

First excitations of the spin 1/2 Heisenberg antiferromagnet on the kagomé lattice

We study the exact low energy spectra of the spin 1/2 Heisenberg antiferromagnet on small samples of the kagomé lattice of up to N=36 sites. In agreement with the conclusions of previous authors, we find that these low energy spectra contradict the hypothesis of Néel type long range order. Certainly, the ground state of this system is a spin liquid, but its properties are rather unusual. The magnetic () excitations are separated from the ground state by a gap. However, this gap is filled with nonmagnetic () excitations. In the thermodynamic limit the spectrum of these nonmagnetic excitations will presumably develop into a gapless continuum adjacent to the ground state. Surprisingly, the eigenstates of samples with an odd number of sites, i.e. samples with an unsaturated spin, exhibit symmetries which could support long range chiral order. We do not know if these states will be true thermodynamic states or only metastable ones. In any case, the low energy properties of the spin 1/2 Heisenberg antiferromagnet on the kagomé lattice clearly distinguish this system from either a short range RVB spin liquid or a standard chiral spin liquid. Presumably they are facets of a generically new state of frustrated two-dimensional quantum antiferromagnets. Received: 27 November 1997 / Accepted: 29 January 1998     

355nm Nd∶YAG激光在H_2中的高效一级斯托克斯转换

对脉冲Nd∶YAG激光(355 nm)在H2和H2∶He-Ar混合气体中的受激拉曼散射(SRS)进行了研究。在0.5 MPa的氢气中,同时测量到从二级反斯托克斯到三级斯托克斯的多波长输出,其总转化效率达88%;而高压下只剩下一级和二级斯托克斯输出,其中二级斯托克斯最大能量转化效率达44%(对应量子效率为63%)。由于高级斯托克斯的竞争,纯氢气中一级斯托克斯的最大能量转换效率不超过43%。通过向3 MPa氢气中掺入2 MPaAr气后,很好地抑制了二级斯托克斯的产生,从而获得了能量转换效率高达71%(对应量子效率为83%)的一级斯托克斯输出。对四波混频和级联受激拉曼散射在氢气多级斯托克斯产生中的作用以及惰性气体对它们的影响进行了讨论。     

Kalman filtering of a space-time Markov random field

In this paper, finite-dimensional recursive filters for space-time Markov random fields are derived. These filters can be used with the expectation maximization (EM) algorithm to yield maximum likelihood estimates of the parameters of the model.     

Effective magnetic properties of a composite material with circular conductive elements

Effective magnetic properties of a composite meta-material consisting of periodically arranged circular conductive elements are studied theoretically. A general expression for the effective bulk permeability is obtained with mutual effects and lattice ordering being taken into account. The resonance frequency of the permeability is found to be strongly dependent on the size and shape of the unit cell. Frequency dispersion of the permeability is studied with special attention paid to the frequency range, where negative values of the permeability are possible. Corresponding recommendations for optimisation of the meta-materials with negative permeability are made. The results are confirmed by numerical simulations of the finite structure behaviour in an external magnetic field. Received 19 April 2002 Published online 31 July 2002     

Comparisons between model equations for long waves

Summary Considered here are model equations for weakly nonlinear and dispersive long waves, which feature general forms of dispersion and pure power nonlinearity. Two variants of such equations are introduced, one of Korteweg-de Vries type and one of regularized long-wave type. It is proven that solutions of the pure initial-value problem for these two types of model equations are the same, to within the order of accuracy attributable to either, on the long time scale during which nonlinear and dispersive effects may accumulate to make an order-one relative difference to the wave profiles.This research was supported in part by the National Science Foundation. A considerable portion of the project was completed while the first author was resident at the Institute for Mathematics and Its Applications, University of Minnesota.     

Global and uniform convergence of subspace correction methods for some convex optimization problems

This paper gives some global and uniform convergence estimates for a class of subspace correction (based on space decomposition) iterative methods applied to some unconstrained convex optimization problems. Some multigrid and domain decomposition methods are also discussed as special examples for solving some nonlinear elliptic boundary value problems.

     

Cosmic dust and our origins

The small solid particles in the space between the stars provide the surfaces for the production of many simple and complex molecules. Processes involving the effects of ultraviolet irradiation of the thin (hundredth micron) mantles are shown to produce a wide range of molecules and ions also seen in comets. Some of the more complex ones inferred from laboratory experiments are expected to play an important role in the origin of life. An outline of the chemical evolution of interstellar dust as observed and as studied in the laboratory is presented. Observations of comets are shown to provide substantial evidence for their being fluffy aggregates of interstellar dust as it was in the protosolar nebula, i.e. the interstellar cloud which collapsed to form the solar system. The theory that comets may have brought the progenitors of life to the earth is summarized.     

Non-universal critical behaviour of heteronuclear dimers on a square lattice––A Monte Carlo study

Monte Carlo simulation within the grand canonical ensemble, the histogram reweighting technique, and finite size scaling analysis are used to explore the phase behaviour of heteronuclear dimers, composed of A and B type atoms, on a square lattice. We have found that for the models with attractive BB and AB nearest-neighbour energy, u_BB=u_AB=−1, and for non-repulsive energy between AA nearest-neighbour sites, u_AA<0, the system belongs to the universality class of the two-dimensional Ising model. However, when u_AA>0, the system exhibits a non-universal critical behaviour. We have evaluated the dependences of the critical point characteristics on the value of u_AA.