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: Janus particles present an important class of building blocks for directional assembly. These are compartmentalized colloids with two different hemispheres. In this work we consider a three-dimensional model of Janus spheres that contain one hydrophobic and one charged hemisphere. Using molecular dynamics simulations, we study the morphology of these particles when confined in a channel-like environment. The interplay between the attractive and repulsive forces on each particle gives rise to a rich phase space where the relative orientation of each particle plays a dominant role in the formation of large-scale clusters. The interest in this system is primarily due to the fact that it could give a better understanding of the mechanisms of the formation of polar membranes. A variety of ordered membranelike morphologies is found consisting of single and multiple connected chain configurations. The helicity of these chains can be chosen by simply changing the salt concentration of the solution. Special attention is given to the formation of Bernal spirals. These helices are composed of regular tetrahedra and are known to exhibit nontrivial translational and rotational symmetry.

Abstract: Using density functional theory, we study the electronic properties of a graphene-silicene bilayer (GSB). A single layer of silicene binds to the graphene layer with adhesion energy of about 25 meV/atom. This adhesion energy between the two layers follows accurately the well-known -1/z(2) dispersion energy as found between two infinite parallel plates. In small flakes of GSB with hydrogenated edges, negative charge is transferred from the graphene layer to the silicene layer, producing a permanent and a switchable polar bilayer, while in an infinite GSB, the negative charge is transferred from the silicene layer to the graphene layer. The graphene-silicene bilayer is a good candidate for a nanocapacitor with piezoelectric capabilities. We found that the permanent dipole of the bilayer can be tuned by an external perpendicular electric field.

Abstract: Hydrogen adsorption on boron and nitrogen doped graphene is investigated in detail by means of first-principles calculations. A comprehensive study is performed of the structural, electronic, and magnetic properties of chemisorbed hydrogen atoms and atom pairs near the dopant sites. The main effect of the substitutional atoms is charge doping which is found to greatly affect the adsorption process by increasing the binding energy at the sites closest to the substitutional species. It is also found that doping does not induce magnetism despite the odd number of electrons per atom introduced by the foreign species, and that it quenches the paramagnetic response of chemisorbed H atoms on graphene. Overall, the effects are similar for B and N doping, with only minor differences in the adsorption energetics due to different sizes of the dopant atoms and the accompanying lattice distortions.

Abstract: The in-plane stiffness (C), Poisson's ratio (nu), Young's modulus and ultimate strength (sigma) along two different crystallographic orientations are calculated for the single layer crystals: MoS2, HfS2 and TiS3 in 1H, 1T and monoclinic phases. We find that MoS2 and HfS2 have isotropic in-plane stiffnesses of 124.24 N m(-1) and 79.86 N m(-1), respectively. While for TiS3 the in-plane stiffness is highly anisotropic due to its monoclinic structure, with C-x = 83.33 N m(-1) and C-y = 133.56 N m(-1) (x and y are parallel to its longer and shorter in-plane lattice vectors.). HfS2 which is in the 1T phase has the smallest anisotropy in its ultimate strength, whereas TiS3 in the monoclinic phase has the largest. Along the armchair direction MoS2 has the largest sigma of 23.48 GPa, whereas along y TiS3 has the largest sigma of 18.32 GPa. We have further analyzed the band gap response of these materials under uniaxial tensile strain, and find that they exhibit different behavior. Along both armchair and zigzag directions, the band gap of MoS2 (HfS2) decreases (increases) as strain increases, and the response is almost isotropic. For TiS3, the band gap decreases when strain is along x, while if strain is along y, the band gap increases first and then decreases beyond a threshold strain value. The different characteristics observed in these sulphides with different structures shed light on the relationship between the structure and properties, which is useful for applications in nanotechnology.

Abstract: Within a tight binding approach we investigate the energy levels of hexagonal and triangular bilayer graphene (BLG) quantum dots (QDs) with zigzag and armchair edges. We study AA- and AB-(Bernal) stacked BLG QDs and obtain the energy levels in both the absence and the presence of a perpendicular electric field (i.e., biased BLG QDs). Our results show that the size dependence of the energy levels is different from that of monolayer graphene QDs. The energy spectrum of AB-stacked BLG QDs with zigzag edges exhibits edge states which spread out into the opened energy gap in the presence of a perpendicular electric field. We found that the behavior of these edges states is different for the hexagonal and triangular geometries. In the case of AA-stacked BLG QDs, the electron and hole energy levels cross each other in both cases of armchair and zigzag edges as the dot size or the applied bias increases.

