Measurement of the TeV atmospheric muon charge ratio with the complete OPERA data set

The OPERA detector, designed to search for νμ → ντ oscillations in the CNGS beam, is located in the underground Gran Sasso laboratory, a privileged location to study TeV-scale cosmic rays. For the analysis here presented, the detector was used to measure the atmospheric muon charge ratio in the TeV region. OPERA collected chargeseparated cosmic ray data between 2008 and 2012. More than 3 million atmospheric muon events were detected and reconstructed, among which about 110000 multiple muon bundles. The charge ratio Rμ ≡ Nμ+/Nμ− was measured separately for single and for multiple muon events. The analysis exploited the inversion of the magnet polarity which was performed on purpose during the 2012 Run. The combination of the two data sets with opposite magnet polarities allowedminimizing systematic uncertainties and reaching an accurate determination of the muon charge ratio. Data were fitted to obtain relevant parameters on the composition of primary cosmic rays and the associated kaon production in the forward fragmentation region. In the surface energy range 1–20 TeV investigated by OPERA, Rμ is well described by a parametric model including only pion and kaon contributions to themuon flux, showing no significant contribution of the prompt component. The energy independence supports the validity of Feynman scaling in the fragmentation region up to 200 TeV/nucleon primary energy.

Detection of back-to-back proton pairs in charged-current neutrino interactions with the ArgoNeuT detector in the NuMI low energy beam line

Short range nucleon-nucleon correlations in nuclei (NN SRC) carry important information on nuclear structure and dynamics. NN SRC have been extensively probed through two-nucleon knock- out reactions in both pion and electron scattering experiments. We report here on the detection of two-nucleon knock-out events from neutrino interactions and discuss their topological features as possibly involving NN SRC content in the target argon nuclei. The ArgoNeuT detector in the Main Injector neutrino beam at Fermilab has recorded a sample of 30 fully reconstructed charged current events where the leading muon is accompanied by a pair of protons at the interaction vertex, 19 of which have both protons above the Fermi momentum of the Ar nucleus. Out of these 19 events, four are found with the two protons in a strictly back-to-back high momenta configuration directly observed in the final state and can be associated to nucleon Resonance pionless mechanisms involving a pre-existing short range correlated np pair in the nucleus. Another fraction (four events) of the remaining 15 events have a reconstructed back-to-back configuration of a np pair in the initial state, a signature compatible with one-body Quasi Elastic interaction on a neutron in a SRC pair. The detection of these two subsamples of the collected (mu- + 2p) events suggests that mechanisms directly involving nucleon-nucleon SRC pairs in the nucleus are active and can be efficiently explored in neutrino-argon interactions with the LAr TPC technology.

AEgIS experiment: Towards antihydrogen beam production for antimatter gravity measurements

AEgIS (Antimatter Experiment: Gravity, Interferometry, Spectroscopy) is an experiment that aims to perform the first direct measurement of the gravitational acceleration g of antihydrogen in the Earth’s field. A cold antihydrogen beam will be produced by charge exchange reaction between cold antiprotons and positronium excited in Rydberg states. Rydberg positronium (with quantum number n between 20 and 30) will be produced by a two steps laser excitation. The antihydrogen beam, after being accelerated by Stark effect, will fly through the gratings of a moir´e deflectometer. The deflection of the horizontal beam due to its free fall will be measured by a position sensitive detector. It is estimated that the detection of about 103 antihydrogen atoms is required to determine the gravitational acceleration with a precision of 1%. In this report an overview of the AEgIS experiment is presented and its current status is described. Details on the production of slow positronium and its excitation with lasers are discussed.

Dispersion relation for hadronic light-by-light scattering: theoretical foundations

