Simulation of laser–plasma interactions and fast-electron transport in inhomogeneous plasma

A new framework is introduced for kinetic simulation of laser–plasma interactions in an inhomogeneous plasma motivated by the goal of performing integrated kinetic simulations of fast-ignition laser fusion. The algorithm addresses the propagation and absorption of an intense electromagnetic wave in an ionized plasma leading to the generation and transport of an energetic electron component. The energetic electrons propagate farther into the plasma to much higher densities where Coulomb collisions become important. The high-density plasma supports an energetic electron current, return currents, self-consistent electric fields associated with maintaining quasi-neutrality, and self-consistent magnetic fields due to the currents. Collisions of the electrons and ions are calculated accurately to track the energetic electrons and model their interactions with the background plasma. Up to a density well above critical density, where the laser electromagnetic field is evanescent, Maxwell’s equations are solved with a conventional particle-based, finite-difference scheme. In the higher-density plasma, Maxwell’s equations are solved using an Ohm’s law neglecting the inertia of the background electrons with the option of omitting the displacement current in Ampere’s law. Particle equations of motion with binary collisions are solved for all electrons and ions throughout the system using weighted particles to resolve the density gradient efficiently. The algorithm is analyzed and demonstrated in simulation examples. The simulation scheme introduced here achieves significantly improved efficiencies.     

Stochastic perturbation of Kerr law optical solitons

The soliton perturbation theory is used to study and analyze the stochastic perturbation of optical solitons, due to Kerr law nonlinearity, in addition to deterministic perturbations of optical solitons that is governed by the nonlinear Schrödingers equation. The Langevin equations are derived and analyzed. The deterministic perturbations that are considered here are of both Hamiltonian as well as of non-Hamiltonian type.     

Coupled out of plane vibrations of spiral beams for micro-scale applications

An analytical method is proposed to calculate the natural frequencies and the corresponding mode shape functions of an Archimedean spiral beam. The deflection of the beam is due to both bending and torsion, which makes the problem coupled in nature. The governing partial differential equations and the boundary conditions are derived using Hamilton’s principle. Two factors make the vibrations of spirals different from oscillations of constant radius arcs. The first is the presence of terms with derivatives of the radius in the governing equations of spirals and the second is the fact that variations of radius of the beam causes the coefficients of the differential equations to be variable. It is demonstrated, using perturbation techniques that the derivative of the radius terms have negligible effect on structure’s dynamics. The spiral is then approximated with many merging constant-radius curved sections joined together to approximate the slow change of radius along the spiral. The equations of motion are formulated in non-dimensional form and the effect of all the key parameters on natural frequencies is presented. Non-dimensional curves are used to summarize the results for clarity. We also solve the governing equations using Rayleigh’s approximate method. The fundamental frequency results of the exact and Rayleigh’s method are in close agreement. This to some extent verifies the exact solutions. The results show that the vibration of spirals is mostly torsional which complicates using the spiral beam as a host for a sensor or energy harvesting device.     

Coarsely grained stochastic Boltzmann equation and its moment equations

The stochastic Boltzmann equation is coarsely grained. The coarsely grained stochastic (CGS) Boltzmann equation has fluctuating terms in its collision term. On the basis of the CGS Boltzmann equation, reduced Grad’s 26 moment equations are derived. Coarsely grained moment equations obtained from the CGS Boltzmann equation show that fluctuating terms remain as nonvanishing terms owing to the nonlinearity in the collision term of the CGS Boltzmann equation. The Navier-Stokes-Fourier law obtained using the CGS Boltzmann equation indicates that the pressure deviator and heat flux include fluctuations of their one-order higher moments.     

Benford’s law in non-equilibrium processes: Droplet collisions case

Benford’s law is investigated for the simulation results generated from non-equilibrium molecular dynamics. A statistic to measure how closely a set of the numbers follows Benford’s law is defined. The simulation data are from the collisions of two nano droplets with different impact velocities. When a non-equilibrium system returns to its equilibrium state, some physical quantities relevant to the non-equilibrium settings follow Benford’s law more closely. The initial settings for the non-equilibrium state can be interpreted as a data fabrication of its corresponding equilibrium state. A connection with the Shannon entropy for the first digit distribution is also discussed.     

