Stencils and problem partitionings: Their influence on the performance of multiple processor systems

Given a discretization stencil, partitioning the problem domain is an important first step for the efficient solution of partial differential equations on multiple processor systems. Partitions are derived that minimize interprocessor communication when the number of processors is known a priori and each domain partition is assigned to a different processor. This partitioning technique uses the stencil structure to select appropriate partition shapes. For square problem domains, it is shown that non-standard partitions (e.g., hexagons) are frequently preferable to the standard square partitions for a variety of commonly used stencils. This investigation is concluded with a formalization of the relationship between partition shape, stencil structure, and architecture, allowing selection of optimal partitions for a variety of parallel systems

Parallel, iterative solution of sparse linear systems: Models and architectures

A model of a general class of asynchronous, iterative solution methods for linear systems is developed. In the model, the system is solved by creating several cooperating tasks that each compute a portion of the solution vector. A data transfer model predicting both the probability that data must be transferred between two tasks and the amount of data to be transferred is presented. This model is used to derive an execution time model for predicting parallel execution time and an optimal number of tasks given the dimension and sparsity of the coefficient matrix and the costs of computation, synchronization, and communication. The suitability of different parallel architectures for solving randomly sparse linear systems is discussed. Based on the complexity of task scheduling, one parallel architecture, based on a broadcast bus, is presented and analyzed

A model of asynchronous iterative algorithms for solving large, sparse, linear systems

Solving large, sparse, linear systems of equations is one of the fundamental problems in large scale scientific and engineering computation. A model of a general class of asynchronous, iterative solution methods for linear systems is developed. In the model, the system is solved by creating several cooperating tasks that each compute a portion of the solution vector. This model is then analyzed to determine the expected intertask data transfer and task computational complexity as functions of the number of tasks. Based on the analysis, recommendations for task partitioning are made. These recommendations are a function of the sparseness of the linear system, its structure (i.e., randomly sparse or banded), and dimension

Sharp eigenvalue enclosures for the perturbed angular Kerr-Newman Dirac operator

A certified strategy for determining sharp intervals of enclosure for the eigenvalues of matrix differential operators with singular coefficients is examined. The strategy relies on computing the second order spectrum relative to subspaces of continuous piecewise linear functions. For smooth perturbations of the angular Kerr-Newman Dirac operator, explicit rates of convergence due to regularity of the eigenfunctions are established. Existing benchmarks are validated and sharpened by several orders of magnitude in the unperturbed setting.Comment: 27 pages, 2 figures, 5 tables. Some errors fixe

The Linear Boltzmann Equation as the Low Density Limit of a Random Schrodinger Equation

We study the evolution of a quantum particle interacting with a random potential in the low density limit (Boltzmann-Grad). The phase space density of the quantum evolution defined through the Husimi function converges weakly to a linear Boltzmann equation with collision kernel given by the full quantum scattering cross section.Comment: 74 pages, 4 figures, (Final version -- typos corrected

A new proof of the analyticity of the electronic density of molecules

We give a new, short proof of the regularity away from the nuclei of the electronic density of a molecule obtained in [1,2]. The new argument is based on the regularity properties of the Coulomb interactions underlined in [3,4] and on well-known elliptic technics. [1] S. Fournais, M. Hoffmann-Ostenhof, T. Hoffmann-Ostenhof, T. Oe stergaard Soerensen: The electron density is smooth away from the nuclei. Comm. Math. Phys. 228, no. 3 (2002), 401-415. [2] S. Fournais, M. Hoffmann-Ostenhof, T. Hoffmann-Ostenhof, T. Oestergaard Soerensen: Analyticity of the density of electronic wave functions. Ark. Mat. 42, no. 1 (2004), 87-106. [3] W. Hunziker: Distortion analyticity and molecular resonances curves. Ann. Inst. H. Poincar\'e, s. A, t. 45, no 4, 339-358 (1986). [4] M. Klein, A. Martinez, R. Seiler, X.P. Wang: On the Born-Oppenheimer expansion for polyatomic molecules. Comm. Math. Phys. 143, no. 3, 607-639 (1992). The paper is published in Letters in Mathematical Physics 93, number 1, pp. 73-83, 2010. The original publication is available at " www.springerlink.com "

Long-Time Dynamics of Variable Coefficient mKdV Solitary Waves

We study the Korteweg-de Vries-type equation dt u=-dx(dx^2 u+f(u)-B(t,x)u), where B is a small and bounded, slowly varying function and f is a nonlinearity. Many variable coefficient KdV-type equations can be rescaled into this equation. We study the long time behaviour of solutions with initial conditions close to a stable, B=0 solitary wave. We prove that for long time intervals, such solutions have the form of the solitary wave, whose centre and scale evolve according to a certain dynamical law involving the function B(t,x), plus an H^1-small fluctuation.Comment: 19 page

Operations and single particle interferometry

Interferometry of single particles with internal degrees of freedom is investigated. We discuss the interference patterns obtained when an internal state evolution device is inserted into one or both the paths of the interferometer. The interference pattern obtained is not uniquely determined by the completely positive maps (CPMs) that describe how the devices evolve the internal state of a particle. By using the concept of gluing of CPMs, we investigate the structure of all possible interference patterns obtainable for given trace preserving internal state CPMs. We discuss what can be inferred about the gluing, given a sufficiently rich set of interference experiments. It is shown that the standard interferometric setup is limited in its abilities to distinguish different gluings. A generalized interferometric setup is introduced with the capacity to distinguish all gluings. We also connect to another approach using the well known fact that channels can be realized using a joint unitary evolution of the system and an ancillary system. We deduce the set of all such unitary `representations' and relate the structure of this set to gluings and interference phenomena.Comment: Journal reference added. Material adde

Effective Hamiltonians for atoms in very strong magnetic fields

We propose three effective Hamiltonians which approximate atoms in very strong homogeneous magnetic fields $B$ modelled by the Pauli Hamiltonian, with fixed total angular momentum with respect to magnetic field axis. All three Hamiltonians describe $N$ electrons and a fixed nucleus where the Coulomb interaction has been replaced by $B$-dependent one-dimensional effective (vector valued) potentials but without magnetic field. Two of them are solvable in at least the one electron case. We briefly sketch how these Hamiltonians can be used to analyse the bottom of the spectrum of such atoms.Comment: 43 page

Zero energy resonance and the logarithmically slow decay of unstable multilevel systems

The long time behavior of the reduced time evolution operator for unstable multilevel systems is studied based on the N-level Friedrichs model in the presence of a zero energy resonance.The latter means the divergence of the resolvent at zero energy. Resorting to the technique developed by Jensen and Kato [Duke Math. J. 46, 583 (1979)], the zero energy resonance of this model is characterized by the zero energy eigenstate that does not belong to the Hilbert space. It is then shown that for some kinds of the rational form factors the logarithmically slow decay of the reduced time evolution operator can be realized.Comment: 31 pages, no figure