Evolutionary algorithm-based analysis of gravitational microlensing lightcurves

A new algorithm developed to perform autonomous fitting of gravitational microlensing lightcurves is presented. The new algorithm is conceptually simple, versatile and robust, and parallelises trivially; it combines features of extant evolutionary algorithms with some novel ones, and fares well on the problem of fitting binary-lens microlensing lightcurves, as well as on a number of other difficult optimisation problems. Success rates in excess of 90% are achieved when fitting synthetic though noisy binary-lens lightcurves, allowing no more than 20 minutes per fit on a desktop computer; this success rate is shown to compare very favourably with that of both a conventional (iterated simplex) algorithm, and a more state-of-the-art, artificial neural network-based approach. As such, this work provides proof of concept for the use of an evolutionary algorithm as the basis for real-time, autonomous modelling of microlensing events. Further work is required to investigate how the algorithm will fare when faced with more complex and realistic microlensing modelling problems; it is, however, argued here that the use of parallel computing platforms, such as inexpensive graphics processing units, should allow fitting times to be constrained to under an hour, even when dealing with complicated microlensing models. In any event, it is hoped that this work might stimulate some interest in evolutionary algorithms, and that the algorithm described here might prove useful for solving microlensing and/or more general model-fitting problems.Comment: 14 pages, 3 figures; accepted for publication in MNRA

FitSKIRT: genetic algorithms to automatically fit dusty galaxies with a Monte Carlo radiative transfer code

We present FitSKIRT, a method to efficiently fit radiative transfer models to UV/optical images of dusty galaxies. These images have the advantage that they have better spatial resolution compared to FIR/submm data. FitSKIRT uses the GAlib genetic algorithm library to optimize the output of the SKIRT Monte Carlo radiative transfer code. Genetic algorithms prove to be a valuable tool in handling the multi- dimensional search space as well as the noise induced by the random nature of the Monte Carlo radiative transfer code. FitSKIRT is tested on artificial images of a simulated edge-on spiral galaxy, where we gradually increase the number of fitted parameters. We find that we can recover all model parameters, even if all 11 model parameters are left unconstrained. Finally, we apply the FitSKIRT code to a V-band image of the edge-on spiral galaxy NGC4013. This galaxy has been modeled previously by other authors using different combinations of radiative transfer codes and optimization methods. Given the different models and techniques and the complexity and degeneracies in the parameter space, we find reasonable agreement between the different models. We conclude that the FitSKIRT method allows comparison between different models and geometries in a quantitative manner and minimizes the need of human intervention and biasing. The high level of automation makes it an ideal tool to use on larger sets of observed data.Comment: 14 pages, 10 figures; accepted for publication in Astronomy and Astrophysic

Advanced Forward Modeling and Inversion of Stokes Profiles Resulting from the Joint Action of the Hanle and Zeeman Effects

A big challenge in solar and stellar physics in the coming years will be to decipher the magnetism of the solar outer atmosphere (chromosphere and corona) along with its dynamic coupling with the magnetic fields of the underlying photosphere. To this end, it is important to develop rigorous diagnostic tools for the physical interpretation of spectropolarimetric observations in suitably chosen spectral lines. Here we present a computer program for the synthesis and inversion of Stokes profiles caused by the joint action of atomic level polarization and the Hanle and Zeeman effects in some spectral lines of diagnostic interest, such as those of the He I 10830 A and D_3 multiplets. It is based on the quantum theory of spectral line polarization, which takes into account all the relevant physical mechanisms and ingredients (optical pumping, atomic level polarization, Zeeman, Paschen-Back and Hanle effects). The influence of radiative transfer on the emergent spectral line radiation is taken into account through a suitable slab model. The user can either calculate the emergent intensity and polarization for any given magnetic field vector or infer the dynamical and magnetic properties from the observed Stokes profiles via an efficient inversion algorithm based on global optimization methods. The reliability of the forward modeling and inversion code presented here is demonstrated through several applications, which range from the inference of the magnetic field vector in solar active regions to determining whether or not it is canopy-like in quiet chromospheric regions. This user-friendly diagnostic tool called "HAZEL" (from HAnle and ZEeman Light) is offered to the astrophysical community, with the hope that it will facilitate new advances in solar and stellar physics.Comment: 62 pages, 19 figures, 3 tables. Accepted for publication in Ap

Stochastic optimization methods for extracting cosmological parameters from CMBR power spectra

