Abstract: The deposition of ultra-thin tungsten films and their related 2D chalcogen compounds on large area dielectric substrates by gas phase reactions is challenging. The lack of nucleation sites complicates the adsorption of W-related precursors and subsequent sulfurization usually requires high temperatures. We propose here a technique in which a thin solid amorphous silicon film is used as reductant for the gas phase precursor WF6 leading to the conversion to metallic W. The selectivity of the W conversion towards the underlying dielectric surfaces is demonstrated. The role of the Si surface preparation, the conversion temperature, and Si thickness on the formation process is investigated. Further, the in situ conversion of the metallic tungsten into thin stoichiometric WS2 is achieved by a cyclic approach based on WF6 and H2S pulses at the moderate temperature of 450 1C, which is much lower than usual oxide sulfurization processes.

Abstract: The generation of high-density plasmas on the surface of porous catalysts is very important for plasma catalysis, as it determines the active surface of the catalyst that is available for the reaction. In this work, we investigate the mechanism of surface and volume plasma streamer formation and propagation near micro-sized pores in dielectric barrier discharges operating in air at atmospheric pressure. A two-dimensional particle-in-cell/ Monte Carlo collision model is used to model the individual kinetic behavior of plasma species. Our calculations indicate that the surface discharge is enhanced on the surface of the catalyst pores compared with the microdischarge inside the catalyst pores. The reason is that the surface ionization wave induces surface charging along the catalyst pore sidewalls, leading to a strong electric field along the pore sidewalls, which in turn further enhances the surface discharge. Therefore, highly concentrated reactive species occur on the surfaces of the catalyst pores, indicating high-density plasmas on the surface of porous catalysts. Indeed, the maximum electron impact excitation and ionization rates occur on the pore surface, indicating the more pronounced production of excited state and electron-ion pairs on the pore surface than inside the pore, which may profoundly affect the plasma catalytic process. Published under license by AIP Publishing.

Abstract: Natural deep eutectic solvents (NADES) are mixtures of naturally derived compounds with a significantly decreased melting point due to the specific interactions among the constituents. NADES have benign properties (low volatility, flammability, toxicity, cost) and tailorable physicochemical properties (by altering the type and molar ratio of constituents), hence they are often considered as a green alternative to common organic solvents. Modeling the relation between their composition and properties is crucial though, both for understanding and predicting their behavior. Several efforts were done to this end, yet this review aims at structuring the present knowledge as an outline for future research. First, we reviewed the key properties of NADES and relate them to their structure based on the available experimental data. Second, we reviewed available modeling methods applicable to NADES. At the molecular level, density functional theory and molecular dynamics allow interpreting density differences and vibrational spectra, and computation of interaction energies. Additionally, properties at the level of the bulk media can be explained and predicted by semi-empirical methods based on ab initio methods (COSMO-RS) and equation of state models (PC-SAFT). Finally, methods based on large datasets are discussed; models based on group contribution methods and machine learning. A combination of bulk media and dataset modeling allows qualitative prediction and interpretation of phase equilibria properties on the one hand, and quantitative prediction of melting point, density, viscosity, surface tension and refractive indices on the other hand. In our view, multiscale modeling, combining the molecular and macroscale methods, will strongly enhance the predictability of NADES properties and their interaction with solutes, yielding truly tailorable solvents to accommodate (bio)chemical reactions.

Abstract: In this paper we employ first-principles calculations to investigate the effect of substitutional Si doping in the amorphous calcium-phosphate (a-HAP) structure on the work of adhesion, integral charge transfer, charge density difference and theoretical tensile strengths between an a-HAP coating and amorphous titanium dioxide (a-TiO2) substrate systemically. Our calculations demonstrate that substitution of a P atom by a Si atom in a-HAP (a-Si-HAP) with the creation of OH-vacancies as charge compensation results in a significant increase of the bonding strength of the coating to the substrate. The work of adhesion of the optimized Si-doped interfaces reaches a value of up to -2.52 J m(-2), which is significantly higher than for the stoichiometric a-HAP/a-TiO2. Charge density difference analysis indicates that the dominant interactions at the interface have significant covalent character, and in particular two Ti-O and three Ca-O bonds are formed for a-Si-HAP/a-TiO2 and one Ti-O and three Ca-O bonds for a-HAP/a-TiO2. From the stress-strain curve, the Young's modulus of a-Si-HAP/a-TiO2 is calculated to be about 25% higher than that of the a-HAP/a-TiO2, and the yielding stress is about 2 times greater than that of the undoped model. Our calculations therefore demonstrate that the presence of Si in the a-HAP structure strongly alters not only the bioactivity and resorption rates, but also the mechanical properties of the a-HAP/a-TiO2 interface. The results presented here provide an important theoretical insight into the nature of the chemical bonding at the a-HAP/a-TiO2 interface, and are particularly significant for the practical medical applications of HAP-based biomaterials.

