Abstract: We use a variational technique to study the many-body phase diagram of electrons and holes in n-doped and p-doped monolayer transition metal dichalcogenides (TMDs). We find a total of four different phases. (i) A fully spin polarized and valley polarized ferromagnetic state. (ii) A state with no global spin polarization but with spin polarization in each valley separately, i.e., spin-valley locking. (iii) A state with spin polarization in one of the valleys and little to no spin polarization in the other valley. (iv) A paramagnetic state with no valley polarization. These phases are separated by first-order phase transitions and are determined by the particle density and the dielectric constant of the substrate. We find that in the presence of a perpendicular magnetic field the four different phases persist. In the case of n-doped MoS2, a fifth phase, which is completely valley polarized but not spin polarized, appears for magnetic fields larger than 7 T and for magnetic fields larger than 23 T completely replaces the second phase.

Abstract: Bilayer hexagonal borophene, which is bound together through pillars, is a novel topological semimetal. Using density functional theory, we investigate its electronic band structure and show that it is a Dirac material which exhibits a nodal line. A tight-binding model was constructed based on the Slater-Koster approach, which accurately models the electronic spectrum. We constructed an effective four-band model Hamiltonian to describe the spectrum near the nodal line. This Hamiltonian can be used as a new platform to study the new properties of nodal line semimetals. We found that the nodal line is created by edge states and is very robust against perturbations and impurities. Breaking symmetries can split the nodal line, but cannot open a gap.

Abstract: Starting from the single-particle Dirac Hamiltonian for charge carriers in monolayer transition metal dichalcogenides (TMDs), we construct a four-band Hamiltonian describing interlayer excitons consisting of an electron in one TMD layer and a hole in the other TMD layer. An expression for the electron-hole interaction potential is derived, taking into account the effect of the dielectric environment above, below, and between the two TMD layers as well as polarization effects in the transition metal layer and in the chalcogen layers of the TMD layers. We calculate the interlayer exciton binding energy and average in-plane interparticle distance for different TMD heterostructures. The effect of different dielectric environments on the exciton binding energy is investigated and a remarkable dependence on the dielectric constant of the barrier between the two layers is found, resulting from competing effects as a function of the in-plane and out-of-plane dielectric constants of the barrier. The polarization effects in the chalcogen layers, which in general reduce the exciton binding energy, can lead to an increase in binding energy in the presence of strong substrate effects by screening the substrate. The excitonic absorbance spectrum is calculated and we show that the interlayer exciton peak depends linearly on a perpendicular electric field, which agrees with recent experimental results.

Abstract: Plasma catalysis is gaining increasing interest for various environmental applications. Catalytic
material can be inserted in different shapes in the plasma, e.g., as pellets, (coated) beads, but also
as honeycomb monolith and 3DFD structures, also called ‘structured catalysts’, which have high
mass and heat transfer properties. In this work, we examine the streamer discharge propagation
and the interaction between plasma and catalysts, inside the channels of such structured catalysts,
by means of a two-dimensional particle-in-cell/Monte Carlo collision model. Our results reveal
that plasma streamers behave differently in various structured catalysts. In case of a honeycomb
structure, the streamers are limited to only one channel, with low or high plasma density when
the channels are parallel or perpendicular to the electrodes, respectively. In contrast, in case of a
3DFD structure, the streamers can distribute to different channels, causing discharge
enhancement due to surface charging on the dielectric walls of the structured catalyst, and
especially giving rise to a broader plasma distribution. The latter should be beneficial for plasma
catalysis applications, as it allows a larger catalyst surface area to be exposed to the plasma.

Abstract: Using the tight-binding approach, we investigate the energy spectrum of square, triangular, and hexagonal MoS2 quantum dots (QDs) in the presence of a perpendicular magnetic field. Novel edge states emerge in MoS2 QDs, which are distributed over the whole edge which we call ring states. The ring states are robust in the presence of spin-orbit coupling (SOC). The corresponding energy levels of the ring states oscillate as a function of the perpendicular magnetic field which are related to Aharonov-Bohm oscillations. Oscillations in the magnetic field dependence of the energy levels and the peaks in the magneto-optical spectrum emerge (disappear) as the ring states are formed (collapsed). The period and the amplitude of the oscillation decrease with the size of the MoS2 QDs.

