First operation of bulk micromegas in low pressure negative ion drift gas mixtures for dark matter searches

A bulk micromegas micropattern charge readout device has, for the first time, been operated at room temperature in low pressure carbon disulphide vapour. This is a key step opening prospects for use of micromegas readout for large volume negative ion time projection chambers (TPCs) without magnets, such as proposed for directional dark matter detectors and other rare event applications. The dependence of the gain on the amplification field, pressure and drift field has been evaluated. For the available gap size of 75 μm a maximum gain of 1300 ± 120 was achieved in 40 torr vapour with an energy resolution of 22% for 5.9 keV ⁵⁵Fe X-rays. From a fit to the data, the Townsend coefficient gas parameters A and B have been derived. Operation has also been successfully achieved in xenon:carbon disulphide blends over a range of partial and total pressures. A gain of 890 ± 130 at an energy resolution of 35% has been recorded for a 1:1 blend at a total pressure of 80 torr. Possible improvements are discussed in the context of operation in directional dark matter TPCs as a replacement for multi-wire proportional counters.     

GEM-based TPC with CCD imaging for directional dark matter detection

The most mature directional dark matter experiments at present all utilize low-pressure gas Time Projection Chamber (TPC) technologies. We discuss some of the challenges for this technology, for which balancing the goal of achieving the best sensitivity with that of cost effective scale-up requires optimization over a large parameter space. Critical for this are the precision measurements of the fundamental properties of both electron and nuclear recoil tracks down to the lowest detectable energies. Such measurements are necessary to provide a benchmark for background discrimination and directional sensitivity that could be used for future optimization studies for directional dark matter experiments. In this paper we describe a small, high resolution, high signal-to-noise GEM-based TPC with a 2D CCD readout designed for this goal. The performance of the detector was characterized using alpha particles, X-rays, gamma-rays, and neutrons, enabling detailed measurements of electron and nuclear recoil tracks. Stable effective gas gains of greater than 1 × 10⁵ were obtained in 100 Torr of pure CF₄ by a cascade of three standard CERN GEMs each with a 140 µm pitch. The high signal-to-noise and sub-millimeter spatial resolution of the GEM amplification and CCD readout, together with low diffusion, allow for excellent background discrimination between electron and nuclear recoils down below ∼10 keVee (∼23 keVr fluorine recoil). Even lower thresholds, necessary for the detection of low mass WIMPs for example, might be achieved by lowering the pressure and utilizing full 3D track reconstruction. These and other paths for improvements are discussed, as are possible fundamental limitations imposed by the physics of energy loss.     

Observation and Study of the CO, CO, HCO and N2H Molecular Lines Towards IRAS 02232+6138

Using the 13.7 m millimeter-wave telescope at the Qinghai Station of Purple Mountain Observatory, we have made observations of ¹³CO, C¹⁸O, HCO⁺ and N₂H⁺ molecular lines towards IRAS 02232+6138. As the excitation density of the probe molecule increases from ¹³CO to HCO⁺, the size of the cloud core associated with IRAS 02232+6138 decreases from 2.40 pc to 0.54 pc, and the virial mass of the cloud core decreases from 2.2 × 10³M to 5.1 × 10²M. A bipolar molecular outflow is found towards IRAS 02232+6138. Using the power function n(r) ∝ r⁻ to fit the spatial density structure of the cloud core, we obtain the power-law index  = 2.3 − 1.2; and we find that, as the probed density increases, the power function becomes more flat. The abundance ratio of ¹³CO to C¹⁸O is 12.4 ± 6.9, comparable with the values 11.8 ± 5.9 for dark clouds and the values 9.0–15.6 for massive cores. The abundance of N₂H⁺ molecules is 3.5 ± 2.5 × 10⁻¹⁰, consistent with the value 1.0 − 5.0 × 10⁻¹⁰ for dark cloud cores and the value 1.2 − 12.8 × 10⁻¹⁰ for massive cores. The abundance of HCO⁺ molecules is 0.9 ± 0.5 × 10⁻⁹, close to the value 1.6 − 2.4 × 10⁻⁹ for massive cores. An increase of HCO⁺ abundance in the outflow region was not found. Combining with the IRAS data, the luminosity-mass ratio of the cloud core is obtained in the range 37–163(L/M). Based on the IRAS luminosity, it is estimated that a main-sequence O7.5 star is probably embedded in the IRAS 02232+6138 cloud core.     

Resolving cosmic structure formation with the Millennium-II Simulation

