A Generalized Nomenclature Scheme for Graphene Pores, Flakes, and Edges, and an Algorithm for Their Generation and Numbering.

In the present study, we generalize our recently proposed nomenclature scheme for porous graphene structures to include graphene flakes and (periodic) edges, i.e., nanographenes and graphene nanoribbons. The proposed nomenclature scheme is a complete scheme that similarly treats all these structures. Beyond this generalization, we study the geometric features of graphene flakes and edges based on ideas from the graph theory, as well as the pore-flake duality. Based on this study, we propose an algorithm for the systematic generation, identification, and numbering of graphene pores, flakes, and edges. The algorithm and the nomenclature scheme can also be used for flakes and edges of similar honeycomb systems.

Evaluating the performance of ReaxFF potentials for sp2 carbon systems (graphene, carbon nanotubes, fullerenes) and a new ReaxFF potential.

We study the performance of eleven reactive force fields (ReaxFF), which can be used to study sp2 carbon systems. Among them a new hybrid ReaxFF is proposed combining two others and introducing two different types of C atoms. The advantages of that potential are discussed. We analyze the behavior of ReaxFFs with respect to 1) the structural and mechanical properties of graphene, its response to strain and phonon dispersion relation; 2) the energetics of (n, 0) and (n, n) carbon nanotubes (CNTs), their mechanical properties and response to strain up to fracture; 3) the energetics of the icosahedral C60 fullerene and the 40 C40 fullerene isomers. Seven of them provide not very realistic predictions for graphene, which made us focusing on the remaining, which provide reasonable results for 1) the structure, energy and phonon band structure of graphene, 2) the energetics of CNTs versus their diameter and 3) the energy of C60 and the trend of the energy of the C40 fullerene isomers versus their pentagon adjacencies, in accordance with density functional theory (DFT) calculations and/or experimental data. Moreover, the predicted fracture strain, ultimate tensile strength and strain values of CNTs are inside the range of experimental values, although overestimated with respect to DFT. However, they underestimate the Young's modulus, overestimate the Poisson's ratio of both graphene and CNTs and they display anomalous behavior of the stress - strain and Poisson's ratio - strain curves, whose origin needs further investigation.

High temperature stability, metallic character and bonding of the Si2BN planar structure.

The family of monolayered Si2BN structures constitute a new class of 2D materials exhibiting metallic character with remarkable stability. Topologically, these structures are very similar to graphene, forming a slightly distorted honeycomb lattice generated by a union of two basic motifs with AA and AB stacking. In the present work we study in detail the structural and electronic properties of these structures in order to understand the factors which are responsible for their structural differences as well as those which are responsible for their metallic behavior and bonding. Their high temperature stability is demonstrated by the calculations of finite temperature phonon modes which show no negative contributions up to and beyond 1000 K. Presence of the negative thermal expansion coefficient, a common feature of one-atom thick 2D structures, is also seen. Comparison of the two motifs reveal the main structural differences to be the differences in their bond angles, which are affected by the third nearest neighbor interactions ofcis-transtype. On the other hand, the electronic properties of these two structures are very similar, including the charge transfers occurring between orbitals and between atoms. Their metallicity is mainly due to thepzorbitals of Si with a minor contribution from thepzorbitals of B, while the contribution from thepzorbitals of N atoms is negligible. There is almost no contributions from the Npzelectrons to the energy states near the Fermi level, and they form a band well below it. I.e., thepzelectrons of N are localized mostly at the N atoms and therefore cannot be considered as mobile electrons of thepzcloud. Moreover, we show that due to the relative positions in the energy axis of the atomic energies of thepzorbitals of B, N and Si atoms, the density of states (DOS) of Si2BN can be considered qualitatively as a combination of the DOS of planar hexagonal BN (h-BN) and hypothetically planar silicene (ph-Si). As a result, the Si2BN behaves electronically at the Fermi level as slightly perturbed ph-Si, having very similar electronic properties as silicene, but with the advantage of having kinetic stability in planar form. As for the bonding, the Si-Si bonds are covalent, while theπback donation mechanism occurs for the B-N bonding, in accordance with the B-N bonding in h-BN.

Atomistic potential for graphene and other sp2 carbon systems.

We introduce a torsional force field for sp2 carbon to augment an in-plane atomistic potential of a previous work [G. Kalosakas et al., J. Appl. Phys., 2013, 113, 134307] so that it is applicable to out-of-plane deformations of graphene and related carbon materials. The introduced force field is fit to reproduce density-functional-theory calculation data of appropriately chosen structures. The aim is to create a force field that is as simple as possible so it can be efficient for large scale atomistic simulations of various sp2 carbon structures without significant loss of accuracy. We show that the complete proposed potential reproduces characteristic properties of fullerenes and carbon nanotubes. In addition, it reproduces very accurately the out-of-plane acoustic and optical modes of graphene's phonon dispersion as well as all phonons with frequencies up to 1000 cm-1.

Structural deformations of two-dimensional planar structures under uniaxial strain: the case of graphene.

In the present work, a method for the study of the structural deformations of two dimensional planar structures under uniaxial strain is presented. The method is based on molecular mechanics using the original stick and spiral model and a modified one which includes second nearest neighbor interactions for bond stretching. As we show, the method allows an accurate prediction of the structural deformations of any two dimensional planar structure as a function of strain, along any strain direction in the elastic regime, if structural deformations are known along specific strain directions, which are used to calculate the stick and spiral model parameters. Our method can be generalized including other strain conditions and not only uniaxial strain. We apply this method to graphene and we test its validity, using results obtained from ab initio density functional theory calculations. What we find is that the original stick and spiral model is not appropriate to describe accurately the structural deformations of graphene in the elastic regime. However, the introduction of second nearest neighbor interactions provides a very accurate description.

