A Self-Similar Dynamics in Viscous Spheres

We study the evolution of radiating and viscous fluid spheres assuming an additional homothetic symmetry on the spherically symmetric space-time. We match a very simple solution to the symmetry equations with the exterior one (Vaidya). We then obtain a system of two ordinary differential equations which rule the dynamics, and find a self-similar collapse which is shear-free and with a barotropic equation of state. Considering a huge set of initial self-similar dynamics states, we work out a model with an acceptable physical behavior.     

Energy Conserving, Liouville, and Symplectic Integrators

In this paper, we construct an integrator that converves volume in phase space. We compare the results obtained using this method and a symplectic integrator. The results of our experiments do not reveal any superiority of the symplectic over strictly volume-preserving integrators. We also investigate the effect of numerically conserving energy in a numerical process by rescaling velocities to keep energy constant at every step. Our results for Henon-Heiles problem show that keeping energy constant in this way destroys ergodicity and forces the solution onto a periodic orbit.     

Inclusive (γ, N), (γ, NN) and (γ, Nπ) reactions in nuclei at intermediate energies

Differential cross sections of nucleons excited in photonuclear reactions in medium and heavy nuclei are studied by considering all relevant reaction mechanisms leading to the excitation of protons or neutrons. We take advantage of previous microscopic studies for the absorption and scattering of photons and photoproduced pions, and implement a simulation code in order to take into account the propagation of the nucleons as well as their collisions with other nucleons in the nuclear medium, which generate secondary excited nucleons. Comparison with experimental data is done. Cross sections for nucleon emission in coincidence with one pion are also calculated, and some coincidence observables are discussed.     

Projectivep-algebras

Presented by R. Freese.     

Canonical methods for Hamiltonian systems: numerical experiments

A Hamiltonian system possesses dynamics (e.g. preservation of volume in phase space and symplectic structure) that call for special numerical integrators, namely canonical methods. Recent research on this aspect have shown that canonical numerical integrators may be needed for Hamiltonian systems. In this paper, we focus on numerical experiments that compare canonical and non-canonical numerical integrators. Test problems are taken from different areas in physical sciences. These experiments help to buttress the claims that canonical numerical integrators give results that mimic the qualitative behavior of the original system and that canonical numerical integrators are suitable for long time integrations. Our experiments indicate that higher-order canonical methods allow for larger timestep than lower-order canonical methods.     

Topochemically Controlled Redox Reactions in the System In/Sn/Cl

The difference in reactivity of the two modifications of InCl in solid state reactions with SnCl₂ is discussed. It is explained on the basis of semi‐empirical and ab initio calculations giving the density of states diagrams of InCl and the energies of possible disproportionation reactions. Their general features are discussed on the basis of a simple bonding picture for open‐packed structures involving inert pair elements. A detailed analysis of the DOS distribution allows to pinpoint the observed redox instability to a specific structural feature of α‐InCl.     

Lagrangian Duality in Set-Valued Optimization

In this paper, we study optimization problems where the objective function and the binding constraints are set-valued maps and the solutions are defined by means of set-relations among all the images sets (Kuroiwa, D. in Takahashi, W., Tanaka, T. (eds.) Nonlinear analysis and convex analysis, pp. 221–228, 1999). We introduce a new dual problem, establish some duality theorems and obtain a Lagrangian multiplier rule of nonlinear type under convexity assumptions. A necessary condition and a sufficient condition for the existence of saddle points are given. The authors thank the two referees for valuable comments and suggestions on early versions of the paper. The research of the first author was partially supported by Ministerio de Educación y Ciencia (Spain) Project MTM2006-02629 and by Junta de Castilla y León (Spain) Project VA027B06.     

Sparse generalized Fourier transforms

Block-diagonalization of sparse equivariant discretization matrices is studied. Such matrices typically arise when partial differential equations that evolve in symmetric geometries are discretized via the finite element method or via finite differences. By considering sparse equivariant matrices as equivariant graphs, we identify a condition for when block-diagonalization via a sparse variant of a generalized Fourier transform (GFT) becomes particularly simple and fast. Characterizations for finite element triangulations of a symmetric domain are given, and formulas for assembling the block-diagonalized matrix directly are presented. It is emphasized that the GFT preserves symmetric (Hermitian) properties of an equivariant matrix. By simulating the heat equation at the surface of a sphere discretized by an icosahedral grid, it is demonstrated that the block-diagonalization is beneficial. The gain is significant for a direct method, and modest for an iterative method. A comparison with a block-diagonalization approach based upon the continuous formulation is made. It is found that the sparse GFT method is an appropriate way to discretize the resulting continuous subsystems, since the spectrum and the symmetry are preserved. AMS subject classification (2000)  43A30, 65T99, 20B25     

Mössbauer studies of the re-entrant spin-glass behaviour of Fe–Al alloys

Fe–Al alloys around the concentration of 30 at. % Al present re-entrant spin-glass behaviour at low temperatures. This behaviour is not completely understood and Mössbauer spectroscopy, combined with other experimental techniques, is useful to describe and explain this behaviour. Results show that the Mössbauer spectra coincide with the magnetic behaviour showed in literature and they can be explained as a magnetic cluster system whose magnetic clusters are getting smaller when the temperature is decreasing. When the temperature is reaching to the spin-glass transition at 92 K the spins in the paramagnetic matrix are moving slower and below this transition the spins are completely frozen.