%Neutron Bands Made of Excited Neutron States of Nuclei on the Lattice Mediated by Protons in Transition-Metal Hydrides (to be

1. Introduction

The nuclear structure of isolated nuclei ^A_ZX has been thoroughly investigated in about sixty years since the discovery of the atomic nucleus in 1911 in order to achieve fundamental understanding in the energy region up to several hundred MeV^1,2). The global features of the exited levels of nucleons and their energy distribution seem to be fairly well described by the Fermi gas mode, while the results have had been mainly confined to light nuclei and a quantitative analysis is plagued with difficulties in the description of the reaction mechanism.¹⁾ This is true even now especially for excited levels with energies very close to the zero level; which corresponds to the neutron level with a binding energy of zero in the nucleus ^A_ZX, or to the state where a neutron and the separated nucleus ^A-1_ZX remain still. (We use this energy standard in this paper unless otherwise stated.)

Therefore, it is interesting to investigate some phenomena that are directly related with the excited levels of nucleons at around zero energy in medium and heavy nuclei.

In this paper, these features of excited states of nuclei in solids are semi-quantitatively investigated on the knowledge of nuclear structures established in nuclear physics and apply them to cold fusion phenomenon (CFP). We use the Fermi gas model for nucleons in a nucleus throughout this work.

2. Excited States of Neutrons and its Density of States in Medium and Heavy Nuclei

It is a common knowledge in nuclear physics that average properties of the excitation spectrum are given by the Fermi gas model as a result of dominance of the particle degrees of freedom over the number of collective modes.¹⁾

In the Fermi gas model, nucleons in a lattice nucleus at a_i is treated as independent particles and their quantum states ψ_{n}(x, a_i) are specified by quantum numbers {n} ≡ (n, l, m, s);

ψ_{n}(x, a_i) =ψ_{nlms}(x – a_i,σ). (1)

The wave function of a neutron in a nucleus ^A_ZX, however, extends far away from the nucleus when the energy E of the state is less than but close to zero and then the wave function outsides the nucleus is approximated by

ψ_{nlms}(x – a_i,σ) = c_i e ^–^η^{|r
– ai|} Y_l,m(θ_i,φ_i)χ_s(σ), (2)

whereη≡η(|E|) is a damping factor of the radial wave function depending on the energy assumed for simplicity to be independent of quantum numbers, and (θ_i,φ_i) are angles measured from the lattice point a_i. In the following treatment, we use the wave function (1) until we need the wave function (2).

The result of the calculation of the total level density for the Fermi gas in a nucleus ^A_ZX is given as:¹⁾

ρ(N,Z,ε)= (6^1/4g₀/12(g₀ε)^5/4)exp((4π²/6)g₀ε)^1/2

(N ≃ Z) (3)

where ε is the excitation energy measured from the ground state level and g₀ is the one-particle level density at the Fermi energy g_F, representing the sum of the proton and neutron level densities

g₀ ≡ g(ε_F) = (3/2)(A/ε_F), (4)

for a case Z = N = A/2. These levels seem very sharp and have fairly long lifetime, which we take as an infinite in the following treatment.

The energy range, where the above formula is applicable, is determined by a relation

ε_F /A ≪ E ≪ε_FA^1/3, (5)

whereε_F ≈ 37 MeV for heavy nuclei.¹⁾ This relation gives an energy range 0.4 – 170 MeV of applicability of the relation (3) for nuclei with mass numbers A ≅ 100.

High density of nuclear levels at high excitation energies, amounts of the order 10⁶ times higher than that corresponding to single-particle motion, has been revealed by densely spaced, sharp resonances in the slow neutron capture reactions and results in formation of the compound nucleus in a nucleus with A ≅ 100.^1,3) The figure 10⁶ will be increased further by several orders when the energy of the slow neutron capture reactions goes down to ≅ 1 eV. In the following discussion, we will take this factor as 10⁹ at its maximum suggested by experimental data for Ag in the range of 2 to 8 MeV³⁾ considering later application to Pd isotopes in the energy range up to 10 MeV.

