Abstract: Motivated by the recent realization of graphene sensors to detect individual gas molecules, we investigate the adsorption of H2O, NH3, CO, NO2, and NO on a graphene substrate using first-principles calculations. The optimal adsorption position and orientation of these molecules on the graphene surface is determined and the adsorption energies are calculated. Molecular doping, i.e., charge transfer between the molecules and the graphene surface, is discussed in light of the density of states and the molecular orbitals of the adsorbates. The efficiency of doping of the different molecules is determined and the influence of their magnetic moment is discussed.

Abstract: The envelope function method traditionally employs a single basis set which, in practice, relates to a single material because the k.p matrix elements are generally only known in a particular basis. In this work, we defined a basis function transformation to alleviate this restriction. The transformation is completely described by the known inter-band momentum matrix elements. The resulting envelope function equation can solve the electronic structure in lattice matched heterostructures without resorting to boundary conditions at the interface between materials, while all unit-cell averaged observables can be calculated as with the standard envelope function formalism. In the case of two coupled bands, this heterostructure formalism is equivalent to the standard formalism while taking position dependent matrix elements. (C) 2015 Elsevier B.V. All rights reserved

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.

Abstract: Using the four-band continuum model we derive a general expression for the infinite-mass boundary condition in bilayer graphene. Applying this new boundary condition we analytically calculate the confined states and the corresponding wave functions in a bilayer graphene quantum dot in the absence and presence of a perpendicular magnetic field. Our results for the energy spectrum show an energy gap between the electron and hole states at small magnetic fields. Furthermore the electron (e) and hole (h) energy levels corresponding to the K and K' valleys exhibit the E-K(e(h)) (m) = E-K'(e(h)) (m) symmetry, where m is the angular momentum quantum number. (C) 2014 Elsevier Ltd. All rights reserved.

Abstract: Molecular-dynamic simulations were performed in order to investigate the melting processes of isotropically confined binary systems. We considered two species of particles, which differ by their amount of electric charge. A Lindemann type of criterion was used to determine the angular melting temperature. We demonstrate that the magic-to-normal cluster transition can evolve in two distinct ways, that is, through a structural phase transition of the first order or via a smooth transition where an increase of the shells' width leads to a continuous decreasing mechanical stability of the system. Moreover, for large systems, we demonstrate that the internal cluster exerts a minor effect on the mechanical stability of the external shell. Furthermore, we show that highly symmetric configurations, such as those found for multiple ring structures, have large mechanical stability, i.e., high angular melting temperature.

Abstract: We extend the unsymmetrized self-consistent-field method (USF) for anharmonic crystals to layered non-Bravais crystals to investigate structural, dynamical and thermodynamic properties of a free-standing graphene monolayer. In this theory, the main anharmonicity of the crystal lattice has been included and the quantum corrections are taken into account in an h-expansion for the one-particle density matrix. The obtained result for the thermal expansion coefficient (TEC) of graphene shows a strong temperature dependence and agrees with experimental results by Bao et al. (Nat. Nanotechnol., 4 (2009) 562). The obtained value of TEC at room temperature (300 K) is -6.4 x 10(- 6) K- 1 and it becomes positive for T > T-alpha = 358K. We find that quantum effects are significant for T < 1000 K. The interatomic distance, effective amplitudes of the graphene lattice vibrations, adiabatic and isothermal bulk moduli, isobaric and isochoric heat capacities are also calculated and their temperature dependences are determined. Copyright (C) EPLA, 2014

