Painleve singularity analysis applied to charged particle dynamics during reconnection

For a plasma in the collisionless regime, test-particle modelling can lend some insight into the macroscopic behavior of the plasma, e.g conductivity and heating. A common example for which this technique is used is a system with electric and magnetic fields given by B = {dollar}\delta y{dollar}cx x + xcx y + {dollar}\gamma{dollar}cx z and E = {dollar}\epsilon{dollar}cx z, where {dollar}\delta{dollar}, {dollar}\gamma{dollar}, and {dollar}\epsilon{dollar} are constant parameters. This model can be used to model plasma behavior near neutral lines, ({dollar}\gamma{dollar} = 0), as well as current sheets ({dollar}\gamma{dollar} = 0, {dollar}\delta{dollar} = 0). The integrability properties of the particle motion in such fields might affect the plasma\u27s macroscopic behavior, and we have asked the question For what values of {dollar}\delta{dollar}, {dollar}\gamma{dollar}, and {dollar}\epsilon{dollar} is the system integrable? to answer this question, we have employed Painleve singularity analysis, which is an examination of the singularity properties of a test particle\u27s equations of motion in the complex time plane. This analysis has identified two field geometries for which the system\u27s particle dynamics are integrable in terms of the second Painleve transcendent: the circular O-line case and the case of the neutral sheet configuration. These geometries yield particle dynamics that are integrable in the Liouville sense (i.e. there exist the proper number of integrals in involution) in an extended phase space which includes the time as a canonical coordinate, and this property is also true for nonzero {dollar}\gamma{dollar}. The singularity property tests also identified a large, dense set of X-line and O-line field geometries that yield dynamics that may possess the weak Painleve property. In the case of the X-line geometries, this result shows little relevance to the physical nature of the system, but the existence of a dense set of elliptical O-line geometries with this property may be related to the fact that for {dollar}\epsilon{dollar} positive, one can construct asymptotic solutions in the limit {dollar}t \to \infty{dollar}

Ultracold Neutral Plasmas

Ultracold neutral plasmas, formed by photoionizing laser-cooled atoms near the ionization threshold, have electron temperatures in the 1-1000 kelvin range and ion temperatures from tens of millikelvin to a few kelvin. They represent a new frontier in the study of neutral plasmas, which traditionally deals with much hotter systems, but they also blur the boundaries of plasma, atomic, condensed matter, and low temperature physics. Modelling these plasmas challenges computational techniques and theories of non-equilibrium systems, so the field has attracted great interest from the theoretical and computational physics communities. By varying laser intensities and wavelengths it is possible to accurately set the initial plasma density and energy, and charged-particle-detection and optical diagnostics allow precise measurements for comparison with theoretical predictions. Recent experiments using optical probes demonstrated that ions in the plasma equilibrate in a strongly coupled fluid phase. Strongly coupled plasmas, in which the electrical interaction energy between charged particles exceeds the average kinetic energy, reverse the traditional energy hierarchy underlying basic plasma concepts such as Debye screening and hydrodynamics. Equilibration in this regime is of particular interest because it involves the establishment of spatial correlations between particles, and it connects to the physics of the interiors of gas-giant planets and inertial confinement fusion devices.Comment: 89 pages, 54 image

Computational modeling of low-density ultracold plasmas

2017 Summer.Includes bibliographical references.In this dissertation I describe a number of different computational investigations which I have undertaken during my time at Colorado State University. Perhaps the most significant of my accomplishments was the development of a general molecular dynamic model that simulates a wide variety of physical phenomena in ultracold plasmas (UCPs). This model formed the basis of most of the numerical investigations discussed in this thesis. The model utilized the massively parallel architecture of GPUs to achieve significant computing speed increases (up to 2 orders of magnitude) above traditional single core computing. This increased computing power allowed for each particle in an actual UCP experimental system to be explicitly modeled in simulations. By using this model, I was able to undertake a number of theoretical investigations into ultracold plasma systems. Chief among these was our lab's investigation of electron center-of-mass damping, in which the molecular dynamics model was an essential tool in interpreting the results of the experiment. Originally, it was assumed that this damping would solely be a function of electron-ion collisions. However, the model was able to identify an additional collisionless damping mechanism that was determined to be significant in the first iteration of our experiment. To mitigate this collisionless damping, the model was used to find a new parameter range where this mechanism was negligible. In this new parameter range, the model was an integral part in verifying the achievement of a record low measured UCP electron temperature of 1.57 ± 0.28K and a record high electron strong coupling parameter, Γ, of 0.35 ± 0.08. Additionally, the model, along with experimental measurements, was used to verify the breakdown of the standard weak coupling approximation for Coulomb collisions. The general molecular dynamics model was also used in other contexts. These included the modeling of both the formation process of ultracold plasmas and the thermalization of the electron component of an ultracold plasma. Our modeling of UCP formation is still in its infancy, and there is still much outstanding work. However, we have already discovered a previously unreported electron heating mechanism that arises from an external electric field being applied during UCP formation. Thermalization modeling showed that the ion density distribution plays a role in the thermalization of electrons in ultracold plasma, a consideration not typically included in plasma modeling. A Gaussian ion density distribution was shown to lead to a slightly faster electron thermalization rate than an equivalent uniform ion density distribution as a result of collisionless effects. Three distinct phases of UCP electron thermalization during formation were identified. Finally, the dissertation will describe additional computational investigations that preceded the general molecular dynamics model. These include simulations of ultracold plasma ion expansion driven by non-neutrality, as well as an investigation into electron evaporation. To test the effects of non-neutrality on ion expansion, a numerical model was developed that used the King model of the electron to describe the electron distribution for an arbitrary charge imbalance. The model found that increased non-neutrality of the plasma led to the rapid expansion of ions on the plasma exterior, which in turn led to a sharp ion cliff-like spatial structure. Additionally, this rapid expansion led to additional cooling of the electron component of the plasma. The evaporation modeling was used to test the underlying assumptions of previously developed analytical expression for charged particle evaporation. The model used Monte Carlo techniques to simulate the collisions and the evaporation process. The model found that neither of the underlying assumption of the charged particle evaporation expressions held true for typical ultracold plasma parameters and provides a route for computations in spite of the breakdown of these two typical assumptions

Electromagnetic Waves

This book is dedicated to various aspects of electromagnetic wave theory and its applications in science and technology. The covered topics include the fundamental physics of electromagnetic waves, theory of electromagnetic wave propagation and scattering, methods of computational analysis, material characterization, electromagnetic properties of plasma, analysis and applications of periodic structures and waveguide components, and finally, the biological effects and medical applications of electromagnetic fields

MS FT-2-2 7 Orthogonal polynomials and quadrature: Theory, computation, and applications

Quadrature rules find many applications in science and engineering. Their analysis is a classical area of applied mathematics and continues to attract considerable attention. This seminar brings together speakers with expertise in a large variety of quadrature rules. It is the aim of the seminar to provide an overview of recent developments in the analysis of quadrature rules. The computation of error estimates and novel applications also are described