Low energy neutron production by inverse β decay in metallic hydride surfaces

It has recently been argued that inverse-β nuclear transmutations might occur at an impressively high rate in a thin layer at a metallic hydride surface under specific conditions. In this note we present a calculation of the transmutation rate, which shows that there is little room for such a remarkable effect.


Introduction and main results
In a number of papers [1][2][3][4][5], A. Widom and collaborators consider the intriguing possibility that low energy neutrons are produced on metallic hydride surfaces by the capture of electrons by protons and the subsequent inverse-β reaction: This process cannot, of course, take place in a hydrogen atom in vacuum; it would have a negative Q-value (Q = M p + m e − M n < 0). The authors advance the hypothesis that mass renormalization from the interaction of electrons in atoms with some electromagnetic external field (e.g. due to a laser source) may increase the effective electron mass, so as to make the reaction (1) possible. The tilde sign on the electron symbol underscores that the electron bound to the proton is a 'dressed' one, with an effective mass denoted by A central element of the analysis of this very intriguing effect is the decay rate of the boundẽp state which is given in Ref. [5] as: 1 Natural units are used throughout, we approximate by unity the cosine of the Cabibbo angle, λ = g A /g V ≈ 1.25 is the axial to vector coupling in neutron β decay. The threshold value for process (1) to occur at all is With the above numerical values and β ≈ 20, the authors find What is surprising in formula (3) is the absence of the factor |ψ(0)| 2 , which gives the probability for finding the electron and the proton at the same point. 2 This factor is a mandatory consequence of the local Fermi lagrangian density: With the electron and proton fields computed at the same point, the factor |ψ(0)| 2 indicates that it must be difficult for the lagrangian to induce transitions from an initial state with, say, the proton localized in our lab and the electron on Mars. We present a calculation of the rate of (1) done in two independent ways. One is to consider the decay of theẽp bound states with spin 0 and 1, relating the amplitude to the basic lagrangian (6) in a way analogous to the quark model calculation of pion decay made time ago by Van Royen and Weisskopf (VRW) [6].
The second is to use the time-honored formula: where σ is the unpolarized cross section for process (1). This is the basic formula in the theory of electron K-capture, also used [7] to derive the decay rate of positronium from the annihilation cross section of an electron-positron pair into photons.
In the following, we give the details of the calculation of Γ . In the end, we comment on the high values of β considered in [1][2][3][4][5].

Decay rates of S-wave bound states
Neutron emission can be described as arising from the decay of the two, S-wave ground states, H 0 and H 1 , with total spin equal to zero and one, respectively. We may describe this decay with two phenomenological parameters, f 0,1 according to: The Fermi constant has been inserted for convenience, H 0 and H i 1 are field operators describing the annihilation of either H state and we describe the creation of the neutrino with the antineutrino field, ν c . The situation is entirely similar to the decay π − → μ − ν c , described by the phenomenological parameter f π . 3 To connect to the VRW formulas, it is convenient to rewrite the Fermi lagrangian so as to have the annihilation operators for the electron and the proton in the same Dirac bilinear (similarly to the quark-antiquark fields in pion decay). This is done by first introducing the positron and antineutrino fields in (6) and then making a Fierz transformation. To wit: 4 Following VRW [6], we define the correctly normalized H state according to r and s are spin indices, φ(r, s) are the Clebsch-Gordan coefficients appropriate to the spin of H , (a p,e ) † are the proton and electron creation operators and f (p) the momentum space wave function, with The parameter f 0 is obtained by comparing (15) with (16), e.g.: Expanding the fields in normal modes, spinors and gamma matrices combine with the appropriate Clebsches, and we obtain (the √ 2 factor comes from the Clebsch, the overall signs of f 0,1 are irrelevant): with the x-space wave function given by Similarly: Direct calculation of the spin-averaged rates gives and the inclusive rate is obtained from as anticipated in (8).

Inclusive rate from the unpolarized cross section
The cross section for e − + p → ν + n is [8,9]: The transmutation rate thus obtained via Eq. (7) confirms the result in (8).

The value of β
Values of β of the order of or even larger than 20 are certainly unusual in condensed matter physics, especially for bound electrons. An estimate of β is given in [1][2][3][4][5] where A is given in terms of the plasma frequency Ω p of the protons: and |u| 2 is the r.m.s. displacement of the protons.

Conclusions
A correct calculation gives a neutron production rate from (1) about 300 times smaller than estimated in [1][2][3][4][5], for the value of the mass renormalization factor β ≈ 20 considered there. In turn, it is questionable that values of β can be realized, in particular for bound electrons, so large as to give rise to useful nuclear transmutation rates. A more detailed analysis of the attainable values of β is needed to obtain more definite conclusions on this interesting phenomenon, should it exist at all.