avatar University of Washington Inc Services


  • Page 1

    INTRODUCTION The Nuclear Physics Laboratory of the University of Washington has for over 40 years supported a broad program of experimental physics research. The current program includes "in-house" research using the local tandem Van de Graaff and superconducting linac accelerators and non-accelerator research in solar neutrino physics at the Sudbury Neutrino Observatory in Canada and at SAGE in Russia, double beta decay and gravitation as well as user-mode research at large accelerator and reactor facilities around the world. The Mass-8 experiment has made significant progress by nearly doubling the amount of data for both 8Li and 8B beta decay. Monitoring of the gas counter performance was added to insure the quality of the data. The new high energy gamma ray beamline and detector setup has been completed and first data have been taken toward obtaining a precise spectrum of photons from the 17 MeV excited states of 8Be. As an important step in interpreting GDR emission spectra, a measurement of pre-equilibrium particle loss prior to compound nuclear formation was made for 200 MeV 18O + 100Mo. The result is, on average, a compound nucleus that is 3 mass units lighter and 20% lower excitation energy than the compound nucleus formed in complete fusion. In a study of sub-barrier fusion of Ca with Ti isotopes it has been shown that quasi-elastic scattering excitation functions can be used to obtain barrier distributions with similar structure to that obtained from fusion data. We have measured the angular distribution of fission fragments in various near- and sub-barrier heavy-ion reactions with actinide targets. We find no evidence of strong peaklike structures in fission fragment anisotropies as a function of center-of-mass energy. An analysis of C+U and O+Th fission fragment anisotropies yields no evidence of an effect due to the entrance channel mass asymmetry. A study of the three reactions C + 235,236,238U indicates a strong influence of the ground state spin of actinide targets on the anisotropy of fission fragments at sub-barrier energies. A coincidence study of the fragmentation of C60 (buckyballs) in collisions with hydrogen has demonstrated that near-binary fragmentation competes with dimer evaporation in producing fragments in the C40 to C50 mass range. Also, a new pathway for producing doubly negative carbon clusters by fragmenting alkali carbide anions has been demonstrated. 9

  • Page 2

    Our development of beryllium targets for a planned experiment has achieved a 9Be target, fabricated in the same manner as will be used eventually to make a 7Be target, with an impurity/Be ratio of 2/1. The Russian-American Gallium Experiment (SAGE) has published the results from a successful 51Cr neutrino source experiment. Analysis of the solar neutrino results is near completion. Work continues at a rapid pace toward completion of the Sudbury Neutrino Observatory. Following four years of planning and development, installation of the acrylic vessel is now over 60% complete, with final completion scheduled for this fall. The UW SNODAQ group has prototyped most of the software elements necessary for the SNO data acquisition, and has made significant progress in the coding of the online data-monitoring system. The code needed for bulk testing of the production front-end electronics is mostly in place. MiniSNO has been used to test the electronics calibration algorithms, and will be used again in the near future to test some of the optical calibration sources. The Neutral-Current Detector project has also made great strides in the last year with production of prototype counters. Preamplifiers, endcaps, cable terminations and delay lines for the proportional counters are being assembled. We have received the majority of the required CVD Ni tubing. The assembly equipment is ready and the production phase will begin soon. We have taken delivery of the remotely operated submersible vehicle to be used to deploy the neutral current detectors in the acrylic vessel. A successful run on the parity nonconserving neutron spin-rotation in 4He was carried out on the NG-6 beamline at the Cold Neutron Research Facility at NIST in Gaithersburg, MD. In December, the emiT experiment, a search for time-reversal violation in neutron beta decay of free neutrons, moved onto this beamline. The emiT system, which includes four beta and four proton detectors segments as well as beam and spin transport systems, observed first data in early January and is scheduled to continue to acquire data through mid-1997. A new algorithm for simulating the effects of multiparticle Bose-Einstein statistics and Coulomb interactions in ultrarelativistic heavy ion collisions at NA49 and STAR has been implemented in a code that is now being tested and will be used to generate a library of Bose-Einstein correlated events simulating collisions at CERN and RHIC. New limits for short range equivalence principle violations have been obtained with the rotating source torsion balance apparatus. A new continuously rotating torsion balance apparatus to test the equivalence principle for a variety of sources that are as exotic as dark matter is under construction. As always, we encourage outside applications for the use of our facilities. As a convenient reference for potential users, the table on the following page lists the vital statistics of our accelerators. For further

  • Page 3

    As always, we encourage outside applications for the use of our facilities. As a convenient reference for potential users, the table on the following page lists the vital statistics of our accelerators. For further information, please write or telephone Professor Derek W. Storm, Director, Nuclear Physics Laboratory, University of Washington, Seattle, WA 98195; (206) 543-4085 (e-mail; STORM@NPL.WASHINGTON.EDU). We close this introduction with a reminder that the articles in this report describe work in progress and are not to be regarded as publications or to be quoted without permission of the authors. In each article the names of the investigators have been listed alphabetically, with the primary author to whom inquiries should be addressed underlined. We thank Richard J. Seymour and Karin M. Hendrickson for their help in producing this report. Kurt Snover Editor Barbara Fulton Assistant Editor

  • Page 4

    TANDEM VAN DE GRAAFF ACCELERATOR A High Voltage Engineering Corporation Model FN purchased in 1966 with NSF funds; operation funded primarily by the U.S. Department of Energy. See W.G. Weitkamp and F.H. Schmidt, "The University of Washington Three Stage Van de Graaff Accelerator," Nucl. Instrum. Meth. 122, 65 (1974). Some Available Energy Analyzed Beams Ion Max. Current Max. Energy Ion Source (particle A) (MeV) 1 H or 2H 50 18 DEIS or 860 3 He or 4He 2 27 Double Charge-Exchange Source 3 He or 4He 30 7.5 Tandem Terminal Source 6 Li or 7Li 1 36 860 11 B 5 54 860 12 C or 13C 10 63 860 * 14N 1 63 DEIS or 860 16 O or 18O 10 72 DEIS or 860 F 10 72 DEIS or 860 * Ca 0.5 99 860 Ni 0.2 99 860 I 0.01 108 860 * Negative ion is the hydride, dihydride, or trihydride. Additional ion species available include the following: Mg, Al, Si, P, S, Cl, Fe, Cu, Ge, Se, Br and Ag. Less common isotopes are generated from enriched material. BOOSTER ACCELERATOR We give in the following table maximum beam energies and expected intensities for several representative ions. See "Status of and Operating Experience with the University of Washington Superconducting Booster Linac," D.W. Storm et al., Nucl. Instrum. Meth. A 287, 247 (1990). Available Energy Analyzed Beams Ion Max. Current Max. Practical (p A) Energy MeV p >1 35 d >1 37 He 0.5 65 Li 0.3 94 C 0.6 170 N 0.03 198 O 0.1 220 Si 0.1 300 35Cl 0.02 358 40Ca 0.001 310 Ni 0.001 395

  • Page 5

    1.0 FUNDAMENTAL SYMMETRIES, WEAK INTERACTIONS AND NUCLEAR ASTROPHYSICS 1.1 Beta delayed alpha spectra from 8 Li and 8 B decays and the neutrino spectrum in 8 B decay E.G. Adelberger, L.D. Carr, J.-M. Casandjian, H.E. Swanson and K.B. Swartz∗ SNO will measure the energy spectrum of 8 B solar neutrinos reaching the earth. If neutrino oscillations occur, the spectrum will be distorted from its original shape and the deviation will contain information on the neutrino mixing parameters. Bahcall et al.1 have pointed out that our ability to predict the undistorted shape of the 8 B neutrinos is limited by our knowledge of the final state continuum fed in 8 B decay. Two kinds of data are useful here, the beta spectrum and the beta-delayed alpha spectrum in 8 B decay. Bahcall et al. showed that existing delayed-alpha data were inconsistent and chose to use a single measurement of the beta spectrum in obtaining a “standard” 8 B spectrum. Because of the difficulty of making accurate beta spectrum shape measurements, we have chosen to make a careful remeasurement of the delayed alpha spectra in 8 Li and 8 B decays, paying special attention to the absolute calibration of the energy scale and to understanding the energy response of detectors. We use the “Mass 8” rotating target and moveable catcher foil apparatus to implant 8 Li (produced by 7 Li( d, p)) and 8 B (produced by 6 Li( 3 He, n )) into 10 microgram/cm 2 C catcher foils. The foils are viewed on opposite sides by a pair of Si telescopes consisting of 75 micron thick E counters followed by 500 micron thick veto detectors. The telescopes have small solid angles (∆Ω / 4π = 2.2 × 10 −3 ) to minimize summing of alphas with the associated beta particle. We can, without breaking vacuum, insert thin 148 Gd and 241Am sources in front of the detectors. In addition we can place thin Al sheets in front of the telescopes to eliminate the alphas and see only betas. Finally we can measure the thickness of the catcher foil by measuring the energy loss of 244Cm alphas passing through the foils - all without breaking vacuum. The detector telescopes are cooled to 0°C, and the electronics (except for preamps) is mounted in a special temperature-controlled rack. Analysis of 8 Li data is in progress. A 8 B data taking run is scheduled. ∗ Yale University, New Haven, CT. 1 J.N. Bahcall et al., Phys. Rev. C 54, 411 (1996). 1

  • Page 6

    1.2 Progress in the Mass-8 β-decay experiment M. Beck, L.D. Carr, H. Lawler, D.W. Storm, H.E. Swanson and J.P.S. van Schagen The Mass-8 isospin triplet presents an ideal system to test CVC and look for the existence of second class currents.2 The weak-decay from 8 Li and 8 B nuclei is being extensively studied in our lab using the Mass-8 apparatus.3 Several hardware improvements were implemented. Air-cooled Ta-slits were installed to stop the beam when it is swept away during the catcher foil arm rotation. Furthermore, the collimator assembly near the target was improved to further reduce the activity landing on the catcher frame. A method has been developed in which the wires of the α -counters are cleaned prior to an experiment by submerging the wire frames in trichloroethane and using ultrasound. The wires are then rinsed with acetone. The software used to obtain an in situ calibration of the β -counters4 was extended to include the possibility to use radioactive β-sources. The program used to fit the β − α angular correlation has also been extended. For each energy E β this angular correlation is given by: W(θ β −α ) = a 0 [1 + a1cos(θ β −α ) + a 2 cos2 (θ β −α )] (1) or an equivalent expression in terms of Legendre polynomials. The new software is able to determine the best values for the coefficients a 0 , a1 and a 2 by a χ 2-minimization procedure either by fitting the angular correlation for each β-energy bin independently or from a global fit in which either or both of a0 and a2 can be constrained to a linear polynomial in E β . We have performed a three week experiment in which the weak decay from 8 Li has been studied. In total 1.5 × 10 +8 events were collected. The preliminary results of the 8 Li experiment are summarized in Fig. 1.2-1 below. In this fit, the angular correlation was fitted using Eq. 1 and the coefficient a2 was constrained to a 2 / E β = constant. We find a 2u / E β = (3.20 ± 0.15) × 10 −3 MeV −1 and a 2d / E β = (2.93 ± 0.15) × 10 −3 MeV −1 , where u,d corresponds to the upstream and downstream α-counter, respectively, with a typical χ 2 / ν of 0.98. These results are in good agreement with a previous measurement5 where a value a 2 / E β = (3.34 ± 0.26) × 10 −3 MeV −1 was obtained. Currently, additional data is being collected for 8 B. Furthermore, the γ-ray decay from the isobaric analog state in 8 Be will be studied in the near future using the extended NaI set-up (see Sections 6.6 and 6.7). a0 a1 a2 1.1 0.0 0.06 0.05 a2 fit (see text) a2 fit (see text) 0.04 1.0 -0.1 0.03 0.02 theory 0.01 0.9 -0.2 0.00 0 5 ✁ 10 15 0 5 ✁ 10 150 5✁ 10 15 Gated Eβ [MeV] Gated Eβ [MeV] ✂ Gated Eβ [MeV] Fig. 1.2-1. Result for the coefficients a 0 , a1 and a 2 for the decay 8 Li→8 Be + e − + ν using Eq. 1. The fit for a 2 has been constrained to a linear polynomial. The solid (open) data involve the downstream (upstream) α -counter. 2 Nuclear Physics Laboratory Annual Report, University of Washington (1996) p. 1. 3 Nuclear Physics Laboratory Annual Report, University of Washington (1995) p. 27. 4 Nuclear Physics Laboratory Annual Report, University of Washington (1996) p. 2. 5 K.B. Swartz, Measurement of the ✄Angular Correlation of Decay in 8 Li , Ph.D. thesis, University of Washington, ☎ ✆ 1996. 2

  • Page 7

    1.3 ν correlation in 32 Ar decay The e-ν E.G. Adelberger, M. Beck, A. Garcia* and H.E. Swanson We are continuing to analyze the results of CERN/ISOLDE experiment IS332. The focus is on the lineshape of the delayed proton group corresponding to the superallowed 0+ → 0+ decay of 32 Ar . Two completely independent Monte Carlo analysis programs have been written - one is based on the CERN PAW package, the other uses Microsoft Fortran on the PC. In order to improve our analysis and minimize errors in the e-ν correlation coefficient a, due to uncertainties in the detector response function, we have made a study of the response function of Si detectors to ionizing fast particles (see Section 1.4). Once this response function is understood we will return to our extraction of the e-ν correlation coefficient. We have an approved proposal for a new measurement of the e-ν correlation in 32 Ar decay (A. Garcia, spokesman). This new measurement will use a superconducting solenoid to “curl up” the beta particles so that they will not reach the proton detector. This will allow us to have larger detector solid angles improving the statistical precision of the result, and reduce systematic uncertainties that arise from the distortion of the superallowed proton peak by summing the betas. Design work on the apparatus is being carried out in collaboration with our Notre Dame colleagues. * Notre Dame University, Notre Dame, IN. 3

  • Page 8

    1.4 Determining the response function of a silicon surface barrier detector E.G. Adelberger, M. Beck, H. Bichsel and H.E. Swanson In order to understand the response function of the proton detector in the 32Ar delayed-proton measurement of the e-ν correlation (Section 1.3) we have calculated the α-particle response function of a Si-detector and compared it to data taken with a mono-energetic α-source. The model includes electronic straggling7,8 in the α-source (thickness Ds) and in the detector deadlayer (thickness Dd), nuclear straggling9 in the active volume of the detector, the source diameter d, and the solid angle between the source and the detector (square of side length l). The ionization efficiency of recoiling nuclei is calculated according to Lindhard et al.3 with the parameterization of Smith and Lewin.10 Ionization statistics are taken into account with W = 3.6 eV/electron-hole-pair and a Fano factor of 0.1. The electronic resolution was measured during the α-source measurement with a precision pulser. Nuclear straggling in the source and in the deadlayer are negligibly small. The electronic straggling functions for α-particles traversing the Si deadlayer and Gd2O3 source were calculated8 from photoabsorption and electron energy loss spectra.11,12 The nuclear straggling function was computed with a Monte Carlo convolution of the nuclear single-collision cross section. After the two straggling functions were combined the result was folded with Gaussian electronic (FWHM = 5.2 keV) and statistical (FWHM = 2.5 keV) resolutions. The calculated lineshape is shown in Fig. 1.4-2. Its width is dominated by the straggling in the source and the deadlayer whereas the low-energy tail is caused by nuclear straggling in the active volume of the detector (Fig. 1.4-1). The experimental lineshape of 3183 keV α-particles in a Si-detector (l=0.82 cm, Dd=58 nm) was measured with a d=0.5 cm 148Gd2O3 source. Ds=29 nm and Dd=58 nm were determined by analyzing spectra taken with the source translated and/or rotated with respect to the detector. Fig. 1.4-2 shows the comparison of the measured lineshape with the calculated response function. We find good agreement between the calculated response function and the measured lineshape. We will now use the model to determine the detector response to protons. Electronic straggling (dotted) Nuclear straggling (solid) Energy axis matched at maximum Calculated area normalized to data Fig. 1.4-1. Contributions to the response function. Fig. 1.4-2. Calculated response function and comparison with the measured lineshape. 7 H.Bichsel and R.P. Saxon, Phys. Rev. A 11, 4, 1286 (1975). 8 H.Bichsel, Rev. Mod. Phys. 60, 3 663 (1988). 9 J.Lindhard et al., K. Dan. Vidensk. Selsk., Mat-Fys. Medd. 33,10 (1963). 10 P.F. Smith and J.D. Lewin, Phys. Rep. 187 203 (1990). 11 P. Schattschneider, F. Foedermayer and D. S. Su, Scanning Microscopy, 2, (1988). 12 E. Hencke, Atomic Data and Nuclear Data Tables 54, (1993). 4

  • Page 9

    1.5 Spin-decomposition of the GT strength in 37 Ca decay E.G. Adelberger, S.E. Darden,† A. Garcia,† W. Haeberli,‡ N.I. Kaloskamis,† E. Miller,† P.A. Quin,‡ B.P. Schwartz‡ and E. Yacoub‡ We13 had previously noted significant discrepancies between the GT strength observed in 37 Ca decay (see references14,15 for more recent experimental results) and shell model calculation using the Universal 2s − 1d Hamiltonian. In an attempt to better understand the discrepancy and provide data that could be used to improve the shell-model Hamiltonian, we made an extensive study of the spins of the final states fed in 37 Ca decay. This work was done by studying the daughter 37 K levels as 36 Ar ( p, p) resonances using the polarized beam at the University of ✝ Wisconsin tandem accelerator. Spins, parities, and widths of 17 37 K levels were determined, and the isospin mixing of the lowest T = 3/2 level was inferred. We found that neither the Universal or the earlier CW shell-model interaction reproduced the results. However, the Universal interaction predictions would be considerably improved if the energy of the J = 5 / 2 + ; T = 1/2 component of the GT giant resonance were lowered by about 2.5 MeV. This result should therefore point the way toward improvements in the shell model residual interaction. As this work has recently been published16 it will not be described in detail here. † Notre Dame University, Notre Dame, IN. ‡ University of Wisconsin, Madison. WI. 13 E.G. Adelberger et al., Phys. Rev. Lett. 67, 3658 (1991). 14 A. Garcia et al., Phys. Rev. C 51, R439 (1995). 15 W. Trinder et al., Phys. Lett. B 349, 267 (1995). 16 N.Kaloskamis et al., Phys. Rev. C 55, 630 (1997). 5

  • Page 10

    1.6 Beta decay of 40 Ti and the efficiency of the ICARUS 40 Ar neutrino detector E.G. Adelberger, R. Anne,† M. Bhattacharya,‡ C. Donzaud,§ A. García,‡ S. Grévy,§ D. Guillemaud- Mueller,§ N.I. Kaloskamis,‡ M. Lewitowicz,† A.C. Mueller,§ F. Pougheon,§ M.G. Saint-Laurent,† O. Sorlin,§ H.E. Swanson and W. Trinder† The proposed large-volume liquid-argon detector ICARUS17 has several advantages over existing solar neutrino detectors because it will detect neutral and charged-current neutrino interactions in a very symmetrical way. ✞ ✟ Neutral-current interactions will be characterized by single-track e ν , ν e events, while charged-current 40 Ar(ν e − ) 40 K ∗ interactions will produce multiple tracks because the J π = 1+ states fed in allowed neutrino e, capture emit several γ rays as they decay to the 40 K J π = 4 − ground state. Therefore the multiplicity and angular distribution of the event will signal its neutral or charged current nature. The neutral-current efficiency of ICARUS can be accurately calculated using electro-weak theory. However, the charged-current efficiency depends on the matrix elements for neutrino-capture transitions on 40 Ar to excited states of 40 K . A recent shell-model calculation18 predicts a capture rate of 6.7 ± 2.5 SNU [1 SNU = 10 −36 captures per target atom per second], where 2.2 SNU arises from the model-independent Fermi cross section and 4.5 SNU is expected from the model-dependent Gamow- Teller transition strengths, B(GT). An empirical calibration of the 40 Ar(ν e, e − ) transition strengths, is therefore essential. We made such a calibration by studying the β + decays of 40 Ti and used isospin symmetry to relate the 40 Ar → 40 K transitions to the strengths of the mirror 40 Ti → 40 Sc transitions. Very little was known about 40 Ti as it lies far from the valley of stability. Detraz et al.19 observed β delayed protons from 40 Ti decay, but their statistics were rather poor. They obtained a half-life of 56 +−18 12 ms and observed four delayed proton branches whose branching ratios summed to 43 ±6 %, implying that about 60% of the decays were followed by γ rays, or by protons that were below their detection threshold. We produced 40 Ti at GANIL by projectile fragmentation of a 82.6 MeV/u 50 Cr beam on a 272.4 mg/cm2 nickel target. The 50 Cr beam was produced in an ECR ion source using isotopically enriched feed material. The momentum analysis of the reaction products was performed using the ALPHA20 spectrometer with a momentum acceptance of 0.6% around Bρ = 715.6 MeV/ec. Fragments of interest were then selected using the LISE3 spectrometer. ✠ ✡ ✠ ✡ Our charged particle detection system consisted of 5 surface-barrier Si detectors, D1 150 µ m , D 2 150 µ m , ✠ ✡ ✠ ✡ ✠ ✡ D 4 500 µ m , D 5 500 µ m and D 6 500 µ m . Each of the first two detectors (D1 and D 2 ) provided energy loss and time-of-flight information which gave two independent identifications of the incoming fragments. The last three detectors (D 4 , D 5 and D 6 ) formed a telescope used for implanting the fragments and detecting their subsequent β decays. A remotely-controlled carbon energy degrader of variable thickness was adjusted so that the 40 Ti ions were implanted in the center of D 5 . Detectors D 4 and D 6 registered the β particles emitted by ions implanted in D 5 . Five “70%-efficient” HPGe detectors, mounted close to D 5 , detected β -delayed γ rays. A preliminary analysis of the data, which contains about 3 × 10 4 40 Ti decays, yields a 40 Ti half-life of 53.7 ± 1.4 ms and we observe proton decays whose branches sum up to close to 100%. Analysis of the data is continuing. † GANIL, BP 50-27, 14021 Caen Cedex, France. ‡ University of Notre Dame, Notre Dame, IN. § Institut de Physique Nucléaire, 91406 Orsay Cedex, France. 17 ICARUS Collaboration, Proposal, ICARUS II, A second-generation proton decay experiment and neutrino observatory at Gran Sasso. 18 W.E. Ormand, P.M. Pizzochero, P.F. Bortignon and R.A. Broglia, Phys. Lett. B 345, 343 (1995). 19 C. Détraz et al., Nucl. Phys. A 519, 529 (1990). 20 M. Lewitowicz et al., Phys. Lett. B 332, 20 (1994). 6

  • Page 11

    1.7 Large helium cell for measurement of the spectrum of photons in 8 Be K.A. Snover, D.W. Storm and J.P.S. van Schagen Part of the Mass-8 test of CVC and search for second class currents involves measuring the distribution of isovector M1 transition strength to the broad 2+ final state at 3 MeV in 8 Be . The initial states are formed by 4 He - 4 He fusion, and the width of the final state is reflected in the photon spectrum. In order to determine the M1 matrix element accurately, the energy distribution of the emitted photons needs to be known. The spectrum can be measured with a detector at 90°, using a long gas cell with the windows shielded from the detector. In earlier measurements it was determined that a significant background was produced by scattered helium which subsequently interacted with the walls of the cell. Consequently we built a large gas cell, long enough that the windows can be well shielded and large enough in diameter that He must be scattered through an angle greater than 45° to hit the part of the walls viewed by the detector. Preliminary design details were presented in last year’s Annual Report.21 At 45° the He has only one-half of the 34-MeV beam energy. The cell walls are lined with Ta, so that subsequent nuclear interactions are suppressed by the Coulomb barrier. We measured the rate for photon production on thick Ta by He of energies from 18 to 24 MeV and fit the energy dependence to an exponential to extrapolate to lower energies. By integrating along the beam, over the He-He scattering cross section, and over the photon production rate, we found that we should expect 4 × 10 −4 photons from the wall for every photon produced in the resonant He-He interaction. Most of these wall- produced photons have energy below 7 MeV, so the problem of background from the wall should be solved. The windows are shielded with 6 inches of lead, which should attenuate the low energy photons by a factor of 1000. The background from the windows can be measured by running with H 2 in the target. We have performed a preliminary measurement of the photon spectrum from the resonant states in 8 Be using this cell and the Seattle NaI detector. The spectrum, which peaks at about 14 MeV, appears to be clean above about 10 MeV. Work is in progress to analyze these data. 21 Nuclear Physics Laboratory Annual Report, University of Washington (1996) p. 57. 7

  • Page 12

    1.8 Target and apparatus development for a planned 7 Be(p, γ ) experiment E.G. Adelberger, N. Bateman,* L. Buchmann,* J.-M. Casandjian, R. Hoffenberg, K.A. Snover, T.D. Steiger, J.F. Wilkerson and A. Zyuzin* We are developing the target fabrication and the experimental apparatus for a measurement of the astrophysically interesting 7 Be( p, γ )8 B reaction in the energy range Ep = 0.3 to 1.0 MeV. The experiment is a collaboration in which the radioactive 7 Be target will be fabricated at TRIUMF and the cross section measurements carried out here at the NPL. Much of the effort of the past year has been focused on target preparation and testing. The previous 22 scheme of Be and Al layers on a Cu backing has been abandoned in favor of a simpler design of a single Be layer on a conductive backing. A series of targets of 5-10 mm diameter of 9 Be have been produced at TRIUMF by a combination of chemistry and vacuum distillation using a procedure similar to that intended for the fabrication of the 7 Be targets. Testing of the targets has been carried out by heavy ion sputtering diagnostic tests in Vancouver, British Columbia, by ( p, γ ) resonance studies and by a search for characteristic γ lines from 2.4 MeV proton bombardment. The 9 Be( p, γ ) reaction at Ep ≅ 1 MeV has been used to characterize the 9 Be thickness of the targets, as well as the purity. The present results of this iterative procedure are targets of the desired thickness with an impurity/Be ratio of about 2/1. Efforts at further improvement are continuing. A beamline (L45) has been identified as suitable for dedication to our experiment. An additional magnetic quadrupole lens has been obtained23 and tested together with some rudimentary optics on L45. Design of the target box and components are now underway. The water-cooled target and rotating arm have been designed, as well as the Faraday cup and the LN2 cold trap. Our goal is to have a working apparatus and preliminary measurements on a low activity 7 Be target by the end of this year. * TRIUMF, Vancouver, BC, Canada. 22 Nuclear Physics Laboratory Annual Report, University of Washington (1996) p 7. 23 Long-term loan from the State University of New York, Stony Brook. 8

  • Page 13

    1.9 Measurement of the PNC spin-rotation of cold neutrons in a liquid helium target E.G. Adelberger, J.H. Gundlach, B.R. Heckel, D.M. Markoff, S.D. Penn and H.E. Swanson We have completed the first measurement of the parity nonconserving spin-rotation of cold neutrons transmitted through a liquid helium target. This parity-violating observable is sensitive to the pion exchange term in the nucleon-nucleon weak interaction. A determination of the pion coupling constant, fπ , will provide information on the neutral current contribution to the N-N weak interaction. Discrepancies remain between the experimental constraints on fπ , set by γ -ray circular polarization measurements in 18 F , and theoretical calculations.24 After tests of the entire apparatus were completed at NPL, the experiment was shipped to the National Institute of Standards and Technology (NIST) in Gaithersburg, Maryland. The apparatus was then reassembled on the Fundamental Physics beam line (NG-6) at the Cold Neutron Research Facility (CNRF) which provided a polarized cold neutron (wavelengths greater than 3.9Å) flux of 2 × 108 neutrons/cm 2 /sec. The spin-rotation apparatus remained on the NG-6 line for three 6-week reactor cycles. The apparatus consists of a crossed polarizer-analyzer system with magnetic field coils to preserve the neutron spin in transition to the low field region, and two target chambers separated along the beam axis by a 180 degree spin-rotating coil, the π -coil. The beam profile was measured at the entrance and exit of the magnetic field transition coils and at the entrance to the detector located immediately downstream of the analyzer. The neutron transmission from the polarizer to the detector was determined to be approximately 3%. The polarization product for this system was measured for various states of target configuration and π -coil settings. During the second and third cycles, we optimized the running of the polarimeter and collected data. A description of the data acquisition system is presented in Section 1.10. To minimize false signals from magnetic field effects, we used two layers of passive magnetic shielding and two active magnetic field trim coils (one producing a uniform field, the other a gradient field) to reduce the residual magnetic fields inside the target region to less than 100 µ Gauss. We have characterized the magnetic field drifts including changes in ambient fields from equipment, most notably the overhead crane. We were able to monitor and minimize magnetic field rotations with the use of our segmented 3 He ionization chamber, described elsewhere in this report (see Section 1.11). The active trim coils were set so that the background angular spin-rotations, on the order of a few milliradians, were identical within 1 mrad in all velocity separated regions of the detector. We use the segmented detector data to separate the geometry and velocity independent parity non-conserving rotation from background rotations. The cryogenic system which supports the liquid helium targets suffered from a number of complications. This included increasingly problematic vacuum cold leaks in the helium stack of the Oxford Instruments cryostat, cold shock created vacuum leaks in the target feedthrough tubes, and damage to the target helium transfer pump from ice entering the liquid helium reservoir. The loss of data acquisition time due to these problems was about 7 weeks. Preliminary analysis of the data indicates a neutron spin-rotation in helium of ( −6.9 ± 7.3) × 10 −7 radians/cm, which is within a factor of 1.2 of statistical noise. Section 1.12 discusses the analysis program and the techniques used to extract the rotation signal from background reactor fluctuations and magnetic field drifts. This result, φ PNC ( n + α ) , in conjunction with the measured longitudinal analyzing power of polarized protons in helium, ☛ A L ( p + α ) ,25 implies a range for the pion coupling constant of 6 ≤ fπ ≤ 55, in units of 3.8 × 10 −8 . (For comparison, ☛ the theoretical range is 0 ≤ fπ ≤ 30 .26) We have successfully shown that the spin-rotation measurement in the n-α system can be used to extract information on fπ . 24 W. Haberli and B.R. Holstein, “Parity Violation in the Nucleon-Nucleon System,” in Symmetries and Fundamental Interactions in Nuclei, W.C. Haxton and E.M. Henley, eds. World Scientific, 1995 p. 58. 25 J. Lang et.al., Phys. Rev. C 34, 1545 (1986). 26 B. Desplanques, J.F. Donoghue, and B.R. Holstein, Ann. Phys. (NY) 172, 100 (1986). 9

  • Page 14

    1.10 NSAC - data acquisition for the neutron spin-rotation experiment H. E. Swanson The data acquisition system and apparatus for the neutron spin-rotation experiment were described in last years annual report.27,28 The program NSAC functions to sequence external fields in the apparatus and to record signals from an array of neutron detectors. It was written in Microsoft Fortran 6.1 with simple real time extensions written in Assembler. NSAC was used throughout 1996 at NIST to obtain the data discussed in the accompanying article in this report (see Section 1.9). The program was upgraded to include the following capabilities: Additional data-taking modes were added to measure the Polarization of the neutron beam and to scan the magnitude of external fields, to tune the polarimeter, and to investigate systematics. The neutron beam polarizations measured for each detector element are written to a file and are used to calculate the spin-rotation angle. Online analysis routines were added to compute and display various combinations of asymmetries in the data. These greatly aided online monitoring of the experiment. Multitasking capability was added to the target cycling program so that analysis could be performed at the same time as the somewhat lengthily process of filling and draining the targets. Non-data-taking modes were added to cycle the mechanical parts of the fill system during cooldown to prevent them from freezing. 27 Nuclear Physics Laboratory Annual Report, University of Washington (1996) p. 8. 28 Nuclear Physics Laboratory Annual Report, University of Washington (1996) p. 55. 10

  • Page 15

    1.11 Design and testing of a 3He ionization chamber for use in measuring the velocity profile of a cold neutron beam E.G. Adelberger, B.R. Heckel, D.M. Markoff, S.D. Penn, and H.E. Swanson We have designed, constructed, and tested a 3He ionization chamber (IC) for use in the Neutron Spin Rotation experiment which ran from July-November 1996 (see Section 1.9). Unlike standard ionization chambers, our detector operates with a low density of 3He and has multiple charge collection regions along the beam axis. This configuration allows us to calculate the beam’s velocity spectrum from the absorption profile in the IC. The detector was originally built and tested in the summer and autumn of 1995. During the past year we have overhauled the gas handling and electrical systems. This final design was tested in July just prior to the experimental run. At standard operating voltage (1200 V) the detector channels registered a dark current of 0.2 pA compared to a typical signal current of 30 nA. Calculated currents differed from measured currents by only 2% where the majority of this error arises from an inability to accurately characterize the neutron beam. Our experiment placed four important restrictions on the detector design. The incident cold neutron beam ( λ N ≈ 5 Å ) was split into two parallel beams separated by 1 cm. We required that the two beams be measured separately with no left-right crosstalk. Since our liquid He targets were frequently cycled, we needed to separately measure the upper and lower beam halves to ensure that the targets had properly filled/drained. Our detector noise had to be small compared to the typical signals of 30 nA. Finally we had to establish collection regions along the beam axis which had similar statistics and maximal velocity separation. To reduce our error in measuring ϕ PNC it was essential to minimize longitudinal magnetic fields which also rotate the neutron spin. Using the fact that α M , the rotation due to the magnetic field, depends on velocity, we required enough velocity separation in the IC to detect the spin rotation arising from a 40 µ Gauss B field, the limit of our static µ−metal shields. The detector consisted of an alternating series of circular HV and grounded sensor plates. The plate separation was 3 cm, and the operating voltage was 1.2 kV. To reduce fringe fields an annular ring at half voltage was positioned midway between the HV and sensor plates. The first plate differed from the rest; it had a plate spacing of 1.5 cm, operated at 600 V and had no half voltage ring. Each sensor plate (5" OD) had a 1" wide outer annulus (grounded) and inner sensing plate (0.8 mm thick, 3" OD) which was divided into four quadrants. The quadrant plates were mounted on a quartz annulus with a separation of 1 mm. The plates defined four sensing regions with depths of 1.5, 3, 6, and 9 cm respectively. The plate stack was surrounded by an Al cylinder and filled with a gas mixture of 3 atm Ar and 0.5 atm 3He. The 3He density was chosen so that the percentage of the cold neutron beam absorbed in each region would be 30, 30, 30, and 10% respectively. The Ar density was chosen so that the proton range would be small compared to 1 cm, the separation of the left and right beams. The detector was tested at the NIST Cold Neutron Research Facility with a full spectrum and a monochromatic beam. The full spectrum beam has a Maxwell distribution peaked at 5 Å with a high energy cut-off at 4.4 Å. Using this distribution, the predicted plate currents precisely matched our measurements, and there was no measurable left-right crosstalk. On the other hand, the monochromatic beam measurement yielded an unexpected result. The monochromatic cold neutron beam is created on the neutron interferometer beamline using Bragg reflection. Our measurements indicated that the assumed monochromatic beam was actually polluted by a second Bragg reflection which accounted for one third of the flux. Our findings were later supported by independent measurements. Currently to measure the beam velocity profile at the CNRF one must scan through the neutron wavelengths using either a Bragg reflection system or a beam chopper for time-of-flight separation. Both processes are time consuming and are susceptible to reactor fluctuations. With our low density 3He detector we have demonstrated that we could essentially take a snapshot of the beam velocity profile, a process which is both easy and independent of beam fluctuations. However since our detector has only four velocity sectors the profile is rather crude. We are currently collaborating with the NIST Neutron Interferometer Group to build a second detector with 10-15 absorption regions to test the practicality of measuring the CNRF beam profile in this manner. 11

  • Page 16

    1.12 Development of the data analysis code, DANSR, for the neutron spin-rotation experiments E.G. Adelberger, B.R. Heckel, D.M. Markoff, S.D. Penn and H.E. Swanson We have developed the code DANSR (Data Analysis for Neutron Spin Rotation) to analyze the data from the Neutron Spin Rotation experiment which was conducted from July - November 1996 (see Section 1.9). The code’s basic functions are to read in the raw data files generated by the acquisition code NSAC (see Section 1.10), and to calculate the PNC spin rotation angle, ϕ PNC , as well as several systematic values such as the magnetic rotation angle, α M and the raw asymmetries. Since ϕ PNC is known to be small (10 −7 rads) compared to the total spin rotation (θ ≈ 10 −3 rads), the main thrust of the code is in developing methods to reduce systematic noise such as reactor fluctuations, magnetic field drifts, magnetic field spikes, and variations in the target conditions. I will explain the current strategies used in the code with the caveat that much of the analysis remains to be done. During data acquisition NSAC repeatedly cycles the experiment through the states for the end coil field (±), data collection subcycles, pi coil field ( π = − 1, 0, 1 ), and target position ( T = 1, 2 ), in order of most to least frequently cycled state. The two end coil states alternately measure the neutron spin projection along ± y• , where the beam direction is along •z and the initial polarization along x• . By switching between target positions in front of and behind the π coil, we alternate the sign (±) of our ϕ PNC signal. For each system state NSAC writes out the integrated charge from the 16 detector channels (see section 1.11). We can denote Si, ±,π ,T as the channel readout for a given system state where we have integrated over subcycles. Let N i,± ,π, T be the total number of neutrons of a given polarization which would be absorbed in an ideal, lossless system on channel i, in state ±, π, T . Then, Si, ±,π ,T = N i, ±,π ,T (1 ± Ρ i,π ,T sin ( θ ) ) 2 where Pi,π ,T is the polarization product. If N +,… = N −,… then the raw asymmetry Ai,π , T = (S + − S− ) (S+ + S− ) = P sin (θ ) , from which we can extract ϕ PNC using quadratic smoothing as described below. Unfortunately because of reactor fluctuations, N +,… ≠ N −,…on the 0.1 % level. To eliminate these fluctuations we assume that for a Left-Right plate pair N L ,+ ,π, T N L , −,π ,T = N R,+ ,π, T N R ,−,π ,T , and we calculate: S R, +,π ,T S L ,−,π , T 2 ℜ π ,T = S R, −,π ,T S L ,+,π , T [ = 1 + 2 PR ,π , T sin (θ R ) − PL ,π ,T sin (θ L ) + O sin (θ ) ] ( ) If we further assume that Pi,π ,T does not depend on target state, then we can define ☞✌✍ ✎✏ ✑ ✒ ✓ T +1 ℜπ , T − 1 2 Pα′ ,π α M + −1 Pϕ′ T ,π ϕ PNC ϕ π′ ,T = = where Pϕ′ ,π = Pϕ′ 1,π + Pϕ′ 2 ,π . Pϕ′ ,π Pϕ′ ,π The combined polarization product and rotation angle due to the magnetic field, Pα′ ,πα M , is a smooth background which is independent of target state. ϕ PNC is then an oscillating signal on a background which slowly varies as the magnetic field drifts. Therefore we can extract ϕ PNC via quadratic smoothing. For a series of target states, labeled j: ϕ PNC,π , j = −1 j+1 ϕ π′ , j − 3ϕ π′ , j+1 + 3ϕ π′ , j+2 − ϕ π′ , j+3 8 ✒ ✓ ✔ ✕ A first pass through our data has yielded a result of ϕ PNC = −6.9 ± 7.3 × 10 −7 rad. This number is still ✖ ✗ quite rough since we have not eliminated corrupted data nor analyzed our systematics. We are also in the process of understanding the dependence of ϕ PNC on several systematics such as α M , A, and the count rate asymmetry ℵi ,π = (Si,π ,1 − Si,π, 2 ) ( Si,π ,1 + Si ,π ,2 ) . 12

  • Page 17

    1.13 Time reversal in neutron beta decay - The emiT experiment S.R. Elliott, H.P. Mumm, A.W. Myers, R.G.H. Robertson, T.D. Steiger, T.D. Van Wechel, D.I. Will and J.F. Wilkerson The emiT experiment is a search for a violation of time-reversal (T) invariance in the beta decay of free neutrons. The experiment utilizes a beam of cold (<10 meV), polarized neutrons from the Cold Neutron Research Facility at the National Institute of Standards and Technology (NIST) in Gaithersburg, MD. A sizable team of scientists has been assembled to perform this experiment from Los Alamos National Laboratory, NIST, the University of California at Berkeley/Lawrence Berkeley National Laboratory, the University of Michigan, the University of Notre Dame and the University of Washington's Nuclear Physics Laboratory. emiT will probe the T-odd triple correlation (between the neutron spin and the momenta of the neutrino and electron decay products) in the neutron beta-decay distribution.29 The coefficient of this correlation, D, will be measured by detecting decay electrons in coincidence with recoil protons while controlling the neutron polarization. Technological advances in neutron polarization and proton detection--in addition to an improved detector geometry-- should allow emiT to attain a sensitivity to D of 3 × 10 −4 . This level of sensitivity represents a factor of five improvement over previous neutron T tests, and may permit restrictions to be placed on several extensions to the Standard Model that allow values of D near 10 −3 . emiT is the first neutron T test to make use of a “supermirror” neutron polarizer. Thus, emiT achieves a polarization of 95±2%, as opposed to the 65-85% polarizations typical in previous experiments. The emiT detector consists of four plastic scintillator paddles for electron detection and four arrays of large-area PIN diodes to detect the protons. The PIN diodes have been extensively tested and shown to be efficient and economical detectors of low- energy (~ 30 keV) protons.30 The eight detector segments are arranged in an alternating octagonal array about the neutron beam so that each segment of one type lies at an angle of 135° relative to two segments of the other type. This geometry takes advantage of the fact that the electron-proton angular distribution is strongly peaked due to the disparate masses of the decay products. When compared to the 90° geometry used in previous experiments, this octagonal geometry results in an increase in signal rate which is the equivalent of a factor of seven increase in neutron beam flux. The emiT experiment was installed on the NG-6 beamline at NIST in November of 1996. Roughly two months were spent carefully characterizing the neutron beam and performing the initial shakedown of the detector. A run commenced in early February 1997 which entailed collecting preliminary data under full running conditions. The detector is operating roughly as expected and electron-proton coincidences are being detected at rates up to 3 Hz. The reactor was down for scheduled maintenance during March, 1997. Hence the bulk of emiT's data collection will occur during April-July, 1997. 29 J.D. Jackson, S.B. Treiman and H.W. Wyld, Jr., Phys Rev. 106, 517 (1957). 30 E.G. Wasserman, Time Reversal Invariance in Polarized Neutron Decay, Ph.D. thesis, Harvard University, 1994. 13

  • Page 18

    1.14 Detection of radio-frequency pulses associated with cosmic ray air showers J.L. Rosner* and J.F. Wilkerson As a result of work in the 1960's and 1970's31,32 some of which has continued beyond then, it is recognized that air showers of energy 1017 eV are accompanied by radio-frequency (RF) pulses, whose polarization and frequency spectrum suggest that they are due mainly to the separation of positive and negative charges of the shower in the Earth's magnetic field.33 The most convincing data have been accumulated in the 30-100 MHz frequency range. However, opinions have differed regarding the strength of the pulses, and atmospheric and ionospheric effects have led to irreproducibility of results. In particular, there may also be pulses associated with cosmic-ray-induced atmospheric discharges.34 A study is being undertaken of the feasibility of equipping the Auger35 array with the ability to detect such pulses. RF pulses may be able to provide auxiliary information about primary composition and shower height. It is possible that the higher energy of the showers to which the array would be sensitive would change the parameters of detection. Before a design for large-scale RF pulse detection can be produced, it has been necessary to retrace some of the steps of the past 30 years by demonstrating the existence of the pulses for 1017 eV showers, and by controlling or monitoring some of the factors which led to their irreproducibility in the past. A prototype detector has been set up at the CASA/MIA36 site in Dugway, Utah. A trigger based on the coincidence of several muon “patches” was set to select ‘large’ showers with a rate of 30-50 per hour . A log-periodic antenna sensitive to 26-170 MHz was mounted about 30 meters to the east of the CASA boundary, at an elevation of about 10 meters to place it just above the lightning protection grid. A digital storage scope was used to register filtered and preamplified RF data on a rolling basis. These data were then captured and stored upon receipt of a large-event trigger. A total of more than 5000 triggers, obtained under various conditions of filtering, preamplification, and noise reduction, have been acquired during runs in September 1996, December 1996, and February 1997. Preliminary results from initial analysis indicates that about 1 in 10 of these triggers may contain information on RF pulses from air showers. Shower sizes, core locations, and incident shower angles are being correlated with time, duration, and intensity of radio bursts. Detailed analysis of these data are in progress. * Enrico Fermi Institute and Department of Physics, University of Chicago, Chicago, IL 60637. 31 P.R. Barker, W.E. Hazen and A.Z. Hendel, Phys. Rev. Lett. 18, 51 (1967); W.E. Hazen ibid., 22, 35 (1969); 24, 476 (1970). 32 H.R. Allan, in Progress in Elementary Particles and Cosmic Ray Physics, v. 10, edited by J.G.Wilson and S.G. Wouthuysen North Holland, Amsterdam 1971) p. 171, and references therein, V.B. Atrashkevich et al., Yad. Fiz. 28, 366 (1978). 33 F.D. Kahn and I. Lerche, Proc. Roy. Soc. A 289, 206 (1966); S.A. Colgate, J. Geophys. Res. 72, 4869 (1972). 34 R.R. Wilson, Phys. Rev. 108, 155 (1967); R.A. Roussel Dupré, A.V. Gurevitch, T. Tunnell, and G.M. Milikh, Los Alamos National Laboratory Report LA 12601, MS, November, 1993. 35 Auger Design Report, March 14, 1997. See http://www-td-auger.fnal.gov:82/. 36 A. Borine et al., Nucl. Instrum. Methods A 346, 322 (1994). 14

  • Page 19

    1.15 Construction of a new high-precision rotating torsion balance E.G. Adelberger, J.H. Gundlach, B.R. Heckel, K. Mauldin, H.E. Swanson and P. Williams We are constructing a continuously rotating torsion balance designed to test the equivalence principle with at least 20 times better sensitivity than our previous work37. Among other things, this will allow us to test whether the 600 km/s velocity of our region of the universe with respect to the rest frame of the cosmic microwave background is predominantly due to ordinary gravitation. The 600 km/s velocity divided by a Hubble time corresponds to an acceleration of about 1.7 × 10 −10 cm / s 2 . Our new device will allow us to compare this acceleration for Be, Al and Ti test bodies with a sensitivity better than 2 × 10 −13 cm / s 2 . This will be sufficient to detect any significant non- gravitational (i.e. equivalence-principle violating) component of the acceleration. The new apparatus is located in the old cyclotron cave. This site has the advantage of being removed from traffic patterns, is thermally stable, and provides a strong topographic source for Yukawa interactions with ranges down to 1 m. The turntable for the new instrument consists of a 25 cm diameter aluminum air bearing powered by a direct- drive eddy-current motor. A Heidenhain 36000-line shaft encoder is directly mounted to the turntable rotor. This system was manufactured by the Professional Instruments Company in Minneapolis. We developed a rotation-rate controller based on a TMS320 digital signal processor that digitizes the sine and cosine outputs of the Heidenhain, computes the turntable angle, and compares it to the desired angle which was set by the digital pulses of a Stanford precision programmable oscillator. This error signal is then used in a digital PID loop to drive the turntable motor. This device is extremely successful, the error signal of the unloaded turntable is less 100 picoradians per root Hz for frequencies below 2.5 Hz. Signals from, and power to, the rotating apparatus are brought out via a 60-contact slipring assembly that is servoed to track the turntable so that the turntable does not provide the torque to rotate the sliprings. The vacuum chamber is completed and will be attached below the turntable. The torsion pendulum will be suspended from a small motorized rotation and vertical translation stage that is mounted on a two-axis gimbaled platform. The autocollimator and its optics will also be suspended from this gimbaled platform. A hermetic thermal shield consisting of double-walled aluminum cylinders that will have temperature- regulated water flowing in them has been completed. A µ -metal shield is mounted inside the shields. We have designed a symmetric pendulum that holds 8 test bodies (4 each of 2 kinds of material - either Be, Al or Ti). The entire pendulum has vanishing q 20 and q 40 mass moments. The barrel-shaped test bodies themselves have vanishing q 20 and q 40 moments, and are located very reproducibly in conical seats in the pendulum. A prototype has been built from aluminum; the final pendulum will be machined from Be. Currently we are testing the turntable with an old torsion balance that is temporarily suspended from a gimbal on the turntable. This setup will enable us to measure gravity gradients in the location of the final apparatus, so that the gravity gradient compensator masses can be rough-cut before the whole system is completed. 37 Y. Su et al., Phys. Rev. D 50, 3614 (1994). 15

  • Page 20

    1.16 Design and testing of autocollimators for the Eötwash experiments E.G. Adelberger, J.H. Gundlach, M.G. Harris, B.R. Heckel, S.D. Penn, G.L. Smith and H.E. Swanson We have designed, built, and tested a new version of autocollimator for use in the Eötwash experiments. This new model features an improved construction and an updated optical system. Our tests on the new autocollimator indicate a fourfold improvement in sensitivity over the old model. The Eötwash autocollimator is used to measure the angular deflection of our torsion pendulum. The autocollimator consists of a diode laser, an optical system, and a position sensitive detector (PSD). The laser is positioned at the focal point of a plano-convex lens. The laser light is reflected through a beam splitter to the lens, where the beam becomes parallel, and is then reflected off a series of mirrors and onto the pendulum’s mirror. The light reflected off the pendulum returns along a similar path but passes through the beam splitter to the PSD. In the old autocollimator we used a standard beam splitter which had two drawbacks. First we lost half our light with each pass through the beam splitter. This 75% light loss was very inefficient, and any spurious light reflected back onto the PSD added noise to the system. In addition half of the return beam was reflected back into the laser which is known to cause laser instability. In the new model we incorporated a polarizing beam splitter (PBS) which had an efficiency of > 98% for splitting the beam along its orthogonal polarizations. We also added a 1/4 waveplate between the PBS and the lens with its optical axis set at 45° to the normal of the optical plane. Thus by correctly aligning the laser polarization, nearly all of its light would be reflected to the lens while nearly all the return beam would be transmitted to the PSD. This eliminated our problems with reflections, back scattering, and light efficiency. The autocollimator design was enhanced in several ways to provide a more accurate and rigid alignment and greater ease of use. The laser holder was altered to include a rotational mount with locking ring to allow to an easy change in polarization without changing laser position. In addition a heat shield was added to guard against temperature variations thereby increasing laser stability. Both the 1/4 waveplate and the lens were mounted in a disk- shaped holder and prevented from slipping by a slightly compressed O-ring. The 1/4 waveplate holder also had a locking ring rotational mount to allow easy and accurate alignment of the optical axis. This locking ring was extended to project slightly from the block so as to rigidly align the block with the spindle piece, which houses the lens and connects the autocollimator to the vacuum can. The PSD holder was also fitted with a rotational mount and locking ring to set the detector angle. The PSD is a 1-D detector so it is important to keep its angular position fixed. However the complete PSD mount is composed of several pieces (mounting block, preamp box, and heat shield) which must sometimes be taken apart. Thus we also devised a way to remove one or all of these components without changing the rotational orientation of the PSD. The detector itself was mounted in a sliding bed with lock screws to allow for the small transverse motion to align the detector center on the optical axis. The PSD mount was also equipped with a scale to measure the detector angle. Between the PSD mount and the block was the stage block, a spacer which during calibration can be replaced by a translation stage. The stage block was modified slightly to allow the stage and its adapter blocks to be inserted as a unit rather than in three pieces as before. Also the stage blocks were equipped with protruding rims to align them on the block. An unfortunate consequence of the new design was that the focused spot on the PSD was too small. For best results the spot must be large compared to the gaps on the PSD. This problem was fixed by increasing the thickness of the stage block which placed the PSD slightly beyond the focal point. A newer version of the autocollimator is currently being designed. This new model will include a longer focal length (from 20 to 30 cm), transverse adjustment for the laser, longitudinal adjustment for the PSD, a more massive laser mount (to better heat sink the laser), a more compact heat shield for the PSD, and the ability to swap positions of the laser mount and detector mount. 16

  • Page 21

    1.17 A short-range test of the Equivalence Principle E.G. Adelberger, J.H. Gundlach, B.R. Heckel, G.L. Smith and H.E. Swanson We developed the Rot-Wash torsion balance to extend our existing tests38 of the equivalence principle (EP) in several ways. 1. to probe Yukawa EP-violating interactions with ranges down to 1 cm 2. to achieve high sensitivity for interactions that couple to charges such as N-Z for which our earth-attractor data was insensitive because the earth contains nearly equal numbers of neutrons and protons. 3. to “fill in” the gap (Yukawa ranges between 10 4 m and 10 6 m) where earth-attractor constants are not available because of the difficulty in calculating the horizontal component of a Yukawa force for these ranges (see reference. 1). As these results should appear in Physical Review Letters39 before this Annual Report is published, we reproduce below the abstract of our publication and refer the reader to that publication for details: We rotated a 3 ton 238 U attractor around a compact torsion balance and compared the accelerations of Cu and Pb toward U. We found that aCu − aPb = ( 0.7 ± 5.7 ) × 10 −13 cm / s 2 compared to the 9.8 × 10 −5 cm / s 2 gravitational acceleration toward the attractor. Our results set new constraints on equivalence-principle violating interactions with Yukawa ranges down to 1 cm and rule out an earlier suggestion of a Yukawa interaction coupled predominantly to N-Z. It is amusing to note that if an acceleration of 5.7 × 10 −13 cm / s 2 were applied to an object beginning at the time of Caesar and continued to the present day, the body would now be moving approximately as fast as the end of a minute hand on a clock. 38 Y. Su et al., Phys. Rev. D 50, 3614 (1994). 39 J.H. Gundlach et al., Phys. Rev. Lett. in press. 17

  • Page 22

    1.18 On the translation of multipoles for a 1/r potential E.G. Adelberger and C. D'Urso We outline an analytic method for calculating with high accuracy inner and outer multipoles of a 1/r potential about an arbitrary point in terms of known multipoles about a given point. This problem is frequently encountered in experimental gravitation where the multipole formalism is used to calculate torques and forces40,41 on test bodies in the fields of attractors with realistic geometries. Many cases of practical interest deal with distributions that can be accurately approximated as a linear superposition of basic geometric shapes. To calculate the inner and outer multipoles of these distributions it is necessary to know the multipoles about the point of interest for each of the superposed shapes. We showed that if the inner and outer multipoles of these simple shapes are known about any one point then analytic solutions can be obtained for the multipoles of a superposed distribution about an arbitrary point. Consider an arbitrary distribution of mass, charge, etc. that can be separated into two disjoint bodies. One body, the “object”, is considered to be in the field produced by the other body, the “source”. Integrating over this distribution, the object has inner and outer multipoles * qlm = ρ ( r )r1Ylm ( r )d 3 r Qlm = ρ s ( r )r -(1+1)Ylm ( r )d 3 r, ✚ ✚ ✘ ✙ ✛ ✙ (1) 0 * and Y are spherical harmonics and the density functions, ρ 0 (r ) and ρs(r ) , correspond to ✚ ✚ respectively, where Ylm lm the object and source bodies, respectively. Assuming that both qlm and Qlm are known about a coordinate origin, O, we now wish to know the corresponding inner and outer multipoles about an arbitrary point, P( r ′,θ ′, φ ′ ). To express the moments about P in terms of the moments about O we use relations found in D.A. Varshalovich et al.42 The inner and outer multipole moments about P, respectively, take the following forms: 4π ( 2L + 1)! q~lm = ∑ r' l' Yl*' m' ( r' )C(l', m', l, m, L, M )δ L,l + l ' qlm , ✙ (2) l',m' l,m ( 2l' +l )!( 2l + 1)! ~ 4π ( 2l )! Qlm = ∑ r' l' Yl*' m' (r' )C(l', m' , l, m, L, M )δ L ,l + l' Qlm , ✙ (3) l',m'l ,m ( 2 L )!( 2l ' +1)! where C is a Clebsch-Gordan coefficient. These expressions provide analytic solutions for the multipole moments about an arbitrary point in terms of known moments about a given point. The equation for q~LM is equivalent to one given in Reference 2. Here we mention selection rules that are useful in understanding the effects on the multipole moments due to perturbations in the position of the objects. For small displacements the leading order terms in Eqs. 2 and 3 have l' = 1. As a result, a small translation of an object with moments qlm and Qlm induces new moments q~l +1,m′ and ~ Ql −1,m′ where m ′ = m if the displacement is along z and m ′ = m ± 1 if the displacement is in the x − y plane. In ✙ second order in the displacement one induces new moments ql =2,m ′′ and Ql −2,m′′ where m ′′ = m for displacements along z and m ′′ = m,m ± 2 for displacements in the x − y plane. ✙ These expressions allow one to calculate, without numerical integration, inner and outer multipoles of complex systems. Perhaps more importantly, one can study the effect of perturbations of systems whose unperturbed multipole moments are known and, for example, to center precisely objects in the field of a source using purely gravitational means. Rotations of the spherical multipoles are easily carried out using standard techniques of angular momentum algebra. We addressed the problem of translations. 40 E.G. Adelberger et al., Phys. Rev. D 42, 3267 (1990). 41 Y. Su et al., Phys. Rev. D 50, 3614 (1994). 42 D.A. Varshalovich, A.N. Moskalev, and V.K. Khersonskii, Quantum Theory of Angular Momentum, World Scientific, 1988, p. 167. 18

  • Page 23

    2.0 NEUTRINO PHYSICS 2.1 The neutral current detector project at the Sudbury Neutrino Observatory (SNO) M.C. Browne, T.H. Burritt, P.J. Doe, C.A. Duba, S.R. Elliott, J.E. Franklin, J.V. Germani, K. M. Heeger, R. Meijer Drees, A. Myers, A.W.P. Poon, R.G.H. Robertson, M.W.E. Smith, T.D. Steiger, T. Van Wechel, J.F. Wilkerson and S. Zsitvay. SNO will detect Cerenkov light emitted from electrons or positrons produced by charged-current neutrino interactions. These measurements will provide a measure of the flux of electron neutrinos from the sun. Neutrinos of any flavor, however, can produce free neutrons in the heavy water by neutral current interactions. Thus the measurement of the neutron production is a measurement of the total flux of neutrinos from the sun. Since solar burning produces only electron neutrinos, a comparison of the total neutrino flux to the electron neutrino flux could provide strong evidence for neutrino oscillations and therefore neutrino mass. The neutral current detectors (NCD's) are 3 He filled proportional counters designed to detect such neutrons. These NCD's are made from CVD Ni tubing and endcaps. A quartz tube forms the high voltage and signal feedthrough to a Cu anode wire. We have received about 300 acceptable 2-meter CVD Ni tubes from Mirotech in Toronto, Canada. The remaining 150 tubes should arrive by June 1997. Recent assay results from LANL of production tubing sampled from standard handling during counter production shows the Th contamination level to be approximately 1 ppt. This level is acceptable in terms of neutral current background. The various parts needed to assemble endcaps are being fabricated. As they are produced, they are being shipped to IJ Research in Santa Ana, California and being assembled. Severe difficulties were encountered in the transition from prototype to production endcaps at IJ Research. Good solder seals between the silica insulator and the nickel parts were elusive. Although prototype endcaps had worked well, the initial production endcaps were not leak tight. A final assembly jig was delivered to IJ from UW which solved these glass-to-metal seal difficulties. Thus after a significant delay, IJ is beginning full production of endcaps and we received 25 finished endcaps in December. As of January 1997, they are metallizing all of the quartz feethroughs and preparing them for radioassay. Beginning in March, full production should be underway. We measured the temperature change that compromised the glass-to-metal seal on the endcaps. Previous estimates indicated that this temperature change could be as small as 40°C. The measurement indicates that it is actually closer to 100°C which is very good news for shipping logistics. The endcap delays led to a decision to store nickel tubes in the Index adit (northern Washington state) under 800 m.w.e. of rock to reduce cosmic-ray activation. The cleanroom in which the counters will be produced is now continuously monitored and is better than Class 1000 when in use. The production factory is ready to use once endcap production resumes. The gas handling system for mixing the counter gas is operational. We have begun making counters and they work quite well. This initial batch of counters is intended to be studied underground at the Waste Isolation Pilot Plant facility in southern New Mexico to verify that they are low in radioactivity. The production cable has arrived from South Bay Cable in Idyllwild, California and looks very good both electrically and mechanically. Samples have been characterized optically by UBC and have been radioassayed by direct counting. All parts for the delay lines have been fabricated and assembled. These have all been tested to verify the impedance and testing for microdischarge characteristics under high voltage will soon begin. Using production cable, we have assembled 6 wet-end terminations, which will soon undergo long-term testing under water. Deep Ocean Engineering has delivered the Mark II ROV to Los Alamos for pool tests this summer and fall. The microdischarge test system is working well and has been installed in the cleanroom. Microdischarge tests have been done for 115 endcap-days at 2200 Volts. Only 10 microdischarge events have been observed. Since the counters will operate at a lower voltage (~1800 V) and most microdischarge events can be identified on an event-by-event basis, this low rate indicates that the endcap electrical breakdown should not result in a significant background for the NC measurement. More details are given elsewhere in this report. Our present production schedule will have all detectors underground in Sudbury by early 1998, ready for deployment that year after a period of cooldown during which cosmogenic 56 Co decays away. 19

  • Page 24

    2.2 The SNO data acquisition system Q.R. Ahmad, J.C. Beck, Y. Chan,* S.R. Elliott, M.A. Howe, M. Isaac.* F. McGirt,† R. Meijer Drees, A.W. Myers, T.D. Steiger, T.D. Van Wechel, J.F. Wilkerson and M. Wilson Our group is responsible for providing the data acquisition (DAQ) system which will read out the signals from SNO's 9557 photomultiplier tubes. The electronics and DAQ system is designed to handle modest background rates in excess of 1 kHz and burst rates in excess of 1 MHz and to have essentially no deadtime, in case of a galactic supernova. The DAQ system is written in an object-oriented programming language (C++), supports the VME bus, and utilizes multiple VME embedded processors (Motorola MVME 167s) for control of the experiment. The SNO DAQ system supports the readout of the 9557 individual channels via 19 custom 9u “SNO Crates” each containing 16 32-channel front end cards and a crate trigger card. Each SNO Crate is controlled via a SNO Crate translator/controller card that connects the SNO Crate backplane to a 6u VME card located in a standard VME crate, termed an interface crate. The interface crate contains additional SNO translator/controller cards as well as VME embedded processors directing readout of each crate. Event information is shipped from these crates via VME-to-VME controllers to VME embedded processors located in the primary VME data crate which handle event building, data recording, and shipment of data to the surface. We are currently developing code using Metrowerks' CodeWarrior C++ compiler and the Symantec TCL 1.1.3 class library. The advances made in the general DAQ framework will be addressed in "Extensions to the SNO DAQ Software" (see Section 2.4). The UW is also responsible for testing the custom front-end electronics as they come off of the production chain, and the progress we've made towards the test-stand is discussed in the section "The SNO Electronics Production Testing System" (see Section 2.3). The last year has seen considerable development on the monitoring side of the DAQ; we have now designed and prototyped the mechanism whereby "live" data can be shipped from a VME-resident cpu to a UNIX workstation running a TCP/IP socket manager called a "Dispatcher", from where it can be further distributed to analysis and visualization programs. The Dispatcher program, developed and in use at CERN, has been installed at NPL and has been extended for use with Macintosh computers. Several prototype analysis and visualization programs have been written, including an HTML interface for automatic updating of a Web page, an event display, an electronics calibration program, and an extended version of the CERN analysis program PAW, tailored for SNO, which can receive events from the Dispatcher and can show updating plots as the data arrives. Our group has also provided code to enable the official SNO offline analysis package (SNOMAN) to be able to hook up to the Dispatcher. *Lawrence Berkeley National Laboratory, Berkeley, CA. † Los Alamos National Laboratory, Los Alamos, NM. 20

  • Page 25

    2.3 The SNO electronics production testing system Q.R. Ahmad, J.C. Beck, M.A. Howe, R. Meijer Drees and J.F. Wilkerson At the time of writing, the design of SNO's data acquisition electronics is complete and the boards have been sent out for manufacture. When this is done (sometime in early April), all components will need to be checked for proper functionality before being used in the detector. Since there are large numbers of certain types of components, SNO is developing an automated system to perform most of this testing. The University of Washington's SNODAQ group is responsible for writing the software that comprises this automated test system. It is designed to test the 350 Front-End Cards (FECs) which acquire all SNO data, and the 350 High-Voltage Cards (HVCs) which connect the FECs to the 9,557 photomultipler tubes of SNO. The system performs basic analog and digital tests, produces a log of all test results, and decides whether the card is working properly or not. If the card passes the ensemble of tests, then another automated routine will adjust the timing of the electronics for proper operation in the detector. If the card fails, it will be sent to another testing station with the testing log for repairs. The software for this automated test system is comprised of routines for two programs that run on the Apple Macintosh computer. The SNODAQ application, written in C++, is the interface program for the SNO electronics. One can perform basic operations to the FEC and HVC, and program more complicated tasks with additional C++ routines. The Userland Frontier scripting language controls the flow of the testing and generates the log files with test results. At the moment, the basic framework of this automated test system is in place. The code for the SNODAQ application has been modified to allow two-way communication between it and Frontier. We have written a group of Frontier scripts that performs a basic set of digital read and write tests on the FEC. These scripts have been useful in debugging the pre-production hardware. Routines which test the ADC pedestal and slopes for all 32 channels of the FEC have recently been added to the SNODAQ code. Soon work will begin on the automated timing setup routines. After that, testing of the electronics will take place at the TRIUMF accelerator facility in Vancouver. All of the software written for this test system will also be useful for multi-card testing and the final assembly of the electronics on-site in Sudbury. 21

  • Page 26

    2.4 Extensions to the SNO DAQ software Q.R. Ahmad, J.C. Beck, Y. Chan,* S.R. Elliott, M.A. Howe, M. Isaac.* F. McGirt,† R. Meijer Drees, A.W. Myers, T.D. Steiger, T.D. Van Wechel, J.F. Wilkerson and M. Wilson The Sudbury Neutrino Observatory will use a data acquisition program that has been under development in a generic form for several years by Frank McGirt and John Wilkerson. This program has been undergoing extensive modifications here at NPL. It now has a more Macintosh-like interface and is much more extensible. The SNO DAQ code is runs on Macintosh computers using C++ compiled with MetroWerks CodeWarrior. The object-oriented nature of C++ allows data acquisition hardware, crate controllers, and other hardware objects to be fully described and encapsulated into software objects which contain a complete interface to that particular piece of hardware. A large number of VME, Camac, and NuBus based hardware modules are supported. In the original program, the job of collecting, processing, and storing data was done with task objects which were written specifically for each experiment to control the interactions between multiple hardware modules. In some cases, these tasks also had self-contained plotting packages. As a result, these task objects were quite complex and were written by people who were intimately familiar with both the hardware and the underlying framework of the code. In an effort to simplify the writing of such tasks, the original program was converted to a data-flow model with a new user-interface that allows all objects to be visually represented on the screen as icons with input/output pads. In this model, each object is as simple and self-contained as possible and can either produce, modify, display, or store data packet objects. The flow of data is set up by drawing lines between input/output pads using the mouse. In this way, data analysis chains (i.e. tasks) can be built or modified in seconds simply by adding/removing objects and drawing or moving connection lines. Data processing objects are simple and reusable. For example, an object that does histogramming or plotting only has to be written once and can then be connected into any analysis chain. While not all of the hardware modules yet support this data-flow model, a number of data process objects have been written. The list includes one and two dimensional histograms, a disk file object, a plotter object, specialized SNO data filters, and even TCP/IP objects to send and receive data over the network. As more of the hardware modules are converted and more processing objects are written, it will become possible to build up non- trivial data collection and analysis chains without writing any code at all. *Lawrence Berkeley National Laboratory, Berkeley, CA. † Los Alamos National Laboratory, Los Alamos, NM. 22

  • Page 27

    2.5 SNO data acquisition software development Q.R. Ahmad, J.C. Beck, Y. Chan,* S.R. Elliott, M.A. Howe, M. Isaac.* F. McGirt,† R. Meijer Drees, A.W. Myers, T.D. Steiger, T.D. Van Wechel, J.F. Wilkerson and M. Wilson The SNO DAQ group is responsible, among other things, for developing the software needed to understand and debug the SNO Front End Cards. Early last year the 32 channel pre-production Front End Cards [FEC32] were delivered to the SNO Electronics group at the University of Pennsylvania (UPenn). Given the necessity to understand the electronics thoroughly in order to write software for it we spent many months at UPenn working closely with the Electronics group. During that time we developed the essential software tools required to interface to the hardware. These tools were developed under our "Object Oriented Software Bus" (OSB) architecture, which has been outlined in last years' Annual Report. The FEC32 resides in a custom SNOBus Crate, which has been designed to minimize the VME bus electronics noise. Two custom made translator hardware modules, XL1 and XL2 respectively, provide the link between the VME bus and SNOBus. The SNOBus Crate consists of a XL2, 16 FEC32s and a Crate Trigger Card. The XL1 resides in the VME Crate. In order to interface to the SNOBus Crate a software module called SNTR[SNO Translator] was developed. This enables the user to interface to any hardware members of the SNOBus Crate. We have used this software module to perform low level debugging of the pre-production cards. We have made the user interface intuitive and have made it flexible enough so that new tools could be added with ease. Recently the electronics group, after debugging the pre-production FEC32, has submitted the FEC32 for final production. The SNTR module will be used in our Test Stand System which is under development to test, verify and characterize all production electronics cards. We have also developed another software module called SMTC [SNO Master Trigger Card] which interfaces to the Master Trigger Card. This is now being used at UPenn to debug the card. Recently we received the Master Trigger Card from UPenn and we are in the midst of integrating it into our system. Very soon we will have all the primary electronics components at NPL and that will no doubt help us shake down the Electronics/DAQ system. In addition to writing software for debugging electronics we have also developed algorithms for taking data, building events and then writing those events to disk. Currently these software tools have been developed on 68K and Power PC stations. In the final system all primary DAQ activities will take place in embedded processors [eCPUs] which reside in VME. We are currently in the process of migrating these tools to eCPUs. *Lawrence Berkeley National Laboratory, Berkeley, CA. † Los Alamos National Laboratory, Los Alamos, NM. 23

  • Page 28

    2.6 Progress on the acrylic vessel for the Sudbury Neutrino Observatory P.J. Doe and the SNO collaborators December, 1996 saw the completion of the upper hemisphere of the acrylic heavy water containment vessel for the SNO. The upper hemisphere and chimney, comprising 60% of the total vessel, is now suspended in position from its 10 Vectran ropes. Construction of the lower hemisphere has begun. This represents a major milestone in the progress to a complete and operating observatory. Achieving this milestone proved to be far more difficult than anticipated; however, valuable lessons were learnt which are being applied to the construction of the lower hemisphere. SNO Institute personnel and physicists are responsible for: 1) positioning and surveying the acrylic panels from which the vessel is built prior to bonding them together; 2) postcuring (heat treating) the bond joints; 3) finishing joints (sanding and polishing) after the panels have been postcured; 4) maintaining cleanliness throughout construction. The company contracted to build the vessel (Reynolds Polymer Technology, Inc.) is responsible for bonding the panels together. The first four of the five rows of acrylic panels comprising the upper hemisphere were installed between November 1995 and April 1996. Problems were encountered while bonding the fifth and last row. These were attributed to unexpected technical difficulties and inadequate quality control of the bonding operations. To overcome this, a dedicated R&D team was formed from members of the collaborating institutes to develop techniques to carry out the necessary repairs to the bonds. In addition the QA/QC (Quality Assurance/Quality Control) program was reviewed and reorganized. For each of the construction activities listed above, an activity leader from the SNO institutes has been appointed. That person is responsible for developing the QA program for his activity and ensuring that it is rigorously adhered to. The activity leaders report to an overall QA activity leader for the vessel. These measures are an important component in ensuring the timely and successful completion of the acrylic vessel. By December 1996 the upper hemisphere and chimney of the vessel were complete and the process of removing all the scaffolding, jigs and fixtures begun. In January 1997, the upper hemisphere was raised 45 feet on its construction platform and suspended in its final position from 10 Vectran ropes. The special jigs and fixtures for constructing the lower hemisphere were then erected on this platform and the bonding of the first row of panels of the lower hemisphere has begun. The schedule calls for the completion of the acrylic vessel by 18 June 1997. Throughout this period, Peter Doe, as group leader for the acrylic vessel will continue to direct the various activity leaders and to function as activity leader for bonding and bonding R&D. Hardy Seifert, Visiting Scientist on the UW faculty from October 1996 to July 1997, has been working with Peter in Sudbury on the SNO bonding R&D team, and many of the SNO group at UW have spent some months on site at Sudbury working to overcome the thorny problems encountered with the vessel construction. 24

  • Page 29

    2.7 Development of a compact 20 MeV gamma-ray source for energy calibration at the Sudbury Neutrino Observatory M.C. Browne, R.J. Komar,* N.P. Kherani,† H.B. Mak,# A.W.P. Poon, R.G.H. Robertson and C.E. Waltham* We are developing a compact 20 MeV gamma-ray source for energy calibration at the Sudbury Neutrino 4 Observatory (SNO). The gamma-rays are produced in the 3 H p, γ He radiative capture reaction. The design and ✁ the operational characteristics of the source can be found in our previous reports.1 We have built the final source to be used in the SNO detector. In Fig. 2.7-1, a photograph of the source is shown. The tritium target used in the source was fabricated at the tritium facility at Ontario Hydro Technologies in Toronto, Canada. The target fabricated was a scandium tritide thin film on a molybdenum substrate. Prior to the target fabrication, the substrate was chemically cleaned and etched to enhance film adhesion. A scandium film of ~10000Å was evaporated. To ensure a good tritium-to-scandium stoichiometric ratio, the scandium film was tritiated in situ. The substrate was heated up to 400°C by an internal heater installed in the evaporation system. Tritium was let into the system and was pumped by the scandium film. Once the scandium film had been saturated by the tritium, the temperature of the substrate was slowly brought down to anneal the film. After the target fabrication and its mounting in the source, the source was brought to Queen's University in Kingston, Ontario for testing. A 12.7cm(Ø) x 7.62 cm BGO crystal with an active cosmic veto was used to detect the 20 MeV gamma-rays. A 12.7cm (Ø) x 5.1 cm liquid scintillator with active cosmic-ray veto was used to monitor the neutron production rate by pulse shape discrimination. The neutrons are produced from 3 H t, nn 4He ✁ ✂ and 3 H d, n 4He reactions. Analyses are being done on the data. Preliminary results of an eight-hour, 25 keV ✁ ✂ beam energy run show that in the energy window of 18 to 20 MeV, we see gamma-rays at 8.7σ above the background. In Fig. 2.7-2, a background subtracted energy spectrum for the eight-hour run is shown. This analysis indicates that the source generates ~1 gamma per second. Neutrons were also observed coming from the source. In the near future, we plan to make a gamma-ray angular distribution measurement with this source here at NPL. 30 3 4 SNO H( p,γ ) He Source 20 Data Counts / 128.3 keV Monte Carlo 10 0 -10 16 18 20 22 24MeV Energy 4 Fig. 2.7-1. A photograph of the 3 H p, γ He source. Fig. 2.7-2. Background subtracted energy spectrum ✁ The ion discharge magnet is removed here in order to of the eight-hour test run. The Monte Carlo spectrum display the source body. is superimposed on the data. * Department of Physics and Astronomy, University of British Columbia, Vancouver, BC, Canada V6T 1Z1. †Ontario Hydro Technologies, 800 Kipling Avenue, Toronto, Ontario, Canada M8Z 5S4. #Department of Physics , Queen's University, Kingston, ON, Canada K7L 3N6 1Nuclear Physics Laboratory Annual Report , University of Washington, (1994) p. 10; Nuclear Physics Laboratory Annual Report, University of Washington, (1995) p 22. 25

  • Page 30

    3.0 NUCLEUS-NUCLEUS REACTIONS 3.1 Preequilibrium particle emission and the GDR in Sn Nuclei M.P. Kelly, J.P. Lestone J.F. Liang, K.A. Snover, A.A. Sonzogni and J.P.S. van Schagen The proper analysis of GDR γ -ray spectra produced in heavy-ion collisions depends on the characterization of the hot nuclei produced in heavy-ion collisions. Our current GDR decay studies1 involve 10 MeV/nucleon projectiles. We report the results of a study of preequilibrium particle emission in 18 O +100 Mo collisions in the region of 10 MeV/nucleon bombarding energy and the effects of this emission on the excitation energy of the fused compound system. Measurements of light charged particles produced in 18 O +100 Mo collisions at 169 MeV and 200 MeV bombarding energy and measurements of evaporation residues for the same reaction from 100 MeV to 217 MeV bombarding energy were performed at the University of Washington Nuclear Physics Laboratory using the FN Tandem Van de Graaff as injector for the Superconducting Linear Accelerator. These studies allow us to accurately determine the average excitation energy and fusion cross section for compound nuclei formed in 18 O +100 Mo collisions. Analysis of particle spectra for 18 O +100 Mo at 200 MeV bombarding energy indicate that on the average, 0.27+0.06 protons and 0.45+0.07 α -particles are lost due to preequilibrium emission. We infer from comparison with model calculations that 0.87 neutrons are lost as well. As a result, the average excitation energy of the fused system is decreased by approximately 20% relative to complete fusion. In addition, evaporation residue measurements show that the default fusion cross section often used in the statistical model code CASCADE2 underpredicts the true (complete + incomplete) fusion cross section by nearly 20% at the same energy. Results of statistical model calculations incorporating the experimentally determined fusion cross sections and excitation energies indicate that both the deduced GDR width and strength are significantly affected by the energy loss due to preequilibrium emission, while the deduced centroid energy remains relatively unchanged. In fact, it was found that if one neglects the effects of preequilibrium emission in 18 O +100 Mo collisions at 200 MeV bombarding energy, one determines both a GDR width and a strength that are 15%-20% lower than those deduced from a proper analysis which accounts for preequilibrium emission. A more detailed account of these studies is being submitted to Physical Review C. 1 Nuclear Physics Laboratory Annual Report, University of Washington (1996) p. 25. 2 F. Puhlhofer, Nucl. Phys. A 280, 267 (1977). 27

  • Page 31

    3.2 Quasi-elastic and fusion barrier distributions for 40 Ca + 46,48,50 Ti J.F. Liang, J.P. Lestone, M.P. Kelly, A.A. Sonzogni and R. Vandenbosch It has recently been suggested that fusion information can be derived from quasi-elastic scattering at near- grazing angles.3 Since fusion excitation functions were measured and barrier distributions derived for 40 Ca + 46,48,50 Ti ,4 an experiment to obtain quasi-elastic excitation functions was performed to compare the different barrier distributions results. The quasi-elastic events were detected with a set of Si detectors, placed at angles close to the grazing angle. Energy and time of flight information was used to distinguish between the different classes of events. The quasi-elastic barrier distributions were obtained by using the relation:1 ✄ dσ qel D qel (E) = − d / dE ✁✂ ☎ . The results can be seen in Fig. 3.2-1, together with the fusion results. A good ✆ dσ Ruth agreement between them is seen. Because of these results and those of references 1 and 3, one concludes that this technique can provide barrier distributions which closely resemble those from fusion data. This can be usefully exploited when trying to probe fusion enhancement with radioactive beams, where the beam flux is too small to be used in fusion cross section measurements at near-barrier energies. Fig. 3.2-1. Fusion and quasi-elastic barrier distributions for 40 Ca + 46,48,50 Ti . 3 H. Timmers et al., Nucl. Phys. A 584, 190 (1995). 4 Nuclear Physics Laboratory Annual Report, University of Washington (1996) p. 27. 3 H. Timmers et al., Instituto Nazionale di Fisica Nucleare, Laboratori Nazionali di Legnano, Annual Report (1995) pp. 88- 90. 28

  • Page 32

    3.3 Projectile mass dependence of fission fragments anisotropies for heavy-ion induced fission on 232 Th at near-barrier energies J.F. Liang, J.P. Lestone, M.P. Kelly, A.A. Sonzogni and R. Vandenbosch Fission fragment angular distributions and angular correlations were measured for 16,18 O + 232 Th at near barrier energies. A brief description of the experimental technique can be found in Section 3.4 of this report. The folding angle distributions were used to separate a transfer-fission component. The percentage of transfer-fission to total fission was similar for both oxygen isotopes (~8-12% at near barrier energies). The anisotropy results can be seen in Fig. 3.3-1, together with those of 19 F + 232 Th ,5 11 B,12 C + 232 Th (see Section 3.4). It can be seen that for a given bombarding energy relative to the Coulomb barrier, the anisotropy increases smoothly with projectile size. No discontinuity in the fission dynamics, due to the mass asymmetry of the system, can be seen in this plot. Fig. 3.3-1. Fission fragment anisotropy as a function of E c.m. / Vb for 11 B, 12 C , 16,18 O and 19 F + 232 Th . 5 N. Majumdar et al., Phys. Rev. C 51, 3109 (1995). 29

  • Page 33

    3.4 Fission fragment anisotropies in the 11 B, 12 C+ 232 Th reactions M.P. Kelly, J.P. Lestone, D.J. Prindle, A.A. Sonzogni and R. Vandenbosch Majumdar et al.6,7 have recently published measurements of the anisotropies of fission fragments following fusion reactions with 11 B, 12 C , 16 O and 19 F projectiles on 232 Th. These measurements show unexpected peaklike structures in the anisotropies as a function of center of mass energy in the near- and sub-barrier energy regions. The 16 O and 19 F measurements are made difficult by the presence of a significant yield of fission following transfer reactions at near- and sub-barrier energies. However, for reactions involving 11 B and 12 C projectiles on Th, the transfer-fission yield is small enough that the anisotropy of fragments following full momentum transfer can be obtained by simply measuring the singles fission fragments at all but the lowest beam energies.1,8 We have measured the angular distribution of fission fragments in the 11 B+ 232 Th and 12 C + 232 Th reactions, using beams from the University of Washington superconducting linac. The target consisted of a 225µg/ cm 2 layer of 232 ThF4 evaporated onto a 100 µg / cm 2 Ni foil. Fission fragments were detected with Si surface barrier telescopes and identified using energy and time-of-flight information. Folding angle distributions were obtained by observing fission fragments in Si telescopes and the complementary fragments in a large area position sensitive Si detector. Our folding angle distributions contain a single peak and confirm that the fission yield in 11 B, 12 C + 232 Th reactions is dominated by fission following complete fusion. Fig. 3.4-1 compares our fission fragment anisotropies to those of other authors. Our anisotropies vary smoothly with center of mass energy and give no hint of the peaklike structures reported in References 1 and 2. In view of our results we feel that the conclusions drawn in Reference 1 should be viewed with caution. Fig. 3.4-1. Fission fragment anisotropies for the 11 B, 12 C + 232 Th reactions. The solid circle shows the present work; (open squared);1,2 (open circles);4 and (open triangles).5 To avoid over crowding of this figure only representative error bars have been shown for the 12 C + 232 Th data of Majundar et al.,2 Karnik et al.9 and Ramamurthy et al..10 The up arrows show the positions of the B+Th and C+Th fusion barriers. 6 N. Majumdar et al., Phys. Rev. Lett. 77, 5027 (1996). 7 N. Majumdar et al., Phys. Rev. C 53, R544 (1996). 8 J.P. Lestone et al., Phys. Rev. C. 55, R16 (1997). 9 A. Karnik et al., Z. Phys. A. 351, 195 (1995). 10 V.S. Ramamurthy et al., Phys. Rev. Lett. 65, 25 (1990). 30

  • Page 34

    3.5 Search for an entrance channel dependence of fission anisotropies M.P. Kelly, J.P. Lestone, A.A. Sonzogni and R. Vandenbosch In recent years several authors claim to have seen an entrance channel dependence of fission anisotropies [A = W(180° ) / W(90° )] for target projectile combinations across the Businaro-Gallone (BG) ridge in the mass asymmetry degree of freedom.11,12 To remove any effects due to the difference in fissility of different compound nuclei, it is necessary to study reactions that form the same compound nuclei with different entrance channel mass asymmetries. At above barrier energies, Vandenbosch et al.13 found no entrance channel dependence of the fission anisotropies for 248 Cf compound nuclei formed in the two reactions 12 C + 236 U and 16 O + 232 Th which span the BG ridge. We present here a study of these two reactions at near- and sub-barrier energies. Fission fragment cross sections, angular distributions and folding angle distributions were measured using 12 C and 16 O ions from the University of Washington superconducting linac. In the 16 O + 232 Th reaction, corrections due to transfer fission were made using our fission fragment folding angle distributions. Our 12 C + 236 U folding angle distributions are consistent with the fission yield being dominated by fusion-fission events. To remove the influence of the differing amounts of angular momentum brought in by the projectiles and the different fusion barriers (VB ) for the 12 C + 236 U and 16 O + 232 Th reactions the quantity 4(A exp − 1)K 20 / < J 2 > has been estimated as a function of the center of mass energy relative to the fusion barrier E c.m. / VB using the present experimental results and those of other studies3,14 (see Fig. 3.5-1). The mean square spin of the compound nuclei < J 2 > was estimated as a function of the center of mass energy using calculations that reproduce our measured fission cross sections. Values of 4(A exp − 1)K 20 / < J 2 >≠ 1 indicate a departure of the experimental anisotropies from the values predicted by the transition state model. Our analysis suggests the anomalous behavior of anisotropies increases smoothly as the center of mass energy drops through the region of the fusion barrier independent of the entrance channel mass asymmetry. Fig. 3.5-1. 4(A exp − 1)K 20 / < J 2 > as a function of E c.m. / VB . † The 16 O + 232 Th point of Reference 4 has been corrected using the transfer fission measurements of Lestone et al.15 11 V.S. Ramamurthy et al., Phys. Rev. Lett. 65, 25 (1990). 12 N. Majumdar et al., Phys. Rev. Lett. 77, 5027 (1996). 13 R. Vandenbosch et al., Phys. Rev. C 54, R977(1996). 14 T. Murakami, et al., Phys. Rev. C 34, 1353 (1986). 15 J.P. Lestone et al., Nucl. Phys. A 509, 178 (1990). 31

  • Page 35

    3.6 Influence of target spin on sub-barrier fission fragment anisotropies M.P. Kelly, J.P. Lestone, A.A. Sonzogni and R. Vandenbosch Liu et al.16 recently observed that the sub-barrier fission anisotropies in the 11 B+ 237 Np reaction are substantially lower than the corresponding anisotropies in the 16 O + 232 Th reaction. It is claimed that this entrance channel effect is due to the Businaro-Gallone ridge in the mass asymmetry degree of freedom. This conclusion is, however, not supported by a similar comparison of 12 C + 236 U and 16 O + 232 Th fission anisotropies (see Section 3.5). Another possible explanation of the large differences in the 11 B+ 237 Np and 16 O + 232 Th anisotropies is the 5/2 spin of 237 Np and the 3/2 spin of the 11 B. At beam energies well above the fusion barrier the angular momentum brought in by the projectile in sufficiently large that the spin of the target and/or projectile can be neglected. At sub-barrier energies where the fusion spins are smaller it is possible for target/projectile spins to have a noticeable effect. To test this possibility we measured the fission anisotropy for the three reactions 12 C + 235,236,238 U , as a function of center of mass energy (see Fig. 3.6-1). Measured folding angle distributions for all three reactions confirm that the fission yields are dominated by fusion-fission events. At above barrier energies the three different 12 C + U reactions have comparable anisotropies. At sub-barrier energies the 12 C + 236,238 U reactions have anisotropies of ~1.6, while the fission anisotropies for the 12 C + ( 7 / 2 spin) 235 U reaction decrease with decreasing center of mass energy to values well below 1.6. This is strong evidence for an influence of target spin on sub-barrier fission anisotropies. The magnitude of our observed sub-barrier target spin effect requires a preferential reaction with the tips of the deformed target nucleus at low energies and some remembrance of the K of the target ground state. Fig. 3.6-1. Fission fragment anisotropies of the three reactions 12 C + 235 U, 236 U and 238 U . The open squares, closed triangles and stars show the present results with 238 U, 236 U and 235 U targets respectively. The solid squares and circles show the 236 U target results of Murakami et al.17 and Vandenbosch. et al.18 The up arrow shows the position of the C+U fusion barrier. 16 Z. Liu et al., Phys. Lett. B 353, 173 (1995). 17 T. Murakami et al., Phys. Rev. C 34, 1353 (1986). 18 R. Vandenbosch et al., Phys. Rev. C 54, R977 (1996). 32

  • Page 36

    4.0 ULTRA-RELATIVISTIC HEAVY IONS 4.1 Event simulation of high-order Bose-Einstein and Coulomb correlations J.G. Cramer and D.D. Weerasundara Last year we reported1 the initial development of a new Monte Carlo event generator program for simulating the high-order multiparticle correlations between pions produced in ultra-relativistic heavy ion collisions. The program uses a new Monte Carlo elimination algorithm that employs the multiparticle correlation formalism developed in a recent paper2 to impose Bose-Einstein and Coulomb correlations of up to 6th order on the π ± particles of a simulated event. Tests of a preliminary version of this code in event simulations showed an unexpected “Bose-Einstein cooling” effect in the generated particle spectrum, with increasing correlation order producing a dramatic decrease of the particle spectrum's inverse-slope or “temperature” parameter. Further, the 2-particle correlations extracted by analyzing the simulated events showed a tendency toward runaway to Bose-Einstein condensation as the correlation order used in the simulation was increased from 2 to 6. Since no such runaway effect is apparent in NA49 data (which should contain true BE and Coulomb correlations up to order 250 or so), there was a puzzling discrepancy between theoretical simulation and experimental data. The root of these problems was recently discovered to be the basic ansatz that the full multiparticle correlation of a “candidate” particle with its nearest neighbors in momentum space is proportional to the probability of acceptance or “wave function collapse” of that particle. This is discussed in another article (see Section 4.2) in which it is shown from quantum mechanical arguments that a new “reduced correlation function” must be used to determine the acceptance or rejection of candidate particles in the rejection algorithm. An improved Monte Carlo event generator program sim_hbt_10 has been developed and is now debugged and running. It includes several significant improvements over the old code: 1) it simulates Bose-Einstein effects in the rejection algorithm using partially-coherent reduced correlation functions generated by Mathematica; 2) it uses Coulomb correlation functions obtained from 6-dimensional Monte-Carlo integration of Coulomb wave functions over the source (see Section 4.3) to take into account the finite source size; 3) when N particles are used in the BE correlation it now includes about 2N-1 particles in the Coulomb correlation; 4) it explicitly includes the effects of the pion “halo” from long-lived resonances decaying well beyond the surface of the primary pion freezout source; and 5) it explicitly includes the effects of mis-identified “contaminant” particles. The first calculations with this new code indicate that the unphysically strong Bose-Einstein effects described above are no longer present. In the new calculations there is no apparent tendency toward BE condensation runaway in the system as the correlation order is increased. We are continuing these investigations. We will soon use the new code to produce a library of simulated events under various conditions of source multiplicity, size, resonance fraction, and contamination that approximate the conditions of experiment NA49 and the planned STAR experiment. 1 Nuclear Physics Laboratory Annual Report, University of Washington (1996) p. 37. 2 J.G. Cramer and K. Kadija, Phys. Rev. 53, 908 (1996). 33

  • Page 37

    4.2 Wave function collapse: Schrödinger's Cat and HBT simulations J.G. Cramer In an ultra-relativistic collision between heavy nuclei leading to the emission and detection of a large number (M) of pi mesons, the quantum mechanics of Bose-Einstein statistics requires that when the measurement of the system is made by the detector (e.g., the NA49 time projection chambers) the full M-pion wave function of the event must collapse to a fully symmetrized M-particle state. On the other hand, in a cascade-type Monte-Carlo simulation that attempts to include Bose-Einstein effects (see Section 4.1) the particles can only be produced and accepted sequentially, one at a time. This is the quantum-mechanical equivalent of making a new measurement each time a particle is produced. Parts (a) and (b) of Fig. 4.2-1 illustrate the difference between these two quantum mechanical situations, with horizontal lines representing measurements and downward arrows representing created particles. This difference leads to a version of the Schrödinger's Cat paradox: repeated measurements of the cat's condition (by repeatedly peeking into the cat box) lead to a different quantum condition from that which would have existed if the cat was left unobserved until the final measurement. In computer simulations we are forced to “peek into the cat box” with the production of each new particle, thereby distorting the quantum result. This qualitative difference between the simulation and the actual quantum condition has two important consequences: 1) it affects the choice of the correlation function used in selecting candidate particles; and 2) it introduces an order-asymmetry in particle treatment that must be “annealed out'' by the destruction as well as creation of particles. (a) Collapse of all particles simultaneously (c) Successive collapses and annihilations in sliding window Measurement Measurement X Measurement X Measurement X (b) Successive particle collapses in cascade simulation Measurement X Measurement X Measurement Measurement X Measurement Measurement X Measurement Measurement X Measurement Measurement . . . . . . Measurement Measurement Measurement Measurement Fig. 4.2-1. Comparison of wave function collapse modes. 34

  • Page 38

    First let us consider the choice of the correlation functions. The multiparticle momentum-space correlation function RM of order M (i.e., the correlation between M identical particles) is defined by the relation: ρ M (p1,... pM ) RM (p1,... p M ) = , (1) ρ 1 (p1 )...ρ 1 ( pM ) where ρ M (p1 ,... pM ) is the inclusive probability density for M particles expressed as a function of the three- momenta pi of the correlated particles, and ρ 1 (pi ) is the single particle probability density of the ith particle. This function can be interpreted as the relative probability of measuring a particular momentum-space configuration (p1 , p2 ,... pM ) of the M particles. If we were able to collapse the full wave function of an M particle event with a single measurement, then RM could be used to select among full M-particle candidate events in a Monte-Carlo rejection algorithm. This, however, is not computationally feasible. In Monte-Carlo simulations, in effect, we “measure” or select candidate particles one at a time in the presence of other pre-existing previously measured particles. In this circumstance it is not appropriate to use RM as the probability of the resulting momentum-space configuration because it include contributions from pre-existing particles that do not “connect” with the candidate particle. The relative probability of a new momentum space configuration with the new particle added is the difference between new M particle configuration and the old M-1 particle configuration. This leads to a “reduced correlation function” SM , as defined by: SM (i, j,k..., M) = 1 + RM (i, j,k..., M) − RM − 1 (j, k,..., M), (2) where i is the index of the candidate particle and j,k ...,M are the indices of its M-1 nearest neighbors. Using the reduced correlation function in selecting among candidate particles insures that the spurious probability enhancements from clusters that are distant or of opposite charge do not influence the correlation generated. We note that the when all particles have the same momentum, the reduced correlation function reaches its maximum value of SM (max) = M! - (M-1)! + 1, which has the values 2, 5, 19, 97, and 601 for M = 2, 3, 4, 5 and 6, respectively. Simulations of Bose Einstein effects should use this distribution function in accepting or rejecting particles. To deal with the second problem mentioned above, that of annealing out the order-asymmetry in particle treatment, we have devised a “sliding window” technique illustrated in part (c) of Fig. 4.2-1. For an event of multiplicity M, we first pre-fill the momentum space with M-1 uncorrelated particles. The generation of the Mth particle and succeeding ones then proceeds by “killing” one old particle as each new particle is accepted. In effect, each new particle is only correlated with particles within a “sliding window” of width M-1. This procedure is repeated for 2M+1 steps, so that a total of 3M particles have been generated. However, only the last M particles generated are retained in the event. We have implemented both of these procedures in a new program and are currently testing them. The initial results are very encouraging and indicate that Monte-Carlo simulation leading to a good approximation of a state with a fully symmetrized M particle wave function is feasible. 35

  • Page 39

    4.3. Finite Coulomb size and finite resolution effects in NA49 mixed-charge pion correlations J.G. Cramer In ultra-relativistic 160 GeV/nucleon heavy ion collisions between lead nuclei at CERN, the size of the pion source produced is large enough (6-8 fm) that the pointlike Coulomb interaction between pions (i.e., a Gamow penetrability) must be modified to take into account the finite size of the source of pions. Therefore, it is important to determine the Coulomb radius of the source before generating the Coulomb corrections used in performing Hanbury-Brown-Twiss analysis on experimental data from experiments like NA49. Partly to investigate this issue, an analysis of the π + π − correlations has been made by the NA49 collaboration from data measured with the NA49 vertex TPCs.3 The π + π − correlation offers the advantages that it is an enhancement at low relative momentum (while the π − π − correlation is a suppression) there is no Bose- Einstein effect between non-identical particles to modify the Coulomb effects, and the attenuation of the Coulomb effect by pions from resonances are about the same in π + π − and π − π − Coulomb correlations. The effect of the finite source size on the Coulomb correlation between charged pions can be calculated by integrating the two-particle Coulomb wave function of the pions over the spatial distribution of the source. The programs pipibig and pipipm were written by Pratt4 to calculate the finite-size Coulomb effects for π − π − and π + π − correlations, respectively. The programs as written, however, did not allow for the uncertainty in momentum determination in experimental data, which smears out and reduces the strength of the Coulomb effects at low relative momentum. We have therefore modified the Pratt programs to include the resolution effects. A sample of the π + π − correlation data of Appelshauser was fitted with the modified pipipm program, varying the Coulomb radius and the resolution width to obtain an optimum fit. The results are shown in Fig. 4.3-1, which plots the ratio of experimental data points to values from the pipipm calculation, with the parameters RCou = 6.7 fm and ∆p = 6.0 MeV. A perfect fit would be a raio of 1.0. As can be seen, the deviations are very small. Fig. 4.3-1. Ratio of π + π − correlation data to calculated finite-size Coulomb correlation. The curve through the points is the Coulomb correlation without resolution smearing. The upper line is the Gamow penetrability. 3 H. Appelshauser, Ph.D. thesis, University of Frankfurt, Germany (1996). 4 S. Pratt, private communication (1992). 36

  • Page 40

    4.4 Phase space filling systematics and superradiance in ultrarelativistic heavy ion collisions J.G. Cramer Collisions of ultra-relativistic heavy ions produce a spatial volume that is for a time densely populated with Bose-Einstein particles, i.e. pions. It is therefore interesting to consider whether the number density in position- momentum phase space is large enough to produce the dramatic phenomena associated with Bose-Einstein statistics: superradiance, stimulated emission, lasing, and Bose-Einstein condensation. Bertsch5 has shown that under the assumptions of a simple but plausible model for heavy ion collisions the filling factor < f > for cells of volume 3 in position-momentum phase space for a given collision event is: ✁ π dNπ − ( c)3 ✁ <f>= (1) 2 dy R' Rs' Ro' T 3 ✂ Here N − / dy is the multiplicity of π − particles per unit rapidity, R' Rs' , and Ro' are the longitudinal, side, and ✂ π out radii of the pion source, respectively, using the NA35 convention for source radii, and T is the characteristic “temperature” or (inverse-slope parameter) of the transverse mass distribution of the pions. Replacing dN − / dy π with the charged particle multiplicity dNch / dy and using the commonly accepted NA44/NA49 convention for source radii ( R' (NA35) = R(NA49)/ 2 ), expression becomes: π dNch ( c)3 ✁ < f >= (2) 8 dy R Rs RoT 3 ✂ Experiments using HBT interferometry on a number of colliding heavy ion systems over a wide range of particle masses and collision energies have revealed certain systematics for the colliding systems: (1) the characteristic temperature T ranges between 150 and 190 MeV; (2) the side and out radii are approximately equal, while the longitudinal radius is perhaps 20% larger; (3) defining the transverse radius Rt as the average of Rs and Ro , over a wide range of multiplicities Rt ≅ 0.9[dNch / dy ]1/ 3 . If we accept this relation and take R / Rt ≅ 1.2, then ✄ the phase space filling factor becomes < f > ≅ 0.24[ c / T ]3 . For a typical mid-range temperature of T = 170 ✁ MeV, < f > ≅ 0.38 . This value scales as [(170 MeV)/T]3 at other temperatures. This value of < f > ≅ 0.38 is an extremely large filling factor. It is to be compared with Bertsch's value for the onset of superradiance which is < f > crit ≅ 0.11 to 0.13. However, < f > ≅ 0.38 is probably an overestimate because not all of the pions come directly from the source. It is estimated that about 25% of the pions observed in a lead-lead collision at CERN come from long-lived resonances that decay well away from the collision point and would not participate in the phenomena of Bose-Einstein statistics. Correcting for this resonance fraction gives < f > ≅ 0.28 , still more than twice the calculated value of < f > crit . This estimate indicates that heavy ion collisions over a wide range of conditions produce Bose-Einstein systems which have number densities leading to superradiance and stimulated emission. We are investigating experimental analysis techniques for observing these effects. 5 G.F. Bertsch, Phys. Rev. Lett. 72, 2349 (1994). 37

  • Page 41

    4.5 Coulomb radii for STAR single event physics? J.G. Cramer The STAR detector at RHIC will begin operation in 1999. One analysis theme of the new detector will be a program of “event-by-event” physics, in which events exhibiting some unusual characteristic are grouped in an ensemble that is then subjected to further analysis, in search of unusual physics. It is expected that each Au+Au collision at RHIC will produce several thousand charged pions, and it remains an open question whether this particle multiplicity is sufficient to gain some indication of the characteristic source size of a single event through correlation analysis. Last year we examined the question of whether two-particle or three-particle Bose-Einstein correlations provide the best indicator of single-event source size6 and showed that two-particle correlations are always statistically superior. Now we are considering a similar comparison of π + π − vs π − π − correlations as indicators of single-event source size. In the RHIC environment it is expected that the size of the source will be at least twice as large as that found in the lead-lead collisions at CERN, i.e., in the range of 12 to 16 fm. For sources of this size, the like-charge Bose-Einstein momentum-space correlation for π − π − and π + π + is forced into the small relative momentum region where it is strongly suppressed by Coulomb repulsion.7 On the other hand, mixed-charge π + π − correlations offer several distinct advantages for estimating single- event source size: 1) The number of opposite-charge pairs is n 2 vs. n(n-1) like-charge pairs, giving slightly better statistics; 2) the π + π − Coulomb correlation is a strong enhancement which is guaranteed to have better statistics than the like-charge correlation at the lowest relative momentum; 3) the π + π − go to well separated parts of the TPC, so there is no need for a low-q cut when the tracks merge and are unresolved; and 4) while the π + π − Coulomb correlation is somewhat suppressed by Coulomb size effects (see Section 4.2), that suppression is less severe than that of the like-charge Bose-Einstein correlation. Fig. 4.5-1 compares the predicted π + π − correlation for Coulomb radii of 0 (Gamow), 6 fm, 12 fm, and 18 fm. We take this as an indication that size estimation from single event π + π − Coulomb correlations may be possible. In the coming year we will use the simulation code sim_hbt_10 (see Section 4.1) to investigate this question further. Fig. 4.5-1. Comparison of π + π − correlations with RCou = 0 (high), 6, 12, and 18 fm (low). 6 Nuclear Physics Laboratory Annual Report, University of Washington (1996) p. 38. 7 J.G. Cramer, Phys. Rev. C 43, 2798 (1991). 38

  • Page 42

    4.6 Application of SControl to the STAR TPC gas system J.G. Cramer and M.A. Howe The time projection chamber (TPC) for the STAR detector is nearing completion at LBL and will be transported by air to Brookhaven for installation at RHIC in October-November of this year. One critical component of the STAR TPC is the gas system, which supplies a mixture of helium and ethane to the TPC and monitors gas purity, pressure, temperature, and other relevant parameters. We are adapting the program SControl8,9 developed as the NA49 slow control monitor and alarm console for use with the STAR gas system. The project involves developing a version of SControl that will operate on a Sun workstation and linking the TPC gas system console, a PC running under Windows 95, the EPICS system which is the NA49 standard for slow control, and SControl as the top-level system communicating with the PC and the overall EPICS system. We plan to test this system during the cosmic ray tests of the STAR TPC, which will be in progress from May-September, 1997. 8 Nuclear Physics Laboratory Annual Report, University of Washington (1996) p. 45. 9 Nuclear Physics Laboratory Annual Report, University of Washington (1995) p. 57. 39

  • Page 43

    4.7 Event-by-event physics at the SPS and RHIC S. J. Bailey,* D.J. Prindle, J.G. Reid, C. Reynolds, T.A. Trainor and D.D. Weerasundara Event-by-event physics in nuclear collisions has been an active field of research for more than a decade. Jet production in p-p and e-e collisions and flow in heavy-ion collisions have been extensively studied at LEP, CSPS, TEVATRON, BEVALAC and AGS among other facilities. Jets and flow are large-amplitude and/or large- scale effects which are detectable despite the rather small particle multiplicities per event from these collision systems. Multiplicities have been limited by small system size (p-p, e-e) or low energy densities (heavy ions). Both jet and flow phenomena are manifestations of symmetry reduction or increased correlation with respect to a nominally thermalized system of produced particles. With the recent availability of a 158 GeV/nucleon lead beam at the CERN SPS one obtains (due to larger system size and increased energy density) substantially increased event multiplicities which can provide sensitivity to smaller-amplitude symmetry reductions over a finite scale interval. This increased statistical power permits event-by-event study of global thermodynamic variables and the possibility to extend this analysis program beyond global or large-scale variables to a scaling analysis approach. This additional analysis capability is needed to explore fully so-called soft or nonperturbative QCD phenomena. NA49 event-by-event analysis has proceeded on two complementary fronts: large-scale or global event- variable analysis -- based primarily on a thermodynamic approach, and scaled correlation analysis of momentum space which attempts to extend the event-by-event concept over a range of scales limited only by event multiplicity. There is also a flow analysis program underway. Examples of global variables are < Pt >, K / π ratios, m t -spectra slope parameters and some measure of the width of the rapidity distribution about the CM (for each species). These quantities all measure in some way the ‘boundedness’ of the phase-space distribution in momentum and flavor composition. One looks for deviations of the distributions of these variables from an event-mixed or other reference population or for correlations among global variables. With the larger collision systems and energies at SPS and RHIC, event multiplicity (statistical power) becomes sufficiently great to extend global-variable analysis by using more detailed measures of the correlation content of the phase-space distribution. At the SPS the total pion multiplicity for 158 GeV/nucleon Pb-Pb collisions is more than 1000, making it worthwhile to pursue a more detailed analysis for this particle species. At RHIC, multiplicities of 5000 or more in the STAR acceptance are expected. The UW group has pioneered a scaled correlation analysis technique for extracting all available information from a multiparticle distribution. The basic measure of this technique is scaled dimension transport, an elaboration of standard entropy-based topological measures on point sets. This technique has been applied to the STAR detector trigger algorithm production, jet-finding and electromagnetic calorimeter analysis, and NA49 event- by-event analysis (see Section 4.8). In the last case a population of anomalous events at the part-per-mil level has been detected in a sample of 200k events, illustrating the sensitivity of this model-independent approach. * Now at Harvard University, Cambridge, MA. 40

  • Page 44

    4.8 Scaled correlation analysis applied to NA49 main TPC data D.D. Weerasundara and T.A, Trainor QCD theory predicts the formation of color deconfined quark matter in the ultrarelativistic heavy ion collisions at the CERN SPS, RHIC and LHC energies. Event-by-event analysis of multiparticle distributions produced in these collisions is widely recognized as a promising way to look for the formation of color deconfinement. We have developed a multiparticle event-by-event scaled correlation analysis based on Renyi entropy10 to extract as much information as possible to search for rare phenomena produced in these collisions. The NA49 experiment provides an ideal test bed to further develop this correlation analysis. As part of the preparation for the large-scale DST production of NA49 data, we have analyzed ~200,000 events of main TPC (MTPC) data during 1996, using the large scale production computing facility (SHIFT) at CERN. The production of this DST sample serves two purposes. First it tests the robustness of the event reconstruction software, helps to develop calibration data bases and refine the DST production scripts. Secondly and most importantly, it provides a sufficiently high-statistics data sample to develop and test various physics analysis software. In order to have a manageable data volume (a few gigabytes) for the physics analysis, we have created from these DSTs, a set of mini-DSTs containing a minimal set of essential information relevant to our correlation analysis, which is being carried out at the UW-NPL HP-UX computer facility (see Section 4.13). We have performed scaled correlation analysis on the transverse mass spectra extracted from particle trajectories reconstructed using the MTPC-only data. For this analysis, we assume all the measured charged particles to be pions. For each event, scaled correlation analysis generates a set of measurements, namely, rank-q Renyi entropy (Sq,q = 1,2,3), scale derivative of entropy: information (Iq) and scale derivative of information: dimension lowering ( ∆dq ) as a function of scale for scales in the range [-4,2]. For the calculation of Iq and, ∆dq , we generate for each run (i.e., a set of ~10k events), a reference entropy Sq as a function of scale by taking the ensemble average of the measured event entropy for that run period. Subsequently, an event space is formed in polar coordinates (r,θ , φ ) by taking ∆dq at three different scale points (x,y,z). Each point in this space represents an event. r is a measure of deviation from the reference and (θ , φ ) represents the “shape” of this deviation. In order to look for events with highly correlated structure, we compare real data to a set of Poisson events with constant mT slope and track multiplicity distribution identical to that from real data. Comparison of Poisson and real data yields a set of anomalous events in real data at the 1 part per mil level. We find that the distribution from this special class of events exhibit a higher yield in the range 0.6 ≤ mT ≤ 1GeV / c 2 , compared to normal events from real data as well as Poisson events. The fact that we have used MTPC-only data limits our ability to select high quality tracks; hence the physics interpretation of the observed signal remains ambiguous. Reconstruction of global tracks that combines information from all four TPCs will allow us to use high quality tracks in our analysis. However, we conclude that from the present results that scaled correlation analysis of multiparticle distributions proven to be a powerful tool to search for rare events in NA49 data. 10 Nuclear Physics Laboratory Annual Report, University of Washington (1996) p. 36. 41

  • Page 45

    4.9 STAR trigger interval review J.G. Reid and T.A. Trainor On June 6-7, 1996 the STAR current trigger effort was reviewed and approved intact. Our contribution to this was primarily the recently developed level-2 algorithms using scaled correlation analysis. However, the level-1 autocorrelation techniques that have been in place since 1995 were also discussed. The STAR trigger relies on data coming from five different sub-detector systems. First, there are vertex position detectors which consists of 48 Cherenkov radiators. Second, there is a central trigger barrel consisting of 240 scintillator slats placed cylindrically about the interaction region covering −1 < η < 1 and 0 ≤ φ ≤ 2π . Third, and most important because of its fine segmentation in η , is the multiwire proportional counter which is about 8000 TPC anode wires covering −1 <|η| < 2 . Fourth, two hadron calorimeters operating collectively as a veto calorimeter. Finally, there is an electromagnetic calorimeter (1200 EM modules) covering η from -1 to 1. Our trigger studies rely most heavily on scaled correlation analysis of data from the CTB and MWCs giving a total coverage in η from -2 to 2. The most significant problem encountered in the level-2 trigger was the trigger processor time constraint. The scaled correlation analysis we have developed requires the data set of interest to be sampled using binning systems of a range of sizes (we call this range the scale window). Also, to avoid aliasing effects we need to resample the data set several times at each bin size. Due to this sampling and resampling of the data, our analysis is very CPU time intensive. Since the time constraint for the level-2 trigger in STAR is 5ms we were forced to develop a subset of the full analysis that maintains most of the discriminating power, but is several orders of magnitude faster. To meet these runtime challenges we performed the full analysis on simulated data and identified a few points in the scale window where we expect most of the data's correlations to be. Rather then looking at the entire scale window we concentrate on three main points. We also turned down the resampling at each scale level, allowing in some noise, but not enough to overwhelm the final results. Finally, after some code optimization we were able to meet the timing requirements. In performing the full scale window analysis on simulated data sets we discovered the negative effects of charge integration in the MWCs (as opposed to hit counting). We found that the charge integration method suffers from Landau fluctuations which degrade the trigger quality significantly. Our analysis showed conclusively that to maintain trigger integrity we must implement hit counting in favor of charge integration in spite of the extra expense. 42

  • Page 46

    4.10 DCC simulations J.G. Reid, T.A. Trainor and D.D. Weerasundara Among the proposed phenomena associated with QCD symmetry restoration the disoriented chiral condensate (DCC) serves as a paradigm. If chiral symmetry restoration occurs within the highest energy-density regions of the collision volume the distinction among pion species is predicted to be lost there. In the subsequent cooling process this symmetric state may decay preferentially into a particular isospin state over some finite volume, resulting in nonstatistical deviations of the neutral pion fraction from its usual value (1/3) in some regions of phase space. The observability of the DCC phenomenon hinges on both the thermodynamic trajectory of the collision process and the observability of the soft-pion component of the particle spectrum. Observable DCC effects seem to depend on significant departure from an equilibrium thermodynamic trajectory (during the cooling phase) according to current model calculations. And manifestation of a DCC in particle spectra is expected to occur only for transverse momenta near or below the pion mass. In general, neutral and charged pions are not observed with the same detector components, there are significant low-momentum detection thresholds, and phase-space acceptance overlap for these two species may be incomplete, giving rise to the possibility of significant systematic error. Nevertheless, theoretical and experimental studies of this phenomenon are being vigorously pursued. An initial event-by-event DCC analysis (CERN experiment WA98) utilizes the total yields of neutral and charged particles for each event, distributed on a 2D scatter plot. This distribution shows a significant linear correlation corresponding to correlated variations of the two species with total event multiplicity. The event-by- event signal of interest is the degree of deviation from the correlation axis. This deviation is equivalent to a variation in the Nγ / N ch ratio, or the neutral particle fraction, as a global event variable. Excessive fluctuation (beyond counting statistics) of this deviation from the correlation axis could be an indication of an anomalous neutral-to-charged pion ratio resulting from DCC formation. No statistically significant DCC signal has been observed in a preliminary analysis. In order to evaluate the optimum analysis procedure for a DCC search, and indeed to determine whether DCCs are observable at all with the available experimental apparatus and conditions, we have undertaken a simulation program using scaled correlation analysis as the diagnostic measure. Fig. 4.10-1 summarizes the preliminary results of this study for the STAR detector at RHIC. While about 5000 particles fall into the acceptance of STAR, of these only about 1400 pions lie in the “soft-pion” region of the p t spectrum (and about 700 neutral pions detectable only by the electromagnetic calorimeter). We have examined the observability of several DCC scenarios having different scaling behavior. The outcome of this study will be a characterization of the sensitivity limits of the STAR detector in terms of DCC scaling structure and amplitude. 43

  • Page 47

    STAR DCC simulation 2100 soft pions - 1400 charged pions into STAR acceptance 100 1 100 1 m r m r 90 0.9 90 0.9 ✠ 80 0.8 80 0.8 ✟ 70 0.7 70 0.7 60 0.6 60 0.6 ✞ 50 0.5 50 0.5 40 0.4 40 0.4 30 0.3 30 0.3 ✝ 20 0.2 20 0.2 10 0.1 10 0.1 0 0 0 0 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 x x x x 12 0.1 12 0.1 number ∆d4 number ∆ d4 0.08 0.08 ✆ 10 10 0.06 0.06 0.04 0.04 8 8 0.02 0.02 6 0 6 0 -0.02 -0.02 4 -0.04 4 -0.04 -0.06 -0.06 2 2 -0.08 -0.08 0 -0.1 0 -0.1 0 0.2 0.4 0.6 0.8 1 -2 -1.5 -1 -0.5 0 0.5 0 0.2 0.4 0.6 0.8 1 -2 -1.5 -1 -0.5 0 0.5 r ☎ charged log(e/L) r ☎ charged log(e/L) no DCC large-scale DCC 100 1 100 1 m r m r 90 0.9 90 0.9 80 0.8 80 0.8 70 0.7 70 0.7 60 0.6 60 0.6 50 0.5 50 0.5 ✌ 40 0.4 40 0.4 ☞ 30 0.3 30 0.3 20 0.2 20 0.2 10 0.1 10 0.1 0 0 0 0 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 ✡ x ✡ x ✡ x x 12 ☛ 0.1 12 ☛ 0.1 number ∆d4 number ∆d4 0.08 0.08 10 10 0.06 0.06 0.04 0.04 8 8 0.02 0.02 6 0 6 0 -0.02 -0.02 4 4 -0.04 -0.04 -0.06 -0.06 2 2 -0.08 -0.08 0 -0.1 0 -0.1 0 0.2 0.4 0.6 0.8 1 -2 -1.5 -1 -0.5 0 0.5 0 0.2 0.4 0.6 0.8 1 -2 -1.5 -1 -0.5 0 0.5 r charged log(e/L) r charged log(e/L) small-scale DCC no DCC - 10 events Fig. 4.10-1. DCC simulations showing charged and neutral pion single-event multiplicity m distributions on an arbitrary parameter x with no DCC structure, large-scale structure and small-scale structure. In each case the neutral pion fraction r is also plotted on x. The frequency distribution on r is plotted, and the corresponding dimension transport ∆d(e ) is shown. At bottom right is an ensemble of 10 events without DCC structure to provide a measure of statistical fluctuations. 44

  • Page 48

    4.11 Human face recognition by scaled correlation analysis Q.R. Ahmad, J.G. Reid and T.A. Trainor We have been developing a very general analysis system for use in the STAR trigger as well as for data analysis in STAR and NA49. Since the analysis is a general method of identifying relative correlations in data sets it can be easily applied to the problems of image analysis and comparisons. One exciting application we have begun to explore is human face recognition. However, since we are using a scaled analysis we can also probe the more difficult problem of identifying smaller scale correlations in the faces in our data set, namely looking for correlations between the similar expressions (i.e., smiling, frowning) as made by different people. This problem is also being approached using neural network methods so we can compare our results to a standard approach, and try to better understand the true nature of our analysis and neural network methods. The data we used in this preliminary study was a set of 24 images taken by CCD camera. There are pictures of five different people in this data set, each person making five different expressions (except one person who only makes four). Our goal in analyzing this data was two-fold. First, we wanted to be able to distinguish between the people in the images, so if we were given another photo of a member of the ensemble we could identify which one it was. Second, we examined the small scale region to see what correlations were apparent. We were very successful in distinguishing between the members of the ensemble. After analyzing the data and forming a phase space we found the data for each individual to lie clustered in a region of the space. Identifying individual's expressions as being similar was inconclusive and we have decided to postpone this analysis until we have a larger data set to work with. Unfortunately, the neural network analysis is still incomplete, so we have yet to do the analysis system comparison. The face recognition analysis could also benefit from a larger data set, and in the future we would like to get a significant amount of data (~1000 images of ~50 people). A system of masking off everything in the image but the person's face would also be beneficial. Some measures were taken to do this in the existing data set, cropping the photos above the neck and having the subjects wear shower caps to make their hair look more uniform, but much more could be done here. Our results from this study were very exciting and we plan to pursue it further in the future. 45

  • Page 49

    4.12 NA49 distortion corrections and system geometry T.A. Trainor, D.D. Weerasundara, and the NA49 Collaboration* System geometry of the NA49 TPC system consists of: 1) an idealized mechanical representation of the system (geometry data base), 2) corrections to this idealization required by mechanical, electronic and other functional departures from the ideal (distortion correction data base and parametrizations), 3) determination of the effective time dependent drift speed,11 4) determination of the large scale alignment of the system with a multitarget tracking procedure12 and 5) confirmation of the accuracy and consistency of the system geometry with global tracking across the enter TPC system. The geometry data base contains an idealized mechanical representation of the TPC system. To complete the conversion from the raw data to the space points in the survey system, a measured system of distortions, represented parametrically or in supplementary data bases, must be removed from the transformed data. Knowledge of these distortions can be inferred from laser tracks, grid pulser calibration data and straight tracks in the TPC system. The distortion system in the horizontal and vertical planes are substantially different. Vertical distortions have three sources: 1) spatial equivalent of time offsets in the pad electronics ranging up to 500 µ m for which a correction data base is obtained with a grid pulser system; 2) mechanical distortion of the top plate and the individual sectors and pad-plane PC boards ranging up to 500 µ m : correction data base obtained from straight track residuals; and 3) large-scale spatial drift speed variations across the TPC volume ranging up to 1 mm near the TPC walls and caused by electric field inhomogeneity arising from top plate deformations: correction data base obtained from tracking residuals, electric field calculations, reconstructed event vertex systematics and laser data. Horizontal distortions have two sources: 1) mechanical misalignments of pad-plane PC boards during gluing operations ranging up to 100 µ m -- correction data base derived from straight track residuals and 2) effective charge transport of up to 1 mm near sector boundaries and electric field inhomgeneities near the top plate - parametric correction derived from reconstructed event vertex systematics. Assuming that the internal representation space of for each subdetector is made distortion free to an acceptable degree, a multitarget alignment procedure13 is performed to determine the positions and the orientations of various subdetectors with respect to each other and to an absolute reference system. The final system geometry has been completely determined and the NA49 experiment has begun the large-scale DST production from 1995 experimental run data. * CERN, Geneva, Switzerland. 11 Nuclear Physics Laboratory Annual Report, University of Washington (1996) p. 42. 12 Nuclear Physics Laboratory Annual Report, University of Washington (1996) p. 43. 13 Nuclear Physics Laboratory Annual Report, University of Washington (1996) p. 43. 46

  • Page 50

    4.13 Large volume NA49 data analysis at UW-NPL C. Reynolds, R.J. Seymour, T.A. Trainor and D.D. Weerasundara URHI group recently upgraded its offline computer system with the addition of three new HP 9000/780 model C180 unix boxes, one 5 slot 2500-XT DLT tape robot, one single slot 2000-XT DLT tape drive and three 9 gigabyte Seagate ST19171WD SCSI disks. Fig. 4.13-1 shows a schematic diagram of the hardware configuration of the new system. These three machines use NFS to share 12 gigabytes of internal disks and the 27 gigabyte logical disk giving a total disk capacity of 39 gigabytes in addition to the existing ~18 gigabyte disk capacity. The three HP machines provide cpu power (3x17.2SPECfp95 and 3x10.8SPECint95) comparable to that of a ~10x200MHz Pentium Pro processor farm. The HP C180, being a 64-bit machine, can handle very large (~100 gigabytes) of disk volumes, the three 9 gigabyte disks are mounted as a single 27 gigabyte logical disk volume which enables us to handle files with large data volumes. This is another advantage over Pentium Pro (32-bit) machines which cannot handle disk volumes with tens of gigabytes. NA49 higher level offline data analysis is one of the several computation-intensive research projects (see Section 4.1) that utilize this system. The NA49 experiment produces ~10 terabytes of raw data (~one million events) per year. A factor of 10 data reduction is achieved at the first level DST production (on the SHIFT computer facility14) yielding one terabyte of reconstructed data. Another factor of 10 in data reduction is achieved by creating miniDSTs which contain selected information relevant to higher level physics analysis. This results in ~100 gigabytes of miniDST data volume that fits onto 10 DLT tapes. We estimate that we could analyze ~1 million miniDST events with the scaled correlation analysis (see Section 4.8) in ~30 days on our local offline computer facility. UW-NPL Offline Production Facility 3x17.2 SPECfp95 = 10 Pentium Pros 200 HP C180 HP C180 HP C180 DLT DLT 5-slot 27 Gbyte 1-slot disk robot drive CERN event spectra DST output selected events Fig. 4.13-1. Schematic diagram of the configuration of this offline production facility. 14 http://wwwinfo.cern.ch/pdp/serv/shift.html. 47

  • View More

Get the full picture and Receive alerts on lawsuits, news articles, publications and more!