Abstract: Plasma is gaining interest for CH4 conversion into higher hydrocarbons and H2. However, the performance in terms of conversion and selectivity toward different hydrocarbons is different for different plasma types, and the underlying mechanisms are not yet fully understood. Therefore, we study here these mechanisms in different plasma sources, by means of a chemical kinetics model. The model is first validated by comparing the calculated conversions and hydrocarbon/H2 selectivities with experimental results in these different plasma types and over a wide range of specific energy input (SEI) values. Our model predicts that vibrational−translational nonequilibrium is negligible in all CH4 plasmas investigated, and instead, thermal conversion is important. Higher gas temperatures also lead to a more selective production of unsaturated hydrocarbons (mainly C2H2) due to neutral dissociation of CH4 and subsequent dehydrogenation processes, while three-body recombination reactions into saturated hydrocarbons (mainly C2H6, but also higher hydrocarbons) are dominant in low temperature plasmas.

Abstract: The surface termination of MXenes greatly determines the electrochemical properties and ion kinetics on their surfaces. So far, hydroxyl-, oxygen-, and fluorine-terminated MXenes have been widely studied for energy storage applications. Recently, sulfur-functionalized MXene structures, which possess low diffusion barriers, have been proposed as candidate materials to enhance battery performance. We performed first-principles calculations on the structural, stability, electrochemical, and ion dynamic properties of Li-adsorbed sulfur-functionalized groups 3B, 4B, 5B, and 6B transition-metal (M)-based MXenes (i.e., M2CS2 with M = Sc, Ti, Zr, Hf, V, Nb, Ta, Cr, Mo, and W). We performed phonon calculations, which indicated that all of the above M2CS2 MXenes, except for Sc, are dynamically stable at T = 0 K. The ground-state structure of each M2CS2 monolayer depends on the type of M atom. For instance, while sulfur prefers to sit at the FCC site on Ti2CS2, it occupies the HCP site of Cr-based MXene. We determined the Li adsorption configurations at different concentrations using the cluster expansion method. The highest maximum open-circuit voltages were computed for the group 4B element (i.e., Ti, Zr, and Hf)-based M2CS2, which are larger than 2.1 V, while their average voltages are approximately 1 V. The maximum voltage for the group 6B element (i.e., Cr, Mo, W)-based M2CS2 is less than 1 V, and the average voltage is less than 0.71 V. We found that S functionalization is helpful for capacity improvements over the O-terminated MXenes. In this respect, the computed storage gravimetric capacity may reach up to 417.4 mAh/g for Ti2CS2 and 404.5 mAh/g for V2CS2. Ta-, Cr-, Mo-, and W-based M2CS2 MXenes show very low capacities, which are less than 100 mAh/g. The Li surface diffusion energy barriers for all of the considered MXenes are less than 0.22 eV, which is favorable for high charging and discharging rates. Finally, ab initio molecular dynamic simulations performed at 400 K and bond-length analysis with respect to Li concentration verify that selected promising systems are robust against thermally induced perturbations that may induce structural transformations or distortions and undesirable Li release.

Abstract: The quantum electrostatic heterostructure (QEH) model allows for efficient computation of the dielectric screening properties of layered van der Waals (vdW)-bonded heterostructures in terms of the dielectric functions of the individual two-dimensional (2D) layers. Here, we extend the QEH model by including (1) contributions to the dielectric function from infrared active phonons in the 2D layers, (2) screening from homogeneous bulk substrates, and (3) intraband screening from free carriers in doped 2D semiconductor layers. We demonstrate the potential of the extended QEH model by calculating the dispersion of coupled phonons in multilayer stacks of hexagonal boron-nitride (hBN), the strong hybridization of plasmons and optical phonons in graphene/hBN heterostructures, the effect of substrate screening on the exciton series of monolayer MoS2, and the properties of hyperbolic plasmons in a doped phosphorene sheet. The new QEH code is distributed as a Python package with a simple command line interface and a comprehensive library of dielectric building blocks for the most common 2D materials, providing an efficient open platform for dielectric modeling of realistic vdW heterostructures.

Abstract: Within a tight-binding approach, we investigate the effect of twisting angle on the energy levels of circular bilayer graphene (BLG) quantum dots (QDs) in both the absence and presence of a perpendicular magnetic field. The QDs are defined by an infinite-mass potential, so that the specific edge effects are not present. In the absence of magnetic field (or when the magnetic length is larger than the moire length), we show that the low-energy states in twisted BLG QDs are completely affected by the formation of moire patterns, with a strong localization at AA-stacked regions. When magnetic field increases, the energy gap of an untwisted BLG QD closes with the edge states, localized at the boundaries between the AA- and AB-stacked spots in a twisted BLG QD. Our observation of the spatial localization of the electrons in twisted BLG QDs can be experimentally probed by low-bias scanning tunneling microscopy measurements.

