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: A gliding arc (GA) plasma has great potential for CO2 conversion into value-added chemicals, because of its high energy efficiency. To improve the application, a 2D/3D fluid model is needed to investigate the CO2 conversion mechanisms in the actual discharge geometry. Therefore, the complex CO2 chemical kinetics description must be reduced due to the huge computational cost associated with 2D/3D models. This paper presents a chemistry reduction method for CO2 plasmas, based on the so-called directed relation graph method. Depending on the defined threshold values, some marginal species are identified. By means of a sensitivity analysis, we can further reduce the chemistry set by removing one by one the marginal species. Based on the socalled flux-sensitivity coupling, we obtain a reduced CO2 kinetics model, consisting of 36 or 15 species (depending on whether the 21 asymmetric mode vibrational states of CO2 are explicitly included or lumped into one group), which is applied to a GA discharge. The results are compared with those predicted with the full chemistry set, and very good agreement is reached. Moreover, the range of validity of the reduced CO2 chemistry set is checked, telling us that this reduced set is suitable for low power GA discharges. Finally, the time and spatial evolution of the CO2 plasma characteristics are presented, based on a 2D model with the reduced kinetics.

Abstract: This study focuses on a magnetically stabilized gliding arc (MGA) plasma. Two fully coupled flow-plasma models (in 3D and 2D) are presented. The 3D model is applied to compare the arc dynamics of the MGA with a traditional gas-driven gliding arc. The 2D model is used for a detailed parametric study on the effect of the external magnetic field. The results show that the relative velocity between the plasma and feed gas is generated due to the Lorentz force, which can increase the plasma-treated gas fraction. The magnetic field also helps to decrease the gas temperature by enhancing heat transfer and to increase the electron number density. This work shows the potential of an external magnetic field to control the gliding arc behavior, for enhanced gas conversion at low gas flow rates.

Abstract: We investigate discharge contraction in a microwave plasma at sub-atmospheric pressure, operating in CO2 and CO2/CH4 mixtures. The rise of the electron number density with plasma contraction intensifies the gas heating in the core of the plasma. This, in turn, initiates fast core-periphery transport and defines the rate of thermal chemistry over plasma chemistry. In this context, power concentration describes the overall mechanism including plasma contraction and chemical kinetics. In a complex chemistry such as dry reforming of methane, transport of reactive species is essential to define the performance of the reactor and achieve the desired outputs. Thus, we couple experimental observations and thermodynamic calculations for model validation and understanding of reactor performance. Adding CH4 alters the thermodynamic properties of the mixture, especially the reactive component of the heat conductivity. The increase in reactive heat conductivity increases the pressure at which plasma contraction occurs, because higher rates of gas heating are required to reach the same temperature. In addition, we suggest that the predominance of heat conduction over convection is a key condition to observe the effect of heat conductivity on gas temperature.

Abstract: Since the onset of pattern transfer technologies for chip manufacturing, various strategies have been developed to circumvent or overcome aspect ratio dependent etching (ARDE). These methods have, however, their own limitations in terms of etch non-idealities, throughput or costs. Moreover, they have mainly been optimized for individual in-device features and die-scale patterns, while occasionally ending up with poor patterning of metrology marks, affecting the alignment and overlay in lithography. Obtaining a better understanding of the underlying mechanisms of ARDE and how to mitigate them therefore remains a relevant challenge to date, for both marks and advanced nodes. In this work, we accordingly assessed the neutral transport mechanisms in ARDE by means of experiments and multiscale modeling for SiO₂etching with CHF₃/Ar and CF₄/Ar plasmas. The experiments revealed a local maximum in the etch rate for an aspect ratio around unity, i.e. the simultaneous occurrence of regular and inverse reactive ion etching lag for a given etch condition. We were able to reproduce this ARDE trend in the simulations without taking into account charging effects and the polymer layer thickness, suggesting shadowing and diffuse reflection of neutrals as the primary underlying mechanisms. Subsequently, we explored four methods with the simulations to regulate ARDE, by varying the incident plasma species fluxes, the amount of polymer deposition, the ion energy and angular distribution and the initial hardmask sidewall angle, for which the latter was found to be promising in particular. Although our study focusses on feature dimensions characteristic to metrology marks and back-end-of-the-line integration, the obtained insights have a broader relevance, e.g. to the patterning of advanced nodes. Additionally, this work supports the insight that physisorption may be more important in plasma etching at room temperature than originally thought, in line with other recent studies, a topic on which we recommend further research.

