Abstract: Motivated by a recent experiment that reported the synthesis of a new 2D material nitrogenated holey graphene (C2N) [Mahmood et al., Nat. Commun., 2015, 6, 6486], the electronic, magnetic, and mechanical properties of nitrogenated (C2N), phosphorated (C2P) and arsenicated (C2As) monolayer holey graphene structures are investigated using first-principles calculations. Our total energy calculations indicate that, similar to the C2N monolayer, the formation of the other two holey structures are also energetically feasible. Calculated cohesive energies for each monolayer show a decreasing trend going from the C2N to C2As structure. Remarkably, all the holey monolayers considered are direct band gap semiconductors. Regarding the mechanical properties (in-plane stiffness and Poisson ratio), we find that C2N has the highest in-plane stiffness and the largest Poisson ratio among the three monolayers. In addition, our calculations reveal that for the C2N, C2P and C2As monolayers, creation of N and P defects changes the semiconducting behavior to a metallic ground state while the inclusion of double H impurities in all holey structures results in magnetic ground states. As an alternative to the experimentally synthesized C2N, C2P and C2As are mechanically stable and flexible semiconductors which are important for potential applications in optoelectronics.

Abstract: Recently it was found that bilayer graphene may exhibit regions with and without van der Waals coupling between the two layers. We show that such structures can exhibit a strong layer selectivity when current flows through the coupled region and that this selectivity can be tuned by means of electrostatic gating. Analysing how this effect depends on the type of bilayer stacking, the potential on the gates and the smoothness of the boundary between the coupled and decoupled regions, we show that nearly perfect layer selectivity is achievable in these systems. This effect can be further used to realise a tunable layer switch.

Abstract: Using the tight-binding approach, we study the influence of a perpendicular magnetic field on the energy levels of hexagonal, triangular, and circular bilayer graphene (BLG) quantum dots (QDs) with zigzag and armchair edges. We obtain the energy levels for AB (Bernal)-stacked BLG QDs in both the absence and the presence of a perpendicular electric field (i.e., biased BLG QDs). We find different regions in the spectrum of biased QDs with respect to the crossing point between the lowest-electron and -hole Landau levels of a biased BLG sheet. Those different regions correspond to electron states that are localized at the center, edge, or corner of the BLG QD. Quantum Hall corner states are found to be absent in circular BLG QDs. The spatial symmetry of the carrier density distribution is related to the symmetry of the confinement potential, the position of zigzag edges, and the presence or absence of interlayer inversion symmetry.

Abstract: Finding novel atomically thin heterostructures and understanding their characteristic properties are critical for developing better nanoscale optoelectronic devices. In this study, we investigate the electronic and optical properties of a GaS-Ca(OH)(2) heterostructure using first-principle calculations. The band gap of the GaS-Ca(OH)(2) heterostructure is significantly reduced when compared to those of the isolated constituent layers. Our calculations showthat the GaS-Ca(OH)(2) heterostructure is a type-II heterojunction which can be used to separate photoinduced charge carriers where electrons are localized in GaS and holes in the Ca(OH)(2) layer. This leads to spatially indirect excitons which are important for solar energy and optoelectronic applications due to their long lifetime. By solving the Bethe-Salpeter equation on top of a single shot GW calculation (G(0)W(0)), the dielectric function and optical oscillator strength of the constituent monolayers and the heterostructure are obtained. The oscillator strength of the optical transition for the GaS monolayer is an order of magnitude larger than the Ca(OH)(2) monolayer. We also found that the calculated optical spectra of different stacking types of the heterostructure show dissimilarities, although their electronic structures are rather similar. This prediction can be used to determine the stacking type of ultrathin heterostructures.

Abstract: We investigate the superconducting states that arise as a consequence of mesoscopic confinement and a multicomponent order parameter in the Ginzburg-Landau model for p-wave superconductivity. Conventional vortices, but also half-quantum vortices and skyrmions, are found as the applied magnetic field and the anisotropy parameters of the Fermi surface are varied. The solutions are well differentiated by a topological charge that for skyrmions is given by the Hopf invariant and for vortices by the circulation of the superconducting velocity. We revealed several unique states combining vortices and skyrmions, their possible reconfiguration with varied magnetic field, as well as temporal and field-induced transitions between vortical and skyrmionic states.

