Abstract: We perform the wavepacket dynamics simulation of a graphene-based device where propagating electron trajectories are tamed by an applied magnetic field toward a normal/superconductor interface. The magnetic field controls the incidence angle of the incoming electronic wavepacket at the interface, which results in the tunable electron-hole ratio in the reflected wave function due to the angular dependence of the Andreev reflection. Here, mapped control of the quasiparticle trajectories by the external magnetic field not only defines an experimental probe for fundamental studies of the Andreev reflection in graphene but also lays the foundation for further development of magnetic focusing devices based on nanoengineered superconducting two-dimensional materials.

Abstract: The effect of the magnetic impurities on the charge transport in a magnetic topological ultra-thin film (MTF) is analytically investigated by applying the semi-classical Boltzmann framework through a modified relaxation-time approximation. Our results for the relaxation time of electrons as well as the charge conductivity of the system exhibit two distinct regimes of transport. We show that the generated charge current in a MTF is always dissipative and anisotropic when both conduction bands are involved in the charge transport. The magnetic impurities induce a chirality selection rule for the transitions of electrons which can be altered by changing the orientation of the magnetic impurities. On the other hand, when a single conduction band participates in the charge transport, the resistivity is isotropic and can be entirely suppressed due to the corresponding chirality selection rule. Our findings propose a method to determine an onset thickness at which a crossover from a three-dimensional magnetic topological insulator to a (two-dimensional) MTF occurs.

Abstract: The physical properties of clays and micas can be controlled by exchanging ions in the crystal lattice. Atomically thin materials can have superior properties in a range of membrane applications, yet the ion-exchange process itself remains largely unexplored in few-layer crystals. Here we use atomic-resolution scanning transmission electron microscopy to study the dynamics of ion exchange and reveal individual ion binding sites in atomically thin and artificially restacked clays and micas. We find that the ion diffusion coefficient for the interlayer space of atomically thin samples is up to 10(4) times larger than in bulk crystals and approaches its value in free water. Samples where no bulk exchange is expected display fast exchange at restacked interfaces, where the exchanged ions arrange in islands with dimensions controlled by the moire superlattice dimensions. We attribute the fast ion diffusion to enhanced interlayer expandability resulting from weaker interlayer binding forces in both atomically thin and restacked materials. This work provides atomic scale insights into ion diffusion in highly confined spaces and suggests strategies to design exfoliated clay membranes with enhanced performance. Layered clays are of interest for membranes and many other applications but their ion-exchange dynamics remain unexplored in atomically thin materials. Here, using electron microscopy, it is found that the ion diffusion for few-layer two-dimensional clays approaches that of free water and that superlattice cation islands can form in twisted and restacked materials.

Abstract: Two-dimensional Transition Metal Dichalcogenides (TMDs) possessing extraordinary physical properties at reduced dimensionality have attracted interest due to their promise in electronic and optical device applications. However, TMD monolayers can show a broad range of different properties depending on their crystal phase; for example, H phases are usually semiconductors, while the T phases are metallic. Thus, controlling phase transitions has become critical for device applications. In this study, the energetically low-lying crystal structures of pristine and Janus TMDs are investigated by using ab initio Nudged Elastic Band and molecular dynamics simulations to provide a general explanation for their phase stability and transition properties. Across all materials investigated, the T phase is found to be the least stable and the H phase is the most stable except for WTe2, while the T' and T '' phases change places according to the TMD material. The transition energy barriers are found to be large enough to hint that even the higher energy phases are unlikely to undergo a phase transition to a more stable phase if they can be achieved except for the least stable T phase, which has zero barrier towards the T ' phase. Indeed, in molecular dynamics simulations the thermodynamically least stable T phase transformed into the T ' phase spontaneously while in general no other phase transition was observed up to 2100 K for the other three phases. Thus, the examined T ', T '' and H phases were shown to be mostly stable and do not readily transform into another phase. Furthermore, so-called mixed phase calculations considered in our study explain the experimentally observed lateral hybrid structures and point out that the coexistence of different phases is strongly stable against phase transitions. Indeed, stable complex structures such as metal-semiconductor-metal architectures, which have immense potential to be used in future device applications, are also possible based on our investigation.

