Studies of phonon anharmonicity in solids

Today our understanding of the vibrational thermodynamics of materials at low temperatures is emerging nicely, based on the harmonic model in which phonons are independent. At high temperatures, however, this understanding must accommodate how phonons interact with other phonons or with other excitations. We shall see that the phonon-phonon interactions give rise to interesting coupling problems, and essentially modify the equilibrium and non-equilibrium properties of materials, e.g., thermodynamic stability, heat capacity, optical properties and thermal transport of materials. Despite its great importance, to date the anharmonic lattice dynamics is poorly understood and most studies on lattice dynamics still rely on the harmonic or quasiharmonic models. There have been very few studies on the pure phonon anharmonicity and phonon-phonon interactions. The work presented in this thesis is devoted to the development of experimental and computational methods on this subject.

Modern inelastic scattering techniques with neutrons or photons are ideal for sorting out the anharmonic contribution. Analysis of the experimental data can generate vibrational spectra of the materials, i.e., their phonon densities of states or phonon dispersion relations. We obtained high quality data from laser Raman spectrometer, Fourier transform infrared spectrometer and inelastic neutron spectrometer. With accurate phonon spectra data, we obtained the energy shifts and lifetime broadenings of the interacting phonons, and the vibrational entropies of different materials. The understanding of them then relies on the development of the fundamental theories and the computational methods.

We developed an efficient post-processor for analyzing the anharmonic vibrations from the molecular dynamics (MD) calculations. Currently, most first principles methods are not capable of dealing with strong anharmonicity, because the interactions of phonons are ignored at finite temperatures. Our method adopts the Fourier transformed velocity autocorrelation method to handle the big data of time-dependent atomic velocities from MD calculations, and efficiently reconstructs the phonon DOS and phonon dispersion relations. Our calculations can reproduce the phonon frequency shifts and lifetime broadenings very well at various temperatures.

To understand non-harmonic interactions in a microscopic way, we have developed a numerical fitting method to analyze the decay channels of phonon-phonon interactions. Based on the quantum perturbation theory of many-body interactions, this method is used to calculate the three-phonon and four-phonon kinematics subject to the conservation of energy and momentum, taking into account the weight of phonon couplings. We can assess the strengths of phonon-phonon interactions of different channels and anharmonic orders with the calculated two-phonon DOS. This method, with high computational efficiency, is a promising direction to advance our understandings of non-harmonic lattice dynamics and thermal transport properties.

These experimental techniques and theoretical methods have been successfully performed in the study of anharmonic behaviors of metal oxides, including rutile and cuprite stuctures, and will be discussed in detail in Chapters 4 to 6. For example, for rutile titanium dioxide (TiO₂), we found that the anomalous anharmonic behavior of the B_1g mode can be explained by the volume effects on quasiharmonic force constants, and by the explicit cubic and quartic anharmonicity. For rutile tin dioxide (SnO₂), the broadening of the B_2g mode with temperature showed an unusual concave downwards curvature. This curvature was caused by a change with temperature in the number of down-conversion decay channels, originating with the wide band gap in the phonon dispersions. For silver oxide (Ag₂O), strong anharmonic effects were found for both phonons and for the negative thermal expansion.

The importance of spin for observing gravitational waves from coalescing compact binaries with LIGO and Virgo

General Relativity predicts the existence of gravitational waves, which carry information about the physical and dynamical properties of their source. One of the many promising sources of gravitational waves observable by ground-based instruments, such as in LIGO and Virgo, is the coalescence of two compact objects (neutron star or black hole). Black holes and neutron stars sometimes form binaries with short orbital periods, radiating so strongly in gravitational waves that they coalesce on astrophysically short timescales. General Relativity gives precise predictions for the form of the signal emitted by these systems. The most recent searches for theses events used waveform models that neglected the effects of black hole and neutron star spin. However, real astrophysical compact objects, especially black holes, are expected to have large spins. We demonstrate here a data analysis infrastructure which achieves an improved sensitivity to spinning compact binaries by the inclusion of spin effects in the template waveforms. This infrastructure is designed for scalable, low-latency data analysis, ideal for rapid electromagnetic followup of gravitational wave events.

