Lyman alpha line formation in starbursting galaxies II. Extremely Thick, Dustless, and Static HI Media

The Lya line transfer in an extremely thick medium of neutral hydrogen is investigated by adopting an accelerating scheme in our Monte Carlo code to skip a large number of core or resonant scatterings. This scheme reduces computing time significantly with no sacrifice in the accuracy of the results. We applied this numerical method to the Lya transfer in a static, uniform, dustless, and plane-parallel medium. Two types of photon sources have been considered, the midplane source and the uniformly distributed sources. The emergent profiles show double peaks and absorption trough at the line-center. We compared our results with the analytic solutions derived by previous researchers, and confirmed that both solutions are in good agreement with each other. We investigated the directionality of the emergent Lya photons and found that limb brightening is observed in slightly thick media while limb darkening appears in extremely thick media. The behavior of the directionality is noted to follow that of the Thomson scattered radiation in electron clouds, because both Lya wing scattering and Thomson scattering share the same Rayleigh scattering phase function. The mean number of wing scatterings just before escape is in exact agreement with the prediction of the diffusion approximation. The Lya photons constituting the inner part of the emergent profiles follow the relationship derived from the diffusion approximation. We present a brief discussion on the application of our results to the formation of Lya broad absorption troughs and P-Cygni type Lya profiles seen in the UV spectra of starburst galaxies.Comment: 24 papges, 12 figures, The revised version submitted to Ap

Solving the Boltzmann Equation on GPU

We show how to accelerate the direct solution of the Boltzmann equation using Graphics Processing Units (GPUs). In order to fully exploit the computational power of the GPU, we choose a method of solution which combines a finite difference discretization of the free-streaming term with a Monte Carlo evaluation of the collision integral. The efficiency of the code is demonstrated by solving the two-dimensional driven cavity flow. Computational results show that it is possible to cut down the computing time of the sequential code of two order of magnitudes. This makes the proposed method of solution a viable alternative to particle simulations for studying unsteady low Mach number flows.Comment: 18 pages, 3 pseudo-codes, 6 figures, 1 tabl

Numerical study of barriers and valleys in the free-energy landscape of spin glasses

We study the problem of glassy relaxations in the presence of an external field in the highly controlled context of a spin-glass simulation. We consider a small spin glass in three dimensions (specifically, a lattice of size L=8, small enough to be equilibrated through a Parallel Tempering simulations at low temperatures, deep in the spin glass phase). After equilibrating the sample, an external field is switched on, and the subsequent dynamics is studied. The field turns out to reduce the relaxation time, but huge statistical fluctuations are found when different samples are compared. After taking care of these fluctuations we find that the expected linear regime is very narrow. Nevertheless, when regarded as a purely numerical method, we find that the external field is extremely effective in reducing the relaxation times.Comment: 22 pages, 10 figures; Published versio

Status and Future Perspectives for Lattice Gauge Theory Calculations to the Exascale and Beyond

In this and a set of companion whitepapers, the USQCD Collaboration lays out a program of science and computing for lattice gauge theory. These whitepapers describe how calculation using lattice QCD (and other gauge theories) can aid the interpretation of ongoing and upcoming experiments in particle and nuclear physics, as well as inspire new ones.Comment: 44 pages. 1 of USQCD whitepapers