GPU accelerated Monte Carlo simulation of the 2D and 3D Ising model

The compute unified device architecture (CUDA) is a programming approach for performing scientific calculations on a graphics processing unit (GPU) as a data-parallel computing device. The programming interface allows to implement algorithms using extensions to standard C language. With continuously increased number of cores in combination with a high memory bandwidth, a recent GPU offers incredible resources for general purpose computing. First, we apply this new technology to Monte Carlo simulations of the two dimensional ferromagnetic square lattice Ising model. By implementing a variant of the checkerboard algorithm, results are obtained up to 60 times faster on the GPU than on a current CPU core. An implementation of the three dimensional ferromagnetic cubic lattice Ising model on a GPU is able to generate results up to 35 times faster than on a current CPU core. As proof of concept we calculate the critical temperature of the 2D and 3D Ising model using finite size scaling techniques. Theoretical results for the 2D Ising model and previous simulation results for the 3D Ising model can be reproduced.     

Solution of a diffusion equation of smoke molecules due to constant point source in gravitational field

From both the Langevin equation, including a gravitational term, the Fokker-Planck equation based on the dynamical behavior of Brownian particles, and equation of smoke molecular diffusion due to a constant point source is introduced and is solved by applying the Laplace transformation with the convolution theorem.The solution is expressed by the complementary error function with a mean pathway of smoke molecules affected by gravity and is proved to be reduced to conventional forms by a certain restriction neglecting gravity without any forced term.     

A roe-average algorithm for a granular-gas model with non-conservative terms

A Roe-average algorithm has been derived for a granular-gas model, proposed by Goldshtein and Shapiro [Goldshtein, Shapiro, Mechanics of collisional motion of granular materials: Part 1. General hydrodynamic equations, J. Fluid Mech. 282 (1995) 75–114], which contains non-conservative terms in the Euler-like hyperbolic governing equations apart from sink terms, which arise from inelastic collision of granules and are present only in the energy equation. The non-conservative terms introduce non-isentropic effects in acoustic-wave propagation within granular media and they also contribute to the Rankine–Hugoniot relations across a discontinuity. A Roe-average algorithm, based on the same granular-gas model, was derived in the literature [V. Kamenetsky, A. Goldshtein, M. Shapiro, D. Degani, Evolution of a shock wave in a granular gas, Phys. Fluids, 12 (2000) 3036–3049] which then required the implementation of a shock-fitting technique at a discontinuity. In the present work, Roe-averaged variables have been obtained from the Rankine–Hugoniot jump relations and the non-conservative terms have been incorporated in the numerical flux formula consistent with upwind principles associated with the granular speed of sound. Results for unsteady one-dimensional granular flows, colliding with a wall, demonstrate the capability of the proposed algorithm to capture strong shocks in addition to flow features not found in molecular gases, such as a fluidized region downstream of the shock and a compacted solid-block region adjacent to the wall.     

Lax representation with spectral parameter on a torus for integrable particle systems

Complete integrability is proved for the most general class of systems of interacting particles on a straight line with the Hamiltonian including elliptic functions of coordinates, depending on seven arbitrary parameters and having the structure defined by the root systems of the classical Lie algebras. The Lax representation for them depends on the spectral parameter given on a complex torus /, where is the lattice of periods of the Jacobi functions dependent on the Hamiltonian parameters. The possibility of constructing explicit solutions to the equations of motion is discussed.     

Time-dependent linear Hamiltonian systems and quantum mechanics

In this Letter, a theorem on time-dependent linear Hamiltonian systems is recalled and its connection with the Schrödinger equation is discussed. The kernel of the evolution operator of such quantum systems is computed. Furthermore, the Lewis and Riesenfeld theory for systems with many degrees of freedom is generalized.     

A high-order super-grid-scale absorbing layer and its application to linear hyperbolic systems

We continue the development of the super-grid-scale model initiated in [T. Colonius, H. Ran, A super-grid-scale model for simulating compressible flow on unbounded domains, J. Comput. Phys. 182 (1) (2002) 191–212] and consider its application to linear hyperbolic systems. The super-grid-scale model consists of two parts: reduction of an unbounded to a bounded domain by a smooth coordinate transformation and a damping of those scales. For linear problems the super-grid scales are analogous to spurious numerical waves. We damp these waves by high-order undivided differences. We compute reflection coefficients for different orders of the damping and find that significant improvements are obtained when high-order damping is used.     

Scattering theory for a model of an atom in a quantized field

We compute the scattering operator for a simplified model of an atom interacting a with a quantized field. The field is restricted to the vacuum and one-particle sectors, and the atom has only two states. We also solve the inverse scattering problem for the same model. The methods used rely on the particular form of the interaction, which is chosen to mimic the interaction between radiation and matter.This work was partially supported by NSF Grant DMS-8922941.     

Big q-Jacobi polynomials,q-Hahn polynomials,and a family of quantum 3-spheres

The big q-Jacobi polynomials and the q-Hahn polynomials are realized as spherical functions on a new quantum SU _q (2)-space which can be regarded as the total space of a family of quantum 3-spheres.     