Abstract: Methane needle-plane discharge has practical application prospect and scientific research significance since methane conversion heavy oil hydrogenation is formed by coupling methane needle-plane discharge with heavy oil hydrogenation, which can achieve high-efficient heavy oil hydrogenation and increase the yields of high value-added light olefins. In this paper, a two-dimensional fluid model is built up for numerically simulating the methane needle-plane discharge plasma at atmospheric pressure. Spatial and axial distributions of electric intensity, electron temperature and particle densities are obtained. Reaction yields are summarized and crucial pathways to produce various kinds of charged and neutral particles are found out. Simulation results indicate that axial evolutions of CH3+ and CH4+ densities, electric intensity and electron temperature are similar and closely related. The CH5+ and C2H5+ densities first increase and then decrease along the axial direction. The CH3 and H densities have nearly identical spatial and axial distributions. Particle density distributions of CH2, C2H4 and C2H5 are obviously different in the area near the cathode but comparatively resemblant in the positive column region. The CH3+ and CH4+ are produced by electron impact ionizations between electrons and CH4. The CH5+ and C2H5+ are respectively generated by molecular impact dissociations between CH3+ and CH4 and between CH4+ and CH4. Electron impact decomposition between electrons and CH4 is a dominated reaction to produce CH3, CH2, CH and H. The reactions between CH2 and CH4 and between electrons and C2H4 are critical pathways to produce C2H4 and C2H2, respectively. In addition, the yields of electron impact decomposition reactions between electrons and CH4 and reactions between CH2 and CH4 account for 52.15% and 47.85% of total yields of H-2 respectively.

Abstract: The dependence of the excitonic photoluminescence (PL) spectrum of monolayer transition metal dichalcogenides (TMDs) on the tilt angle of an applied magnetic field is studied. Starting from a four-band Hamiltonian we construct a theory which quantitatively reproduces the available experimental PL spectra for perpendicular and in-plane magnetic fields. In the presence of a tilted magnetic field, we demonstrate that the dark exciton PL peaks brighten due to the in-plane component of the magnetic field and split for light with different circular polarizations as a consequence of the perpendicular component of the magnetic field. This splitting is more than twice as large as the splitting of the bright exciton peaks in tungsten-based TMDs. We propose an experimental setup that will allow for accessing the predicted splitting of the dark exciton peaks in the PL spectrum.

Abstract: Defects are almost inevitably present in realistic materials and defective materials are expected to exhibit very different properties than their nondefective (perfect) counterparts. Here, using a combination of the tight-binding approach and the scattering matrix formalism, we investigate the electronic transport properties of defective phosphorene nanoribbons (PNRs) containing atomic vacancies. We find that for both armchair PNRs (APNRs) and zigzag PNRs (ZPNRs), single vacancies can create quasilocalized states, which can affect their conductance. With increasing vacancy concentration, three different transport regimes are identified: ballistic, diffusive, and Anderson localized ones. In particular, ZPNRs that are known to be metallic due to the presence of edge states become semiconducting: edge conductance vanishes and transport gap appears due to Anderson localization. Moreover, we find that for a fixed vacancy concentration, both APNRs and ZPNRs of narrower width and/or longer length are more sensitive to vacancy disorder than their wider and/or shorter counterparts, and that for the same ribbon length and width, ZPNRs are more sensitive to vacancy disorder than APNRs.

Abstract: The permeation of water between neighboring graphene oxide (GO) flakes, i.e., 2D nanochannels, are investigated using a simple model for the GO membrane. We simulate the hydrophilic behavior of nanocapillaries and study the effect of surface charge on the dynamical properties of water flow and the influence of Na+ and Cl- ions on water permeation. Our approach is based on extensive equilibrium molecular dynamics simulations to obtain a better understanding of water permeation through charged nanochannels in the presence of ions. We found significant change in the slippage dynamics of confined water such as a profound increase in viscosity/slip length with increasing charges over the surface. The slip length decreases one order of magnitude (i.e., 1/30) with increasing density of surface charge, while it increases by a factor of 2 with ion concentration. We found that commensurability induced by nanoconfinement plays an important role on the intrinsic dynamical properties of water.

