Experimental deuteron momentum distributions with reduced final state interactions

This dissertation presents a study of the D( e, e′p)n reaction carried out at the Thomas Jefferson National Accelerator Facility (Jefferson Lab) for a set of fixed values of four-momentum transfer Q 2 = 2.1 and 0.8 (GeV/c)2 and for missing momenta pm ranging from pm = 0.03 to pm = 0.65 GeV/c. The analysis resulted in the determination of absolute D(e,e′ p)n cross sections as a function of the recoiling neutron momentum and it's scattering angle with respect to the momentum transfer [vector] q. The angular distribution was compared to various modern theoretical predictions that also included final state interactions. The data confirmed the theoretical prediction of a strong anisotropy of final state interaction contributions at Q2 of 2.1 (GeV/c)2 while at the lower Q2 value, the anisotropy was much less pronounced. At Q2 of 0.8 (GeV/c)2, theories show a large disagreement with the experimental results. The experimental momentum distribution of the bound proton inside the deuteron has been determined for the first time at a set of fixed neutron recoil angles. The momentum distribution is directly related to the ground state wave function of the deuteron in momentum space. The high momentum part of this wave function plays a crucial role in understanding the short-range part of the nucleon-nucleon force. At Q2 = 2.1 (GeV/c)2, the momentum distribution determined at small neutron recoil angles is much less affected by FSI compared to a recoil angle of 75°. In contrast, at Q2 = 0.8 (GeV/c)2 there seems to be no region with reduced FSI for larger missing momenta. Besides the statistical errors, systematic errors of about 5–6 % were included in the final results in order to account for normalization uncertainties and uncertainties in the determi- nation of kinematic veriables. The measurements were carried out using an electron beam energy of 2.8 and 4.7 GeV with beam currents between 10 to 100 &mgr; A. The scattered electrons and the ejected protons originated from a 15cm long liquid deuterium target, and were detected in conicidence with the two high resolution spectrometers of Hall A at Jefferson Lab.^

Electronic excitations in solids studied using inelastic x-ray scattering

Inelastic x-ray scattering can be used to study the electronic structure of matter. The x rays scattered from the target both induce and carry information on the electronic excitations taking place in the system. These excitations are the manifestations of the electronic structure and the physics governing the many-body system. This work presents results of non-resonant inelastic x-ray scattering experiments on a range of materials including metallic, insulating and semiconducting compounds as well as an organic polymer. The experiments were carried out at the National Synchrotron Light Source, USA and at the European Synchrotron Radiation Facility, France. The momentum transfer dependence of the experimental valence- and core-electron excitation spectra is compared with the results of theoretical first principles computations that incorporate the electron-hole interaction. A recently developed method for analyzing the momentum transfer dependence of core-electron excitation spectra is studied in detail. This method is based on real space multiple scattering calculations and is used to extract the angular symmetry components of the local unoccupied density of final states.

Secondary Perturbation Effects in Keplerian Accretion Disks: Elliptical Instability

Origin of turbulence in cold accretion disks, particularly in 3D, which is expected to be hydrodynamic but not magnetohydrodynamic, is a big puzzle. While the flow must exhibit some turbulence in support of the transfer of mass inward and angular momentum outward, according to the linear perturbation theory it should always be stable. We demonstrate that the 3D secondary disturbance to the primarily perturbed disk which exhibits elliptical vortices into the system solves the problem. This result is essentially applicable to the outer region of accretion disks in active galactic nuclei where the gas is significantly cold and neutral in charge and the magnetic Reynolds number is smaller than 10^4.

Basic equations for suspension flows with interphase mass-transfer