3. Effective Potential for the Super-nuclear Interaction between Neutrons in Adjacent Lattice Nuclei of Metal Hydrides and Deuterides

In the transition-metal hydrides MeH_x, on the other hand, the crystal structure is dependent on the concentration x of hydrogen isotopes which can be introduced into the crystal lattice of the metal Me continuously until a definite limit and kept stably there (occluded).^4,5) We confine our investigation to crystals of stoichiometric compounds PdH for our object in the following treatment. In this compound, hydrogen atoms occluded in the crystal are ionized and occupy octahedral interstices having six Pd atoms each as nearest neighbors on the crystallographic axes half way of the lattice constant a. The lattice constant a of the compound PdH_x depends on the composition and that of PdH is a little larger than that of Pd crystal 3.89 Å. In the following treatment, however, we ignore the dependence of a on the composition x and use the value for Pd crystal as for the compound PdH.

Dynamical behavior of the proton occluded in transition-metal hydrides is described as a harmonic oscillator in its ground and lower excited states. The wave function, φ_p(R–b_j,σ), of a proton in a state specified by quantum numbers p ≡ (n_p,l,m,s_p) at an interstice b_j can have finite probability density at nearby lattice point at a_i, a nearest neighbor of b_j, especially when the proton is in its excited states. If we ignore mutual interaction of Z protons on different interstices, the total proton wave function may be expressed as a product of wave functions on the interstices (neglecting anti-symmetrization),

Φ_{p_α_}(X₁, X₂, ,Xz,) = Π_jφ{p_j}(R_j –b_j,σ_j), (6)

where {p_α} ≡{p₁, p₂, , p_z}.

The overlapping of the proton wave functionφ{p_j}(R_j – b_j) on the interstice b_j with a nucleon (neutron) wave functionψ{n}(r – a_i), Eq. (1), of an adjacent lattice nucleus at a_i results in the proton-neutron interaction through the nuclear force. The nuclear interaction is expressed by a potential whose form is taken, for example, as the square-well type;

V_s(r – R) = – V^(s)₀, (|r – R| < b) (7)

= 0, (|r – R > b)

where V^(s)₀ ≅ 3.5 MeV and b ≅ 2.2 ×10^–13 cm.⁶⁾ The choice of this potential out of several possible types does not make a large difference to the result for low energy phenomena we are considering in this paper.

This interaction pulls two neutron states in different lattice nuclei into coupling as shown below that we will call the "super-nuclear interaction." In the following investigation, we concentrate on excited neutrons in lattice nuclei than protons, which needs more energy to be raised to the excited levels with the same energy than neutrons due to the fact Z ≪ N. (In Pd, Z = 46 and N = 56 – 64.)

Let us consider a neutron in an excited state {n} of one of lattice nuclei. The regularity of the crystal lattice determines the coefficients of the linear combination as required by the Bloch's theorem.⁷⁾ Then in a periodic potential of lattice nuclei, a neutron in an excited state {n} of a lattice nucleus at a_i should be expressed by a Bloch function (omitting the spin part)

ψ_k(r) = Σ_ie^i(kai)ψ_{n}(r – a_i). (8)

Therefore, the total wave function of the system composed of a neutron Bloch waveψ_k(r) and z occluded protons in the state {p_α} = {p₁, p₂, , p_z,} at interstices is expressed as (omitting spin parts)

Ψ_k,{p}(r;R₁, R_{2 ,}R_z)

=ψ_k(r) Φ_{p_α_}(X₁, X₂, ,Xz,) . (9)

The total energy E_k,{p_α_} of this system in the second-order perturbation approximation is expressed as follows taking the square well potential for the nuclear interaction:

E_k,{p_α_} = E_{n,p_α_}