Abstract: Magnetically coupled superconductorferromagnet hybrids offer advanced routes for nanoscale control of superconductivity. Magnetotransport characteristics and scanning tunneling microscopy images of vortex structures in superconductorferromagnet hybrids reveal rich superconducting phase diagrams. Focusing on a particular combination of a ferromagnet with a well-ordered periodic magnetic domain structure with alternating out-of-plane component of magnetization, and a small coherence length superconductor, we find directed nucleation of superconductivity above the domain wall boundaries. We show that near the superconductor-normal state phase boundary the superconductivity is localized in narrow mesoscopic channels. In order to explore the Abrikosov flux line ordering in F/S hybrids, we use a combination of scanning tunneling microscopy and GinzburgLandau simulations. The magnetic stripe domain structure induces periodic local magnetic induction in the superconductor, creating a series of pinninganti-pinning channels for externally added magnetic flux quanta. Such laterally confined Abrikosov vortices form quasi-1D arrays (chains). The transitions between multichain states occur through propagation of kinks at the intermediate fields. At high fields we show that the system becomes nonlinear due to a change in both the number of vortices and the confining potential. In F/S/F hybrids we demonstrate the evolution of the anisotropic conductivity in the superconductor that is magnetically coupled with two adjacent ferromagnetic layers. Stripe magnetic domain structures in both F-layers are aligned under each other, resulting in a directional superconducting order parameter in the superconducting layer. The conductance anisotropy strongly depends on the period of the magnetic domains and the strength of the local magnetization. The anisotropic conductivity of up to three orders of magnitude can be achieved with a spatial critical temperature modulation of 5% of T c. Induced anisotropic properties in the F/S and F/S/F hybrids have a potential for future application in switching and nonvolatile memory elements operating at low temperatures.

Abstract: Themean-field Hubbard model is used to investigate the formation of the antiferromagnetic phase in hexagonal graphene rings with inner zigzag edges. The outer edge of the ring was taken to be either zigzag or armchair, and we found that both types of structures can have a larger antiferromagnetic interaction as compared with hexagonal dots. This difference could be partially ascribed to the larger number of zigzag edges per unit area in rings than in dots. Furthermore, edge states localized on the inner ring edge are found to hybridize differently than the edge states of dots, which results in important differences in the magnetism of graphene rings and dots. The largest staggered magnetization is found when the outer edge has a zigzag shape. However, narrow rings with armchair outer edge are found to have larger staggered magnetization than zigzag hexagons. The edge defects are shown to have the least effect on magnetization when the outer ring edge is armchair shaped. DOI: 10.1103/PhysRevB.87.085434

Abstract: The technique of beam optics is applied to the description of the wave function of Dirac electrons. This approach is illustrated by considering electron transmission through simple nonhomogeneous structures, such as flat and bent p-n junctions and superlattices. We found that a convex p-n junction compresses the beam waist, while a concave interface widens it without loosing its focusing properties. At a flat p-n junction the waist of the transmitted Gaussian beam can be narrowed or widened, depending on the angle of incidence. A general condition is derived for the occurrence of beam collimation in a superlattice which is less stringent than previous discussed.

Abstract: In this study we investigate the effect of device topology on the ballistic current in n-channel metal-oxide-semiconductor field-effect transistors. Comparison of the nanoscale planar and double-gate devices reveals that, down to a certain thickness of the double gate film, the ballistic current flowing in the double gate device is twice as large compared to its planar counterpart. On the other hand, further thinning of the film beyond this threshold is found to change noticeably the confinement and transport characteristics, which are strongly depending on the film material and the surface orientation. For double gate Ge and Si devices there exists a critical film thickness below which the transverse gate field is no longer effectively screened by the inversion layer electron gas and mutual inversion of the two gates is turned on. In the case of GaAs and other similar IIIV compounds, a decrease in the film thickness may drastically change the occupation of the L-valleys and therefore amend the transport properties. The simulation results show that, in both cases, the ballistic current and the transconductance are considerably enhanced.

Abstract: First-principles density functional theory (DFT) calculations were carried out to investigate the structural and electronic properties of beryllium (Be) doped and Be and boron (B) co-doped graphene systems. We observed that not only the concentration of impurity atoms is important to tune the band-gap to some desired level, but also the specific substitution sites play a key role. In our system, which consists of 32 atoms, a maximum of 4Be and, in the co-doped state, 2Be and 3B atom substitutions are investigated. Both dopants are electron deficient relative to C atoms and cause the Fermi level to shift downward (p-type doping). A maximum band gap of 1.44 eV can be achieved on incorporation of 4Be atoms. The introduction of Be is more sensitive in terms of geometry and stability than B. However, in opening the energy gap, Be is more effective than B and N (nitrogen). Our results offer the possibility to modify the band-gap of graphene sufficiently for utilization in diverse electronic device applications.