Abstract: The dielectric properties of confined water is of fundamental interest and is still controversial. For water confined in channels with height smaller than h = 8 angstrom, we found a commensurability effect and an extraordinary decrease in the out-of-plane dielectric constant down to the limit of the dielectric constant of optical water. Spatial resolved polarization density data obtained from molecular dynamics simulations are found to be antisymmetric across the channel and are used as input in a mean-field model for the dielectric constant as a function of the height of the channel for h > 15 angstrom. Our results are in excellent agreement with a recent experiment [L. Fumagalli et al., Science 360, 1339 (2018)].

Abstract: Non-thermal plasma (NTP), an ionized gas generated at ambient pressure and temperature, has been an emerging technology for medical applications. Through controlled delivery of reactive oxygen and nitrogen species (ROS/RNS), NTP can elicit hormetic cellular responses, thus stimulating broad therapeutic effects. To enable clinical translation of the promising preclinical research into NTP therapy, a deeper understanding of NTP interactions with clinical substrates is profoundly needed. Since NTP-generated ROS/RNS will inevitably interact with blood in several clinical contexts, understanding their stability in this system is crucial. In this study, two medically relevant NTP delivery modalities were used to assess the stability of NTP-generated ROS/RNS in three aqueous solutions with increasing organic complexities: phosphate-buffered saline (PBS), blood plasma (BP), and processed whole blood. NTP-generated RNS collectively (NO2−, ONOO−), H2O2, and ONOO− exclusively were analyzed over time. We demonstrated that NTP-generated RNS and H2O2 were stable in PBS but scavenged by different components of the blood. While RNS remained stable in BP after initial scavenging effects, it was completely reduced in processed whole blood. On the other hand, H2O2 was completely scavenged in both liquids over time. Our previously developed luminescent probe europium(III) was used for precision measurement of ONOO− concentration. NTP-generated ONOO− was detected in all three liquids for up to at least 30 seconds, thus highlighting its therapeutic potential. Based on our results, we discussed the necessary considerations to choose the most optimal NTP modality for delivery of ROS/RNS to and via blood in the clinical context.

Abstract: Non-thermal plasma (NTP) therapy has been emerging as a promising cancer treatment strategy, and recently, its ability to locally induce immunogenic cancer cell death is being unraveled. We hypothesized that the chemical species produced by NTP reduce immunosuppressive surface proteins and checkpoints that are overexpressed on cancerous cells. Here, 3D in vitro tumor models, an in vivo mouse model, and molecular dynamics simulations are used to investigate the effect of NTP on CD47, a key innate immune checkpoint. CD47 is immediately modulated after NTP treatment and simulations reveal the potential oxidized salt-bridges responsible for conformational changes. Umbrella sampling simulations of CD47 with its receptor, signal-regulatory protein alpha (SIRPα), demonstrate that the induced-conformational changes reduce its binding affinity. Taken together, this work provides new insight into fundamental, chemical NTP-cancer cell interaction mechanisms and a previously overlooked advantage of present NTP cancer therapy: reducing immunosuppressive signals on the surface of cancer cells.

Abstract: Glioblastoma multiforme (GBM) is the most frequent and aggressive primary malignant brain tumor in adults. Current standard radiotherapy and adjuvant chemotherapy with the alkylating agent temozolomide (TMZ) yield poor clinical outcome. This is due to the stem-like properties of tumor cells and genetic abnormalities in GBM, which contribute to resistance to TMZ and progression. In this study, we used cold atmospheric plasma (CAP) to enhance the sensitivity to TMZ through inhibition of antioxidant signaling (linked to TMZ resistance). We demonstrate that CAP indeed enhances the cytotoxicity of TMZ by targeting the antioxidant specific glutathione (GSH)/glutathione peroxidase 4 (GPX4) signaling. We optimized the threshold concentration of TMZ on five different GBM cell lines (U251, LN18, LN229, U87-MG and T98G). We combined TMZ with CAP and tested it on both TMZ-sensitive (U251, LN18 and LN229) and TMZ-resistant (U87-MG and T98G) cell lines using two-dimensional cell cultures. Subsequently, we used a three-dimensional spheroid model for the U251 (TMZ-sensitive) and U87-MG and T98G (TMZ-resistant) cells. The sensitivity of TMZ was enhanced, i.e., higher cytotoxicity and spheroid shrinkage was obtained when TMZ and CAP were administered together. We attribute the anticancer properties to the release of intracellular reactive oxygen species, through inhibiting the GSH/GPX4 antioxidant machinery, which can lead to DNA damage. Overall, our findings suggest that the combination of CAP with TMZ is a promising combination therapy to enhance the efficacy of TMZ towards the treatment of GBM spheroids.

