Abstract: We have investigated the plasma processes in an atmospheric pressure glow discharge (APGD) in He used for analytical spectrometry by means of fluid and Monte Carlo (MC) simulations. Typical results include the potential and electric field distributions in the plasma, the density profiles of the various plasma species throughout the discharge, the mean electron energy, as well as the rates of the various collision processes in the plasma, and the relative importance of the different production and loss rates for the various species. The similarities and differences with low-pressure glow discharges are discussed. The main differences are a very small cathode dark space region and a large positive column as well as the dominant role of molecular ions. Some characteristic features of the APGD, such as the occurrence of the different spatial zones in the discharge, are illustrated, with links to experimental observations.

Abstract: Both even- and odd-numbered neutral carbon clusters Cn (n=210) are systematically studied using the energy minimization method and the modified Brenner potential for the carbon-carbon interactions. Many stable configurations were found, and several new isomers are predicted. For the lowest energy stable configurations we obtained their binding energies and bond lengths. We found that for n5 the linear isomer is the most stable one while for n>5 the monocyclic isomer becomes the most stable. The latter was found to be regular for all studied clusters. The dependence of the binding energy for linear and cyclic clusters versus the cluster size n (n=210) is found to be in good agreement with several previous calculations, in particular with ab initio calculations as well as with experimental data for n=25.

Abstract: Using reactive molecular dynamics simulations, the melting behavior of nickelcarbon nanoclusters is examined. The phase diagrams of icosahedral and Wulff polyhedron clusters are determined using both the Lindemann index and the potential energy. Formulae are derived for calculating the equilibrium constants and the solid and liquid fractions during a phase transition, allowing more rational determination of the melting temperature with respect to the arbitrary Lindemann value. These results give more insight into the properties of nickelcarbon nanoclusters in general and can specifically be very useful for a better understanding of the synthesis of carbon nanotubes using the catalytic chemical vapor deposition method.

Abstract: Five geothermal waters from the Mons area (southern Belgium) have been studied: one natural hot spring at Stambruges, one stagnant warm water from the inclined tunnels at Baudour, and three thermal waters from the drillholes at St. Ghislain, Ghlin and Douvrain, originating from the carbonate/anhydrite-bearing Visean strata, at depths of ca. 2600, 1550 and 1300 m, respectively. Multielement chemical analysis of the filtered water and its suspended matter > 0.4 μm) was carried out by instrumental neutron activation. Temperature in depth, calculated using the silica (chalcedony) chemical geothermometer, ranged from 75 to 88°C, in good agreement with experimentally determined values. Na/K and Na/K/Ca geothermometers yieilded erratic results, as expected from the geological environment in the aquifer. From the analytical data it can be calculated that the thermal waters of St. Ghislain, Ghlin and Douvrain are not only saturated with respect to chalcedony, but also to anhydrite, calcite, fluorite, barite, strontianite, and possibly zinc silicate, iron (III) hydroxide or siderite, albite, microcline, gibbsite and kaolinite. They are oversaturated with respect to muscovite. Data are also presented for the other thermal waters, and a cold spring water (Claire Fontaine, Stambruges). The similar trace-element composition of the thermal waters can be explained by percolation of the water in the same distant recharge zone, from where it descends, becomes heated at depth and rises along collapse breccia, and locally (Baudour, Stambruges) along fissures. The uptake of higher amounts of Ca, Mg, Sr and sulfate in St. Ghislain and Ghlin, as compared to Douvrain and Baudour is correlated with the boundary between the non-dissolved and dissolved evaporitic zones. This boundary is situated between St. Ghislain and Douvrain, and is roughly parallel with the direction of the groundwater flow (WNW).

Abstract: A two-dimensional model is presented that describes the behavior of argon metastable atoms, copper atoms, and copper ions in an argon direct. current glow discharge, in the standard cell of the VG9000 glow discharge mass spectrometer for analyzing flat samples. The model is combined with a previously developed model for the electrons, argon ions, and atoms in the same cell to obtain an overall picture of the glow discharge, The results of the present model comprise the number densities of the described plasma species, the relative contributions of different production and loss processes for the argon metastable atoms, the thermalization profile of the sputtered copper atoms, the relative importance of the different ionization mechanisms for the copper atoms, the ionization degree of copper, the copper ion-to-argon ion density ratio, and the relative roles of copper ions, argon ions, and atoms in the sputtering process. All these quantities are calculated for a range of voltages and pressures, Moreover, since the sticking coefficient of copper atoms on solid surfaces is not well-known in the literature, the influence of this parameter on the results is briefly discussed.