Abstract: The compositional influence of Cr and Y on the microstructure of Mg―Cr―O, and Mg―Y―O films synthesized by reactive magnetron sputtering has been investigated by transmission electron microscopy, X-ray diffraction and molecular dynamics simulations. A decrease in crystallinity is observed in these films as the M (Cr or Y) content is increased. It is found that M forms a solid solution with MgO for metal ratios up to ~ 70% and ~ 50% for Cr and Y respectively. Above ~ 70% Cr metal ratio the Mg―Cr―O films are found to be completely amorphous. The Mg―Y―O films are composed of Mg(Y)O and Y2O3 nano crystallites, up to ~ 50% Y metal ratio. Above this ratio, only Y2O3 nano crystallites are found. The preferential < 111> MgO grain alignment is strongly affected by the increase in M content. For M metal ratios up to ~ 50%, there is a selective promotion of the < 100> MgO grain alignments and a decline in the < 111> grain alignments.

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: Microdischarges are an important type of plasma discharge that possess several unique characteristics, such as the presence of a stable glow discharge, high plasma density and intense excimer radiation, leading to several potential applications. The intense and controllable gas heating within the extremely small dimensions of microdischarges has been exploited in microthruster technologies by incorporating a micro-nozzle to generate the thrust. This kind of microthruster has a significantly improved specific impulse performance compared to conventional cold gas thrusters, and can meet the requirements arising from the emerging development and application of micro-spacecraft. In this paper, we performed a self-consistent 2D particle-in-cell simulation, with a Monte Carlo collision model, of a microdischarge operating in a prototype micro-plasma thruster with a hollow cylinder geometry and a divergent micro-nozzle. The model takes into account the thermionic electron emission including the Schottky effect, the secondary electron emission due to cathode bombardment by the plasma ions, several different collision processes, and a non-uniform argon background gas density in the cathode-anode gap. Results in the high-pressure (several hundreds of Torr), high-current (mA) operating regime showing the behavior of the plasma density, potential distribution, and energy flux towards the hollow cathode and anode are presented and discussed. In addition, the results of simulations showing the effect of different argon gas pressures, cathode material work function and discharge voltage on the operation of the microdischarge thruster are presented. Our calculated properties are compared with experimental data under similar conditions and qualitative and quantitative agreements are reached.

Abstract: The effects of hole injection in amorphous indium-gallium-zinc-oxide (a-IGZO) are analyzed by means of first-principles calculations. The injection of holes in the valence band tail states leads to their capture as a polaron, with high self-trapping energies (from 0.44 to 1.15 eV). Once formed, they mediate the formation of peroxides and remain localized close to the hole injection source due to the presence of a large diffusion energy barrier (of at least 0.6 eV). Their diffusion mechanism can be mediated by the presence of hydrogen. The capture of these holes is correlated with the low off-current observed for a-IGZO transistors, as well as with the difficulty to obtain a p-type conductivity. The results further support the formation of peroxides as being the root cause of Negative Bias Illumination Stress (NBIS). The strong self-trapping substantially reduces the injection of holes from the contact and limits the creation of peroxides from a direct hole injection. In the presence of light, the concentration of holes substantially rises and mediates the creation of peroxides, responsible for NBIS. Published by AIP Publishing.