Abstract: The hunt for high-temperature superfluidity has received new impetus from the discovery of atomically thin stable materials. Electron-hole superfluidity in coupled MoSe2-WSe2 monolayers is investigated using a mean-field multiband model that includes band splitting caused by strong spin-orbit coupling. This splitting leads to a large energy misalignment of the electron and hole bands which is strongly modified by interchanging the doping of the monolayers. The choice of doping determines if the superfluidity is tunable from one to two components.

Abstract: Two-dimensional graphene-based van der Waals heterostructures have received considerable interest because of their intriguing characteristics compared with the constituent single-layer two-dimensional materials. Here, we investigate the interfacial characteristics, Schottky contact, and optical performance of graphene/Ga2SSe van der Waals (vdW) heterostructure using first-principles calculations. The effects of stacking patterns, electric gating, and interlayer coupling on the interfacial properties of graphene/Ga2SSe heterostructures are also examined. Our results demonstrate that the Dirac cone of graphene is well preserved at the F point in all stacking patterns due to the weak vdW interactions, which keep the heterostructures feasible such that they can be obtained in further experiments. Moreover, depending on the stacking patterns, a small band gap of about 13-17 meV opens in graphene and has a high carrier mobility, indicating that the graphene/Ga2SSe heterostructures are potential candidates for future high-speed nanoelectronic applications. In the ground state, the graphene/Ga2SSe heterostructures form an n-type Schottky contact. The transformation from an n-type to a p-type Schottky contact or to an Ohmic contact can be forced by electric gating or by varying the interlayer coupling. Our findings could provide physical guidance for designing controllable Schottky nanodevices with high electronic and optical performances.

Abstract: Atomic scale defects in semiconductors enable their technological applications and realization of different quantum states. Using scanning tunneling microscopy and spectroscopy complemented by ab initio calculations we determine the nature of defects in the anisotropic van der Waals layered semiconductor ReS2. We demonstrate the in-plane anisotropy of the lattice by directly visualizing chains of rhenium atoms forming diamond-shaped clusters. Using scanning tunneling spectroscopy we measure the semiconducting gap in the density of states. We reveal the presence of lattice defects and by comparison of their topographic and spectroscopic signatures with ab initio calculations we determine their origin as oxygen atoms absorbed at lattice point defect sites. These results provide an atomic-scale view into the semiconducting transition metal dichalcogenides, paving the way toward understanding and engineering their properties.

Abstract: By using density functional theory (DFT) calculations, we predict a puckered, dynamically stable Janus single-layer black arsenic-phosphorus (b-AsP), which is composed of two different atomic sublayers, arsenic and phosphorus atoms. The calculated phonon spectrum reveals that Janus single-layer b-AsP is dynamically stable with either pure or coupled optical phonon branches arising from As and P atoms. The calculated Raman spectrum indicates that due to the relatively strong P-P bonds, As atoms have no contribution to the highfrequency optical vibrations. In addition, the orientation-dependent isovolume heat capacity reveals anisotropic contributions of LA and TA phonon branches to the low-temperature thermal properties. Unlike pristine single layers of b-As and b-P, Janus single-layer b-AsP exhibits additional out-of-plane asymmetry which leads to important consequences for its electronic, optical, and elastic properties. In contrast to single-layer b-As, Janus single-layer b-AsP is found to possess a direct band gap dominated by the P atoms. Moreover, real and imaginary parts of the dynamical dielectric function, including excitonic effects, reveal the highly anisotropic optical feature of the Janus single-layer. A tight-binding (TB) model is also presented for Janus single-layer b-AsP, and it is shown that, with up to seven nearest hoppings, the TB model reproduces well the DFT band structure in the low-energy region around the band gap. This TB model can be used in combination with the Green's function approach to study, e.g., quantum transport in finite systems based on Janus single-layer b-AsP. Furthermore, the linear-elastic properties of Janus single-layer b-AsP are investigated, and the orientation-dependent in-plane stiffness and Poisson ratio are calculated. It is found that the Janus single layer exhibits strong in-plane anisotropy in its Poisson ratio much larger than that of single-layer b-P. This Janus single layer is relevant for promising applications in optical dichroism and anisotropic nanoelasticity.

Abstract: High tunability of two-dimensional magnetic materials (by strain, gating, heterostructuring, or otherwise) provides unique conditions for studying versatile magnetic properties and controlling emergent magnetic phases. Expanding the scope of achievable magnetic phenomena in such materials is important for both fundamental and technological advances. Here we perform atomistic spin-dynamics simulations to explore the (chiral) magnetic phases of atomic monolayers in the limit of suppressed first-neighbors exchange interaction. We report the rich phase diagram of exotic magnetic configurations, obtained for both square and honeycomb lattice symmetries, comprising coexistence of ferromagnetic and antiferromagnetic spin cycloids, as well as multiple types of magnetic skyrmions. We perform a minimum-energy path analysis for the skyrmion collapse to evaluate the stability of such topological objects and reveal that magnetic monolayers could be good candidates to host the antiferromagnetic skyrmions that are experimentally evasive to date.

