T-matrix method for electromagnetic scattering from scatterers with complex structure

We describe a T-matrix program for light scattering calculations from particles with complex structure. The code treats the cases of homogeneous, layered and composite scatterers. These results are combined with basic results concerning the scattering by inhomogeneous scatterers and aggregates to apply to more general types of scatterers. Some numerical simulations are presented.     

The discrete Gaussian chain with 1/r n interactions: Exact results

We show that the discrete Gaussian chain with interactionV(r) = 1/(r ²–1/4) is self-dual. At the dual temperaturek _B T = 1 we calculate the height-height correlation function and find that the system is rough. A duality relation is established for the temperature-dependent correlation function exponent. We also consider interactionsV(r)–1/r ⁿ and show that absence of a phase transition for 2 <n < 3 implies absence of a phase transition for 1 <n < 2. All these results have their counterparts in a linear system of charges interacting through a potential which is asymptotically logarithmic (forn = 2) or power-law-like (forn 2.On leave of absence from Chemistry Laboratory III, Universitetsparken 5, 2100 ©, Copenhagen, Denmark.     

Linearization of the interaction principle: Analytic Jacobians in the “Radiant” model

In this paper we present a new linearization of the Radiant radiative transfer model. Radiant uses discrete ordinates for solving the radiative transfer equation in a multiply-scattering anisotropic medium with solar and thermal sources, but employs the adding method (interaction principle) for the stacking of reflection and transmission matrices in a multilayer atmosphere. For the linearization, we show that the entire radiation field is analytically differentiable with respect to any surface or atmospheric parameter for which we require Jacobians (derivatives of the radiance field). Derivatives of the discrete ordinate solutions are based on existing methods developed for the LIDORT radiative transfer models. Linearization of the interaction principle is completely new and constitutes the major theme of the paper. We discuss the application of the Radiant model and its linearization in the Level 2 algorithm for the retrieval of columns of carbon dioxide as the main target of the Orbiting Carbon Observatory (OCO) mission.     

Improved Lower Bounds for the Critical Probability of Oriented Bond Percolation in Two Dimensions

We present a coupled decreasing sequence of random walks on Z that dominate the edge process of oriented bond percolation in two dimensions. Using the concept of random walk in a strip, we describe an algorithm that generates an increasing sequence of lower bounds that converges to the critical probability of oriented percolation p_c. From the 7th term on, these lower bounds improve upon 0.6298, the best rigorous lower bound at present, establishing 0.63328 as a rigorous lower bound for p_c. Finally, a Monte Carlo simulation technique is presented; the use thereof establishes 0.64450 as a non-rigorous five-digit-precision (lower) estimate for p_c. Mathematics Subject Classification (1991): 60K35 Supported by CNPq (grant N.301637/91-1). Supported by a grant from CNPq.     

GPU accelerated simulations of 3D deterministic particle transport using discrete ordinates method

Graphics Processing Unit (GPU), originally developed for real-time, high-definition 3D graphics in computer games, now provides great faculty in solving scientific applications. The basis of particle transport simulation is the time-dependent, multi-group, inhomogeneous Boltzmann transport equation. The numerical solution to the Boltzmann equation involves the discrete ordinates (S_n) method and the procedure of source iteration. In this paper, we present a GPU accelerated simulation of one energy group time-independent deterministic discrete ordinates particle transport in 3D Cartesian geometry (Sweep3D). The performance of the GPU simulations are reported with the simulations of vacuum boundary condition. The discussion of the relative advantages and disadvantages of the GPU implementation, the simulation on multi GPUs, the programming effort and code portability are also reported. The results show that the overall performance speedup of one NVIDIA Tesla M2050 GPU ranges from 2.56 compared with one Intel Xeon X5670 chip to 8.14 compared with one Intel Core Q6600 chip for no flux fixup. The simulation with flux fixup on one M2050 is 1.23 times faster than on one X5670.     

光学元件双转子抛光技术中类高斯型去除函数建模与逼近

在Preston假设基础上,研究并推导了光学元件抛光过程中"花瓣"型截面磨头在双转子运动形式下的去除函数.提出了"花瓣"型截面磨头的离散点弧长表达方法,并讨论了该方法的适用范围;基于离散点弧长法,推导了"花瓣"型截面磨头在双转子运动形式下的去除函数表达式;根据可调步长式曲线逼近原理研究了类高斯型目标去除函数的逼近算法;并得到了类余弦型去除函数的双转子抛光模参量.从而证明了双转子抛光技术可以得到类高斯型去除函数.     

Impulsive generalized function synchronization of complex dynamical networks

This Letter investigates generalized function synchronization of continuous and discrete complex networks by impulsive control. By constructing the reasonable corresponding impulsively controlled response networks, some criteria and corollaries are derived for the generalized function synchronization between the impulsively controlled complex networks, continuous and discrete networks are both included. Furthermore, the generalized linear synchronization and nonlinear synchronization are respectively illustrated by several examples. All the numerical simulations demonstrate the correctness of the theoretical results.     

A hidden analytic structure of the Rabi model

The Rabi model describes the simplest interaction between a cavity mode with a frequency ω_c${\omega }_c$ and a two-level system with a resonance frequency ω₀${\omega }_0$. It is shown here that the spectrum of the Rabi model coincides with the support of the discrete Stieltjes integral measure in the orthogonality relations of recently introduced orthogonal polynomials. The exactly solvable limit of the Rabi model corresponding to Δ=ω₀/(2ω_c)=0$\Delta ={\omega }_0/\left(2{\omega }_c\right)=0$, which describes a displaced harmonic oscillator, is characterized by the discrete Charlier polynomials in normalized energy ?$?$, which are orthogonal on an equidistant lattice. A non-zero value of Δ$\Delta$ leads to non-classical discrete orthogonal polynomials ?_k(?)${?}_k\left(?\right)$ and induces a deformation of the underlying equidistant lattice. The results provide a basis for a novel analytic method of solving the Rabi model. The number of ca. 1350 calculable energy levels per parity subspace obtained in double precision (cca 16 digits) by an elementary stepping algorithm is up to two orders of magnitude higher than is possible to obtain by Braak’s solution. Any first n$n$ eigenvalues of the Rabi model arranged in increasing order can be determined as zeros of ?_N(?)${?}_N\left(?\right)$ of at least the degree N=n+n_t$N=n+n_t$. The value of n_t>0$n_t>0$, which is slowly increasing with n$n$, depends on the required precision. For instance, n_t?26$n_t?26$ for n=1000$n=1000$ and dimensionless interaction constant κ=0.2$\kappa =0.2$, if double precision is required. Given that the sequence of the l$l$th zeros x_nl$x_{nl}$’s of ?_n(?)${?}_n\left(?\right)$’s defines a monotonically decreasing discrete flow with increasing n$n$, the Rabi model is indistinguishable from an algebraically solvable model in any finite precision. Although we can rigorously prove our results only for dimensionless interaction constant κ<1$\kappa <1$, numerics and exactly solvable example suggest that the main conclusions remain to be valid also for κ≥1$\kappa \ge 1$.     

Long range effect of ion irradiation on diffusion

Argon ion irradiation enhanced diffusion between Cu and Ni was found beyond the radiation effect zone in classical theory. The enhanced diffusion effect attenuates along with the distance between the diffusion zone and the radiated surface. An interpretation based on irradiation induced discrete breathers was employed as a possible mechanism.     

Alternating direction implicit schemes for the two-dimensional fractional sub-diffusion equation

New numerical techniques are presented for the solution of a two-dimensional anomalous sub-diffusion equation with time fractional derivative. In these methods, standard central difference approximation is used for the spatial discretization, and, for the time stepping, two new alternating direction implicit (ADI) schemes based on the L₁ approximation and backward Euler method are considered. The two ADI schemes are constructed by adding two different small terms, which are different from standard ADI methods. The solvability, unconditional stability and H¹ norm convergence are proved. Numerical results are presented to support our theoretical analysis and indicate the efficiency of both methods.