Abstract: The DFT-1/2 method is a band-gap correction with GW precision at a density functional theory (DFT) computational cost. The method was also extended to correct the gap between defect levels, allowing for the calculation of optical transitions. However, this method fails when the atomic character of the occupied and unoccupied defect levels is similar as we illustrate by two examples, the tetrahedral hydrogen interstitial and the negatively charged vacancy in diamond. We solve this problem by decoupling the effect of the occupied and unoccupied defect levels and call this the decoupled DFT-1/2 method for defects.

Abstract: This study focuses on optimizing the composition and firing temperature of porcelain tiles using statistical analysis techniques. A full factorial design, including model adequacy checking, analysis of variance, Pareto charts, interaction plots, regression model, and response optimizer is employed. The key factors were the Seger ratios of SiO2/Al2O3, Na2O/K2O, MgO/CaO, and firing temperature. The response variables investigated were bulk density, water absorption, linear shrinkage, coefficient of thermal expansion (at 500 degrees C), and strength. The statistical analysis revealed highly significant results, which were further validated, confirming their reliability for practical use in the production of porcelain tiles. The study demonstrated the effectiveness of utilizing Seger formulas and properties of typical raw materials to accurately predict the final properties of ceramic tiles. By employing SiO2/Al2O3 = 5.2, Na2O/K2O = 1.50, MgO/CaO = 3.0, and firing temperature of 1180 degrees C, optimized properties, such as maximum strength, maximum bulk density, and minimum water absorption, was achieved with a composite desirability of 0.9821.

Abstract: Based on first-principles calculations, we study the structural stability of the H and T phases of monolayer MoS2 upon Li doping. Our calculations demonstrate that it is possible to stabilize a distorted T phase of MoS2 over the H phase through adsorption of Li atoms on the MoS2 surface. Through molecular dynamics and phonon calculations, we show that the T phase of MoS2 is dynamically unstable and undergoes considerable distortions. The type of distortion depends on the concentration of adsorbed Li atoms and changes from zigzag-like to diamond-like when increasing the Li doping. There exists a substantial energy barrier to transform the stable H phase to the distorted T phases, which is considerably reduced by increasing the concentration of Li atoms. We show that it is necessary that the Li atoms adsorb on both sides of the MoS2 monolayer to reduce the barrier sufficiently. Two processes are examined that allow for such two-sided adsorption, namely, penetration through the MoS2 layer and diffusion over the MoS2 surface. We show that while there is only a small barrier of 0.24 eV for surface diffusion, the amount of energy needed to pass through a pure MoS2 layer is of the order of similar or equal to 2 eV. However, when the MoS2 layer is covered with Li atoms the amount of energy that Li atoms should gain to penetrate the layer is drastically reduced and penetration becomes feasible.

Abstract: The structural, electronic, and magnetic properties of pristine, defective, and oxidized monolayer TiS3 are investigated using first-principles calculations in the framework of density functional theory. We found that a single layer of TiS3 is a direct band gap semiconductor, and the bonding nature of the crystal is fundamentally different from other transition metal chalcogenides. The negatively charged surfaces of single layer TiS3 makes this crystal a promising material for lubrication applications. The formation energies of possible vacancies, i.e. S, Ti, TiS, and double S, are investigated via total energy optimization calculations. We found that the formation of a single S vacancy was the most likely one among the considered vacancy types. While a single S vacancy results in a nonmagnetic, semiconducting character with an enhanced band gap, other vacancy types induce metallic behavior with spin polarization of 0.3-0.8 mu(B). The reactivity of pristine and defective TiS3 crystals against oxidation was investigated using conjugate gradient calculations where we considered the interaction with atomic O, O-2, and O-3. While O-2 has the lowest binding energy with 0.05-0.07 eV, O-3 forms strong bonds stable even at moderate temperatures. The strong interaction (3.9-4.0 eV) between atomic O and TiS3 results in dissociative adsorption of some O-containing molecules. In addition, the presence of S-vacancies enhances the reactivity of the surface with atomic O, whereas it had a negative effect on the reactivity with O-2 and O-3 molecules.

