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: The lowest-energy configurations of Cn(n≤55) clusters are obtained using the energy-minimization technique with the conjugate gradient method where a modified Brenner potential is invoked to describe the carbon and hydrocarbon interaction. We found that the ground-state configuration consists of a single ring for small number of C atoms and multiring structures are found with increasing n, which can be in planar, bowl-like or caplike form. Contrary to previous predictions, the binding energy Eb does not show even-odd oscillations and only small jumps are found in the Eb(n) curve as a consequence of specific types of edges or equivalently the number of secondary atoms. We found that hydrogenation of the edge atoms may change the ground-state configuration of the nanocluster. In both cases we determined the magic clusters. Special attention is paid to trigonal and hexagonal shaped carbon clusters and to clusters having a graphenelike configuration. Trigonal clusters are never the ground state while hexagonal-shaped clusters are only the ground state when they have zigzag edges.

Abstract: We propose that low-strain Si/O superlattices can be constructed by connecting reconstructed Si{001} surfaces by Si-O-Si bridges. Ab initio calculations show that our models are energetically more favorable than all the models proposed so far. The part of our Si/O superlattice model is experimentally accessible just by oxidizing a Si( 001) substrate. To complete our Si/O superlattice model, we propose a three-step method. We also explore the potential of our Si/O superlattice models for new materials used in future Si electronics. We find that the location of the channel where the carriers travel can be controlled between the interfaces and the Si layers by the insertion of O atoms into the Si-Si dimers. By revealing the origins of the interface electron and hole states, we find that similar interface states should be easily achieved for Si slabs and Si substrates. Interestingly, the interface electrons and holes have small effective masses in the direction parallel to the channel and large effective masses in the direction normal to the channel, which makes the Si/O superlattices attractive to be used for channel materials. We also find that the valley splitting of Si is enhanced by the formation of the Si/O/Si interfaces, which is ideal for developing Si-based qubits. Our findings open new perspectives to design and control the electronic properties of Si.

Abstract: Both even- and odd-numbered neutral carbon clusters Cn (n=210) are systematically studied using the energy minimization method and the modified Brenner potential for the carbon-carbon interactions. Many stable configurations were found, and several new isomers are predicted. For the lowest energy stable configurations we obtained their binding energies and bond lengths. We found that for n5 the linear isomer is the most stable one while for n>5 the monocyclic isomer becomes the most stable. The latter was found to be regular for all studied clusters. The dependence of the binding energy for linear and cyclic clusters versus the cluster size n (n=210) is found to be in good agreement with several previous calculations, in particular with ab initio calculations as well as with experimental data for n=25.

Abstract: To investigate the impact of both the CoFe ferromagnetic layer thickness and the capping paramagnetic layer on the tunnel magnetoresistance (TMR), we performed first-principles simulations on epitaxial magnetic tunnel junctions contacted with either CoFe or Ta paramagnetic capping layers. We observed a strong oscillation of the TMR amplitude with respect to the thickness of the ferromagnetic layer. The TMR is found to be amplified whenever the MgO spin tunnel barrier is thickened. Quantization of the electronic structure of the ferromagnetic layers is found to be at the origin of this oscillatory behavior. Metals such as Ta contacting the magnetic layer are found to enhance the amplitude of the oscillations due to the occurrence of an interface dipole. The latter drives the band alignment and tunes the nature of the spin channels that are active during the tunneling process. Subsequently, the regular transmission spin channels are modulated in the magnetic tunnel junction stack and other complex ones are being activated.

Abstract: Amorphous semiconductors are usually characterized by a low charge carrier mobility, essentially related to their lack of long-range order. The development of such material with higher charge carrier mobility is hence challenging. Part of the issue comes from the difficulty encountered by first-principles simulations to evaluate concepts such as the electron effective mass for disordered systems since the absence of periodicity induced by the disorder precludes the use of common concepts derived from condensed matter physics. In this paper, we propose a methodology based on first-principles simulations that partially solves this problem, by quantifying the degree of delocalization of a wave function and of the connectivity between the atomic sites within this electronic state. We validate the robustness of the proposed formalism on crystalline and molecular systems and extend the insights gained to disordered/amorphous InGaZnO4 and Si. We also explore the properties of p-type oxide semiconductor candidates recently reported to have a low effective mass in their crystalline phases [G. Hautier et al., Nat. Commun. 4, 2292 (2013)]. Although in their amorphous phase none of the candidates present a valence band with delocalization properties matching those found in the conduction band of amorphous InGaZnO4, three of the seven analyzed materials show some potential. The most promising candidate, K2Sn2O3, is expected to possess in its amorphous phase a slightly higher hole mobility than the electron mobility in amorphous silicon.

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: A two-dimensional self-consistent fluid model combined with the full set of Maxwell equations is developed to investigate an argon capacitively coupled plasma, focusing on the electromagnetic effects on the discharge characteristics at various discharge conditions. The results indicate that there exist distinct differences in plasma characteristics calculated with the so-called electrostatic model (i.e., without taking into account the electromagnetic effects) and the electromagnetic model (which includes the electromagnetic effects), especially at very high frequencies. Indeed, when the excitation source is in the high frequency regime and the electromagnetic effects are taken into account, the plasma density increases significantly and meanwhile the ionization rate evolves to a very different distribution when the electromagnetic effects are dominant. Furthermore, the dependence of the plasma characteristics on the voltage and pressure is also investigated, at constant frequency. It is observed that when the voltage is low, the difference between these two models becomes more obvious than at higher voltages. As the pressure increases, the plasma density profiles obtained from the electromagnetic model smoothly shift from edge-peaked over uniform to a broad maximum in the center. In addition, the edge effect becomes less pronounced with increasing frequency and pressure, and the skin effect rather than the standing-wave effect becomes dominant when the voltage is high.

