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.

Monte Carlo simulation of cluster growth in surface defects induced by the tip of a scanning tunnelling microscope

First steps are taken to model the electrochemical deposition of metals in nanometer-sized cavities. In the present work, the electrochemical deposition of Cu atoms in nanometer-sized holes dug on Au(111) is investigated through Monte Carlo simulations using the embedded atom method to represent particle interactions. By sweeping the chemical potential of Cu, a cluster is allowed to grow within the hole rising four atomic layers above the surface. Its lateral extension remains confined to the area defined by the borders of the original defect. (C) 2004 Elsevier B.V. All rights reserved.

The impact of virtual simulation in palliative radiotherapy for non-small-cell lung cancer.

Background and purpose: Radiotherapy is widely used to palliate local symptoms in non-small-cell lung cancer. Using conventional X-ray simulation, it is often difficult to accurately localize the extent of the tumour. We report a randomized, double blind trial comparing target localization with conventional and virtual simulation.Methods: Eighty-six patients underwent both conventional and virtual simulation. The conventional simulator films were compared with digitally reconstructed radiographs (DRRs) produced from the computed tomography (CT) data. The treatment fields defined by the clinicians using each modality were compared in terms of field area, position and the implications for target coverage.Results: Comparing fields defined by each study arm, there was a major mis-match in coverage between fields in 66.2% of cases, and a complete match in only 5.2% of cases. In 82.4% of cases, conventional simulator fields were larger (mean 24.5+/-5.1% (95% confidence interval)) than CT-localized fields, potentially contributing to a mean target under-coverage of 16.4+/-3.5% and normal tissue over-coverage of 25.4+/-4.2%.Conclusions: CT localization and virtual simulation allow more accurate definition of the target volume. This could enable a reduction in geographical misses, while also reducing treatment-related toxicity.

Electron collisions with iron peak elements: a computational grand challenge

This paper discusses the calculation of electron impact collision strengths and effective collision strengths for iron peak elements of importance in the analysis of many astronomical and laboratory spectra. It commences with a brief overview of R-matrix theory which is the basis of computer programs which have been widely used to calculate the relevant atomic data used in this analysis. A summary is then given of calculations carried out over the last 20 y for electron collisions with Fe II. The grand challenge, represented by the calculation of accurate collision strengths and effective collision strengths for this ion, is then discussed. A new parallel R-matrix program PRMAT, which is being developed to meet this challenge, is then described and results of recent calculations, using this program to determine optically forbidden transitions in e- – Ni IV on a Cray T3E-1200 parallel supercomputer, are presented. The implications of this e- – Ni IV calculation for the determination of accurate data from an isoelectronic e- – Fe II calculation are discussed and finally some future directions of research are reviewed.

Dynamical simulation of inelastic quantum transport

A method for correlated quantum electron-ion dynamics is combined with a method for electronic open boundaries to simulate in real time the heating, and eventual equilibration at an elevated vibrational energy, of a quantum ion under current flow in an atomic wire, together with the response of the current to the ionic heating. The method can also be used to extract inelastic current voltage corrections under steady-state conditions. However, in its present form the open-boundary method contains an approximation that limits the resolution of current-voltage features. The results of the simulations are tested against analytical results from scattering theory. Directions for the improvement of the method are summarized at the end.

Ion Association in [bmim][PF6]/Naphthalene Mixtures: An Experimental and Computational Study

Mixtures of room temperature ionic liquids (IL) with neutral organic molecules provide a valuable testing ground to investigate the interplay of the ionic and molecular-dipolar state in dense Coulomb systems at near ambient conditions. In the present study, the viscosity eta and the ionic conductivity a of 1-n-butyl-3-methylimidazolium hexafluorophosphate ([bmim][PF6])/naphthalene mixtures at T = 80 degrees C have been measured at 10 stoichiometries spanning the composition range from pure naphthalene to pure [bmim][PF6]. The viscosity grows nearly monotonically with increasing IL mole fraction (x), whereas the conductivity per ion displays a clear peak at x approximate to 15%. The origin of this maximum has been investigated using molecular dynamics simulations based on a classical force field. Snapshots of the simulated samples show that the conductivity maximum is due to the gradual transition in the IL component from an ionic state at high x to a dipolar fluid made of neutral ion pairs at low x. At concentrations x <0.20 the ion pairs condense into molecular-thin filaments bound by dipolar forces and extending in between nanometric droplets of IL. These results are confirmed and complemented by the computation of dynamic and transport properties in [bmim][PF6]/naphthalene mixtures at low IL concentration.

Spectral simulation of laser ablated magnesium plasmas

We describe a collisional-radiative equilibrium model for predicting the optical emission spectrum of low-temperature magnesium plasmas, specifically those created by laser ablation. In the model, levels are populated by a balance of collisional and radiative rates. We include Stark widths of lines and trapping of radiation in the calculations. By use of this model we discuss various issues of importance in spectral analysis of laser ablated plasma plumes, such as the partial local thermodynamic equilibrium approximation, line trapping and time dependence.

A computational model for the cooling phase of injection moulding

This paper discusses the approaches and techniques used to build a realistic numerical model to analyse the cooling phase of the injection moulding process. The procedures employed to select an appropriate mesh and the boundary and initial conditions for the problem are discussed and justified. The final model is validated using direct comparisons with experimental results generated in an earlier study. The model is shown to be a useful tool for further studies aimed at optimising the cooling phase of the injection moulding process. Using the numerical model provides additional information relating to changes in conditions throughout the process, which otherwise could not be deduced or assessed experimentally. These results, and other benefits related to the use of the model, are also discussed in the paper. © 2007 Elsevier B.V. All rights reserved.