Energy levels of small titanium oxide clusters obtained from SCF calculations

The energy levels of small titanium oxide clusters [(TiO₂)₂, and (TiO₂H)₃, and (TiO₂H)₂] have been calculated using ab initio SCF methods. Both crystal and relaxed geometries have been considered. Systematic changes in the valencelevel structure resulting from geometry relaxation are found, which may be related to band-gap changes experimentally observed for small titanium oxide particles. In addition, a Ti? OH local surface state is found to be well described within a limited cluster model. © 1994 John Wiley & Sons, Inc.     

A semi-analytical numerical method for time-dependent radiative transfer problems in slab geometry with coherent isotropic scattering

We describe a semi-analytical numerical method for coherent isotropic scattering time-dependent radiative transfer problems in slab geometry. This numerical method is based on a combination of two classes of numerical methods: the spectral methods and the Laplace transform (LTS_N) methods applied to the radiative transfer equation in the discrete ordinates (S_N) formulation. The basic idea is to use the essence of the spectral methods and expand the intensity of radiation in a truncated series of Laguerre polynomials in the time variable and then solve recursively the resulting set of “time-independent” S_N problems by using the LTS_N method. We show some numerical experiments for a typical model problem.     

Study of the charge-exchange reactions\pi ^ - p \to (\pi ^0 ,\eta ,\eta \prime )n at 63 GeV

Muon electron pairs were detected in an Al multiplate spark chamber, exposed to a neutrino beam from the CERN PS. The leptons were not accompanied by other particles, except occasionally by protons. The background came mainly from muon associated π⁰ production, with one decay gamma lost. It was determined empirically, together with the small contribution from υ _e reactions. For electron energies above 2 GeV the background is 5.7±1.5 events, whereas 18 (μe)-candidates have been observed. Hence the effect is established, with a rate of about 10^?4 as compared to the muonic reactions above 3 GeV. Charm creation as the origin of this (μe)-production process is excluded; heavy neutral lepton production does not fit the kinematics observed. Instead the events are compatible with the two-body decay of an object with variable invariant mass of order 1 GeV, possibly resulting from axion interactions.     

An algebraic characterization ofB-convergent Runge-Kutta methods

Summary In the analysis of discretization methods for stiff intial value problems, stability questions have received most part of the attention in the past.B-stability and the equivalent criterion algebraic stability are well known concepts for Runge-Kutta methods applied to dissipative problems. However, for the derivation ofB-convergence results — error bounds which are not affected by stiffness — it is not sufficient in many cases to requireB-stability alone. In this paper, necessary and sufficient conditions forB-convergence are determined.This paper was written while J. Schneid was visiting the Centre for Mathematics and Computer Science with an Erwin-Schrödinger stipend from the Fonds zur Förderung der wissenschaftlichen Forschung     

Mass Power Spectrum in a Universe Dominated by the Chaplygin Gas

The mass power spectrum for a Universe dominated by the Chaplygin gas is evaluated numerically from scales of the order of the Hubble horizon to 100 Mpc. The results are compared with a pure baryonic Universe and a cosmological constant model. In all three cases, the spectrum increases with k, the wavenumber of the perturbations. The slope of the spectrum is higher for the baryonic model and smaller for the cosmological constant model, the Chaplygin gas interpolating these two models. The results are analyzed in terms of the sound velocity of the Chaplygin gas and the moment the Universe begins to accelerate.     

Stability of a Model of Relativistic Quantum Electrodynamics

The relativistic “no pair” model of quantum electrodynamics uses the Dirac operator, D(A) for the electron dynamics together with the usual self-energy of the quantized ultraviolet cutoff electromagnetic field A– in the Coulomb gauge. There are no positrons because the electron wave functions are constrained to lie in the positive spectral subspace of some Dirac operator, D, but the model is defined for any number, N, of electrons, and hence describes a true many-body system. In addition to the electrons there are a number, K, of fixed nuclei with charges ≤Z. If the fields are not quantized but are classical, it was shown earlier that such a model is always unstable (the ground state energy E=−∞) if one uses the customary D(0) to define the electron space, but is stable (E > − const.(N+K)) if one uses D(A) itself (provided the fine structure constant α and Z are not too large). This result is extended to quantized fields here, and stability is proved for α= 1/137 and Z≤ 42. This formulation of QED is somewhat unusual because it means that the electron Hilbert space is inextricably linked to the photon Fock space. But such a linkage appears to better describe the real world of photons and electrons. Received: 8 September 2001 / Accepted: 18 March 2002     

VIBRATION OF TWISTED AND CURVED CYLINDRICAL PANELS WITH VARIABLE THICKNESS

A numerical procedure, with an exact strain-displacement relationship of twisted and curved cylindrical panels having variable thickness derived by considering the Green strain tensor on general shell theory, is presented using the principle of virtual work and the Rayleigh-Ritz method with algebraic polynomials as in-plane and transverse displacement functions. The accuracy and applicability of the procedure are verified by comparing the present results with previous experimental and theoretical results for several panels. The effects of variation ratio of thickness in chordwise and lengthwise directions, twist, and curvature both in two directions aforementioned on vibrations of cylindrical panels are studied in detail, and typical vibration mode shapes are plotted to demonstrate the effects.     

Ermakov Approach for Minisuperspace Oscillators

The Wheeler–DeWitt equation of arbitrary Hartle–Hawking factor ordering for several minisuperspace universe models, such as the pure gravity Friedmann–Robertson–Walker and Taub ones, is mapped onto the dynamics of corresponding classical oscillators. The latter ones are studied by the classical Ermakov invariant method, which is a natural approach in this context. For the more realistic case of a minimally coupled massive scalar field, one can study, within the same type of approach, the corresponding squeezing features as a possible means of describing cosmological evolution. Finally, we comment on the analogy with the accelerator physics.