Abstract: Non-equilibrium inductively coupled plasmas (ICPs) operating in hydrogen are of significant interest for applications including large-area materials processing. Increasing control of spatial gas heating, which drives the formation of neutral species density gradients and the rate of gas-temperature-dependent reactions, is critical. In this study, we use 2D fluid-kinetic simulations with the Hybrid Plasma Equipment Model to investigate the spatially resolved production of atomic hydrogen in a low-pressure planar ICP operating in pure hydrogen (10-20 Pa or 0.075-0.15 Torr, 300 W). The reaction set incorporates self-consistent calculation of the spatially resolved gas temperature and 14 vibrationally excited states. We find that the formation of neutral-gas density gradients, which result from spatially non-uniform electrical power deposition at constant pressure, can drive significant variations in the vibrational distribution function and density of atomic hydrogen when gas heating is spatially resolved. This highlights the significance of spatial gas heating on the production of reactive species in relatively high-power-density plasma processing sources.

Abstract: It is well established that lipid aldehydes (LAs) are able to increase the permeability of cell membranes and induce their rupture. However, it is not yet clear how LAs are distributed in phase-separated membranes (PSMs), which are responsible for the transport of selected molecules and intracellular signaling. Thus, we investigate here the distribution of LAs in a PSM by coarse-grained molecular dynamics simulations. Our results reveal that LAs derived from mono-unsaturated lipids tend to accumulate at the interface between the liquid-ordered/liquiddisordered domains, whereas those derived from poly-unsaturated lipids remain in the liquid-disordered domain. These results are important for understanding the effects caused by oxidized lipids in membrane structure, properties and organization.

Abstract: This work proposes to use core-shell structured spheres to evaluate whether it allows to individually optimize bulk and surface effects of a packing material, in order to optimize conversion and energy efficiency. Different core-shell materials have been prepared by spray coating, using dense spheres (as core) and powders (as shell) of SiO2, Al2O3, and BaTiO3. The materials are investigated for their performance in CO2 dissociation and compared against a benchmark consisting of a packed-bed reactor with the pure dense spheres, as well as an empty reactor. The results in terms of CO2 conversion and energy efficiency show various interactions between the core and shell material, depending on their combination. Al2O3 was found as the best core material under the applied conditions here, followed by BaTiO3 and SiO2, in agreement with their behaviour for the pure spheres. Applying a thin shell layer on the cores showed equal performance between the different shell materials. Increasing the layer thickness shifts this behaviour, and strong combination effects were observed depending on the specific material. Therefore, this method of core-shell spheres has the potential to allow tuning of the packing properties more closely to the application by designing an optimal combination of core and shell.

