Abstract: The Covid era has made us aware of the need for resilient, self-sufficient, and local production. We are likely willing to pay an extra price for that quality. Ammonia (NH3) synthesis accounts for 2 % of global energy production and is an important point of attention for the development of green energy technologies. Therefore, we propose a thermally integrated process for H2 production and NH3 synthesis using plasma technology, and we evaluate its techno-economic performance and CO2 footprint by life cycle assessment (LCA). The key is to integrate energy-wise a high-temperature plasma (HTP) process, with a (low-temperature) non-thermal plasma (NTP) process and to envision their joint economic potential. This particularly means raising the temperature of the NTP process, which is typically below 100 ◦ C, taking advantage of the heat released from the HTP process. For that purpose, we proposed the integrated process and conducted chemical kinetics simulations in the NTP section to determine the thermodynamically feasible operating window of this novel combined plasma process. The results suggest that an NH3 yield of 2.2 mol% can be attained at 302 ◦ C at an energy yield of 1.1 g NH3/kWh. Cost calculations show that the economic performance is far from commercial, mainly because of the too low energy yield of the NTP process. However, when we base our costs on the best literature value and plausible future scenarios for the NTP energy yield, we reach a cost prediction below 452 $/tonne NH3, which is competitive with conventional small-scale Haber-Bosch NH3 synthesis for distributed production. In addition, we demonstrate that biogas can be used as feed, thus allowing the proposed integrated reactor concept to be part of a biogas-to-ammonia circular concept. Moreover, by LCA we demonstrate the environmental benefits of the pro­posed plant, which could cut by half the carbon emissions when supplied by photovoltaic electricity, and even invert the carbon balance when supplied by wind power due to the avoided emissions of the carbon black credits.

Abstract: The treatment of a well plate by an atmospheric pressure plasma jet is common for<italic>in vitro</italic>plasma medicine research. Here, reactive species are largely produced through the mixing of the jet effluent with the surrounding atmosphere. This mixing can be influenced not only by the ambient conditions, but also by the geometry of the treated well. To limit this influence and control the atmosphere, a shielding gas is sometimes applied. However, the interplay between the gas shield and the well geometry has not been investigated. In this work, we developed a 2D-axisymmetric computational fluid dynamics model of the kINPen plasma jet, to study the mixing of the jet effluent with the surrounding atmosphere, with and without gas shield. Our computational and experimental results show that the choice of well type can have a significant influence on the effluent conditions, as well as on the effectiveness of the gas shield. Furthermore, the geometry of the shielding gas device can substantially influence the mixing as well. Our results provide a deeper understanding of how the choice of setup geometry can influence the plasma treatment, even when all other operating parameters are unchanged.

Abstract: Catalysis science and technology increased efforts recently to progress beyond conventional “thermal” catalysis and face the challenges of net-zero emissions and electrification of production. Nevertheless, a better gaps and opportunities analysis is necessary. This review analyses four emerging areas of unconventional or less- conventional catalysis which share the common aspect of using directly renewable energy sources: (i) plasma catalysis, (ii) catalysis for flow chemistry and process intensification, (iii) application of electromagnetic (EM) fields to modulate catalytic activity and (iv) nanoscale generation at the catalyst interface of a strong local EM by plasmonic effect. Plasma catalysis has demonstrated synergistic effects, where the outcome is higher than the sum of both processes alone. Still, the underlying mechanisms are complex, and synergy is not always obtained. There is a crucial need for a better understanding to (i) design catalysts tailored to the plasma environment, (ii) design plasma reactors with optimal transport of plasma species to the catalyst surface, and (iii) tune the plasma conditions so they work in optimal synergy with the catalyst. Microfluidic reactors (flow chemistry) is another emerging sector leading to the intensification of catalytic syntheses, particularly in organic chemistry. New unconventional catalysts must be designed to exploit in full the novel possibilities. With a focus on (a) continuous-flow photocatalysis, (b) electrochemical flow catalysis, (c) microwave flow catalysis and (d) ultra­ sound flow activation, a series of examples are discussed, with also indications on scale-up and process indus­ trialisation. The third area discussed regards the effect on catalytic performances of applying oriented EM fields spanning several orders of magnitude. Under well-defined conditions, gas breakdown and, in some cases, plasma formation generates activated gas phase species. The EM field-driven chemical conversion processes depend further on structured electric/magnetic catalysts, which shape the EM field in strength and direction. Different effects influencing chemical conversion have been reported, including reduced activation energy, surface charging, hot spot generation, and selective local heating. The last topic discussed is complementary to the third, focusing on the possibility of tuning the photo- and electro-catalytic properties by creating a strong localised electrical field with a plasmonic effect. The novel possibilities of hot carriers generated by the plasmonic effect are also discussed. This review thus aims to stimulate the reader to make new, creative catalysis to address the challenges of reaching a carbon-neutral world.

