Strains and Jets in Black Hole Fields

We study the behaviour of an initially spherical bunch of particles emitted along trajectories parallel to the symmetry axis of a Kerr black hole. We show that, under suitable conditions, curvature and inertial strains compete to generate jet-like structures.Comment: To appear in the Proceedings of the Spanish Relativity Meeting 2007 held in Tenerife (Spain) 3 Figure

Quotients of the Dwork pencil

In this paper we investigate the geometry of the Dwork pencil in any dimension. More specifically, we study the automorphism group G of the generic fiber of the pencil over the complex projective line, and the quotients of it by various subgroups of G. In particular, we compute the Hodge numbers of these quotients via orbifold cohomology

Quasiseparable Hessenberg reduction of real diagonal plus low rank matrices and applications

We present a novel algorithm to perform the Hessenberg reduction of an $n\times n$ matrix $A$ of the form $A = D + UV^*$ where $D$ is diagonal with real entries and $U$ and $V$ are $n\times k$ matrices with $k\le n$. The algorithm has a cost of $O(n^2k)$ arithmetic operations and is based on the quasiseparable matrix technology. Applications are shown to solving polynomial eigenvalue problems and some numerical experiments are reported in order to analyze the stability of the approac

Solving polynomial eigenvalue problems by means of the Ehrlich-Aberth method

Given the $n\times n$ matrix polynomial $P(x)=\sum_{i=0}^kP_i x^i$, we consider the associated polynomial eigenvalue problem. This problem, viewed in terms of computing the roots of the scalar polynomial $\det P(x)$, is treated in polynomial form rather than in matrix form by means of the Ehrlich-Aberth iteration. The main computational issues are discussed, namely, the choice of the starting approximations needed to start the Ehrlich-Aberth iteration, the computation of the Newton correction, the halting criterion, and the treatment of eigenvalues at infinity. We arrive at an effective implementation which provides more accurate approximations to the eigenvalues with respect to the methods based on the QZ algorithm. The case of polynomials having special structures, like palindromic, Hamiltonian, symplectic, etc., where the eigenvalues have special symmetries in the complex plane, is considered. A general way to adapt the Ehrlich-Aberth iteration to structured matrix polynomial is introduced. Numerical experiments which confirm the effectiveness of this approach are reported.Comment: Submitted to Linear Algebra App

Spinning test particles and clock effect in Kerr spacetime

We study the motion of spinning test particles in Kerr spacetime using the Mathisson-Papapetrou equations; we impose different supplementary conditions among the well known Corinaldesi-Papapetrou, Pirani and Tulczyjew's and analyze their physical implications in order to decide which is the most natural to use. We find that if the particle's center of mass world line, namely the one chosen for the multipole reduction, is a spatially circular orbit (sustained by the tidal forces due to the spin) then the generalized momentum $P$ of the test particle is also tangent to a spatially circular orbit intersecting the center of mass line at a point. There exists one such orbit for each point of the center of mass line where they intersect; although fictitious, these orbits are essential to define the properties of the spinning particle along its physical motion. In the small spin limit, the particle's orbit is almost a geodesic and the difference of its angular velocity with respect to the geodesic value can be of arbitrary sign, corresponding to the spin-up and spin-down possible alignment along the z-axis. We also find that the choice of the supplementary conditions leads to clock effects of substantially different magnitude. In fact, for co-rotating and counter-rotating particles having the same spin magnitude and orientation, the gravitomagnetic clock effect induced by the background metric can be magnified or inhibited and even suppressed by the contribution of the individual particle's spin. Quite surprisingly this contribution can be itself made vanishing leading to a clock effect undistiguishable from that of non spinning particles. The results of our analysis can be observationally tested.Comment: IOP macros, eps figures n. 12, to appear on Classical and Quantum Gravity, 200

Kerr metric, static observers and Fermi coordinates

The coordinate transformation which maps the Kerr metric written in standard Boyer-Lindquist coordinates to its corresponding form adapted to the natural local coordinates of an observer at rest at a fixed position in the equatorial plane, i.e., Fermi coordinates for the neighborhood of a static observer world line, is derived and discussed in a way which extends to any uniformly circularly orbiting observer there.Comment: 15 page latex iopart class documen

On Functions of quasi Toeplitz matrices

Let $a(z)=\sum_{i\in\mathbb Z}a_iz^i$ be a complex valued continuous function, defined for $|z|=1$, such that $\sum_{i=-\infty}^{+\infty}|ia_i|<\infty$. Consider the semi-infinite Toeplitz matrix $T(a)=(t_{i,j})_{i,j\in\mathbb Z^+}$ associated with the symbol $a(z)$ such that $t_{i,j}=a_{j-i}$. A quasi-Toeplitz matrix associated with the continuous symbol $a(z)$ is a matrix of the form $A=T(a)+E$ where $E=(e_{i,j})$, $\sum_{i,j\in\mathbb Z^+}|e_{i,j}|<\infty$, and is called a CQT-matrix. Given a function $f(x)$ and a CQT matrix $M$, we provide conditions under which $f(M)$ is well defined and is a CQT matrix. Moreover, we introduce a parametrization of CQT matrices and algorithms for the computation of $f(M)$. We treat the case where $f(x)$ is assigned in terms of power series and the case where $f(x)$ is defined in terms of a Cauchy integral. This analysis is applied also to finite matrices which can be written as the sum of a Toeplitz matrix and of a low rank correction