Theoretical and numerical analysis on multispectral bioluminescence tomography

Recently, molecular imaging has been rapidly developed to studyphysiological and pathological processes in vivo at the cellularand molecular levels. Among molecular imaging modalities, opticalimaging has attracted a major attention for its unique advantages.In this paper, we establish a mathematical framework for multispectralbioluminescence tomography (BLT) that allows simultaneous studiesof multiple optical reporters. We show solution existence, uniquenessand continuous dependence on data as well as the limiting behaviourswhen the regularization parameter approaches zero or when thepenalty parameter approaches infinity. Then, we propose twonumerical schemes for multispectral BLT and derive error estimatesfor the corresponding solutions.     

OnQ-order andR-order of convergence

We give sufficient conditions for a sequence to have theQ-order and/or theR-order of convergence greater than one. If an additional condition is satisfied, then the sequence has an exactQ-order of convergence. We show that our results are sharp and we compare them with older results.This work was supported in part by the National Science Foundation under Grant No. DMS-85-03365. The author wishes to thank J. E. Dennis and R. A. Tapia for helpful comments, and the referee for pointing out a number of typographical and mathematical errors in the original version of this paper.     

Uniform convergence of Galerkin‐multigrid methods for nonconforming finite elements for second‐order problems with less than full elliptic regularity

In this article we prove uniform convergence estimates for the recently developed Galerkin‐multigrid methods for nonconforming finite elements for second‐order problems with less than full elliptic regularity. These multigrid methods are defined in terms of the “Galerkin approach,” where quadratic forms over coarse grids are constructed using the quadratic form on the finest grid and iterated coarse‐to‐fine intergrid transfer operators. Previously, uniform estimates were obtained for problems with full elliptic regularity, whereas these estimates are derived with less than full elliptic regularity here. Applications to the nonconforming P₁, rotated Q₁, and Wilson finite elements are analyzed. The result applies to the mixed method based on finite elements that are equivalent to these nonconforming elements. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 203–217, 2002; DOI 10.1002/num.10004     

The modified characteristic finite difference fractional steps method for the coupled system of fluid dynamics in porous media and its analysis

For the coupled system of multilayer fluid dynamics in porous media, the modified characteristic finite difference fractional steps method applicable to parallel arithmetic is put forward and two‐dimensional and three‐dimensional schemes are used to form a complete set. Some techniques, such as calculus of variations, energy method, piecewise biquadratic interpolation, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates are adopted. Optimal order estimates in L² norm are derived to determine the error in the approximate solution. This method has already been applied to the numerical simulation of multilayer fluid dynamics in porous media. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 665–681, 2003.     

The convergence of the Rayleigh-Ritz Method in quantum chemistry

Two sufficient criteria for the convergence of the Rayleigh-Ritz Method (RRM) with respect to the eigenvalues (E-convergence) of non-relativistic electronic Hamiltonians of molecules are discussed and compared. Moreover, a necessary and sufficient criterion is given. By example (Sect. 9) it is shown that the L ²-completeness of the basis is not sufficient to guarantee E-convergence. The convergence of the wave functions in different norms (-convergence) is also investigated. In particular, sufficient conditions for the one-particle basis functions (orbitals) are given, such that a CI calculation in this basis is both E- and -convergent.     

Grouped coordinate minimization using Newton's method for inexact minimization in one vector coordinate

Letf(x,y) be a function of the vector variablesx R ⁿ andy R ^m. The grouped (variable) coordinate minimization (GCM) method for minimizingf consists of alternating exact minimizations in either of the two vector variables, while holding the other fixed at the most recent value. This scheme is known to be locally,q-linearly convergent, and is most useful in certain types of statistical and pattern recognition problems where the necessary coordinate minimizers are available explicitly. In some important cases, the exact minimizer in one of the vector variables is not explicitly available, so that an iterative technique such as Newton's method must be employed. The main result proved here shows that a single iteration of Newton's method solves the coordinate minimization problem sufficiently well to preserve the overall rate of convergence of the GCM sequence.The authors are indebted to Professor R. A. Tapia for his help in improving this paper.     

Rare numbers

Suppose thatX ₁,X ₂,... is a sequence of iid random variables taking values inZ ⁺. Consider the random sequenceA(X)(X ₁,X ₂,...). LetY _n be the number of integers which appear exactly once in the firstn terms ofA(X). We investigate the limit behavior ofn ^–(1–) Y _n for [0, 1].     

Approximations to Isentropic Planar Magneto-Hydrodynamics Equations by Relaxed Euler-Type Systems

In this paper, the authors consider an approximation to the isentropic planar Magneto-hydrodynamics (MHD for short) equations by a kind of relaxed Euler-type system. The approximation is based on the generalization of the Maxwell law for nonNewtonian fluids together with the Maxwell correction for the Amp`ere law, hence the approximate system becomes a first-order quasilinear symmetrizable hyperbolic systems with partial dissipation. They establish the global-in-time smooth solutions to the approximate Euler-type equations in a small neighbourhood of constant equilibrium states and obtain the global-in-time convergence towards the isentropic planar MHD equations. In addition, they also establish the global-in-time error estimates of the limit based on stream function techniques and energy estimates for error variables.     

无穷积分余项与其被积函数比值的收敛阶

对一类函数的无穷积分余项与该函数的比值得到当x趋于无穷大时的收敛阶,这类函数是幂函数与指数函数的乘积函数,并将其应用到Mittag-Leffler函数.同时考虑了对应的级数情形.