In this paper we make a further step towards a dispersive description of the hadronic light-by-light (HLbL) tensor, which should ultimately lead to a data-driven evaluation of its contribution to (g − 2) μ . We first provide a Lorentz decomposition of the HLbL tensor performed according to the general recipe by Bardeen, Tung, and Tarrach, generalizing and extending our previous approach, which was constructed in terms of a basis of helicity amplitudes. Such a tensor decomposition has several advantages: the role of gauge invariance and crossing symmetry becomes fully transparent; the scalar coefficient functions are free of kinematic singularities and zeros, and thus fulfill a Mandelstam double-dispersive representation; and the explicit relation for the HLbL contribution to (g − 2) μ in terms of the coefficient functions simplifies substantially. We demonstrate explicitly that the dispersive approach defines both the pion-pole and the pion-loop contribution unambiguously and in a model-independent way. The pion loop, dispersively defined as pion-box topology, is proven to coincide exactly with the one-loop scalar QED amplitude, multiplied by the appropriate pion vector form factors.

Progress in Simulations with Twisted Mass Fermions at the Physical Point

In this contribution, results from Nf = 2 lattice QCD simulations at one lattice spacing using twisted mass fermions with a clover term at the physical pion mass are presented. The mass splitting between charged and neutral pions (including the disconnected contribution) is shown to be around 20(20) MeV. Further, a first measurement using the clover twisted mass action of the average momentum fraction of the pion is given. Finally, an analysis of pseudoscalar meson masses and decay constants is presented involving linear interpolations in strange and charm quark masses. Matching to meson mass ratios allows the calculation of quark mass ratios: ms=ml = 27:63(13), mc=ml = 339:6(2:2) and mc=ms = 12:29(10). From this mass matching the quantities fK = 153:9(7:5) MeV, fD = 219(11) MeV, fDs = 255(12) MeV and MDs = 1894(93) MeV are determined without the application of finite volume or discretization artefact corrections and with errors dominated by a preliminary estimate of the lattice spacing.

Analysis of Synthetic Diamond Wafer Interferograms Using a Parallelized Simulated Annealing Algorithm

Diamonds are known for both their beauty and their durability. Jefferson National Lab in Newport News, VA has found a way to utilize the diamond's strength to view the beauty of the inside of the atomic nucleus with the hopes of finding exotic forms of matter. By firing very fast electrons at a diamond sheet no thicker than a human hair, high energy particles of light known as photons are produced with a high degree of polarization that can illuminate the constituents of the nucleus known as quarks. The University of Connecticut Nuclear Physics group has responsibility for crafting these extremely thin, high quality diamond wafers. These wafers must be cut from larger stones that are about the size of a human finger, and then carefully machined down to the final thickness. The thinning of these diamonds is extremely challenging, as the diamond's greatest strength also becomes its greatest weakness. The Connecticut Nuclear Physics group has developed a novel technique to assist industrial partners in assessing the quality of the final machining steps, using a technique based on laser interferometry. The images of the diamond surface produced by the interferometer encode the thickness and shape of the diamond surface in a complex way that requires detailed analysis to extract. We have developed a novel software application to analyze these images based on the method of simulated annealing. Being able to image the surface of these diamonds without requiring costly X-ray diffraction measurements allows rapid feedback to the industrial partners as they refine their thinning techniques. Thus, by utilizing a material found to be beautiful by many, the beauty of nature can be brought more clearly into view.

Development and Characterization of an in vitro Four-species Anaerobic Dental Biofilm Model

