Restricted rotation dynamics and far-infrared spectra of liquid water governed by elastic interactions

A semi-phenomenological model describing wideband dielectric and far-infrared spectra of liquid water was proposed recently by the same authors [J. Mol. Struct. 606 (2002) 9], where a small dipole-moment component changing harmonically with time determines a weak absorption band (termed here the R-band) centred at the wavenumber v similar to 200 cm(-1). In the present work, a rough molecular theory of the R-band based on the concept of elastic interactions is given. Stretching and bending of hydrogen bonds cause restricted rotation (RR) of a polar water molecule in terms of a dimer comprising the H- bonded molecules. Analytical expression for the RR frequency nu(str) is derived as a function of the RR amplitude, geometrical parameters and force constants. The density g(nu(str)) of frequency distribution is shown to be centred in the R-band. The spectrum of the dipolar auto-correlation function calculated for this structural-dynamical model is found. A composite model comprising two intermolecular potentials is proposed, which yields for water a good description of the experimental wideband (from 0 to 1000 cm(- 1)) spectra of complex permittivity and of absorption coefficient. The presented interpretation of these spectra is based on a concept that water presents a two-component solution, with components differing by the types of molecular rotation. (C) 2003 Elsevier B.V. All rights reserved.

Specific and unspecific interactions in polar fluids in view of wideband dielectric far-infrared spectra

A simple molecular analytical theory of dielectric relaxation in strongly polar fluids is considered in terms of a semi- phenomenological approach. Theoretical spectra epsilon(v), a(v) of complex permittivity and absorption coefficient are fully determined by a form of intermolecular potential well, in which a dipole reorients. In a recent publication by VI. Gaiduk, O.F. Nielsen, and T.S. Perova [J. Molliq 95 (1002) 1-25] the wideband spectra of liquid H2O and D2O were described in terms of a composite model comprising the rectangular and the cosine squared potential wells. Much better results are achieved in this work, where the rectangular well is replaced by a well with a rounded bottom termed the hat-curved well. The spectrum of the auto-correlation function (ACF) is calculated for such a potential. The proposed theory of a composite model, comprising hat-curved and parabolic wells, is applied for liquid water. This model is capable for describing the Debye relaxation region, the second relaxation region in the submillimeter wavelength range, and the far infra-red (FIR) e(v), a(v) spectra, where an intense librational band and an additional weak band are placed, respectively, near 700 cm(-1) and 200 cm(-1). The latter band reflects the features of so-called specific (viz. directly related to H-bonds) interactions and the former band reflects the features of unspecific interactions. The physical mechanisms connected with these types of interactions are discussed in terms of two relevant types of water structure (types of molecular rotation). The proposed theory is also applied to a non-associated liquid in terms of one hat-curved potential well. (C) 2004 Elsevier B.V. All rights reserved.

Statistical Monitoring of Dynamic Multivariate Processes: Part I: Modeling Auto- and Cross-Correlation

The work in this paper is of particular significance since it considers the problem of modelling cross- and auto-correlation in statistical process monitoring. The presence of both types of correlation can lead to fault insensitivity or false alarms, although in published literature to date, only autocorrelation has been broadly considered. The proposed method, which uses a Kalman innovation model, effectively removes both correlations. The paper (and Part 2 [2]) has emerged from work supported by EPSRC grant GR/S84354/01 and is of direct relevance to problems in several application areas including chemical, electrical, and mechanical process monitoring.

On the structure and dynamics of ionic liquids

The structure and dynamics of the ionic liquid 1-ethyl-3-methylimidazolium nitrate is studied by molecular dynamics simulations. We find long-range spatial correlations between the ions and a three-dimensional local structure that reflects the asymmetry of the cations. The main contribution to the configurational energy comes from the electrostatic interactions which leads to charge-ordering effects. Radial screening and threedimensional distribution of charge are also analyzed. The motion of a single ion is studied via velocity and reorientational correlation functions. It is found that ions "rattle" in a long-lived cage, while the orientational structure relaxes on a time scale longer than 200 ps. As in a supercooled liquid, the mean square displacements reveal a subdiffusive dynamics. In addition, the presence of dynamic heterogeneities can be detected by analyzing the non-Gaussian behavior of the van Hove correlation function and the spatial arrangement of the most mobile ions. The short-time collective dynamics is also studied through the electric current time correlation function.

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.