Abstract: In order to study quantum effects in a two-dimensional crystal lattice of a free-standing monolayer graphene, we have performed both path-integral Monte Carlo (PIMC) and classical Monte Carlo (MC) simulations for temperatures up to 2000 K. The REBO potential is used for the interatomic interaction. The total energy, interatomic distance, root-mean-square displacement of the atom vibrations, and the free energy of the graphene layer are calculated. The obtained lattice vibrational energy per atom from the classical MC simulation is very close to the energy of a three-dimensional harmonic oscillator 3k(B)T. The PIMC simulation shows that quantum effects due to zero-point vibrations are significant for temperatures T < 1000 K. The quantum contribution to the lattice vibrational energy becomes larger than that of the classical lattice for T < 400 K. The lattice expansion due to the zero-point motion causes an increase of 0.53% in the lattice parameter. A minimum in the lattice parameter appears at T similar or equal to 500 K. Quantum effects on the atomic vibration amplitude of the graphene lattice and its free energy are investigated.

Abstract: The thermal expansion alpha(T) in layered crystals is of fundamental and technological interest. As suggested by I. M. Lifshitz in 1952, in thin solid films (crystalline membranes) a negative contribution to alpha(T) is due to anharmonic couplings between in-plane stretching modes and out-of-plane bending (flexural modes). Genuine in-plane anharmonicities give a positive contribution to alpha(T). The competition between these two effects can lead to a change of sign (crossover) from a negative value of alpha(T) in a temperature (T) range T <= T-alpha to a positive value of alpha(T) for T > T-alpha in layered crystals. Here, we present an analytical lattice dynamical theory of these phenomena for a two-dimensional (2D) hexagonal crystal. We start from a Hamiltonian that comprises anharmonic terms of third and fourth order in the lattice displacements. The in-plane and out-of-plane contributions to the thermal expansion are studied as functions of T for crystals of different sizes. Besides, renormalization of the flexural mode frequencies plays a crucial role in determining the crossover temperature T-alpha. Numerical examples are given for graphene where the anharmonic couplings are determined from experiments. The theory is applicable to other layer crystals wherever the anharmonic couplings are known. (C) 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim

Abstract: We study ballistic transport of Dirac fermions in silicene through arrays of barriers, of width d, in the presence of an exchange field M and a tunable potential of height U or depth-U. The spin-and valley-resolved conductances as functions of U or M, exhibit resonances away from the Dirac point (DP) and close to it a pronounced dip that becomes a gap when a critical electric field E-z is applied. This gap widens by increasing the number of barriers and can be used to realize electric field-controlled switching of the current. The spin p(s) and valley p(v) polarizations of the current near the DP increase with Ez or M and can reach 100% for certain of their values. These field ranges widen significantly by increasing the number of barriers. Also, ps and pv oscillate nearly periodically with the separation between barriers or wells and can be inverted by reversing M.

Abstract: We study the effect of mechanical deformations on the elec- tronic properties of hexagonal flakes of bilayer graphene. The behavior of electrons induced by triaxial strain can be de- scribed by an effective pseudo-magnetic field which is homo- geneous in the center of the flake. We find that in-plane strain, applied to both layers equally, can break the layer symmetry leading to different behavior in the top and bottom layers of graphene. At low energy, just one of the layers feels
the pseudo-magnetic field: the zero-energy pseudo-Landau level is missing in the second layer, thus creating a gap be- tween the lowest non-zero levels. While the layer asymmetry is most significant at zero energy, interaction with the edges of the flake extends the effect to higher pseudo-Landau lev- els. The behavior of the top and bottom layers may be re- versed by rotating the triaxial strain by 60°.

Abstract: We theoretically investigate the electronic transport properties of two closely spaced L-shaped semiconductor quantum wires, for different configurations of the output channel widths as well as the distance between the wires. Within the effective-mass approximation, we solve the time-dependent Schrodinger equation using the split-operator technique that allows us to calculate the transmission probability, the total probability current, the conductance, and the wave function scattering between the energy subbands. We determine the maximum distance between the quantum wires below which a relevant non-zero transmission is still found. The transmission probability and the conductance show a strong dependence on the width of the output channel for small distances between the wires. (C) 2015 AIP Publishing LLC.

Abstract: We use the anisotropic time-dependent Ginzburg-Landau theory to investigate the effect of a square array of out-of-plane magnetic dots on the dynamics of Josephson vortices (fluxons) in artificial stacks of superconducting-normal-superconducting (SNS) Josephson junctions in the presence of external DC and AC currents. Periodic pinning due to the magnetic dots distorts the triangular lattice of fluxons and results in the appearance of commensurability features in the current-voltage characteristics of the system. For the larger values of the magnetization, additional peaks appear in the voltage-time characteristics of the system due to the creation and annihilation of vortex-antivortex pairs. Peculiar changes in the response of the system to the applied current is found resulting in a “superradiant” vortex-flow state at large current values, where a rectangular lattice of moving vortices is formed. Synchronizing the motion of fluxons by adding a small ac component to the biasing dc current is realized. However, we found that synchronization becomes difficult for large magnetization of the dots due to the formation of vortex-antivortex pairs.