Abstract: Atomistic simulation methods for the quantification of free energies are in wide use. These methods operate by sampling the probability density of a system along a small set of suitable collective variables (CVs), which is, in turn, expressed in the form of a free energy surface (FES). This definition of the FES can capture the relative stability of metastable states but not that of the transition state because the barrier height is not invariant to the choice of CVs. Free energy barriers therefore cannot be consistently computed from the FES. Here, we present a simple approach to calculate the gauge correction necessary to eliminate this inconsistency. Using our procedure, the standard FES as well as its gauge-corrected counterpart can be obtained by reweighing the same simulated trajectory at little additional cost. We apply the method to a number of systems—a particle solvated in a Lennard-Jones fluid, a Diels–Alder reaction, and crystallization of liquid sodium—to demonstrate its ability to produce consistent free energy barriers that correctly capture the kinetics of chemical or physical transformations, and discuss the additional demands it puts on the chosen CVs. Because the FES can be converged at relatively short (sub-ns) time scales, a free energy-based description of reaction kinetics is a particularly attractive option to study chemical processes at more expensive quantum mechanical levels of theory.

Abstract: In this letter, CMOS-compatible Ni/HfO2/TiN resistive random access memory stacks demonstrated attractive unipolar switching properties, showing >10(3) endurance and long retention at 150 degrees C. The Ni bottom electrode (BE) improved the switching yield over the NiSiPt BE. To better understand the unipolar forming mechanism, ab initio simulation and time of flight-secondary ion mass spectroscopy were utilized. Compared to the NiSiPt BE, Ni BE gives larger Ni diffusion in the HfO2 and lower formation enthalpy of Ni2+ species during electrical forming. Both the electrical and physical results supported a Ni-injection mechanism for the filament formation. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.3695078]

Abstract: Superimposition of dual frequencies (DFs) is one of the methods used for controlling plasma distribution in an inductively coupled plasma (ICP) source. The effects of a superimposed DF on the argon plasma characteristics have been investigated using a two-dimensional self-consistent fluid model. When both currents are fixed at 6A, the plasma density drops with decrease in one of the source frequencies due to less efficient heating and the plasma uniformity improves significantly. Moreover, for ICP operated with superimposed DFs (i.e., 4.52MHz/13.56MHz and 2.26MHz/13.56MHz), the current source exhibits the same period as the low frequency (LF) component, and the plasma density is higher than that obtained at a single frequency (i.e., 4.52 and 2.26MHz) with the same total current of 12A. However, at superimposed current frequencies of 6.78MHz/13.56MHz, the plasma density is lower than that obtained at a single frequency of 6.78MHz due to the weaker negative azimuthal electric field between two positive maxima during one period of 6.78MHz. When the superimposed DF ICP operates at 2.26 and 13.56MHz, the rapid oscillations of the induced electric field become weaker during one period of 2.26MHz as the current ratio of 2.26MHz/13.56MHz rises from 24A/7 A to 30A/1 A, and the plasma density drops with the current ratio due to weakened electron heating. The uniformity of plasma increases due to sufficient diffusion under the low-density condition.