We present the Millennium-II Simulation (MS-II), a very large N -body simulation of dark matter evolution in the concordance Λ cold dark matter (ΛCDM) cosmology. The MS-II assumes the same cosmological parameters and uses the same particle number and output data structure as the original Millennium Simulation (MS), but was carried out in a periodic cube one-fifth the size  (100  h ⁻¹ Mpc)  with five times better spatial resolution (a Plummer equivalent softening of  1.0  h ⁻¹ kpc  ) and with 125 times better mass resolution (a particle mass of  6.9 × 10⁶  h ⁻¹ M_⊙  ). By comparing results at MS and MS-II resolution, we demonstrate excellent convergence in dark matter statistics such as the halo mass function, the subhalo abundance distribution, the mass dependence of halo formation times, the linear and non-linear autocorrelations and power spectra, and halo assembly bias. Together, the two simulations provide precise results for such statistics over an unprecedented range of scales, from haloes similar to those hosting Local Group dwarf spheroidal galaxies to haloes corresponding to the richest galaxy clusters. The 'Milky Way' haloes of the Aquarius Project were selected from a lower resolution version of the MS-II and were then resimulated at much higher resolution. As a result, they are present in the MS-II along with thousands of other similar mass haloes. A comparison of their assembly histories in the MS-II and in resimulations of 1000 times better resolution shows detailed agreement over a factor of 100 in mass growth. We publicly release halo catalogues and assembly trees for the MS-II in the same format within the same archive as those already released for the MS.     

Measurement of Na and I recoil discrimination and detection efficiency in a cooled ‘UVIS’ NaI crystal for dark matter searches

The UVIS dark matter detector, proposed by Spooner and Smith (Phys. Lett. B 314 (1993) 430), consists of < 1000 ppm Tl doped NaI scintillator operated at 100-200 K, with measurement of UV and visible scintillation components used to discriminate between electrons (gammas) and Na or I recoils. Presented here are results of measurements of the gamma/nuclear recoil discrimination power of a Kyropolous NaI (0.5 ppm Tl) UVIS test detector operated at 160 K, with monoenergetic neutrons used to induce nuclear recoils via elastic scattering. Defining the statistical gain factor for N events as , the coefficient C_m was measured as a function of photoelectron pulse height. At 40 photoelectrons C_m was found to be 0.5 (corresponding to 90% rejection of gammas with 20% loss of neutron events). The scintillation efficiency for Na and I recoils relative to electrons was also measured at 160 K and found to be 35 ± 5% for Na recoils and 10 ± 2% for iodine.     

Dark matter haloes within clusters

We examine the properties of dark matter haloes within a rich galaxy cluster using a high-resolution simulation that captures the cosmological context of a cold dark matter universe. The mass and force resolution permit the resolution of 150 haloes with circular velocities larger than 80 km s⁻¹ within the cluster virial radius of 2 Mpc (with Hubble constant H ₀ = 50 km s⁻¹ Mpc⁻¹). This enables an unprecedented study of the statistical properties of a large sample of dark matter haloes evolving in a dense environment. The cumulative fraction of mass attached to these haloes varies from close to zero per cent at 200 kpc to 13 per cent at the virial radius. Even at this resolution the overmerging problem persists; haloes that pass within 100–200 kpc of the cluster centre are tidally disrupted. Additional substructure is lost at earlier epochs within the massive progenitor haloes. The median ratio of apocentric to pericentric radii is 6:1, so that the orbital distribution is close to isotropic, circular orbits are rare and radial orbits are common. The orbits of haloes are unbiased with respect to both position within the cluster and the orbits of the smooth dark matter background, and no velocity bias is detected. The tidal radii of surviving haloes are generally well-fitted using the simple analytic prediction applied to their orbital pericentres. Haloes within clusters have higher concentrations than those in the field. Within the cluster, halo density profiles can be modified by tidal forces and individual encounters with other haloes that cause significant mass loss —'galaxy harassment'. Mergers between haloes do not occur inside the cluster virial radius.     

Gamma ray spectrum of the crab nebula in the multi TeV region

The THEMISTOCLE array of 18 Cherenkov detectors which has a 3 TeV gamma energy threshold, has detected a signal from the Crab nebula at a 5.8 standard deviation level. Information on the energy spectrum is obtained in the range 3–15 TeV. The integrated flux can be fitted with the form, Φ (> E) = (3.7 ± 0.5) × 10^-12 (E/5)-^{−1.5 ± 0.20} cm⁻² s⁻¹ (E in TeV) compatible with the extrapolation of results at lower energies. The Crab signal is used to measure the angular resolution of the multi-telescope technique. The value obtained is 2.3 mr (0.15°) in agreement with the results of simulations, and confirms the interest of this new method for multi-TeV gamma-ray detection.     

