modeling excess heat in nih peter hagelstein 1 and irfan chaudhary 2 1 massachusetts institute of...

39
Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology, Lahore

Upload: philip-byrd

Post on 31-Dec-2015

217 views

Category:

Documents


1 download

TRANSCRIPT

Page 1: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Modeling excess heat in NiH

Peter Hagelstein1 and Irfan Chaudhary2

1Massachusetts Institute of Technology2University of Engineering and Technology,

Lahore

Page 2: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

NiH

Page 3: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Before looking at expts…

•Excess heat seen in NiH experiments

•Effect first reported in electrolysis experiments by

Mills and Kneizys (1991)

•Excess heat in gas loading experiments reported by

Piantelli et al (1994)

•NiH is not PdD

•Differences are important

Page 4: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

NiH lattice structure (fcc)

NiH

Page 5: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Phase diagram

Note: 400 Mpa = 3948 atm

Page 6: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Pressure vs loading

D/Pd

0.0 0.2 0.4 0.6 0.8 1.0

p (a

tm)

10-1

100

101

102

103

104

100 MPa

1 GPa

Page 7: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

X-ray diffraction data in electrochemical loading

No evidence for intermediate a values for loadings in the - phase region. Observed only is the change in volume occupied by -phase NiH.

Juskenas et al, Electrochimica Acta 43 1903 (1998)

Page 8: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Recall electron density near vacancy in Pd…

s [111] (Angstrom)

-3 -2 -1 0 1 2 3

(e

/Ang

stro

m3 )

0.00

0.05

0.10

0.15

0.20

0.25

O

T

V

D2

Page 9: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

[100] Displacement of D2 in PdD Supercell Monovacancy

L Dechiaro, Quantum Espresso DFT calculation

Page 10: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Electron density even higher in Ni (which is why

H doesn’t load well). But electron density reduced near a vacancy, and can form H2. Ni is

closest analog of Pd for H2/D2 formation near

vacancy. Expect issues for vacancy creation and

HD molecule formation to be similar.

Page 11: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Vacancies made more readily

Page 12: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Ni phonon modes

D. A. Dimitrov et al Phys. Rev. B 60 6204 (1999)

Page 13: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Issues with the development of optical phonon

modes in gas loading since H concentration is low

Page 14: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Donor-receiver model

HD

3He

Phonon mode

AZ*

AZ

Donor system Receiver system

Page 15: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

D/H

0.00 0.02 0.04 0.06 0.08

Pxs

(W)

0.2

0.3

0.4

0.5

0.6

0.7

200 mW

500 mWPin = 1000 mW

More Pxs with D added

M. R. Swartz, G. M. Verner, and A. H. Frank, Proc. ICCF9 p. 335 (2002).

Page 16: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

More issues

•HD/3He transition fine for donor•Reduced mass smaller than for D2/4He system

•So tunneling is orders of magnitude larger

•Deuterium natural abundance is 1/6240 of hydrogen

•E = 5.49 MeV, so need to exchange few quanta

1 2

1 1 1 2

2 3D

DD H HD H

MM M

m m

Page 17: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Take away message

•Excess heat seen in NiH

•Electrochemical systems, gas systems

•Harder to load

•Easier to make vacancies

•HD formation good

•Donor-receiver model happy

•Is some D in H•Larger interaction matrix element since Gamow factor smaller

Page 18: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Piantelli experiment

Page 19: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Piantelli experiment

S Focardi, R Habel, and F Piantelli, Il Nuovo Cimento, 107A 163 (1994)

10 cm 9 cm 0.5 cm diameter

Ni rod:

Page 20: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Calibration

Page 21: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Data showing Pxs

Input power: 140? WExcess power: 20 W

Page 22: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

T vs Pin for Pxs = 0, 20, 50W

Page 23: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

1998 Piantelli experiment

S Focardi, V Gabbani, V Montalbano, F Piantelli, S Veronesi, Il Nuovo Cimento 111 1233 (1998)

Page 24: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Calibration curves

T1

T3

T2

T4

Page 25: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Excitation of the sample

T1 = 381.7T1 = 467.4

Approx: Pin = 60 W Pxs = 20 W

Page 26: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Pressure-composition isotherms

•NiH is like PdH…1 atm = 105 Pa

Page 27: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Solubility of H at low pressure

Page 28: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Uptake of H2 after several loading cycles

Page 29: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Thinking about result

•Bulk Ni does not load much (5-20x10-5) near 1 atm

•Need O(6000 atm) to pressure-load bulk Ni with H2

•But some loading observed nonetheless in Piantelli expts

•Number of H atoms absorbed is several times O(3x1021)

•Number of Ni atoms in sample (1994) is O(4x1023)

•Loading in Cammarota replication from H2 absorption is NiH0.2

•Must be (non-bulk) special sites (defects or impurities)

•But not enough impurities in Cammarota version!

Page 30: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

X-ray diffraction data in electrochemical loading

No evidence for intermediate a values for loadings in the - phase region. Observed only is the change in volume occupied by -phase NiH.

Juskenas et al, Electrochimica Acta 43 1903 (1998)

Page 31: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Diffusion of H

•Diffusion of H in Ni is much slower

than in Pd

•Elevated temperature D in NiH is

similar to D in PdD at 300 K

D = D0 e-E/kT

D0 = 7.04x10-3 cm2/sec

E = 409 meV

D(300K) = 9.5x10-10 cm2/sec (NiH)

D(300K) = 5.5x10-7 cm2/sec (PdD)

Page 32: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

PdD diffusion model at 300 K

D(oct)/Pd(location)

0.0 0.2 0.4 0.6 0.8 1.0

DD c

m2 /s

ec

10-8

10-7

10-6

-

Page 33: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Excitation of the sample

T1 = 381.7T1 = 467.4

Approx: Pin = 60 W Pxs = 20 W

H outgassing with temperature rise

Page 34: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

-emission events

S Focardi and F Piantelli, “Produzione de energia e reazioni nucleari in sistemi NiHa 400 C” (2000)

Page 35: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

New elements in Miley expts

G. H. Miley and P. J. Shrestha, Slides presented at ICCF12 (2005).

Page 36: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

New elements in Piantelli expt

Page 37: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

•Low H loading in NiH gas systems

•Not enough H for good optical phonon mode

•NiH systems so far probably work based on acoustic

mode excitation

•Ni then participates strongly in vibrations

•Coupling of energy through nuclear excited states in

Ni

•Some long-lived ones will have fission decay

pathways

•Lattice-induced fission produces new elements

•Eats up significant amount of produced energy

Page 38: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Wish list for Santa

•Measurements of 3He in gas phase

•And correlation with energy (is Q-value 5.5 MeV?)

•Quantitative assay of Ni disintegration products

•More experiments with HD/3He branch in Pd codep

expts

•Two-laser experiments in Ni to shed light on acoustic

mode operation

•Like to have NiX alloy with lower electron density

Page 39: Modeling excess heat in NiH Peter Hagelstein 1 and Irfan Chaudhary 2 1 Massachusetts Institute of Technology 2 University of Engineering and Technology,

Take away message

•NiH excess heat discovered by Mills and Kneizys

(1991)

•Gas loaded excess heat by Piantelli group (1994)

•H doesn’t load well in Ni

•Indirect evidence for vacancy generation in

Piantelli experiment

•Excess power stimulated by H flux

•High temperature operation favors 3He diffusion

•Excess heat stays on after flux initiation

(consistent with nuclear energy going into

acoustic phonon modes

•Acoustic mode operation connected with

transmutation