The reconstruction of the CMBR power spectrum from a map represents a major computational challenge to which much effort has been applied. However, once the power spectrum has been recovered there still remains the problem of extracting cosmological parameters from it. Doing this involves optimizing a complicated function in a many dimensional parameter space. Therefore efficient algorithms are necessary in order to make this feasible. We have tested several different types of algorithms and found that the technique known as simulated annealing is very effective for this purpose. It is shown that simulated annealing is able to extract the correct cosmological parameters from a set of simulated power spectra, but even with such fast optimization algorithms, a substantial computational effort is needed.Comment: 7 pages revtex, 3 figures, to appear in PR

A low pre-infall mass for the Carina dwarf galaxy from disequilibrium modelling

Dark matter only simulations of galaxy formation predict many more subhalos around a Milky Way like galaxy than the number of observed satellites. Proposed solutions require the satellites to inhabit dark matter halos with masses between one to ten billion solar masses at the time they fell into the Milky Way. Here we use a modelling approach, independent of cosmological simulations, to obtain a preinfall mass of 360 (+380,-230) million solar masses for one of the Milky Way's satellites: Carina. This determination of a low halo mass for Carina can be accommodated within the standard model only if galaxy formation becomes stochastic in halos below ten billion solar masses. Otherwise Carina, the eighth most luminous Milky Way dwarf, would be expected to inhabit a significantly more massive halo. The implication of this is that a population of "dark dwarfs" should orbit the Milky Way: halos devoid of stars and yet more massive than many of their visible counterparts.Comment: 5 pages, 3 figures, 1 table, and supplementary material availabl

Quantum Limitations on the Storage and Transmission of Information

Information must take up space, must weigh, and its flux must be limited. Quantum limits on communication and information storage leading to these conclusions are here described. Quantum channel capacity theory is reviewed for both steady state and burst communication. An analytic approximation is given for the maximum signal information possible with occupation number signal states as a function of mean signal energy. A theorem guaranteeing that these states are optimal for communication is proved. A heuristic "proof" of the linear bound on communication is given, followed by rigorous proofs for signals with specified mean energy, and for signals with given energy budget. And systems of many parallel quantum channels are shown to obey the linear bound for a natural channel architecture. The time--energy uncertainty principle is reformulated in information language by means of the linear bound. The quantum bound on information storage capacity of quantum mechanical and quantum field devices is reviewed. A simplified version of the analytic proof for the bound is given for the latter case. Solitons as information caches are discussed, as is information storage in one dimensional systems. The influence of signal self--gravitation on communication is considerd. Finally, it is shown that acceleration of a receiver acts to block information transfer.Comment: Published relatively inaccessible review on a perennially interesting subject. Plain TeX, 47 pages, 5 jpg figures (not embedded

Discovery of a high-redshift Einstein ring

We report the discovery of a partial Einstein ring of radius 1.48arcsec produced by a massive (and seemingly isolated) elliptical galaxy. The spectroscopic follow-up at the VLT reveals a 2L* galaxy at z=0.986, which is lensing a post-starburst galaxy at z=3.773. This unique configuration yields a very precise measure of the mass of the lens within the Einstein radius, (8.3e11 +- 0.4)/h70 Msolar. The fundamental plane relation indicates an evolution rate of d [log (M/L)B] / dz = -0.57+-0.04, similar to other massive ellipticals at this redshift. The source galaxy shows strong interstellar absorption lines indicative of large gas-phase metallicities, with fading stellar populations after a burst. Higher resolution spectra and imaging will allow the detailed study of an unbiased representative of the galaxy population when the universe was just 12% of its current age.Comment: 5 pages, 3 figures, accepted in A&A Le

The distribution of interstellar dust in CALIFA edge-on galaxies via oligochromatic radiative transfer fitting

We investigate the amount and spatial distribution of interstellar dust in edge-on spiral galaxies, using detailed radiative transfer modeling of a homogeneous sample of 12 galaxies selected from the CALIFA survey. Our automated fitting routine, FitSKIRT, was first validated against artificial data. This is done by simultaneously reproducing the SDSS $g$-, $r$-, $i$- and $z$-band observations of a toy model in order to combine the information present in the different bands. We show that this combined, oligochromatic fitting, has clear advantages over standard monochromatic fitting especially regarding constraints on the dust properties. We model all galaxies in our sample using a three-component model, consisting of a double exponential disc to describe the stellar and dust discs and using a S\'ersic profile to describe the central bulge. The full model contains 19 free parameters, and we are able to constrain all these parameters to a satisfactory level of accuracy without human intervention or strong boundary conditions. Apart from two galaxies, the entire sample can be accurately reproduced by our model. We find that the dust disc is about 75% more extended but only half as high as the stellar disc. The average face-on optical depth in the V-band is $0.76$ and the spread of $0.60$ within our sample is quite substantial, which indicates that some spiral galaxies are relatively opaque even when seen face-on.Comment: 18 pages, 6 figures, 4 tables, Accepted for publication in MNRA