Confronting models on cosmic ray interactions with particle physics at LHC energies

Inelastic pp collisions are dominated by soft (low momentum transfer) physics where perturbative QCD cannot be fully applied. A deep understanding of both soft and semi-hard processes is crucial for predictions of minimum bias and underlying events of the now coming on line pp Large Hadron Collider (LHC). Moreover, the interaction of cosmic ray particles entering in the atmosphere is extremely sensitive to these soft processes and consequently cannot be formulated from first principles. Because of this, air shower analyses strongly rely on hadronic interaction models, which extrapolate collider data several orders of magnitude. A comparative study of Monte Carlo simulations of pp collisions (at the LHC center-of-mass energy ~ 14 TeV) using the most popular hadronic interaction models for ultrahigh energy cosmic ray (SIBYLL and QGSJET) and for collider physics (the PYTHIA multiparton model) is presented. The most relevant distributions are studied including those observables from diffractive events with the aim of discriminating between the different models.Comment: 8 pages revtex, 8 figures, added reference

Influence of shower fluctuations and primary composition on studies of the shower longitudinal development

We study the influence of shower fluctuations, and the possible presence of different nuclear species in the primary cosmic ray spectrum, on the experimental determination of both shower energy and the proton air inelastic cross section from studies of the longitudinal development of atmospheric showers in fluorescence experiments. We investigate the potential of track length integral and shower size at maximum as estimators of shower energy. We find that at very high energy (~10^19-10^20 eV) the error of the total energy assignment is dominated by the dependence on the hadronic interaction model, and is of the order of 5%. At lower energy (~10^17-10^18 eV), the uncertainty of the energy determination due to the limited knowledge of the primary cosmic ray composition is more important. The distribution of depth of shower maximum is discussed as a measure of the proton-air cross section. Uncertainties in a possible experimental measurement of this cross section introduced by intrinsic shower fluctuations, the model of hadronic interactions, and the unknown mixture of primary nuclei in the cosmic radiation are numerically evaluated.Comment: 12 pages, 11 figures, 4 table

Hazard ranking method for populations exposed to arsenic in private water supplies: relation to bedrock geology

Approximately one million people in the UK are served by private water supplies (PWS) where main municipal water supply system connection is not practical or where PWS is the preferred option. Chronic exposure to contaminants in PWS may have adverse effects on health. South West England is an area with elevated arsenic concentrations in groundwater and over 9000 domestic dwellings here are supplied by PWS. There remains uncertainty as to the extent of the population exposed to arsenic (As), and the factors predicting such exposure. We describe a hazard assessment model based on simplified geology with the potential to predict exposure to As in PWS. Households with a recorded PWS in Cornwall were recruited to take part in a water sampling programme from 2011 to 2013. Bedrock geologies were aggregated and classified into nine Simplified Bedrock Geological Categories (SBGC), plus a cross-cutting “mineralized” area. PWS were sampled by random selection within SBGCs and some 508 households volunteered for the study. Transformations of the data were explored to estimate the distribution of As concentrations for PWS by SBGC. Using the distribution per SBGC, we predict the proportion of dwellings that would be affected by high concentrations and rank the geologies according to hazard. Within most SBGCs, As concentrations were found to have log-normal distributions. Across these areas, the proportion of dwellings predicted to have drinking water over the prescribed concentration value (PCV) for As ranged from 0% to 20%. From these results, a pilot predictive model was developed calculating the proportion of PWS above the PCV for As and hazard ranking supports local decision making and prioritization. With further development and testing, this can help local authorities predict the number of dwellings that might fail the PCV for As, based on bedrock geology. The model presented here for Cornwall could be applied in areas with similar geologies. Application of the method requires independent validation and further groundwater-derived PWS sampling on other geological formation

Multiparton Interactions in Photoproduction at HERA

