Abstract: A theoretical and experimental study of the work function of few-layer graphene is reported. The influence of the number of layers on the work function is investigated in the presence of a substrate, a molecular dipole layer, and combinations of the two. The work function of few-layer graphene is almost independent of the number of layers with only a difference between monolayer and multilayer graphene of about 60 meV. In the presence of a charge-donating substrate the charge distribution is found to decay exponentially away from the substrate and this is directly reflected in the work function of few-layer graphene. A dipole layer changes the work function only when placed in between the substrate and few-layer graphene through a change of the charge transfer between the two.

Abstract: Trions and biexcitons in anisotropic two-dimensional materials are investigated within an effective mass theory. Explicit results are obtained for phosphorene and arsenene, materials that share features such as a direct quasiparticle gap and anisotropic conduction and valence bands. Trions are predicted to have remarkably high binding energies and an elongated electron-hole structure with a preference for alignment along the armchair direction, where the effective masses are lower. We find that biexciton binding energies are also notably large, especially for monolayer phosphorene, where they are found to be twice as large as those for typical monolayer transition metal dichalcogenides.

Abstract: Electronic transport properties of a graphene-silicene bilayer system are studied using density-functional theory in combination with the nonequilibrium Green's function formalism. Depending on the energy of the electrons, the transmission can be larger in this system as compared to the sum of the transmissions of separated graphene and silicene monolayers. This effect is related to the increased electron density of states in the bilayer sample. At some energies, the electronic states become localized in one of the layers, resulting in the suppression of the electron transmission. The effect of an applied voltage on the transmission becomes more pronounced in the layered sample as compared to graphene due to the larger variation of the electrostatic potential profile. Our findings will be useful when creating hybrid nanoscale devices where enhanced transport properties will be desirable. (C) 2015 AIP Publishing LLC.

Abstract: Both even- and odd-numbered neutral carbon clusters Cn (n=210) are systematically studied using the energy minimization method and the modified Brenner potential for the carbon-carbon interactions. Many stable configurations were found, and several new isomers are predicted. For the lowest energy stable configurations we obtained their binding energies and bond lengths. We found that for n5 the linear isomer is the most stable one while for n>5 the monocyclic isomer becomes the most stable. The latter was found to be regular for all studied clusters. The dependence of the binding energy for linear and cyclic clusters versus the cluster size n (n=210) is found to be in good agreement with several previous calculations, in particular with ab initio calculations as well as with experimental data for n=25.

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: The self rotation of a graphene flake over graphite is controlled by the size, initial misalignment and temperature. Using both ab initio calculations and molecular dynamics simulations, we investigate annealing effects on the self rotation of a graphene flake on a graphene substrate. The energy barriers for rotation and drift of a graphene flake over graphene is found to be smaller than 25 meV/atom which is comparable to thermal energy. We found that small flakes (of about similar to 4 nm) are more sensitive to temperature and initial misorientation angles than larger one (beyond 10 nm). The initial stacking configuration of the flake is found to be important for its dynamics and time evolution of misalignment. Large flakes, which are initially in the AA-or AB-stacking state with small misorientation angle, rotate and end up in the AB-stacking configuration. However small flakes can they stay in an incommensurate state specially when the initial misorientation angle is larger than 2 degrees. Our results are in agreement with recent experiments.

Abstract: Different allotropes of phosphorene are possible of which black and blue phosphorus are the most stable. While blue phosphorus has isotropic properties, black phosphorus is strongly anisotropic in its electronic and optical properties due to its anisotropic crystal structure. In this work, we systematically investigated the lattice thermal properties of black and blue phosphorene by using first-principles calculations based on the quasiharmonic approximation approach. Similar to the optoelectronic and electronic properties, we predict that black phosphorene has highly anisotropic thermal properties, in contrast to the blue phase. The linear thermal expansion coefficients along the zigzag and armchair direction differ up to 20% in black phosphorene. The armchair direction of black phosphorene is more expandable as compared to the zigzag direction and the biaxial expansion of blue phosphorene under finite temperature. Our comparative analysis reveals that the inclusion of finite-temperature effects makes the blue phase thermodynamically more stable over the black phase above 135 K.

