Computationally efficient algorithms for the two-dimensional Kolmogorov-Smirnov test

Goodness-of-fit statistics measure the compatibility of random samples against some theoretical or reference probability distribution function. The classical one-dimensional Kolmogorov-Smirnov test is a non-parametric statistic for comparing two empirical distributions which defines the largest absolute difference between the two cumulative distribution functions as a measure of disagreement. Adapting this test to more than one dimension is a challenge because there are 2^d-1 independent ways of ordering a cumulative distribution function in d dimensions. We discuss Peacock's version of the Kolmogorov-Smirnov test for two-dimensional data sets which computes the differences between cumulative distribution functions in 4n^2 quadrants. We also examine Fasano and Franceschini's variation of Peacock's test, Cooke's algorithm for Peacock's test, and ROOT's version of the two-dimensional Kolmogorov-Smirnov test. We establish a lower-bound limit on the work for computing Peacock's test of Omega(n^2.lg(n)), introducing optimal algorithms for both this and Fasano and Franceschini's test, and show that Cooke's algorithm is not a faithful implementation of Peacock's test. We also discuss and evaluate parallel algorithms for Peacock's test

Exploring the inner region of Type 1 AGNs with the Keck interferometer

The exploration of extragalactic objects with long-baseline interferometers in the near-infrared has been very limited. Here we report successful observations with the Keck interferometer at K-band (2.2 um) for four Type 1 AGNs, namely NGC4151, Mrk231, NGC4051, and the QSO IRAS13349+2438 at z=0.108. For the latter three objects, these are the first long-baseline interferometric measurements in the infrared. We detect high visibilities (V^2 ~ 0.8-0.9) for all the four objects, including NGC4151 for which we confirm the high V^2 level measured by Swain et al.(2003). We marginally detect a decrease of V^2 with increasing baseline lengths for NGC4151, although over a very limited range, where the decrease and absolute V^2 are well fitted with a ring model of radius 0.45+/-0.04 mas (0.039+/-0.003 pc). Strikingly, this matches independent radius measurements from optical--infrared reverberations that are thought to be probing the dust sublimation radius. We also show that the effective radius of the other objects, obtained from the same ring model, is either roughly equal to or slightly larger than the reverberation radius as a function of AGN luminosity. This suggests that we are indeed partially resolving the dust sublimation region. The ratio of the effective ring radius to the reverberation radius might also give us an approximate probe for the radial structure of the inner accreting material in each object. This should be scrutinized with further observations.Comment: accepted for publication in A&A Letter

Quantum Coherence of Relic Neutrinos

We argue that in at least a portion of the history of the universe the relic background neutrinos are spatially-extended, coherent superpositions of mass states. We show that an appropriate quantum mechanical treatment affects the neutrino mass values derived from cosmological data. The coherence scale of these neutrino flavor wavepackets can be an appreciable fraction of the causal horizon size, raising the possibility of spacetime curvature-induced decoherence.Comment: 4 pages, 4 figures; matches publication in PR

Mapping the radial structure of AGN tori

We present mid-IR interferometric observations of 6 type 1 AGNs at multiple baseline lengths of 27--130m, reaching high angular resolutions up to lambda/B~0.02 arcseconds. For two of the targets, we have simultaneous near-IR interferometric measurements as well. The multiple baseline data directly probe the radial distribution of the material on sub-pc scales. Within our sample, which is small but spans over ~2.5 orders of magnitudes in the UV/optical luminosity L of the central engine, the radial distribution clearly and systematically changes with luminosity. First, we show that the brightness distribution at a given mid-IR wavelength seems to be rather well described by a power law, which makes a simple Gaussian or ring size estimation quite inadequate. Here we instead use a half-light radius R_1/2 as a representative size. We then find that the higher luminosity objects become more compact in normalized half-light radii R_1/2 /R_in in the mid-IR, where R_in is the dust sublimation radius empirically given by the L^1/2 fit of the near-IR reverberation radii. This means that, contrary to previous studies, the physical mid-IR emission size (e.g. in pc) is not proportional to L^1/2, but increases with L much more slowly, or in fact, nearly constant at 13 micron. Combining the size information with the total flux specta, we infer that the radial surface density distribution of the heated dust grains changes from a steep ~r^-1 structure in high luminosity objects to a shallower ~r^0 structure in those of lower luminosity. The inward dust temperature distribution does not seem to smoothly reach the sublimation temperature -- on the innermost scale of ~R_in, a relatively low temperature core seems to co-exist with a slightly distinct brightness concentration emitting roughly at the sublimation temperature.Comment: accepted for publication in A&

