Abstract: Recently, Sr-doped LaCrO3 has been experimentally introduced as a new p-type transparent conducting oxide. It is demonstrated that substituting Sr for La results in inducing p-type conductivity in LaCrO3. Performing first principles calculations we study the electronic structure and formation energy of various point defects in LaCrO3. Our results for the formation energies show that in addition to Sr, two more divalent defects, Ca and Ba, substituting for La in LaCrO3, behave as shallow acceptors in line with previous experimental reports. We further demonstrate that under oxygen-poor growth conditions, these shallow acceptors will be compensated by intrinsic donor-like defects (an oxygen vacancy and Cr on an oxygen site), but in the oxygen-rich growth regime the shallow acceptors have the lowest formation energies between all considered defects and will lead to p-type conductivity.

Abstract: The vast biomedical potential of cold atmospheric pressure plasmas (CAPs) is governed by the formation of reactive species. These biologically active species are formed upon the interaction of CAPs with the surroundings. In biological milieu, water plays an essential role. The development of biomedical CAPs thus requires understanding of the sources of the reactive species in aqueous media exposed to the plasma. This is especially important in case of the COST RF plasma jet, which is developed as a reference microplasma system. In this work, we investigated the formation of the OH radicals, H atoms and H2O2 in aqueous solutions exposed to the COST plasma jet. This was done by combining experimental and modelling approaches. The liquid phase species were analysed using UV-Vis spectroscopy and spin trapping with hydrogen isotopes and electron paramagnetic resonance (EPR) spectroscopy. The discrimination between the species formed from the liquid phase and the gas phase molecules was performed by EPR and 1H-NMR analyses of the liquid samples. The concentrations of the reactive species in the gas phase plasma were obtained using a zero-dimensional (0D) chemical kinetics computational model. A three-dimensional (3D) fluid dynamics model was developed to provide information on the induced humidity in the plasma effluent. The comparison of the experimentally obtained trends for the formation of the species as a function of the feed gas and effluent humidity with the modelling results suggest that all reactive species detected in our system are mostly formed in the gas phase plasma inside the COST jet, with minor amounts arising from the plasma effluent humidity.

Abstract: A number of recent computational material design studies based on density functional theory (DFT) calculations have put forward a new class of materials with electrically switchable chemical characteristics that can be exploited in the development of tunable gas storage and electrocatalytic applications. We find systematic flaws in almost every computational study of gas adsorption on polarized or charged surfaces, stemming from an improper and unreproducible treatment of periodicity, leading to very large errors of up to 3 eV in some cases. Two simple corrective procedures that lead to consistent results are proposed, constituting a crucial course correction to the research in the field.

Abstract: The thermodynamic aspects of various 2D materials are explored using Density Functional Theory (DFT). Various metal chalcogenides (MX2, M = metal, chalcogen X = S, Se, Te) are investigated with respect to their interaction and stability under different ambient conditions met in the integration process of a transistor device. Their interaction with high- dielectrics is also addressed, in order to assess their possible integration in Complementary Metal Oxide Semiconductor (CMOS) field effect transistors. 2D materials show promise for high performance nanoelectronic devices, but the presence of defects (vacancies, grain boundaries,...) can significantly impact their electronic properties. To assess the impact of defects, their enthalpies of formation and their signature levels in the density of states have been studied. We find, consistently with literature reports, that chalcogen vacancies are the most likely source of defects. It is shown that while pristine 2D materials are in general stable whenever set in contact with different ambient atmospheres, the presence of defective sites affects the electronic properties of the 2D materials to varying degrees. We observe that all the 2D materials studied in the present work show strong reactivity towards radical oxygen plasma treatments while reactivity towards other common gas phase chemical such as O-2 and H2O and groups present at the high- surface varies significantly between species. While energy band-gaps, effective masses and contact resistivities are key criteria in selection of 2D materials for scaled CMOS and tunneling based devices, the phase and ambient stabilities might also play a very important role in the development of reliable nanoelectronic applications.