Abstract: By means of replica exchange molecular dynamics simulations we investigate how the length of a silk-like, alternating diblock oligopeptide influences its secondary and quaternary structure. We carry out simulations for two protein sizes consisting of three and five blocks, and study the stability of a single protein, a dimer, a trimer and a tetramer. Initial configurations of our simulations are β-roll and β-sheet structures. We find that for the triblock the secondary and quaternary structures upto and including the tetramer are unstable: the proteins melt into random coil structures and the aggregates disassemble either completely or partially. We attribute this to the competition between conformational entropy of the proteins and the formation of hydrogen bonds and hydrophobic interactions between proteins. This is confirmed by our simulations on the pentablock proteins, where we find that, as the number of monomers in the aggregate increases, individual monomers form more hydrogen bonds whereas their solvent accessible surface area decreases. For the pentablock β-sheet protein, the monomer and the dimer melt as well, although for the β-roll protein only the monomer melts. For both trimers and tetramers remain stable. Apparently, for these the entropy loss of forming β-rolls and β-sheets is compensated for in the free-energy gain due to the hydrogen-bonding and hydrophobic interactions. We also find that the middle monomers in the trimers and tetramers are conformationally much more stable than the ones on the top and the bottom. Interestingly, the latter are more stable on the tetramer than on the trimer, suggesting that as the number of monomers increases protein-protein interactions cooperatively stabilize the assembly.
According to our simulations, the β-roll and β-sheet aggregates must be approximately equally
stable.

Abstract: In the quest for miniaturizing magnetic devices, the effects of boundaries and surfaces become increasingly important. Here we show how the recently predicted boundary-induced Dzyaloshinskii-Moriya interaction (DMI) affects the magnetization of ferromagnetic films with a C-infinity v symmetry and a perpendicular magnetic anisotropy. For an otherwise uniformly magnetized film, we find a surface twist when the magnetization in the bulk is canted by an in-plane external field. This twist at the surfaces caused by the boundary-induced DMI differs from the common canting caused by internal DMI observed at the edges of a chiral magnet. Furthermore, we find that the surface twist due to the boundary-induced DMI strongly affects the width of the domain wall at the surfaces. We also find that the skyrmion radius increases in the depth of the film, with the average size of the skyrmion increasing with boundary-induced DMI. This increase suggests that the boundary-induced DMI contributes to the stability of the skyrmion.

Abstract: Previous studies have shown that positron-annihilation spectroscopy is a highly sensitive probe of the electronic structure and surface composition of ligand-capped semiconductor quantum dots (QDs) embedded in thin films. The nature of the associated positron state, however, whether the positron is confined inside the QDs or localized at their surfaces, has so far remained unresolved. Our positron-annihilation lifetime spectroscopy studies of CdSe QDs reveal the presence of a strong lifetime component in the narrow range of 358–371 ps, indicating abundant trapping and annihilation of positrons at the surfaces of the QDs. Furthermore, our ab initio calculations of the positron wave function and lifetime employing a recent formulation of the weighted density approximation demonstrate the presence of a positron surface state and predict positron lifetimes close to experimental values. Our study thus resolves the long-standing question regarding the nature of the positron state in semiconductor QDs and opens the way to extract quantitative information on surface composition and ligand-surface interactions of colloidal semiconductor QDs through highly sensitive positron-annihilation techniques.

Abstract: We have discovered two novel types of planar defects that appear in heteroepitaxial YBa2Cu3O7-delta(YBCO123) thin films, grown by pulsed-laser deposition (PLD) either with or without a La2/3Ca1/3MnO3 (LCMO) overlayer, using the combination of highangle annular dark-field scanning transmission electron microscopy (HAADF-STEM) imaging and electron energy loss spectroscopy (EELS) mapping for unambiguous identification. These planar lattice defects are based on the intergrowth of either a BaO plane between two CuO chains or multiple Y-O layers between two CuO2 planes, resulting in non-stoichiometric layer sequences that could directly impact the high-Tc superconductivity. (C) 2018 Author(s).

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: 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.