Impact of laser-accelerated micron-size projectile on dense plasma

The impact of a laser-accelerated micron-size projectile on a dense plasma target is studied using two-dimensional particle-in-cell simulations. The projectile is first accelerated by an ultraintense laser. It then impinges on the dense plasma target and merges with the latter. Part of the kinetic energy of the laser-accelerated ions in the projectile is deposited in the fused target, and an extremely high concentration of plasma ions with a mean kinetic energy needed for fusion reaction is induced. The interaction is thus useful for laser-driven impact fusion and as a compact neutron source.

Studies in nuclear reactor dynamics : I. The accuracy of point kinetics. II. The effect of delayed neutrons on the spectrum of the group-diffusion operator

This thesis is a theoretical work on the space-time dynamic behavior of a nuclear reactor without feedback. Diffusion theory with G-energy groups is used.

In the first part the accuracy of the point kinetics (lumped-parameter description) model is examined. The fundamental approximation of this model is the splitting of the neutron density into a product of a known function of space and an unknown function of time; then the properties of the system can be averaged in space through the use of appropriate weighting functions; as a result a set of ordinary differential equations is obtained for the description of time behavior. It is clear that changes of the shape of the neutron-density distribution due to space-dependent perturbations are neglected. This results to an error in the eigenvalues and it is to this error that bounds are derived. This is done by using the method of weighted residuals to reduce the original eigenvalue problem to that of a real asymmetric matrix. Then Gershgorin-type theorems .are used to find discs in the complex plane in which the eigenvalues are contained. The radii of the discs depend on the perturbation in a simple manner.

In the second part the effect of delayed neutrons on the eigenvalues of the group-diffusion operator is examined. The delayed neutrons cause a shifting of the prompt-neutron eigenvalue s and the appearance of the delayed eigenvalues. Using a simple perturbation method this shifting is calculated and the delayed eigenvalues are predicted with good accuracy.

An investigation of techniques for the measurement and interpretation of cosmic ray isotopic abundances

An instrument, the Caltech High Energy Isotope Spectrometer Telescope (HEIST), has been developed to measure isotopic abundances of cosmic ray nuclei in the charge range 3 ≤ Z ≤ 28 and the energy range between 30 and 800 MeV/nuc by employing an energy loss -- residual energy technique. Measurements of particle trajectories and energy losses are made using a multiwire proportional counter hodoscope and a stack of CsI(TI) crystal scintillators, respectively. A detailed analysis has been made of the mass resolution capabilities of this instrument.

Landau fluctuations set a fundamental limit on the attainable mass resolution, which for this instrument ranges between ~.07 AMU for z~3 and ~.2 AMU for z~2b. Contributions to the mass resolution due to uncertainties in measuring the path-length and energy losses of the detected particles are shown to degrade the overall mass resolution to between ~.1 AMU (z~3) and ~.3 AMU (z~2b).

A formalism, based on the leaky box model of cosmic ray propagation, is developed for obtaining isotopic abundance ratios at the cosmic ray sources from abundances measured in local interstellar space for elements having three or more stable isotopes, one of which is believed to be absent at the cosmic ray sources. This purely secondary isotope is used as a tracer of secondary production during propagation. This technique is illustrated for the isotopes of the elements O, Ne, S, Ar and Ca.

The uncertainties in the derived source ratios due to errors in fragmentation and total inelastic cross sections, in observed spectral shapes, and in measured abundances are evaluated. It is shown that the dominant sources of uncertainty are uncorrelated errors in the fragmentation cross sections and statistical uncertainties in measuring local interstellar abundances.

These results are applied to estimate the extent to which uncertainties must be reduced in order to distinguish between cosmic ray production in a solar-like environment and in various environments with greater neutron enrichments.

