Abstract: The interaction between cold atmospheric pressure plasma and liquids is receiving increasing attention for various applications. In particular, the use of plasma-treated liquids (PTL) for biomedical applications is of growing importance, in particular for sterilization and cancer treatment. However, insight into the
underlying mechanisms of plasma–liquid interactions is still scarce. Here, we present a 2D fluid dynamics model for the interaction between a plasma jet and liquid water. Our results indicate that the formed reactive species originate from either the gas phase (with further solvation) or are formed at the liquid interface. A clear increase in the aqueous density of H2O2, HNO2/NO2- and NO3-
is observed as a function of time, while the densities of O3, HO2/O2- and ONOOH/ONOO- are found to quickly reach a maximum due to chemical reactions in solution. The trends observed in our model correlate well with experimental observations from the literature.

Abstract: The interaction of silicene, the silicon counterpart of graphene, with (0001) ZnS surfaces is investigated theoretically, using first-principles simulations. The charge transfer occurring at the silicene/(0001) ZnS interface leads to the opening of an indirect energy band gap of about 0.7 eV in silicene. Remarkably, the nature (indirect or direct) and magnitude of the energy band gap of silicene can be controlled by an external electric field: the energy gap is predicted to become direct for electric fields larger than about 0.5 V angstrom(-1), and the direct energy gap decreases approximately linearly with the applied electric field. The predicted electric field tunable energy band gap of the silicene/(0001) ZnS interface is very promising for its potential use in nanoelectronic devices.

Abstract: The laser ablation-induced plasma was used as a composition-con trolled source for ion implantation in Si crystals. Then, laser ablation in combination with inductively coupled plasma mass spectrometry was used for the elemental depth profiling of the implanted samples. Monte Carlo simulations permitted us to conclude that a depth resolution of tens of nm would be necessary to define the shape of the implantation profiles, as is obtained using XPS and RBS, whereas a hundred nm depth resolution is sufficient to determine the total implanted dose. The detection power of LA-ICP-MS would routinely allow rapid analytical control on the trace level implanted dose. Nevertheless, this technique is limited in terms of depth profiling resolution due to pulse mixing and signal tailing induced during the aerosol transport. Raw signal processing procedures were developed for the minimization of shapeline dispersion, deconvolution of pulse mixing and more appropriate assessment of the implanted profiles. Shapeline dispersion could be corrected for by determining the signal waning constant and implementing this information for a non-affine alibi transformation of the LA-ICP-MS signal traces. Pulse mixing deconvolution was attained with an algorithm that considered accumulated signal intensity due to pulse-on-pulse stacking, i.e., the latest pulse on top of all antecedent individual pulses' exponential tails proportionally.