Abstract: The effect of proximity-induced interactions such as Rashba spin-orbit coupling (SOC) and exchange interaction on the electronic and optical properties of n-type monolayer (ML) MoS2 is investigated. We predict and demonstrate that the Rashba SOC can induce an in-plane spin splitting with terahertz (THz) energy, while the exchange interaction lifts the energy degeneracy in different valleys. Thus, spin polarization can be achieved in an n-type ML MoS2 and valley Hall or optical Hall effect can be observed using linearly polarized THz radiation. In such a case, the transverse optical conductivity sigma(xy) (omega) results from spin-flip transition within spin-split conduction bands and from the fact that contributions from electrons with different spin orientations in different valleys can no longer be canceled out. Interestingly, we find that for fixed effective Zeeman field (or exchange interaction) the lowest spin-split conduction band in ML MoS2 can be tuned from one in the K valley to another one in the K' valley by varying the Rashba parameter lambda(R). Therefore, by changing lambda(R) we can turn the sign of the spin polarization and Im sigma(xy) (omega) from positive to negative. Moreover, we find that the dominant contribution of the selection rules to sigma(xx)(omega) is from electrons in the K valley and to sigma(xy) (omega) is from electrons in the K' valley. These important and interesting theoretical findings can be helpful to experimental observation of the optical Hall effect in valleytronic systems using linearly polarized THz radiation fields.

Abstract: We report the direct measurement of the topological skyrmion energy barrier through a hysteresis of the skyrmion lattice in the chiral magnet MnSi. Measurements were made using small-angle neutron scattering with a custom-built resistive coil to allow for high-precision minor hysteresis loops. The experimental data were analyzed using an adapted Preisach model to quantify the energy barrier for skyrmion formation and corroborated by the minimum-energy path analysis based on atomistic spin simulations. We reveal that the skyrmion lattice in MnSi forms from the conical phase progressively in small domains, each of which consisting of hundreds of skyrmions, and with an activation barrier of several eV.

Abstract: Using first-principle calculations, we systematically study strain tuning of the electronic properties of two- dimensional graphitic carbon nitride nanosheets with empirical formula CnNm. We found the following: (i) the ferromagnetic ordered state in the metal-free systems (n, m) = (4,3), (10,9), and (14,12) remains stable in the presence of strain of about 6%. However, the system (9,7) loses its ferromagnetic ordering when increasing strain. This is due to the presence of topological defects in the (9,7) system, which eliminates the asymmetry between spin up and spin down of the p(z) orbitals when strain is applied. (ii) By applying uniaxial strain, a band gap opens in systems which are initially gapless. (iii) In semiconducting systems which have an initial gap of about 1 eV, the band gap is closed with applying uniaxial strain.

Abstract: The strain fields of periodically buckled graphene induce a periodic pseudomagnetic field (PMF) that modifies the electronic band structure. From the geometry, amplitude, and period of the periodic pseudomagnetic field, we determine the necessary conditions to access the regime of correlated phases by examining the band flattening. As compared to twisted bilayer graphene the proposed system has the advantages that (1) only a single layer of graphene is needed, (2) one is not limited to hexagonal superlattices, and (3) narrower flat bandwidth and larger separation between flat bands can be induced. We, therefore, propose that periodically strained graphene single layers can become a platform for the exploration of exotic many-body phases.

Abstract: Cold atmospheric pressure plasmas have proven to provide an alternative treatment of cancer by targeting tumorous cells while leaving their healthy counterparts unharmed. However, the underlying mechanisms of the plasma–cell interactions are not yet fully understood. Reactive oxygen species, and in particular hydroxyl radicals (OH), are known to play a crucial role in plasma driven apoptosis of
malignant cells. In this paper we investigate the interaction of OH radicals, as well as H2O2 molecules and HO2 radicals, with DNA by means of reactive molecular dynamics simulations using the ReaxFF force field. Our results provide atomic-scale insight into the dynamics of oxidative stress on DNA caused by the OH radicals, while H2O2 molecules appear not reactive within the considered timescale. Among the observed processes are the formation of 8-OH-adduct radicals, forming the first stages towards the formation of 8-oxoGua and 8-oxoAde, H-abstraction reactions of the amines, and the partial opening of loose DNA ends in aqueous solution.

Abstract: We analyzed pattern formation and identified various morphologies in a system of particles interacting through a non-monotonic potential with a competing range interaction characterized by a repulsive core (r < r(c)) and an attractive tail (r > r(c)), using molecular-dynamics simulations. Depending on parameters, the interaction potential models the inter-particle interaction in various physical systems ranging from atoms, molecules and colloids to vortices in low kappa type-II superconductors and in recently discovered 'type-1.5' superconductors. We constructed a 'morphology diagram' in the plane 'critical radius r(c)-density n' and proposed a new approach to characterizing the different types of patterns. Namely, we elaborated a set of quantitative criteria in order to identify the different pattern types, using the radial distribution function (RDF), the local density function and the occupation factor.