Abstract: We present a zero-dimensional plasma kinetics model, including both surface and gas phase kinetics, to determine the role of vibrationally excited states in plasma-catalytic ammonia synthesis. We defined a new method to systematically capture the conditions of dielectric barrier discharges (DBDs), including those found in packed bed DBDs. We included the spatial and temporal nature of such discharges by special consideration of the number of micro-discharges in the model. We introduce a parameter that assigns only a part of the plasma power to the microdischarges, to scale the model conditions from filamentary to uniform plasma. Because of the spatial and temporal behaviour of the micro-discharges, not all micro-discharges occurring in the plasma reactor during a certain gas residence time are affecting the molecules. The fraction of power considered in the model ranges from 0.005 %, for filamentary plasma, to 100 %, for uniform plasma. If vibrational excitation is included in the plasma chemistry, these different conditions, however, yield an ammonia density that is only varying within one order of magnitude. At only 0.05 % of the power put into the uniform plasma component, a model neglecting vibrational excitation clearly does not result in adequate amounts of ammonia. Thus, our new model, which accounts for the concept in which not all the power is deposited by the micro-discharges, but some part may also be distributed in between them, suggests that vibrational kinetic processes are really important in (packed bed) DBDs. Indeed, vibrational excitation takes place in both the uniform plasma between the micro-discharges and in the strong micro-discharges, and is responsible for an increased N2 dissociation rate. This is shown here for plasma-catalytic ammonia synthesis, but might also be valid for other gas conversion applications.

Abstract: Plasma has been studied for several years to convert CO2 into value-added products. If CO2 could be converted in the presence of H2O as a cheap H-source for making syngas and oxygenates, it would mimic natural photosynthesis. However, CO2/H2O plasmas have not yet been extensively studied, not by experiments, and certainly not computationally. Therefore, we present here a kinetic modelling study to obtain a greater understanding of the vibrational kinetics of a CO2/H2O microwave plasma. For this purpose, we first created an electron impact cross section set for H2O, using a swarm-derived method. We added the new cross section set and CO2/H2O-related chemistry to a pure CO2 model. While it was expected that H2O addition mainly causes quenching of the CO2 asymmetric mode vibrational levels due to the additional CO2/H2O vibrational-translational relaxation, our model shows that the modifications in the vibrational kinetics are mainly induced by the strong electron dissociative attachment to H2O molecules, causing a reduction in electron density, and the corresponding changes in the input of energy into the CO2 vibrational levels by electron impact processes.

Abstract: Plasma catalysis is gaining increasing interest for various applications, but the underlying mechanisms are still far from understood. Hence, more fundamental research is needed to understand these mechanisms. This can be obtained by both modelling and experiments. This foundations paper describes the fundamental insights in plasma catalysis, as well as efforts to gain more insights by modelling and experiments. Furthermore, it discusses the state-of-the-art of the major plasma catalysis applications, as well as successes and challenges of technology transfer of these applications.

Abstract: Vibrational excitation represents an efficient channel to drive the dissociation of CO₂in a non-thermal plasma. Its viability is investigated in low-pressure pulsed discharges, with the intention of selectively exciting the asymmetric stretching mode, leading to stepwise excitation up to the dissociation limit of the molecule. Gas heating is crucial for the attainability of this process, since the efficiency of vibration–translation (V–T) relaxation strongly depends on temperature, creating a feedback mechanism that can ultimately thermalize the discharge. Indeed, recent experiments demonstrated that the timeframe of V–T non-equilibrium is limited to a few milliseconds at ca. 6 mbar, and shrinks to the<italic>μ</italic>s-scale at 100 mbar. With the aim of backtracking the origin of gas heating in pure CO₂plasma, we perform a kinetic study to describe the energy transfers under typical non-thermal plasma conditions. The validation of our kinetic scheme with pulsed glow discharge experiments enables to depict the gas heating dynamics. In particular, we pinpoint the role of vibration–vibration–translation relaxation in redistributing the energy from asymmetric to symmetric levels of CO₂, and the importance of collisional quenching of CO₂electronic states in triggering the heating feedback mechanism in the sub-millisecond scale. This latter finding represents a novelty for the modelling of low-pressure pulsed discharges and we suggest that more attention should be paid to it in future studies. Additionally, O atoms convert vibrational energy into heat, speeding up the feedback loop. The efficiency of these heating pathways, even at relatively low gas temperature and pressure, underpins the lifetime of V–T non-equilibrium and suggests a redefinition of the optimal conditions to exploit the ‘ladder-climbing’ mechanism in CO₂discharges.

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: We have developed a comprehensive kinetic model to study the O atom kinetics in an O₂plasma and its afterglow. By adopting a pseudo-1D plug-flow formalism within the kinetic model, our aim is to assess how far the O atoms travel in the plasma afterglow, evaluating its potential as a source of O atoms for post-plasma gas conversion applications. Since we could not find experimental data for pure O₂plasma at atmospheric pressure, we first validated our model at low pressure (1–10 Torr) where very good experimental data are available. Good agreement between our model and experiments was achieved for the reduced electric field, gas temperature and the densities of the dominant neutral species, i.e. O₂(a), O₂(b) and O. Subsequently, we confirmed that the chemistry set is consistent with thermodynamic equilibrium calculations at atmospheric pressure. Finally, we investigated the O atom densities in the O₂plasma and its afterglow, for which we considered a microwave O₂plasma torch, operating at a pressure between 0.1 and 1 atm, for a flow rate of 20 slm and an specific energy input of 1656 kJ mol⁻¹. Our results show that for both pressure conditions, a high dissociation degree of ca. 92% is reached within the discharge. However, the O atoms travel much further in the plasma afterglow for<italic>p</italic>= 0.1 atm (9.7 cm) than for<italic>p</italic>= 1 atm (1.4 cm), attributed to the longer lifetime (3.8 ms at 0.1 atm vs 1.8 ms at 1 atm) resulting from slower three-body recombination kinetics, as well as a higher volumetric flow rate.

