Optimized implementation of the Lanczos method for magnetic systems

Numerically exact investigations of interacting spin systems provide a major tool for an understanding of their magnetic properties. For medium size systems the approximate Lanczos diagonalization is the most common method. In this article we suggest two improvements: efficient basis coding in subspaces and simple restructuring for openMP parallelization.Comment: 9 pages, 2 figues, submitted to Journal of Computational Physic

Model estimation and identification of manual controller objectives in complex tracking tasks

A methodology is presented for estimating the parameters in an optimal control structural model of the manual controller from experimental data on complex, multiinput/multioutput tracking tasks. Special attention is devoted to estimating the appropriate objective function for the task, as this is considered key in understanding the objectives and strategy of the manual controller. The technique is applied to data from single input/single output as well as multi input/multi outpuut experiments, and results discussed

Creation and Manipulation of Anyons in the Kitaev Model

We analyze the effect of local spin operators in the Kitaev model on the honeycomb lattice. We show, in perturbation around the isolated-dimer limit, that they create Abelian anyons together with fermionic excitations which are likely to play a role in experiments. We derive the explicit form of the operators creating and moving Abelian anyons without creating fermions and show that it involves multi-spin operations. Finally, the important experimental constraints stemming from our results are discussed.Comment: 4 pages, 3 figures, published versio

Theory of magnetization plateaux in the Shastry-Sutherland model

Using perturbative continuous unitary transformations, we determine the long-range interactions between triplets in the Shastry-Sutherland model, and we show that an unexpected structure develops at low magnetization with plateaux progressively appearing at 2/9, 1/6, 1/9 and 2/15 upon increasing the inter-dimer coupling. A critical comparison with previous approaches is included. Implications for the compound SrCu$_2$(BO$_3$)$_2$ are also discussed: we reproduce the magnetization profile around localized triplets revealed by NMR, we predict the presence of a 1/6 plateau, and we suggest that residual interactions beyond the Shastry-Sutherland model are responsible for the other plateaux below 1/3.Comment: 5 pages, 6 figure

Newtonian Limit of Conformal Gravity

We study the weak-field limit of the static spherically symmetric solution of the locally conformally invariant theory advocated in the recent past by Mannheim and Kazanas as an alternative to Einstein's General Relativity. In contrast with the previous works, we consider the physically relevant case where the scalar field that breaks conformal symmetry and generates fermion masses is nonzero. In the physical gauge, in which this scalar field is constant in space-time, the solution reproduces the weak-field limit of the Schwarzschild--(anti)DeSitter solution modified by an additional term that, depending on the sign of the Weyl term in the action, is either oscillatory or exponential as a function of the radial distance. Such behavior reflects the presence of, correspondingly, either a tachion or a massive ghost in the spectrum, which is a serious drawback of the theory under discussion.Comment: 9 pages, comments and references added; the version to be published in Phys. Rev.

Navigation and guidance analysis for a Mars mission Interim study report

Error propagation program simulating earth based tracking for navigation and guidance analysis of Mars missio

Bounding and approximating parabolas for the spectrum of Heisenberg spin systems

We prove that for a wide class of quantum spin systems with isotropic Heisenberg coupling the energy eigenvalues which belong to a total spin quantum number S have upper and lower bounds depending at most quadratically on S. The only assumption adopted is that the mean coupling strength of any spin w.r.t. its neighbours is constant for all N spins. The coefficients of the bounding parabolas are given in terms of special eigenvalues of the N times N coupling matrix which are usually easily evaluated. In addition we show that the bounding parabolas, if properly shifted, provide very good approximations of the true boundaries of the spectrum. We present numerical examples of frustrated rings, a cube, and an icosahedron.Comment: 8 pages, 3 figures. Submitted to Europhysics Letter

First-order layering and critical wetting transitions in non-additive hard sphere mixtures

Using fundamental-measure density functional theory we investigate entropic wetting in an asymmetric binary mixture of hard spheres with positive non-additivity. We consider a general planar hard wall, where preferential adsorption is induced by a difference in closest approach of the different species and the wall. Close to bulk fluid-fluid coexistence the phase rich in the minority component adsorbs either through a series of first-order layering transitions, where an increasing number of liquid layers adsorbs sequentially, or via a critical wetting transition, where a thick film grows continuously.Comment: 4 pages, 4 figure