Abstract: Plasma medicine, or the biomedical application of cold atmospheric plasma (CAP), is an expanding field within plasma research. CAP has demonstrated remarkable versatility in diverse biological applications, including cancer treatment, wound healing, microorganism inactivation, and skin disease therapy. However, the precise mechanisms underlying the effects of CAP remain incompletely understood. The therapeutic effects of CAP are largely attributed to the generation of reactive oxygen and nitrogen species (RONS), which play a crucial role in the biological responses induced by CAP. Specifically, RONS produced during CAP treatment have the ability to chemically modify cell membranes and membrane proteins, causing nitro-oxidative stress, thereby leading to changes in membrane permeability and disruption of cellular processes. To gain atomic-level insights into these interactions, non-reactive molecular dynamics (MD) simulations have emerged as a valuable tool. These simulations facilitate the examination of larger-scale system dynamics, including protein-protein and protein-membrane interactions. In this comprehensive review, we focus on the applications of non-reactive MD simulations in studying the effects of CAP on cellular components and interactions at the atomic level, providing a detailed overview of the potential of CAP in medicine. We also review the results of other MD studies that are not related to plasma medicine but explore the effects of nitro-oxidative stress on cellular components and are therefore important for a broader understanding of the underlying processes.

Abstract: A two-dimensional hybrid Monte Carlofluid model, incorporating a full-wave solution of Maxwell's equations, is employed to describe the behavior of high frequency (HF) and very high frequency capacitively coupled plasmas (CCPs), operating both at single frequency (SF) and dual frequency (DF) in a CF4/O2 gas mixture. First, the authors investigate the plasma composition, and the simulations reveal that besides CF4 and O2, also COF2, CF3, and CO2 are important neutral species, and CF+3 and F− are the most important positive and negative ions. Second, by comparing the results of the model with and without taking into account the electromagnetic effects for a SF CCP, it is clear that the electromagnetic effects are important, both at 27 and 60 MHz, because they affect the absolute values of the calculation results and also (to some extent) the spatial profiles, which accordingly affects the uniformity in plasma processing. In order to improve the plasma radial uniformity, which is important for the etch process, a low frequency (LF) source is added to the discharge. Therefore, in the major part of the paper, the plasma uniformity is investigated for both SF and DF CCPs, operating at a HF of 27 and 60 MHz and a LF of 2 MHz. For this purpose, the authors measure the etch rates as a function of position on the wafer in a wide range of LF powers, and the authors compare them with the calculated fluxes toward the wafer of the plasma species playing a role in the etch process, to explain the trends in the measured etch rate profiles. It is found that at a HF of 60 MHz, the uniformity of the etch rate is effectively improved by adding a LF power of 2 MHz and 300 W, while its absolute value increases by about 50%, thus a high etch rate with a uniform distribution is observed under this condition.

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: Self-organized striated structure has been observed experimentally and numerically in CF4 plasmas in radio-frequency capacitively coupled plasmas recently (Liu et al 2016 Phys. Rev. Lett. 116 255002). In this work, the striated structure is investigated in a capacitively coupled oxygen discharge with the introduction of the effect from the secondary electron emission, based on a particle-in-cell/Monte Carlo collision model. As we know, the transport of positive and negative ions plays a key role in the formation of striations in electronegative gases, for which, the electronegativity needs to be large enough. As the secondary electron emission increases, electrons in the sheaths gradually contribute more ionization to the discharge. Meanwhile, the increase of the electron density, especially in the plasma bulk, leads to an increased electrical conductivity and a reduced bulk electric field, which would shield the ions' mobility. These changes result in enlarged striation gaps. And then, with more emitted electrons, obvious disruption of the striations is observed accompanied with a transition of electron heating mode. Due to the weakened field, the impact ionization in the plasma bulk is attenuated, compared with the enhanced ionization caused by secondary electrons. This would lead to the electron heating mode transition from striated (STR) mode to gamma-mode. Besides, our investigation further reveals that gamma-mode is more likely to dominate the discharge under high gas pressures or driving voltages.