In the case of suspension flows, the rate of interphase momentum transfer M(k) and that of interphase energy transfer E(k), which were expressed as a sum of infinite discontinuities by Ishii, have been reduced to the sum of several terms which have concise physical significance. M(k) is composed of the following terms: (i) the momentum carried by the interphase mass transfer; (ii) the interphase drag force due to the relative motion between phases; (iii) the interphase force produced by the concentration gradient of the dispersed phase in a pressure field. And E(k) is composed of the following four terms, that is, the energy carried by the interphase mass transfer, the work produced by the interphase forces of the second and third parts above, and the heat transfer between phases. It is concluded from the results that (i) the term, (-alpha-k-nabla-p), which is related to the pressure gradient in the momentum equation, can be derived from the basic conservation laws without introducing the "shared-pressure presumption"; (ii) the mean velocity of the action point of the interphase drag is the mean velocity of the interface displacement, upsilonBAR-i. It is approximately equal to the mean velocity of the dispersed phase, upsilonBAR-d. Hence the work terms produced by the drag forces are f(dc) . upsilonBAR-d, and f(cd) . upsilonBAR-d, respectively, with upsilonBAR-i not being replaced by the mean velocity of the continuous phase, upsilonBAR-c; (iii) by analogy, the terms of the momentum transfer due to phase change are upsilonBAR-d-GAMMA-c, and upsilonBAR-d-GAMMA-d, respectively; (iv) since the transformation between explicit heat and latent heat occurs in the process of phase change, the algebraic sum of the heat transfer between phases is not equal to zero. Q(ic) and Q(id) are composed of the explicit heat and latent heat, so that the sum Q(ic) + Q(id)) is equal to zero.

Non-linear theory of a positive-column in a magnetic-field

A nonlinear theory of an intermediate pressure discharge column in a magnetic field is presented. Motion of the neutral gas is considered. The continuity and momentum transfer equations for charged particles and neutral particles are solved by numerical methods. The main result obtained is that the rotating velocities of ionic gas and neutral gas are approximately equal. Bohm's criterion and potential inversion in the presence of neutral gas motion are also discussed.

The calculation of photoelectron angular distributions with jet-like structure from scattering theory

Photoelectron angular distributions produced in above-threshold ionization (ATI) are analysed using a nonperturbative scattering theory. The numerical results are in good qualitative agreement with recent measurements. Our study shows that the origin of the jet-like structure arises from the inherent properties of the ATI process and not from the angular momentum of either the initial or the excited states of the atom.

Photoelectron angular distributions from above threshold ionization of hydrogen atoms in strong laser fields

We apply a scattering theory of nonperturbative quantum electrodynamics to study the photoelectron angular distributions (PADs) of a hydrogen atom irradiated by linearly polarized laser light. The calculated PADs show main lobes and jetlike structure. Previous experimental studies reveal that in a set of above-threshold-ionization peaks when the absorbed-photon number increases by one, the jet number also increases by one. Our study confirms this experimental observation. Our calculations further predict that in some cases three more jets may appear with just one-more-photon absorption. With consideration of laser-frequency change, one less jet may also appear with one-more-photon absorption. The jetlike structure of PADs is due to the maxima of generalized phased Bessel functions, not an indication of the quantum number of photoelectron angular momentum states.

Anomaly in photoelectron angular distributions

The recently observed anomaly in photoelectron angular distributions (PADs), the disappearance of the main lobes of PADs which should be usually in the direction of laser polarization, is reinterpreted as a minimum of generalized Bessel functions in the laser-polarization direction with the theory of nonperturbative quantum electrodynamics. The reinterpretation has no artificial fitting parameters and explains more features of the experimentally observed PADs, in contrast to the existing interpretation in which the anomaly is interpreted as a quantum interference of angular momentum partial waves. Some hierarchy anomalies are predicted for further experimental observations.

Toward Understanding Astrophysical Phenomena

Fast radio bursts (FRBs), a novel type of radio pulse, whose physics is not yet understood at all. Only a handful of FRBs had been detected when we started this project. Taking account of the scant observations, we put physical constraints on FRBs. We excluded proposals of a galactic origin for their extraordinarily high dispersion measures (DM), in particular stellar coronas and HII regions. Therefore our work supports an extragalactic origin for FRBs. We show that the resolved scattering tail of FRB 110220 is unlikely to be due to propagation through the intergalactic plasma. Instead the scattering is probably caused by the interstellar medium in the FRB's host galaxy, and indicates that this burst sits in the central region of that galaxy. Pulse durations of order $\ms$ constrain source sizes of FRBs implying enormous brightness temperatures and thus coherent emission. Electric fields near FRBs at cosmological distances would be so strong that they could accelerate free electrons from rest to relativistic energies in a single wave period. When we worked on FRBs, it was unclear whether they were genuine astronomical signals as distinct from `perytons', clearly terrestrial radio bursts, sharing some common properties with FRBs. Recently, in April 2015, astronomers discovered that perytons were emitted by microwave ovens. Radio chirps similar to FRBs were emitted when their doors opened while they were still heating. Evidence for the astronomical nature of FRBs has strengthened since our paper was published. Some bursts have been found to show linear and circular polarizations and Faraday rotation of the linear polarization has also been detected. I hope to resume working on FRBs in the near future. But after we completed our FRB paper, I decided to pause this project because of the lack of observational constraints.