Abstract: The low-pressure hollow cathode discharge made of a hollow circular tube and an anode is a type of simple structure discharge system. In particular, under the arc discharge mode, hollow cathodes have high plasma density and energy density with a wide range of adaptability of pressure and current. Low-pressure hollow cathode arc (HCA) discharges have been widely used as plasma sources in various fields such as manufacturing, vacuum welding, and aerospace since the 1960s. Despite the early experimental and applied researches on low-pressure HCA discharges, the basic theoretical study was relatively lagged much behind, resulting in many unanswered questions, such as the optimal discharge operating parameters, the power deposition inside the cathode, the causes of plasma instability, and how to effectively reduce cathode erosion and so on. Due to the special discharge structure of the hollow cathode, it is difficult to make an accurate experimental diagnosis, so a reasonable numerical simulation is an effective study method. However, up to now, there is still a lack of complete and effective numerical models which can evaluate various physical fields in the low-pressure hollow cathode discharges. To address the above problems and difficulties, a comprehensive and self-consistent 2D multi-physical coupling numerical model based on a commercial program of finite element method, the COMSOL Multiphysics, was provided in this paper. The model involves plasma transport, arc flow and heat transfer, and cathode thermal equilibrium, and can consider the effect of an applied magnetic field. The processes of secondary electron emission, thermal-field electron emission, ions and backflow high-energy electrons bombardment, and thermal radiation from the cathode surface are considered in the cathode thermal equilibrium process. Based on the above background, this paper works from the following aspects： In Chapter 1, the basic concepts of low-pressure HCA discharge including the hollow cathode effect, the basic characteristics, and operation modes were introduced firstly; Secondly, the application fields, development history, and overseas and domestic research status of hollow cathode discharge were reviewed; finally, the problems were presented and the research background was explained, and the research purpose of this paper was clarified. In Chapter 2, a complete and self-consistent numerical model of low-pressure hollow cathode discharge was proposed based on the fundamental theory and assumptions, and the set of control equations and boundary conditions in the model were elaborated. In addition, the electron energy distribution function, the collision processes, the solving tools of this model, and calculation schemes were introduced in detail. Finally, a validation example was given to test the rationality and applicability of the numerical model. In Chapter 3, the fundamental plasma properties of low-pressure hollow cathode arcs were investigated. Firstly, the ion Joule heating effect was studied. The results showed that the temperature distributions of the arc and cathode are only able to approach the experimental measurements after considering the ion Joule heating, which shows that the Joule heating of ions is crucial for the heating of the arc plasma. Secondly, by comparing the radial distribution of electron and ion density inside the cathode, the structure of the cathode sheath could be simulated well using this model. Finally, it was shown that the thermal radiation from the cathode surface is an important cooling mechanism of the cathode and only under higher surface emissivity can balance the larger heat flow given by the plasma to the cathode, and the temperature distribution of the cathode shows a non-monotonic increasing trend and is consistent with the profile of experimental measurement so that the so-called active zone is formed. In Chapter 4, the power deposition in the low-pressure HCA was studied in simulation. Two main aspects were considered: the power deposition into particles (both electrons and heavy particles) and the power deposition onto the cathode. It was found that the deposited power into particles increases with the rise of discharge current, but there is no effect on the total power deposition onto the cathode. In high-density plasmas, Coulomb collisions between electrons and ions also become very important, especially since a portion of the deposition energy on heavy particles comes mainly from the energy transfer from electrons to ions. It was also found that regardless of external parameters, half of the power deposition onto the cathode always comes from the particle contribution, while the other half is the net contribution of heat transfer and cathode radiation. The HCA model also allows the simulation of multiple discharge modes for low-pressure HCA discharges over a wide range of gas flow rates. It was also shown that the discharge operating conditions and the external magnetic field can change the distribution of the particle flow on the cathode wall. In Chapter 5, the ion sputtering erosion process on the cathode was simulated by coupling the HCA numerical model with the moving grid technique. The results showed that the ion sputtering erosion on the cathode depends on the ion flux and the plasma potential near the cathode wall and that their distribution and magnitude jointly determine the erosion morphology of the cathode. It was also found that the location of the most severe erosion on the cathode is located in the region of the densest ion flux on the cathode wall, rather than in the longitudinal correspondence with the central region of the internal positive column (IPC). The external magnetic fields can mitigate the cathode erosion and reduce the erosion depth, but stronger magnetic fields lead to a concentration of current density at the cathode tip, which can enhance erosion slightly at the cathode outlet end. Finally, the conclusions and innovation highlights were summarized, and prospects for future work were discussed.

Abstract: The mechanism of interaction of cold nonequilibrium plasma jets with mammalian cells in physiologic liquid is reported. The major biological active species produced by an argon RF plasma jet responsible for cell viability reduction are analyzed by experimental results obtained through physical, biological, and chemical diagnostics. This is complemented with chemical kinetics modeling of the plasma source to assess the dominant reactive gas phase species. Different plasma chemistries are obtained by changing the feed gas composition of the cold argon based RF plasma jet from argon, humidified argon (0.27%), to argon/oxygen (1%) and argon/air (1%) at constant power. A minimal consensus physiologic liquid was used, providing isotonic and isohydric conditions and nutrients but is devoid of scavengers or serum constituents. While argon and humidified argon plasma led to the creation of hydrogen peroxide dominated action on the mammalian cells, argonoxygen and argonair plasma created a very different biological action and was characterized by trace amounts of hydrogen peroxide only. In particular, for the argonoxygen (1%), the authors observed a strong negative effect on mammalian cell proliferation and metabolism. This effect was distance dependent and showed a half life time of 30 min in a scavenger free physiologic buffer. Neither catalase and mannitol nor superoxide dismutase could rescue the cell proliferation rate. The strong distance dependency of the effect as well as the low water solubility rules out a major role for ozone and singlet oxygen but suggests a dominant role of atomic oxygen. Experimental results suggest that O reacts with chloride, yielding Cl2 − or ClO−. These chlorine species have a limited lifetime under physiologic conditions and therefore show a strong time dependent biological activity. The outcomes are compared with an argon MHz plasma jet (kinpen) to assess the differences between these (at least seemingly) similar plasma sources.