In--Flight (K−,pK^-,p) Reactions for the Formation of Kaonic Atoms and Kaonic Nuclei in Green function method

We study theoretically the kaonic atom and kaonic nucleus formations in the in--flight ($K^-,p$) reactions using the Green function method, which is suited to evaluate formation rates both of stable and unstable bound systems. We consider $^{12}$C and $^{16}$O as the targets and calculate the spectra of the ($K^-,p$) reactions. We conclude that a no peak structure due to kaonic nucleus formation is expected in the reaction spectra calculated with the chiral unitary kaon--nucleus optical potential. In the spectra with the phenomenological deep kaon--nucleus potential, we may have possibilities to observe some structures due to kaonic nucleus states. For all cases, we have peaks due to the kaonic atom formations in the reaction spectra.Comment: 10 pages, 9 figures, newly calculated results added, revisions and updated references, to appear in Physical Review

Connection between rotation and miscibility in a two-component Bose-Einstein condensate

A two-component Bose-Einstein condensate rotating in a toroidal trap is investigated. The topological constraint depends on the density distribution of each component along the circumference of the torus, and therefore the quantization condition on the circulation can be controlled by changing the miscibility using the Feshbach resonance. We find that the system exhibits a variety of dynamics depending on the initial angular momentum when the miscibility is changed.Comment: 6 pages, 7 figure

Formation of Deeply Bound Kaonic Atoms in (K^-,N) Reactions

We study theoretically the (K^-,N) reactions for the formation of the deeply bound kaonic atoms, which were predicted to be quasi--stable with narrow widths, using the Green function method. We consider various cases with different target nuclei and energies systematically and find the clear signals in the theoretical spectra for all cases considered in this article. The signals show very interesting structures, such as the $RESONANCE DIP$ instead of the resonance peak. We discuss the origins of the interesting structures and possibilities to get new information on the existence of the kaonic nuclei from the spectra of the atomic state formations.Comment: 11 pages, 9 figure

Self-Consistent Velocity Dependent Effective Interactions

The theory of self-consistent effective interactions in nuclei is extended for a system with a velocity dependent mean potential. By means of the field coupling method, we present a general prescription to derive effective interactions which are consistent with the mean potential. For a deformed system with the conventional pairing field, the velocity dependent effective interactions are derived as the multipole pairing interactions in doubly-stretched coordinates. They are applied to the microscopic analysis of the giant dipole resonances (GDR's) of ${}^{148,154}Sm$, the first excited $2^+$ states of Sn isotopes and the first excited $3^-$ states of Mo isotopes. It is clarified that the interactions play crucial roles in describing the splitting and structure of GDR peaks, in restoring the energy weighted sum rule, and in reducing the values of $B(E\lambda)$.Comment: 35 pages, RevTeX, 7 figures (available upon request), to appear in Phys.Rev.

A dust-parallax distance of 19 megaparsecs to the supermassive black hole in NGC 4151

The active galaxy NGC 4151 has a crucial role as one of only two active galactic nuclei for which black hole mass measurements based on emission line reverberation mapping can be calibrated against other dynamical methods. Unfortunately, effective calibration requires an accurate distance to NGC 4151, which is currently not available. Recently reported distances range from 4 to 29 megaparsecs (Mpc). Strong peculiar motions make a redshift-based distance very uncertain, and the geometry of the galaxy and its nucleus prohibit accurate measurements using other techniques. Here we report a dust-parallax distance to NGC 4151 of $D_A = 19.0^{+2.4}_{-2.6}$ Mpc. The measurement is based on an adaptation of a geometric method proposed previously using the emission line regions of active galaxies. Since this region is too small for current imaging capabilities, we use instead the ratio of the physical-to-angular sizes of the more extended hot dust emission as determined from time-delays and infrared interferometry. This new distance leads to an approximately 1.4-fold increase in the dynamical black hole mass, implying a corresponding correction to emission line reverberation masses of black holes if they are calibrated against the two objects with additional dynamical masses.Comment: Authors' version of a letter published in Nature (27 November 2014); 8 pages, 5 figures, 1 tabl