Abstract: The mechanisms of plasma in medicine are broadly attributed to plasma-derived reactive oxygen and nitrogen species (RONS). In order to exert any intracellular effects, these plasma-derived RONS must first traverse a major barrier in the cell membrane. The cell membrane lipid composition, and thereby the magnitude of this barrier, is highly variable between cells depending on type and state (e.g. it is widely accepted that healthy and cancerous cells have different membrane lipid compositions). In this study, we investigate how plasma-derived RONS interactions with lipid membrane components can potentially be exploited in the future for treatment of diseases. We couple phospholipid vesicle experiments, used as simple cell models, with molecular dynamics (MD) simulations of the lipid membrane to provide new insights into how the interplay between phospholipids and cholesterol may influence the response of healthy and diseased cell membranes to plasma-derived RONS. We focus on the (i) lipid tail saturation degree, (ii) lipid head group type, and (iii) membrane cholesterol fraction. Using encapsulated molecular probes, we study the influence of the above membrane components on the ingress of RONS into the vesicles, and subsequent DNA damage. Our results indicate that all of the above membrane components can enhance or suppress RONS uptake, depending on their relative concentration within the membrane. Further, we show that higher RONS uptake into the vesicles does not always correlate with increased DNA damage, which is attributed to ROS reactivity and lifetime. The MD simulations indicate the multifactorial chemical and physical processes at play, including (i) lipid oxidation, (ii) lipid packing, and (iii) lipid rafts formation. The methods and findings presented here provide a platform of knowledge that could be leveraged in the development of therapies relying on the action of plasma, in which the cell membrane and oxidative stress response in cells is targeted.

Abstract: While spectroscopic data on small hydrocarbons in interstellar media in combination with crossed molecular beam (CMB) experiments have provided a wealth of information on astrochemically relevant species, much of the underlying mechanistic pathways of their formation remain elusive. Therefore, in this work, the chemical reaction mechanisms of C(³P_J) + C₆H₆and C⁺(²P) + C₆H₆systems using the quantum mechanical molecular dynamics (QMMD) technique at the PBE0-D3(BJ) level of theory is investigated, mimicking a CMB experiment. Both the dynamics of the reactions as well as the electronic structure for the purpose of the reaction network are evaluated. The method is validated for the first reaction by comparison to the available experimental data. The reaction scheme for the C(³P_J) + C₆H₆system covers the literature data,<italic>e.g.</italic>the major products are the 1,2-didehydrocycloheptatrienyl radical (C₇H₅) and benzocyclopropenyl radical (C₆H₅–CH), and it reveals the existence of less common pathways for the first time. The chemistry of the C⁺(²P_J) + C₆H₆system is found to be much richer, and we have found that this is because of more exothermic reactions in this system in comparison to those in the C(³P_J) + C₆H₆system. Moreover, using the QMMD simulation, a number of reaction paths have been revealed that produce three distinct classes of reaction products with different ring sizes. All in all, at all the collision energies and orientations, the major product is the heptagon molecular ion for the ionic system. It is also revealed that the collision orientation has a dominant effect on the reaction products in both systems, while the collision energy mostly affects the charged system. These simulations both prove the applicability of this approach to simulate crossed molecular beams, and provide fundamental information on reactions relevant for the interstellar medium.

Abstract: Electron acceleration from the interaction of an intense short-pulse laser with low density plasma is considered. The relation between direct electron acceleration within the laser pulse and that in the wake is investigated analytically. The magnitude and location of the ponderomotive-force-caused charge separation field with respect to that of the pulse determine the relative effectiveness of the two acceleration mechanisms. It is shown that there is an optimum condition for acceleration in the wake. Electron acceleration within the pulse dominates as the pulse becomes sufficiently short, and the latter directly drives and even traps the electrons. The latter can reach ultrahigh energies and can be extracted by impinging the pulse on a solid target. (C) 2003 American Institute of Physics.

