Computing the lowest eigenvalues of the Fermion matrix by subspace iterations

Subspace iterations are used to minimise a generalised Ritz functional of a large, sparse Hermitean matrix. In this way, the lowest $m$ eigenvalues are determined. Tests with $1 \leq m \leq 32$ demonstrate that the computational cost (no. of matrix multiplies) does not increase substantially with $m$. This implies that, as compared to the case of a $m=1$, the additional eigenvalues are obtained for free.Comment: Talk presented at LATTICE96(algorithms), 3 pages, 2 Postscript figures, uses epsf.sty, espcrc2.st

The Height of a Giraffe

A minor modification of the arguments of Press and Lightman leads to an estimate of the height of the tallest running, breathing organism on a habitable planet as the Bohr radius multiplied by the three-tenths power of the ratio of the electrical to gravitational forces between two protons (rather than the one-quarter power that Press got for the largest animal that would not break in falling over, after making an assumption of unreasonable brittleness). My new estimate gives a height of about 3.6 meters rather than Press's original estimate of about 2.6 cm. It also implies that the number of atoms in the tallest runner is very roughly of the order of the nine-tenths power of the ratio of the electrical to gravitational forces between two protons, which is about 3 x 10^32.Comment: 12 pages, LaTe

Dipole Excitation of Dipositronium

The energy interval between the ground and the P-wave excited states of the recently discovered positronium molecule Ps_2 is evaluated, including the relativistic and the leading logarithmic radiative corrections, E_P-E_S = 0.181 586 7(8) a.u. The P-state, decaying usually via annihilation, is found to decay into the ground state by an electric dipole transition 19 percent of the time. Anticipated observation of this transition will provide insight into this exotic system.Comment: 5 page

Ab initio mass tensor molecular dynamics

Mass tensor molecular dynamics was first introduced by Bennett [J. Comput. Phys. 19, 267 (1975)] for efficient sampling of phase space through the use of generalized atomic masses. Here, we show how to apply this method to ab initio molecular dynamics simulations with minimal computational overhead. Test calculations on liquid water show a threefold reduction in computational effort without making the fixed geometry approximation. We also present a simple recipe for estimating the optimal atomic masses using only the first derivatives of the potential energy.Comment: 19 pages, 5 figure

Effects of Line-tying on Resistive Tearing Instability in Slab Geometry

The effects of line-tying on resistive tearing instability in slab geometry is studied within the framework of reduced magnetohydrodynamics (RMHD).\citep{KadomtsevP1974,Strauss1976} It is found that line-tying has a stabilizing effect. The tearing mode is stabilized when the system length $L$ is shorter than a critical length $L_{c}$, which is independent of the resistivity $\eta$. When $L$ is not too much longer than $L_{c}$, the growthrate $\gamma$ is proportional to $\eta$ . When $L$ is sufficiently long, the tearing mode scaling $\gamma\sim\eta^{3/5}$ is recovered. The transition from $\gamma\sim\eta$ to $\gamma\sim\eta^{3/5}$ occurs at a transition length $L_{t}\sim\eta^{-2/5}$.Comment: Correct a typ

Relativistic time dilatation and the spectrum of electrons emitted by 33 TeV lead ions penetrating thin foils

We study the energy distribution of ultrarelativistic electrons produced when a beam of 33 TeV Pb$^{81+}$(1s) ions penetrates a thin Al foil. We show that, because of a prominent role of the excitations of the ions inside the foil which becomes possible due to the relativistic time dilatation, the width of this distribution can be much narrower compared to the case when the ions interact with rarefied gaseous targets. We also show that a very similar shape of the energy distribution may arise when 33 TeV Pb$^{82+}$ ions penetrate a thin Au foil. These results shed some light on the origin of the very narrow electron energy distributions observed experimentally about a decade ago.Comment: Four pages, two figure

The locality of the square-root method for improved staggered quarks

We study the effects of improvement on the locality of square-rooted staggered Dirac operators in lattice QCD simulations. We find the localisation lengths of the improved operators (FAT7TAD and ASQTAD) to be very similar to that of the one-link operator studied by Bunk et al., being at least the Compton wavelength of the lightest particle in the theory, even in the continuum limit. We conclude that improvement has no effect. We discuss the implications of this result for the locality of the nth-rooted fermion determinant used to reduce the number of sea quark flavours, and for possible staggered valence quark formulations

Quantum dynamics of an Ising spin-chain in a random transverse field

We consider an Ising spin-chain in a random transverse magnetic field and compute the zero temperature wave vector and frequency dependent dynamic structure factor numerically by using Jordan-Wigner transformation. Two types of distributions of magnetic fields are introduced. For a rectangular distribution, a dispersing branch is observed, and disorder tends to broaden the dispersion peak and close the excitation gap. For a binary distribution, a non-dispersing branch at almost zero energy is recovered. We discuss the relationship of our work to the neutron scattering measurement in $\mathrm{LiHoF_4}$.Comment: 4 pages and 6 eps figures; minor clarifications were made; the text was shortened to add an additional figur

Magnification Ratio of the Fluctuating Light in Gravitational Lens 0957+561

Radio observations establish the B/A magnification ratio of gravitational lens 0957+561 at about 0.75. Yet, for more than 15 years, the optical magnfication ratio has been between 0.9 and 1.12. The accepted explanation is microlensing of the optical source. However, this explanation is mildly discordant with (i) the relative constancy of the optical ratio, and (ii) recent data indicating possible non-achromaticity in the ratio. To study these issues, we develop a statistical formalism for separately measuring, in a unified manner, the magnification ratio of the fluctuating and constant parts of the light curve. Applying the formalism to the published data of Kundi\'c et al. (1997), we find that the magnification ratios of fluctuating parts in both the g and r colors agrees with the magnification ratio of the constant part in g-band, and tends to disagree with the r-band value. One explanation could be about 0.1 mag of consistently unsubtracted r light from the lensing galaxy G1, which seems unlikely. Another could be that 0957+561 is approaching a caustic in the microlensing pattern.Comment: 12 pages including 1 PostScript figur

Stable resonances and signal propagation in a chaotic network of coupled units

We apply the linear response theory developed in \cite{Ruelle} to analyze how a periodic signal of weak amplitude, superimposed upon a chaotic background, is transmitted in a network of non linearly interacting units. We numerically compute the complex susceptibility and show the existence of specific poles (stable resonances) corresponding to the response to perturbations transverse to the attractor. Contrary to the poles of correlation functions they depend on the pair emitting/receiving units. This dynamic differentiation, induced by non linearities, exhibits the different ability that units have to transmit a signal in this network.Comment: 10 pages, 3 figures, to appear in Phys. rev.