The pulsar triple system, J0733+1715, has its orbital parameters fitted to high accuracy owing to the precise timing of the central $\ms$ pulsar. The two orbits are highly hierarchical, namely $P_{\mathrm{orb,1}}\ll P_{\mathrm{orb,2}}$, where 1 and 2 label the inner and outer white dwarf (WD) companions respectively. Moreover, their orbital planes almost coincide, providing a unique opportunity to study secular interaction associated purely with eccentricity beyond the solar system. Secular interaction only involves effect averaged over many orbits. Thus each companion can be represented by an elliptical wire with its mass distributed inversely proportional to its local orbital speed. Generally there exists a mutual torque, which vanishes only when their apsidal lines are parallel or anti-parallel. To maintain either mode, the eccentricity ratio, $e_1/e_2$, must be of the proper value, so that both apsidal lines precess together. For J0733+1715, $e_1\ll e_2$ for the parallel mode, while $e_1\gg e_2$ for the anti-parallel one. We show that the former precesses $\sim 10$ times slower than the latter. Currently the system is dominated by the parallel mode. Although only a little anti-parallel mode survives, both eccentricities especially $e_1$ oscillate on $\sim 10^3\yr$ timescale. Detectable changes would occur within $\sim 1\yr$. We demonstrate that the anti-parallel mode gets damped $\sim 10^4$ times faster than its parallel brother by any dissipative process diminishing $e_1$. If it is the tidal damping in the inner WD, we proceed to estimate its tidal quantity parameter ($Q$) to be $\sim 10^6$, which was poorly constrained by observations. However, tidal damping may also happen during the preceding low-mass X-ray binary (LMXB) phase or hydrogen thermal nuclear flashes. But, in both cases, the inner companion fills its Roche lobe and probably suffers mass/angular momentum loss, which might cause $e_1$ to grow rather than decay.

Several pairs of solar system satellites occupy mean motion resonances (MMRs). We divide these into two groups according to their proximity to exact resonance. Proximity is measured by the existence of a separatrix in phase space. MMRs between Io-Europa, Europa-Ganymede and Enceladus-Dione are too distant from exact resonance for a separatrix to appear. A separatrix is present only in the phase spaces of the Mimas-Tethys and Titan-Hyperion MMRs and their resonant arguments are the only ones to exhibit substantial librations. When a separatrix is present, tidal damping of eccentricity or inclination excites overstable librations that can lead to passage through resonance on the damping timescale. However, after investigation, we conclude that the librations in the Mimas-Tethys and Titan-Hyperion MMRs are fossils and do not result from overstability.

Rubble piles are common in the solar system. Monolithic elements touch their neighbors in small localized areas. Voids occupy a significant fraction of the volume. In a fluid-free environment, heat cannot conduct through voids; only radiation can transfer energy across them. We model the effective thermal conductivity of a rubble pile and show that it is proportional the square root of the pressure, $P$, for $P\leq \epsy^3\mu$ where $\epsy$ is the material's yield strain and $\mu$ its shear modulus. Our model provides an excellent fit to the depth dependence of the thermal conductivity in the top $140\,\mathrm{cm}$ of the lunar regolith. It also offers an explanation for the low thermal inertias of rocky asteroids and icy satellites. Lastly, we discuss how rubble piles slow down the cooling of small bodies such as asteroids.

Electromagnetic (EM) follow-up observations of gravitational wave (GW) events will help shed light on the nature of the sources, and more can be learned if the EM follow-ups can start as soon as the GW event becomes observable. In this paper, we propose a computationally efficient time-domain algorithm capable of detecting gravitational waves (GWs) from coalescing binaries of compact objects with nearly zero time delay. In case when the signal is strong enough, our algorithm also has the flexibility to trigger EM observation {\it before} the merger. The key to the efficiency of our algorithm arises from the use of chains of so-called Infinite Impulse Response (IIR) filters, which filter time-series data recursively. Computational cost is further reduced by a template interpolation technique that requires filtering to be done only for a much coarser template bank than otherwise required to sufficiently recover optimal signal-to-noise ratio. Towards future detectors with sensitivity extending to lower frequencies, our algorithm's computational cost is shown to increase rather insignificantly compared to the conventional time-domain correlation method. Moreover, at latencies of less than hundreds to thousands of seconds, this method is expected to be computationally more efficient than the straightforward frequency-domain method.