A Cerenkov-ΔΕ-Cerenkov detector for high energy cosmic ray isotopes and an accelerator study of ^(40)Ar & ^(56)Fe fragmentation

This thesis has two major parts. The first part of the thesis will describe a high energy cosmic ray detector -- the High Energy Isotope Spectrometer Telescope (HEIST). HEIST is a large area (0.25 m²sr) balloon-borne isotope spectrometer designed to make high-resolution measurements of isotopes in the element range from neon to nickel (10 ≤ Z ≤ 28) at energies of about 2 GeV/nucleon. The instrument consists of a stack of 12 NaI(Tl) scintilla tors, two Cerenkov counters, and two plastic scintillators. Each of the 2-cm thick NaI disks is viewed by six 1.5-inch photomultipliers whose combined outputs measure the energy deposition in that layer. In addition, the six outputs from each disk are compared to determine the position at which incident nuclei traverse each layer to an accuracy of ~2 mm. The Cerenkov counters, which measure particle velocity, are each viewed by twelve 5-inch photomultipliers using light integration boxes.

HEIST-2 determines the mass of individual nuclei by measuring both the change in the Lorentz factor (Δγ) that results from traversing the NaI stack, and the energy loss (ΔΕ) in the stack. Since the total energy of an isotope is given by Ε = γM, the mass M can be determined by M = ΔΕ/Δγ. The instrument is designed to achieve a typical mass resolution of 0.2 amu.

The second part of this thesis presents an experimental measurement of the isotopic composition of the fragments from the breakup of high energy ⁴⁰Ar and ⁵⁶Fe nuclei. Cosmic ray composition studies rely heavily on semi-empirical estimates of the cross-sections for the nuclear fragmentation reactions which alter the composition during propagation through the interstellar medium. Experimentally measured yields of isotopes from the fragmentation of ⁴⁰Ar and ⁵⁶Fe are compared with calculated yields based on semi-empirical cross-section formulae. There are two sets of measurements. The first set of measurements, made at the Lawrence Berkeley Laboratory Bevalac using a beam of 287 MeV/nucleon ⁴⁰Ar incident on a CH₂ target, achieves excellent mass resolution (σ_m ≤ 0.2 amu) for isotopes of Mg through K using a Si(Li) detector telescope. The second set of measurements, also made at the Lawrence Berkeley Laboratory Bevalac, using a beam of 583 MeV/nucleon ⁵⁶FeFe incident on a CH₂ target, resolved Cr, Mn, and Fe fragments with a typical mass resolution of ~ 0.25 amu, through the use of the Heavy Isotope Spectrometer Telescope (HIST) which was later carried into space on ISEE-3 in 1978. The general agreement between calculation and experiment is good, but some significant differences are reported here.

相位恢复算法在量子关联衍射成像中的应用研究

随着研究工作的逐步深入,目前已经利用经典热光源实现了关联衍射成像,使得该技术有望在X射线以及中子衍射成像等方面得到广泛应用。在实验利用非相干光得到物体无透镜傅里叶变换频谱的基础上,采用误差消除与输入输出恢复算法,并结合过采样理论,实现了实验所用物体透射率函数的恢复。分别得到了纯振幅物体的振幅分布函数与纯相位物体的相位分布函数。此外,还讨论了实验所得傅里叶变换频谱的噪声等因素对图像恢复结果的影响。

Metallaboratrane Facilitated E‒H Bond Activation and Hydrogenation Catalysis

