EasyVVUQ: A library for verification, validation and uncertainty quantification in high performance computing

EasyVVUQ is an open source Python library (https://github.com/UCL-CCS/EasyVVUQ) designed to facilitate verification, validation and uncertainty quantification (VVUQ) for a wide variety of simulations. The goal of EasyVVUQ is to make it as easy as possible to implement advanced VVUQ techniques for existing application codes or workflows. Our aim is to expose these features in an accessible way for users of scientific software, in particular for simulation codes running on high performance computers

Winter Wheat Yield Response to Plant Density as a Function of Yield Environment and Tillering Potential: A Review and Field Studies

Wheat (Triticum aestivum L.) grain yield response to plant density is inconsistent, and the mechanisms driving this response are unclear. A better understanding of the factors governing this relationship could improve plant density recommendations according to specific environmental and genetics characteristics. Therefore, the aims of this paper were to: i) execute a synthesis-analysis of existing literature related to yield-plant density relationship to provide an indication of the need for different agronomic optimum plant density (AOPD) in different yield environments (YEs), and ii) explore a data set of field research studies conducted in Kansas (USA) on yield response to plant density to determine the AOPD at different YEs, evaluate the effect of tillering potential (TP) on the AOPD, and explain changes in AOPD via variations in wheat yield components. Major findings of this study are: i) the synthesis-analysis portrayed new insights of differences in AOPD at varying YEs, reducing the AOPD as the attainable yield increases (with AOPD moving from 397 pl m-2 for the low YE to 191 pl m-2 for the high YE); ii) the field dataset confirmed the trend observed in the synthesis-analysis but expanded on the physiological mechanisms underpinning the yield response to plant density for wheat, mainly highlighting the following points: a) high TP reduces the AOPD mainly in high and low YEs, b) at canopy-scale, both final number of heads and kernels per square meter were the main factors improving yield response to plant density under high TP, c) under varying YEs, at per-plant-scale, a compensation between heads per plant and kernels per head was the main factor contributing to yield with different TP.Fil: Bastos, Leonardo M.. Kansas State University; Estados UnidosFil: Carciochi, Walter Daniel. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Mar del Plata; Argentina. Kansas State University; Estados UnidosFil: Lollato, Romulo P.. Kansas State University; Estados UnidosFil: Jaenisch, Brent R.. Kansas State University; Estados UnidosFil: Rezende, Caio R.. Kansas State University; Estados UnidosFil: Schwalbert, Rai. Kansas State University; Estados UnidosFil: Vara Prasad, P.V.. Kansas State University; Estados UnidosFil: Zhang, Guorong. Kansas State University; Estados UnidosFil: Fritz, Allan K.. Kansas State University; Estados UnidosFil: Foster, Chris. John Deer; Estados UnidosFil: Wright, Yancy. John Deer; Estados UnidosFil: Young, Steven. John Deer; Estados UnidosFil: Bradley, Pauley. John Deer; Estados UnidosFil: Ciampitti, Ignacio Antonio. Kansas State University; Estados Unido

Electron Scattering From High-Momentum Neutrons in Deuterium

