Ab initio simulation of charged slabs at constant chemical potential

We present a practical scheme for performing ab initio supercell calculations of charged slabs at constant electron chemical potential mu, rather than at constant number of electrons N-e. To this end, we define the chemical potential relative to a plane (or "reference electrode") at a finite distance from the slab (the distance should reflect the particular geometry of the situation being modeled). To avoid a net charge in the supercell, and thus make possible a standard supercell calculation, we restore the electroneutrality of the periodically repeated unit by means of a compensating charge, whose contribution to the total energy and potential is subtracted afterwards. The "constant mu" mode enables one to perform supercell calculation on slabs, where the slab is kept at a fixed potential relative to the reference electrode. We expect this to be useful in modeling many experimental situations, especially in electro-chemistry. (C) 2001 American Institute of Physics.

Angular distribution for the elastic scattering of electrons from Ar+(3s^2 3p^5 ^2P) above the first inelastic threshold

The measured angular differential cross section (DCS) for the elastic scattering of electrons from Ar+(3s2 3p5 2P) at the collision energy of 16 eV is presented. By solving the Hartree-Fock equations, we calculate the corresponding theoretical DCS including the coupling between the orbital angular momenta and spin of the incident electron and those of the target ion and also relaxation effects. Since the collision energy is above one inelastic threshold for the transition 3s2 3p5 2P–3s 3p6 2S, we consider the effects on the DCS of inelastic absorption processes and elastic resonances. The measurements deviate significantly from the Rutherford cross section over the full angular range observed, especially in the region of a deep minimum centered at approximately 75°. Our theory and an uncoupled, unrelaxed method using a local, spherically symmetric potential by Manson [Phys. Rev. 182, 97 (1969)] both reproduce the overall shape of the measured DCS, although the coupled Hartree-Fock approach describes the depth of the minimum more accurately. The minimum is shallower in the present theory owing to our lower average value for the d-wave non-Coulomb phase shift s2, which is due to the high sensitivity of s2 to the different scattering potentials used in the two models. The present measurements and calculations therefore show the importance of including coupling and relaxation effects when accurately modeling electron-ion collisions. The phase shifts obtained by fitting to the measurements are compared with the values of Manson and the present method.

Differential cross section measurements for elastic scattering of electrons from Ar2+ and Xe2+

A crossed-beams energy-loss spectrometer has been used to investigate angular distributions for electron scattering from Ar2+ and Xe2+ ions, at a collision energy of 16 eV. For Ar2+ the measurements are compared with the predictions of a partial waves calculation based on a semi-empirical potential, where it is shown that the interference term governs the position of the observed minimum in the angular distribution.

Elastic scattering of electrons from argon ions

A crossed-beams energy-loss spectrometer has been used to investigate angular distributions for electron scattering from Ar2+ and Ar3+ ions, at a collision energy of 16 eV. Results are compared with the predictions of a partial waves calculation based on a semi-empirical potential, and with the classical Rutherford formula.

Dynamical Equations For The Period Vectors In A Periodic System Under Constant External Stress

The purpose of this paper is to derive the dynamical equations for the period vectors of a periodic system under constant external stress. The explicit starting point is Newton’s second law applied to halves of the system. Later statistics over indistinguishable translated states and forces associated with transport of momentum are applied to the resulting dynamical equations. In the final expressions, the period vectors are driven by the imbalance between internal and external stresses. The internal stress is shown to have both full interaction and kinetic-energy terms.

Plasma boundary sheath in the afterglow of a pulsed inductively coupled RF plasma

The sheath dynamics in the afterglow of a pulsed inductively coupled plasma, operated in hydrogen, is investigated. It is found that the sheath potential does not fully collapse in the early post-discharge. Time resolved measurements of the positive ion flux in a hydrogen plasma, using a mass resolved ion energy analyser, reveal that a constant 2 eV mean ion energy persists for several hundred micro-seconds in the afterglow. The presence of a finite sheath potential is explained by super-elastic collisions between vibrationally excited hydrogen molecules and electrons in the afterglow, leading to an electron temperature of about 0.5 eV. Plasma density decay times measured using both the mass resolved energy analyser and a Langmuir probe are in good agreement. Vibrational temperatures measured using optical emission spectroscopy support the theory of electron heating through super-elastic collisions with vibrationally excited hydrogen molecules. Measurements are also supported by numerical simulations and modelling results.

Evolution of elastic x-ray scattering in laser-shocked warm dense lithium

We have studied the dynamics of warm dense Li with near-elastic x-ray scattering. Li foils were heated and compressed using shock waves driven by 4-ns-long laser pulses. Separate 1-ns-long laser pulses were used to generate a bright source of 2.96 keV Cl Ly-alpha photons for x-ray scattering, and the spectrum of scattered photons was recorded at a scattering angle of 120 degrees using a highly oriented pyrolytic graphite crystal operated in the von Hamos geometry. A variable delay between the heater and backlighter laser beams measured the scattering time evolution. Comparison with radiation-hydrodynamics simulations shows that the plasma is highly coupled during the first several nanoseconds, then relaxes to a moderate coupling state at later times. Near-elastic scattering amplitudes have been successfully simulated using the screened one-component plasma model. Our main finding is that the near-elastic scattering amplitudes are quite sensitive to the mean ionization state Z and by extension to the choice of ionization model in the radiation-hydrodynamics simulations used to predict plasma properties within the shocked Li.