Abstract: The modelling of a gliding arc discharge (GAD) is studied by means of the quasineutral (QN) plasma modelling approach. The model is first evaluated for reliability and proper description of a gliding arc discharge at atmospheric pressure, by comparing with a more elaborate non-quasineutral (NQN) plasma model in two different geometries – a 2D axisymmetric and a Cartesian geometry. The NQN model is considered as a reference, since it provides a continuous self-consistent plasma description, including the near electrode regions. In general, the results of the QN model agree very well with those obtained from the NQN model. The small differences between both models are attributed to the approximations in the derivation of the QN model. The use of the QN model provides a substantial reduction of the computation time compared to the NQN model, which is crucial for the development of more complex models in three dimensions or with complicated chemistries. The latter is illustrated for (i) a reverse vortex flow(RVF) GAD in argon, and (ii) a GAD in CO2. The RVF discharge is modelled in three dimensions and the effect of the turbulent heat transport on the plasma and gas characteristics is
discussed. The GAD model in CO2 is in a 1D geometry with axial symmetry and provides results for the time evolution of the electron, gas and vibrational temperature of CO2, as well as for the molar fractions of the different species.

Abstract: The nucleation and growth of titanium dioxide (TiO2) and zinc oxide (ZnO) thin films on Fe2O3 (hematite), Al2O3 (á-alumina) and SiO2 (á-quartz) are studied by molecular dynamics simulations. The results show the formation of a strong interface region between the substrate and the film in the six systems studied here. A combination of polycrystalline and amorphous phases are observed in the TiO2 films grown on the three substrates. ZnO deposition on the Fe2O3 and Al2O3 crystals yields a monocrystalline film growth. The ZnO film deposited on the SiO2 crystal exhibits less crystallinity. The simulation results are compared with experimental results available in the literature.

Abstract: The origin of the ferroelectric polarization switching in orthorhombic HfO2 has been investigated by first principles calculations. The phenomenon can be regarded as being the coordinated displacement of four O ions in the orthorhombic unit cell, which can lead to a saturated polarization as high as 53 mu C/cm(2). We show the correlation between the computed polarization reversal barrier and the experimental coercive fields. (C) 2014 AIP Publishing LLC.

Abstract: The peculiarities of sputtering processes at 0.5-5 keV Ne grazing ion bombardment of Si(001) and SiC(001) surfaces and their possible application for the surface modification have been studied by computer simulation. Sputtering yields in the primary knock-on recoil atoms regime versus the initial energy of incident ions (E(0) = 0.5-5 keV) and angle of incidence (psi = 0-30 degrees) counted from a target surface have been calculated. Comparative studies of layer-by-layer sputtering for Si(001) and SiC(001) surfaces versus the initial energy of incident ions as well as an effective sputtering and sputtering threshold are discussed.

Abstract: The phase transition processes induced by ultrashort, 100 fs pulsed laser irradiation of Au, Cu, and Ni are studied by means of a combined atomistic-continuum approach. A moderately low absorbed laser fluence range, from 200 to 600 J/m2 is considered to study phase transitions by means of a local and a nonlocal order parameter. At low laser fluences, the occurrence of layer-by-layer evaporation has been observed, which suggests a direct solid to vapor transition. The calculated amount of molten material remains very limited under the conditions studied, especially for Ni. Therefore, our results show that a kinetic equation that describes a direct solid to vapor transition might be the best approach to model laser-induced phase transitions by continuum models. Furthermore, the results provide more insight into the applicability of analytical superheating theories that were implemented in continuum models and help the understanding of nonequilibrium phase transitions.

Abstract: The physical processes in an Ar/O2 magnetron discharge used for the reactive sputter deposition of TiOx thin films were simulated with a 2d3v particle-in-cell/Monte Carlo collisions (PIC/MCC) model. The plasma species taken into account are electrons, Ar+ ions, fast Arf atoms, metastable Arm* atoms, Ti+ ions, Ti atoms, O+ ions, O2+ ions, O− ions and O atoms. This model accounts for plasmatarget interactions, such as secondary electron emission and target sputtering, and the effects of target poisoning. Furthermore, the deposition process is described by an analytical surface model. The influence of the O2/Ar gas ratio on the plasma potential and on the species densities and fluxes is investigated. Among others, it is shown that a higher O2 pressure causes the region of positive plasma potential and the O− density to be more spread, and the latter to decrease. On the other hand, the deposition rates of Ti and O are not much affected by the O2/Ar proportion. Indeed, the predicted stoichiometry of the deposited TiOx film approaches x=2 for nearly all the investigated O2/Ar proportions.