Abstract: Amorphous platinum clusters supported on porous carbon have been envisaged for high-performance fuel cell electrodes. For this application, it is crucial to control the morphology of the Pt layer and the Ptsubstrate interaction to maximize activity and stability. We thus investigate the morphology evolution during Pt cluster growth on a porous carbon substrate employing atomic scale molecular dynamics simulations. The simulations are based on the Pt-C interaction potential using parameters derived from density functional theory and are found to yield a Pt cluster morphology similar to that observed in low loaded fuel cell electrodes prepared by plasma sputtering. Moreover, the simulations show amorphous Pt cluster growth in agreement with X-ray diffraction and transmission electron microscopy experiments on high performance low Pt content (10 μgPt cm−2) loaded fuel cell electrodes and provide a fundamental insight in the cluster growth mechanism.

Abstract: A study of ns-laser ablation is presented, which focuses on the transient behavior of the physical processes that act in and above a copper sample. A dimensionless multiphase collisional radiative model describes the interplay between the ablation, collisional, and radiative mechanisms. Calculations are done for a 6 ns-Nd:YAG laser pulse operating at 532 nm and fluences up to 15 J/cm2. Temporal intensity profiles as well as transmissivities are in good agreement with experimental results. It is found that volumetric ablation mechanisms and photo-processes both play an essential role in the onset of ns-laser induced breakdown.

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: Porous materials are commonly used in microelectronics, as they can meet the demand for continuously shrinking electronic feature dimensions. However, they are facing severe challenges in plasma etching, due to plasma induced damage. In this paper, we present both the plasma characteristics and surface processing during the etching of porous materials. We explain how the damage occurs in the porous material during plasma etching for a wide range of chuck temperatures and the responsible mechanism for plasma damage-free etching at cryogenic temperature, by a combination of experiments and numerical modeling.

Abstract: The unique capabilities of microsecond dwell time (DT) single-particle inductively coupled plasma mass spectrometry (spICPMS) were utilized to characterize the cloud of ions generated from the introduction of suspensions of gold nanoparticles (AuNPs) into the plasma. A set of narrowly distributed particles with diameters ranging from 15.4 to 100.1 nm was synthesized and characterized according to established protocols. Statistically significant numbers of the short transient spICPMS events were evaluated by using 50 μs DT for their summed intensity, maximum intensity, and duration, of which all three were found to depend on the particle diameter. The summed intensity increases from 10 to 1661 counts and the maximum intensity from 6 to 309 counts for AuNPs with diameters from 15.4 to 83.2 nm. The event duration rises from 322 to 1007 μs upon increasing AuNP diameter. These numbers represent a comprehensive set of key data points of the ion clouds generated in ICPMS from AuNPs. The extension of event duration is of high interest to appoint the maximum possible particle number concentration at which separation of consecutive events in spICPMS can still be achieved. Moreover, the combined evaluation of all above-mentioned ion cloud characteristics can explain the regularly observed prolonged single-particle events. The transport and ionization behavior of AuNPs in the ICP was also computationally modeled to gain insight into the size-dependent signal generation. The simulated data reveals that the plasma temperature, and therefore the point of ionization of the particles, is the same for all diameters. However, the maximum number density of Au+, as well as the extent of the ion cloud, depends on the particle diameter, in agreement with the experimental data, and it provides an adequate explanation for the observed ion cloud characteristics.

Abstract: Water permeation between stacked layers of hBN sheets forming 2D nanochannels is investigated using large-scale ab initio-quality molecular dynamics simulations. A high-dimensional neural network potential trained on density-functional theory calculations is employed. We simulate water in van der Waals nanocapillaries and study the impact of nanometric confinement on the structure and dynamics of water using both equilibrium and nonequilibrium methods. At an interlayer distance of 10.2 A confinement induces a first-order phase transition resulting in a well-defined AA-stacked bilayer of hexagonal ice. In contrast, for h < 9 A, the 2D water monolayer consists of a mixture of different locally ordered patterns of squares, pentagons, and hexagons. We found a significant change in the transport properties of confined water, particularly for monolayer water where the water-solid friction coefficient decreases to half and the diffusion coefficient increases by a factor of 4 as compared to bulk water. Accordingly, the slip-velocity is found to increase under confinement and we found that the overall permeation is dominated by monolayer water adjacent to the hBN membranes at extreme confinements. We conclude that monolayer water in addition to bilayer ice has a major contribution to water transport through 2D nanochannels.