Abstract: A fluid model is self-consistently established to investigate the harmonic effects in an inductively coupled plasma, where the electromagnetic field is solved by the finite difference time domain technique. The spatiotemporal distribution of harmonic current density, harmonic potential, and other plasma quantities, such as radio frequency power deposition, plasma density, and electron temperature, have been investigated. Distinct differences in current density have been observed when calculated with and without Lorentz force, which indicates that the nonlinear Lorentz force plays an important role in the harmonic effects, especially at low frequencies. Moreover, the even harmonics are larger than the odd harmonics both in the current density and the potential. Finally, the dependence of various plasma quantities with and without the Lorentz force on various driving frequencies is also examined. It is shown that the deposited power density decreases and the depth of penetration increases slightly because of the Lorentz force. The electron density increases distinctly while the electron temperature remains almost the same when the Lorentz force is taken into account.

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: Based on a rational classification of defects in amorphous materials, we propose a simplified model to describe intrinsic defects and hydrogen impurities in amorphous indium gallium zinc oxide (a-IGZO). The proposed approach consists of organizing defects into two categories: point defects, generating structural anomalies such as metal-metal or oxygen-oxygen bonds, and defects emerging from changes in the material stoichiometry, such as vacancies and interstitial atoms. Based on first-principles simulations, it is argued that the defects originating from the second group always act as perfect donors or perfect acceptors. This classification simplifies and rationalizes the nature of defects in amorphous phases. In a-IGZO, the most important point defects are metal-metal bonds (or small metal clusters) and peroxides (O-O single bonds). Electrons are captured by metal-metal bonds and released by the formation of peroxides. The presence of hydrogen can lead to two additional types of defects: metal-hydrogen defects, acting as acceptors, and oxygen-hydrogen defects, acting as donors. The impact of these defects is linked to different instabilities observed in a-IGZO. Specifically, the diffusion of hydrogen and oxygen is connected to positive-and negative-bias stresses, while negative-bias illumination stress originates from the formation of peroxides.

Abstract: We performed a detailed investigation of the oxidation of the technologically relevant GaAs(001)-beta 2(2x4) surface via density functional calculations. The purpose is to gain insights on the atomistic mechanisms and local bondings that underlie the degradation of the surface properties once exposed to oxygen. The study comprises the adsorption of single O atoms, through the sampling of several adsorption sites, and the subsequent formation of the O adsorbate at increasing coverage by taking into account multiple-atom adsorption. Based on the evaluation of the energetics and the structural properties of the atomistic models generated, the results here reported delineate a consistent picture of the initial stage of the surface oxidation: (i) at low coverage, in the limit of single O insertions, oxygen is incorporated on the surface forming a twofold-bridging Ga-O-As bond; (ii) at increasing coverage, as multiple O atoms are involved, this is accompanied by the formation of a threefold-coordinated bond (with two Ga and one As atoms); (iii) the latter has important implications regarding the electronic properties of the adsorbate since this O bonding may result in the formation of As dangling bonds. Moreover, a clear trend of increased energy gain for the incorporation of neighboring O atoms compared to single O insertions indicates that the formation of oxide clusters is favored over a regime of uniform oxidation. Our findings provide a detailed description of the O bonding and stress the importance of modeling the adsorption of multiple O atoms for an accurate description of the surface oxidation.

Abstract: Monte Carlo (MC) methods have a long-standing history as partners of molecular dynamics (MD) to simulate the evolution of materials at the atomic scale. Among these techniques, the uniform-acceptance force-bias Monte Carlo (UFMC) method [ G. Dereli Mol. Simul. 8 351 (1992)] has recently attracted attention [ M. Timonova et al. Phys. Rev. B 81 144107 (2010)] thanks to its apparent capacity of being able to simulate physical processes in a reduced number of iterations compared to classical MD methods. The origin of this efficiency remains, however, unclear. In this work we derive a UFMC method starting from basic thermodynamic principles, which leads to an intuitive and unambiguous formalism. The approach includes a statistically relevant time step per Monte Carlo iteration, showing a significant speed-up compared to MD simulations. This time-stamped force-bias Monte Carlo (tfMC) formalism is tested on both simple one-dimensional and three-dimensional systems. Both test-cases give excellent results in agreement with analytical solutions and literature reports. The inclusion of a time scale, the simplicity of the method, and the enhancement of the time step compared to classical MD methods make this method very appealing for studying the dynamics of many-particle systems.