Abstract: The formation process of a microdischarge (MD) in both μm- and nm-sized catalyst pores is simulated by a two-dimensional particle-in-cell/Monte Carlo collision model. A parallel-plate dielectric barrier discharge configuration in filamentary mode is considered in ambient air. The discharge is powered by a high voltage pulse. Our calculations reveal that a streamer can penetrate into the surface features of a porous catalyst and MDs can be formed inside both μm- and nm-sized pores, yielding ionization inside the pore. For the μm-sized pores, the ionization mainly occurs inside the pore, while for the nm-sized pores the ionization is strongest near and inside the pore. Thus, enhanced discharges near and inside the mesoporous catalyst are observed. Indeed, the maximum values of the electric field, ionization rate and electron density occur near and inside the pore. The maximum electric field and electron density inside the pore first increase when the pore size rises from 4 nm to 10 nm, and then they decrease for the 100 nm pore, due to
a more pronounced surface discharge for the smaller pores. However, the ionization rate is highest for the 100 nm pore due to the largest effective ionization region.

Abstract: A one-dimensional, implicit particle-in-cell Monte Carlo collision model is used to simulate the plasma kinetic properties at a steady state in a parallel-plate direct current argon glow microdischarge under various operating conditions, such as driving voltage (301000 V) and gap size (101000 µm) at atmospheric pressure. First, a comparison between rf and dc modes is shown for the same pressure, driving voltage and gap spacing. Furthermore, the effect of gap size scaling (in the range of 101000 µm) on the breakdown voltage, peak electron density and peak electron current density at the breakdown voltage is examined. The breakdown voltage is lower than 150 V in all gaps considered. The microdischarge is found to have a neutral bulk plasma region and a cathode sheath region with size varying with the applied voltage and the discharge gap. In our calculations, the electron and ion densities are of the order of 10181023 m−3, which is in the glow discharge limit, as the ionization degree is lower than 1% . The electron energy distribution function shows a two-energy group distribution at a gap of 10 µm and a three-energy group distribution at larger gaps such as 200 µm and 1000 µm, emphasizing the importance of the gap spacing in dc microdischarges.

Abstract: Precipitation is accompanied by the formation and migration of heterophase interfaces. Using the combined approach of advanced imaging and atomistic simulations, we studied the precipitate-matrix interfaces in various aluminium alloy systems, aiming to resolve their detailed atomic structures and illuminate their role in phase transformations.

Abstract: Nanoparticle modified electrodes constitute an attractive way to tailor-make efficient carbon dioxide (CO2) reduction catalysts. However, the restructuring and sintering processes of nanoparticles under electrochemical reaction conditions not only impedes the widespread application of nanoparticle catalysts, but also misleads the interpretation of the selectivity of the nanocatalysts. Here, we colloidally synthesized metallic copper (Cu) and silver (Ag) nanoparticles with a narrow size distribution (<10%) and utilized them in electrochemical CO2 reduction reactions. Monometallic Cu and Ag nanoparticle electrodes showed severe nanoparticle sintering already at low overpotential of -0.8 V vs. RHE, as evidenced by ex situ SEM investigations, and potential-dependent variations in product selectivity that resemble bulk Cu (14% for ethylene at -1.3 V vs. RHE) and Ag (69% for carbon monoxide at -1.0 V vs. RHE). However, by co-deposition of Cu and Ag nanoparticles, a nanoparticle stabilization effect was observed between Cu and Ag, and the sintering process was greatly suppressed at CO2 reducing potentials (-0.8 V vs. RHE). Furthermore, by varying the Cu/Ag nanoparticle ratio, the CO2 reduction reaction (CO2RR) selectivity towards methane (maximum of 20.6% for dense Cu-2.5-Ag-1 electrodes) and C-2 products (maximum of 15.7% for dense Cu-1-Ag-1 electrodes) can be tuned, which is attributed to a synergistic effect between neighbouring Ag and Cu nanoparticles. We attribute the stabilization of the nanoparticles to the positive enthalpies of Cu-Ag solid solutions, which prevents the dissolution-redeposition induced particle growth under CO2RR conditions. The observed nanoparticle stabilization effect enables the design and fabrication of active CO2 reduction nanocatalysts with high durability.