Abstract: Molecular dynamics simulations using improved TersoffBrenner potential parameters were performed to investigate Cl+ etching of a {2×1} reconstructed Si(100) surface. Steady-state Si etching accompanying the Cl coverage of the surface is observed. Furthermore, a steady-state chlorinated reaction layer is formed. The thickness of this reaction layer is found to increase with increasing energy. The stoichiometry of SiClx species in the reaction layer is found to be SiCl:SiCl2:SiCl3 = 1.0:0.14:0.008 at 50 eV. These results are in excellent agreement with available experimental data. While elemental Si products are created by physical sputtering, most SiClx (0<x<4) etch products are produced by chemical-enhanced physical sputtering.

Abstract: In metallized film capacitors, there exists an air gap of about 0.2 μm between the films, with a pressure ranging generally from 130 atm. Because of the created potential difference between the two films, a microdischarge is formed in this gap. In this paper, we use an implicit particle-in-cell Monte Carlo collision simulation method to study the discharge properties in this direct-current microdischarge with 0.2 μm gap in a range of different voltages and pressures. The discharge process is significantly different from a conventional high pressure discharge. Indeed, the high electric field due to the small gap sustains the discharge by field emission. At low applied voltage (~15 V), only the electrons are generated by field emission, while both electrons and ions are generated as a stable glow discharge at medium applied voltage (~50 V). At still higher applied voltage (~100 V), the number of electrons and ions rapidly multiplies, the electric field reverses, and the discharge changes from a glow to an arc regime.

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: The limitation in time and length scale is a major issue of molecular dynamics (MD) simulation. Although several methods have been developed to extend the MD time scale, their performance usually deteriorates with increasing system size. Therefore, an acceleration method which is applicable to large systems is required to bridge the gap between the MD simulations and target phenomena. In this study, an accelerated MD method for large system is developed based on the collective variable-driven hyperdynamics (CVHD) method [K.M. Bal and E.C. Neyts, 2015]. The key idea is to run CVHD in parallel with rate control and accelerate multiple possible events simultaneously. Using this novel method, carbon diffusion in bcc-iron bicrystal with grain boundary is examined as an application for practical materials. Carbon atoms reaching at the grain boundary are trapped whereas carbon atoms in the bulk region diffuse randomly, and both dynamic regimes can be simultaneously accelerated with the parallel CVHD technique.

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: Piezoelectricity is considered to be one of the key functionalities in biomaterials to boost bone tissue regeneration, however, integrating biocompatibility, biodegradability and 3D structure with pronounced piezoresponse remains a material challenge. Herein, novel hybrid biocompatible 3D scaffolds based on biodegradable poly(3-hydroxybutyrate) (PHB) and reduced graphene oxide (rGO) flakes have been developed. Nanoscale insights revealed a more homogenous distribution and superior surface potential values of PHB fibers (33 +/- 29 mV) with increasing rGO content up to 1.0 wt% (314 +/- 31 mV). The maximum effective piezoresponse was detected at 0.7 wt% rGO content, demonstrating 2.5 and 1.7 times higher out-of-plane and in-plane values, respectively, than that for pure PHB fibers. The rGO addition led to enhanced zigzag chain formation between paired lamellae in PHB fibers. In contrast, a further increase in rGO content reduced the alpha-crystal size and prevented zigzag chain conformation. A corresponding model explaining structural and molecular changes caused by rGO addition in electrospun PHB fibers is proposed. In addition, finite element analysis revealed a negligible vertical piezoresponse compared to lateral piezoresponse in uniaxially oriented PHB fibers based on alpha-phase (P2(1)2(1)2(1) space group). Thus, the present study demonstrates promising results for the development of biodegradable hybrid 3D scaffolds with an enhanced piezoresponse for various tissue engineering applications.