Analytical Elastic Stiffness Model for 3D Woven Orthogonal Interlock Composites

This research presents the development of an analytical model to predict the elastic stiffness performance of orthogonal interlock bound 3D woven composites as a consequence of altering the weaving parameters and constituent material types. The present approach formulates expressions at the micro level with the aim of calculating more representative volume fractions of a group of elements to the layer. The rationale in representing the volume fractions within the unit cell more accurately was to improve the elastic stiffness predictions compared to existing analytical modelling approaches. The models developed in this work show good agreement between experimental data and improvement on existing predicted values by models published in literature.

A heuristic model for the simulation of the deformation of elastic and spongy material for virtual reality applications

This paper presents a practical algorithm for the simulation of interactive deformation in a 3D polygonal mesh model. The algorithm combines the conventional simulation of deformation using a spring-mass-damping model, solved by explicit numerical integration, with a set of heuristics to describe certain features of the transient behaviour, to increase the speed and stability of solution. In particular, this algorithm was designed to be used in the simulation of synthetic environments where it is necessary to model realistically, in real time, the effect on non-rigid surfaces being touched, pushed, pulled or squashed. Such objects can be solid or hollow, and have plastic, elastic or fabric-like properties. The algorithm is presented in an integrated form including collision detection and adaptive refinement so that it may be used in a self-contained way as part of a simulation loop to include human interface devices that capture data and render a realistic stereoscopic image in real time. The algorithm is designed to be used with polygonal mesh models representing complex topology, such as the human anatomy in a virtual-surgery training simulator. The paper evaluates the model behaviour qualitatively and then concludes with some examples of the use of the algorithm.

Dielectric constant characterization using a numerical method for the microstrip ring resonator

A new method of dielectric-constant measurement is developed. The dielectric constant epsilon(r) RF/microwave substrate is extracted by combining the microstrip ring resonator measurement with Ansoft HFSS electromagnetic simulation software. The developed method has two advantages: (i) characterization of dielectric constant versus multiple frequency points, and (ii) compatibility with electronics design automation (EDA) software tools. This characterization method can reduce the design cycle of microwave circuits and devices. (C) 2004 Wiley Periodicals, Inc.

How can cells sense the elasticity of a substrate? An analysis using a cell tensegrity model

A eukaryotic cell attaches and spreads on substrates, whether it is the extracellular matrix naturally produced by the cell itself, or artificial materials, such as tissue-engineered scaffolds. Attachment and spreading require the cell to apply forces in the nN range to the substrate via adhesion sites, and these forces are balanced by the elastic response of the substrate. This mechanical interaction is one determinant of cell morphology and, ultimately, cell phenotype. In this paper we use a finite element model of a cell, with a tensegrity structure to model the cytoskeleton of actin filaments and microtubules, to explore the way cells sense the stiffness of the substrate and thereby adapt to it. To support the computational results, an analytical 1D model is developed for comparison. We find that (i) the tensegrity hypothesis of the cytoskeleton is sufficient to explain the matrix-elasticity sensing, (ii) cell sensitivity is not constant but has a bell-shaped distribution over the physiological matrix-elasticity range, and (iii) the position of the sensitivity peak over the matrix-elasticity range depends on the cytoskeletal structure and in particular on the F-actin organisation. Our model suggests that F-actin reorganisation observed in mesenchymal stem cells (MSCs) in response to change of matrix elasticity is a structural-remodelling process that shifts the sensitivity peak towards the new value of matrix elasticity. This finding discloses a potential regulatory role of scaffold stiffness for cell differentiation.

On visco-elastic modelling of polyethylene terephthalate behaviour during multiaxial elongations slightly over the glass transition temperature

The mechanical response of Polyethylene Terephthalate (PET) in elongation is strongly dependent on temperature, strain and strain rate. Near the glass transition temperature Tg, the stress-strain curve presents a strain softening effect vs strain rate but a strain hardening effect vs strain under conditions of large deformations. The main goal of this work is to propose a viscoelastic model to predict the PET behaviour when subjected to large deformations and to determine the material properties from the experimental data. To represent the non–linear effects, an elastic part depending on the elastic equivalent strain and a non-Newtonian viscous part depending on both viscous equivalent strain rate and cumulated viscous strain are tested. The model parameters can then be accurately obtained trough a comparison with the experimental uniaxial and biaxial tests. The in?uence of the temperature on the viscous part is also modelled and an evaluation of the adiabatic self heating of the specimen is compared to experimental results.