Green function for the diffusion limit of one-dimensional continuous time random walks

It is shown how the fractional probability density diffusion equation for the diffusion limit of one-dimensional continuous time random walks may be derived from a generalized Markovian Chapman-Kolmogorov equation. The non-Markovian behaviour is incorporated into the Markovian Chapman-Kolmogorov equation by postulating a Levy like distribution of waiting times as a kernel. The Chapman-Kolmogorov equation so generalised then takes on the form of a convolution integral. The dependence on the initial conditions typical of a non-Markovian process is treated by adding a time dependent term involving the survival probability to the convolution integral. In the diffusion limit these two assumptions about the past history of the process are sufficient to reproduce anomalous diffusion and relaxation behaviour of the Cole-Cole type. The Green function in the diffusion limit is calculated using the fact that the characteristic function is the Mittag-Leffler function. Fourier inversion of the characteristic function yields the Green function in terms of a Wright function. The moments of the distribution function are evaluated from the Mittag-Leffler function using the properties of characteristic functions and a relation between the powers of the second moment and higher order even moments is derived. (C) 2004 Elsevier B.V. All rights reserved.

Escape times for rigid Brownian rotators in a bistable potential from the time evolution of the Green function and the characteristic time of the probability evolution

The greatest relaxation time for an assembly of three- dimensional rigid rotators in an axially symmetric bistable potential is obtained exactly in terms of continued fractions as a sum of the zero frequency decay functions (averages of the Legendre polynomials) of the system. This is accomplished by studying the entire time evolution of the Green function (transition probability) by expanding the time dependent distribution as a Fourier series and proceeding to the zero frequency limit of the Laplace transform of that distribution. The procedure is entirely analogous to the calculation of the characteristic time of the probability evolution (the integral of the configuration space probability density function with respect to the position co-ordinate) for a particle undergoing translational diffusion in a potential; a concept originally used by Malakhov and Pankratov (Physica A 229 (1996) 109). This procedure allowed them to obtain exact solutions of the Kramers one-dimensional translational escape rate problem for piecewise parabolic potentials. The solution was accomplished by posing the problem in terms of the appropriate Sturm-Liouville equation which could be solved in terms of the parabolic cylinder functions. The method (as applied to rotational problems and posed in terms of recurrence relations for the decay functions, i.e., the Brinkman approach c.f. Blomberg, Physica A 86 (1977) 49, as opposed to the Sturm-Liouville one) demonstrates clearly that the greatest relaxation time unlike the integral relaxation time which is governed by a single decay function (albeit coupled to all the others in non-linear fashion via the underlying recurrence relation) is governed by a sum of decay functions. The method is easily generalized to multidimensional state spaces by matrix continued fraction methods allowing one to treat non-axially symmetric potentials, where the distribution function is governed by two state variables. (C) 2001 Elsevier Science B.V. All rights reserved.

Spectroscopic measurements of phase-resolved electron energy distribution functions in RF-excited discharges

The reliable measurement of the electron energy distribution function (EEDF) of plasmas is one of the most important subjects of plasma diagnostics, because this piece of information is the key to understand basic discharge mechanisms. Specific problems arise in the case of RF-excited plasmas, since the properties of electrons are subject to changes on a nanosecond time scale and show pronounced spatial anisotropy. We report on a novel spectroscopic method for phase- and space-resolved measurements of the electron energy distribution function of energetic (> 12 eV) electrons in RF discharges. These electrons dominate excitation and ionization processes and are therefore of particular interest. The technique is based on time-dependent measurements during the RF cycle of excited-state populations of rare gases admixed in small fractions. These measurements yield � in combination with an analytical model � detailed information on the excitation processes. Phase-resolved optical emission spectroscopy allows us to overcome the difficulties connected with the very low densities (107�109 cm�3) and the transient character of the electrons in the sheath region. The EEDF of electrons accelerated in the sheath region can be described by a shifted Maxwellian with a drift velocity component in direction of the electric field. The method yields the high-energy tail of the EEDF on an absolute scale. The applicability of the method is demonstrated at a capacitively coupled RF discharge in hydrogen.

Analysis procedure for calculation of electron energy distribution functions from incoherent Thomson scattering spectra

Incoherent Thomson scattering (ITS) provides a nonintrusive diagnostic for the determination of one-dimensional (1D) electron velocity distribution in plasmas. When the ITS spectrum is Gaussian its interpretation as a three-dimensional (3D) Maxwellian velocity distribution is straightforward. For more complex ITS line shapes derivation of the corresponding 3D velocity distribution and electron energy probability distribution function is more difficult. This article reviews current techniques and proposes an approach to making the transformation between a 1D velocity distribution and the corresponding 3D energy distribution. Previous approaches have either transformed the ITS spectra directly from a 1D distribution to a 3D or fitted two Gaussians assuming a Maxwellian or bi-Maxwellian distribution. Here, the measured ITS spectrum transformed into a 1D velocity distribution and the probability of finding a particle with speed within 0 and given value v is calculated. The differentiation of this probability function is shown to be the normalized electron velocity distribution function. (C) 2003 American Institute of Physics.