The vacuum Regge Trajectory in conventional field theory

The effect on the scattering amplitude of the existence of a pole in the angular momentum plane near J = 1 in the channel with the quantum numbers of the vacuum is calculated. This is then compared with a fourth order calculation of the scattering of neutral vector mesons from a fermion pair field in the limit of large momentum transfer. The presence of the third double spectral function in the perturbation amplitude complicates the identification of pole trajectory parameters, and the limitations of previous methods of treating this are discussed. A gauge invariant scheme for extracting the contribution of the vacuum trajectory is presented which gives agreement with unitarity predictions, but further calculations must be done to determine the position and slope of the trajectory at s = 0. The residual portion of the amplitude is compared with the Gribov singularity.

Particle fluid mechanics in shear flows, acoustic waves, and shock waves

Three different categories of flow problems of a fluid containing small particles are being considered here. They are: (i) a fluid containing small, non-reacting particles (Parts I and II); (ii) a fluid containing reacting particles (Parts III and IV); and (iii) a fluid containing particles of two distinct sizes with collisions between two groups of particles (Part V).

Part I

A numerical solution is obtained for a fluid containing small particles flowing over an infinite disc rotating at a constant angular velocity. It is a boundary layer type flow, and the boundary layer thickness for the mixture is estimated. For large Reynolds number, the solution suggests the boundary layer approximation of a fluid-particle mixture by assuming W = W_p. The error introduced is consistent with the Prandtl’s boundary layer approximation. Outside the boundary layer, the flow field has to satisfy the “inviscid equation” in which the viscous stress terms are absent while the drag force between the particle cloud and the fluid is still important. Increase of particle concentration reduces the boundary layer thickness and the amount of mixture being transported outwardly is reduced. A new parameter, β = 1/Ω τ_v, is introduced which is also proportional to μ. The secondary flow of the particle cloud depends very much on β. For small values of β, the particle cloud velocity attains its maximum value on the surface of the disc, and for infinitely large values of β, both the radial and axial particle velocity components vanish on the surface of the disc.

Part II

The “inviscid” equation for a gas-particle mixture is linearized to describe the flow over a wavy wall. Corresponding to the Prandtl-Glauert equation for pure gas, a fourth order partial differential equation in terms of the velocity potential ϕ is obtained for the mixture. The solution is obtained for the flow over a periodic wavy wall. For equilibrium flows where λ_v and λ_T approach zero and frozen flows in which λ_v and λ_T become infinitely large, the flow problem is basically similar to that obtained by Ackeret for a pure gas. For finite values of λ_v and λ_T, all quantities except v are not in phase with the wavy wall. Thus the drag coefficient C_D is present even in the subsonic case, and similarly, all quantities decay exponentially for supersonic flows. The phase shift and the attenuation factor increase for increasing particle concentration.

Part III

Using the boundary layer approximation, the initial development of the combustion zone between the laminar mixing of two parallel streams of oxidizing agent and small, solid, combustible particles suspended in an inert gas is investigated. For the special case when the two streams are moving at the same speed, a Green’s function exists for the differential equations describing first order gas temperature and oxidizer concentration. Solutions in terms of error functions and exponential integrals are obtained. Reactions occur within a relatively thin region of the order of λ_D. Thus, it seems advantageous in the general study of two-dimensional laminar flame problems to introduce a chemical boundary layer of thickness λ_D within which reactions take place. Outside this chemical boundary layer, the flow field corresponds to the ordinary fluid dynamics without chemical reaction.

Part IV

The shock wave structure in a condensing medium of small liquid droplets suspended in a homogeneous gas-vapor mixture consists of the conventional compressive wave followed by a relaxation region in which the particle cloud and gas mixture attain momentum and thermal equilibrium. Immediately following the compressive wave, the partial pressure corresponding to the vapor concentration in the gas mixture is higher than the vapor pressure of the liquid droplets and condensation sets in. Farther downstream of the shock, evaporation appears when the particle temperature is raised by the hot surrounding gas mixture. The thickness of the condensation region depends very much on the latent heat. For relatively high latent heat, the condensation zone is small compared with Ʌ_D.