Abstract: The plasma behavior in a parallel-plate dielectric barrier discharge (DBD) is simulated by a two-dimensional particle-in-cell/Monte Carlo collision model, comparing for the first time an unpacked (empty) DBD with a packed bed DBD, i.e., a DBD filled with dielectric spheres in the gas gap. The calculations are performed in air, at atmospheric pressure. The discharge is powered by a pulse with a voltage amplitude of −20 kV. When comparing the packed and unpacked DBD reactors with the same dielectric barriers, it is clear that the presence of the dielectric packing leads to a transition in discharge behavior from a combination of negative streamers and unlimited surface streamers on the bottom dielectric surface to a combination of predominant positive streamers and limited surface discharges on the dielectric surfaces of the beads and plates. Furthermore, in the packed bed DBD, the electric field is locally enhanced inside the dielectric material, near the contact points between the beads and the plates, and therefore also in the plasma between the packing beads and between a bead and the dielectric wall, leading to values of $4\times {10}

Abstract: High density polyethylene surfaces were exposed to the atmospheric post-discharge of a radiofrequency plasma torch supplied in helium and oxygen. Dynamic water contact angle measurements were performed to evaluate changes in surface hydrophilicity and angle resolved x-ray photoelectron spectroscopy was carried out to identify the functional groups responsible for wettability changes and to study their subsurface depth profiles, up to 9 nm in depth. The reactions leading to the formation of CO, C = O and OC = O groups were simulated by molecular dynamics. These simulations demonstrate that impinging oxygen atoms do not react immediately upon impact but rather remain at or close to the surface before eventually reacting. The simulations also explain the release of gaseous species in the ambient environment as well as the ejection of low molecular weight oxidized materials from the surface.

Abstract: Monte Carlo (MC) methods have a long-standing history as partners of molecular dynamics (MD) to simulate the evolution of materials at the atomic scale. Among these techniques, the uniform-acceptance force-bias Monte Carlo (UFMC) method [ G. Dereli Mol. Simul. 8 351 (1992)] has recently attracted attention [ M. Timonova et al. Phys. Rev. B 81 144107 (2010)] thanks to its apparent capacity of being able to simulate physical processes in a reduced number of iterations compared to classical MD methods. The origin of this efficiency remains, however, unclear. In this work we derive a UFMC method starting from basic thermodynamic principles, which leads to an intuitive and unambiguous formalism. The approach includes a statistically relevant time step per Monte Carlo iteration, showing a significant speed-up compared to MD simulations. This time-stamped force-bias Monte Carlo (tfMC) formalism is tested on both simple one-dimensional and three-dimensional systems. Both test-cases give excellent results in agreement with analytical solutions and literature reports. The inclusion of a time scale, the simplicity of the method, and the enhancement of the time step compared to classical MD methods make this method very appealing for studying the dynamics of many-particle systems.

Abstract: In order to increase the understanding of the underlying processes in organic solvent nanofiltration (OSN), a study has been undertaken aimed at clarifying the solvent flux behaviour of ceramic nanofiltration and ultrafiltration membranes. Ceramic membranes were chosen for their non-swelling character. Pure water and a variation of 11 different organic solvents were measured on a series of different ceramic membranes with pore-size diameters ranging from 0.9 nm up to 100 nm. To avoid any historical effects, each flux measurement was carried out on a new membrane. The flux results were analysed in a phenomenological way, and a common very simple linear relationship was observed between the product of flux and viscosity of the solvent, and the total Hansen solubility parameter of the solvent. The linear relationship was found for all membranes, independent of the membrane pore size and the membrane material. The slope of the linear relationship was found to depend exponentially on the pore-size diameter and on the polarity of the membrane surface. This result emphasizes the importance of viscosity in the solvent transport, but also of the polarity difference between membrane surface and solvent. The very simple flux model deduced, allows a straightforward prediction of the flux of any solvent or solvent mixture, once the water flux of the membrane is known. At the high pore-size end, the phenomenological model naturally transforms into the viscous-flow or pore-flow behaviour as required. A tentative physical explanation of the model takes into account the presence and extension of a water layer adsorbed to the total pore surface of these membranes. This work also shows that the water flux of a hydrophilic membrane gives a good indication of its molecular weight cut-off (MWCO), and therefore of its separation performance in water. (C)0 2013 Elsevier B.V. All rights reserved.