Far-Field Compression for Fast Kernel Summation Methods in High Dimensions

We consider fast kernel summations in high dimensions: given a large set of points in $d$ dimensions (with $d \gg 3$) and a pair-potential function (the {\em kernel} function), we compute a weighted sum of all pairwise kernel interactions for each point in the set. Direct summation is equivalent to a (dense) matrix-vector multiplication and scales quadratically with the number of points. Fast kernel summation algorithms reduce this cost to log-linear or linear complexity. Treecodes and Fast Multipole Methods (FMMs) deliver tremendous speedups by constructing approximate representations of interactions of points that are far from each other. In algebraic terms, these representations correspond to low-rank approximations of blocks of the overall interaction matrix. Existing approaches require an excessive number of kernel evaluations with increasing $d$ and number of points in the dataset. To address this issue, we use a randomized algebraic approach in which we first sample the rows of a block and then construct its approximate, low-rank interpolative decomposition. We examine the feasibility of this approach theoretically and experimentally. We provide a new theoretical result showing a tighter bound on the reconstruction error from uniformly sampling rows than the existing state-of-the-art. We demonstrate that our sampling approach is competitive with existing (but prohibitively expensive) methods from the literature. We also construct kernel matrices for the Laplacian, Gaussian, and polynomial kernels -- all commonly used in physics and data analysis. We explore the numerical properties of blocks of these matrices, and show that they are amenable to our approach. Depending on the data set, our randomized algorithm can successfully compute low rank approximations in high dimensions. We report results for data sets with ambient dimensions from four to 1,000.Comment: 43 pages, 21 figure

Sporadic occurrence of hemorrhagic colitis associated with Escherichia coli O157:H7 in Newfoundland

During a 9-month period in 1984, 113 fecal samples obtained from 92 patients with diarrheal illness were cultured for Escherichia coli serotype O157:H7 to determine the occurrence of this agent in diarrheal illness in Newfoundland. E. coli O157:H7 was isolated in almost pure culture from 12 stool specimens obtained from 7 (15%) of the 47 patients who had grossly bloody diarrhea; none of the 12 yielded any of the usual enteric pathogens. The agent was not isolated from the stool specimens obtained from the remaining 45 patients, who did not have bloody diarrhea. All seven patients whose specimens were positive for E. coli O157:H7 had clinical manifestations typical of hemorrhagic colitis, but the syndrome was clinically suspected and a specific test requested in only two cases. The seven cases were not clustered geographically or temporally, and no common exposure was identified. To determine whether hamburger meat was the source of E. coli O157:H7, 66 samples obtained from various retail outlets were tested; none were found to be positive. Hemorrhagic colitis may be a common disease, and E. coli O157:H7 should be considered as an agent in bloody diarrheal illness

Statistical Description of Hydrodynamic Processes in Ionic Melts with taking into account Polarization Effects

Statistical description of hydrodynamic processes for ionic melts is proposed with taking into account polarization effects caused by the deformation of external ionic shells. This description is carried out by means of the Zubarev nonequilibrium statistical operator method, appropriate for investigations of both strong and weak nonequilibrium processes. The nonequilibrium statistical operator and the generalized hydrodynamic equations that take into account polarization processes are received for ionic-polarization model of ionic molten salts when the nonequilibrium averaged values of densities of ions number, their momentum, dipole momentum and total energy are chosen for the reduced description parameters. A spectrum of collective excitations is investigated within the viscoelastic approximation for ion-polarization model of ionic melts.Comment: 24 pages, RevTex4.1-format, no figure

Sharp crossover from composite fermionization to phase separation in mesoscopic mixtures of ultracold bosons

We show that a two-component mixture of a few repulsively interacting ultracold atoms in a one-dimensional trap possesses very different quantum regimes and that the crossover between them can be induced by tuning the interactions in one of the species. In the composite fermionization regime, where the interactions between both components are large, none of the species show large occupation of any natural orbital. Our results show that by increasing the interaction in one of the species, one can reach the phase-separated regime. In this regime, the weakly interacting component stays at the center of the trap and becomes almost fully phase coherent, while the strongly interacting component is displaced to the edges of the trap. The crossover is sharp, as observed in the in the energy and the in the largest occupation of a natural orbital of the weakly interacting species. Such a transition is a purely mesoscopic effect which disappears for large atom numbers.Comment: 5 pages, 3 figure

Correlations in a two-dimensional Bose gas with long range interactions

We study the correlations of two-dimensional dipolar excitons in coupled quantum wells with a dipole -- dipole repulsive interaction. We show that at low concentrations, the Bose degeneracy of the excitons is accompanied by strong multi-particle correlations and the system behaves as a Bose liquid. At high concentration the particles interaction suppresses quantum coherence and the system behaves similar to a classical liquid down to a temperature lower than typical for a Bose gas. We evaluate the interaction energy per particle and the resulting blue shift of the exciton luminescence that is a direct tool to measure the correlations. This theory can apply to other systems of bosons with extended interaction.Comment: 11 pages including 2 figure

Quantum correlations and spatial localization in one-dimensional ultracold bosonic mixtures

We present the complete phase diagram for one-dimensional binary mixtures of bosonic ultracold atomic gases in a harmonic trap. We obtain exact results with direct numerical diagonalization for small number of atoms, which permits us to quantify quantum many-body correlations. The quantum Monte Carlo method is used to calculate energies and density profiles for larger system sizes. We study the system properties for a wide range of interaction parameters. For the extreme values of these parameters, different correlation limits can be identified, where the correlations are either weak or strong. We investigate in detail how the correlation evolve between the limits. For balanced mixtures in the number of atoms in each species, the transition between the different limits involves sophisticated changes in the one- and two-body correlations. Particularly, we quantify the entanglement between the two components by means of the von Neumann entropy. We show that the limits equally exist when the number of atoms is increased, for balanced mixtures. Also, the changes in the correlations along the transitions among these limits are qualitatively similar. We also show that, for imbalanced mixtures, the same limits with similar transitions exist. Finally, for strongly imbalanced systems, only two limits survive, i.e., a miscible limit and a phase-separated one, resembling those expected with a mean-field approach.Comment: 18 pages, 8 figure