Abstract: CO2 dissociation and its subsequent conversion into added-value chemicals is a promising strategy for recycling CO2 gas into reusable products. One of the possible methods is direct plasma-induced dissociation. In this work we study the efficiency of CO2 dissociation in pulsed atmospheric-pressure gas discharge between two conducting electrodes by a 0-D numerical plasma model. The purpose of the study is to provide results on the optimal conditions of CO2 conversion with respect to the energy efficiency and dissociation by varying the maximum power density value and the pulse length. The power density is directly related to the discharge current and the reduced electric field in the discharge. We consider pulse lengths in the range from hundreds of nanosecond up to milliseconds. The results obtained show that the dissociation degree and energy efficiency are sensitive to the pulse length (duration) and the power density, so that a considerable improvement of the discharge performance can be achieved by fine-tuning these parameters. The study is intended to provide guidance in designing an experimental set-up and a power supply with the characteristics necessary to achieve optimal conversion.

Abstract: Targeting the redox balance of malignant cells via the delivery of high oxidative stress unlocks a potential therapeutic strategy against glioblastoma (GBM). We investigated a novel reactive oxygen species (ROS)-inducing combination treatment strategy, by increasing exogenous ROS via cold atmospheric plasma and inhibiting the endogenous protective antioxidant system via auranofin (AF), a thioredoxin reductase 1 (TrxR) inhibitor. The sequential combination treatment of AF and cold atmospheric plasma-treated PBS (pPBS), or AF and direct plasma application, resulted in a synergistic response in 2D and 3D GBM cell cultures, respectively. Differences in the baseline protein levels related to the antioxidant systems explained the cell-line-dependent sensitivity towards the combination treatment. The highest decrease of TrxR activity and GSH levels was observed after combination treatment of AF and pPBS when compared to AF and pPBS monotherapies. This combination also led to the highest accumulation of intracellular ROS. We confirmed a ROS-mediated response to the combination of AF and pPBS, which was able to induce distinct cell death mechanisms. On the one hand, an increase in caspase-3/7 activity, with an increase in the proportion of annexin V positive cells, indicates the induction of apoptosis in the GBM cells. On the other hand, lipid peroxidation and inhibition of cell death through an iron chelator suggest the involvement of ferroptosis in the GBM cell lines. Both cell death mechanisms induced by the combination of AF and pPBS resulted in a significant increase in danger signals (ecto-calreticulin, ATP and HMGB1) and dendritic cell maturation, indicating a potential increase in immunogenicity, although the phagocytotic capacity of dendritic cells was inhibited by AF. In vivo, sequential combination treatment of AF and cold atmospheric plasma both reduced tumor growth kinetics and prolonged survival in GBM-bearing mice. Thus, our study provides a novel therapeutic strategy for GBM to enhance the efficacy of oxidative stress-inducing therapy through a combination of AF and cold atmospheric plasma.