Long-term Variability Properties of Mkn501

Optical photometric observations at the I and R wavebands were carried out towards Mkn501 using the 1.56 meter optical telescope of Shanghai Astronomical Observatory at Sheshan. Combining our new observations with the published historical data, we have obtained the light curves of Mkn501 at the optical, infrared and radio wavebands with a time coverage of nearly 80 years. The relationship between the light variability and the color index is discussed, it is found that a strong correlation exists between the color indices (B − V) and (B − R), with the correlation coefficient reaching r = 0.73. The correlations of multi-band light variabilities are analyzed by the DCF method, it is found that certain positive correlations of B-band light variability with the 4.8 GHz and infrared light variabilities exist. And the spectral analysis on the B-band light curve with the CLEANest method indicates that the light curve of Mkn501 contains probably two periodical components of (10.06 ± 0.04) yr and (21.60 ± 0.17) yr.     

Voids in a ΛCDM universe

We study the formation and evolution of voids in the dark matter distribution using various simulations of the popular Λ cold dark matter cosmogony. We identify voids by requiring them to be regions of space with a mean overdensity of −0.8 or less – roughly the equivalent of using a spherical overdensity group finder for haloes. Each of the simulations contains thousands of voids. The distribution of void sizes in the different simulations shows good agreement when differences in particle and grid resolution are accounted for. Voids very clearly correspond to minima in the smoothed initial density field. Apart from a very weak dependence on the mass resolution, the rescaled mass profiles of voids in the different simulations agree remarkably well. We find a universal void mass profile of the form  ρ(< r )/ρ( r _eff) ∝ exp[( r / r _eff)^α]  , where r _eff is the effective radius of a void and  α∼ 2  . The mass function of haloes in voids is steeper than that of haloes that populate denser regions. In addition, the abundances of void haloes seem to evolve somewhat more strongly between redshifts ∼1 and 0 than the global abundances of haloes.     

Does the Fornax dwarf spheroidal have a central cusp or core?

The dark matter dominated Fornax dwarf spheroidal has five globular clusters orbiting at ∼1 kpc from its centre. In a cuspy cold dark matter halo the globulars would sink to the centre from their current positions within a few Gyr, presenting a puzzle as to why they survive undigested at the present epoch. We show that a solution to this timing problem is to adopt a cored dark matter halo. We use numerical simulations and analytic calculations to show that, under these conditions, the sinking time becomes many Hubble times; the globulars effectively stall at the dark matter core radius. We conclude that the Fornax dwarf spheroidal has a shallow inner density profile with a core radius constrained by the observed positions of its globular clusters. If the phase space density of the core is primordial then it implies a warm dark matter particle and gives an upper limit to its mass of ∼0.5 keV, consistent with that required to significantly alleviate the substructure problem.     

Properties of the infrared dark cloud G79.2+0.38

The MSX infrared dark cloud G79.2+0.38 has been observed over a 11′×′ region simultaneously in the J=1-0 rotational transition lines of the ¹²CO and its isotopic molecules ¹³CO and ¹⁸CO. The dense molecular cores defined by the C¹⁸O line are found to be associated with the two high-extinction patches shown in the MSX A-band image. The two dense cores have the column density N (H₂) (5 – 12) × 10²² cm⁻² and the mean number density n (3 ± 1) × 10⁴ cm⁻³. Their sizes are 1.7 and 1.2 pc in ¹³CO(1-0) line, 1.2 and 0.6 pc in C¹⁸O(1-0) line, respectively. The masses of these cloud cores are estimated to be in the range from 2 × 10² to 2 × 10³ M. The profile of radial mean density of the cloud core can be described by the exponential function ¯n(p) p^{−0.34±0.02}. Compared with the cases of typical optical dark clouds, the abundances of the CO isotopic molecules ¹³CO and C¹⁸O in this MSX infrared dark cloud appear to be depleted by a factor of 4–11, but at present there is no evidence for any obvious variation of the relative abundance ratio X_13/18 between ¹³CO and C¹⁸O with the column density.     

The HNCO ring in the Sgr B2 region