For solid particles suspended initially in an inert gas, the relaxation zone immediately following the compression wave consists of a region where the particle temperature is first being raised to its melting point. When the particles are totally melted as the particle temperature is further increased, evaporation of the particles also plays a role.

The equilibrium condition downstream of the shock can be calculated and is independent of the model of the particle-gas mixture interaction.

Part V

For a gas containing particles of two distinct sizes and satisfying certain conditions, momentum transfer due to collisions between the two groups of particles can be taken into consideration using the classical elastic spherical ball model. Both in the relatively simple problem of normal shock wave and the perturbation solutions for the nozzle flow, the transfer of momentum due to collisions which decreases the velocity difference between the two groups of particles is clearly demonstrated. The difference in temperature as compared with the collisionless case is quite negligible.

Part I: Total yield measurements for the ^(21)Ne(α,n)^)(24)Mg reaction. Part II: A study of the ^(12)C(He,p)^(14)N reaction

PART I

The total cross-section for the reaction ²¹Ne(α, n)²⁴Mg has been measured in the energy range 1.49 Mev ≤ E_cm ≤ 2.6 Mev. The cross-section factor, S(O), for this reaction has been determined, by means of an optical model calculation, to be in the range 1.52 x 10¹² mb-Mev to 2.67 x 10¹² mb-Mev, for interaction radii in the range 5.0 fm to 6.6 fm. With S(O) ≈ 2 x 10¹² mb-Mev, the reaction ²¹Ne(α, n)²⁴Mg can produce a large enough neutron flux to be a significant astrophysical source of neutrons.

PART II

The reaction¹²C(³He, p)¹⁴N has been studied over the energy range 12 Mev ≤ E_lab ≤ 18 Mev. Angular distributions of the proton groups leading to the lowest seven levels in ¹⁴N were obtained.

Distorted wave calculations, based on two-nucleon transfer theory, were performed, and were found to be reliable for obtaining the value of the orbital angular momentum transferred. The present work shows that such calculations do not yield unambiguous values for the spectroscopic factors.

Formation of superheavy nuclei in cold fusion reactions

Within the concept of the dinuclear system (DNS), a dynamical model is proposed for describing the formation of superheavy nuclei in complete fusion reactions by incorporating the coupling of the relative motion to the nucleon transfer process. The capture of two heavy colliding nuclei, the formation of the compound nucleus, and the de-excitation process are calculated by using an empirical coupled channel model, solving a master equation numerically and applying statistical theory, respectively. Evaporation residue excitation functions in cold fusion reactions are investigated systematically and compared with available experimental data. Maximal production cross sections of superheavy nuclei in cold fusion reactions with stable neutron-rich projectiles are obtained. Isotopic trends in the production of the superheavy elements Z=110, 112, 114, 116, 118, and 120 are analyzed systematically. Optimal combinations and the corresponding excitation energies are proposed.

Production of Superheavy Nuclei in Massive fusion reactions

A master equation is constructed to treat the nucleon transfer process in heavy ion fusion reactions to form superheavy nucleus. The relative motion concerning the energy, the angular momentum and the fragment deformation relaxations is explicitly treated to couple with the diffusion process. The nucleon transition probabilities, which are derived microscopically, are thus time dependent. The calculated evaporation residue cross-sections for both cold and hot fusion are in good agreement with the known experimental data.

Competing fusion and quasifission reaction mechanisms in the production of superheavy nuclei

Within the framework of a dinuclear system model, a new master equation is constructed and solved, which includes the relative distance of nuclei as a new dynamical variable in addition to the mass asymmetry variable so that the nucleon transfer, which leads to fusion and the evolution of the relative distance, which leads to quasifission (QF) are treated simultaneously in a consistent way. The QF mass yields and evaporation residual cross sections to produce superheavy nuclei are systematically investigated under this framework. The results fit the experimental data well. It is shown that the Kramers formula gives results of QF, which agree with those by our diffusion treatment, only if the QF barrier is high enough. Otherwise some large discrepancies occur.