The E‒H bond activation chemistry of tris-phosophino-iron and -cobalt metallaboratranes is discussed. The ferraboratrane complex (TPB)Fe(N₂) heterolytically activates H‒H and the C‒H bonds of formaldehyde and arylacetylenes across an Fe‒B bond. In particular, H‒H bond cleavage at (TPB)Fe(N₂) is reversible and affords the iron-hydride-borohydride complex (TPB)(μ‒H)Fe(L)(H) (L = H₂, N₂). (TPB)(μ‒H)Fe(L)(H) and (TPB)Fe(N₂) are competent olefin and arylacetylene hydrogenation catalysts. Stoichiometric studies indicate that the B‒H unit is capable of acting as a hydride shuttle in the hydrogenation of olefin and arylacetylene substrates. The heterolytic cleavage of H₂ by the (TPB)Fe system is distinct from the previously reported (TPB)Co(H₂) complex, where H₂ coordinates as a non-classical H₂ adduct based on X-ray, spectroscopic, and reactivity data. The non-classical H₂ ligand in (TPB)Co(H₂) is confirmed in this work by single crystal neutron diffraction, which unequivocally shows an intact H‒H bond of 0.83 Å in the solid state. The neutron structure also shows that the H₂ ligand is localized at two orientations on cobalt trans to the boron. This localization in the solid state contrasts with the results from ENDOR spectroscopy that show that the H₂ ligand freely rotates about the Co‒H₂ axis in frozen solution. Finally, the (TPB)Fe system, as well as related tris-phosphino-iron complexes that contain a different apical ligand unit (Si, PhB, C, and N) in place of the boron in (TPB)Fe, were studied for CO₂ hydrogenation chemistry. The (TPB)Fe system is not catalytically competent, while the silicon, borate, carbon variants, (SiP^R₃)Fe, (PhBP^iPr₃)Fe, and (CP^iPr₃)Fe, respectively, are catalysts for the hydrogenation of CO₂ to formate and methylformate. The hydricity of the CO₂ reactive species in the silatrane system (SiP^iPr₃)Fe(N₂)(H) has been experimentally estimated.

Exploring the Dynamic Radio Sky: The Search for Slow Transients with the VLA

While synoptic surveys in the optical and at high energies have revealed a rich discovery phase space of slow transients, a similar yield is still awaited in the radio. Majority of the past blind surveys, carried out with radio interferometers, have suffered from a low yield of slow transients, ambiguous transient classifications, and contamination by false positives. The newly-refurbished Karl G. Jansky Array (Jansky VLA) offers wider bandwidths for accurate RFI excision as well as substantially-improved sensitivity and survey speed compared with the old VLA. The Jansky VLA thus eliminates the pitfalls of interferometric transient search by facilitating sensitive, wide-field, and near-real-time radio surveys and enabling a systematic exploration of the dynamic radio sky. This thesis aims at carrying out blind Jansky VLA surveys for characterizing the radio variable and transient sources at frequencies of a few GHz and on timescales between days and years. Through joint radio and optical surveys, the thesis addresses outstanding questions pertaining to the rates of slow radio transients (e.g. radio supernovae, tidal disruption events, binary neutron star mergers, stellar flares, etc.), the false-positive foreground relevant for the radio and optical counterpart searches of gravitational wave sources, and the beaming factor of gamma-ray bursts. The need for rapid processing of the Jansky VLA data and near-real-time radio transient search has enabled the development of state-of-the-art software infrastructure. This thesis has successfully demonstrated the Jansky VLA as a powerful transient search instrument, and it serves as a pathfinder for the transient surveys planned for the SKA-mid pathfinder facilities, viz. ASKAP, MeerKAT, and WSRT/Apertif.

A molecular dynamics study of the microscopic properties of simple dense fluids

The microscopic properties of a two-dimensional model dense fluid of Lennard-Jones disks have been studied using the so-called "molecular dynamics" method. Analyses of the computer-generated simulation data in terms of "conventional" thermodynamic and distribution functions verify the physical validity of the model and the simulation technique.

The radial distribution functions g(r) computed from the simulation data exhibit several subsidiary features rather similar to those appearing in some of the g(r) functions obtained by X-ray and thermal neutron diffraction measurements on real simple liquids. In the case of the model fluid, these "anomalous" features are thought to reflect the existence of two or more alternative configurations for local ordering.

