Abstract: In solids, charge polarity can one-to-one correspond to spin polarity phenomenologically, e.g., ferroelectricity/ferromagnetism, antiferroelectricity/antiferromagnetism, and even dipole-vortex/magnetic-vortex, but ferrielectricity/ferrimagnetism kept telling a disparate story in microscopic level. Since the definition of a charge dipole involves more than one ion, there may be multiple choices for a dipole unit, which makes most ferrielectric orders equivalent to ferroelectric ones, i.e., this ferrielectricity is not necessary to be a real independent branch of polarity. In this work, by using the spherical aberration-corrected scanning transmission electron microscope, we visualize a nontrivial ferrielectric structural evolution in BaFe2Se3, in which the development of two polar sub-lattices is out-of-sync, for which we term it as irreducible ferrielectricity. Such irreducible ferrielectricity leads to a non-monotonic behavior for the temperature-dependent polarization, and even a compensation point in the ordered state. Our finding unambiguously distinguishes ferrielectrics from ferroelectrics in solids.

Abstract: Understanding the molecular-level mechanisms of phase transformation in solids is of fundamental interest for functional materials such as zeolites. Two-dimensional (2D) zeolites, when used as shape-selective catalysts, can offer improved access to the catalytically active sites and a shortened diffusion length in comparison with their 3D analogues. However, few materials are known to maintain both their intralayer microporosity and structure during calcination for organic structure-directing agent (SDA) removal. Herein we report that PST-9, a new 2D zeolite which has been synthesized via the multiple inorganic cation approach and fulfills the requirements for true layered zeolites, can be transformed into the small-pore zeolite EU-12 under its crystallization conditions through the single-layer folding process, but not through the traditional dissolution/recrystallization route. We also show that zeolite crystal growth pathway can differ according to the type of organic SDAs employed.

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: In traditional capacitively coupled plasmas, the discharge can be described by an electrostatic model, in which the Poisson equation is employed to determine the electrostatic electric field. However, current plasma reactors are much larger and driven at a much higher frequency. If the excitation wavelength k in the plasma becomes comparable to the electrode radius, and the plasma skin depth d becomes comparable to the electrode spacing, the electromagnetic (EM) effects will become significant and compromise the plasma uniformity. In this regime, capacitive discharges have to be described by an EM model, i.e., the full set of Maxwells equations should be solved to address the EM effects. This paper gives an overview of the theory, simulation and experiments that have recently been carried out to understand these effects, which cause major uniformity problems in plasma processing for microelectronics and flat panel display industries. Furthermore, some methods for improving the plasma uniformity are also described and compared.

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: Oriented Co nanoparticles were obtained by Co ion implantation in crystalline ZnO thin films grown by pulsed laser deposition. Transmission electron microscopy revealed the presence of elliptically shaped Co precipitates with nanometer size, which are embedded in the ZnO thin films, resulting in anisotropic magnetic behavior. The low-temperature resistance of the Co-implanted ZnO thin films follows the Efros-Shklovskii type variable-range-hopping. Large negative magnetoresistance (MR) exceeding 10% is observed in a magnetic field of 1 T at 2.5K and the negative MR survives up to 250K (0.3%). The negative MR reveals hysteresis as well as anisotropy that correlate well with the magnetic properties, clearly demonstrating the presence of spin-dependent transport. (C) 2013 AIP Publishing LLC.