Abstract: Upgrading ethane with CO2 as a soft oxidant represents a desirable means of obtaining oxygenated hydrocarbons. This reaction is not thermodynamically feasible under mild conditions and has not been previously achieved as a one-step process. Nonthermal plasma was implemented as an alternative means of supplying energy to overcome activation barriers, leading to the production of alcohols, aldehydes, and acids as well as C1−C5+ hydrocarbons under ambient pressure, with a maximum total oxygenate selectivity of 12%. A plasma chemical kinetic computational model was developed and found to be in good agreement with the experimental trends. Results from this study illustrate the potential to use plasma for the direct synthesis of value-added alcohols, acids, and aldehydes from ethane and CO2 under mild conditions.

Abstract: Oxidative stress-inducing therapies are characterized as a specific treatment that involves the production of reactive oxygen and nitrogen species (RONS) by external or internal sources. To protect cells against oxidative stress, cells have evolved a strong antioxidant defense system to either prevent RONS formation or scavenge them. The maintenance of the redox balance ensures signal transduction, development, cell proliferation, regulation of the mechanisms of cell death, among others. Oxidative stress can beneficially be used to treat several diseases such as neurodegenerative disorders, heart disease, cancer, and other diseases by regulating the antioxidant system. Understanding the mechanisms of various endogenous antioxidant systems can increase the therapeutic efficacy of oxidative stress-based therapies, leading to clinical success in medical treatment. This review deals with the recent novel findings of various cellular endogenous antioxidant responses behind oxidative stress, highlighting their implication in various human diseases, such as ulcers, skin pathologies, oncology, and viral infections such as SARS-CoV-2.

Abstract: Biological membranes are under constant attack of reactive oxygen and nitrogen species (RONS), which may lead to a complex mixture of nitro-oxidized lipids that are responsible for structural and dynamic changes on the membrane. Because of that, nitro-oxidized lipids are also associated with several tumors and inflammatory and neurodegenerative diseases. Moreover, lipid oxidation may induce membrane phase-separated domains, which also drastically affect the membrane function. Evidence suggests that domain interfaces are “hot spots” for pore formation, but the underlying mechanisms remain elusive. There is an urgent need for an improved understanding of oxidation-induced phase separation on membrane properties. Likewise, the molecular structure at domain interfaces still needs to be elucidated. To evaluate the effect of lipid nitro-oxidation on the permeability of single-phase (homogeneous) and phase-separated (heterogeneous) phospholipid bilayers (PLBs), we performed atomistic molecular dynamics (MD) simulations using: (1) single-phase PLBs composed of several isomers of nitrated and/or oxidized lipids; (2) phase-separated PLBs composed of coexisting liquid ordered (Lo) and liquid disordered (Ld) domains, where the Ld domain is composed of non-oxidized and/or oxidized lipids. Our results show that nitrated lipids increase the membrane permeability of single-phase PLBs by three-fold compared to oxidized lipids. In addition, we show that oxidized lipids in the presence of nitrated lipids decrease the membrane permeability, suggesting an interaction between nitrated and oxidized lipids. Overall, the permeability of single-phase and phase-separated PLBs was comparable, and the presence of oxidized lipids increases the membrane permeability only in single-phase PLBs. Despite the latter, the presence of only 1.5% of lipid aldehydes at the Lo/Ld domain interfaces of phase-separated PLBs was able to increase the membrane permeability. In consequence of this, we also performed coarse-grained MD simulations to evaluate whether lipid aldehydes have a preference to accumulate at the interface between Lo/Ld domains. Our results show that lipid aldehydes derived from mono-unsaturated lipids accumulate at the interface, but those derived from poly-unsaturated lipids remain in the Ld domain. This study is of interest for photodynamic therapy and plasma medicine for cancer treatment, to understand the effects caused by RONS in cell membranes.