Towards a physics of Internet traffic in a geographic network

A set of equations from a biased random walk are shown to describe the time-based Gaussian distributions of Internet traffic relative to the Earth’s time zones. The Internet is an example of a more general physical problem dealing with motion near the speed of light relative to different time frames of reference. The second order differential equation (DE) takes the form of ‘time diffusion’ near the speed of light or alternatively considered as a complex variable with real time and imaginary longitudinal components. Congestion waves are generated by peak global traffic from different time zones following the Earth’s revolution. The DE is divided into space and time operators for discussion and each component solution, including constants, is illustrated using data from a global network compiled by the Stanford Linear Accelerator Centre (SLAC). Indices of global and regional phase congestion for the monitoring sites are calculated from standardised regressions from the Earth’s rotation. There is also a J-curve limit to transferring information by the Internet and this is expressed as an inequality underpinned by the speed of light with examples from US and European traffic. The research returns to an often little known theme of Isaac Newton’s: mixing physics with geography. In our case, the equations define trajectories of information packets travelling near the speed of light, navigating within networks and between longitudes, relative to the Earth’s rotation.     

Estimating parameters in stochastic systems: A variational Bayesian approach

This work is concerned with approximate inference in dynamical systems, from a variational Bayesian perspective. When modelling real world dynamical systems, stochastic differential equations appear as a natural choice, mainly because of their ability to model the noise of the system by adding a variation of some stochastic process to the deterministic dynamics. Hence, inference in such processes has drawn much attention. Here a new extended framework is derived that is based on a local polynomial approximation of a recently proposed variational Bayesian algorithm. The paper begins by showing that the new extension of this variational algorithm can be used for state estimation (smoothing) and converges to the original algorithm. However, the main focus is on estimating the (hyper-) parameters of these systems (i.e. drift parameters and diffusion coefficients). The new approach is validated on a range of different systems which vary in dimensionality and non-linearity. These are the Ornstein-Uhlenbeck process, the exact likelihood of which can be computed analytically, the univariate and highly non-linear, stochastic double well and the multivariate chaotic stochastic Lorenz ’63 (3D model). As a special case the algorithm is also applied to the 40 dimensional stochastic Lorenz ’96 system. In our investigation we compare this new approach with a variety of other well known methods, such as the hybrid Monte Carlo, dual unscented Kalman filter, full weak-constraint 4D-Var algorithm and analyse empirically their asymptotic behaviour as a function of observation density or length of time window increases. In particular we show that we are able to estimate parameters in both the drift (deterministic) and the diffusion (stochastic) part of the model evolution equations using our new methods.     

Low-pressure RF discharge in the free-flight regime

The self-consistent equations system for low-pressure RF discharge in the free-flight regime is formulated. The expressions for the electron energy diffusion coefficient due to electron-neutral collisions and to the electron collisions with the plasma-space charge moving boundary (stochastic heating) are derived. If the electron-neutral elastic collisions frequency exceeds the inelastic one, the conventional two-term approximation for the electron distribution function (EDF) can be generalized, and the space-time-averaged electron kinetic equation can be reduced to the one-dimensional energy diffusion one. The fast electrons attached to the electrode surface can also be accounted for in this equation. It is shown that in the cases of (a) spatially uniform ion profile, (b) for frequencies that are small compared with the electron bounce frequency, and (c) for frequencies exceeding the electron plasma one in the sheath, the stochastic heating vanishes     

A path integral method for coarse-graining noise in stochastic differential equations with multiple time scales

We present a new path integral method to analyze stochastically perturbed ordinary differential equations with multiple time scales. The objective of this method is to derive from the original system a new stochastic differential equation describing the system’s evolution on slow time scales. For this purpose, we start from the corresponding path integral representation of the stochastic system and apply a multi-scale expansion to the associated path integral kernel of the corresponding Lagrangian. As a concrete example, we apply this expansion to a system that arises in the study of random dispersion fluctuations in dispersion-managed fiber-optic communications. Moreover, we show that, for this particular example, the new path integration method yields the same result at leading order as an asymptotic expansion of the associated Fokker-Planck equation.     