Thermodynamic properties and atomistic structure of the dry amorphous silica surface from a reactive force field model

A force field model of the Keating type supplemented by rules to break, form, and interchange bonds is applied to investigate thermodynamic and structural properties of the amorphous SiO2 surface. A simulated quench from the liquid phase has been carried out for a silica sample made of 3888 silicon and 7776 oxygen atoms arranged on a slab similar to 40 angstrom thick, periodically repeated along two directions. The quench results into an amorphous sample, exposing two parallel square surfaces of similar to 42 nm(2) area each. Thermal averages computed during the quench allow us to determine the surface thermodynamic properties as a function of temperature. The surface tension turns out to be gamma=310 +/- 20 erg/cm(2) at room temperature and gamma=270 +/- 30 at T=2000 K, in fair agreement with available experimental estimates. The entropy contribution Ts-s to the surface tension is relatively low at all temperatures, representing at most similar to 20% of the surface energy. Almost without exceptions, Si atoms are fourfold coordinated and oxygen atoms are twofold coordinated. Twofold and threefold rings appear only at low concentration and are preferentially found in proximity of the surface. Above the glass temperature T-g=1660 +/- 50 K, the mobility of surface atoms is, as expected, slightly higher than that of bulk atoms. The computation of the height-height correlation function shows that the silica surface is rough in the equilibrium and undercooled liquid phase, becoming smooth below the glass temperature T-g.

Improved Principal Component Monitoring of Large-Scale Processes

This is the first paper that shows and theoretically analyses that the presence of auto-correlation can produce considerable alterations in the Type I and Type II errors in univariate and multivariate statistical control charts. To remove this undesired effect, linear inverse ARMA filter are employed and the application studies in this paper show that false alarms (increased Type I errors) and an insensitive monitoring statistics (increased Type II errors) were eliminated.

The role of Mittag-Leffler functions in anomalous relaxation

The inertia-corrected Debye model of rotational Brownian motion of polar molecules was generalized by Coffey et al. [Phys. Rev. E, 65, 32 102 (2002)] to describe fractional dynamics and anomalous rotational diffusion. The linear- response theory of the normalized complex susceptibility was given in terms of a Laplace transform and as a function of frequency. The angular-velocity correlation function was parametrized via fractal Mittag-Leffler functions. Here we apply the latter method and complex-contour integral- representation methods to determine the original time-dependent amplitude as an inverse Laplace transform using both analytical and numerical approaches, as appropriate. (C) 2004 Elsevier B.V. All rights reserved.

Thermal and quantum fluctuations in chains of ultracold polar molecules

Ultracold polar molecules, in highly anisotropic traps and interacting via a repulsive dipolar potential, may form one-dimensional chains at high densities. According to classical theory, at low temperatures there exists a critical value of the density at which a second-order phase transition from a linear to a zigzag chain occurs. We study the effect of thermal and quantum fluctuations on these self-organized structures using classical and quantum Monte Carlo methods, by means of which we evaluate the pair correlation function and the static structure factor. Depending on the parameters, these functions exhibit properties typical of a crystalline or of a liquid system. We compare the thermal and the quantum results, identifying analogies and differences. Finally, we discuss experimental parameter regimes where the effects of quantum fluctuations on the linear-zigzag transition can be observed.

Ramsey interferometry with a spin embedded in a Coulomb chain