Graphical display techniques have been used extensively to provide some intuitive insight into the various microscopic phenomena occurring in the model. For example, "snapshots" of the instantaneous system configurations for different times show that the "excess" area allotted to the fluid is collected into relatively large, irregular, and surprisingly persistent "holes". Plots of the particle trajectories over intervals of 2.0 to 6.0 x 10^-12 sec indicate that the mechanism for diffusion in the dense model fluid is "cooperative" in nature, and that extensive diffusive migration is generally restricted to groups of particles in the vicinity of a hole.

A quantitative analysis of diffusion in the model fluid shows that the cooperative mechanism is not inconsistent with the statistical predictions of existing theories of singlet, or self-diffusion in liquids. The relative diffusion of proximate particles is, however, found to be retarded by short-range dynamic correlations associated with the cooperative mechanism--a result of some importance from the standpoint of bimolecular reaction kinetics in solution.

A new, semi-empirical treatment for relative diffusion in liquids is developed, and is shown to reproduce the relative diffusion phenomena observed in the model fluid quite accurately. When incorporated into the standard Smoluchowski theory of diffusion-controlled reaction kinetics, the more exact treatment of relative diffusion is found to lower the predicted rate of reaction appreciably.

Finally, an entirely new approach to an understanding of the liquid state is suggested. Our experience in dealing with the simulation data--and especially, graphical displays of the simulation data--has led us to conclude that many of the more frustrating scientific problems involving the liquid state would be simplified considerably, were it possible to describe the microscopic structures characteristic of liquids in a concise and precise manner. To this end, we propose that the development of a formal language of partially-ordered structures be investigated.

Hyperfine interactions and nuclear structure

An automatic experimental apparatus for perturbed angular correlation measurements, capable of incorporating Ge(Li) detectors as well as scintillation counters, has been constructed.

The gamma-gamma perturbed angular correlation technique has been used to measure magnetic dipole moments of several nuclear excited states in the osmium transition region. In addition, the hyperfine magnetic fields, experienced by nuclei of 'impurity' atoms embedded in ferromagnetic host lattices, have been determined for several '4d' and '5d' impurity atoms.

The following magnetic dipole moments were obtained in the osmium transition region μ_2⁺(¹⁹⁰Os) = 0.54 ± 0.06 nm μ_4⁺(¹⁹⁰Os) = 0.88 ± 0.48 nm μ_2⁺(¹⁹²Os) = 0.56 ± 0.08 nm μ_2⁺(¹⁹²Pt) = 0.56 ± 0.06 nm μ_2^+’(¹⁹²Pt) = 0.62 ± 0.14 nm.

These results are discussed in terms of three collective nuclear models; the cranking model, the rotation-vibration model and the pairing-plus-quadrupole model. The measurements are found to be in satisfactory agreement with collective descriptions of low lying nuclear states in this region.

The following hyperfine magnetic fields of 'impurities' in ferromagnetic hosts were determined; H_int(Cd Ni) = - (64.0 ± 0.8)kG H_int(Hg Fe) = - (440 ± 105)kG H_int(Hg Co) = - (370 ± 78)kG H_int(Hg Ni) = - (86 ± 22)kG H_int(Tl Fe) = - (185 ± 70)kG H_int(Tl Co) = - (90 ± 35)kG H_int(Ra Fe) = - (105 ± 20)kG H_int(Ra Co) = - (80 ± 16)kG H_int(Ra Ni) = - (30 ± 10)kG, where in H_int(AB); A is the impurity atom embedded in the host lattice B. No quantitative theory is available for comparison. However, these results are found to obey the general systematics displayed by these fields. Several mechanisms which may be responsible for the appearance of these fields are mentioned.

Finally, a theoretical expression for time-differential perturbed angular correlation measurement, which duplicates experimental conditions is developed and its importance in data analysis is discussed.

A search for fractionally charged particles in cosmic rays and a theoretical interpretation of the results