+Σ_k_’,i,i’,jexp(–i(ka_i – k’a_i’))v_np(ii’j), (10)

v_np(ii'j)

=Σ_p'(<np;ij|V|n'p';ij><n'p';i'j|V|np;i'j>)/(E_{n',p'} – E_{n,p}),

=Σ_{p’}_≠_{p}P∫dEρ×

(<np;ij|V|n'p';ij><n'p';i'j|V|np;i'j>)/(E+ε_p'p), (11)

E_{n,p_α_}= E_{n}^(p)₊Σ_jε_pj, V(r) = V_s(r), (12)

<np;ij|V|n'p';ij> = ∬drdR_jψ*_{n}(r – a_i)φ*{p}(R_j –b_j,)

×V_s(r – R) ψ_{n’}(r – a_i)φ_p’(R_j –b_j,), (13)

where summations over i and i' in (10) are only over the nearest neighbor lattice points a_i and a_i’ of an interstice b_j, ρ_n(E) is a density of states for neutron quantum states, ε_p'p ≡ε_p’ –ε_p, and E ≡E_{n'} – E_{n}. Further, the summation over {p’} reduces to a factor, (n_p+1)(n_p+2), the degeneracy of the energyε_np. E_{n}^(p) is an energy of a neutron in an excited stateψ_{n’}(r – a_i) in a lattice nucleus at a_i when occluded protons are in states {p_α}, and ε_pj in (12) is an energy of a proton in a stateφ_pj(R_j –b_j,) at an interstice b_j. We ignore, however, p-dependence of E_{n}^(p) hereafter in this work.

For the neutron wave function (1) in the Fermi gas model, we can describe wave functions ψ_{n’}(r – a_i) by those determined in the nuclear harmonic oscillator potential in a nucleus to calculate matrix elements (13) in the above equation (11):

ψ_nlms(r,θ,φ,σ)

= R_nl(r)Y_l^m(θ,φ)χ_s(σ), (|m| ≦ l) (14)

E_nlms = (n + 3/2)(h/2π)ω_n+ Δε_lms (15)

whereΔε_lms expresses the l·s and other coupling energies taken symbolically into consideration to distinguish energies of the states with the same n and different l, m, and s, ω_n is the circular frequency of the harmonic oscillator and Y_l^m(θ,φ) are the spherical harmonics.

In nuclei of palladium isotopes, we can use an excited neutron state 2f_7/2 as shown by shell model calculation with a Woods-Saxon potential¹⁾ for the order of magnitude estimation of (14):

ψ_2f
7/2,_s (r,θ,φ,σ)

= R₅₃(z)Y₃^m(θ,φ)χ_s(σ), (|m|≦ 3) (16)

R₅₃(z) = C_n(32/210)^1/2z^3/2(1 – (2/9)z)e^–z/2, (17)

C_n = 2(8α_n³/π)^1/4, z = 2α_nr², α_n = πm_nω_n/h,

where m_n is the mass of the neutron and ω_n
=41/A^1/3 MeV.⁸⁾

For the interstitial proton wave functionsφ_p’(R_j –b_j,) in PdH, on the other hand, we can use a wave functionφ_1d(R,Θ,Φ) in a lattice harmonic oscillator potential centered at an interstice determined by diffusion data;⁹⁾

φ_p’(R_j) =φ_nplms’(R, Θ,Φ,σ_p)

= ξ_npl(R)Y_lm (Θ,Φ)χ_s(σ_p), (|m|≦l) (18)

ε_n_plm= 2π(n_p + 3/2)hω_p, (19)

φ_1d’(R, Θ,Φ) =ξ_1d(Z)Y₂₀(Θ,Φ), (n = 2) (20)

ξ_1d(Z) = C_p(4/15)^1/2Zexp(–Z/2), (21)

C_p = 2(8α_p³/π)^1/4, Z = 2α_pR², α_p = {m_pπω_p/h}, ω_p = (K/m_p)^1/2,

or by Hermite polynomials H_n(ξ);¹⁰⁾

φ_p’(R_j –b_j, σ_p) = u_nx(x) u_ny(y) u_nz(z) χ_s(σ_p), (22)

u_nx(x) = N_nH_n(αx)exp(–(1/2)α²x²), (23)

α⁴= 4π²m_pK/h², N_n = (α/π ^1/22ⁿn!))^1/2.

where R = (R, Θ,Φ), n_p is an integer, l ≦ n_p and |m| ≦l, ε_nlm is the proton energy of the stateφ_nlm’(R), ω_p = (K/m_p)^1/2, m_p is the mass of the proton, K is the force constant, and n_i(i = x, y, z) are integers.

The proton wave functions thus determined include already effects of screening by itinerant electrons and electrons bound in atoms, and also the effect of Coulomb repulsion by lattice nuclei.

The analysis based on the diffusion data⁹⁾ showed that appropriate wave functions for a proton in the NbH is that with n = 2 in the above equation and the corresponding force constant K is given as

K_H = 1.44 ×10¹⁹ eV/m² (NbH) . (24)

We use this value for PdH to make an order of magnitude estimation in this paper.

A concrete expression of the matrix element (14) for PdH is expressed as follows using wave functions (15), (19), and others:

<2f_{7/2}1d;ij|V|2p_{3/2}2s;ij>

= – ∬drdR_jR₅₃(z_i)Y_3,0(θ_i,φ_i)χ_1d(Z_j)Y_2,0(Θ_j,Φ_j)

×V_s(r – R_j)R₅₁(z_i)Y_1,0(θ_i,φ_i )χ_2s(Z_j) Y_0,0(Θ_j,Φ_j), (25)

z_i = 2α_n|r – a_i|², Z_j = 2α_p|R_j– b_j|²,

where a_i is a nearest neighbor lattice site of an interstice b_j, K = K_H inα_p in Eq.(22), and (θ_i,φ_i) and (Θ_j,Φ_j) are angles measured from origins at a_i and b_j, respectively.

To estimate an order of magnitude of the effective potential v_np(ii'j) (11), we utilize the property of the densely spaced excited states explained before and ignore selection rules associated with single configurations. Furthermore, we put the numerator of (11) as a constant and take it as the value of the matrix element (25) for PdH.

Then, the order of magnitude of the effective potential v_np(ii'j) given in Eq.(11) is estimated as follows: the proton wave functionφ_p’(R) is slowly varying in the range of the nuclear force, and the nuclear wave functionψ_n(r) is approximated by a delta-function. Then, an order of magnitude of the matrix elements <np;ij|V|n'p';ij> is given as

|<np;ij|V|n'p';ij>|

≅ ∫ψ_n(r)* ψ_n(r)dr<V>φ_p’(R) *φ_p’(R)Ω (26)

≅ 1×{4/3}πr₀³×|u₂(x_N)|²|u₀(0)|²|u₀(0)|²

= 3.2 ×10^–14eV, (27)

where Ωis the volume of the Pd nucleus, <V> = |V₀^(s)|= 3.5 MeV (Eq.(7)), φ_p’(R) is taken as u₂(x)u₀(y)u₀(z) and x_{N =}1.95Å is the position of the lattice nucleus measured from the interstice.

Putting this value (30) into Eq.(11), we can estimate the effective potential v_np(ii'j) as a function of the principal value of the integration appeared in that equation, assuming the insensitiveness of the matrix elements to the energy:

v_np(ii'j) ≈ 1×10^–27 eV²I,

I ≡ P∫(ρ_n(E)/E)dE. (28)

We can estimate the approximate value of the integral I, taking following valuesρ_n(E) ≈10⁹keV^–1, δε≈ 10^–9 keV, and Δε≈1 keV on the assumption that single particle energy level difference is ≈ 1 keV and the level density increases to 10⁹times larger than that of single particle motion:

I ≈ (ρ_n(ε)/δε) Δε = 10¹⁵ eV^–1.

v_np(ii'j) ≈ 1 ×10^–12 ( eV). (29)

4. Tight-Binding Neutron Bands in Metal Hydrides and Deuterides

The effective super-nuclear interaction energy obtained above is used to calculate band structure of neutron energy in transition-metal hydrides that is originated in the excited states of neutrons in lattice nuclei and mediated by occluded hydrogen isotopes.

To show briefly crystal-structure dependence of the bandwidth, we will make a simplification of the super-nuclear interaction (11) between adjacent nuclei at a_i and a_i’ assuming that it depends only on the magnitude of the vector a_l ≡ a_i – a_i’.

Then, we can rewrite the total energy (10) and have energy spectrum of the neutron Bloch waves in the face centered cubic (fcc) lattice (a is the side of the lattice cube);⁷⁾

E = E_{n,p_α_} – α- 2×4γ(cos(1/4)k_ya cos(1/4)k_za + cos(1/4)k_za cos(1/4)k_xa + cos(1/4)k_xa cos(1/4)k_ya)

– 2γ(cosk_xa + cosk_ya + cosk_za) (fcc) (30)

E_{n,p_α_} = E_{n} + Σ_jε_pj,

– α = v_np(0), –γ = v_np(ii'j), (31)

The factor 2 in the third term on the right comes from the fact that nearest neighbor lattice nuclei are mediated by two protons at different interstices while next nearest ones are by only one. A characteristic of this energy band formation is the contributions from nearest neighbors ((0, ±a/2, ±a/2) etc.) and also from next nearest neighbors ((±a,0,0) etc.) to the k-dependent terms.

The neutron energy bands originating in the excited states of lattice nuclei are located below zero energy in contrast to those originating in free neutron states above zero worked out in a previous paper.¹¹⁾ The former could be called neutron valence band and the latter neutron conduction band to distinguish them in the following discussion of the nuclear reactions in solids.

Using the value of v_np(ii'j) given in (29), we obtain a semi-quantitative estimation of the valence band width Δ from Eq.(30):

Δ = 24 v_np(ii'j) ≈ 10^–8 (meV) (PdH). (32)

Thus, it is concluded that the matrix elements (25) should be 10⁵ times larger than the values estimated in (27) to substantially keep the neutron bands below zero which was determined to form in solids with a width Δ ≥ 25 meV that is not destroyed by the thermal motion of ions at room temperature. This is realized only when the neutron wave function (1) extends out as the wave function (2) from a lattice nucleus to regions where a wave function of the occluded proton (23) has a larger value by a factor 10⁵ than that at the lattice nuclei. The main term of the proton wave function relevant to this behavior is the exponential factor exp{–α²x²/2} in (23) and it gives this value at x₀ = 1.43 Å from an interstice (or 0.52 Å from a lattice point). If this behavior is coupled with an extension of the neutron wave function (2), then the neutron-proton interaction can contribute to formation of a neutron valence band with a width of Δ ≿ 25 meV.

From a point of view of the isolated nucleus treated in conventional nuclear physics, this is an unconceivable situation. While, the extension of a neutron wave function (2) far away to 0.52 Å= 5.2×10^–9 cm over the nuclear extent range of r₀ = 10^–13 cm, i.e. 10⁴times longer than r₀, is not absurd in the situation we are considering here.

As was shown by numerical calculation in a previous paper,¹¹⁾ the energy of thermal neutrons interacting with lattice nuclei by attractive nuclear force is pulled down below zero; the states of propagating waves then become quasi-localized states around lattice nuclei with a damping factor depending on the strength of the attractive interaction. The same situation is also realized from opposite direction as a limit of highest bound states as shown in Eq.(2). We consider here an s-type wave function for the state, for simplicity:

ψ_η(r- a_i) = c_iexp( –iη|r – ai|). (33)

To extend the neutron wave function to the range of λ= 5.2 ×10^–9 cm referred above, the decay constant of the state η(|E|) = 1/λ should be 1.9 ×10⁸ cm^–1 and this corresponds to an energy E:

|E| = (h²/8π²m_n)η(|E|)² = 7.4 (eV) (34)

below zero, where m_n = 1.67×10^–24 g is the neutron mass. In other words, the excited states of isolated lattice nuclei with energies of from zero to 7 eV can participate to the neutron valence band, or the neutron bands below zero, in transition-metal hydrides considered above.

If the state has less energy, i.e. far from zero, and the extension of the state is less than 5.2×10^–9 cm, the band state fails to be substantially formed even in PdH and neutrons are essentially in single particle states in isolated lattice nuclei.

5. Discussion

When there are many neutrons in a neutron band, there appear interesting features of neutron's behavior at boundaries that reflect neutrons back into the crystal; "local coherence" of neutron Bloch waves, and therefore, high densities of neutrons (neutron liquid) appear there.¹²⁾ High-density neutrons in the boundary region¹³⁾or in neutron star matters ^14,15) induce formation of "neutron drops" (or clusters of many neutrons and a few protons and corresponding electrons) in a thin neutron background. These neutron liquid (NL) and neutron drops (ND) in a thin neutron background interact with nuclei to produce new nuclear effects in the boundary region.

Scenario of the CFP will be written down as follows. The background thermal neutrons in ambience trapped in a sample of the transition-metal hydrides or deuterides are in a neutron conduction band. Their density at boundary region becomes high due to the local coherence but may be not so large and not enough to form neutron drops. The neutrons in the band, however, can reacts with nuclei in the boundary region and the reactions are the trigger reactions.^16,17)The nuclear products of the trigger reactions induce breeding reactions resulting in multiplication of the number of neutrons in the conduction band and also excitation of neutrons in lattice nuclei.

The latter effect makes possible formation of neutron valence bands (NVB) in the CF matter we are now considering. The density of neutrons in the NVB will be very large enough to form neutron liquid (NL) and neutron drops (ND) in the boundary region. The ND thus formed may be in a lattice (a Coulomb lattice) with smaller lattice constants coexisting with original crystal lattice of the transition metal with larger (≈ 10² times) lattice constants. This is a new state of solids not noticed and not observed until CFP was detected.

The NL and/or ND thus formed can give or exchange nucleons with lattice nuclei and/or with nuclei of minor elements in the boundary layer. Nuclear reactions investigated in nuclear physics in 20th century were mainly those occur in free space except rare cases of n-p cluster formations in the neutron star matter.^14,15) The nuclear reactions in surface layers where are lattice nuclei, occluded hydrogen isotopes, and high density neutrons (NL and ND) should be distinguished from those occurring in free space and treated with similar cautions to the n-p cluster formation in the neutron star matter.¹⁵⁾

The fundamental differences related with nuclear transmutations observed in CFP are possibilities of 1) nucleon exchange between NL (and/or ND) and lattice nuclei (and/or minor nuclei) and 2) energy exchange between NL (and/or ND) and nuclides in excited states in the CF matter. The former gives a possibility to generate new nuclides with largely different mass and proton numbers from lattice or minor nuclei in the CF matter and the latter gives a possibility to stabilize excited states of nuclides without emission of γ– rays.

In our treatment of experimental data sets in CFP,^{11,12,16–18)} we have applied the TNCF model to various events only using reactions where occurs absorption of a neutron by a nucleus followed by β– or α– decay or by fission to explain various products with successful results. The nuclear transmutations, however, have shown large changes of mass numbers up to several tens in the experiments showing NT_F^19–23) and recent experiment of NT_A^23–27) which needs possibility to absorb large number of neutrons or sometimes the n-p clusters simultaneously. The formation of NL and/or the neutron drops (ND) gives natural explanation of these absorptions.

As we have seen in this paper, CFP is a wide spread phenomenon including excess heat generation, three types of NT, production of light elements, ³H and ⁴He, emissions of neutrons, gammas and X-rays with various energies up to about 10 MeV, and decay-time shortenings^16,28-30) occurring in complex systems composed of transition-metal hydrides and deuterides and others at about room temperature in ambient radiation.

The events with large variety from nuclear transmutations to emissions of light particles and γ-rays are evidences of nuclear reactions occurring in surface layers of CF materials, especially transition-metal hydrides and deuterides, intermittently and sporadically. Investigating this phenomenon, we could figure out physics of CFP as neutron physics in crystals occluding hydrogen isotopes; formation of two types of neutron bands, neutron liquid and neutron drops in surface layers where appears the local coherence of neutron Bloch waves.

Knowing physics of CFP, we can explore various applications of this phenomenon ranging from production of new nuclides, remediation of hazardous radioactivity, and production of thermal energy although limited by our poor imagination at present. Really, world of application of this phenomenon will be wider exceeding our present imagination.

Acknowledgment

The author would like to express his thanks to John Dash, Makoto Takeo and Jon Warner for valuable discussions during this work. This work is supported by a grant from the New York Community Trust and Professional Development Fund for Part-time Faculty of Portland State University.

References

[1] A. Bohr and B.R. Mottelson, Nuclear Structure I, Benjamin, New York, 1969.

[2] A. Bohr and B.R. Mottelson, Nuclear Structure II, Benjamin, New York, 1975.

[3] K. Tsukada, S. Tanaka, M. Maruyama and Y. Tomita, "Energy Dependence of the Nuclear Level Density below the Neutron Binding Energy" Nuclear Physics, 78, 369 (1966).

[4] R. Bau ed., Transition Metal Hydrides, American Chemical Society, Washington, D.C., 1978.

[5] G. Alefeld and J. Voelkl ed. Hydrogen in Metals I, Springer-Verlag, Berlin, 1978.

[6] J.M. Blatt and V.F. Weisskopf, Theoretical Nuclear Physics, Chapter II, John-Wiley & Sons, New York, 1952.

[7] For instance, N.F. Mott and H. Jones, The Theory of the Properties of Metals and Alloys, Dover, New York, 1958.

[8] W.F. Hornyak, Nuclear Structure, Chapter IV, Academic Press, New York, 1975.

[9] J.A. Sussmann and Y. Weissman, "Application of the Quantum Theory of Diffusion to H and D in Niobium" Phys. Stat. Sol. B53, 419 (1972).

[10] For instance, L. Pauling and E.B. Wilson, Introduction to Quantum Mechanics, McGraw-Hill Book Co., New York, 1935.

[11] H. Kozima, "Neutron Band in Solids", J. Phys. Soc. Japan 67, 3310 (1998).

[12] H. Kozima, K. Arai, M. Fujii, H. Kudoh, K. Yoshimoto and K. Kaki, "Nuclear Reactions in Surface Layers of Deuterium-Loaded Solids" Fusion Technol. 36, 337 (1999).

[13] H. Kozima, "Neutron Drop: Condensation of Neutrons in Metal Hydrides and Deuterides", Fusion Technol. 37, 253 (2000).

[14] G. Baym, H.A. Bethe and C.J. Pethick, "Neutron Star Matter," Nuclear Physics A175, 225 (1971).

[15] J.W. Negele and D. Vautherin, "Neutron Star Matter at Sub-nuclear Densities" Nuclear Physics A207, 298 (1973)

[16] H. Kozima, Discovery of the Cold Fusion Phenomenon - Evolution of the, Solid State-Nuclear Physics and the Energy Crisis in 21st Century, Ohtake Shuppan KK., Tokyo, Japan, 1998.

[17] H. Kozima, K. Kaki and M. Ohta, "Anomalous Phenomenon in Solids Described by the TNCF Model", Fusion Technol. 33, 52 (1998).

[18] H. Kozima, M. Ohta, M. Fujii, K. Arai and H. Kudoh, "Possible Explanation of ⁴He Production in a Pd/D₂ System by the TNCF Model" Fusion Science and Technology 40, 86 (2001).

[19] J.O'M. Bockris and Z. Minevski, "Two Zones of "Impurities" Observed after Prolonged Electrolysis of Deuterium on Palladium", Infinite Energy Nos. 5 & 6, 67 (1995-96). (NT_F)

[20] T. Mizuno, T. Akimoto, T. Ohmori and M. Enyo, "Confirmation of the Changes of Isotopic Distribution for the Elements on Palladium Cathode after Strong Electrolysis in D₂O Solution", Int. J. Soc. of Materials Engin. for Resources 6-1, 45 (1998). (NT_F)

[21] T. Ohmori, M. Enyo, T. Mizuno, Y. Nodasaka and H. Minagawa, "Transmutation in the Electrolysis of Light Water - Excess Energy and Iron Production in a Gold Electrode", Fusion Technol. 31, 210 (1997). (NT_F)

[22] G.H. Miley, G. Narne, M.J. Williams, J.A. Patterson, J. Nix, D. Cravens and H. Hora, "Quantitative Observation of Transmutation Products Occurring in Thin-Film Coated Microspheres during Electrolysis", Proc. ICCF6, p.629. (NT_F)

[23] H. Yamada, S. Narita, Y. Fujii, T. Sato, S. Sasaki and T. Omori, "Production of Ba and Several Anomalous Elements in Pd under Light Water Electrolysis" Proc. ICCF9 (to be published); Abstracts of ICCF9, p.123 (2002). (NT_Fand NT_A; Pd → Ba, Pb)

[24] S. Miguet and J. Dash, "Microanalysis of Palladium after Electrolysis in Heavy Water", Proceedings of 1st Low Energy Nuclear Reactions Conference, College Station, Texas, p. 23 (1995) (NT_A; Pd → Cd).

[25) R. Kopecek and J. Dash, "Excess Heat and Unexpected Elements from Electrolysis of Heavy Water with Titanium Cathodes", Proceedings of 2nd Low Energy Nuclear Reactions Conference, College Station, Texas, p. 46 (1996) (NT_A; Ti → Cr).

[26] J. Warner and J. Dash, "Heat Production during the Electrolysis of D₂O with Titanium Cathodes", Conference Proceedings 70 (Proceedings of 8th International Conference on Cold Fusion, Lerici, Italy), p.161 (2000) (NT_A; Ti → Cr).

[27] Y. Iwamura, M. Sakano and T. Itoh, "Elemental analysis of Pd Complexes: Effects of D₂ Gas Permeation" Jpn. J. Appl. Phys. 41, 4642 (2002) (NT_A; Cs → Pr, Sr → Mo)

[28] J. Dash, I. Savvatimova, S. Frantz, E. Weis and H. Kozima, "Effects of Glow Discharge with Hydrogen Isotope Plasmas on Radioactivity of Uranium", Proc. ICCF9 (2002) (to be published). (Decay-time shortening of ²³⁸₉₂U)

[29] I.V. Goryachev, "Abnormal Results of Experimenting with Excited Substances and Interpretation of the Discovered Effects within the Frames of the Model of Collective Interactions", Proc. ICCF9 (2002) (to be published). (Decay-time shortening of radioactive nuclides)

[30] J. Dash, I. Savvatimova, G. Goddard, S. Frantz, E. Weis and H. Kozima, "Effects of Hydrogen Isotope on Radioactivity of Uranium" Proc. ICENES2002 (to be published). (Decay-time shortening of ²³⁸₉₂U)