Dental caries is the most common chronic disease worldwide. It is characterized by the demineralization of tooth enamel caused by acid produced by cariogenic dental bacteria growing on tooth surfaces, termed bacterial biofilms. Cariogenesis is a complex biological process that is influence by multiple factors and is not attributed to a sole causative agent. Instead, caries is associated with multispecies microbial biofilm communities composed of some bacterial species that directly influence the development of a caries lesion and other species that are seemingly benign but must contribute to the community in an uncharacterized way. Clinical analysis of dental caries and its microbial populations is challenging due to many factors including low sensitivity of clinical measurement tools, variability in saliva chemistry, and variation in the microbiota. Our laboratory has developed an in vitro anaerobic biofilm model for dental carries to facilitate both clinical and basic research-based analyses of the multispecies dynamics and individual factors that contribute to cariogenicity. The rational for development of this system was to improve upon the current models that lack key elements. This model places an emphasis on physiological relevance and ease of maintenance and reproducibility. The uniqueness of the model is based on integrating four critical elements: 1) a biofilm community composed of four distinct and representative species typically associated with dental caries, 2) a semi-defined synthetic growth medium designed to mimic saliva, 3) physiologically relevant biofilm growth substrates, and 4) a novel biofilm reactor device designed to facilitate the maintenance and analysis. Specifically, human tooth sections or hydroxyapatite discs embedded into poly(methyl methacrylate) (PMMA) discs are incubated for an initial 24 hr in a static inverted removable substrate (SIRS) biofilm reactor at 37°C under anaerobic conditions in artificial saliva (CAMM) without sucrose in the presence of 1 X 106 cells/ml of each Actinomyces odontolyticus, Fusobacterium nucleatum, Streptococcus mutans, and Veillonella dispar. During days 2 and 3 the samples are maintained continually in CAMM with various exposures to 0.2% sucrose; all of the discs are transferred into fresh medium every 24 hr. To validate that this model is an appropriate in vitro representation of a caries-associated multispecies biofilm, research aims were designed to test the following overarching hypothesis: an in vitro anaerobic biofilm composed of four species (S. mutans, V. dispar, A. odontolyticus, and F. nucleatum) will form a stable biofilm with a community profile that changes in response to environmental conditions and exhibits a cariogenic potential. For these experiments the biofilms as described above were exposed on days 2 and 3 to either CAMM lacking sucrose (no sucrose), CAMM with 0.2% sucrose (constant sucrose), or were transferred twice a day for 1 hr each time into 0.2% sucrose (intermittent sucrose). Four types of analysis were performed: 1) fluorescence microscopy of biofilms stained with Syto 9 and hexidium idodine to determine the biofilm architecture, 2) quantitative PCR (qPCR) to determine the cell number of each species per cm2, 3) vertical scanning interferometry (VSI) to determine the cariogenic potential of the biofilms, and 4) tomographic pH imaging using radiometric fluorescence microscopy after exposure to pH sensitive nanoparticles to measure the micro-environmental pH. The qualitative and quantitative results reveal the expected dynamics of the community profile when exposed to different sucrose conditions and the cariogenic potential of this in vitro four-species anaerobic biofilm model, thus confirming its usefulness for future analysis of primary and secondary dental caries.

Digital elevation model (DEM) of the ice sheet in western Dronning Maud Land, Antarctica

In this paper, a new digital elevation model (DEM) is derived for the ice sheet in western Dronning Maud Land, Antarctica. It is based on differential interferometric synthetic aperture radar (SAR) from the European Remote Sensing 1/2 (ERS-1/2) satellites, in combination with ICESat's Geoscience Laser Altimeter System (GLAS). A DEM mosaic is compiled out of 116 scenes from the ERS-1 ice phase in 1994 and the ERS-1/2 tandem mission between 1996 and 1997 with the GLAS data acquired in 2003 that served as ground control. Using three different SAR processors, uncertainties in phase stability and baseline model, resulting in height errors of up to 20 m, are exemplified. Atmospheric influences at the same order of magnitude are demonstrated, and corresponding scenes are excluded. For validation of the DEM mosaic, covering an area of about 130,000 km**2 on a 50-m grid, independent ICESat heights (2004-2007), ground-based kinematic GPS (2005), and airborne laser scanner data (ALS, 2007) are used. Excluding small areas with low phase coherence, the DEM differs in mean and standard deviation by 0.5 +/- 10.1, 1.1 +/- 6.4, and 3.1 +/- 4.0 m from ICESat, GPS, and ALS, respectively. The excluded data points may deviate by more than 50 m. In order to suppress the spatially variable noise below a 5-m threshold, 18% of the DEM area is selectively averaged to a final product at varying horizontal spatial resolution. Apart from mountainous areas, the new DEM outperforms other currently available DEMs and may serve as a benchmark for future elevation models such as from the TanDEM-X mission to spatially monitor ice sheet elevation.

TanDEM-X ice thickness dataset for Wilkins Ice Shelf, Antarctic Peninsula, with link to GeoTIFF