Multisoliton solutions and energy sharing collisions in coupled nonlinear Schrödinger equations with focusing,defocusing and mixed type nonlinearities

Bright and bright-dark type multisoliton solutions of the integrable N-coupled nonlinear Schrödinger (CNLS) equations with focusing, defocusing and mixed type nonlinearities are obtained by using Hirota’s bilinearization method. Particularly, for the bright soliton case, we present the Gram type determinant form of the n-soliton solution (n:arbitrary) for both focusing and mixed type nonlinearities and explicitly prove that the determinant form indeed satisfies the corresponding bilinear equations. Based on this, we also write down the multisoliton form for the mixed (bright-dark) type solitons. For the focusing and mixed type nonlinearities with vanishing boundary conditions the pure bright solitons exhibit different kinds of nontrivial shape changing/energy sharing collisions characterized by intensity redistribution, amplitude dependent phase-shift and change in relative separation distances. Due to nonvanishing boundary conditions the mixed N-CNLS system can admit coupled bright-dark solitons. Here we show that the bright solitons exhibit nontrivial energy sharing collision only if they are spread up in two or more components, while the dark solitons appearing in the remaining components undergo mere standard elastic collisions. Energy sharing collisions lead to exciting applications such as collision based optical computing and soliton amplification. Finally, we briefly discuss the energy sharing collision properties of the solitons of the (2+1) dimensional long wave-short wave resonance interaction (LSRI) system.     

Generalizing Planck’s distribution by using the Carati–Galgani model of molecular collisions

Classical systems of coupled harmonic oscillators are studied using the Carati–Galgani model. We investigate the consequences for Einstein’s conjecture by considering that the exchange of energy in molecular collisions follows the Lévy type statistics. We develop a generalization of Planck’s distribution admitting that there are analogous relations in the equilibrium quantum statistical mechanics of the relations found using the nonequilibrium classical statistical mechanics approach. The generalization of Planck’s law based on the nonextensive statistical mechanics formalism is compatible with our analysis.     

Modelling of linear wave propagation in spatial fluid filled pipe systems consisting of elastic curved and straight elements

This paper is concerned with the theoretical analysis of time harmonic dynamics of compound elastic pipes with and without internal fluid loading. Compound pipes are assembled as a sequence of segments, each of which has a constant curvature. As a prerequisite for the wave propagation analysis, dispersion equations are solved, Green’s matrices are formulated and Somigliana’s identities are derived for an isolated curved segment. The governing equations of wave motion of a compound pipe are obtained as an ensemble of the boundary integral equations for individual segments and the continuity conditions at their interfaces. The proposed methodology is validated in several benchmark problems and then applied for analysis of the periodicity effects. The results obtained for piping systems with a variable number of identical curved segments are put into the context of the classical Floquet theory. Brief parametric studies suggest that the curved inserts can be employed as a tool for the passive control of wave propagation in fluid-filled pipes, and their stop band characteristics may be tailored to reach desirable attenuation levels in prescribed frequency ranges.     

Effects of the spectral inhomogeneity on the optical properties in the four-wave mixing signal

Effects of the inhomogeneous broadening on the optical responses of a two-level system in the four-wave mixing signal (FWM) were studied in this work. In this model, we introduced in the optical Bloch equations a stochastic Bohr frequency, generated as a consequence of the solute-solvent molecular collisions. These expressions are averaged in the statistical ensemble over all realizations of the random variables. Using the convolution theorem, it is possible to incorporate the saturation effects of the electromagnetic field treated at all orders of perturbation. For the averages, we have considered two distributions of linewidth. From the analytical expressions, for both linear and nonlinear responses, numerical calculations were carried out to obtain surfaces as a function of the magnitude of the pump field and the detuning factor. Saturation effect is shown in the linear and the coupled nonlinear responses. In this work, we have already considered the inhomogeneous spectral line and studied the cases in which all the fields propagate through the medium.     

Understanding the accuracy of Nanbu’s numerical Coulomb collision operator