Abstract: 200 years ago, humanity started the industrial revolution by discovering fossil fuels, which lead to unprecedented technological advancements. However it has become alarmingly clear that the major environmental concerns associated with fossil fuels require a short-term transition from a carbon-based energy economy to a sustainable one based on green electricity. A key step concerning this transition exists in developing electricity-driven alternatives for chemical processes that rely on fossil fuels as a raw material. A technology that is gaining increasing interest to achieve this, is plasma technology. Using plasmas to induce chemical reactions by selectively heating electrons in a gas has already delivered promising results for gas conversion applications like CO2 conversion and N2 fixation, but plasma reactors still require optimization to be considered industrially competitive to existing fossil fuel-based processes and emerging other electricity-based technologies. In this thesis I develop computational models to describe plasma reactors and identify key mechanisms in three different plasma reactors for three different gas conversion applications, i.e. N2 fixation, combined CO2-CH4 conversion and CO2 splitting. I first developed models to describe a new rotating gliding arc (GA) reactor operating in two arc modes, which, as revealed by my model, are characterized by distinct plasma chemistry pathways. Subsequently, my colleague and I study the quenching effect of an effusion nozzle to this rotating GA reactor, reaching the best results to date for N2 fixation into NOx at atmospheric pressure, i.e., NOx concentrations up to 5.9%, at an energy cost down to 2.1 MJ/mol. Afterwards, I investigate the possible improvement of N2 admixtures in plasma-based CO2 and CH4 conversion, as significant amounts of N2 are often found in industrial CO2 waste streams, and gas separations are financially costly. Through combining my models with the experiment from a fellow PhD student, we reveal that moderate amounts of N2 (i.e. around 20%) increase both the electron density and the gas temperature to yield an overall energy cost reduction of 21%. Finally, I model quenching nozzles for plasma-based CO2 conversion in a microwave reactor, to explain the enhancements in CO2 conversion that were demonstrated in experiments. Through computational modelling I reveal that the nozzle introduces fast gas quenching resulting in the suppression of recombination reactions, which have more impact at low flow rates, where recombination is the most limiting factor in the conversion process.

Abstract: Natural deep eutectic solvents (NADES) show great promise as media for enzymatic reactions in areas where (bio)compatibility with natural or medicinal products is a must. While in theory they can be tailored to the intended reaction to ensure optimized yields, the knowledge to date is predominantly empirical, with some mechanistic reports providing a fragmented view at best. Therefore, it is not easy to explain experimental observations, let alone make predictions. The aim of this study was to develop a structured, holistic understanding of the effects of NADES media on enzymatic reactions, distinguishing between effects on solubility, solvation, viscosity, inhibition and denaturation. Experimental and computational chemistry methods were combined to separately study the interactions between enzyme, substrate, and NADES as reaction media. The initial enzyme activity and the final conversion of vinyl laurate transesterification by immobilized Candida antarctica lipase were studied experimentally. The direct effect of NADES on the same enzyme was modeled by molecular dynamics simulation. The effect of solubility was studied by both experimental and computational methods. To predict the solubility and viscosity of NADES, data-driven models were developed by combining group contribution and machine learning methods, based on the accumulated experimental knowledge on NADES found in the literature. Finally, the composed relationships and prediction models were applied to the practical example of deacetylation of mannosylerythritol lipids (MELs). The experimental findings show that the chosen NADES system has a significant effect on both the apparent initial activity and the final conversion. However, in the simulations, the enzyme retains its original structure; moreover, NADES has an additional stabilizing effect on the enzyme. In addition, changes in the molar ratio of the compounds in NADES do not show a significant effect on the stability of the enzyme. These results indicate that the main effect of NADES on the reaction is mainly related to the substrate-solvent interactions (solvation energy) and the viscosity of the system. On the other hand, the experimental results only confirmed the significance of solvation, viscosity did not show a clear correlation with the studied reaction parameters. The machine learning models built for solubility and viscosity gave quantitative predictions of these properties. The accumulated knowledge was used to optimize the yield in the deacetylation reaction of MELs. The combination of these methods provides fundamental knowledge about the effect of NADES on biocatalysis, but the results are also applicable to other uses of NADES.