The dataset shows the ice thickness over Wilkins Ice Shelf, Antarctic Peninsula derived from TanDEM-X Interferometry. The data has been acquired between June and August 2012. The TanDEM-X heights have been linked to CryoSAT-2 heights (V. Helm) from the respective time stamp. Elevations have been transformed from WGS84 ellipsoidal heights to the EGM2008 geoid. The ice shelf thickness was estimated assuming hydrostatic equilibrium and a mean ice density of 915 kg/m³.

Accumulation, velocities, fluxes and mass balance of glacial tributaries to the Lambert Glacier-Amery Ice Shelf system, Antarctica

By incorporating recently available remote sensing data, we investigated the mass balance for all individual tributary glacial basins of the Lambert Glacier-Amery Ice Shelf system, East Antarctica. On the basis of the ice flow information derived from SAR interferometry and ICESat laser altimetry, we have determined the spatial configuration of eight tributary drainage basins of the Lambert-Amery glacial system. By combining the coherence information from SAR interferometry and the texture information from SAR and MODIS images, we have interpreted and refined the grounding line position. We calculated ice volume flux of each tributary glacial basin based on the ice velocity field derived from Radarsat three-pass interferometry together with ice thickness data interpolated from Australian and Russian airborne radio echo sounding (RES) surveys and inferred from ICESat laser altimetry data. Our analysis reveals that three tributary basins have a significant net positive imbalance, while five other subbasins are slightly positive or close to zero balance. Overall, in contrast to previous studies, we find that the grounded ice in Lambert Glacier-Amery Ice Shelf system has a positive mass imbalance of 22.9 ± 4.4 Gt/a. The net basal melting for the entire Amery Ice Shelf is estimated to be 27.0 ± 7.0 Gt/a. The melting rate decreases rapidly from the grounding zone to the ice shelf front. Significant basal refreezing is detected in the downstream section of the ice shelf. The mass balance estimates for both the grounded ice sheet and the ice shelf mass differ substantially from other recent estimates.

Geodetic and satellite derived ice surface velocities of Vestfonna ice cap

During 2007 we launched a geodetic campaign on the Svalbard ice cap Vestfonna in order to estimate the velocity field of the ice cap. This was done within the frame of the IPY project KINNVIKA. We present here the velocity measurements derived from our campaigns 2007-2010 and compare the geodetic measurements against InSAR velocity fields from satellite platforms from 1995/96 and 2008. We find the spatial distribution of ice speeds from the InSAR is in good agreement within the uncertainty limits with our geodetic measurements. We observe no clear indication of seasonal ice speed differences, but we find a speed-up of the outlet glacier Franklinbreen between the InSAR campaigns, and speculate the outlet is having a surge phase.

Atmospheric, oceanic and hydrological polar motion excitation mechanisms based on a combination of geometric and gravimetric space observations

The goal of our study is to determine accurate time series of geophysical Earth rotation excitations to learn more about global dynamic processes in the Earth system. For this purpose, we developed an adjustment model which allows to combine precise observations from space geodetic observation systems, such as Satellite Laser Ranging (SLR), Global Navigation Satellite Systems (GNSS), Very Long Baseline Interferometry (VLBI), Doppler Orbit determination and Radiopositioning Integrated on Satellite (DORIS), satellite altimetry and satellite gravimetry in order to separate geophysical excitation mechanisms of Earth rotation. Three polar motion time series are applied to derive the polar motion excitation functions (integral effect). Furthermore we use five time variable gravity field solutions from Gravity Recovery and Climate Experiment (GRACE) to determine not only the integral mass effect but also the oceanic and hydrological mass effects by applying suitable filter techniques and a land-ocean mask. For comparison the integral mass effect is also derived from degree 2 potential coefficients that are estimated from SLR observations. The oceanic mass effect is also determined from sea level anomalies observed by satellite altimetry by reducing the steric sea level anomalies derived from temperature and salinity fields of the oceans. Due to the combination of all geodetic estimated excitations the weaknesses of the individual processing strategies can be reduced and the technique-specific strengths can be accounted for. The formal errors of the adjusted geodetic solutions are smaller than the RMS differences of the geophysical model solutions. The improved excitation time series can be used to improve the geophysical modeling.