Abstract: We apply an implicit particle-in-cell Monte-Carlo (PIC-MC) method to study a radio-frequency argon microdischarge at steady state in the glow discharge limit, in which the microdischarge is sustained by secondary electron emission from the electrodes. The plasma density, electron energy distribution function (EEDF), and electron temperature are calculated in a wide range of operating conditions, including driving voltage, microdischarge gap, and pressure. Also, the effect of gap size scaling (in the range of 50-1000 μm) on the plasma sustaining voltage and peak electron density at atmospheric pressure is examined, which has not been explored before. In our simulations, three different EEDFs, i.e., a so-called three temperature hybrid mode, a two temperature α mode, and a two temperature γ mode distribution, are identified at different gaps and voltages. The maximum sustaining voltage to avoid a transition from the glow mode to an arc is predicted, as well as the minimum sustaining voltage for a steady glow discharge. Our calculations elucidate that secondary electrons play an essential role in sustaining the discharge, and as a result the relationship between breakdown voltage and gap spacing is far away from the Paschen law at atmospheric pressure.

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: In this paper, a self-consistent fluid model is developed focusing on the plasma parameters in capacitively coupled 20% N 2-80% Ar discharges. Measurements of ion density are performed with the help of a floating double probe, and the emission intensities from Ar(4p) and N 2 ( B ) transitions are detected by an optical emission spectroscopy to estimate their relative densities. The consistency between the numerical and experimental results confirms the reliability of the simulation. Then the plasma characteristics, specifically the reaction mechanisms of active species, are analyzed under various voltages. The increasing voltage leads to a monotonous increase in species density, whereas a less homogeneous radial distribution is observed at a higher voltage. Due to the high concentration of Ar gas, Ar + becomes the main ion, followed by the N 2 +</mml:msubsup> ion. Besides the electron impact ionization of neutrals, the charge transfer processes of Ar +/ N 2 and N 2 +</mml:msubsup>/Ar are found to have an impact on the ionic species. The results indicate that adopting the lower charge transfer reaction rate coefficients weakens the Ar + ion density and yields a higher N 2 +</mml:msubsup> ion density. However, the effect on the species spatial distributions and other species densities is limited. As for the excited-state species, the electron impact excitation of background gases remains overwhelming in the formation of Ar(4p), N 2 ( B ), and N 2 ( a ' ), whereas the <mml:msub> N 2 ( A ) molecules are mainly formed by the decay of <mml:msub> N 2 ( B ). In addition, the dissociation of <mml:msub> N 2 collided by excited-state Ar atoms dominates the N generation, which are mostly depleted to produce N + ions.

Abstract: A two-dimensional self-consistent fluid model coupled with the full set of Maxwell equations is established to investigate the phase-shift effect on the transient behaviour of electrodynamics and power deposition in a hydrogen capacitively coupled plasma. The effect has been examined at 13.56 MHz and 100 MHz, respectively, because of the different phase-shift modulation when the electromagnetic effects are dominant. The results indicate that the spatiotemporal distributions of the plasma characteristics obtained for various phase-shift cases are obviously different both in shape and especially in absolute values. Indeed, when the phase difference varies from 0 to π, there is an increase in the electron flux, thus the power deposition becomes more pronounced. At the frequency of 13.56 MHz, the axial electron flux in the bulk plasma becomes uniform along the z-axis, and the radial electron flux exhibits two peaks within one period at the reverse-phase case, whereas the oscillation is less pronounced at the in-phase case. Furthermore, in the very high frequency discharge, the radial electron flux is alternately positive and negative with four peaks during one period, and the ionization mainly occurs in the sheath region, due to the prominent power deposition there at a phase difference equal to π.

Abstract: A two-dimensional self-consistent fluid model combined with the full set of Maxwell equations is employed to investigate an Ar/CF4 capacitively coupled plasma, focusing on the phase-shift effect on the plasma characteristics at various frequencies and gas mixture ratios. When the discharge is sustained by a single frequency at 13.56 MHz in an Ar/CF4 mixture with a ratio of 0.9/0.1, no obvious difference is detected between the electron densities obtained in the so-called electrostatic model (with only the static electric fields taken into account) and the electromagnetic model (which includes the electromagnetic effects). However, as the frequency increases to 60 and 100 MHz, the difference becomes distinct, due to the significant influence of the electromagnetic effects. The phase-shift effect on the plasma radial uniformity has also been investigated in a dual frequency discharge, i.e. when the top driven source is switched on with a phase difference phiv ranging from 0 to π, in the frequency range 13.56100 MHz. At low concentration of CF4 (10%), Ar+ ions are the major positive ions in the entire range of frequencies. When the frequency is low, i.e. 13.56 MHz, the Ar+ density exhibits an off-axis peak at phiv = 0 due to the edge effect, and a better uniformity caused by the phase-shift modulation is obtained at phiv = π. At 60 MHz, the Ar+ density varies from edge-peaked at phiv = 0 to uniform (i.e. at phiv = 0.53π), and finally at phiv = π, a broad maximum is observed at the centre due to the standing-wave effect. As the frequency increases to 100 MHz, the best radial uniformity is reached at 0.25π, and the maximum moves again towards the radial wall in the reverse-phase case (phiv = π) due to the dominant skin effect. When the frequency is fixed at 100 MHz, the phase-shift control shows a different behaviour at a high concentration of CF4. For instance, the ${\rm CF}_3

Abstract: Compared with dust-free plasmas, the existence of dust particles in plasmas may greatly influence the plasma properties. such as the plasma density, electron temperature, sheath properties, electron energy distribution function (EEDF) as well as the heating mechanism. In this work, a 1D hybrid fluid/MC model has been developed to investigate the interaction between dust and plasma in a low-pressure silane discharge sustained in a radio frequency capacitively coupled plasma, in which we assume spherical dust particles with a given radius are generated by taking the sum of the production rate of Si2H4- and Si2H5- as the nucleation rate. From our simulation, the plasma may experience definite perturbation by dust particles with a certain radius (more than 50nm) with an increase in electron temperature first, which further induces a rapid rise in the positive and negative ion densities. Then, the densities begin to decline due to the gradual lack of sufficient seed electrons. In addition, as the dust radius increases, the high energy tails of the EEDFs will be enhanced for discharge maintenance, accompanied by a decline in the population of low-energy electrons in comparison with those of pristine plasma. Furthermore, an obvious bulk heating is observed apart from the a-mode and local field reversal heating. This may contribute to the enhanced bulk electric field (also called the drift field) as a result of electron depletion via the dust. In addition, large-sized dust particles that accumulate near the sheaths tend to form two stable density peaks with their positions largely influenced by the time-averaged sheath thickness. A detailed study of the effects of the external parameters, including pressure, voltage and frequency, on the spatial distribution of dust particles is also conducted.

Abstract: A systematic investigation of voltage effect on the plasma parameters, especially the species densities and chemical reaction mechanisms, in the capacitive N-2 discharges is performed by employing a two-dimensional self-consistent fluid model. The validity of the numerical model is first demonstrated by the qualitative agreement of the calculated and experimental results. Then, the densities, production mechanisms, and loss mechanisms of species from simulation are examined at various voltages. It is found that all the species densities increase monotonically with the voltage, whereas their spatial profiles at lower voltages are quite different from those at higher voltages. The electrons and Nthorn 2 ions are mainly generated by the electron impact ionization of N-2 gas, while the Nthorn ions, whose density is one or two orders of magnitude lower, are mostly formed by the ionization of N atoms. The electron impact dissociation of N-2 gas dominates the generation of N atoms, which are mostly destroyed for the Nthorn ion production. As for the excited N-2 levels, the level conversion processes play a very important role in their production and depletion mechanisms, except for the electron impact excitation of the ground state N-2 molecules. Published under license by AIP Publishing.

Abstract: Plasma generation inside catalyst pores is of utmost importance for plasma catalysis, as the existence of plasma species inside the pores affects the active surface area of the catalyst available to the plasma species for catalytic reactions. In this paper, the electric field enhancement, and thus the plasma production inside catalyst pores with different pore shapes is studied with a two-dimensional fluid model. The results indicate that the electric field will be significantly enhanced near tip-like structures. In a conical pore with small opening, the strongest electric field appears at the opening and bottom corners of the pore, giving rise to a prominent ionization rate throughout the pore. For a cylindrical pore, the electric field is only enhanced at the bottom corners of the pore, with lower absolute value, and thus the ionization rate inside the pore is only slightly enhanced. Finally, in a conical pore with large opening, the electric field is characterized by a maximum at the bottom of the pore, yielding a similar behavior for the ionization rate. These results demonstrate that the shape of the pore has a significantly influence on the electric field enhancement, and thus modifies the plasma properties.

Abstract: Plasma medicine has attracted tremendous interest in a variety of medical conditions, ranging from wound healing to antimicrobial applications, even in cancer treatment, through the interactions of cold atmospheric plasma (CAP) and various biological tissues directly or indirectly. The underlying mechanisms of CAP treatment are still poorly understood although the oxidative effects of CAP with amino acids, peptides, and proteins have been explored experimentally. In this study, machine learning (ML) technology is introduced to efficiently unveil the interaction mechanisms of amino acids and reactive oxygen species (ROS) in seconds based on the data obtained from the reactive molecular dynamics (MD) simulations, which are performed to probe the interaction of five types of amino acids with various ROS on the timescale of hundreds of picoseconds but with the huge computational load of several days. The oxidative reactions typically start with H-abstraction, and the details of the breaking and formation of chemical bonds are revealed; the modification types, such as nitrosylation, hydroxylation, and carbonylation, can be observed. The dose effects of ROS are also investigated by varying the number of ROS in the simulation box, indicating agreement with the experimental observation. To overcome the limits of timescales and the size of molecular systems in reactive MD simulations, a deep neural network (DNN) with five hidden layers is constructed according to the reaction data and employed to predict the type of oxidative modification and the probability of occurrence only in seconds as the dose of ROS varies. The well-trained DNN can effectively and accurately predict the oxidative processes and productions, which greatly improves the computational efficiency by almost ten orders of magnitude compared with the reactive MD simulation. This study shows the great potential of ML technology to efficiently unveil the underpinning mechanisms in plasma medicine based on the data from reactive MD simulations or experimental measurements. In this study, since reactive molecular dynamics simulation can currently only describe interactions between a few hundred atoms in a few hundred picoseconds, deep neural networks (DNN) are introduced to enhance the simulation results by predicting more data efficiently. image