Single ionization of helium by antiproton impact

We study the ionization of helium (fie) in collision with antiprotons (p) in the energy range from 10 keV to 1000 keV. We adopt a semiclassical single center close coupling approach in which the wave function for the electron is expanded in a B-spline basis centered on the nucleus of the He atom, The calculations are performed using two different models: the independent particle (IP) model and the one-electron (OE) approximation. The interaction between the active electron and the rest of the atom, i.e. passive electron and nucleus, is represented by a model potential. The results obtained are compared with experimental data as well as with existing theoretical calculations. (c) 2005 Published by Elsevier B.V.

A variational approach to the relaxation of single domain magnetic particles based on Brown's model

Brown's model for the relaxation of the magnetization of a single domain ferromagnetic particle is considered. This model results in the Fokker-Planck equation of the process. The solution of this equation in the cases of most interest is non- trivial. The probability density of orientations of the magnetization in the Fokker-Planck equation can be expanded in terms of an infinite set of eigenfunctions and their corresponding eigenvalues where these obey a Sturm-Liouville type equation. A variational principle is applied to the solution of this equation in the case of an axially symmetric potential. The first (non-zero) eigenvalue, corresponding to the largest time constant, is considered. From this we obtain two new results. Firstly, an approximate minimising trial function is obtained which allows calculation of a rigorous upper bound. Secondly, a new upper bound formula is derived based on the Euler-Lagrange condition. This leads to very accurate calculation of the eigenvalue but also, interestingly, from this, use of the simplest trial function yields an equivalent result to the correlation time of Coffey et at. and the integral relaxation time of Garanin. (C) 2004 Elsevier B.V. All rights reserved.

Interpolation formula between very low and intermediate-to-high damping Kramers escape rates for single-domain ferromagnetic particles

It is shown that the Mel'nikov-Meshkov formalism for bridging the very low damping (VLD) and intermediate-to-high damping (IHD) Kramers escape rates as a function of the dissipation parameter for mechanical particles may be extended to the rotational Brownian motion of magnetic dipole moments of single-domain ferromagnetic particles in nonaxially symmetric potentials of the magnetocrystalline anisotropy so that both regimes of damping, occur. The procedure is illustrated by considering the particular nonaxially symmetric problem of superparamagnetic particles possessing uniaxial anisotropy subject to an external uniform field applied at an angle to the easy axis of magnetization. Here the Mel'nikov-Meshkov treatment is found to be in good agreement with an exact calculation of the smallest eigenvalue of Brown's Fokker-Planck equation, provided the external field is large enough to ensure significant departure from axial symmetry, so that the VLD and IHD formulas for escape rates of magnetic dipoles for nonaxially symmetric potentials are valid.

Microscopic models for dielectric relaxation in disordered systems

It is shown how the Debye rotational diffusion model of dielectric relaxation of polar molecules (which may be described in microscopic fashion as the diffusion limit of a discrete time random walk on the surface of the unit sphere) may be extended to yield the empirical Havriliak-Negami (HN) equation of anomalous dielectric relaxation from a microscopic model based on a kinetic equation just as in the Debye model. This kinetic equation is obtained by means of a generalization of the noninertial Fokker-Planck equation of conventional Brownian motion (generally known as the Smoluchowski equation) to fractional kinetics governed by the HN relaxation mechanism. For the simple case of noninteracting dipoles it may be solved by Fourier transform techniques to yield the Green function and the complex dielectric susceptibility corresponding to the HN anomalous relaxation mechanism.

Competitive adsorption of reactive dyes from solution: equilibrium isotherm studies in single and multisolute systems

Experimental data of the adsorption of reactive dyestuffs onto Filtrasorb 400 activated carbon (FS400) were determined in an equilibrium isotherm study. As most industrial wastewater contains more than one pollutant, an investigation into the effect of multisolute systems (using the unhydrolysed form of the reactive dyes) on the adsorption capacity was undertaken. Equilibrium isotherm models were employed to describe the adsorption capacities of single, binary and ternary dye solutions. The results of these analyses showed that adsorption of reactive dyes from single and multisolute systems can be successfully described by Langmuir, and Redlich–Peterson equilibrium isotherm models. Experimental data indicated that competitive adsorption for active sites on the carbon surface results in a reduction in the overall uptake capacity of the reactive dyes investigated.

Solving the advection-dispersion equation using the discrete puff particle method