Abstract: We have grown epitaxial Cr-doped V2O3 thin films with Cr concentrations between 0% and 20% on (0001)-Al2O3 by oxygen-assisted molecular beam epitaxy. For the highly doped samples (>3%), a regular and monotonous increase of the resistance with decreasing temperature is measured. Strikingly, in the low doping samples (between 1% and 3%), a collapse of the insulating state is observed with a reduction of the low temperature resistivity by up to 5 orders of magnitude. A vacuum annealing at high temperature of the films recovers the low temperature insulating state for doping levels below 3% and increases the room temperature resistivity towards the values of Cr-doped V2O3 single crystals. It is well-know that oxygen excess stabilizes a metallic state in V2O3 single crystals. Hence, we propose that Cr doping promotes oxygen excess in our films during deposition, leading to the collapse of the low temperature insulating state at low Cr concentrations. These results suggest that slightly Cr-doped V2O3 films can be interesting candidates for field effect devices. (C) 2015 AIP Publishing LLC.

Abstract: The properties of a new class of phosphorene nanotubes (PNT) are investigated by performing first-principles calculations. We demonstrate that it is advantageous to use blue phosphorene in order to make small nanotubes and propose a way to create low-energy PNTs by the inclusion of defect lines. Five different types of defect lines are investigated and incorporated in various combinations. The resulting defect-induced faceted PNTs have negligible bending stresses which leads to a reduction in the formation energy with respect to round PNTs. Our armchair faceted PNTs have similar formation energies than the recently proposed multiphase faceted PNTs, but they have a larger variety of possible structures. Our zigzag faceted PNTs have lower formation energies than round tubes and multiphase faceted nanotubes. The electronic properties of the defect-induced faceted PNTs are determined by the defect lines which control the band gap and the shape of the electronic states at the band edges. These band gaps increase with the radius of the nanotubes and converge to those of isolated defect lines.

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: The structural transitions of the ground state of a system of repulsively interacting particles confined in a quasi-one-dimensional channel, and the effect of the interparticle interaction as well as the functional form of the confinement potential on those transitions are investigated. Although the nonsequential ordering of transitions (non-SOT), i.e., the 1 – 2 – 4 – 3 – 4 – 5 – 6 – ... sequence of chain configurations with increasing density, is widely robust as predicted in a number of theoretical studies, the sequential ordering of transitions (SOT), i.e., the 1 – 2 – 3 – 4 – 5 – 6 – ... chain, is found as the ground state for long-ranged interparticle interaction and hard-wall-like confinement potentials. We found an energy barrier between every two different phases around its transition point, which plays an important role in the preference of the system to follow either a SOT or a non-SOT. However, that preferential transition requires also the stability of the phases during the transition. Additionally, we analyze the effect of a small structural disorder on the transition between the two phases around its transition point. Our results show that a small deformation of the triangular structure changes dramatically the picture of the transition between two phases, removing in a considerable region the non-SOT in the system. This feature could explain the fact that the non-SOT is, up to now, not observed in experimental systems, and suggests a more advanced experimental setup to detect the non-SOT.

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

Abstract: Superconductors, ideally diamagnetic when in the Meissner state, can also exhibit paramagnetic behavior due to trapped magnetic flux. In the absence of pinning such paramagnetic response is weak, and ceases with increasing sample thickness. Here we show that in multiband superconductors paramagnetic response can be observed even in slab geometries, and can be far larger than any previous estimate – even multiply larger than the diamagnetic Meissner response for the same applied magnetic field. We link the appearance of this giant paramagnetic response to the broad crossover between conventional Type-I and Type-II superconductors, where Abrikosov vortices interact non-monotonically and multibody effects become important, causing unique flux configurations and their locking in the presence of surfaces.

Abstract: A comparative study is made of the propagation of wave packets and fronts in three different meta-media, i.e. graphene, left-handed media (LHM) and transmission lines, using one-dimensional models. It is shown that a potential step in graphene influences only the frequency of the electronic wave, i.e., the particular spectrum branch (electron or hole) to which the wave belongs to, while the envelop function (the wave front or packet form) remains unchanged. Although the model for a vacuum and LHM interface is similar to that of the potential step in graphene, the solutions are quite different due to differences in the chirality of the waves. Comparing the propagation of wave fronts and packets in a standard transmission line and its meta-analog we demonstrate that the propagating packets in the meta-line are much more deformed as compared to the standard one, including broadening, asymmetry and even the appearance of fast moving precursors. This influence is seen not only in the case of packets with steep fronts but in soft Gaussian packets as well.

Abstract: A Hg(Cd)Te insulator heterostructure hosts a two-dimensional electron system that can simulate the physics of Dirac fermions with only a single valley. We investigate the magnetotransport properties of this structure and show that, unlike most two-dimensional crystals with spin and valley coupled levels, the Shubnikov-de Haas oscillations exhibit a high spin polarization in the absence of any valley degree of freedom. This effect can be observed using magnetospectroscopy measurements for quantum well thicknesses corresponding to either the topologically trivial or quantum spin Hall phases. The pseudospin texture of the electrons near the Fermi level is also studied and we show that a tunable pseudospin-polarized quantum Hall liquid can only be observed for thicknesses corresponding to the inverted regime.