Abstract: Using first principles calculations we investigate the structural and electronic properties of interfaces between fluorographane and MoS2. Unsymmetrical functionalization of graphene with H and F results in an intrinsic dipole moment perpendicular to the plane of the buckled graphene skeleton. Depending on the orientation of this dipole moment, the electronic properties of a physically absorbed MoS2 monolayer can be switched from n-to p-type or vice versa. We show that one can realize vanishing n-type/p-type Schottky barrier heights when contacting MoS2 to fluorographane. By applying a perpendicular electric field, the size of the Schottky barrier and the degree of doping can be tuned. Our calculations indicate that a fluorographane monolayer is a promising candidate for bipolar doping of MoS2, which is vital in the design of novel technological applications based on two-dimensional materials.

Abstract: The structural properties of partially hydrogenated and fluorinated graphene with different percentages of H/F atoms are investigated using molecular dynamics simulations based on reactive force field (ReaxFF) potentials. We found that the roughness of graphene varies with the percentage (p) of H or F and in both cases is maximal around p = 50%. Similar results were obtained for partially oxidized graphene. The two-dimensional area size of partially fluorinated and hydrogenated graphene exhibits a local minimum around p = 35% coverage. The lattice thermal contraction in partially functionalized graphene is found to be one order of magnitude larger than that of fully covered graphene. We also show that the armchair structure for graphene oxide (similar to the structure of fully hydrogenated and fluorinated graphene) is unstable. Our results show that the structure of partially functionalized graphene changes nontrivially with the C : H and C : F ratio as well as with temperature.

Abstract: Piezoelectricity is a unique material property that allows one to convert mechanical energy into electrical one or vice versa. Transition metal dichalcogenides (TMDC) and transition metal dioxides (TMDO) are expected to have great potential for piezoelectric device applications due to their noncentrosymmetric and two-dimensional crystal structure. A detailed theoretical investigation of the piezoelectric stress (e 11 ) and piezoelectric strain (d 11 ) coefficients of single layer TMDCs and TMDOs with chemical formula MX 2 (where M= Cr, Mo, W, Ti, Zr, Hf, Sn and X = O, S, Se, Te) is presented by using first-principles calculations based on density func- tional theory. We predict that not only the Mo- and W-based members of this family but also the other materials with M= Cr, Ti, Zr and Sn exhibit highly promising piezoelectric properties. CrTe 2 has the largest e 11 and d 11 coefficients among the group VI elements (i.e., Cr, Mo, and W). In addition, the relaxed-ion e 11 and d 11 coefficients of SnS 2 are almost the same as those of CrTe 2 . Furthermore, TiO 2 and ZrO 2 pose comparable or even larger e 11 coefficients as compared to Mo- and W-based TMDCs and TMDOs. Our calculations reveal that TMDC and TMDO structures are strong candidates for future atomically thin piezoelectric applications such as transducers, sensors, and energy harvesting devices due to their piezoelectric coefficients that are comparable (even larger) to currently used bulk piezoelectric materials.

Abstract: The electron states in axially symmetric quantum wires are computed by means of the effective-mass Schrodinger equation, which is written in cylindrical coordinates phi, rho, and z. We show that a direct discretization of the Schrodinger equation by central finite differences leads to a non-symmetric Hamiltonian matrix. Because diagonalization of such matrices is more complex it is advantageous to transform it in a symmetric form. This can be done by the Liouville-like transformation proposed by Rizea et al. (2008), which replaces the wave function psi(rho) with the function F(rho) = psi(rho)root rho and transforms the Hamiltonian accordingly. Even though a symmetric Hamiltonian matrix is produced by this procedure, the computed wave functions are found to be inaccurate near the origin, and the accuracy of the energy levels is not very high. In order to improve on this, we devised a finite-difference scheme which discretizes the Schrodinger equation in the first step, and then applies the Liouville-like transformation to the difference equation. Such a procedure gives a symmetric Hamiltonian matrix, resulting in an accuracy comparable to the one obtained with the finite element method. The superior efficiency of the new finite-difference scheme (FDM) is demonstrated for a few p-dependent one-dimensional potentials which are usually employed to model the electron states in free-standing and core shell quantum wires. The new scheme is compared with the other FDM schemes for solving the effective-mass Schrodinger equation, and is found to deliver energy levels with much smaller numerical error for all the analyzed potentials. It also gives more accurate results than the scheme of Rizea et al., except for the ground state of an infinite rectangular potential in freestanding quantum wires. Moreover, the PT symmetry is invoked to explain similarities and differences between the considered FDM schemes. (C) 2015 Elsevier B.V. All rights reserved.