Challenges of Profile Likelihood Evaluation in Multi-Dimensional SUSY Scans

Statistical inference of the fundamental parameters of supersymmetric theories is a challenging and active endeavor. Several sophisticated algorithms have been employed to this end. While Markov-Chain Monte Carlo (MCMC) and nested sampling techniques are geared towards Bayesian inference, they have also been used to estimate frequentist confidence intervals based on the profile likelihood ratio. We investigate the performance and appropriate configuration of MultiNest, a nested sampling based algorithm, when used for profile likelihood-based analyses both on toy models and on the parameter space of the Constrained MSSM. We find that while the standard configuration is appropriate for an accurate reconstruction of the Bayesian posterior, the profile likelihood is poorly approximated. We identify a more appropriate MultiNest configuration for profile likelihood analyses, which gives an excellent exploration of the profile likelihood (albeit at a larger computational cost), including the identification of the global maximum likelihood value. We conclude that with the appropriate configuration MultiNest is a suitable tool for profile likelihood studies, indicating previous claims to the contrary are not well founded.Comment: 21 pages, 9 figures, 1 table; minor changes following referee report. Matches version accepted by JHE

Dark Matter in 3D

We discuss the relevance of directional detection experiments in the post-discovery era and propose a method to extract the local dark matter phase space distribution from directional data. The first feature of this method is a parameterization of the dark matter distribution function in terms of integrals of motion, which can be analytically extended to infer properties of the global distribution if certain equilibrium conditions hold. The second feature of our method is a decomposition of the distribution function in moments of a model independent basis, with minimal reliance on the ansatz for its functional form. We illustrate our method using the Via Lactea II N-body simulation as well as an analytical model for the dark matter halo. We conclude that O(1000) events are necessary to measure deviations from the Standard Halo Model and constrain or measure the presence of anisotropies.Comment: 36 pages, 13 figure

Deep Sequencing of Pyrethroid-Resistant Bed Bugs Reveals Multiple Mechanisms of Resistance within a Single Population

A frightening resurgence of bed bug infestations has occurred over the last 10 years in the U.S. and current chemical methods have been inadequate for controlling this pest due to widespread insecticide resistance. Little is known about the mechanisms of resistance present in U.S. bed bug populations, making it extremely difficult to develop intelligent strategies for their control. We have identified bed bugs collected in Richmond, VA which exhibit both kdr-type (L925I) and metabolic resistance to pyrethroid insecticides. Using LD50 bioassays, we determined that resistance ratios for Richmond strain bed bugs were ∼5200-fold to the insecticide deltamethrin. To identify metabolic genes potentially involved in the detoxification of pyrethroids, we performed deep-sequencing of the adult bed bug transcriptome, obtaining more than 2.5 million reads on the 454 titanium platform. Following assembly, analysis of newly identified gene transcripts in both Harlan (susceptible) and Richmond (resistant) bed bugs revealed several candidate cytochrome P450 and carboxylesterase genes which were significantly over-expressed in the resistant strain, consistent with the idea of increased metabolic resistance. These data will accelerate efforts to understand the biochemical basis for insecticide resistance in bed bugs, and provide molecular markers to assist in the surveillance of metabolic resistance

The Herpesvirus Associated Ubiquitin Specific Protease, USP7, Is a Negative Regulator of PML Proteins and PML Nuclear Bodies

The PML tumor suppressor is the founding component of the multiprotein nuclear structures known as PML nuclear bodies (PML-NBs), which control several cellular functions including apoptosis and antiviral effects. The ubiquitin specific protease USP7 (also called HAUSP) is known to associate with PML-NBs and to be a tight binding partner of two herpesvirus proteins that disrupt PML NBs. Here we investigated whether USP7 itself regulates PML-NBs. Silencing of USP7 was found to increase the number of PML-NBs, to increase the levels of PML protein and to inhibit PML polyubiquitylation in nasopharyngeal carcinoma cells. This effect of USP7 was independent of p53 as PML loss was observed in p53-null cells. PML-NBs disruption was induced by USP7 overexpression independently of its catalytic activity and was induced by either of the protein interaction domains of USP7, each of which localized to PML-NBs. USP7 also disrupted NBs formed from some single PML isoforms, most notably isoforms I and IV. CK2α and RNF4, which are known regulators of PML, were dispensable for USP7-associated PML-NB disruption. The results are consistent with a novel model of PML regulation where a deubiquitylase disrupts PML-NBs through recruitment of another cellular protein(s) to PML NBs, independently of its catalytic activity

Atmospheric retrieval of exoplanets