Abstract: The physics of plasmas containing positive and negative ions is discussed with special attention to the recently produced pair-ion plasma containing ions of equal mass and opposite charge. The effects of the density gradient in the direction perpendicular to the ambient magnetic field vector are discussed. The possible presence of electrons is discussed in the context of plasma modes propagating at an angle with respect to the magnetic field vector. It is shown that the electron plasma mode may become a backward mode in the presence of a density gradient, and this behavior may be controlled either by the electron number density or the mode number in the perpendicular direction. In plasmas with hot electrons an instability may develop, driven by the combination of electron collisions and the density gradient, and in the regime of a sound ions' response. In the case of a pure pair-ion plasma, for lower frequencies and for parameters close to those used in the recent experiments, the perturbed ions may feel the effects of the magnetic field. In this case the plasma mode also becomes backward, resembling features of an experimentally observed but yet unexplained backward mode. (C) 2008 American Institute of Physics.

Abstract: The plasma-based dry reforming in a dielectric barrier discharge (DBD) reactor is important to achieve sustainable goals, but many challenges remain. For example, the conversion and energy yield of DBD reactors are relatively low, and the catalysts or packing materials used in existing studies cannot improve them, possibly due to the unsuitable properties and structures of catalysts or packing materials for plasma processes. In order to study the effect of catalyst structure on plasma-based dry reforming, a controllable synthesis of the catalyst supports or templates was explored. In Chapter 2, an initially immiscible synthesis method was proposed to synthesize uniform silica spheres, which can replace the organic solvent-based Stöber method to successfully synthesize silica particles with the same size ranges as the original Stöber process without addition of organic solvents. Using the silica spheres as templates, 3D porous Cu and CuO catalysts with different pore sizes were synthesized in Chapter 3 to study the effect of catalyst pore size on the plasma-catalytic dry reforming. In most cases, the smaller the pore size, the higher the conversion of CH4 and CO2 due to the reaction of radicals and ions formed in the plasma. An exception are the samples synthesized from 1 μm silica, which show better performance due to the electric field enhancement for pore sizes close to the Debye length. Besides the pore size, the particle diameter of the catalyst or packing is also one of the important factors affecting the interaction between plasma and catalyst. In Chapter 4, SiO2 spheres (with or without supported metal) were used to study the effect of different support particle sizes on plasma-based dry reforming. We found that a uniform SiO2 packing improves the conversion of plasma-based dry reforming. The conversion of plasma-based dry reforming first increases and then decreases with increasing particle size, due to the balance between the promoting and hindering effect of the particle filling on the plasma discharge. Chapter 5 is to improve the design of the DBD reactor itself, in order to try to increase its low energy yield. Some stainless steel rings were put over the inner electrode rod of the DBD reactor. The presence of rings increases the local electric field, the displaced charges and the discharge fraction, and also makes the discharge more stable and with more uniform intensity. The placement of the rings improves the performance of the reactor at 30 W supplied power.

Abstract: The possibility of confinement effects inside a carbon nanotube provides new application opportunities, e.g., growth of novel carbon nanostructures. However, the understanding the precise role of catalystfeedstock in the nanostructure synthesis is still elusive. In our simulation-based study, we investigate the Ni-catalyzed growth mechanism of encapsulated carbon nanostructures, viz. double-wall carbon nanotube and graphene nanoribbon, from carbon and hydrocarbon growth precursors, respectively. Specifically, we find that the tube and ribbon growth is determined by a catalyst-vs-feedstock competition effect. We compare our results, i.e., growth mechanism and structure morphology with all available theoretical and experimental data. Our calculations show that all encapsulated nanostructures contain metal (catalyst) atoms and such structures are less stable than their pure counterparts. Therefore, we study the purification mechanism of these structures. In general, this study opens a possible route to the controllable synthesis of tubular and planar carbon nanostructures for today’s nanotechnology.

Abstract: The presence of iodide ions during the growth of gold nanorods strongly affects the shape of the final products, which is proposed to be due to selective iodide adsorption on certain crystallographic facets. Therefore, a detailed structural and morphological characterization of the starting rods is crucial toward understanding this effect. Electron tomography is used to determine the crystallographic indices of the lateral facets of gold nanorods, as well as those present at the tips. Based on this information, density functional theory calculations are used to determine the surface and interface energies of the observed facets and provide insight into the relationship between the amount of iodide ions in the growth solution and the final morphology of anisotropic gold nanoparticles.