Abstract: Thermal fluctuations of single layer hydrogenated graphene (graphane) are investigated using large scale atomistic simulations. By analyzing the mean square value of the height fluctuations < h(2)> and the height-height correlation function H(q) for different system sizes and temperatures, we show that hydrogenated graphene is an unrippled system in contrast to graphene. The height fluctuations are bounded, which is confirmed by a H(q) tending to a constant in the long wavelength limit instead of showing the characteristic scaling law q(4-eta)(eta similar or equal to 0.85) predicted by membrane theory. This unexpected behavior persists up to temperatures of at least 900 K and is a consequence of the fact that in graphane the thermal energy can be accommodated by in-plane bending modes, i.e., modes involving C-C-C bond angles in the buckled carbon layer, instead of leading to significant out-of-plane fluctuations that occur in graphene.

Abstract: Using atomistic simulations we investigate the thermodynamical properties of a single atomic layer of hexagonal boron nitride (h-BN). The thermal induced ripples, heat capacity, and thermal lattice expansion of large scale h-BN sheets are determined and compared to those found for graphene (GE) for temperatures up to 1000 K. By analyzing the mean-square height fluctuations < h(2)> and the height-height correlation function H(q) we found that the h-BN sheet is a less stiff material as compared to graphene. The bending rigidity of h-BN (i) is about 16% smaller than the one of GE at room temperature (300 K), and (ii) increases with temperature as in GE. The difference in stiffness between h-BN and GE results in unequal responses to external uniaxial and shear stress and different buckling transitions. In contrast to a GE sheet, the buckling transition of a h-BN sheet depends strongly on the direction of the applied compression. The molar heat capacity, thermal-expansion coefficient, and Gruneisen parameter are estimated to be 25.2 J mol(-1) K-1, 7.2 x 10(-6) K-1, and 0.89, respectively.

Abstract: Starting from the simplified linear combination of atomic orbitals method in combination with first-principles calculations, we construct a tight-binding (TB) model in the two-centre approximation for borophene and hydrogenated borophene (borophane). The Slater and Koster approach is applied to calculate the TB Hamiltonian of these systems. We obtain expressions for the Hamiltonian and overlap matrix elements between different orbitals for the different atoms and present the SK coefficients in a nonorthogonal basis set. An anisotropic Dirac cone is found in the band structure of borophane. We derive a Dirac low-energy Hamiltonian and compare the Fermi velocities with that of graphene.

Abstract: We present the Tight-Binding Studio (TB Studio) software package that calculates the different parameters of a tight-binding Hamiltonian from a set of Bloch energy bands obtained from first principle theories such as density functional theory, Hartree-Fock calculations or semi-empirical band-structure theory. This will be helpful for scientists who are interested in studying electronic and optical properties of structures using Green's function theory within the tight-binding approximation. TB Studio is a cross-platform application written in C++ with a graphical user interface design that is user-friendly and easy to work with. This software is powered by Linear Algebra Package C interface library for solving the eigenvalue problems and the standard high performance OpenGL graphic library for real time plotting. TB Studio and its examples together with the tutorials are available for download from tight-binding.com. Program summary Program Title: Tight-Binding Studio Program Files doi:http://dx.doi.org/10.17632/j6x5mwzm2d.1 Licensing provisions: LGPL Programming language: C++ External routines: BLAS, LAPACK, LAPACKE, wxWidgets, OpenGL, MathGL Nature of problem: Obtaining Tight-Binding Hamiltonian from a set of Bloch energy bands obtained from first-principles calculations. Solution method: Starting from the simplified LCAO method, a tight-binding model in the two-center approximation is constructed. The Slater and Koster (SK) approach is used to calculate the parameters of the TB Hamiltonian. By using non-linear fitting approaches the optimal values of the SK parameters are obtained such that the TB energy eigenvalues are as close as possible to those from first-principles calculations. We obtain the expression for the Hamiltonian and the overlap matrix elements between the different orbitals of the different atoms in an orthogonal or non-orthogonal basis set. (C) 2020 Elsevier B.V. All rights reserved.