Abstract: Low-temperature plasma has received extensive attention due to its promising application prospects in the field of air pollutants degradation and energy conversion. To fulfill the need for particular applications, constructing stable plasma sources and investigating the interaction mechanisms between plasma and substances have been hot research topics. This thesis reports the diagnosis and improvement of plasma sources, diagnosis of the active species in plasma and a modeling study of chemical kinetics processes. The main research contents are as follows: In Chapter 3, a diffuse sine AC dielectric barrier discharge (DBD) is successfully obtained by optimizing the electrode structure. It is found that using double-layer dielectric plates can limit the discharge current intensity and significantly improve the discharge uniformity. The electrical characteristics and gas temperature with different operating time show that the discharge stability is also improved by using double-layer dielectric plates. In Chapter 4, nanosecond pulses are employed to generate diffuse DBD plasmas. Three main discharge stages are distinguished by ICCD images, i.e., the streamer breakdown from the needle tip to the plate electrode, the regime transition from streamer to diffuse plasma, and the propagation of surface discharge on the plate electrode surface. The chapter reveales that in nanosecond pulsed discharges the vibrational temperature of N2 increases with the discharge duration, while the rotational temperature mainly stays constant, which means electron energy is transferred into the vibrational levels, but gas heating is not obvious during the discharge pulse. In Chapter 5, both sine AC DBD and nanosecond pulsed DBD, studied in Chapter 2 and 3, are used for formaldehyde degradation. It is found that nanosecond pulsed DBD has more homogenous characteristics, better stability, and lower plasma gas temperature. Moreover, the energy consumption of nanosecond pulsed DBD is much lower than that of AC DBD. In Chapter 6, a 0D chemical kinetics model is developed to investigate the underlying plasma chemistry of methane dry reforming in a nanosecond pulsed discharge. An overview of the dominant reaction pathways of CO2 and CH4 conversion into the major products is given. Furthermore, most of the CO2 molecules are populated into vibrational states during the pulse. Hence, the vibrational states of CO2 play an important role in its dissociation process. In general, this PhD thesis contributes to a better insight in the mechanisms of sinusoidal AC DBD and nanosecond pulsed DBD plasmas and their applications, i.e., decomposition of formaldehyde and dry reforming of methane.

