On the forbidden gap of finite graphene nanoribbons

The electronic structure of isolated finite graphene nanoribbons is investigated by solving, at the Hartree-Fock (HF) level, the Pariser, Parr and Pople (PPP) many-body Hamiltonian. The study is mainly focused on 7-AGNR and 13-AGNR (Armchair Graphene Nano-Ribbons), whose electronic structures have been recently experimentally investigated. Only paramagnetic solutions are considered. The characteristics of the forbidden gap are studied as a function of the ribbon length. For a 7-AGNR, the gap monotonically decreases from a maximum value of ~6.5 eV for short nanoribbons to a very small value of ~0.12 eV for the longer calculated systems. Gap edges are defined by molecular orbitals that are spatially localized near the nanoribbon extremes, that is, near both zig-zag edges. On the other hand, two delocalized orbitals define a much larger gap of about 5 eV. Conductance measurements report a somewhat smaller gap of ~3 eV. The small real gap lies in the middle of the one given by extended states and has been observed by STM and reproduced by DFT calculations. On the other hand, the length dependence of the gap is not monotonous for a 13-AGNR. It decreases initially but sharply increases for lengths beyond 30 Å remaining almost constant thereafter at a value of ~2.1 eV. Two additional states localized at the nanoribbon extremes show up at energies 0.31 eV below the HOMO (Highest Occupied Molecular Orbital) and above the LUMO (Lowest Unoccupied Molecular Orbital). These numbers compare favorably with those recently obtained by means of STS for a 13-AGNR sustained by a gold surface, namely 1.4 eV for the energy gap and 0.4 eV for the position of localized band edges. We show that the important differences between 7- and 13-AGNR should be ascribed to the charge rearrangement near the zig-zag edges obtained in our calculations for ribbons longer than 30 Å, a feature that does not show up for a 7-AGNR no matter its length.

Quantum theory of spin waves in finite chiral spin chains

We calculate the effect of spin waves on the properties of finite-size spin chains with a chiral spin ground state observed on biatomic Fe chains deposited on iridium(001). The system is described with a Heisenberg model supplemented with a Dzyaloshinskii-Moriya coupling and a uniaxial single ion anisotropy that presents a chiral spin ground state. Spin waves are studied using the Holstein-Primakoff boson representation of spin operators. Both the renormalized ground state and the elementary excitations are found by means of Bogoliubov transformation, as a function of the two variables that can be controlled experimentally, the applied magnetic field and the chain length. Three main results are found. First, because of the noncollinear nature of the classical ground state, there is a significant zero-point reduction of the ground-state magnetization of the spin spiral. Second, there is a critical external field from which the ground state changes from chiral spin ground state to collinear ferromagnetic order. The character of the two lowest-energy spin waves changes from edge modes to confined bulk modes over this critical field. Third, in the spin-spiral state, the spin-wave spectrum exhibits oscillatory behavior as function of the chain length with the same period of the spin helix.

A Finite Element Numerical Algorithm for Modelling and Data Fitting in Complex Systems

Numerical modelling methodologies are important by their application to engineering and scientific problems, because there are processes where analytical mathematical expressions cannot be obtained to model them. When the only available information is a set of experimental values for the variables that determine the state of the system, the modelling problem is equivalent to determining the hyper-surface that best fits the data. This paper presents a methodology based on the Galerkin formulation of the finite elements method to obtain representations of relationships that are defined a priori, between a set of variables: y = z(x1, x2,...., xd). These representations are generated from the values of the variables in the experimental data. The approximation, piecewise, is an element of a Sobolev space and has derivatives defined in a general sense into this space. The using of this approach results in the need of inverting a linear system with a structure that allows a fast solver algorithm. The algorithm can be used in a variety of fields, being a multidisciplinary tool. The validity of the methodology is studied considering two real applications: a problem in hydrodynamics and a problem of engineering related to fluids, heat and transport in an energy generation plant. Also a test of the predictive capacity of the methodology is performed using a cross-validation method.