The GENGA Code: Gravitational Encounters in N-body simulations with GPU Acceleration

We describe an open source GPU implementation of a hybrid symplectic N-body integrator, GENGA (Gravitational ENcounters with Gpu Acceleration), designed to integrate planet and planetesimal dynamics in the late stage of planet formation and stability analyses of planetary systems. GENGA uses a hybrid symplectic integrator to handle close encounters with very good energy conservation, which is essential in long-term planetary system integration. We extended the second order hybrid integration scheme to higher orders. The GENGA code supports three simulation modes: Integration of up to 2048 massive bodies, integration with up to a million test particles, or parallel integration of a large number of individual planetary systems. We compare the results of GENGA to Mercury and pkdgrav2 in respect of energy conservation and performance, and find that the energy conservation of GENGA is comparable to Mercury and around two orders of magnitude better than pkdgrav2. GENGA runs up to 30 times faster than Mercury and up to eight times faster than pkdgrav2. GENGA is written in CUDA C and runs on all NVIDIA GPUs with compute capability of at least 2.0.Comment: Accepted by ApJ. 18 pages, 17 figures, 4 table

Statistics of Lead Changes in Popularity-Driven Systems

We study statistical properties of the highest degree, or most popular, nodes in growing networks. We show that the number of lead changes increases logarithmically with network size N, independent of the details of the growth mechanism. The probability that the first node retains the lead approaches a finite constant for popularity-driven growth, and decays as N^{-phi}(ln N)^{-1/2}, with phi=0.08607..., for growth with no popularity bias.Comment: 4 pages, 4 figures, 2 column revtex format. Minor changes in response to referee comments. For publication in PR

Filament eruption connected to protospheric activity

Two cases of activation of filaments that occured in regions of intense magnetic activity was studied. The simultaneous observations from Debrecen Observatory (white light and H alpha filtergram), and from Meudon Observatory (magnetogram, MSDP dopplergram and intensity maps in H alpha) gave a complementary set of data from which can be produced evidence of the influence of the photospheric magnetic field on the destabilization process of the filaments. On June 22, 1980, the eruption of the filament is associated with the motion of pores, which are manifestations of emerging flux knots. On September 3, 1980, the twisting motions in the filament are associated to the birth of a pore in its neighborhood. These observations are discussed

Quantum Algorithm for the Collision Problem

In this note, we give a quantum algorithm that finds collisions in arbitrary r-to-one functions after only O((N/r)^(1/3)) expected evaluations of the function. Assuming the function is given by a black box, this is more efficient than the best possible classical algorithm, even allowing probabilism. We also give a similar algorithm for finding claws in pairs of functions. Furthermore, we exhibit a space-time tradeoff for our technique. Our approach uses Grover's quantum searching algorithm in a novel way.Comment: 8 pages, LaTeX2

Singular components of spectral measures for ergodic Jacobi matrices

For ergodic 1d Jacobi operators we prove that the random singular components of any spectral measure are almost surely mutually disjoint as long as one restricts to the set of positive Lyapunov exponent. In the context of extended Harper's equation this yields the first rigorous proof of the Thouless' formula for the Lyapunov exponent in the dual regions.Comment: to appear in the Journal of Mathematical Physics, vol 52 (2011