19 69ApJ. . .155. .587R -2 -2 2 Earth andMars,intheprotonalpha-particle fluxes(O’GallagherandSimpson dicated bythemeasurementsofaninterplanetary gradient,betweentheorbitsof path lengthsdeducedfromthesestudiesrange afewgcmtomorethan10 by cosmicraysatlowenergiesisseriouslycomplicated byanunknownresidualsolar cm, andtheymayornotbeenergydependent. Eventhoughthesevaluesareof ma (1958),FosterandMulvey(1963),BadhwarDanielDahanayake, gested, therefore,thattheseisotopesareofsecondaryorigin,producedbynuclearinter- nificantly largerthanthatexpectedfromthestudyofuniversalabundance B (ShapiroandSilverberg1968),thedetermination oftheamountmaterialtraversed the sameorderaspathlengthdeducedfrom observedabundancesofLi,Be,and Sreenivasan (1968),andMeyer,Hagge,McDonald (1968).Themeancosmic-ray Kaplon, andLavakare(1964),Kuzhevskii(1966),Fanco-workers(1966a,b), and calculatedintensitiesofdeuteriumhelium-3coulddeterminethemeanamount modulation. Fichtel andReames(1966),BadhwarKaplon (1966),Biswas,Ramadurai,and of materialtraversedbycosmicraysininterstellarspace. actions ofcosmicrayswiththeinterstellargasandthatcomparisonobserved Biswas elal.1967;DennisFreierandWaddington1968).Ithasbeensug- elements (Fanelal.1966a;Fan,Gloeckler,andSimpson1966b;O’Dellet1966; planetary fieldmodulationatsolarminimumwhichabove600MVisoftheformexp(—ri/Rß)with unique determinationofboththecosmic-raypathlengthandresidualinterplanetaryfieldmodulation The AstrophysicalJournal,Vol.155,February1969 nuclei. Thisdeterminationdoesnotdependonanyassumptionsregardingeitherthesourcespectraor at solarminimummaybemadefromacomparisonofthecalculatedandmeasuredintensitiesthesetwo and helium-3bycosmic-rayinteractionsintheGalaxyiscalculatedconsideringenergy-dependent 7] —0.35±0.15BV,whereRandßarerigidityvelocity,respectively. stood intermsofanenergy-independent,cosmic-raypathlength4±1gcm-andaresidualinter- tensities attheEarth,weconcludethatwithinexperimentaluncertaintiesalldatacanbeunder- cross-sections, interactionkinematics,anddemodulatedcosmic-rayspectra.Theresultingfluxatthe the unmodulatedproton-to-alphaparticleratioofprimarycosmicrays.Theproductiondeuterium modulation. Fromacomparisonofthesecalculationswiththemeasureddeuteriumandhelium-3in- Earth isobtainedbytakingintoaccountleakagefromtheGalaxy,ionizationlosses,nuclearbreakup,and 3 COSMIC-RAY DEUTERIUMANDHELIUM-3OFSECONDARYORIGIN © American Astronomical Society •Provided bytheNASA Astrophysics DataSystem The existenceofsucharesidualmodulation,active evenatsolarminimum,isin- Assuming thatcosmic-raydeuteronsandhelium-3nucleiareofsecondaryorigin,weshowa Previous studiesofdandHeproductionweremadebyHayakawa,Ito,Terashi- The abundanceofdeuteriumandhelium-3inthecosmicrayswasfoundtobesig- * NAS-NASAPostdoctoralResident ResearchAssociate. f Researchsupportedbythe National ScienceFoundationundergrantGP-849, Institute ofGeophysicsandPlanetaryPhysics,UniversityCalifornia,LosAngeles AND THERESIDUALMODULATIONOFCOSMICRAYS NASA, GoddardSpaceFlightCenter,Greenbelt,Maryland Received April19,1968;revisedJuly22,1968 Richard E.LiNGENFELTERf Reuven Ramaty* I. INTRODUCTION ABSTRACT AND 587 19 69ApJ. . .155. .587R 3 3 3 3 3 -1-3 production (RamatyandLingenfelter1968),wesuggestedthatbecauseofthesecondary total residualmodulation,which,shortofdirectmeasurementsoverlargedistancesin kinematics andthesecondaryenergydistributionsofvariousinteractionsproducing calculated bytakingintoaccounttheenergydependenceofcross-sections, cosmic-ray protons,alphaparticles,andmediumnucleiwiththeinterstellargas,are charge ratiosorspectralshapesoftheprimarycosmicrayseitherininterstellarspace the relativemodulationoveradistanceofroughlyanastronomicalunitandnot lating theresultantsecondaryequilibriumfluxeswithsamemodulatingfunctionas served primarycosmic-rayfluxesattheEarth,thenbyusingthesedemodulatedintensi- deuterons andhelium-3nuclei.Usingthesespectra,weevaluatetheequilibriumd modulation, thistechniquereliessolelyontheassumptionthatobserveddandHe study ofdandHeinthecosmicraysallowsanessentiallyindependentdetermination origin ofthesenucleiandbecausetheirdifferentcharge-to-massratio,acomparative the solarsystem,remainsessentiallyundetermined. 588 tion intheinterplanetarymagneticfield.Unlikepreviousattemptstoestimatetotal of boththeamountmaterialtraversedbycosmicraysandtheirtotalresidualmodula- length andaresidualsolarmodulationwhichisbothvelocityrigiditydependent. fluxes attheEarthcanbebestunderstoodintermsofanenergy-independentpath used intheinitialdemodulationprocess.BycomparingresultantdandHespectra up. Theresidualsolarmodulationistakenintoaccountbyfirstdemodulatingtheob- He fluxesbyconsideringleakagefromtheGalaxy,ionizationlosses,andnuclearbreak- with theobservations,itisdemonstratedthatmeasureddeuteronandhelium-3 and conclusions. at theirsources.Inthepresentpaperwegiveadetaileddiscussionofthesecalculations are ofsecondaryoriginanddoesnotdependonanyarbitraryassumptionaboutthe 1967; O’Gallagher1967).Themeasurementofthisgradient,however,determinesonly ties tocalculatetheproductionofdandHeininterstellarspace,finallybymodu- helium-3, andtritium,whichdecaysinto arelistedinTable1.Eachofthese acting cosmic-raynuclei;nisthenumberoftarget nucleipergramofinterstellarma- stellar materialmaybewrittenas per nucleonindEaroundE. particles, respectively;jistheequilibriumintensity ininterstellarspaceoftheinter- nuclei intheratioIrlO:!)(SuessandUrey1956; Cameron1959). interactions mayalsoincludethegenerationofpions. cleus, s.Theinteractionswhichwehaveconsidered fortheproductionofdeuterium, cosmic-ray nucleusofenergypernucleonEwillproduce asecondarynucleusofenergy stellar material;Eandaretheenergiespernucleonofprimarysecondary where qistheproductionspectrumofsecondarynucleipersecondgraminter- in theevaluationofequation (1). production cross-sections, andthekinematicssecondaryenergyspectra tobeused terial; cistheinteractioncross-section;andF(E,E )dE istheprobabilitythataprimary 8s 8 s s8 © American Astronomical Society •Provided bytheNASA Astrophysics DataSystem The productionspectraofdeuteriumandhelium-3,resultingfromtheinteraction In apreliminaryreportofratherextensivecalculationsdeuteriumandhelium-3 The rateofproductionsecondarynucleibycosmic-rayinteractionsintheinter- The productionissummedoverallinteractions, i, whichproducethesecondarynu- The interstellarmaterialisassumedtobecomposed ofhydrogen,helium,andCNO We shallnowprovidedetailed discussionsoftheprimarycosmic-ray intensities, REUVEN RAMATYANDRICHARDE.LINGENFELTER (£) =^fdEjmm^FiiEß,),(1) ?ss II. DEUTERONANDHELIUM-3PRODUCTION i Vol. 155 19 69ApJ. . .155. .587R ured inconsiderabledetail fromthresholdtoafewBeVbySchulz(1952), Crawford values obtainedfromthetemporalandspatialvariations mentionedabove,wehave parameter R.Sincethereisnogeneralagreement astothevalueofRo,andsince With thisfunctionalformfor/(R),themodulatingfunctionmentionedabovecanbe modulating parameterswhichdefinethetotalmodulation maybedifferentfromthe explained bysuchamodulationbuthavededuced differentvaluesforthemodulating with RoessentiallytimedependentandoftheorderseveralhundredMVin1965. treated bothr¡andRoasfreeparameterstobedetermined fromthepresentstudy. lagher andSimpson(1967)O’Gallagher havealsoshownthatthemeasured Jokipii (1966,1968)andColeman(1968)demonstratedthatsuchaformis irregularities. ItwasfirstsuggestedbyDorman(1963)that/(R)couldbeapproximated temporal andspatialvariationsofcosmicraysin theinterplanetarymediumcanbe rigidity dependingonthedistributionofmagneticirregularities.Subsequently, by f(R)=RoforR
O^i OLLO
INCIDENT KINETIC ENERGY (MEV/NUCLEON) Fig. 5.—Ranges of the incident kinetic energies per nucleon which can contribute to a given secondary energy for the indicated deuteron-producing reactions. These ranges are determined from the require- ments that cos 0* < 1 and that the deuteron Lorentz factor in the center of mass be less than or equal to its maximum value determined from the conservation of energy and momentum.
Fig. 6.—Angular distribution in the center of mass of deuterons from the reaction He4(/>,d)He3 for various incident proton energies.
© American Astronomical Society • Provided by the NASA Astrophysics Data System 19 69ApJ. . .155. .587R 4 34 A3 1 1 4 + -1 3+ 43 Lorentz-factor distributionsinthecenterofmassarepeakedbackwarddirection exp (—Fi/Fo),whereEisthetritonkineticenergypernucleonand1/Foequalto uncertainties, thetritonspectrumfrommultibodyptbreakupofHeisindependent equation (18),whereyand*aredeterminedbyequations(5)(8)withMi= exponential distributioninthecenterofmassangle,0*.Weassumed,therefore,thatfor of theincidentneutronenergyandcanbefittedbyasimpleexponentialform: shown inFigure7,normalizedtoaunitintegral.Ascanbeseen,withintheexperimental in thereactionsH^a^HeandHe(/>,He)d. F(7,7) forbothdeuteronsandhelium-3nucleiisdeterminedfromequations(9) where cosdodependsontheincidentLorentzfactor7.Inënergyregionof the multibodyreactions,He(p,p')npd,dd,ptandwHe,probabilitydistribution mass. Assumingthatthecharacteristicenergy,Fo,isindependentofincident Innés (1957)forincidentneutronsof90and300MeV,respectively.Thesedataare in proton-heliumreactionsareejectedasfromthecollisionoffreenucleons.Because In termsofthecharacteristic energy,E,definedabove,cos0isgivenby laboratory distributionfunctionF(y,y)becomes the exponentialdistributioninkineticenergypernucleonmaybetransformedintoan and aroundLorentzfactorswhichareaboutequaltothefactorofcenter tial with1/Fo=0.063(MeVnucleon“)“. 0.27 (MeVnucleon“)“.Similarly,thedeuteronspectrummaybefittedbyanexponen- and oftritonsfromptbreakupHeweremeasuredbyTannenwald(1953) to thereactionH(ÿ,d)7r,andthereforedistributionfunctionsF{y,y^)aregivenby and weassumedthisenergydependenceforallvaluesof7considered. 30-90 MeVnucleon^,cos0ocanbewellfittedby to beproducedpredominantlywithlowenergies.Thiswillhappeniftheangularand 2. IthasbeenshownbyRiddifordandWilliams(1960),however,thatpionsproduced and dtTT*alsobecomesignificant,ascanbeseenfromthecross-sectionsshowninFigure where inpainteractions,theparametereequals —1 and(7/)=y*for< function F(7,7s*,cos0*)isgivenby of mass. M =MeMd,andMr. this, wehaveassumedthat,kinematically,thedHe7r°anddt7rreactionsaresimilar 598 REUVENRAMATYANDRICHARDE.LINGENFELTERVol.155 (19). Theparametereis+1inthereactionsHe(ÿ,d)HeandH^ajHe)^—1 t cm S 0O s mccz (7s*) =7cfor>7c«,andinapinteractions,e =+1,andthemeanLorentzfactor tpyv (7s*) isthesameasinpainteractionsfor totalavailableenergyinthecenter 4 3 © American Astronomical Society •Provided bytheNASA Astrophysics DataSystem These spectraindicatethatinthelaboratoryframemultinucleonfragmentstend The energyspectrainthelaboratoryframeofdeuteronsfrompndbreakupHe Above thepion-productionthresholdmorecomplexreactionsleadingtoJHe7r° Since thisreactionresultsintwobodiesthefinalstate,distributionfunction Using equations(4)and(9)withbothy*y* replacedby(73*),wefindthatthe m F(y,7s) =exp[- 2 P(y,7s*, cos0*) X {cos0o[l—exp(—2/cos0o)][(7cD((7/) ~DP}• cos 0o 2 cos 0o[(7c-1)((7s*YDP cos 0o=exp[9.38(7—1)]>(20) l7s}) exp [—(1—ecos0*)/cos0],*/*\\ ^ Eo cos 0o[l-exp(-2/cos0)]^’ O O 2 ^[(7c -D«7.*>DP ' «(t. —7c(7/» COS 0oJ/ 1/ (21) (23) (22) 19 69ApJ. . .155. .587R -1 _1 3 3 4 4 3 be fittedbyanexponential characterizedby1/E=0.14(MeVnucleon) . Assuming by 190MeVprotonsoncarbonareshowninFigure 8,normalizedforunitintegral.As by equations(22)and(23). energy, thed,t,andHe distributions for^CNOreactionsarethesameas thosegiven can beseen,thesespectrahavethesameformas those showninFigure7,andcanalso that thetritondistribution alsohasthisformandthatEisindependent ofincident d, t,andHeproducedinaareactionswithwithout pionproduction: breakup ofHe. equations (4)-(9),withMi=m,we derivethedistributionfunctionfor 0 0 t8r sumed arbitrarilythatthebreakupproductsofeachHenucleuseffectivelycluster ters” havinganisotropicangulardistributioninthecenter-of-masssystem.Thus,from together allowingustotreattheinteractionasatwo-bodyprocesswithtwo“clus- the multibodyreactionsnpd,dd,ptand/¿Hewhenpionsareproduced. Since Eoisindependentoftheincidentenergy,energydependencecos0ode- No. 2,1969COSMIC-RAYDEUTERIUMANDHELIUM-3599 termined solelyby7andy. C8 3 © American Astronomical Society •Provided bytheNASA Astrophysics DataSystem The measurementsofBailey(1956)thespectra of deuteronsandhelium-3produced Fig. 7.—Energydistributioninthelaboratoryframeofdeuterons andtritonsfromthemultibody There arenomeasurementsofd,t,orHeproductionbyaareactions.Wehaveas- We havealsoassumedthesamedistributionfunction,givenbyequation(22),for 0 10203040506070 DEUTERON KINETICENERGY,E.(Mev/Nucleon) -1 F(y,y) =(t-l)• (24) s d 19 69ApJ. . .155. .587R 43 -1 -1 -1 + reactions, anuncertainty ofthisorderisintroducedinthecalculationfrom thissource because ofthelackexperimental dataonthecross-sectionsandkinematics ofthese probability bythemultibodybreakupreactions of cosmic-rayhelium-4onhydrogen. alone. It shouldalsobenotedthattheaareactions,under theassumptionsmadeabove,con- tribute approximately20percentofthetotaldeuteron productionatallenergies,and, H(He, d)Heandbycosmic-rayprotonsinthemultibody breakupofhelium;deuterons with equalprobabilitybycosmic-rayhelium-4 on hydrogeninthepickupreaction helium-3-production spectra,respectively,forrj= 350MVandRo=0.Ascanbeseen in theenergyrangefromabout70to200MeVnucleon areproducedmainlyinthe from Figure9,deuteronsofenergieslessthanabout 70MeVnucleonareproduced of energygreaterthanabout200MeVnucleon areproducedwithroughlyequal carbon interactions. (§>§) resonancereactionofcosmic-rayprotonson hydrogen,H(ÿ,d)7r;anddeuterons using thedemodulatedcosmic-rayspectra,productioncross-sections,andsecondary with increasingt?butisquiteinsensitivetovariationsinR.Therelativecontribution energy distributionsdiscussedabove.Wehaveperformedthesecalculationsforavariety nuclei asafunctionofsecondarykineticenergypernucleoncannowbecalculatedby modulating parameters.Forillustration,weshowinFigures9and10thedeuteron- of modulatingparametersrjandi?o,ingeneraltheabsoluteproductionincreases of thevariousproductionreactions,however,isnotstronglydependentonassumed 600 0 © American Astronomical Society •Provided bytheNASA Astrophysics DataSystem Fig. 8.—Energydistributionsinthelaboratoryframeofdeuteronsandhelium-3nucleifromproton- The rateofproductionpergraminterstellarmaterialdeuteronsandhelium-3 REUVEN RAMATYANDRICHARDE.LINGENFELTER Deuteron andHelium-3KineticEnergy,E(Mev/Nucleon) $ d) ProductionSpectra Vol. 155 19 69ApJ. . .155. .587R 43 3 4 -1 where uisthesecondary cosmic-rayequilibriumdensity;qisthesecondary production rate calculatedabove,per gramofinterstellarmaterial;pisthedensity of interstellar reaction H^He,He)d,asarethelow-energydeuterons,whilehelium-3nucleiof account leakagefromtheGalaxy,ionizationlosses, andnuclearbreakupofthesecond- nuclei attheEarth.Theequilibriumdensitiesof thesenucleiintheGalaxycanbeob- Rq ^0. tained bysolvingasteady-statecontinuityequation inenergyspacewhichtakesinto in theGalaxyandthen,takingintoaccountmodulation, thefluxofdandHe calculated above,wemaynowdeterminetheequilibrium densityofthesesecondaries Ro —0. in thereactionH^He,t)2p,tritiumdecayingtohelium-3.Theareactionsalso higher energiesareproducedpredominantlybythebreakupofcosmic-rayhelium-4 nucleon areproducedprincipallybycosmic-rayhelium-4onhydrogeninthepickup produce about20percentofthehelium-3andintroduceacorrespondinguncertainty aries. in Table1.Thesespectrawerecomputedforademodulatedcosmic-rayintensitywithrj=350MVand No. 2,1969 in Table1.Thesespectrawerecomputedforademodulatedcosmic-rayintensitywith77=350MVand in thetotalproduction. 8 s © American Astronomical Society •Provided bytheNASA Astrophysics DataSystem Fig. 10.—Helium-3-productionspectrapergramofinterstellarmaterialforthevariousmodesshown From Figure10wealsoseethathelium-3nucleiofenergieslessthanabout40MeV From theratesofdeuteriumandhelium-3production intheinterstellarmedium, Fig. 9.—Deuteron-productionspectrapergramofinterstellarmaterialforthevariousmodesshown III. DEUTERONANDHELIUM-3FLUXESATTHEEARTH COSMIC-RAY DEUTERIUMANDHELIUM-3 Fig 910 + ^x7 slb°-J-w-®' u .vd[dE“I s 601 (25 19 69ApJ. . .155. .587R -3 t,T 2 -1 \dE/dt\ loss duetoionization;and\d/visthemeanlifetimeagainstnuclearbreakup.Thesolu- material ingcm;tlistheleakagelifetimefromGalaxy;dE/dtrateofenergy was producedatatimetinthepast.Intermsofthisprobability,numberparticles where ristheeffectivelifetimeagainstleakageandnuclearbreakupgivenby satisfy thefollowingrelationship: where theenergiesEand£',timeintervaltenergylossratedE/dthaveto of energyindEaroundEwhichwereproducedthetimeintervaldt/,with per unittime,/>(/),suchthatP{t)dtistheprobabilityanobservedparticle(at/=0) Equation (28)canbeintegratedbyassumingthattherateofparticleproduction energies indE'around£',isgivenby to equation(26)givenabove,canbeobtainedbyintroducingaprobabilitydistribution 602 tion ofthisequationcanbewrittenas diffusion causedbyscatteringoffmagneticirregularities. Theleakagelifetimetx, formly distributedinatrappingvolumeofradius a,fromwhichcosmicraysescapeby e~. Suchaformisvalidforthesuddenlossesresulting fromnuclearbreakup,andit expression, however,dependsonavarietyoffactors,suchasthespatialdistribution pq{E) istimeindependent.BychangingthevariableofintegrationfromttoE',wefind dE/dt duetoionizationcanbereplacedbypvdE/dR, whereR(E)arechargedparticle mean freepath.Forsuchanexponentialdistribution, equation(30)becomesequivalent mentioned above,maythenbewrittenasrz,= a/(2Xz,fl) whereXlisthescattering, might beagoodapproximationtophysicalsituation inwhichthesourcesareuni- properties ofthetrappingvolumeparticles. tribution isreplacedbythemoregeneralexpression,P{t).Thefunctionalformofthis Equation (30)issimilartoequation(26)givenabove,exceptthattheexponentialdis- that theequilibriumdensityu(E)isgivenby ranges givenbyBarkasandBerger(1964).The deuteron-breakup cross-sectionvaries the sources,natureofparticlepropagationfromsourcestoEarth,and with energybutsincewe findthatthelifetimerisdeterminedprincipally byleakage, we havetakenittobea constantof170mbwhichcorrespondstotheaverage value to equation(26). over theenergyrange10-100 MeVnucleon(VanOers1963;Davisonetal. 1963).The 8 © American Astronomical Society •Provided bytheNASA Astrophysics DataSystem A moregeneralexpressionforwhichunderanappropriatesimplificationreduces In thepresentstudywehaveassumedasimple exponentialdistribution,P(t)= We nowconsiderthenumericalevaluationofequation (26).Theenergylossrate REUVEN RAMATYANDRICHARDE.LINGENFELTER “•® ■-Rwi“>>(-/- u{E) = s du{E) =pqE')^P(t)dt,(28) s \dE/dt\ P t = fdE'q(E')p(f a E \ E 4 \dE/dt\‘ dE" r I dE/dt/’ dE" \ Vol. 155 (26) (29) (27) (30) r-p¿¡ coLO LO No. 2, 1969 COSMIC-RAY DEUTERIUM AND HELIUM-3 603
helium-3-breakup cross-section is much smaller yet and hence was neglected. The break- l LOO^i up mean path length, Xd, is given by X¿ = {(TdNap)~ , where Na is Avogadro’s number. OL Introducing a mean path length against leakage, x = pvri, we can rewrite equation (26) as
m = ^ / iE,.(E) exp [- /(i + „N) § dE'] , (31)
where ys(E) is the secondary cosmic-ray equilibrium intensity in interstellar space. We have evaluated equation (31) for a variety of production spectra, qs(-E), corresponding to a range of values of the modulating parameters rj and Ro, defined above. We have
DEUTERON KINETIC ENERGY (MEV/NUCLEON) HELIUM-3 KINETIC ENERGY (MEV/NUCLEON ) Fig 11 Fig 12 Fig. 11.—Observed intensity of deuterons, together with the calculated deuteron intensities at the Earth for various values of interstellar path lengths and modulating parameters. Fig. 12.—Observed intensity of helium-3 nuclei, together with the calculated helium-3 intensities at the Earth for various values of interstellar path lengths and modulating parameters.
then obtained the d and He3 intensities at the Earth by modulating the resultant equilib- rium spectra in interstellar space with the same modulating function as that used in the demodulation process. The results of these calculations, for 7?0 = 0 and for several values of the mean path length x and the modulating parameter rj, are shown in Figures 11 and 12 for deuterons and helium-3 nuclei, respectively. The net effect of the process of demodulation and the subsequent modulation, de- scribed above, is a decrease of the secondary intensity, because secondaries of a given energy are produced by primaries of higher energy which are therefore less modulated than the secondaries which they produce. This can be seen in Figures 11 and 12 from a comparison of curves A and C, for zero modulation and modulation characterized by 7] = 0.35 BV. Also shown in Figures 11 and 12 are measurements at the Earth of deuter- ons and He3 nuclei.
© American Astronomical Society • Provided by the NASA Astrophysics Data System 19 69ApJ. . .155. .587R -2 3 3 3 3 3-1 3 3 -2 -2 3-2 23 z 3 3 -1 -13 -1 fitting curveCandthelimiting curvesBandDinFigures1112. form exp(—r)/Rß)characterized byrj=0.35±0.15BV.Thiscanbeseen bythebest- which bestfittheobservedintensities.Within experimentaluncertainties,wefind an energy-independentpath length,x=4±1gcm,andaresidualmodulation ofthe that theavailabledataon boththedeuteronsandhelium-3nucleicanbe best fittedby and arangeofvaluesxrj,therebydetermining thevaluesoftheseparameters of about600MV,wehavecalculatedthedand He spectraattheEarthforR=0, spectrum oftheresidualsolarmodulation. parison offuturemeasurementsthe¿/Heratio withthecalculationsdescribedin is about1.3,whereasthemeasuredratioatsameenergypernucleonperhapsas velocity aswellrigiditydependent.Becauseof this, wewishtosuggestthatthecom- residual modulationbutalsothatabovearigidity ofabout600MVthismodulationis dependent modulationwouldincreasethecalculated ¿/Heratio.Thus,thecomparison have ahigherrigidityandhencearelessmodulatedthanHenuclei,suchrigidity- calculated ¿/Heratioat,forexample,50MeVneucleonzeroresidualmodulation resented byapowerlawinRwithvaryingspectralindex,rangingfromabout—0.5 planetary fieldvariationsindicatesthatthefunction/(J?),definedabove,canberep- the presentpapercoulddeterminemoreaccurately thechargedependenceandenergy great as6.Suchadifferencecouldresultfrommodulationoftheformexp(—rj/Rß), deuteron andhelium-3data.AscanbeseenfromcurvesAinFigures1112,the at lowrigiditiesto—2highrigidities.Withintheexperimentalerrors,thiswouldbe According toJokipiiandColeman(1968)theobservedpowerspectrumofinter- as discussedabove,therearesignificantdisagreementsamongthevariousexperimenters by theobservedgradientincosmic-rayintensitybetweenorbitsofEarth as givenbyequation(2)forR>Ro.Sinceatthesameenergypernucleon,douterons 600 MV.Sucharigidity-dependentmodulationisindeeddemandedbytheobserved consistent withthemodulatingfunctiongivenbyequation(2)jRooforder of thedandHemeasurementswithcurveAsuggests notonlythenecessityofsome sidual interplanetaryfieldmodulationatsolarminimum.Thisconclusionissupported leakage pathlength.Ontheotherhand,comparisonofcurvesB,C,andDwith as totheenergydependenceofobservedcosmic-rayspatialandtemporalvariations. and Mars(O’GallagherSimpson1967;O’Gallagher1967). d andHedatasuggeststhatthesespectracanbeunderstoodeasilyintermsofare- clude thatthemeasuredspectraareprobablynotresultofanenergy-dependent lated equilibriumisdeterminedprincipallybyionizationlosses.Ascanbeseenfrom nuclei becomemuchsmallerthanthe4gcmleakagepathlength,andthuscalcu- intensities (O’Delle/a/.1966;Biswase/aZ.1967;DennisFreierandWadding- energy, xwouldhavetobeabout0.1gcm.Thisrequiresaverystrongvariationin where therangeofanHenucleusis1gcm,thisdiscrepancyamountstomorethan the calculatedequilibriumfluxisdeterminedprincipallybyleakagefromGalaxy, length ofabout4gcm“andzeromodulation,thecalculateddHeintensitiesat the pathlength(oforderxE)whichseemsratherunlikely.Wethereforecon- and theintensityisdirectlyproportionaltox.AtlowerenergiesrangesofdHe an orderofmagnitude,andthusintoaccountfortheobservedintensityatthis curve Aattheselowerenergies,thecalculateddandHeintensitiesaremuchgreater against ionizationlossesisgreaterthanthevaluesofxinvestigated.Attheseenergies ton 1968).Asensitivedeterminationofthemeanleakagepathlength,x,canbemadeonly sidual modulation.AscanbeseenfromcurveA,wefindthatforacosmic-raypath 604 than themeasuredvalues(Fanetal.1966a,b).Forexample,at50MeVnucleon, at thesehigherenergies(>150MeVnucleon),wheretherangeofdandHenuclei energies greaterthanabout200MeVnucleonareingoodagreementwiththemeasured 0 © American Astronomical Society •Provided bytheNASA Astrophysics DataSystem Since thepresentcalculationsdonotdependon processes occurringbelowrigidities The rigiditydependenceofthisresidualmodulation,however,isnotwellknown,and In comparingthesedatawithourcalculations,wefirstconsiderthecaseofzerore- REUVEN RAMATYANDRICHARDE.LINGENFELTER Vol. 155 19 69ApJ. . .155. .587R -1 -2 3 3 particle spectrumliesabove thiscurve.Ifsignificant,changeofspectral shapemay modulated protonandalpha-particle spectraareconsistentwithacosmic-ray source calculations cannotbemeaningfully applied.Above200MeVnucleon boththede- per nucleoncorrespondingtoarigidityof600MV, belowwhichtheresultsofpresent which isapowerlawin total energy(curveIII)whereasbelowthis thealpha- solar-minimum protonandalpha-particlespectra withthemodulatingfunctiongiven respectively, for#=4gcm.Forcomparison wehavedemodulatedtheobserved from ahomogeneouslydistributedgalacticsource withagivenenergyspectrum.Such are alsoshowninFigure15.Thesedemodulatedcurves extenddowntokineticenergies rigidity, andtotalenergypernucleonareshown in Figure15bycurves/,//,andIII, spectra resultingfromsourceswhicharepower lawsinkineticenergypernucleon, the interstellarequilibriumintensityofprimary protons andalphaparticlesresulting above (77=0.35BV),andtheresultantinterstellar protonandalpha-particlespectra at theEarthforvariousvaluesofinterstellarpathlengths andmodulatingparameters. the Earthforvariousvaluesofinterstellarpathlengthsandmodulatingparameters. rigidity dependentandisoftheformexp(—rj/Rß). ¿/He ratioscanbeexplainedbyanenergy-independent,path-lengthdistributioncoupled However, ascanbeseenfromthesefigures,ourcalculationsshowthatbothHe/Heand and residualinterplanetarymodulationontheprimarycosmic-raysourcespectra.Using with aresidualsolarmodulation,whichaboveabout600MVisbothvelocityand Figures 13and14.Biswasetal.(1968)haverecentlyconcludedthattheenergydepend- No. 2,1969 the cosmic-raypropagationmodelassumedinprecedingcalculations,wecompute ence oftheHe/Heratiocanbeaccountedforonlybyanenergy-dependentpathlength. 3 3 © American Astronomical Society •Provided bytheNASA Astrophysics DataSystem Fig. 14.—Observedratioofhelium-3toheliumnuclei,togetherwiththecalculatedHe/Heratios Fig. 13.—Observedratioofdouteronstoheliumnuclei,togetherwiththecalculated¿/Heratiosat The ¿/HeandHe/Heratiosderivedfromthepresentcalculationsareshownin Finally weshallexaminetheimplicationsofinterstellarcosmic-raypathlength COSMIC-RAY DEUTERIUMANDHELIUM-3 Fig 13Fig.14 605 19 69ApJ. . .155. .587R 2 interactions ininterstellarspace. reflect eitherasimilarchangeinthespectraofcosmic-raysourcesormaybecaused by avarietyofeffectssuchassolarorstellaraccelerationandproductionnuclear 606 energy pernucleon,respectively,forx—cm-. represent theequilibriumspectrawhichwouldresultfromhomogeneouslydistributedcosmic-ray Biswas, S.,Lavakare,P.J.,Ramadurai, S.,andSreenivasan,N.1967,Proc.IndianAcad. Sei.,65,104. Benveniste, J.,andCork,B. 1953, Phys.Reo.,89,422. .19666,/.Geophys.Res.,71,1771. Balasubrahmanyan, V.K.,Hagge,D.E.,Ludwig,G.H., andMcDonald,F.B.1966a,Proc.Internat, Bailey, L.E.1956,ReprintU.C.R.L.-3334,UniversityofCalifornia RadiationLab.,Berkeley. Badhwar, G.D.,andDaniel,R.1963,Progr.Theoret.Phys., 30,615. sources havingenergyspectrawhicharepowerlawsinkineticpernucleon,rigidity,andtotal Barkas, W.F.,andBerger,M.J.1964,NASASP-3013. Baldoni, B.,Focardi,S.,Hromadnik,H.,Monari,L.,Saporetti, F.,Femino,S.,Mezzanares,Berto- Balasubrahmanyan, V.K.,Boldt,E.,andPalmeira,R.A. 1967,/.Geophys.Res.,72,27. Badhwar, G.D.,andKaplon,M.F.1966,/.Geophys.Res., 71,5175. the observedsolar-minimumspectrawithmodulatingfunctionexp(—0.35/jR).Curves/,//,andIII © American Astronomical Society •Provided bytheNASA Astrophysics DataSystem lini, E.,andGialanella,G.1962,NuovoCimento,26,1376. Fig. 15.—^Interstellarequilibriumintensityofprotonsandalphaparticlesobtainedbydemodulating Conf CosmicRays,London,1,427. REUVEN RAMATYANDRICHARDE.LINGENFELTER REFERENCES Vol. 155 19 69ApJ. . .155. .587R Rosenfeld, A.H.1954,Phys.Rev.,96,139. Schwarzschild, A.,andZupancic, C.1963,Phys.Rev.,129,854. Schulz, A.H.1952,U.C.R.L.Rept.No.1756,Universityof CaliforniaRadiationLab,Berkeley. Riddiford, L.,andWilliams,A.W.1960,Proc.Roy.Soc.London, A,257,316. Ramaty, R.,andLingenfelter,R.E.1968,CanadianJ.Phys., 46,S627. Pickup, E.,Robinson,D.K.,andSalant,E.O.1962,Phys. Rev.,125,2091. Ormes, J.F.,andWebber,W.R.1966,Proc.Internat.Conf. CosmicRays,London,1,349. O’Gallagher, J.,andSimpson,J.A.1967,Ap.147,819. O’Gallagher, J.1967,Ap.J.,150,675. Oers, W.T.H.van.1963,unpublishedPh.D.,thesis.University ofAmsterdam. O’Dell, F.W.,Shapiro,N.M.,Silberberg,R.,andStiller, B.1966,Proc.Internat.ConfCosmicRays, Selove, W.,andTeem,J.M. 1958, Phys.Rev.,112,1658. Parker, E.N.1958,Phys.Rev.,110,1445. Neganov, B.S.,andParfenov,L.1958,Zh.Eksp.iTeoret. Fiz.,34,767. Legg, J.C.1963,Phys.Rev.,129,272. Jokipii, J.R.,andColeman,P.1968,Geophys.Res.,73,549. Innés, W.H.1957,U.C.R.L.Rept.No.8040,UniversityofCaliforniaRadiationLab.,Berkeley. Nadi, M.,andRiad,F.1964,Nucl.Phys.,50,33. Moulthrop, P.H.1955,Phys.Rev.,99,1509. Meyer, J.P.,Hagge,D.E.,andMcDonald,F.B.1968,CanadianPhys.,46,S503. Kozadaev, M.S.,Kulyukin,M.,Sulyaev,R.Filippov,A.I.,Shcherbakov,Yu1960,Zh. .1968,CanadianJ.Phys.,46,S950. Jokipii, J.R.1966,Ap.J.,146,480. Fireman, E.L.,andRowland,F.S.1955,Phys.Rev.,97,780. .1968,ibid.,73,4261. Kuzhevskii, B.M.1966,Yader.Fiz.,4,130. Kavaloski, C.D.,Bassani,G.,andHintz,N.M.1963,Phys.Rev.,132,813. Honda, M.,andLal,D.1960,Phys.Rev.,118,1618. Hofmann, D.J.,andWinckler,J.R.1966,Phys.Rev.Letters,16,109. Hayakawa, S.,Ito,K.,andTerashima,Y.1958,Progr.Theoret.Phys.,Suppl.,6,1. Fichtel, C.E.,andReames,D.V.1966,Phys.Rev.,149,995. Hayakawa, S.,Horikawa,N.,Kajikawa,R.,Kikuchi,K.,Kobayakawa,H.,Matsuda,Nagata, Freier, P.S.,andWaddington,C.J.1965,Geophys.Res.,70,5753. Foster, F.,andMulvey,J.H.1963,NuovoCimento,27,93. Fields, T.H.,Fox,J.G.,Kane,A.,Stallwood,R.andSutton,B.1954,Phys.Rev.,95,638. Fickinger, W.J.,Pickup,E.,Robinson,D.K.,andSalant,E.O.1962,Phys.Rev.,125,2082. Fan, C.Y.,Gloeckler,G.,andSimpson,J.A.1966Ô,Phys.Rev.Letters,17,329. Fan, C.Y.,Gloeckler,G.,Hsieh,K.C.,andSimpson,J.A.1966a,Phys.Rev.Letters,16,813. Durbin, B.,Loar,H.,andSteinberger,J.1951,Phys.Rev.,84,581. Dorman, L.1.1963,ProgressinElem.ParticleandCosmic-RayPhys.,Vol.7(Amsterdam:North-Holland Bunch, S.M.,Forster,H.H.,andKim,C.1964,Nucl.Phys,53,241. Brun, C.,Leforth,M.,andTarrago,X.1962,/.Phys.Rad.,23,167. Biswas, S,Ramadurai,S.,andSreenivasan,N1968,CanadianJ.Phys46,S593. Hart, E.L.,Louttit,R.I.,Luers,D.,Morris,T.W.,Willis,W.J,andYamamoto,S.1962,Phys.Rev., Davison, M.,Hopkins,H.W.K.,Lyons,L.,andShaw,D.F.1963,Phys.Letters,3,358. Guzhavin, V.M.,Kliger,G.K.,Kolganov,Z.,Lebedev.A.V.,Marish,K.S,Prokoshkin,Yu.D., Gloeckler, G.,andJokipii,J.R.1966,Phys.Rev.Letters,17,203. Dennis, B.R.,Badhwar,G.D.,Deney,C.L.,andKaplon,M.F.1967,Bull.Am.Phys.Soc.,12,583. Cameron, A.G.W.1959,Ap.J.,129,676. No. 2,1969 Dahanayake, C.,Kaplon,M.F.,andLevakare,P.J.1964,Geophys.Res.,69,3681. Crawford, F.S.,1953,U.C.R.L.Rept.No.1756,UniversityofCaliforniaRadiationLab.,Berkeley. Comstock, G.M.,Fan,C.Y.,andSimpson,J.A.1966,Ap.J.,146,51. Chapman, K.R.,Jones,T.W.,Khan,Q.H.,McKee,J.S.C.,Ray,H.B.vander,andTanimura,Y. Cairns, D.J.,Griffith,T.C.,Lush,GJ,Metheringham,A.andThomas,R.H.1964,Nucl.Phys.,60, Currie, L.A.,Libby,W.F.,andWolfgang,R.1956,Phys.Rev.,101,1557. Currie, L.A.1959,Phys.Rev.,114,878. © American Astronomical Society •Provided bytheNASA Astrophysics DataSystem London, 1,412. Eksp. iTeoret.Fiz.,38,708. 369. and Sumi,Y.1964,Phys.Letters,8,330. Publishing Co.),p.1. Smolyankin, V.T.,Sokolov,A.P.,Soroko,L.M.,andWachuang,T.1964,Zh.Eksp.iTeor.Fiz.46, 1964, Phys.Letters,11,253. 126, 747. 1245. COSMIC-RAY DEUTERIUMANDHELIUM-3 607 19 69ApJ. . .155. .587R 608 Stevenson, M.L.1953,U.C.R.L.Rept.No.2188,UniversityofCaliforniaRadiationLab.,Berkeley. Smith, G.A.,Courant,A.H.,Fowler,E.C.,Kraybill,Sandweiss,J.,andTaft,H.1961,Phys.Rev., Shapiro, M.M.,andSilberberg,R.1968,CanadianJ.Phys,46,S561. Tannenwald, P.0.1953,Phys.Rev.,89,508. Stadler, H.L.1954,Phys.Rev.,96,496. Smith, W.R.,andIvash,E.V.1962,Phys.Rev.,128,1175. Waddington, C.J.,andFreier,P.S.1966,Proc.Internat.Conf.CosmicRays,London,1,345. Suess, H.E.,andUrey,C.1956,Rev.Mod.Phys.,28,53. Wickersham, A.F.1957,Phys.Rev.,107,1050. © American Astronomical Society •Provided bytheNASA Astrophysics DataSystem 123, 2160. REUVEN RAMATYANDRICHARDE.LINGENFELTER © 1969.TheUniversityofChicago.Allrightsreserved.PrintedinU.SA.