Abstract: Despite the natural abundance and promising properties of Si, there are few examples of crystalline Si-based catalysts. Here, the authors report an epitaxial growth method to construct Co single atoms on Si for light driven CO2 reduction to syngas. Improving the dispersion of active sites simultaneous with the efficient harvest of photons is a key priority for photocatalysis. Crystalline silicon is abundant on Earth and has a suitable bandgap. However, silicon-based photocatalysts combined with metal elements has proved challenging due to silicon's rigid crystal structure and high formation energy. Here we report a solid-state chemistry that produces crystalline silicon with well-dispersed Co atoms. Isolated Co sites in silicon are obtained through the in-situ formation of CoSi2 intermediate nanodomains that function as seeds, leading to the production of Co-incorporating silicon nanocrystals at the CoSi2/Si epitaxial interface. As a result, cobalt-on-silicon single-atom catalysts achieve an external quantum efficiency of 10% for CO2-to-syngas conversion, with CO and H-2 yields of 4.7 mol g((Co))(-1) and 4.4 mol g((Co))(-1), respectively. Moreover, the H-2/CO ratio is tunable between 0.8 and 2. This photocatalyst also achieves a corresponding turnover number of 2 x 10(4) for visible-light-driven CO2 reduction over 6 h, which is over ten times higher than previously reported single-atom photocatalysts.

Abstract: Direct hydroxylation via the functionalization of tertiary benzylic C(sp3)-H bond is of great significance for obtaining tertiary alcohols which find wide applications in pharmaceuticals as well as in fine chemical industries. However, current synthetic procedures use toxic reagents and therefore, the development of a sustainable strategy for the synthesis of tertiary benzyl alcohols is highly desirable. To solve this problem, herein, we report a metal-free
heterogeneous photocatalyst to synthesize the hydroxylated products using oxygen as the key reagent. Various benzylic substrates were employed into our mild reaction conditions to afford the desirable products in good to excellent yields. More importantly, gram-scale reaction was achieved via harvesting direct solar energy and exhibited high quantity of the product. The high stability of the catalyst was proved via recycling the catalyst and spectroscopic analyses. Finally, a possible mechanism was proposed based on the EPR and other experimental
evidence.

Abstract: Cryogenic electronmicroscopy can provide high-resolution reconstructionsof macromolecules embedded in a thin layer of ice from which atomicmodels can be built de novo. However, the interactionbetween the ionizing electron beam and the sample results in beam-inducedmotion and image distortion, which limit the attainable resolutions.Sample charging is one contributing factor of beam-induced motionsand image distortions, which is normally alleviated by including partof the supporting conducting film within the beam-exposed region.However, routine data collection schemes avoid strategies wherebythe beam is not in contact with the supporting film, whose rationaleis not fully understood. Here we characterize electrostatic chargingof vitreous samples, both in imaging and in diffraction mode. We mitigatesample charging by depositing a single layer of conductive grapheneon top of regular EM grids. We obtained high-resolution single-particleanalysis (SPA) reconstructions at 2 & ANGS; when the electron beamonly irradiates the middle of the hole on graphene-coated grids, usingdata collection schemes that previously failed to produce sub 3 & ANGS;reconstructions without the graphene layer. We also observe that theSPA data obtained with the graphene-coated grids exhibit a higher b factor and reduced particle movement compared to dataobtained without the graphene layer. This mitigation of charging couldhave broad implications for various EM techniques, including SPA andcryotomography, and for the study of radiation damage and the developmentof future sample carriers. Furthermore, it may facilitate the explorationof more dose-efficient, scanning transmission EM based SPA techniques.