The high energy photoproduction of jets is being observed at the ep collider, HERA. It may be that the HERA centre-of-mass energy is sufficiently large that the production of more than one pair of jets per ep collision becomes possible, owing to the large number density of the probed gluons. We construct a Monte Carlo model of such multiparton interactions and study their effects on a wide range of physical observables. The conclusion is that multiple interactions could have very significant effects upon the photoproduction final state and that this would for example make extractions of the gluon density in the photon rather difficult. Total rates for the production of many (i.e. > 2) jets could provide direct evidence for the presence of multiple interactions, although parton showering and hadronization significantly affect low transverse energy jets.Comment: 21 pages, 8 figures include

Quark initiated coherent diffractive production of muon pair and W boson at hadron colliders

The large transverse momentum muon pair and W boson productions in the quark initiated coherent diffractive processes at hadron colliders are discussed under the framework of the two-gluon exchange parametrization of the Pomeron model. In this approach, the production cross sections are related to the small-x off-diagonal gluon distribution and the large-x quark distribution in the proton (antiproton). By approximating the off-diagonal gluon distribution by the usual gluon distribution function, we estimate the production rates of these processes at the Fermilab Tevatron.Comment: 11pages, 6 PS figures, to appear in PR

New physics, the cosmic ray spectrum knee, and pppp cross section measurements

We explore the possibility that a new physics interaction can provide an explanation for the knee just above $10^6$ GeV in the cosmic ray spectrum. We model the new physics modifications to the total proton-proton cross section with an incoherent term that allows for missing energy above the scale of new physics. We add the constraint that the new physics must also be consistent with published $pp$ cross section measurements, using cosmic ray observations, an order of magnitude and more above the knee. We find that the rise in cross section required at energies above the knee is radical. The increase in cross section suggests that it may be more appropriate to treat the scattering process in the black disc limit at such high energies. In this case there may be no clean separation between the standard model and new physics contributions to the total cross section. We model the missing energy in this limit and find a good fit to the Tibet III cosmic ray flux data. We comment on testing the new physics proposal for the cosmic ray knee at the Large Hadron Collider.Comment: 17 pages, 4 figure

On the fourth-order accurate compact ADI scheme for solving the unsteady Nonlinear Coupled Burgers' Equations

The two-dimensional unsteady coupled Burgers' equations with moderate to severe gradients, are solved numerically using higher-order accurate finite difference schemes; namely the fourth-order accurate compact ADI scheme, and the fourth-order accurate Du Fort Frankel scheme. The question of numerical stability and convergence are presented. Comparisons are made between the present schemes in terms of accuracy and computational efficiency for solving problems with severe internal and boundary gradients. The present study shows that the fourth-order compact ADI scheme is stable and efficient

Diffractive light quark jet production at hadron colliders in the two-gluon exchange model

Massless quark and antiquark jet production at large transverse momentum in the coherent diffractive processes at hadron colliders is calculated in the two-gluon exchange parametrization of the Pomeron model. We use the helicity amplitude method to calculate the cross section formula. We find that for the light quark jet production the diffractive process is related to the differential off-diagonal gluon distribution function in the proton. We estimate the production rate for this process at the Fermilab Tevatron by approximating the off-diagonal gluon distribution function by the usual diagonal gluon distribution in the proton. And we find that the cross sections for the diffractive light quark jet production and the charm quark jet production are in the same order of magnitude. We also use the helicity amplitude method to calculate the diffractive charm jet production at hadron colliders, by which we reproduce the leading logarithmic approximation result of this process we previously calculated.Comment: 15 pages, 4 PS figures, Revte

Total photoproduction cross-section at very high energy

In this paper we apply to photoproduction total cross-section a model we have proposed for purely hadronic processes and which is based on QCD mini-jets and soft gluon re-summation. We compare the predictions of our model with the HERA data as well as with other models. For cosmic rays, our model predicts substantially higher cross-sections at TeV energies than models based on factorization but lower than models based on mini-jets alone, without soft gluons. We discuss the origin of this difference.Comment: 13 pages, 9 figures. Accepted for publication in EPJC. Changes concern added references, clarifications of the Soft Gluon Resummation method used in the paper, and other changes requested by the Journal referee which do not change the results of the original versio