An array of two spark chambers and six trays of plastic scintillation counters was used to search for unaccompanied fractionally charged particles in cosmic rays near sea level. No acceptable events were found with energy losses by ionization between 0.04 and 0.7 that of unit-charged minimum-ionizing particles. New 90%-confidence upper limits were thereby established for the fluxes of fractionally charged particles in cosmic rays, namely, (1.04 ± 0.07)x10^-10 and (2.03 ± 0.16)x10^-10 cm^-2sr^-1sec^-1 for minimum-ionizing particles with charges 1/3 and 2/3, respectively.

In order to be certain that the spark chambers could have functioned for the low levels of ionization expected from particles with small fractional charges, tests were conducted to estimate the efficiency of the chambers as they had been used in this experiment. These tests showed that the spark-chamber system with the track-selection criteria used might have been over 99% efficient for the entire range of energy losses considered.

Lower limits were then obtained for the mass of a quark by considering the above flux limits and a particular model for the production of quarks in cosmic rays. In this model, which is one involving the multi-peripheral Regge hypothesis, the production cross section and a corresponding mass limit are critically dependent on the Regge trajectory assigned to a quark. If quarks are "elementary'' with a flat trajectory, the mass of a quark can be expected to be at least 6 ± 2 BeV/c². If quarks have a trajectory with unit slope, just as the existing hadrons do, the mass of a quark might be as small as 1.3 ± 0.2 BeV/c². For a trajectory with unit slope and a mass larger than a couple of BeV/c², the production cross section may be so low that quarks might never be observed in nature.

Excitation energies and decay properties of T = 3/2 states in ^(17)O, ^(17)F and ^(21)Na

The two lowest T = 3/2 levels in ²¹Na have been studied in the ¹⁹F(³He, n), ²⁰Ne (p,p) and ²⁰Ne (p,p’) reactions, and their excitation energies, spins, parities and widths have been determined. In a separate investigation, branching ratios were measured for the isospin-nonconserving particle decays of the lowest T = 3/2 levels in ¹⁷O and ¹⁷F to the ground state and first two excited states of ¹⁶O, by studying the ¹⁵N(³He,n) ¹⁷F*(p) ¹⁶O and ¹⁸O(³He, α)¹⁷O*(n) ¹⁶O reactions.

The ¹⁹F(³He,n) ²¹Na reaction was studied at incident energies between 4.2 and 5.9 MeV using a pulsed-beam neutron-time-of-flight spectrometer. Two T = 3/2 levels were identified at excitation energies of 8.99 ± 0.05 MeV (J > ½) and 9.22 ± 0.015 MeV (J ^π = ½⁺, Γ ˂ 40 keV). The spins and parities were determined by a comparison of the measured angular distributions with the results of DWBA calculations.

These two levels were also obsesrved as isospin-forbidden resonances in the ²⁰Ne(p,p) and ²⁰Ne(p,p’) reactions. Excitation energies were measured and spins, parities, and widths were determined from a single level dispersion theory analysis. The following results were obtained:

E_x = 8.973 ± 0.007 MeV, J ^π = 5/2 ⁺ or 3/2⁺, Γ ≤ 1.2 keV,

^Γp_o = 0.1 ± 0.05 keV; E_x = 9.217 ± 0.007 MeV, J^π = ½ ⁺,

Γ = 2.3 ± 0.5 keV, ^Γp_o = 1.1 ± 0.3 keV.

Isospin assignments were made on the basis of excitation energies, spins, parities, and widths.

Branching ratios for the isospin-nonconserving proton decays of the 11.20 MeV, T = 3/2 level in ¹⁷F were measured by the ¹⁵N(³He,n) ¹⁷ F*(p) ¹⁶O reaction to be 0.088 ± 0.016 to the ground state of ¹⁶O and 0.22 ± 0.04 to the unresolved 6.05 and 6.13 MeV levels of ¹⁶O. Branching ratios for the neutron decays of the analogous T = 3/2 level, at 11.08 MeV in ¹⁷O, were measured by the ¹⁶O(³He, α)¹⁷O*(n)¹⁶O reaction to be 0.91 ± 0.15 to the ground state of ¹⁶O and 0.05 ± 0.02 to the unresolved 6.05 and 6.13 MeV states. By comparing the ratios of reduced widths for the mirror decays, the form of the isospin impurity in the T = 3/2 levels is shown to depend on T_z.