Abstract: A weakly coupled system of two crossed graphene nanoribbons exhibits direct tunneling due to the overlap of the wavefunctions of both ribbons. We apply the Bardeen transfer Hamiltonian formalism, using atomistic band structure calculations to account for the effect of the atomic structure on the tunneling process. The strong quantum-size confinement of the nanoribbons is mirrored by the one-dimensional character of the electronic structure, resulting in properties that differ significantly from the case of inter-layer tunneling, where tunneling occurs between bulk two-dimensional graphene sheets. The current-voltage characteristics of the inter-ribbon tunneling structures exhibit resonance, as well as stepwise increases in current. Both features are caused by the energetic alignment of one-dimensional peaks in the density-of-states of the ribbons. Resonant tunneling occurs if the sign of the curvature of the coupled energy bands is equal, whereas a step-like increase in the current occurs if the signs are opposite. Changing the doping modulates the onset-voltage of the effects as well as their magnitude. Doping through electrostatic gating makes these structures promising for application towards steep slope switching devices. Using the atomistic empirical pseudopotentials based Bardeen transfer Hamiltonian method, inter-ribbon tunneling can be studied for the whole range of two-dimensional materials, such as transition metal dichalcogenides. The effects of resonance and of step-like increases in the current we observe in graphene ribbons are also expected in ribbons made from these alternative two-dimensional materials, because these effects are manifestations of the one-dimensional character of the density-of-states. Published by AIP Publishing.

Abstract: Using advanced Ginzburg-Landau simulations, we study the superconducting state of a thin superconducting film under a ferromagnetic layer, separated by an insulating oxide, in applied external magnetic field and electric current. The taken uniaxial ferromagnet is organized into a series of parallel domains with alternating polarization of out-of-plane magnetization, sufficiently strong to induce vortex-antivortex pairs in the underlying superconductor in absence of other magnetic field. We show the organization of such vortex-antivortex matter into rich configurations, some of which are not matching the periodicity of the ferromagnetic film. The variety of possible configurations is enhanced by applied homogeneous magnetic field, where additional vortices in the superconductor may lower the energy of the system by either annihilating the present antivortices under negative ferromagnetic domains or by lowering their own energy after positioning under positive ferromagnetic domains. As a consequence, both the vortex-antivortex reordering in increasing external field and the evolution of the energy of the system are highly nontrivial. Finally, we reveal the very interesting effects of applied dc electric current on the vortex-antivortex configurations, since resulting Lorentzian force has opposite direction for vortices and antivortices, while direction of the applied current with respect to ferromagnetic domains is of crucial importance for the interaction of the applied and the Meissner current, as well as the consequent vortex-antivortex dynamics-both of which are reflected in the anisotropic critical current of the system. (C) 2016 AIP Publishing LLC.

Abstract: The feasibility of Veselago lensing in graphene with a p-n junction is investigated numerically for realistic injection leads. Two different set-ups with two narrow leads are considered with absorbing or reflecting side edges. This allows us to separately determine the influence of scattering on electron focusing for the edges and the p-n interface. Both semiclassical and tight-binding simulations show a distinctive peak in the transmission probability that is attributed to the Veselago lensing effect. We investigate the robustness of this peak on the width of the injector, the position of the p-n interface, and different gate potential profiles. Furthermore, the influence of scattering by both short- and long-range impurities is considered.