A flexible, mass-conservative numerical technique for solving the advection-dispersion equation for miscible contaminant transport is presented. The method combines features of puff transport models from air pollution studies with features from the random walk particle method used in water resources studies, providing a deterministic time-marching algorithm which is independent of the grid Peclet number and scales from one to higher dimensions simply. The concentration field is discretised into a number of particles, each of which is treated as a point release which advects and disperses over the time interval. The dispersed puff is itself discretised into a spatial distribution of particles whose masses can be pre-calculated. Concentration within the simulation domain is then calculated from the mass distribution as an average over some small volume. Comparison with analytical solutions for a one-dimensional fixed-duration concentration pulse and for two-dimensional transport in an axisymmetric flow field indicate that the algorithm performs well. For a given level of accuracy the new method has lower computation times than the random walk particle method.

Comparison of optimised isotherm models for basic dye adsorption by kudzu

This study assesses the use of dried (5% w/w moisture) kudzu (Peuraria lobata ohwi) as an adsorbent medium for the removal of two basic dyes, Basic Yellow 21 and Basic Red 22, from aqueous solutions. The extent of adsorption was measured through equilibrium sorption isotherms for the single component systems. Equilibrium was achieved after 21 days. The experimental isotherm data were analysed using Langmuir, Freundlich, Redlich-Peterson, Temkin and Toth isotherm equations. A detailed error analysis was undertaken to investigate the effect of using different error criteria for the determination of the single component isotherm parameters. The performance of the kudzu was compared with an activated carbon (Chemviron F-400). Kudzu was found to be an effective adsorbent for basic dye colour removal, though its capacity for colour removal was not as high as an activated carbon, the potential appeared to exist to use it as an alternative to activated carbon where carbon cost was prohibitive. (C) 2002 Elsevier Science Ltd. All rights reserved.

Mixed-dimensional coupling in finite element models

In many finite element analysis models it would be desirable to combine reduced- or lower-dimensional element types with higher-dimensional element types in a single model. In order to achieve compatibility of displacements and stress equilibrium at the junction or interface between the differing element types, it is important in such cases to integrate into the analysis some scheme for coupling the element types. A novel and effective scheme for establishing compatibility and equilibrium at the dimensional interface is described and its merits and capabilities are demonstrated. Copyright (C) 2000 John Wiley & Sons, Ltd.

Low SNR Capacity for MIMO Rician and Rayleigh-Product Fading Channels with Single Co-channel Interferer and Noise

This paper studies the ergodic capacity of multiple-input multiple-output (MIMO) systems with a single co-channel interferer in the low signal-to-noise-ratio (SNR) regime. Two MIMO models namely Rician and Rayleigh-product channels are investigated. Exact analytical expressions for the minimum energy per information bit, Eb/N0min, and wideband slope, S0, are derived for both channels. Our results show that the minimum energy per information bit is the same for both channels while their wideband slopes differ significantly. Further, the impact of the numbers of transmit and receive antennas, the Rician K factor, the channel mean matrix and the interference-to-noise-ratio (INR) on the capacity, is addressed. Results indicate that interference degrades the capacity by increasing the required minimum energy per information bit and reducing the wideband slope. Simulation results validate our analytical results.

Numerical simulations of fluid-structure interactions in single-reed mouthpieces

Most single-reed woodwind instrument models rely on a quasistationary approximation to describe the relationship between the volume flow and. the pressure difference across the reed channel. Semiempirical models based on the quasistationary approximation are very useful in explaining the fundamental characteristics of this family of instruments such as self-sustained oscillations and threshold of blowing pressure. However, they fail at explaining more complex phenomena associated with the fluid-structure interaction during dynamic flow regimes, such as the transient and steady-state behavior of the system as a function. of the mouthpiece geometry. Previous studies have discussed the accuracy of the quasistationary approximation but the amount of literature on the subject is sparse, mainly due to the difficulties involved in the measurement of dynamic flows in channels with an oscillating reed. In this paper, a numerical technique based on the lattice Boltzmann method and a finite difference scheme is proposed in order to investigate the characteristics of fully coupled fluid-structure interaction in single-reed mouthpieces with different channel configurations. Results obtained for a stationary simulation with a static reed agree very well with those predicted by the literature based on the quasistationary approximation. However, simulations carried out for a dynamic regime with dn oscillating reed show that the phenomenon associated with flow detachment and reattachment diverges considerably frorn the theoretical assumptions. Furthermore, in the case of long reed channels, the results obtained for the vena contracta factor are in significant disagreement with those predicted by theory. For short channels, the assumption of constant vena contracta was found to be valid for only 40% of the duty cycle. (c) 2007 Acoustical Society of America.

Combined R-matrix eigenstate basis set and finite-difference propagation method for the time-dependent Schrodinger equation: The one-electron case