Graphene allotropes under extreme uniaxial strain: an ab initio theoretical study.

Using density functional theory calculations, we study the response of three representative graphene allotropes (two pentaheptites and octagraphene) as well as graphene, to uniaxial strain up to their fracture limit. Those allotropes can be seen as distorted graphene structures formed upon periodically arranged Stone-Walles transformations. We calculate their mechanical properties (Young's modulus, Poisson's ratio, speed of sound, ultimate tensile strength and the corresponding strain), and we describe the pathways of their fracture. Finally, we study strain as a factor for the conversion of graphene into those allotropes upon Stone-Walles transformations. For specific sets of Stone-Walles transformations leading to an allotrope, we determine the strain directions and the corresponding minimum strain value, for which the allotrope is more favorable energetically than graphene. We find that the minimum strain values which favor those conversions are of the order of 9-13%. Moreover, we find that the energy barriers for the Stone-Walles transformations decrease dramatically under strain, however, they remain prohibitive for structural transitions. Thus, strain alone cannot provide a synthetic route to these allotropes, but could be a part of composite procedures for this purpose.

Successive spin polarizations underlying a new magnetic coupling contribution in diluted magnetic semiconductors.

We propose a new type of magnetic coupling (MC) that is found in diluted magnetic semiconductors (DMSs). The origin of this is found to be the result of charge transfer processes followed by successive spin polarizations (SSPs) along successive cation-anion segments which include the impurities. The basic process underlying the SSP-based MC (SSP-MC) is the sharing of a single spin orbital by two neighboring impurities. As such, it can be considered as a localized double exchange as it is not mediated by free carriers. SSP-MC can be either ferromagnetic (SSP-FMC) or antiferromagnetic (SSP-AFMC) and, as demonstrated here, the SSP-FMC can be significantly enhanced via codoping; it can act in competition with superexchange and/or double and/or p-d exchange interactions. While the SSP-MC is not directly related to the magnitude of the magnetic moments of the impurities, it depends strongly on the energy difference of the host and impurity d-band centers, the difference of their electronegativities and rather weakly on the coupling interactions between them as well as between the cations and their mediating anions. The validity of the proposed SSP-MC as a new type of magnetic coupling is demonstrated by ab initio results for DMSs, namely ZnO, GaN, GaP, TiO2 and MoS2 monodoped (with Co, Cu and Mn) and codoped (with Co-Cu-Co and Mn-Cu-Mn).

Topologically protected conduction state at carbon foam surfaces: an ab initio study.

We report results of ab initio electronic structure and quantum conductance calculations indicating the emergence of conduction at the surface of semiconducting carbon foams. The occurrence of new conduction states is intimately linked to the topology of the surface and not limited to foams of elemental carbon. Our interpretation based on rehybridization theory indicates that conduction in the foam derives from first- and second-neighbor interactions between p∥ orbitals lying in the surface plane, which are related to p⊥ orbitals of graphene. The topologically protected conducting state occurs on bare and hydrogen-terminated foam surfaces and is thus unrelated to dangling bonds. Our results for carbon foam indicate that the conductance behavior may be further significantly modified by surface patterning.

Uncovering the FUTREX-6100XL prediction equation for the percentage body fat.

Based on the near infra-red (NIR) interactance method, the FUTREX company has developed a series of instruments, for the estimation of the body fat percentage (%BF). %BF is estimated through prediction equations incorporated in the instruments, which for the newest models (FUTREX-6100XL and FUTREX-6100A/ZL) are proprietary and they are not published anywhere. This missing knowledge may lead to several misunderstandings and confusion and degrades those instruments to 'black boxes'. The present work uncovers and presents the prediction equation of FUTREX-6100/XL and discusses the contribution of each term of that equation to the %BF. Furthermore, this study presents the method used, which can be used to uncover equations incorporated in other instruments. This method is based on the idea of firstly uncovering the dependence of the equation on each parameter separately and then combining those dependencies to uncover the unknown equation.

Nanomechanical energy storage in twisted nanotube ropes.

We determine the deformation energetics and energy density of twisted carbon nanotubes and nanotube ropes that effectively constitute a torsional spring. Using ab initio and parametrized density functional calculations, we identify structural changes in these systems and determine their elastic limits. The deformation energy of twisted nanotube ropes contains contributions associated not only with twisting but also with stretching, bending, and compression of individual nanotubes. We quantify these energy contributions and show that their relative role changes with the number of nanotubes in the rope.

Applicability of the Hunjan-Ramaswamy global optimization method.

Int. J. Mol. Sci. 3, 30 (2002)]; Phys. Rev. E 66, 046704 (2002)]] a method for global optimization, according to which the global minimum of a potential V(f) can be found, if a potential V(i) (with a known global minimum) is transformed adiabatically in time to V(f) , with the use of a switching function of time g (t) , which interpolates between 0 and 1, and lies in the [0,1] interval. In the present work, the method is examined in detail. With the use of a very simple one-dimensional hypothetical potential, it is shown that the potential transformation may not always be followed by a global minimum transformation, which indicates that the method may not always be safely applied in determining the global minimum. An attempt to improve the method is made by allowing the switching function g (t) to take values outside the [0,1] interval. This improved method is shown to succeed in three different realistic problems, for which the original method fails.