Exoplanetary atmospheric retrieval refers to the inference of atmospheric properties of an exoplanet given an observed spectrum. The atmospheric properties include the chemical compositions, temperature profiles, clouds/hazes, and energy circulation. These properties, in turn, can provide key insights into the atmospheric physicochemical processes of exoplanets as well as their formation mechanisms. Major advancements in atmospheric retrieval have been made in the last decade, thanks to a combination of state-of-the-art spectroscopic observations and advanced atmospheric modeling and statistical inference methods. These developments have already resulted in key constraints on the atmospheric H2O abundances, temperature profiles, and other properties for several exoplanets. Upcoming facilities such as the JWST will further advance this area. The present chapter is a pedagogical review of this exciting frontier of exoplanetary science. The principles of atmospheric retrievals of exoplanets are discussed in detail, including parametric models and statistical inference methods, along with a review of key results in the field. Some of the main challenges in retrievals with current observations are discussed along with new directions and the future landscape

Detection of titanium oxide in the atmosphere of a hot Jupiter.

As an exoplanet transits its host star, some of the light from the star is absorbed by the atoms and molecules in the planet's atmosphere, causing the planet to seem bigger; plotting the planet's observed size as a function of the wavelength of the light produces a transmission spectrum. Measuring the tiny variations in the transmission spectrum, together with atmospheric modelling, then gives clues to the properties of the exoplanet's atmosphere. Chemical species composed of light elements-such as hydrogen, oxygen, carbon, sodium and potassium-have in this way been detected in the atmospheres of several hot giant exoplanets, but molecules composed of heavier elements have thus far proved elusive. Nonetheless, it has been predicted that metal oxides such as titanium oxide (TiO) and vanadium oxide occur in the observable regions of the very hottest exoplanetary atmospheres, causing thermal inversions on the dayside. Here we report the detection of TiO in the atmosphere of the hot-Jupiter planet WASP-19b. Our combined spectrum, with its wide spectral coverage, reveals the presence of TiO (to a confidence level of 7.7σ), a strongly scattering haze (7.4σ) and sodium (3.4σ), and confirms the presence of water (7.9σ) in the atmosphere

All-sky search for gravitational-wave bursts in the second joint LIGO-Virgo run

We present results from a search for gravitational-wave bursts in the data collected by the LIGO and Virgo detectors between July 7, 2009 and October 20, 2010: data are analyzed when at least two of the three LIGO-Virgo detectors are in coincident operation, with a total observation time of 207 days. The analysis searches for transients of duration < 1 s over the frequency band 64-5000 Hz, without other assumptions on the signal waveform, polarization, direction or occurrence time. All identified events are consistent with the expected accidental background. We set frequentist upper limits on the rate of gravitational-wave bursts by combining this search with the previous LIGO-Virgo search on the data collected between November 2005 and October 2007. The upper limit on the rate of strong gravitational-wave bursts at the Earth is 1.3 events per year at 90% confidence. We also present upper limits on source rate density per year and Mpc^3 for sample populations of standard-candle sources. As in the previous joint run, typical sensitivities of the search in terms of the root-sum-squared strain amplitude for these waveforms lie in the range 5 10^-22 Hz^-1/2 to 1 10^-20 Hz^-1/2. The combination of the two joint runs entails the most sensitive all-sky search for generic gravitational-wave bursts and synthesizes the results achieved by the initial generation of interferometric detectors.Comment: 15 pages, 7 figures: data for plots and archived public version at https://dcc.ligo.org/cgi-bin/DocDB/ShowDocument?docid=70814&version=19, see also the public announcement at http://www.ligo.org/science/Publication-S6BurstAllSky

Neutrino tomography of Earth

Cosmic-ray interactions with the atmosphere produce a flux of neutrinos in all directions with energies extending above the TeV scale. The Earth is not a fully transparent medium for neutrinos with energies above a few TeV, as the neutrino–nucleon cross-section is large enough to make the absorption probability non-negligible. Since absorption depends on energy and distance travelled, studying the distribution of the TeV atmospheric neutrinos passing through the Earth offers an opportunity to infer its density profile. This has never been done, however, due to the lack of relevant data. Here we perform a neutrino-based tomography of the Earth using actual data—one-year of through-going muon atmospheric neutrino data collected by the IceCube telescope. Using only weak interactions, in a way that is completely independent of gravitational measurements, we are able to determine the mass of the Earth and its core, its moment of inertia, and to establish that the core is denser than the mantle. Our results demonstrate the feasibility of this approach to study the Earth’s internal structure, which is complementary to traditional geophysics methods. Neutrino tomography could become more competitive as soon as more statistics is available, provided that the sources of systematic uncertainties are fully under control