Comment on “Mathematical and physical aspects of Kappa velocity distribution” [ Phys. Plasmas 14, 110702 (2007) ]

A recent paper [L.-N. Hau and W.-Z. Fu, Phys. Plasmas 14, 110702 (2007)] deals with certain mathematical and physical properties of the kappa distribution. We comment on the authors' use of a form of distribution function that is different from the "standard" form of the kappa distribution, and hence their results, inter alia for an expansion of the distribution function and for the associated number density in an electrostatic potential, do not fully reflect the dependence on kappa that would be associated with the conventional kappa distribution. We note that their definition of the kappa distribution function is also different from a modified distribution based on the notion of nonextensive entropy.

Improved Nonlinear PCA for Process Monitoring Using Support Vector Data Description

Nonlinear principal component analysis (PCA) based on neural networks has drawn significant attention as a monitoring tool for complex nonlinear processes, but there remains a difficulty with determining the optimal network topology. This paper exploits the advantages of the Fast Recursive Algorithm, where the number of nodes, the location of centres, and the weights between the hidden layer and the output layer can be identified simultaneously for the radial basis function (RBF) networks. The topology problem for the nonlinear PCA based on neural networks can thus be solved. Another problem with nonlinear PCA is that the derived nonlinear scores may not be statistically independent or follow a simple parametric distribution. This hinders its applications in process monitoring since the simplicity of applying predetermined probability distribution functions is lost. This paper proposes the use of a support vector data description and shows that transforming the nonlinear principal components into a feature space allows a simple statistical inference. Results from both simulated and industrial data confirm the efficacy of the proposed method for solving nonlinear principal component problems, compared with linear PCA and kernel PCA.

Computational Verification of Two Universal Relations for Simple Ionic Liquids. Kinetic Properties of a Model 2:1 Molten Salt

Two semianalytical relations [Nature, 1996, 381, 137 and Phys. Rev. Lett. 2001, 87, 245901] predicting dynamical coefficients of simple liquids on the basis of structural properties have been tested by extensive molecular dynamics simulations for an idealized 2:1 model molten salt. In agreement with previous simulation studies, our results support the validity of the relation expressing the self-diffusion coefficient as a Function of the radial distribution functions for all thermodynamic conditions such that the system is in the ionic (ie., fully dissociated) liquid state. Deviations are apparent for high-density samples in the amorphous state and in the low-density, low-temperature range, when ions condense into AB(2) molecules. A similar relation predicting the ionic conductivity is only partially validated by our data. The simulation results, covering 210 distinct thermodynamic states, represent an extended database to tune and validate semianalytical theories of dynamical properties and provide a baseline for the interpretation of properties of more complex systems such as the room-temperature ionic liquids.

ELECTRON-ENERGY DISTRIBUTION-FUNCTIONS AND NEGATIVE-ION CONCENTRATIONS IN TANDEM AND HYBRID MULTICUSP NEGATIVE HYDROGEN-ION SOURCES

The second derivative of a Langmuir probe characteristic is used to establish the electron energy distribution function (EEDF) in both a tandem and hybrid multicusp H- ion source. Moveable probes are used to establish the spatial variation of the EEDF. The negative ion density is measured by laser induced photo-detachment. In the case of the hybrid source the EEDF consists of a cold Maxwellian in the central region of the source; the electron temperature increases with increasing discharge current (rising from 0.3 eV at 1 A to 1.2 eV at 50 A when the pressure is 0.4 Pa). A hot-electron tail exists in the EEDF of the driver region adjacent to each filament which is shown to consist of a distinct group of primary electrons at low pressure (0.08 Pa) but becomes degraded mainly through inelastic collisions at higher pressures (0.27 Pa). The tandem source, on the other hand, has a single driver region which extends throughout the central region. The primary electron confinement times are much longer so that even at the lowest pressure considered (0.07 Pa) the primaries are degraded. In both cases the measured EEDF at specific locations and values of discharge operating parameters are used to establish the rate coefficients for the processes of importance in H- production and destruction.

Distribution Network Voltage Support Using Sensitivity-based Dispatch of Distributed Generation

This paper is concerned with the voltage and reactive power issues surrounding the connection of Distributed Generation (DG) on the low-voltage (LV) distribution network. The presented system-wide voltage control algorithm consists of three stages. Firstly available reactive power reserves are utilized. Then, if required, DG active power output is curtailed. Finally, curtailment of non-critical site demand is considered. The control methodology is tested on a variant of the 13-bus IEEE Node Radial Distribution Test Feeder. The presented control algorithm demonstrated that the distribution system operator (DSO) can maintain voltage levels within a desired statutory range by dispatching reactive power from DG or network devices. The practical application of the control strategy is discussed.