We investigate the accuracy of and assumptions underlying the numerical binary Monte Carlo collision operator due to Nanbu [K. Nanbu, Phys. Rev. E 55 (1997) 4642]. The numerical experiments that resulted in the parameterization of the collision kernel used in Nanbu’s operator are argued to be an approximate realization of the Coulomb–Lorentz pitch-angle scattering process, for which an analytical solution for the collision kernel is available. It is demonstrated empirically that Nanbu’s collision operator quite accurately recovers the effects of Coulomb–Lorentz pitch-angle collisions, or processes that approximate these (such interspecies Coulomb collisions with very small mass ratio) even for very large values of the collisional time step. An investigation of the analytical solution shows that Nanbu’s parameterized kernel is highly accurate for small values of the normalized collision time step, but loses some of its accuracy for larger values of the time step. Careful numerical and analytical investigations are presented, which show that the time dependence of the relaxation of a temperature anisotropy by Coulomb–Lorentz collisions has a richer structure than previously thought, and is not accurately represented by an exponential decay with a single decay rate. Finally, a practical collision algorithm is proposed that for small-mass-ratio interspecies Coulomb collisions improves on the accuracy of Nanbu’s algorithm.     

Statistical Dynamics of Solitons in Optical Fibers

The soliton perturbation theory is used to study and analyze the stochastic perturbation of optical solitons in addition to deterministic perturbations of optical solitons that are governed by the nonlinear Schro¨dinger's equation. The Langevin equations are derived and analyzed. The deterministic perturbations that are considered here are of both Hamiltonian as well as of non-Hamiltonian type.     

Suppression and enhancement of impurity scattering in a Bose-Einstein condensate

Impurity atoms propagating at variable velocities through a trapped Bose-Einstein condensate were produced using a stimulated Raman transition. The redistribution of momentum by collisions between the impurity atoms and the stationary condensate was observed in a time-of-flight analysis. The collisional cross section was dramatically reduced when the impurity velocity was reduced below the condensate speed of sound, in agreement with the Landau criterion for superfluidity. For large numbers of impurity atoms, we observed an enhancement of atomic collisions due to bosonic stimulation. This enhancement is analogous to optical super-radiance.     

Sequential reconstruction of driving-forces from nonlinear nonstationary dynamics

This paper describes a functional analysis-based method for the estimation of driving-forces from nonlinear dynamic systems. The driving-forces account for the perturbation inputs induced by the external environment or the secular variations in the internal variables of the system. The proposed algorithm is applicable to the problems for which there is too little or no prior knowledge to build a rigorous mathematical model of the unknown dynamics. We derive the estimator conditioned on the differentiability of the unknown system’s mapping, and smoothness of the driving-force. The proposed algorithm is an adaptive sequential realization of the blind prediction error method, where the basic idea is to predict the observables, and retrieve the driving-force from the prediction error. Our realization of this idea is embodied by predicting the observables one-step into the future using a bank of echo state networks (ESN) in an online fashion, and then extracting the raw estimates from the prediction error and smoothing these estimates in two adaptive filtering stages. The adaptive nature of the algorithm enables to retrieve both slowly and rapidly varying driving-forces accurately, which are illustrated by simulations. Logistic and Moran-Ricker maps are studied in controlled experiments, exemplifying chaotic state and stochastic measurement models. The algorithm is also applied to the estimation of a driving-force from another nonlinear dynamic system that is stochastic in both state and measurement equations. The results are judged by the posterior Cramer-Rao lower bounds. The method is finally put into test on a real-world application; extracting sun’s magnetic flux from the sunspot time series.     

Beat Heating of Collisional and Magnetised Plasma

The heating of a plasma by stimulating plasma electrons as well as plasma ions with two anti-parallel electromagnetic waves under the influence of a uniform static magnetic field is studied using Maxwell's equations and equations of motion. A formula for the power absorption per unit volume of the plasma is derived and effects of collisions and magnitude and orientation of the magnetic field on the beat heating are examined numerically. It is observed that the average power absorption in the absence of ion-neutral collisions in the plasma barely exceeds unity in the units of pure Langmuir mode excitation where as in the presence of ion-neutral collisions the power absorption immediately shoots up to a very high value.