We report results from an experiment measuring the semi-inclusive reaction $d(e,e'p_s)$ where the proton $p_s$ is moving at a large angle relative to the momentum transfer. If we assume that the proton was a spectator to the reaction taking place on the neutron in deuterium, the initial state of that neutron can be inferred. This method, known as spectator tagging, can be used to study electron scattering from high-momentum (off-shell) neutrons in deuterium. The data were taken with a 5.765 GeV electron beam on a deuterium target in Jefferson Laboratory's Hall B, using the CLAS detector. A reduced cross section was extracted for different values of final-state missing mass $W^{*}$, backward proton momentum $\vec{p}_{s}$ and momentum transfer $Q^{2}$. The data are compared to a simple PWIA spectator model. A strong enhancement in the data observed at transverse kinematics is not reproduced by the PWIA model. This enhancement can likely be associated with the contribution of final state interactions (FSI) that were not incorporated into the model. A ``bound neutron structure function'' $F_{2n}^{eff}$ was extracted as a function of $W^{*}$ and the scaling variable $x^{*}$ at extreme backward kinematics, where effects of FSI appear to be smaller. For $p_{s}>400$ MeV/c, where the neutron is far off-shell, the model overestimates the value of $F_{2n}^{eff}$ in the region of $x^{*}$ between 0.25 and 0.6. A modification of the bound neutron structure function is one of possible effects that can cause the observed deviation.Comment: 33 pages RevTeX, 9 figures, to be submitted to Phys. Rev. C. Fixed 1 Referenc

Observation of exclusive DVCS in polarized electron beam asymmetry measurements

We report the first results of the beam spin asymmetry measured in the reaction e + p -> e + p + gamma at a beam energy of 4.25 GeV. A large asymmetry with a sin(phi) modulation is observed, as predicted for the interference term of Deeply Virtual Compton Scattering and the Bethe-Heitler process. The amplitude of this modulation is alpha = 0.202 +/- 0.028. In leading-order and leading-twist pQCD, the alpha is directly proportional to the imaginary part of the DVCS amplitude.Comment: 6 pages, 5 figure

A Bayesian analysis of pentaquark signals from CLAS data

We examine the results of two measurements by the CLAS collaboration, one of which claimed evidence for a $\Theta^{+}$ pentaquark, whilst the other found no such evidence. The unique feature of these two experiments was that they were performed with the same experimental setup. Using a Bayesian analysis we find that the results of the two experiments are in fact compatible with each other, but that the first measurement did not contain sufficient information to determine unambiguously the existence of a $\Theta^{+}$. Further, we suggest a means by which the existence of a new candidate particle can be tested in a rigorous manner.Comment: 5 pages, 3 figure

Complete measurement of three-body photodisintegration of 3He for photon energies between 0.35 and 1.55 GeV

The three-body photodisintegration of 3He has been measured with the CLAS detector at Jefferson Lab, using tagged photons of energies between 0.35 GeV and 1.55 GeV. The large acceptance of the spectrometer allowed us for the first time to cover a wide momentum and angular range for the two outgoing protons. Three kinematic regions dominated by either two- or three-body contributions have been distinguished and analyzed. The measured cross sections have been compared with results of a theoretical model, which, in certain kinematic ranges, have been found to be in reasonable agreement with the data.Comment: 22 pages, 25 eps figures, 2 tables, submitted to PRC. Modifications: removed 2 figures, improvements on others, a few minor modifications to the tex

A Kinematically Complete Measurement of the Proton Structure Function F2 in the Resonance Region and Evaluation of Its Moments

We measured the inclusive electron-proton cross section in the nucleon resonance region (W < 2.5 GeV) at momentum transfers Q**2 below 4.5 (GeV/c)**2 with the CLAS detector. The large acceptance of CLAS allowed for the first time the measurement of the cross section in a large, contiguous two-dimensional range of Q**2 and x, making it possible to perform an integration of the data at fixed Q**2 over the whole significant x-interval. From these data we extracted the structure function F2 and, by including other world data, we studied the Q**2 evolution of its moments, Mn(Q**2), in order to estimate higher twist contributions. The small statistical and systematic uncertainties of the CLAS data allow a precise extraction of the higher twists and demand significant improvements in theoretical predictions for a meaningful comparison with new experimental results.Comment: revtex4 18 pp., 12 figure

eta-prime photoproduction on the proton for photon energies from 1.527 to 2.227 GeV

Differential cross sections for the reaction gamma p -> eta-prime p have been measured with the CLAS spectrometer and a tagged photon beam with energies from 1.527 to 2.227 GeV. The results reported here possess much greater accuracy than previous measurements. Analyses of these data indicate for the first time the coupling of the etaprime N channel to both the S_11(1535) and P_11(1710) resonances, known to couple strongly to the eta N channel in photoproduction on the proton, and the importance of j=3/2 resonances in the process.Comment: 6 pages, 3 figure

Measurement of the Deuteron Structure Function F2 in the Resonance Region and Evaluation of Its Moments

Inclusive electron scattering off the deuteron has been measured to extract the deuteron structure function F2 with the CEBAF Large Acceptance Spectrometer (CLAS) at the Thomas Jefferson National Accelerator Facility. The measurement covers the entire resonance region from the quasi-elastic peak up to the invariant mass of the final-state hadronic system W~2.7 GeV with four-momentum transfers Q2 from 0.4 to 6 (GeV/c)^2. These data are complementary to previous measurements of the proton structure function F2 and cover a similar two-dimensional region of Q2 and Bjorken variable x. Determination of the deuteron F2 over a large x interval including the quasi-elastic peak as a function of Q2, together with the other world data, permit a direct evaluation of the structure function moments for the first time. By fitting the Q2 evolution of these moments with an OPE-based twist expansion we have obtained a separation of the leading twist and higher twist terms. The observed Q2 behaviour of the higher twist contribution suggests a partial cancellation of different higher twists entering into the expansion with opposite signs. This cancellation, found also in the proton moments, is a manifestation of the "duality" phenomenon in the F2 structure function