A search technique for global optimization in a chaotic environment

We describe a new algorithm which uses the trajectories of a discrete dynamical system to sample the domain of an unconstrained objective function in search of global minima. The algorithm is unusually adept at avoiding nonoptimal local minima and successfully converging to a global minimum. Trajectories generated by the algorithm for objective functions with many local minima exhibit chaotic behavior, in the sense that they are extremely sensitive to changes in initial conditions and system parameters. In this context, chaos seems to have a beneficial effect: failure to converge to a global minimum from a given initial point can often be rectified by making arbitrarily small changes in the system parameters.     

On a class of operators related to second-order differential equations

In the present paper, we find a class of linear homogeneous differential equations of order n + 1 (n > 1) whose fundamental system of solutions is constructed from the fundamental system of solutions of a second-order differential equation. The spectral properties of differential operators generated by these differential expressions are investigated. In particular, sufficient conditions are obtained for the coefficients of a second-order differential equation under which the case of maximal deficiency indices is realized. Dedicated to the memory of B. M. Levitan     

High pressure X-ray diffraction study of CdAl2Se4 and Raman study of AAl2Se4 (A=Hg, Zn) and CdAl2X4 (X=Se, S)

We report the results of an X-ray diffraction study of CdAl₂Se₄ and of Raman studies of HgAl₂Se₄ and ZnAl₂Se₄ at room temperature, and of CdAl₂S₄ and CdAl₂Se₄ at 80 K at high pressure. The ambient pressure phase of CdAl₂Se₄ is stable up to a pressure of 9.1 GPa above which a phase transition to a disordered rock salt phase is observed. A fit of the volume pressure data to a Birch-Murnaghan type equation of state yields a bulk modulus of 52.1 GPa. The relative volume change at the phase transition at ∼9 GPa is about 10%. The analysis of the Raman data of HgAl₂Se₄ and ZnAl₂Se₄ reveals a general trend observed for different defect chalcopyrite materials. The line widths of the Raman peaks change at intermediate pressures between 4 and 6 GPa as an indication of the pressure induced two stage order-disorder transition observed in these materials. In addition, we include results of a low temperature Raman study of CdAl₂S₄ and CdAl₂Se₄, which shows a very weak temperature dependence of the Raman-active phonon modes.     

Smart Geodesic Tracing in GRworkbench

We have developed a new tool for numerical work in General Relativity: GRworkbench. We discuss how GRworkbench's implementation of a numerically-amenable analogue to Differential Geometry facilitates the development of robust and chart-independent numerical algorithms. We consider, as an example, geodesic tracing on two charts covering the exterior Schwarzschild space-time.     

On the Reality of the Eigenvalues for a Class of -Symmetric Oscillators

 We study the eigenvalue problem with the boundary conditions that decays to zero as z tends to infinity along the rays , where is a real polynomial and . We prove that if for some we have for all , then the eigenvalues are all positive real. We then sharpen this to a larger class of polynomial potentials. In particular, this implies that the eigenvalues are all positive real for the potentials when with , and with the boundary conditions that decays to zero as z tends to infinity along the positive and negative real axes. This verifies a conjecture of Bessis and Zinn-Justin. Received: 16 January 2002 / Accepted: 1 May 2002 Published online: 6 August 2002     

Eigenvalue distributions of random unitary matrices

 Let U be an n × n random matrix chosen from Haar measure on the unitary group. For a fixed arc of the unit circle, let X be the number of eigenvalues of M which lie in the specified arc. We study this random variable as the dimension n grows, using the connection between Toeplitz matrices and random unitary matrices, and show that (X -E [X])/(\Var (X))^1/2 is asymptotically normally distributed. In addition, we show that for several fixed arcs I ₁ , ..., I _m , the corresponding random variables are jointly normal in the large n limit. Received: 15 November 2000 / Revised version: 27 September 2001 / Published online: 17 May 2002     

A primer on perfect simulation for spatial point processes

This primer provides a self-contained exposition of the case where spatial birth-and-death processes are used for perfect simulation of locally stable point processes. Particularly, a simple dominating coupling from the past (CFTP) algorithm and the CFTP algorithms introduced in [13], [14], and [5] are studied. Some empirical results for the algorithms are discussed. Received: 30 June 2002     

Every diassociative A-loop is Moufang

An A-loop is a loop in which every inner mapping is an automorphism. A problem which had been open since 1956 is settled by showing that every diassociative A-loop is Moufang.

     

Three-dimensional radiative transfer in an anisotropically scattering, semi-infinite medium: generalized reflection function

Three-dimensional radiative transfer in an anisotropic scattering medium exposed to spatially varying, collimated radiation is studied. The generalized reflection function for a semi-infinite medium with a very general scattering phase function is the focus of this investigation. An integral transform is used to reduce the three-dimensional transport equation to a one-dimensional form, and a modified Ambarzumian's method is applied to formulate a nonlinear integral equation for the generalized reflection function. The integration is over both the polar and azimuthal angles; hence, the integral equation is said to be in the double-integral form. The double-integral, reflection function formulation can handle a variety of anisotropic phase functions and does not require an expansion of the phase function in a Legendre polynomial series. Complicated kernel transformations of previous single-integral studies are eliminated. Single and double scattering approximations are developed. Numerical results are presented for a Rayleigh phase function to illustrate the computational characteristics of the method and are compared to results obtained with the single-integral method. Agreement between the two approaches is excellent; however, as the transform variable increases beyond five the number of quadrature points required for the double-integral method to produce accurate solutions significantly increases. A new interpolation scheme produces accurate results when the transform variable is large.