On the Multivariate Gamma-Gamma (ΓΓ) Distribution with Arbitrary Correlation and Applications in Wireless Communications

The statistical properties of the multivariate GammaGamma (ΓΓ) distribution with arbitrary correlation have remained unknown. In this paper, we provide analytical expressions for the joint probability density function (PDF), cumulative distribution function (CDF) and moment generation function of the multivariate ΓΓ distribution with arbitrary correlation. Furthermore, we present novel approximating expressions for the PDF and CDF of the su m of ΓΓ random variables with arbitrary correlation. Based on this statistical analysis, we investigate the performance of radio frequency and optical wireless communication systems. It is noteworthy that the presented expressions include several previous results in the literature as special cases.

A Hybrid Forward Algorithm for RBF Neural Network Construction

This paper proposes a novel hybrid forward algorithm (HFA) for the construction of radial basis function (RBF) neural networks with tunable nodes. The main objective is to efficiently and effectively produce a parsimonious RBF neural network that generalizes well. In this study, it is achieved through simultaneous network structure determination and parameter optimization on the continuous parameter space. This is a mixed integer hard problem and the proposed HFA tackles this problem using an integrated analytic framework, leading to significantly improved network performance and reduced memory usage for the network construction. The computational complexity analysis confirms the efficiency of the proposed algorithm, and the simulation results demonstrate its effectiveness

Solvation effects on equilibria: Triazoles and N-methyl piperidinol

We have performed calculations of the solvation effects on a number of equilibrium constants in water using a recently proposed hybrid quantum classical scheme in which the liquid environment is modelled using classical solvent molecules and the solute electronic structure is computed using modern quantum chemical methods. The liquid phase space is sampled from a fully classical simulation. We find that solvation effects on both triazole tautomeric equilibrium constants and piperidinol conformational equilibrium constants can be interpreted in terms of subtle differences in the local environment which can be seen in probability densities and radial distribution functions. Lower level calculations were performed for comparison and we conclude that the solvation thermodynamics can be predicted from a good classical model of solvent and solute molecules, but the implicit models that we tried are less successful.

Development of complex classical force fields through force matching to ab initio data: Application to a room-temperature ionic liquid

Recent experimental neutron diffraction data and ab initio molecular dynamics simulation of the ionic liquid dimethylimidazolium chloride ([dmim]Cl) have provided a structural description of the system at the molecular level. However, partial radial distribution functions calculated from the latter, when compared to previous classical simulation results, highlight some limitations in the structural description offered by force fieldbased simulations. With the availability of ab initio data it is possible to improve the classical description of [dmim]Cl by using the force matching approach, and the strategy for fitting complex force fields in their original functional form is discussed. A self-consistent optimization method for the generation of classical potentials of general functional form is presented and applied, and a force field that better reproduces the observed first principles forces is obtained. When used in simulation, it predicts structural data which reproduces more faithfully that observed in the ab initio studies. Some possible refinements to the technique, its application, and the general suitability of common potential energy functions used within many ionic liquid force fields are discussed.

Energy analysis of hyperthermal hydrogen atoms generated through surface neutralisation of ions

Hydrogen ions (H+, H-2(+) and H-3(+)) are produced in a magnetically confined inductively coupled radio frequency plasma. Ions are accelerated in the plasma boundary sheath potential, of several hundred volts, in front of a biased metal electrode immersed in the plasma. Backscattered hyperthermal hydrogen atoms are investigated by optical emission spectroscopy and an energy-resolved mass spectrometer. Ionisation of fast neutrals through electron stripping of atoms in the plasma allows energy analysis of the resulting ions. Thereby, the energy distribution function of the hyperthermal atoms can be deduced. The energy spectra can be explained as a superposition of individual spectra of the various ion species. The measured spectra also shows contributions of negative ions created at the electrode surface. In addition to experimental measurements, simulations of the neutral flux of backscattered atoms are carried out.

Prospects of phase resolved optical emission spectroscopy as a powerful diagnostic tool for RF-discharges

Phase resolved optical emission spectroscopy (PROES) bears considerable potential for diagnostics of RF discharges that give detailed insight of spatial and temporal variations of excitation processes. Based on phase and space resolved measurements of the population dynamics of excited states several diagnostic techniques have been developed. Results for a hydrogen capacitively coupled RF (CCRF) discharge are discussed as an example. The gas temperature, the degree of dissociation and the temporally and spatially resolved electron energy distribution function (EEDF) of energetic electrons (>12eV) are measured. Furthermore, the pulsed electron impact excitation during the field reversal phase, typical for hydrogen CCRF discharges, is exploited for measurements of atomic and molecular data like lifetimes of excited states, coefficients for radiationless collisional de-excitation (quenching coefficients), and cascading processes from higher electronic states.