Abstract: The addition of potassium thiocyanate (KSCN) to the FAPbBr3 structure and subsequent post-treatment of nanocrystals (NCs) lead to high quantum confinement, resulting in a photoluminescent quantum yield (PLQY) approaching unity and microsecond decay times. This synergistic approach demonstrated exceptional stability under humid conditions, retaining 70% of the PLQY for over a month, while the untreated NCs degrade within 24 h. Additionally, the devices incorporating the post-treated NCs displayed 1.5% external quantum efficiency (EQE), a 5-fold improvement over untreated devices. These results provide promising opportunities for the use of perovskites in moisture-stable optoelectronics.

Abstract: Pt nanoparticles play an essential role in a wide variety of catalytic reactions. The activity of the particles strongly depends on their three-dimensional (3D) structure and exposed facets, as well as on the reactive environment. High-resolution electron microscopy has often been used to characterize nanoparticle catalysts but unfortunately most observations so far have been either performed in vacuum and/or using conventional (2D) in situ microscopy. The latter however does not provide direct 3D morphological information. We have implemented a quantitative methodology to measure variations of the 3D atomic structure of nanoparticles under the flow of a selected gas. We were thereby able to quantify refaceting of Pt nanoparticles with atomic resolution during various oxidation−reduction cycles. In a H2 environment, a more faceted surface morphology of the particles was observed with {100} and {111} planes being dominant. On the other hand, in O2 the percentage of {100} and {111} facets decreased and a significant increase of higher order facets was found, resulting in a more rounded morphology. This methodology opens up new opportunities toward in situ characterization of catalytic nanoparticles because for the first time it enables one to directly measure 3D morphology variations at the atomic scale in a specific gaseous reaction environment.

Abstract: Plasma catalysis is gaining increasing interest for various environmental applications, but the crucial question is whether plasma can be created inside catalyst pores and under which conditions. In practice, various catalytic support materials are used, with various dielectric constants. We investigate here the influence of the dielectric constant on the plasma properties inside catalyst pores and in the sheath in front of the pores, for various pore sizes. The calculations are performed by a two-dimensional fluid model for an atmospheric pressure dielectric barrier discharge in helium. The electron impact ionization rate, electron temperature, electron and ion density, as well as the potential distribution and surface charge density, are analyzed for a better understanding of the discharge behavior inside catalyst pores. The results indicate that, in a 100 μm pore, the electron impact ionization in the pore, which is characteristic for the plasma generation inside the pore, is greatly enhanced for dielectric constants below 300. Smaller pore sizes only yield enhanced ionization for smaller dielectric constants, i.e., up to εr = 200, 150, and 50 for pore sizes of 50, 30, and 10 μm. Thus, the most common catalyst supports, i.e., Al2O3 and SiO2, which have dielectric constants around εr = 8−11 and 4.2, respectively, should allow more easily that microdischarges can be formed inside catalyst pores, even for smaller pore sizes. On the other hand, ferroelectric materials with dielectric constants above 300 never seem to yield plasma enhancement inside catalyst pores, not even for 100 μm pore sizes. Furthermore, it is clear that the dielectric constant of the material has a large effect on the extent of plasma enhancement inside the catalyst pores, especially in the range between εr = 4 and εr = 200. The obtained results are explained in detail based on the surface charge density at the pore walls,
and the potential distribution and electron temperature inside and above the pores. The results obtained with this model are
important for plasma catalysis, as the production plasma species in catalyst pores might affect the catalyst properties, and thus
improve the applications of plasma catalysis.

Abstract: Plasma catalysis is promising for various environmental, energy and chemical synthesis applications, but the underlying mechanisms are far from understood. Modeling can help to obtain a better insight in these mechanisms. Some burning questions relate to the plasma behavior inside packed bed reactors and whether plasma can penetrate into catalyst pores. In this paper, we try to provide answers to these questions, by means of both fluid modeling and particle-in-cell/Monte Carlo collision simulations. We present a short overview of recent findings obtained in our group by means of modeling, i.e., the enhanced electric field near the contact points and the streamer propagation through the packing in packed bed reactors, as well as the plasma behavior in catalyst pores, to determine the minimum pore size in which plasma streamers can penetrate.