In this work we present the theoretical framework for the solution of the time-dependent Schrödinger equation (TDSE) of atomic and molecular systems under strong electromagnetic fields with the configuration space of the electron’s coordinates separated over two regions; that is, regions I and II. In region I the solution of the TDSE is obtained by an R-matrix basis set representation of the time-dependent wave function. In region II a grid representation of the wave function is considered and propagation in space and time is obtained through the finite-difference method. With this, a combination of basis set and grid methods is put forward for tackling multiregion time-dependent problems. In both regions, a high-order explicit scheme is employed for the time propagation. While, in a purely hydrogenic system no approximation is involved due to this separation, in multielectron systems the validity and the usefulness of the present method relies on the basic assumption of R-matrix theory, namely, that beyond a certain distance (encompassing region I) a single ejected electron is distinguishable from the other electrons of the multielectron system and evolves there (region II) effectively as a one-electron system. The method is developed in detail for single active electron systems and applied to the exemplar case of the hydrogen atom in an intense laser field.

Dynamic pricing models for used products in remanufacturing with lost-sales and uncertain quality

In this paper, we investigate the remanufacturing problem of pricing single-class used products (cores) in the face of random price-dependent returns and random demand. Specifically, we propose a dynamic pricing policy for the cores and then model the problem as a continuous-time Markov decision process. Our models are designed to address three objectives: finite horizon total cost minimization, infinite horizon discounted cost, and average cost minimization. Besides proving optimal policy uniqueness and establishing monotonicity results for the infinite horizon problem, we also characterize the structures of the optimal policies, which can greatly simplify the computational procedure. Finally, we use computational examples to assess the impacts of specific parameters on optimal price and reveal the benefits of a dynamic pricing policy. © 2013 Elsevier B.V. All rights reserved.

A method to reconstruct patient-specific proximal femur surface models from planar pre-operative radiographs

Three-dimensional reconstruction from volumetric medical images (e.g. CT, MRI) is a well-established technology used in patient-specific modelling. However, there are many cases where only 2D (planar) images may be available, e.g. if radiation dose must be limited or if retrospective data is being used from periods when 3D data was not available. This study aims to address such cases by proposing an automated method to create 3D surface models from planar radiographs. The method consists of (i) contour extraction from the radiograph using an Active Contour (Snake) algorithm, (ii) selection of a closest matching 3D model from a library of generic models, and (iii) warping the selected generic model to improve correlation with the extracted contour.

This method proved to be fully automated, rapid and robust on a given set of radiographs. Measured mean surface distance error values were low when comparing models reconstructed from matching pairs of CT scans and planar X-rays (2.57–3.74 mm) and within ranges of similar studies. Benefits of the method are that it requires a single radiographic image to perform the surface reconstruction task and it is fully automated. Mechanical simulations of loaded bone with different levels of reconstruction accuracy showed that an error in predicted strain fields grows proportionally to the error level in geometric precision. In conclusion, models generated by the proposed technique are deemed acceptable to perform realistic patient-specific simulations when 3D data sources are unavailable.

Performance models of storage contention in cloud environments

We propose simple models to predict the performance degradation of disk requests due to storage device contention in consolidated virtualized environments. Model parameters can be deduced from measurements obtained inside Virtual Machines (VMs) from a system where a single VM accesses a remote storage server. The parameterized model can then be used to predict the effect of storage contention when multiple VMs are consolidated on the same server. We first propose a trace-driven approach that evaluates a queueing network with fair share scheduling using simulation. The model parameters consider Virtual Machine Monitor level disk access optimizations and rely on a calibration technique. We further present a measurement-based approach that allows a distinct characterization of read/write performance attributes. In particular, we define simple linear prediction models for I/O request mean response times, throughputs and read/write mixes, as well as a simulation model for predicting response time distributions. We found our models to be effective in predicting such quantities across a range of synthetic and emulated application workloads. 

Domain annihilation due to temperature and thickness gradients in single-crystal BaTiO3

The manner in which 90? ferroelectric-ferroelastic domains respond to changes in temperature has been mapped in BaTiO3 single crystals using atomic force microscopy. Domain periodicity remains unaltered until approximately 2 ? C below the Curie temperature (TC ), whereupon domains coarsened dramatically. This behavior was successfully rationalized by considering the temperature dependence of the parameters associated with standard models of ferroelastic domain formation. However, while successful in describing the expected radical increase in equilibrium period with temperature, the model did not predict the unusual mechanism by which domain coarsening occurred; this was not continuous at a local level but instead involved discrete domain annihilation events. Subsequent insights from a combination of free energy analysis for the system and further experimental data from an analogous situation, in which domain period increases with increasing crystal thickness, suggested that domain annihilation is inevitable whenever a component of the relevant gradient that affects domain period is orientated parallel to the domain walls. Consistent with this thesis, we note that, for the observations presented herein, the thermal gradient possessed a significant component parallel to the domain walls. We suggest that domain annihilation is a general feature of domain structures in gradient fields.