The large-scale structure associated with the 2′N HNCO peak in Sgr B2 [Minh, Y.C., Haikala, L., Hjalmarson, Å., Irvine, W.M., 1998. ApJ 498, 261 (Paper I)] has been investigated. A ring-like morphology of the HNCO emission has been found; this structure may be colliding with the Principal Cloud of Sgr B2. This “HNCO Ring” appears to be centered at (l,b) = (0.7°,−0.07°), with a radius of 5 pc and a total mass of 1.0 × 10⁵ to 1.6 × 10⁶ M. The expansion velocity of the Ring is estimated to be 30–40 km s⁻¹, which gives an expansion time scale of 1.5 × 10⁵ year. The morphology suggests that collision between the Ring and the Principal Cloud may be triggering the massive star formation in the Sgr B2 cloud sequentially, with the latest star formation taking place at the 2′N position. The chemistry related to HNCO is not certain yet, but if it forms mainly via reaction with the evaporated OCN⁻ from icy grain mantles, the observed enhancement of the HNCO abundance can be understood as resulting from shocks associated with the collision between the Principal Cloud and the expanding HNCO Ring.     

The tightening of wide binaries in dSph galaxies through dynamical friction as a test of the dark matter hypothesis

We estimate the time-scales for orbital decay of wide binaries embedded within dark matter haloes, due to dynamical friction against the dark matter particles. We derive analytical scalings for this decay and calibrate and test them through the extensive use of N -body simulations, which accurately confirm the predicted temporal evolution. For density and velocity dispersion parameters as inferred for the dark matter haloes of local dSph galaxies, we show that the decay time-scales become shorter than the ages of the dSph stellar populations for binary stars composed of  1 M_⊙  stars, for initial separations larger than 0.1 pc. Such wide binaries are conspicuous and have been well measured in the solar neighbourhood. The prediction of the dark matter hypothesis is that they should now be absent from stellar populations embedded within low velocity dispersion, high-density dark mater haloes, as currently inferred for the local dSph galaxies, having since evolved into tighter binaries. Relevant empirical determinations of this will become technically feasible in the near future, and could provide evidence to discriminate between dark matter particle haloes or modified gravitational theories, to account for the high dispersion velocities measured for stars in local dSph galaxies.     

Rydberg matter in space: low-density condensed dark matter

Here Rydberg matter is proposed as a candidate for the missing dark matter or dark baryonic matter in the Universe. Spectroscopic and other experimental studies give valuable information on the properties of Rydberg matter, especially its very weak interaction with light caused by the very small overlap with low states, and because of the necessary two-electron transitions even for disturbed matter. Recently, the unidentified infrared (UIR) bands have been shown to agree well with calculations and experiments on Rydberg matter. This is the reason for the present, somewhat speculative, proposal that dark matter has, at least partially, the form of Rydberg matter. The UIR bands have also been observed directly in emission from Rydberg matter in the laboratory. The unique space-filling properties of Rydberg matter are described: a hydrogen atom in this matter occupies a volume  5×10¹²  times larger than in its ground state or in a hydrogen molecule.     

The Bullet Cluster 1E0657-558 evidence shows modified gravity in the absence of dark matter

A detailed analysis of the 2006 November 15 data release X-ray surface density Σ-map and the strong and weak gravitational lensing convergence κ-map for the Bullet Cluster 1E0657-558 is performed and the results are compared with the predictions of a modified gravity (MOG) and dark matter. Our surface density Σ-model is computed using a King β-model density, and a mass profile of the main cluster and an isothermal temperature profile are determined by the MOG. We find that the main cluster thermal profile is nearly isothermal. The MOG prediction of the isothermal temperature of the main cluster is   T = 15.5 ± 3.9 keV  , in good agreement with the experimental value   T = 14.8^+2.0_−1.7 keV  . Excellent fits to the 2D convergence κ-map data are obtained without non-baryonic dark matter, accounting for the 8σ spatial offset between the Σ-map and the κ-map reported in Clowe et al. The MOG prediction for the κ-map results in two baryonic components distributed across the Bullet Cluster 1E0657-558 with averaged mass fraction of 83 per cent intracluster medium (ICM) gas and 17 per cent galaxies. Conversely, the Newtonian dark matter κ-model has on average 76 per cent dark matter (neglecting the indeterminant contribution due to the galaxies) and 24 per cent ICM gas for a baryon to dark matter mass fraction of 0.32, a statistically significant result when compared to the predicted Λ-cold dark matter cosmological baryon mass fraction of 0.176^+0.019_−0.012.     

On the nature of superoutbursts in dwarf novae

We determine a crucial feature of the dark halo density distribution from the fact that the luminous matter dominates the gravitational potential at about one disc scalelength R _d, but at the optical edge     the dark matter has already become the main component of the galaxy density. From the kinematics of 137 spirals we find that the dark matter halo density profiles are self-similar at least out to R _opt and show core radii much larger than the corresponding disc scalelengths. The luminous regions of spirals consist of stellar discs embedded in dark haloes with roughly constant density. This invariant dark matter profile is very difficult to reconcile with the fundamental properties of the density distribution of cold dark matter haloes. With respect to previous work, the present evidence is obtained by means of a robust method and for a large and complete sample of normal spirals.     