Abstract: Efficient quantum mechanical simulation of tunnel field-effect transistors (TFETs) is indispensable to allow for an optimal configuration identification. We therefore present a full-zone 15-band quantum mechanical solver based on the envelope function formalism and employing a spectral method to reduce computational complexity and handle spurious solutions. We demonstrate the versatility of the solver by simulating a 40 nm wide In0.53Ga0.47As lineTFET and comparing it to p-n-i-n configurations with various pocket and body thicknesses. We find that the lineTFET performance is not degraded compared to semi-classical simulations. Furthermore, we show that a suitably optimized p-n-i-n TFET can obtain similar performance to the lineTFET. (C) 2015 AIP Publishing LLC.

Abstract: In this study, we present a theoretical investigation of structural, electronic, and mechanical properties of pentagonal monolayers of carbon (p-graphene), boron nitride (p-B2N4 and p-B4N2), and silver azide (p-AgN3) by performing state-of-the-art first principles calculations. Our total energy calculations suggest feasible formation of monolayer crystal structures composed entirely of pentagons. In addition, electronic band dispersion calculations indicate that while p-graphene and p-AgN3 are semiconductors with indirect bandgaps, p-BN structures display metallic behavior. We also investigate the mechanical properties (in-plane stiffness and the Poisson's ratio) of four different pentagonal structures under uniaxial strain. p-graphene is found to have the highest stiffness value and the corresponding Poisson's ratio is found to be negative. Similarly, p-B2N4 and p-B4N2 have negative Poisson's ratio values. On the other hand, the p-AgN3 has a large and positive Poisson's ratio. In dynamical stability tests based on calculated phonon spectra of these pentagonal monolayers, we find that only p-graphene and p-B2N4 are stable, but p-AgN3 and p-B4N2 are vulnerable against vibrational excitations.

Abstract: Two-dimensional (2D) carbon nitride materials play an important role in energy-harvesting, energy-storage and environmental applications. Recently, a new carbon nitride, 2D polyaniline (C3N) was proposed [PNAS 113 (2016) 7414-7419]. Based on the structure model of this C3N monolayer, we propose two new carbon nitride monolayers, named dumbbell (DB) C4N-I and C4N-II. Using first-principles calculations, we systematically study the structure, stability, and band structure of these two materials. In contrast to other carbon nitride monolayers, the orbital hybridization of the C/N atoms in the DB C4N monolayers is sp(3). Remarkably, the band structures of the two DB C4N monolayers have a Dirac cone at the K point and their Fermi velocities (2.6/2.4 x 10(5) m/s) are comparable to that of graphene. This makes them promising materials for applications in high-speed electronic devices. Using a tight-binding model, we explain the origin of the Dirac cone. (C) 2017 Elsevier Ltd. All rights reserved.

Abstract: A first-principles study of the structural and electronic properties of carbon impurities in CuInSe2 and CuGaSe2 is presented. Carbon is present in organic molecules in the precursor solutions used in non-vacuum growth methods for CuInSe2 and CuGaSe2 based photovoltaic cells. These growth methods make more efficient use of material, time, and energy than traditional vacuum methods. The formation energies of several carbon impurities are calculated using the hybrid HSE06 functional. C Cu acts as a shallow donor, CIn and interstitial C yield deep donor levels in CuInSe2, while in CuGaSe2 CGa and interstitial C act as deep amphoteric defects. So, these defects reduce the majority carrier (hole) concentration in p-type CuInSe2 and CuGaSe2 by compensating the acceptor levels. The deep defects are likely to act as recombination centers for the photogenerated charge carriers and are thus detrimental for the performance of the photovoltaic cells. On the other hand, the formation energies of the carbon impurities are high, even under C-rich growth conditions. Thus, few C impurities will form in CuInSe2 and CuGaSe2 in thermodynamic equilibrium. However, the deposition of the precursor solution in non-vacuum growth methods presents conditions far from thermodynamic equilibrium. In this case, our calculations show that C impurities formed in non-equilibrium tend to segregate from CuInSe2 and CuGaSe2 by approaching thermodynamic equilibrium, e.g., via thorough annealing.