Limits on MeV Dark Matter from the Effective Number of Neutrinos

Thermal dark matter that couples more strongly to electrons and photons than to neutrinos will heat the electron-photon plasma relative to the neutrino background if it becomes nonrelativistic after the neutrinos decouple from the thermal background. This results in a reduction in N_eff below the standard-model value, a result strongly disfavored by current CMB observations. Taking conservative lower bounds on N_eff and on the decoupling temperature of the neutrinos, we derive a bound on the dark matter particle mass of m_\chi > 3-9 MeV, depending on the spin and statistics of the particle. For p-wave annihilation, our limit on the dark matter particle mass is stronger than the limit derived from distortions to the CMB fluctuation spectrum produced by annihilations near the epoch of recombination.Comment: 5 pages, 1 figure, discussion added, references added and updated, labels added to figure, to appear in Phys. Rev.

Post-hoc derivation of SOHO Michelson doppler imager flat fields

<p><b>Context:</b> The SOHO satellite now offers a unique perspective on the Sun as it is the only space-based instrument that can provide large, high-resolution data sets over an entire 11-year solar cycle. This unique property enables detailed studies of long-term variations in the Sun. One significant problem when looking for such changes is determining what component of any variation is due to deterioration of the instrument and what is due to the Sun itself. One of the key parameters that changes over time is the apparent sensitivity of individual pixels in the CCD array. This can change considerably as a result of optics damage, radiation damage, and aging of the sensor itself. In addition to reducing the sensitivity of the telescope over time, this damage significantly changes the uniformity of the flat field of the instrument, a property that is very hard to recalibrate in space. For procedures such as feature tracking and intensity analysis, this can cause significant errors.</p> <p><b>Aims:</b> We present a method for deriving high-precision flat fields for high-resolution MDI continuum data, using analysis of existing continuum and magnetogram data sets.</p> <p><b>Methods:</b> A flat field is constructed using a large set (1000-4000 frames) of cospatial magnetogram and continuum data. The magnetogram data is used to identify and mask out magnetically active regions on the continuum data, allowing systematic biases to be avoided. This flat field can then be used to correct individual continuum images from a similar time.</p> <p><b>Results:</b> This method allows us to reduce the residual flat field error by around a factor 6-30, depending on the area considered, enough to significantly change the results from correlation-tracking analysis. One significant advantage of this method is that it can be done retrospectively using archived data, without requiring any special satellite operations.</p&gt

A method for the estimation of p-mode parameters from averaged solar oscillation power spectra

A new fitting methodology is presented which is equally well suited for the estimation of low-, medium-, and high-degree mode parameters from $m$-averaged solar oscillation power spectra of widely differing spectral resolution. This method, which we call the "Windowed, MuLTiple-Peak, averaged spectrum", or WMLTP Method, constructs a theoretical profile by convolving the weighted sum of the profiles of the modes appearing in the fitting box with the power spectrum of the window function of the observing run using weights from a leakage matrix that takes into account both observational and physical effects, such as the distortion of modes by solar latitudinal differential rotation. We demonstrate that the WMLTP Method makes substantial improvements in the inferences of the properties of the solar oscillations in comparison with a previous method that employed a single profile to represent each spectral peak. We also present an inversion for the internal solar structure which is based upon 6,366 modes that we have computed using the WMLTP method on the 66-day long 2010 SOHO/MDI Dynamics Run. To improve both the numerical stability and reliability of the inversion we developed a new procedure for the identification and correction of outliers in a frequency data set. We present evidence for a pronounced departure of the sound speed in the outer half of the solar convection zone and in the subsurface shear layer from the radial sound speed profile contained in Model~S of Christensen-Dalsgaard and his collaborators that existed in the rising phase of Solar Cycle~24 during mid-2010

On the mass of atoms in molecules: Beyond the Born-Oppenheimer approximation

Describing the dynamics of nuclei in molecules requires a potential energy surface, which is traditionally provided by the Born-Oppenheimer or adiabatic approximation. However, we also need to assign masses to the nuclei. There, the Born-Oppenheimer picture does not account for the inertia of the electrons and only bare nuclear masses are considered. Nowadays, experimental accuracy challenges the theoretical predictions of rotational and vibrational spectra and requires to include the participation of electrons in the internal motion of the molecule. More than 80 years after the original work of Born and Oppenheimer, this issue still is not solved in general. Here, we present a theoretical and numerical framework to address this problem in a general and rigorous way. Starting from the exact factorization of the electron-nuclear wave function, we include electronic effects beyond the Born-Oppenheimer regime in a perturbative way via position-dependent corrections to the bare nuclear masses. This maintains an adiabatic-like point of view: the nuclear degrees of freedom feel the presence of the electrons via a single potential energy surface, whereas the inertia of electrons is accounted for and the total mass of the system is recovered. This constitutes a general framework for describing the mass acquired by slow degrees of freedom due to the inertia of light, bounded particles. We illustrate it with a model of proton transfer, where the light particle is the proton, and with corrections to the vibrational spectra of molecules. Inclusion of the light particle inertia allows to gain orders of magnitude in accuracy

Balltracking: an highly efficient method for tracking flow fields

We present a method for tracking solar photospheric flows that is highly efficient, and demonstrate it using high resolution MDI continuum images. The method involves making a surface from the photospheric granulation data, and allowing many small floating tracers or balls to be moved around by the evolving granulation pattern. The results are tested against synthesised granulation with known flow fields and compared to the results produced by Local Correlation tracking (LCT). The results from this new method have similar accuracy to those produced by LCT. We also investigate the maximum spatial and temporal resolution of the velocity field that it is possible to extract, based on the statistical properties of the granulation data. We conclude that both methods produce results that are close to the maximum resolution possible from granulation data. The code runs very significantly faster than our similarly optimised LCT code, making real time applications on large data sets possible. The tracking method is not limited to photospheric flows, and will also work on any velocity field where there are visible moving features of known scale length

Cosmic String Formation from Correlated Fields

We simulate the formation of cosmic strings at the zeros of a complex Gaussian field with a power spectrum $P(k) \propto k^n$, specifically addressing the issue of the fraction of length in infinite strings. We make two improvements over previous simulations: we include a non-zero random background field in our box to simulate the effect of long-wavelength modes, and we examine the effects of smoothing the field on small scales. The inclusion of the background field significantly reduces the fraction of length in infinite strings for $n < -2$. Our results are consistent with the possibility that infinite strings disappear at some $n = n_c$ in the range $-3 \le n_c < -2.2$, although we cannot rule out $n_c = -3$, in which case infinite strings would disappear only at the point where the mean string density goes to zero. We present an analytic argument which suggests the latter case. Smoothing on small scales eliminates closed loops on the order of the lattice cell size and leads to a ``lattice-free" estimate of the infinite string fraction. As expected, this fraction depends on the type of window function used for smoothing.Comment: 24 pages, latex, 10 figures, submitted to Phys Rev