Abstract: A theoretical study has been conducted of ultrashort pulsed laser induced electron emission from an aluminum surface. Electron emission fluxes retrieved from the commonly employed Fowler-DuBridge theory were compared to fluxes based on a laser-induced non-equilibrium electron distribution. As a result, the two-and three-photon photoelectron emission parameters for the Fowler-DuBridge theory have been approximated. We observe that at regimes where photoemission is important, laser-induced electron emission evolves in a more smooth manner than predicted by the Fowler-DuBridge theory. The importance of the actual electron distribution decreases at higher laser fluences, whereas the contribution of thermionic emission increases. Furthermore, the influence of a space charge effect on electron emission was evaluated by a one dimensional particle-in-cell model. Depending on the fluences, the space charge reduces the electron emission by several orders of magnitude. The influence of the electron emission flux profiles on the effective electron emission was found to be negligible. However, a non-equilibrium electron velocity distribution increases the effective electron emission significantly. Our results show that it is essential to consider the non-equilibrium electron distribution as well as the space charge effect for the description of laser-induced photoemission. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4729071]

Abstract: The role of the Cu atoms sputtered from the cathode material in a cylindrical hollow cathode discharge (HCD) and the corresponding Cu+ ions are studied with a self-consistent model based on the principle of Monte Carlo (MC) and fluid simulations. In order to obtain a more realistic view of the discharge processes, this model is coupled with other submodels, which describe the behavior of electrons, fast Ar atoms, Ar+ ions, and Ar metastable atoms, also based on the principles of MC and fluid simulations. Typical results are, among others, the thermalization profile of the Cu atoms, the fast Cu atom, the thermal Cu atom and Cu+ ion fluxes and densities, and the energy distribution of the Cu+ ions. It was found that the contribution of the Ar+ ions to the sputtering was the most significant, followed by the fast Ar atoms. At the cathode bottom, there was no net sputtered flux but a net amount of redeposition. Throughout the discharge volume, at all the conditions investigated, the largest concentration of Cu atoms was found in the lower half of the HCD, close to the bottom. Penning ionization was found the main ionization mechanism for the Cu atoms. The ionization degree of copper atoms was found to be in the same order as for the argon atoms (10-4). (c) 2005 American Institute of Physics.

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

Abstract: We report on the thin film resistivity of several platinum-group metals (Ru, Pd, Ir, and Pt). Platinum-group thin films show comparable or lower resistivities than Cu for film thicknesses below about 5 nm due to a weaker thickness dependence of the resistivity. Based on experimentally determined mean linear distances between grain boundaries as well as ab initio calculations of the electron mean free path, the data for Ru, Ir, and Cu were modeled within the semiclassical Mayadas-Shatzkes model [Phys. Rev. B 1, 1382 (1970)] to assess the combined contributions of surface and grain boundary scattering to the resistivity. For Ru, the modeling results indicated that surface scattering was strongly dependent on the surrounding material with nearly specular scattering at interfaces with SiO2 or air but with diffuse scattering at interfaces with TaN. The dependence of the thin film resistivity on the mean free path is also discussed within the Mayadas-Shatzkes model in consideration of the experimental findings. Published by AIP Publishing.

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: We examined the vacancy trapping proficiency of Sn and P atoms in germanium using positron annihilation spectroscopy measurements, sensitive to the open-volume defects. Epitaxial Ge1 xSnx films were grown by chemical vapor deposition with different P concentrations in the 3: 0 1019-1: 5 1020 cm 3 range. We corroborate our findings with first principles simulations. Codoping of Ge with a Sn concentration of up to 9% is not an efficient method to suppress the free vacancy concentration and the formation of larger phosphorus-vacancy complexes. Experimental results confirm an increase in the number of P atoms around the monovacancy with P-doping, leading to dopant deactivation in epitaxial germanium-tin layers with similar Sn content. Vice versa, no impact on the improvement of maximum achieved P activation in Ge with increasing Sn-doping has been observed. Theoretical calculations also confirm that Pn-V (vacancy) complexes are energetically more stable than the corresponding SnmPn-V and Snm-V defect structures with the same number of alien atoms (Sn or P) around the monovacancy. he strong attraction of vacancies to the phosphorus atoms remains the dominant dopant deactivation mechanism in Ge as well as in Ge1 xSnx. Published under license by AIP Publishing.

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: The unique physical properties of carbyne, a novel carbon nanostructure, have attracted considerable interest in modern nanotechnology. While carbyne synthesis has been accomplished successfully using diverse techniques, the underlying mechanisms governing the carbon monoxide-dependent catalytic synthesis of endohedral carbyne remain poorly understood. In this simulation-based study, we investigate the synthesis of endohedral carbyne from carbon and carbon monoxide radicals in the presence of a nickel catalyst inside double-walled carbon nanotubes with a (5,5)@(10,10) structure. The outcome of our investigation demonstrates that the incorporation of the carbon atom within the Ni-n@(5,5)@(10,10) model system initiates the formation of an elongated carbon chain. In contrast, upon the introduction of carbon monoxide radicals, the growth of the carbyne chain is inhibited as a result of the oxidation of endohedral nickel clusters by oxygen atoms after the initial steps of nucleation. Our findings align with prior theoretical, simulation, and experimental investigations, reinforcing their consistency and providing valuable insights into the synthesis of carbyne-based nanodevices that hold promising potential for future advancements in nanotechnology.

Abstract: Lipid oxidation is associated with several inflammatory and neurodegenerative diseases, but many questions to unravel its effects on biomembranes are still open due to the complexity of the topic. For instance, recent studies indicated that phase-separated domains can have a significant effect on membrane function. It is reported that domain interfaces are “hot spots” for pore formation, but the underlying mechanisms and the effect of oxidation-induced phase separation on membranes remain elusive. Thus, to evaluate the permeability of the membrane coexisting of liquid-ordered (Lo) and liquid-disordered (Ld) domains, we performed atomistic molecular dynamics simulations. Specifically, we studied the membrane permeability of nonoxidized or oxidized homogeneous membranes (single-phase) and at the Lo/Ld domain interfaces of heterogeneous membranes, where the Ld domain is composed of either oxidized or nonoxidized lipids. Our simulation results reveal that the addition of only 1.5% of lipid aldehyde molecules at the Lo/Ld domain interfaces of heterogeneous membranes increases the membrane permeability, whereas their addition at homogeneous membranes does not have any effect. This study is of interest for a better understanding of cancer treatment methods based on oxidative stress (causing among others lipid oxidation), such as plasma medicine and photodynamic therapy.

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

Abstract: The evaluation of nucleation rates from molecular dynamics trajectories is hampered by the slow nucleation time scale and impact of finite size effects. Here, we show that accurate nucleation rates can be obtained in a very general fashion relying only on the free energy barrier, transition state theory, and a simple dynamical correction for diffusive recrossing. In this setup, the time scale problem is overcome by using enhanced sampling methods, in casu metadynamics, whereas the impact of finite size effects can be naturally circumvented by reconstructing the free energy surface from an appropriate ensemble. Approximations from classical nucleation theory are avoided. We demonstrate the accuracy of the approach by calculating macroscopic rates of droplet nucleation from argon vapor, spanning 16 orders of magnitude and in excellent agreement with literature results, all from simulations of very small (512 atom) systems.