Cosmological Evolution of Supergiant Star-Forming Clouds

In an exploration of the birthplaces of globular clusters, we present a careful examination of the formation of self-gravitating gas clouds within assembling dark matter haloes in a hierarchical cosmological model. Our high-resolution smoothed particle hydrodynamical simulations are designed to determine whether or not hypothesized supergiant molecular clouds (SGMCs) form and, if they do, to determine their physical properties and mass spectra. It was suggested in earlier work that clouds with a median mass of several 10^8 M_sun are expected to assemble during the formation of a galaxy, and that globular clusters form within these SGMCs. Our simulations show that clouds with the predicted properties are indeed produced as smaller clouds collide and agglomerate within the merging dark matter haloes of our cosmological model. We find that the mass spectrum of these clouds obeys the same power-law form observed for globular clusters, molecular clouds, and their internal clumps in galaxies, and predicted for the supergiant clouds in which globular clusters may form. We follow the evolution and physical properties of gas clouds within small dark matter haloes up to z = 1, after which prolific star formation is expected to occur. Finally, we discuss how our results may lead to more physically motivated "rules" for star formation in cosmological simulations of galaxy formation.Comment: Accepted to The Astrophysical Journal; 17 pages, 8 figure

Full-sky maps for gravitational lensing of the CMB

We use the large cosmological Millennium Simulation (MS) to construct the first all-sky maps of the lensing potential and the deflection angle, aiming at gravitational lensing of the CMB, with the goal of properly including small-scale non-linearities and non-Gaussianity. Exploiting the Born approximation, we implement a map-making procedure based on direct ray-tracing through the gravitational potential of the MS. We stack the simulation box in redshift shells up to $z\sim 11$, producing continuous all-sky maps with arcminute angular resolution. A randomization scheme avoids repetition of structures along the line of sight and structures larger than the MS box size are added to supply the missing contribution of large-scale (LS) structures to the lensing signal. The angular power spectra of the projected lensing potential and the deflection-angle modulus agree quite well with semi-analytic estimates on scales down to a few arcminutes, while we find a slight excess of power on small scales, which we interpret as being due to non-linear clustering in the MS. Our map-making procedure, combined with the LS adding technique, is ideally suited for studying lensing of CMB anisotropies, for analyzing cross-correlations with foreground structures, or other secondary CMB anisotropies such as the Rees-Sciama effect.Comment: LaTeX file, 10 pages, MNRAS in press, scales larger than the Millennium Simulation box size semi-analytically added, maps changed, references added, typos correcte

Dark Matter Direct Detection with Non-Maxwellian Velocity Structure

The velocity distribution function of dark matter particles is expected to show significant departures from a Maxwell-Boltzmann distribution. This can have profound effects on the predicted dark matter - nucleon scattering rates in direct detection experiments, especially for dark matter models in which the scattering is sensitive to the high velocity tail of the distribution, such as inelastic dark matter (iDM) or light (few GeV) dark matter (LDM), and for experiments that require high energy recoil events, such as many directionally sensitive experiments. Here we determine the velocity distribution functions from two of the highest resolution numerical simulations of Galactic dark matter structure (Via Lactea II and GHALO), and study the effects for these scenarios. For directional detection, we find that the observed departures from Maxwell-Boltzmann increase the contrast of the signal and change the typical direction of incoming DM particles. For iDM, the expected signals at direct detection experiments are changed dramatically: the annual modulation can be enhanced by more than a factor two, and the relative rates of DAMA compared to CDMS can change by an order of magnitude, while those compared to CRESST can change by a factor of two. The spectrum of the signal can also change dramatically, with many features arising due to substructure. For LDM the spectral effects are smaller, but changes do arise that improve the compatibility with existing experiments. We find that the phase of the modulation can depend upon energy, which would help discriminate against background should it be found.Comment: 34 pages, 16 figures, submitted to JCAP. Tables of g(v_min), the integral of f(v)/v from v_min to infinity, derived from our simulations, are available for download at http://astro.berkeley.edu/~mqk/dmdd

A Nonlinear Force-Free Magnetic Field Approximation Suitable for Fast Forward-Fitting to Coronal Loops. II. Numeric Code and Tests

Based on a second-order approximation of nonlinear force-free magnetic field solutions in terms of uniformly twisted field lines derived in Paper I, we develop here a numeric code that is capable to forward-fit such analytical solutions to arbitrary magnetogram (or vector magnetograph) data combined with (stereoscopically triangulated) coronal loop 3D coordinates. We test the code here by forward-fitting to six potential field and six nonpotential field cases simulated with our analytical model, as well as by forward-fitting to an exactly force-free solution of the Low and Lou (1990) model. The forward-fitting tests demonstrate: (i) a satisfactory convergence behavior (with typical misalignment angles of $\mu \approx 1^\circ-10^\circ$), (ii) relatively fast computation times (from seconds to a few minutes), and (iii) the high fidelity of retrieved force-free $\alpha$-parameters ($\alpha_{\rm fit}/\alpha_{\rm model} \approx 0.9-1.0$ for simulations and $\alpha_{\rm fit}/\alpha_{\rm model} \approx 0.7\pm0.3$ for the Low and Lou model). The salient feature of this numeric code is the relatively fast computation of a quasi-forcefree magnetic field, which closely matches the geometry of coronal loops in active regions, and complements the existing {\sl nonlinear force-free field (NLFFF)} codes based on photospheric magnetograms without coronal constraints.Comment: Solar PHysics, (in press), 25 pages, 11 figure