The Right Mutation Strength for Multi-Valued Decision Variables

The most common representation in evolutionary computation are bit strings. This is ideal to model binary decision variables, but less useful for variables taking more values. With very little theoretical work existing on how to use evolutionary algorithms for such optimization problems, we study the run time of simple evolutionary algorithms on some OneMax-like functions defined over $\Omega = \{0, 1, \dots, r-1\}^n$. More precisely, we regard a variety of problem classes requesting the component-wise minimization of the distance to an unknown target vector $z \in \Omega$. For such problems we see a crucial difference in how we extend the standard-bit mutation operator to these multi-valued domains. While it is natural to select each position of the solution vector to be changed independently with probability $1/n$, there are various ways to then change such a position. If we change each selected position to a random value different from the original one, we obtain an expected run time of $\Theta(nr \log n)$. If we change each selected position by either $+1$ or $-1$ (random choice), the optimization time reduces to $\Theta(nr + n\log n)$. If we use a random mutation strength $i \in \{0,1,\ldots,r-1\}^n$ with probability inversely proportional to $i$ and change the selected position by either $+i$ or $-i$ (random choice), then the optimization time becomes $\Theta(n \log(r)(\log(n)+\log(r)))$, bringing down the dependence on $r$ from linear to polylogarithmic. One of our results depends on a new variant of the lower bounding multiplicative drift theorem.Comment: an extended abstract of this work is to appear at GECCO 201

A system for online beam emittance measurements and proton beam characterization

A system for online measurement of the transverse beam emittance was developed. It is named $^{4}$PrOB$\varepsilon$aM (4-Profiler Online Beam Emittance Measurement) and was conceived to measure the emittance in a fast and efficient way using the multiple beam profiler method. The core of the system is constituted by four consecutive UniBEaM profilers, which are based on silica fibers passing across the beam. The $^{4}$PrOB$\varepsilon$aM system was deployed for characterization studies of the 18~MeV proton beam produced by the IBA Cyclone 18 MeV cyclotron at Bern University Hospital (Inselspital). The machine serves daily radioisotope production and multi-disciplinary research, which is carried out with a specifically conceived Beam Transport Line (BTL). The transverse RMS beam emittance of the cyclotron was measured as a function of several machine parameters, such as the magnetic field, RF peak voltage, and azimuthal angle of the stripper. The beam emittance was also measured using the method based on the quadrupole strength variation. The results obtained with both techniques were compared and a good agreement was found. In order to characterize the longitudinal dynamics, the proton energy distribution was measured. For this purpose, a method was developed based on aluminum absorbers of different thicknesses, a UniBEaM detector, and a Faraday cup. The results were an input for a simulation of the BTL developed in the MAD-X software. This tool allows machine parameters to be tuned online and the beam characteristics to be optimized for specific applications.Comment: published in Journal of Instrumentatio

Inference of the Cold Dark Matter substructure mass function at z=0.2 using strong gravitational lenses

We present the results of a search for galaxy substructures in a sample of 11 gravitational lens galaxies from the Sloan Lens ACS Survey. We find no significant detection of mass clumps, except for a luminous satellite in the system SDSS J0956+5110. We use these non-detections, in combination with a previous detection in the system SDSS J0946+1006, to derive constraints on the substructure mass function in massive early-type host galaxies with an average redshift z ~ 0.2 and an average velocity dispersion of 270 km/s. We perform a Bayesian inference on the substructure mass function, within a median region of about 32 kpc squared around the Einstein radius (~4.2 kpc). We infer a mean projected substructure mass fraction $f = 0.0076^{+0.0208}_{-0.0052}$ at the 68 percent confidence level and a substructure mass function slope $\alpha$ < 2.93 at the 95 percent confidence level for a uniform prior probability density on alpha. For a Gaussian prior based on Cold Dark Matter (CDM) simulations, we infer $f = 0 .0064^{+0.0080}_{-0.0042}$ and a slope of $\alpha$ = 1.90$^{+0.098}_{-0.098}$ at the 68 percent confidence level. Since only one substructure was detected in the full sample, we have little information on the mass function slope, which is therefore poorly constrained (i.e. the Bayes factor shows no positive preference for any of the two models).The inferred fraction is consistent with the expectations from CDM simulations and with inference from flux ratio anomalies at the 68 percent confidence level.Comment: Accepted for publication on MNRAS, some typos corrected and some important references adde

Luminous Satellites II: Spatial Distribution, Luminosity Function and Cosmic Evolution

We infer the normalization and the radial and angular distributions of the number density of satellites of massive galaxies ($\log_{10}[M_{h}^*/M\odot]>10.5$) between redshifts 0.1 and 0.8 as a function of host stellar mass, redshift, morphology and satellite luminosity. Exploiting the depth and resolution of the COSMOS HST images, we detect satellites up to eight magnitudes fainter than the host galaxies and as close as 0.3 (1.4) arcseconds (kpc). Describing the number density profile of satellite galaxies to be a projected power law such that P(R)\propto R^{\rpower}, we find \rpower=-1.1\pm 0.3. We find no dependency of \rpower on host stellar mass, redshift, morphology or satellite luminosity. Satellites of early-type hosts have angular distributions that are more flattened than the host light profile and are aligned with its major axis. No significant average alignment is detected for satellites of late-type hosts. The number of satellites within a fixed magnitude contrast from a host galaxy is dependent on its stellar mass, with more massive galaxies hosting significantly more satellites. Furthermore, high-mass late-type hosts have significantly fewer satellites than early-type galaxies of the same stellar mass, likely a result of environmental differences. No significant evolution in the number of satellites per host is detected. The cumulative luminosity function of satellites is qualitatively in good agreement with that predicted using subhalo abundance matching techniques. However, there are significant residual discrepancies in the absolute normalization, suggesting that properties other than the host galaxy luminosity or stellar mass determine the number of satellites.Comment: 23 pages, 12 figures, Accepted for publication in the Astrophysical Journa