Abstract: Sulfur hexafluoride (SF6), as a man-made gas, is widely used in power industry, semiconductor industry and metal-processing industry. However, SF6 is a greenhouse gas and its global warming potential is 23500 times that of CO2. Besides, SF6 is very stable, with a lifetime in the atmosphere for more than one thousand years. Under natural conditions, only the ultraviolet light can make it slowly decomposed. Thus, the emission of SF6 has a great threat to the environment. In recent years, with the development of our national economy, the use of SF6 increased dramatically. And 90% of the SF6 emissions come from the power industry. In the meantime, the emission of SF6 exists a ‘hysteresis effect’, as many of the SF6-gas insulation equipment will retire in next decades, the emission of SF6 may increase sharply, and this may put great pressure on the environment. Therefore, it’s necessary to make efforts in controlling and treating the SF6 emission. Among the SF6 abatement technologies, the non-thermal plasma(NTP) represented by the dielectric barrier discharge(DBD) can effectively degrade SF6 and is suitable for large-scale industry applications. However, its energy efficiency still gets room for improvement and this kind of method has a defect that it’s hard to regulate the degradation by-products. Therefore, this paper proposed the combination of the packed bed reactor and the DBD technology to form a packed DBD discharge system for SF6 degradation, so that to further improve the energy efficiency and regulate the selectivity of by-products. By experiment and simulation research, the following innovations have been achieved: (1) Based on the packed bed DBD platform, the power parameter and gas-phase parameters of SF6 degradation were studied. It was found that the discharge process was significantly enhanced with the addition of packing particles, and the discharge energy efficiency was improved. The increase of input voltage can obviously increase the degradation rate, but reduces the energy efficiency. The increase of SF6 initial concentration and gas flow rate can improve the energy efficiency, but reduce the degradation rate. Therefore, both degradation rate and energy efficiency should be considered in deciding basic experimental conditions. (2) Active gases, such as O2, H2O and NH3, could effectively promote the degradation rate of SF6, and changed the product selectivity. In our packed bed DBD system, O2 and H2O have the optimal concentration conditions, which are 2% and 1%, respectively. The addition of O2 can promote the generation of S-O-F products, and inhibit the selectivity of SO2, while the addition of H2O had the opposite effects. In addition, the synergistic degradation of NH3 and SF6 will produce solid products, such as NH3HF, NH4HF2 and elemental S. For gaseous products, the increase of NH3 will lead to the generation of SO2 in the final degradation products and inhibit the generation of S-O-F products. (3) Different kinds of packing materials have great impacts on the degradation system in the discharge parameters, degradation rate and energy efficiency, as well as the products distribution. In the experiment, we compared the degradation results in three systems: glass beads packing, γ-Al2O3 packing and no-packing system. The packing of glass beads effectively improved the discharge voltage amplitude and discharge power, while had a limited effect on the equivalent capacitance of the dielectric. Besides, γ-Al2O3 packing had little effect on voltage amplitude, but obviously increased the equivalent capacitance of the dielectric. Furthermore, the degradation rate and energy efficiency in γ-Al2O3 system was higher than that of glass bead system. For products selectivity, γ-Al2O3 system was more desirable, where S-O-F type of product selectivity was suppressed and the SO2 selectivity increased significantly. By contrast, the glass beads system hardly affected the product selectivity. This results are presumably due to the relatively high dielectric constant of γ-Al2O3 particles and γ-Al2O3 itself may act as a reactant or a catalyst participating in the degradation reactions. (4) The size and status of the packing particles also have significant effects on the degradation process. The systems packed with 1, 2 and 4mm γ-Al2O3 particles for SF6 degradation were compared, and the 2mm system had the best performance, which may because the 2mm system had a good balance between the active contact area and the gas residence time. In addition, the packing pellets suffered from a hydration process slightly reduced the discharge parameters in the γ-Al2O3 packing system and significantly reduced the degradation rate was, which may because the H2O molecules pre-occupied the active sites on the γ-Al2O3 surface and reduced the discharge process. (5) Based on density functional theory (DFT), the degradation process of SF6 in the packed bed DBD system was studied at atomic scale. It was found that the SF6 can occur a physical adsorption at AlⅢ active sites on γ-Al2O3 surface. The activation barrier for the first degradation step of SF6 on γ-Al2O3 surface is much lower than in gas phase, which proved that the SF6 molecule is activated on the γ-Al2O3 surface. In addition, the plasma may affect the γ-Al2O3 surface to generate excess electrons or external electric fields. This two effects can change the adsorbed SF6 molecules from physical adsorption to chemisorption, together with an obvious stretching of S-F bonds, indicating that the plasma surface effects prmote the activation and decomposition of SF6 molecules. Furthermore, the stepwise degradation process of SF6 on γ-Al2O3 surface were investigated. The influence of radicals produced by plasma on the degradation process was analyzed. It was found that via Eley–Rideal (ER) reactions, high-energy radicals could effectively reduce the activation barriers and promote the surface reactions. Finally, the degradation mechanism of SF6 molecules in the packed bed plasma system was summarized, which may provide a theoretical basis for the study of harmless degradation of SF6. Keywords: SF6; Packed Bed DBD; Discharge Parameters; Products Analysis; Degradation Mechanism