We show that the statistical properties of a Coulomb crystal can be measured by means of a standard interferometric procedure performed on the spin of one ion in the chain. The ion spin, constituted by two internal levels of the ion, couples to the crystal modes via spatial displacement induced by photon absorption. The loss of contrast in the interferometric signal allows one to measure the autocorrelation function of the crystal observables. Close to the critical point, where the chain undergoes a second-order phase transition to a zigzag structure, the signal gives the behavior of the correlation function at the critical point.

Nonequilibrium critical scaling from quantum thermodynamics

The emerging field of quantum thermodynamics is contributing important results and insights into archetypal many-body problems, including quantum phase transitions. Still, the question whether out-of-equilibrium quantities, such as fluctuations of work, exhibit critical scaling after a sudden quench in a closed system has remained elusive. Here, we take a novel approach to the problem by studying a quench across an impurity quantum critical point. By performing density matrix renormalization group computations on the two-impurity Kondo model, we are able to establish that the irreversible work produced in a quench exhibits finite-size scaling at quantum criticality. This scaling faithfully predicts the equilibrium critical exponents for the crossover length and the order parameter of the model, and, moreover, implies a new exponent for the rescaled irreversible work. By connecting the irreversible work to the two-impurity spin correlation function, our findings can be tested experimentally.

Exchange-correlation energy and potential as approximate functionals of occupied and virtual Kohn-Sham orbitals: Application to dissociating H-2

The standard local density approximation and generalized gradient approximations fail to properly describe the dissociation of an electron pair bond, yielding large errors (on the order of 50 kcal/mol) at long bond distances. To remedy this failure, a self-consistent Kohn-Sham (KS) method is proposed with the exchange-correlation (xc) energy and potential depending on both occupied and virtual KS orbitals. The xc energy functional of Buijse and Baerends [Mol. Phys. 100, 401 (2002); Phys. Rev. Lett. 87, 133004 (2001)] is employed, which, based on an ansatz for the xc-hole amplitude, is able to reproduce the important dynamical and nondynamical effects of Coulomb correlation through the efficient use of virtual orbitals. Self-consistent calculations require the corresponding xc potential to be obtained, to which end the optimized effective potential (OEP) method is used within the common energy denominator approximation for the static orbital Green's function. The problem of the asymptotic divergence of the xc potential of the OEP when a finite number of virtual orbitals is used is addressed. The self-consistent calculations reproduce very well the entire H-2 potential curve, describing correctly the gradual buildup of strong left-right correlation in stretched H-2. (C) 2003 American Institute of Physics.

TIME-RESOLVED OPTICAL-EMISSION AND ELECTRON-ENERGY DISTRIBUTION FUNCTION MEASUREMENTS IN RF PLASMAS

Comparisons between experimentally measured time-dependent electron energy distribution functions and optical emission intensities are reported for low-frequency (100 and 400 kHz) radio-frequency driven discharges in argon. The electron energy distribution functions were measured with a time-resolved Langmuir probe system. Time-resolved optical emissions of argon resonance lines at 687.1 and 750.4 nm were determined by photon-counting methods. Known ground-state and metastable-state excitation cross sections were used along with the measured electron energy distribution functions to calculate the time dependence of the optical emission intensity. It was found that a calculation using only the ground-state cross sections gave the best agreement with the time dependence of the measured optical emission. Time-dependent electron density, electron temperature, and plasma potential measurements are also reported.

TIME-RESOLVED ELECTRON-ENERGY DISTRIBUTION FUNCTION MEASUREMENTS IN A PULSED MAGNETIC MULTIPOLE HYDROGEN DISCHARGE

The time evolution of measured plasma parameters, including the electron energy distribution function (EEDF), in the discharge and post-discharge regime of a pulsed hydrogen magnetic multipole plasma is presented. The time necessary for the plasma to reach equilibrium has been established as 160-mu-s. The present results clarify the mechanisms which initiate the discharge. The decay rates of the charged-particle density and energy in the post-discharge have been measured. These measurements indicate that particle transport to the wall is the dominant loss mechanism for both charged-particle density and energy. The time-resolved EEDF is found to be non-Maxwellian in the discharge and Maxwellian in the late post-discharge.