Threshold Electrodisintegration of ^3He

Cross sections were measured for the near-threshold electrodisintegration of ^3He at momentum transfer values of q=2.4, 4.4, and 4.7 fm^{-1}. From these and prior measurements the transverse and longitudinal response functions R_T and R_L were deduced. Comparisons are made against previously published and new non-relativistic A=3 calculations using the best available NN potentials. In general, for q<2 fm^{-1} these calculations accurately predict the threshold electrodisintegration of ^3He. Agreement at increasing q demands consideration of two-body terms, but discrepancies still appear at the highest momentum transfers probed, perhaps due to the neglect of relativistic dynamics, or to the underestimation of high-momentum wave-function components.Comment: 9 pages, 7 figures, 1 table, REVTEX4, submitted to Physical Review

Automated derivation of the adjoint of high-level transient finite element programs

In this paper we demonstrate a new technique for deriving discrete adjoint and tangent linear models of finite element models. The technique is significantly more efficient and automatic than standard algorithmic differentiation techniques. The approach relies on a high-level symbolic representation of the forward problem. In contrast to developing a model directly in Fortran or C++, high-level systems allow the developer to express the variational problems to be solved in near-mathematical notation. As such, these systems have a key advantage: since the mathematical structure of the problem is preserved, they are more amenable to automated analysis and manipulation. The framework introduced here is implemented in a freely available software package named dolfin-adjoint, based on the FEniCS Project. Our approach to automated adjoint derivation relies on run-time annotation of the temporal structure of the model, and employs the FEniCS finite element form compiler to automatically generate the low-level code for the derived models. The approach requires only trivial changes to a large class of forward models, including complicated time-dependent nonlinear models. The adjoint model automatically employs optimal checkpointing schemes to mitigate storage requirements for nonlinear models, without any user management or intervention. Furthermore, both the tangent linear and adjoint models naturally work in parallel, without any need to differentiate through calls to MPI or to parse OpenMP directives. The generality, applicability and efficiency of the approach are demonstrated with examples from a wide range of scientific applications

Photo- and Electro-Disintegration of 3He at Threshold and pd Radiative Capture

The present work reports results for: pd radiative capture observables measured at center-of-mass (c.m.) energies in the range 0--100 keV and at 2 MeV by the TUNL and Wisconsin groups, respectively; contributions to the Gerasimov-Drell-Hearn (GDH) integral in 3He from the two- up to the three-body breakup thresholds, compared to experimental determinations by the TUNL group in this threshold region; longitudinal, transverse, and interference response functions measured in inclusive polarized electron scattering off polarized 3He at excitation energies below the threshold for breakup into ppn, compared to unpolarized longitudinal and transverse data from the Saskatoon group. The calculations are based on a realistic Hamiltonian with two- and three-nucleon interactions and a realistic current operator, including one- and two-body components. The theoretical predictions obtained by including only one-body currents are in violent disagreement with data. These differences between theory and experiment are, to a large extent, removed when two-body currents are taken into account, although some rather large discrepancies remain in the c.m. energy range 0--100 keV, particularly for the pd differential cross section and tensor analyzing power at small angles, and contributions to the GDH integral. A rather detailed analysis indicates that these discrepancies have, in large part, a common origin, and can be traced back to an excess strength obtained in the theoretical calculation of the E1 reduced matrix element associated with the pd channel having L,S,J=1,1/2,3/2. It is suggested that this lack of E1 strength observed experimentally might have implications for the nuclear interaction at very low energies. Finally, the validity of the long-wavelength approximation for electric dipole transitions is discussed.Comment: 47 pages RevTex file, 10 PostScript figures, submitted to Phys. Rev.

Model Calculations for the Two-Fragment Electro-Disintegration of 4^4He

Differential cross sections for the electro-disintegration process $e + {^4He} \longrightarrow {^3H}+ p + e'$ are calculated, using a model in which the final state interaction is included by means of a nucleon-nucleus (3+1) potential constructed via Marchenko inversion. The required bound-state wave functions are calculated within the integrodifferential equation approach (IDEA). In our model the important condition that the initial bound state and the final scattering state are orthogonal is fulfilled. The sensitivity of the cross section to the input $p{^3H}$ interaction in certain kinematical regions is investigated. The approach adopted could be useful in reactions involving few cluster systems where effective interactions are not well known and exact methods are presently unavailable. Although, our Plane-Wave Impulse Approximation results exhibit, similarly to other calculations, a dip in the five-fold differential cross-section around a missing momentum of $\sim 450 MeV/c$, it is argued that this is an artifact of the omission of re-scattering four-nucleon processes.Comment: 16 pages, 6 figures, accepted for publication by Phys.Rev.

Variational methods

International audienceThis contribution presents derivative-based methods for local sensitivity analysis, called Variational Sensitivity Analysis (VSA). If one defines an output called the response function, its sensitivity to inputs variations around a nominal value can be studied using derivative (gradient) information. The main issue of VSA is then to provide an efficient way of computing gradients. This contribution first presents the theoretical grounds of VSA: framework and problem statement, tangent and adjoint methods. Then it covers pratical means to compute derivatives, from naive to more sophisticated approaches, discussing their various 2 merits. Finally, applications of VSA are reviewed and some examples are presented, covering various applications fields: oceanography, glaciology, meteorology