Abstract: The low-pressure hollow cathode discharge made of a hollow circular tube and an anode is a type of simple structure discharge system. In particular, under the arc discharge mode, hollow cathodes have high plasma density and energy density with a wide range of adaptability of pressure and current. Low-pressure hollow cathode arc (HCA) discharges have been widely used as plasma sources in various fields such as manufacturing, vacuum welding, and aerospace since the 1960s. Despite the early experimental and applied researches on low-pressure HCA discharges, the basic theoretical study was relatively lagged much behind, resulting in many unanswered questions, such as the optimal discharge operating parameters, the power deposition inside the cathode, the causes of plasma instability, and how to effectively reduce cathode erosion and so on. Due to the special discharge structure of the hollow cathode, it is difficult to make an accurate experimental diagnosis, so a reasonable numerical simulation is an effective study method. However, up to now, there is still a lack of complete and effective numerical models which can evaluate various physical fields in the low-pressure hollow cathode discharges. To address the above problems and difficulties, a comprehensive and self-consistent 2D multi-physical coupling numerical model based on a commercial program of finite element method, the COMSOL Multiphysics, was provided in this paper. The model involves plasma transport, arc flow and heat transfer, and cathode thermal equilibrium, and can consider the effect of an applied magnetic field. The processes of secondary electron emission, thermal-field electron emission, ions and backflow high-energy electrons bombardment, and thermal radiation from the cathode surface are considered in the cathode thermal equilibrium process. Based on the above background, this paper works from the following aspects： In Chapter 1, the basic concepts of low-pressure HCA discharge including the hollow cathode effect, the basic characteristics, and operation modes were introduced firstly; Secondly, the application fields, development history, and overseas and domestic research status of hollow cathode discharge were reviewed; finally, the problems were presented and the research background was explained, and the research purpose of this paper was clarified. In Chapter 2, a complete and self-consistent numerical model of low-pressure hollow cathode discharge was proposed based on the fundamental theory and assumptions, and the set of control equations and boundary conditions in the model were elaborated. In addition, the electron energy distribution function, the collision processes, the solving tools of this model, and calculation schemes were introduced in detail. Finally, a validation example was given to test the rationality and applicability of the numerical model. In Chapter 3, the fundamental plasma properties of low-pressure hollow cathode arcs were investigated. Firstly, the ion Joule heating effect was studied. The results showed that the temperature distributions of the arc and cathode are only able to approach the experimental measurements after considering the ion Joule heating, which shows that the Joule heating of ions is crucial for the heating of the arc plasma. Secondly, by comparing the radial distribution of electron and ion density inside the cathode, the structure of the cathode sheath could be simulated well using this model. Finally, it was shown that the thermal radiation from the cathode surface is an important cooling mechanism of the cathode and only under higher surface emissivity can balance the larger heat flow given by the plasma to the cathode, and the temperature distribution of the cathode shows a non-monotonic increasing trend and is consistent with the profile of experimental measurement so that the so-called active zone is formed. In Chapter 4, the power deposition in the low-pressure HCA was studied in simulation. Two main aspects were considered: the power deposition into particles (both electrons and heavy particles) and the power deposition onto the cathode. It was found that the deposited power into particles increases with the rise of discharge current, but there is no effect on the total power deposition onto the cathode. In high-density plasmas, Coulomb collisions between electrons and ions also become very important, especially since a portion of the deposition energy on heavy particles comes mainly from the energy transfer from electrons to ions. It was also found that regardless of external parameters, half of the power deposition onto the cathode always comes from the particle contribution, while the other half is the net contribution of heat transfer and cathode radiation. The HCA model also allows the simulation of multiple discharge modes for low-pressure HCA discharges over a wide range of gas flow rates. It was also shown that the discharge operating conditions and the external magnetic field can change the distribution of the particle flow on the cathode wall. In Chapter 5, the ion sputtering erosion process on the cathode was simulated by coupling the HCA numerical model with the moving grid technique. The results showed that the ion sputtering erosion on the cathode depends on the ion flux and the plasma potential near the cathode wall and that their distribution and magnitude jointly determine the erosion morphology of the cathode. It was also found that the location of the most severe erosion on the cathode is located in the region of the densest ion flux on the cathode wall, rather than in the longitudinal correspondence with the central region of the internal positive column (IPC). The external magnetic fields can mitigate the cathode erosion and reduce the erosion depth, but stronger magnetic fields lead to a concentration of current density at the cathode tip, which can enhance erosion slightly at the cathode outlet end. Finally, the conclusions and innovation highlights were summarized, and prospects for future work were discussed.