Creation of collective many-body states and single photons from two-dimensional Rydberg lattice gases

The creation of collective many-body quantum states from a two-dimensional lattice gas of atoms is studied. Our approach relies on the van-der-Waals interaction that is present between alkali metal atoms when laser excited to high-lying Rydberg s-states. We focus on a regime in which the laser driving is strong compared to the interaction between Rydberg atoms. Here energetically low-lying many-particle states can be calculated approximately from a quadratic Hamiltonian. The potential usefulness of these states as a resource for the creation of deterministic single-photon sources is illustrated. The properties of these photon states are determined from the interplay between the particular geometry of the lattice and the interatomic spacing.Comment: 12 pages, 8 figure

A GPU-based hyperbolic SVD algorithm

A one-sided Jacobi hyperbolic singular value decomposition (HSVD) algorithm, using a massively parallel graphics processing unit (GPU), is developed. The algorithm also serves as the final stage of solving a symmetric indefinite eigenvalue problem. Numerical testing demonstrates the gains in speed and accuracy over sequential and MPI-parallelized variants of similar Jacobi-type HSVD algorithms. Finally, possibilities of hybrid CPU--GPU parallelism are discussed.Comment: Accepted for publication in BIT Numerical Mathematic

Identification and correction of systematic error in high-throughput sequence data

A feature common to all DNA sequencing technologies is the presence of base-call errors in the sequenced reads. The implications of such errors are application specific, ranging from minor informatics nuisances to major problems affecting biological inferences. Recently developed “next-gen” sequencing technologies have greatly reduced the cost of sequencing, but have been shown to be more error prone than previous technologies. Both position specific (depending on the location in the read) and sequence specific (depending on the sequence in the read) errors have been identified in Illumina and Life Technology sequencing platforms. We describe a new type of _systematic_ error that manifests as statistically unlikely accumulations of errors at specific genome (or transcriptome) locations. We characterize and describe systematic errors using overlapping paired reads form high-coverage data. We show that such errors occur in approximately 1 in 1000 base pairs, and that quality scores at systematic error sites do not account for the extent of errors. We identify motifs that are frequent at systematic error sites, and describe a classifier that distinguishes heterozygous sites from systematic error. Our classifier is designed to accommodate data from experiments in which the allele frequencies at heterozygous sites are not necessarily 0.5 (such as in the case of RNA-Seq). Systematic errors can easily be mistaken for heterozygous sites in individuals, or for SNPs in population analyses. Systematic errors are particularly problematic in low coverage experiments, or in estimates of allele-specific expression from RNA-Seq data. Our characterization of systematic error has allowed us to develop a program, called SysCall, for identifying and correcting such errors. We conclude that correction of systematic errors is important to consider in the design and interpretation of high-throughput sequencing experiments

Spectral asymmetry and Riemannian geometry. III

In Parts I and II of this paper ((4),(5)) we studied the 'spectral asymmetry' of certain elliptic self-adjoint operators arising in Riemannian geometry. More precisely, for any elliptic self-adjoint operator A on a compact manifold we defined ηA(s)=Σλ+0signλ|λ|-8, where λ runs over the eigenvalues of A. For the particular operators of interest in Riemannian geometry we showed that ηA(s) had an analytic continuation to the whole complex s-plane, with simple poles, and that s=0 was not a pole. The real number ηA(0), which is a measure of 'spectral asymmetry', was studied in detail particularly in relation to representations of the fundamental group

Novel Modifications of Parallel Jacobi Algorithms

We describe two main classes of one-sided trigonometric and hyperbolic Jacobi-type algorithms for computing eigenvalues and eigenvectors of Hermitian matrices. These types of algorithms exhibit significant advantages over many other eigenvalue algorithms. If the matrices permit, both types of algorithms compute the eigenvalues and eigenvectors with high relative accuracy. We present novel parallelization techniques for both trigonometric and hyperbolic classes of algorithms, as well as some new ideas on how pivoting in each cycle of the algorithm can improve the speed of the parallel one-sided algorithms. These parallelization approaches are applicable to both distributed-memory and shared-memory machines. The numerical testing performed indicates that the hyperbolic algorithms may be superior to the trigonometric ones, although, in theory, the latter seem more natural.Comment: Accepted for publication in Numerical Algorithm

New representation of orbital motion with arbitrary angular momenta

A new formulation is presented for a variational calculation of $N$-body systems on a correlated Gaussian basis with arbitrary angular momenta. The rotational motion of the system is described with a single spherical harmonic of the total angular momentum $L$, and thereby needs no explicit coupling of partial waves between particles. A simple generating function for the correlated Gaussian is exploited to derive the matrix elements. The formulation is applied to various Coulomb three-body systems such as $e^-e^-e^+, tt\mu, td\mu$, and $\alpha e^-e^-$ up to $L=4$ in order to show its usefulness and versatility. A stochastic selection of the basis functions gives good results for various angular momentum states.Comment: Revte

Intensity interferometry of single x-ray pulses from a synchrotron storage ring

We report on measurements of second-order intensity correlations at the high brilliance storage ring PETRA III using a prototype of the newly developed Adaptive Gain Integrating Pixel Detector (AGIPD). The detector recorded individual synchrotron radiation pulses with an x-ray photon energy of 14.4 keV and repetition rate of about 5 MHz. The second-order intensity correlation function was measured simultaneously at different spatial separations that allowed to determine the transverse coherence length at these x-ray energies. The measured values are in a good agreement with theoretical simulations based on the Gaussian Schell-model.Comment: 16 pages, 6 figures, 42 reference

A note on monopole moduli spaces

We discuss the structure of the framed moduli space of Bogomolny monopoles for arbitrary symmetry breaking and extend the definition of its stratification to the case of arbitrary compact Lie groups. We show that each stratum is a union of submanifolds for which we conjecture that the natural $L^2$ metric is hyperKahler. The dimensions of the strata and of these submanifolds are calculated, and it is found that for the latter, the dimension is always a multiple of four.Comment: 17 pages, LaTe

Universal time-evolution of a Rydberg lattice gas with perfect blockade

We investigate the dynamics of a strongly interacting spin system that is motivated by current experimental realizations of strongly interacting Rydberg gases in lattices. In particular we are interested in the temporal evolution of quantities such as the density of Rydberg atoms and density-density correlations when the system is initialized in a fully polarized state without Rydberg excitations. We show that in the thermodynamic limit the expectation values of these observables converge at least logarithmically to universal functions and outline a method to obtain these functions. We prove that a finite one-dimensional system follows this universal behavior up to a given time. The length of this universal time period depends on the actual system size. This shows that already the study of small systems allows to make precise predictions about the thermodynamic limit provided that the observation time is sufficiently short. We discuss this for various observables and for systems with different dimensions, interaction ranges and boundary conditions.Comment: 16 pages, 3 figure