Abstract: Nowadays, lasers are commonly applied in spectrochemical analysis methods, for sampling, plasma formation or a combination of both. Despite the numerous investigations that have been performed on these applications, the underlying processes are still insufficiently understood. In order to fasten progress in the field and in honor of the lifework of professor Rick Russo, we here provide a brief overview of the fundamental mechanisms in lasermatter interaction as proposed in literature, and throw the spotlight on some aspects that have not received much attention yet. For an organized discussion, we choose laser ablation, laser desorption and the associated gaseous plasma formation as the central processes in this perspective article, based on a classification of the laserbased spectrochemical analysis techniques and the corresponding laser-matter interaction regimes. First, we put the looking glass over the excitation and thermalization mechanisms in the laser-irradiated condensed phase, for which we propose the so-called multi-plasma model. This novel model can be understood as an extension of the well-known two-temperature model, featuring multiple thermodynamic dimensions, each of which corresponds to a quasi-particle type. Next, the focus is placed on the mass transfer and ionization mechanisms, after which we shortly highlight the possible role of anisotropic and magnetic effects in the laser-excited material. We hope this perspective article motivates more fundamental research on laser-matter interaction, as a continuation of the lifework of Rick Russo.

Abstract: Plasma catalysis is gaining increasing interest for various environmental applications. Catalytic
material can be inserted in different shapes in the plasma, e.g., as pellets, (coated) beads, but also
as honeycomb monolith and 3DFD structures, also called ‘structured catalysts’, which have high
mass and heat transfer properties. In this work, we examine the streamer discharge propagation
and the interaction between plasma and catalysts, inside the channels of such structured catalysts,
by means of a two-dimensional particle-in-cell/Monte Carlo collision model. Our results reveal
that plasma streamers behave differently in various structured catalysts. In case of a honeycomb
structure, the streamers are limited to only one channel, with low or high plasma density when
the channels are parallel or perpendicular to the electrodes, respectively. In contrast, in case of a
3DFD structure, the streamers can distribute to different channels, causing discharge
enhancement due to surface charging on the dielectric walls of the structured catalyst, and
especially giving rise to a broader plasma distribution. The latter should be beneficial for plasma
catalysis applications, as it allows a larger catalyst surface area to be exposed to the plasma.

Abstract: We perform first-principles calculations to investigate the preferred positions of B and P dopants, both neutral and in their preferred charge state, in Si and Si/SiO2 core-shell nanowires (NWs). In order to understand the observed trends in the formation energy, we isolate the different effects that determine these formation energies. By making the distinction between the unrelaxed and the relaxed formation energy, we separate the impact of the relaxation from that of the chemical environment. The unrelaxed formation energies are determined by three effects: (i) the effect of strain caused by size mismatch between the dopant and the host atoms, (ii) the local position of the band edges, and (iii) a screening effect. In the case of the SiNW (Si/SiO2 NW), these effects result in an increase of the formation energy away from the center (interface). The effect of relaxation depends on the relative size mismatch between the dopant and host atoms. A large size mismatch causes substantial relaxation that reduces the formation energy considerably, with the relaxation being more pronounced towards the edge of the wires. These effects explain the surface segregation of the B dopants in a SiNW, since the atomic relaxation induces a continuous drop of the formation energy towards the edge. However, for the P dopants, the formation energy starts to rise when moving from the center but drops to a minimum just next to the surface, indicating a different type of behavior. It also explains that the preferential location for B dopants in Si/SiO2 core-shell NWs is inside the oxide shell just next to the interface, whereas the P dopants prefer the positions next to the interface inside the Si core, which is in agreement with recent experiments. These preferred locations have an important impact on the electronic properties of these core-shell NWs. Our simulations indicate the possibility of hole gas formation when B segregates into the oxide shell.