Influence of the solar activity on the Indian Monsoon rainfall

We use 130 years data for studying correlative effects due to solar cycle and activity phenomena on the occurrence of the Indian Monsoon rainfall. We compute the correlation coefficients and significance of correlation coefficients for the seasonal and the annual data. We find that: (i) for the whole years 1871–2000, the spring and southwest monsoon rainfall variabilities have significant positive correlations with the sunspot activity during the corresponding period, (ii) the FFT and the wavelet analyses of the southwest monsoon rainfall variability show the periods 2.7, 16 and 22 year, respectively (similar to the periods found in sunspot occurrence data) and, (iii) there is a long-term trend indicating a gradual decrease of occurrence of rainfall variability by nearly 2.3 ± 1.3 mm/year and increase of sunspot activity by nearly 3.9 ± 1.5 sunspots/year compared to the activity of previous solar cycle.

We speculate in this study a possible physical connection between the occurrence of the rainfall variability and the sunspot activity, and the flux of galactic cosmic rays. Owing to long-term positive and significant correlation of the spring and southwest monsoon rainfall variabilities with the sunspot activity, it is suggested that solar activity may be included as one of the crucial parameter in modeling and predicting the Indian monsoon rainfall.     

Is an 11 eV sterile neutrino consistent with clusters, the cosmic microwave background and modified Newtonian dynamics?

In this paper, we show that if a single sterile neutrino exists such that     , it can serendipitously solve all outstanding issues of the Modified Newtonian Dynamics. We focus on fitting the angular power spectrum of the cosmic microwave background (CMB) in detail which is possible using a flat Universe with     and the usual baryonic and dark energy components. One cannot match the CMB if there is more than one massive sterile neutrino, nor with three active neutrinos of 2 eV. This model has the same expansion history as the Λ cold dark matter  (ΛCDM)  model and only differs at the galactic scale, where the modified dynamics outperform  ΛCDM  comprehensively. We discuss how an 11 eV sterile neutrino can explain the dark matter of galaxy clusters without influencing individual galaxies and potentially match the matter power spectrum.     

Geminid Meteor shower activity 2003–2005 as observed by Gadanki radar

The distribution of meteor signals reflected from a backscatter radar is considered according to their duration. This duration time (T) is used to classify the meteor echoes and to calculate the mass index (S) of different meteoroids of shower plus sporadic background. Observational data on particle size distribution of the Geminid meteor shower are very scarce, particularly at low latitudes. In this paper the observational data from Gadanki radar (13.46°N, 79.18°E) have been used to determine the particle size distribution and the number density of meteoroids inside the stream of the Geminid meteor shower. The mean variation of meteor number density across the stream has been determined for three echo duration classes, T<0.4, T=0.4–1 and T>1 s. We are more interested in the appearance of echoes of various durations and therefore meteors of various masses in order to understand more on the filamentary structure of the stream. It is observed that the faint particle flux peaks earlier than the larger particles. We found a decreasing trend in the mass index values from the day of peak activity to the next observation days. The mass index profile was found to be U-shaped with a minimum value near the time of peak activity. The observed minimum s values are 1.64±0.05 and 1.65±0.04 in the years 2003 and 2005, respectively. The activity of the shower indicates the mass segregation of meteoroids inside the stream. Our results are best comparable with the “scissors” structure model of the meteoroid stream formation of Ryabova [2007. Mathematical modeling of the Geminid meteoroid stream. Mon. Not. R. Astron. Soc. 375, 1371–1380] by considering the asteroid 3200 Phaethon as an extinct comet.     

L 43: the late stages of a molecular outflow

Our new 21-arcsec resolution CO J  = 2 → 1 map of the L 43 dark cloud shows a poorly collimated molecular outflow, with little evidence for wings at velocities 10 km s⁻¹. The outflow appears not to be currently driven by a jet: its structure can instead be modelled as a slowly expanding shell. The shell may be compressed either by a wide-angled wind catching up with an existing shell (as in the case of planetary nebulæ), or by the thermal pressure of a hot low-emissivity medium interior to the shell. The outflow is most probably in a late stage of evolution, and appears to be in the process of blowing away its molecular cloud. We also present a 45-arcsec resolution CO J  = 1 → 0 map of the whole molecular cloud, showing that the outflow structure is clearly visible even in the integrated intensity of this low excitation line, and suggesting that rapid mapping may prove useful as a way of finding regions of outflow activity. We also examine the immediate surroundings of the driving source with 450 μm imaging: this confirms that the outflow has already evacuated a bay in the vicinity of the young stellar object.