Synchronization of switched system and application in communication

In this Letter, synchronization of switched system is investigated based on Lyapunov method. A sufficient condition is derived to ensure the synchronization between two switched systems, and a new communication scheme is also proposed based on this. Furthermore, some secure analysis works, such as return map attack and moving average error attack, are also given to show the security of the proposed scheme. Finally, simulation examples are constructed to verify the theoretical analysis and its application for communication.     

Geometric quantization of an OSp(1/2) coadjoint orbit

In this Letter, we show how the complete geometric quantization extends to specific supersymplectic supermanifolds. More precisely, we extend this procedure to OSp(1/2)-coadjoint orbits, which are graded extensions of elliptic Sp(2, )-coadjoint orbits. Our approach exploits results obtained in a previous work, where the notion of a super-Kähler supermanifold was defined, and the former orbits were shown to be nontrivial examples of such a notion. As their underlying Kähler manifolds, these supermanifolds carry a natural (super-Kähler) polarization, a crucial notion that was so far lacking. Geometric quantization leads here to a nontrivial representation of osp(1/2), which is realized in a space of square integrable holomorphic sections of a super-Hermitian complex line bundle sheaf-with-connection over the homogenous space OSp(1/2)/U(1).     

The star exponential and path integrals on compact groups II: The weyl character formula

This is a continuation of the work begun by Cadavid and Nakashima inLett. Math. Phys. 23, 111–115 (1991). An expression for the Weyl Character Formula is obtained in terms of the star-product path integral; and the relationship between the star-product path integral and the path integral developed on coadjoint orbits is established.     

A model of enhanced STM current due to semiconductor optical absorption

We present a simple model for the change in tunneling current between a semiconductor surface and a metal tip under spectroscopic illumination in a scanning tunneling microscope. This model predicts a sharp increase in the tunneling current due to the increase in the conduction band carrier density when the photon energy exceeds the optical band gap. The tunneling current for a large diffusion length has a more pronounced onset than for a small length. Our model should provide, when combined with experiments, a method of determining localized effective stoichiometry, and therefore provides a localized alternative to the use of optical absorption measurements. Our theoretical tunneling current versus photon energy curves are in good qualitative agreement with the existing experimentally measured curves for Si, GaAs, and InP obtained by Qian and Wessels. In addition, we have examined the effects of temperature, surface recombination velocity, and degeneracy on our theoretical results for the Hg_1−xCd_xTe, Hg_1−xZn_xTe, and Hg_1−xZn_xSe ternary narrow gap semiconductor systems.     

On the dynamics and the temperature states of the spin-boson model

Recently, the finite temperature states of the spin-boson model were studied. Here with different techniques in the case of level-splitting zero, the explicit expressions for both the dynamics and the equilibrium states are derived, namely Trotter's product formula and analyticity methods from perturbation theory are used.     

Geometric realizations,curvature decompositions,and Weyl manifolds

We show that any Weyl curvature model can be geometrically realized by a Weyl manifold.     

An improved algorithm and its parallel implementation for solving a general blood-tissue transport and metabolism model

Fast algorithms for simulating mathematical models of coupled blood-tissue transport and metabolism are critical for the analysis of data on transport and reaction in tissues. Here, by combining the method of characteristics with the standard grid discretization technique, a novel algorithm is introduced for solving a general blood-tissue transport and metabolism model governed by a large system of one-dimensional semilinear first order partial differential equations. The key part of the algorithm is to approximate the model as a group of independent ordinary differential equation (ODE) systems such that each ODE system has the same size as the model and can be integrated independently. Thus the method can be easily implemented in parallel on a large-scale multiprocessor computer. The accuracy of the algorithm is demonstrated for solving a simple blood-tissue exchange model introduced by Sangren and Sheppard [W.C. Sangren, C.W. Sheppard, A mathematical derivation of the exchange of a labeled substance between a liquid flowing in a vessel and an external compartment, Bull. Math. Biophys. 15 (1953) 387–394], which has an analytical solution. Numerical experiments made on a distributed-memory parallel computer (an HP Linux cluster) and a shared-memory parallel computer (a SGI Origin 2000) demonstrate the parallel efficiency of the algorithm.     

On a ‘super’ version of Lie's third fundamental theorem

Invoking the techniques of nonstandard analysis, we associate a super Lie group in the sense of M. Rothstein to any finite-dimensional super Lie module over an arbitrary graded-commutative Banach algebra, which obeys duality (in the sense that it may be recaptured from its dual coalgebra).     

A superresolution optical disk with a combination of optical path filtering and saturable absorption techniques

Calculations have been done for the modulation transfer function of a superresolution optical disk with a combination of optical path filtering and saturable absorption techniques. The results reveal that optical path filtering alone can only equalize, to a certain degree, the spatial frequency response of the opticaldisk system, but it can not attain any response beyond the cutoff frequency. With an additional saturable absorption layer in-between the substrate and the information layer, it is found that the optical response both under and beyond the cutoff frequency can be greatly improved.