Um método SN híbrido direto para cálculos de sistemas combustível-moderador em geometria unidimensional

Desenvolvemos nesta dissertação um método híbrido direto para o cálculo do fator de desvantagem e descrição da distribuição do fluxo de nêutrons em sistemas combustível-moderador. Na modelagem matemática, utilizamos a equação de transporte de Boltzmann independente do tempo, considerando espalhamento linearmente anisotrópico no modelo monoenergético e espalhamento isotrópico no modelo multigrupo, na formulação de ordenadas discretas (SN), em geometria unidimensional. Desenvolvemos nesta dissertação um método híbrido direto para o cálculo do fator de desvantagem e descrição da distribuição do fluxo de nêutrons em sistemas combustível-moderador. Na modelagem matemática, utilizamos a equação de transporte de Boltzmann independente do tempo, considerando espalhamento linearmente anisotrópico no modelo monoenergético e espalhamento isotrópico no modelo multigrupo, na formulação de ordenadas discretas (SN), em geometria unidimensional. Descrevemos uma análise espectral das equações de ordenadas discretas (SN)a um grupo e a dois grupos de energia, onde seguimos uma analogia com o método de Case. Utilizamos, neste método, quadraturas angulares diferentes no combustível (NC) e no moderador (NM), onde em geral assumimos que NC > NM . Condições de continuidade especiais que acoplam os fluxos angulares que emergem do combustível (moderador) e incidem no moderador (combustível), foram utilizadas com base na equivalência entre as equações SN e PN-1, o que caracteriza a propriedade híbrida do modelo proposto. Sendo um método híbrido direto, utilizamos as NC + NM equações lineares e algébricas constituídas pelas (NC + NM)/2 condições de contorno reflexivas e (NC + NM)/2 condições de continuidade para determinarmos as NC + NM constantes. Com essas constantes podemos calcular os valores dos fluxos angulares e dos fluxos escalares em qualquer ponto do domínio. Apresentamos resultados numéricos para ilustrar a eficiência e a precisão do método proposto.

Primary cosmic-ray positrons and negatrons in 1968 at energies between 11 and 204 MeV

The cosmic-ray positron and negatron spectra between 11 and 204 MeV have been measured in a series of 3 high-altitude balloon flights launched from Fort Churchill, Manitoba, on July 16, July 21, and July 29, 1968. The detector system consisted of a magnetic spectrometer utilizing a 1000-gauss permanent magnet, scintillation counters, and a lucite Čerenkov counter.

Launches were timed so that the ascent through the 100 g/cm² level of residual atmosphere occurred after the evening geomagnetic cutoff transition. Data gathered during ascent are used to correct for the contribution of atmospheric secondary electrons to the flux measured at float altitude. All flights floated near 2.4 g/cm².

A pronounced morning intensity increase was observed in each flight. We present daytime positron and negatron data which support the interpretation of the diurnal flux variation as a change in the local geomagnetic cutoff. A large diurnal variation was observed in the count rate of positrons and negatrons with magnetic rigidities less than 11 MV and is evidence that the nighttime cutoff was well below this value.

Using nighttime data we derive extraterrestrial positron and negatron spectra. The positron-to-total-electron ratio which we measure indicates that the interstellar secondary, or collision, source contributes ≾50 percent of the electron flux within this energy interval. By comparing our measured positron spectrum with the positron spectrum calculated for the collision source we derive the absolute solar modulation for positrons in 1968. Assuming negligible energy loss during modulation, we derive the total interstellar electron spectrum as well as the spectrum of directly accelerated, or primary, electrons. We examine the effect of adiabatic deceleration and find that many of the conclusions regarding the interstellar electron spectrum are not significantly altered for an assumed energy loss of up to 50 percent of the original energy.