Abstract: Herein, the efficacy of WOx-promoted CeO2-SiO2 and CeO2-ZrO2 mixed oxide catalysts in the solvent-free selective oxidation of benzyl alcohol to benzaldehyde using molecular oxygen as an oxidant is reported. We evaluated the effects of the oxidant and catalyst concentration, reaction duration, and temperature on the reaction with an aim to optimize the reaction conditions. The as-prepared CeO2, CeO2-ZrO2, CeO2-SiO2, WOx/CeO2, WOx/CeO2-ZrO2, and WOx/CeO2-SiO2 catalysts were characterized by X-ray diffraction (XRD), N-2 adsorption-desorption, Raman spectroscopy, temperature-programmed desorption of ammonia (TPD-NH3), X-ray photoelectron spectroscopy (XPS), and transmission electron microscopy (TEM). These characterisation results indicated that the WOx/CeO2-SiO2 catalyst possessed improved physicochemical (i.e., structural, textural, and acidic) properties owing to the strong interactivity between WOx and CeO2-SiO2. A higher number of Ce3+ ions (I-u '''/I-Total) were created with the WOx/CeO2-SiO2 catalyst than those with the other catalysts in this work, indicating the generation of a high number of oxygen vacancies. The WOx/CeO2-SiO2 catalyst exhibited a high conversion of benzyl alcohol (>99%) and a high selectivity (100%) toward benzaldehyde compared to the other promoted catalysts (i.e., WOx/CeO2 and WOx/CeO2-ZrO2), which is attributed to the smaller particle size of the WOx and CeO2 and their high specific surface area, more significant number of acidic sites, and superior number of oxygen vacancies. The WOx/CeO2-SiO2 catalyst could be quickly recovered and utilized at least five times without suffering any appreciable activity loss.

Abstract: The overexpression of voltage dependent anion channels (VDACs), particularly VDAC1, in cancer cells compared to normal cells, plays a crucial role in cancer cell metabolism, apoptosis regulation, and energy homeostasis. In this study, we used molecular dynamics (MD) simulations to investigate the effect of a low level of VDAC1 oxidation (induced e.g., by cold atmospheric plasma (CAP)) on the pyruvate (Pyr) uptake by VDAC1. Inhibiting Pyr uptake through VDAC1 can suppress cancer cell proliferation. Our primary target was to study the translocation of Pyr across the native and oxidized forms of hVDAC1, the human VDAC1. Specifically, we employed MD simulations to analyze the hVDAC1 structure by modifying certain cysteine residues to cysteic acids and methionine residues to methionine sulfoxides, which allowed us to investigate the effect of oxidation. Our results showed that the free energy barrier for Pyr translocation through the native and oxidized channel was approximately 4.3 +/- 0.7 kJ mol-1 and 10.8 +/- 1.8 kJ mol-1, respectively. An increase in barrier results in a decrease in rate of Pyr permeation through the oxidized channel. Thus, our results indicate that low levels of CAP oxidation reduce Pyr translocation, resulting in decreased cancer cell proliferation. Therefore, low levels of oxidation are likely sufficient to treat cancer cells given the inhibition of Pyr uptake.

Abstract: Electrification and carbon capture technologies are essential for achieving net-zero emissions in the chemical sector. A crucial strategy involves converting captured CO2 into CO, a valuable chemical feedstock. This study evaluates the feasibility of two innovative methods: plasma activation and electrolysis, using clean electricity and captured CO2. Specifically, it compares a gliding arc plasma reactor with an embedded novel carbon bed system to a modern zero-gap type low-temperature electrolyser. The plasma method stood out with an energy cost of 19.5 GJ per tonne CO, marking a 43% reduction compared to electrolysis and conventional methods. CO production costs for plasma- and electrolysis-based plants were $671 and $962 per tonne, respectively. However, due to high uncertainty regarding electrolyser costs, the CO production costs in electrolysis-based plants may actually range from $570 to $1392 per tonne. The carbon bed system in the plasma method was a key factor in facilitating additional CO generation from O-2 and enhancing CO2 conversion, contributing to its cost-effectiveness. Challenges for electrolysis included high costs of equipment and low current densities. Addressing these limitations could significantly decrease production costs, but challenges arise from the mutual relationship between intrinsic parameters, such as CO2 conversion, CO2 input flow, or energy cost. In a future scenario with affordable feedstocks and equipment, costs could drop below $500 per tonne for both methods. While this may be more challenging for electrolysis due to complexity and expensive catalysts, plasma-based CO production appears more viable and competitive.

Abstract: Asymmetric dye molecules have unusual optical and electronic properties1, 2, 3. For instance, they show a strong second-order nonlinear optical (NLO) response that has attracted great interest for potential applications in electro-optic modulators for optical telecommunications and in wavelength conversion of lasers2, 3. However, the strong Coulombic interaction between the large dipole moments of these molecules favours a pairwise antiparallel alignment that cancels out the NLO response when incorporated into bulk materials. Here, we show that by including an elongated dipolar dye (p,p′-dimethylaminonitrostilbene, DANS, a prototypical asymmetric dye with a strong NLO response4) inside single-walled carbon nanotubes (SWCNTs)5, 6, an ideal head-to-tail alignment in which all electric dipoles point in the same sense is naturally created. We have applied this concept to synthesize solution-processible DANS-filled SWCNTs that show an extremely large total dipole moment and static hyperpolarizability (β0 = 9,800 × 10−30 e.s.u.), resulting from the coherent alignment of arrays of ∼70 DANS molecules.

Abstract: The behavior of hydrocarbon species at step edges of diamond terraces is investigated by means of combined molecular dynamics−Metropolis Monte Carlo simulations. The results show that the formation of ballas-like diamond films (like UNCD) and well-faceted diamond films (like NCD) can be related to the gas phase concentrations of CxHy in a new manner: Species that have high concentrations above the growing UNCD films suppress the extension of step edges through defect formation. The species that are present above the growing NCD film, however, enhance the extension of diamond terraces, which is believed to result in well-faceted diamond films. Furthermore, it is shown that, during UNCD growth, CxHy species with x ≥ 2 play an important role, in contrast to the currently adopted CVD diamond growth mechanism. Finally, the probabilities for the extension of the diamond (100) terrace are much higher than those for the diamond (111) terrace, which is in full agreement with the experimental observation that diamond (100) facets are more favored than diamond (111) facets during CVD diamond growth.

Abstract: It is known that film structure may change dramatically with the extent of surface diffusion during the film growth process. In the present work, surface diffusion, induced thermally or activated by energetic impacts, is investigated theoretically under conditions appropriate for magnetron sputter-deposition of MgAlO thin films with varying stoichiometry. The distribution of surface diffusion energy barriers available to the system was determined for each stoichiometry, which allowed assessing in a qualitative way how much surface diffusion will take place on the time scale available between deposition events. The activation energy barriers increase with the Al concentration in the film, and therefore, the surface diffusion rates in the time frame of typical deposition rates drop, which can explain the decrease in crystallinity in the film structure and the transition to amorphous structure. The deposition process and the immediate surface diffusion enhanced by the energetic adatoms are simulated by means of a molecular dynamics model. The longer-time thermal surface diffusion and the energy landscape are studied by the temperature accelerated dynamics method, applied in an approximate way. The surface diffusion enhanced by the energetic impacts appears to be very important for the film structure in the low-temperature deposition regime.

Abstract: van der Waals (vdW) forces play an important role in the adsorption of molecules on the surface of solids. However, the choice of the most suitable vdW functional for different systems is an essential problem which must be addressed for different systems. The lack of a systematic study on the performance of the vdW functionals in the bulk and adsorption properties of metal-oxides motivated us to examine different vdW approaches and compute the bulk and molecular adsorption properties of alpha-Cr2O3, alpha-Fe2O3, and alpha-Al2O3. For the bulk properties, we compared our results for the heat of formation, cohesive energy, lattice parameters and bond distances between the different vdW functionals and available experimental data. Next we studied the adsorption of benzene and CH3 molecules on top of different oxide surfaces. We employed different approximations to exchange and correlation within DFT, namely, the Perdew-Burke-Ernzerhof (PBE) GGA, (PBE)+U, and vdW density functionals [ DFT(vdW-DF/DF2/optPBE/optB86b/optB88)+U] as well as DFT-D2/D3(+U) methods of Grimme for the bulk calculations and optB86b-vdW(+U) and DFT-D2(+U) for the adsorption energy calculations. Our results highlight the importance of vdW interactions not only in the adsorption of molecules, but importantly also for the bulk properties. Although the vdW contribution in the adsorption of CH3 (as a chemisorption interaction) is less important compared to the adsorption of benzene (as a physisorption interaction), this contribution is not negligible. Also adsorption of benzene on ferryl/chromyl terminated surfaces shows an important chemisorption contribution in which the vdW interactions become less significant.