Abstract: The lowest-energy configurations of Cn(n≤55) clusters are obtained using the energy-minimization technique with the conjugate gradient method where a modified Brenner potential is invoked to describe the carbon and hydrocarbon interaction. We found that the ground-state configuration consists of a single ring for small number of C atoms and multiring structures are found with increasing n, which can be in planar, bowl-like or caplike form. Contrary to previous predictions, the binding energy Eb does not show even-odd oscillations and only small jumps are found in the Eb(n) curve as a consequence of specific types of edges or equivalently the number of secondary atoms. We found that hydrogenation of the edge atoms may change the ground-state configuration of the nanocluster. In both cases we determined the magic clusters. Special attention is paid to trigonal and hexagonal shaped carbon clusters and to clusters having a graphenelike configuration. Trigonal clusters are never the ground state while hexagonal-shaped clusters are only the ground state when they have zigzag edges.

Abstract: We propose that low-strain Si/O superlattices can be constructed by connecting reconstructed Si{001} surfaces by Si-O-Si bridges. Ab initio calculations show that our models are energetically more favorable than all the models proposed so far. The part of our Si/O superlattice model is experimentally accessible just by oxidizing a Si( 001) substrate. To complete our Si/O superlattice model, we propose a three-step method. We also explore the potential of our Si/O superlattice models for new materials used in future Si electronics. We find that the location of the channel where the carriers travel can be controlled between the interfaces and the Si layers by the insertion of O atoms into the Si-Si dimers. By revealing the origins of the interface electron and hole states, we find that similar interface states should be easily achieved for Si slabs and Si substrates. Interestingly, the interface electrons and holes have small effective masses in the direction parallel to the channel and large effective masses in the direction normal to the channel, which makes the Si/O superlattices attractive to be used for channel materials. We also find that the valley splitting of Si is enhanced by the formation of the Si/O/Si interfaces, which is ideal for developing Si-based qubits. Our findings open new perspectives to design and control the electronic properties of Si.

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: To investigate the impact of both the CoFe ferromagnetic layer thickness and the capping paramagnetic layer on the tunnel magnetoresistance (TMR), we performed first-principles simulations on epitaxial magnetic tunnel junctions contacted with either CoFe or Ta paramagnetic capping layers. We observed a strong oscillation of the TMR amplitude with respect to the thickness of the ferromagnetic layer. The TMR is found to be amplified whenever the MgO spin tunnel barrier is thickened. Quantization of the electronic structure of the ferromagnetic layers is found to be at the origin of this oscillatory behavior. Metals such as Ta contacting the magnetic layer are found to enhance the amplitude of the oscillations due to the occurrence of an interface dipole. The latter drives the band alignment and tunes the nature of the spin channels that are active during the tunneling process. Subsequently, the regular transmission spin channels are modulated in the magnetic tunnel junction stack and other complex ones are being activated.

Abstract: Amorphous semiconductors are usually characterized by a low charge carrier mobility, essentially related to their lack of long-range order. The development of such material with higher charge carrier mobility is hence challenging. Part of the issue comes from the difficulty encountered by first-principles simulations to evaluate concepts such as the electron effective mass for disordered systems since the absence of periodicity induced by the disorder precludes the use of common concepts derived from condensed matter physics. In this paper, we propose a methodology based on first-principles simulations that partially solves this problem, by quantifying the degree of delocalization of a wave function and of the connectivity between the atomic sites within this electronic state. We validate the robustness of the proposed formalism on crystalline and molecular systems and extend the insights gained to disordered/amorphous InGaZnO4 and Si. We also explore the properties of p-type oxide semiconductor candidates recently reported to have a low effective mass in their crystalline phases [G. Hautier et al., Nat. Commun. 4, 2292 (2013)]. Although in their amorphous phase none of the candidates present a valence band with delocalization properties matching those found in the conduction band of amorphous InGaZnO4, three of the seven analyzed materials show some potential. The most promising candidate, K2Sn2O3, is expected to possess in its amorphous phase a slightly higher hole mobility than the electron mobility in amorphous silicon.

Abstract: Biological membranes are under constant attack of free radicals, which may lead to lipid nitro-oxidation, pro­ ducing a complex mixture of nitro-oxidized lipids that are responsible for structural and dynamic changes on the membrane. Despite the latter, nitro-oxidized lipids are also associated with several inflammatory and neuro­ degenerative diseases, the underlying mechanisms of which remain elusive. We perform atomistic molecular dynamics simulations using several isomers of nitro-oxidized lipids to study their effect on the structure and permeability of the membrane, as well as the interaction between the mixture of these products in the phos­pholipid membrane environment. Our results show that the stereo- and positional isomers have a stronger effect on the properties of the membrane composed of oxidized lipids compared to that containing nitrated lipids. Nevertheless, nitrated lipids lead to three-fold increase in water permeability compared to oxidized lipids. In addition, we show that in a membrane consisting of combined nitro-oxidized lipid products, the presence of oxidized lipids protects the membrane from transient pores. Is well stablished that plasma application and photodynamic therapy produces a number of oxidative species used to kill cancer cells, through membrane damage induced by nitro-oxidative stress. This study is important to elucidate the mechanisms and the molecular level properties involving the reactive species produced during that cancer therapies.

Abstract: This Special Issue collects three reviews, eight articles, and two communications related to the design of catalysts for environmental applications, such as the transformation of several pollutants into harmless or valuable products [...]

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.