Abstract: The physics of plasmas containing positive and negative ions is discussed with special attention to the recently produced pair-ion plasma containing ions of equal mass and opposite charge. The effects of the density gradient in the direction perpendicular to the ambient magnetic field vector are discussed. The possible presence of electrons is discussed in the context of plasma modes propagating at an angle with respect to the magnetic field vector. It is shown that the electron plasma mode may become a backward mode in the presence of a density gradient, and this behavior may be controlled either by the electron number density or the mode number in the perpendicular direction. In plasmas with hot electrons an instability may develop, driven by the combination of electron collisions and the density gradient, and in the regime of a sound ions' response. In the case of a pure pair-ion plasma, for lower frequencies and for parameters close to those used in the recent experiments, the perturbed ions may feel the effects of the magnetic field. In this case the plasma mode also becomes backward, resembling features of an experimentally observed but yet unexplained backward mode. (C) 2008 American Institute of Physics.

Abstract: Particle-in-cell/Monte Carlo collision (PIC/MCC) simulations are performed to investigate the asymmetric secondary electron emission (SEE) effects when electrons strike two different material electrodes in low pressure capacitively coupled plasmas (CCPs). To describe the electron-surface interactions, a realistic model, considering the primary electron impact energy and angle, as well as the corresponding surface property-dependent secondary electron yields, is employed in PIC/MCC simulations. In this model, three kinds of electrons emitted from the surface are considered: (i) elastically reflected electrons, (ii) inelastically backscattered electrons, and (iii) electron induced secondary electrons (SEs, i.e., delta-electrons). Here, we examined the effects of electron-surface interactions on the ionization dynamics and plasma characteristics of an argon discharge. The discharge is driven by a voltage source of 13.56MHz with amplitudes in the range of 200-2000V. The grounded electrode material is copper (Cu) for all cases, while the powered electrode material is either Cu or silicon dioxide (SiO2). The simulations reveal that the electron impact-induced SEE is an essential process at low pressures, especially at high voltages. Different electrode materials result in an asymmetric response of SEE. Depending on the instantaneous local sheath potential and the phase of the SEE, these SEs either are reflected by the opposite sheath or strike the electrode surface, where they can induce delta-electrons upon their residual energies. It is shown that highly energetic delta-electrons contribute significantly to the ionization rate and a self-bias forms when the powered electrode material is assumed to be made of SiO2. Complex dynamics is observed due to the multiple electron-surface interaction processes and asymmetric yields of SEs in CCPs.

Abstract: Superimposition of dual frequencies (DFs) is one of the methods used for controlling plasma distribution in an inductively coupled plasma (ICP) source. The effects of a superimposed DF on the argon plasma characteristics have been investigated using a two-dimensional self-consistent fluid model. When both currents are fixed at 6A, the plasma density drops with decrease in one of the source frequencies due to less efficient heating and the plasma uniformity improves significantly. Moreover, for ICP operated with superimposed DFs (i.e., 4.52MHz/13.56MHz and 2.26MHz/13.56MHz), the current source exhibits the same period as the low frequency (LF) component, and the plasma density is higher than that obtained at a single frequency (i.e., 4.52 and 2.26MHz) with the same total current of 12A. However, at superimposed current frequencies of 6.78MHz/13.56MHz, the plasma density is lower than that obtained at a single frequency of 6.78MHz due to the weaker negative azimuthal electric field between two positive maxima during one period of 6.78MHz. When the superimposed DF ICP operates at 2.26 and 13.56MHz, the rapid oscillations of the induced electric field become weaker during one period of 2.26MHz as the current ratio of 2.26MHz/13.56MHz rises from 24A/7 A to 30A/1 A, and the plasma density drops with the current ratio due to weakened electron heating. The uniformity of plasma increases due to sufficient diffusion under the low-density condition.

Abstract: Gold nanocluster properties exhibit unique size-dependence. In this contribution, we employ reactive molecular dynamics simulations to calculate the size- and temperature-dependent surface energies, strain energies and atomic displacements for icosahedral, cuboctahedral, truncated octahedral and decahedral Au-nanoclusters. The calculations demonstrate that the surface energy decreases with increasing cluster size at 0 K but increases with size at higher temperatures. The calculated melting curves as a function of cluster size demonstrate the Gibbs-Thomson effect. Atomic displacements and strain are found to strongly depend on the cluster size and both are found to increase with increasing cluster size. These results are of importance for understanding the size-and temperature-dependent surface processes on gold nanoclusters.