Light Nuclei solving Auger puzzles. The Cen-A imprint

Ultra High Energy Cosmic Rays (UHECR) map at 60 EeV have been found recently by AUGER group spreading anisotropy signatures in the sky. The result have been interpreted as a manifestation of AGN sources ejecting protons at GZK edges mostly from Super-galactic Plane. The result is surprising due to the absence of much nearer Virgo cluster. Moreover, early GZK cut off in the spectra may be better reconcile with light nuclei (than with protons). In addition a large group (nearly a dozen) of events cluster suspiciously along Cen-A. Finally, proton UHECR composition nature is in sharp disagreement with earlier AUGER claim of a heavy nuclei dominance at 40 EeV. Therefore we interpret here the signals as mostly UHECR light nuclei (He, Be, B, C, O), very possibly mostly the lightest (He,Be) ones, ejected from nearest AGN Cen-A, UHECR smeared by galactic magnetic fields, whose random vertical bending is overlapping with super-galactic arm. The eventual AUGER misunderstanding took place because of such a rare coincidence between the Super Galactic Plane (arm) and the smeared (randomized) signals from Cen-A, bent orthogonally to the Galactic fields. Our derivation verify the consistence of the random smearing angles for He, Be and B, C, O, in reasonable agreement with the AUGER main group events around Cen-A. Only few other rare events are spread elsewhere. The most collimated from Cen-A are the lightest. The most spread the heavier. Consequently Cen-A is the best candidate UHE neutrino tau observable by HEAT and AMIGA as enhanced AUGER array at tens-hundred PeV energy. This model maybe soon tested by new events clustering around the Cen-A and by composition imprint study.Comment: 4 pages, 5 figures

Observations of Ultra-High Energy Cosmic Rays

The status of measurements of the arrival directions, mass composition and energy spectrum of cosmic rays above 3 x 10^18 eV (3 EeV) is reviewed using reports presented at the 29th International Cosmic Ray Conference held in Pune, India, in August 2005. The paper is based on a plenary talk given at the TAUP2005 meeting in Zaragoza, 10 - 14 September 2005.Comment: 7 pages and two figure

A Three-Point Cosmic Ray Anisotropy Method

The two-point angular correlation function is a traditional method used to search for deviations from expectations of isotropy. In this paper we develop and explore a statistically descriptive three-point method with the intended application being the search for deviations from isotropy in the highest energy cosmic rays. We compare the sensitivity of a two-point method and a "shape-strength" method for a variety of Monte-Carlo simulated anisotropic signals. Studies are done with anisotropic source signals diluted by an isotropic background. Type I and II errors for rejecting the hypothesis of isotropic cosmic ray arrival directions are evaluated for four different event sample sizes: 27, 40, 60 and 80 events, consistent with near term data expectations from the Pierre Auger Observatory. In all cases the ability to reject the isotropic hypothesis improves with event size and with the fraction of anisotropic signal. While ~40 event data sets should be sufficient for reliable identification of anisotropy in cases of rather extreme (highly anisotropic) data, much larger data sets are suggested for reliable identification of more subtle anisotropies. The shape-strength method consistently performs better than the two point method and can be easily adapted to an arbitrary experimental exposure on the celestial sphere.Comment: Fixed PDF erro

A method to suppress dielectric breakdowns in liquid argon ionization detectors for cathode to ground distances of several millimeters

We present a method to reach electric field intensity as high as 400 kV/cm in liquid argon for cathode-ground distances of several millimeters. This can be achieved by suppressing field emission from the cathode, overcoming limitations that we reported earlier

On the Electric Breakdown in Liquid Argon at Centimeter Scale

We present a study on the dependence of electric breakdown discharge properties on electrode geometry and the breakdown field in liquid argon near its boiling point. The measurements were performed with a spherical cathode and a planar anode at distances ranging from 0.1 mm to 10.0 mm. A detailed study of the time evolution of the breakdown volt-ampere characteristics was performed for the first time. It revealed a slow streamer development phase in the discharge. The results of a spectroscopic study of the visible light emission of the breakdowns complement the measurements. The light emission from the initial phase of the discharge is attributed to electro-luminescence of liquid argon following a current of drifting electrons. These results contribute to set benchmarks for breakdown-safe design of ionization detectors, such as Liquid Argon Time Projection Chambers (LAr TPC).Comment: Minor revision according to editor report. 17 pages, 15 figures, 2 tables. Turboencabulato

Prospects for GMRT to Observe Radio Waves from UHE Particles Interacting with the Moon

Ultra high energy (UHE) particles of cosmic origin impact the lunar regolith and produce radio signals through Askaryan effect, signals that can be detected by Earth based radio telescopes. We calculate the expected sensitivity for observation of such events at the Giant Metrewave Radio Telescope (GMRT), both for UHE cosmic rays (CR) and UHE neutrino interactions. We find that for 30 days of observation time a significant number of detectable events is expected above $10^{20}$ eV for UHECR or neutrino fluxes close to the current limits. Null detection over a period of 30 days will lower the experimental bounds on UHE particle fluxes by magnitudes competitive to both present and future experiments at the very highest energies.Comment: 21 pages, 9 figure