Tunable gap in stable arsenene nanoribbons opens the door to electronic applications.

Arsenic has been predicted to present significantly more diverse 2D phases than other elemental compounds like graphene. While practical applications must be based on finite arsenene samples, like nanoribbons, theory has so far focused on the infinite sheet. Our ab initio simulations show the clear contrast between the properties of arsenene nanoribbons and those of the monolayer, ranging from phase stability to electronic structure. We include nanoribbons derived from the buckled, puckered and square/octagon structures of bulk arsenene. The flexibility afforded by different parent structures, widths and edge passivations leads to a rich variety of semiconducting structures with tunable gaps.

What will freestanding borophene nanoribbons look like? An analysis of their possible structures, magnetism and transport properties.

We report a density-functional theory study of the stability and electronic structure of two recently proposed borophene sheets with Pmmn symmetry and nonzero thickness. We then investigate nanoribbons (BNRs) derived from these nanostructures, with particular attention to technologically relevant properties like magnetism and electronic transport. We consider two perpendicular directions for the edges of the stripes as well as different lateral widths. We show that the Pmmn8 sheet, with 8 atoms in its unit cell and generated by two interpenetrating lattices, has a larger binding energy than the Pmmn2 sheet, with only 2 atoms per unit cell. We also use their phonon spectra to show that the mechanical stability of the Pmmn8 sheet is superior to that of the Pmmn2 sheet. Nanoribbons derived from Pmmn8 are not only more stable than those derived from Pmmn2, but also more interesting from the technological point of view. We find a rich variety of magnetic solutions, depending on the borophene "mother structure", edge orientation, width and, in the case of Pmmn8-derived BNRs, the sublattice of edge atoms. We show that one can build BNRs with magnetic moment in both, one or none of the edges, as well as with parallel or antiparallel magnetic coupling between the edges when magnetic; moreover, their electronic character can be semiconducting, metallic or half-metallic, creating a perfect spin valve at low bias. These different behaviors are reflected in their densities of states, spin density and electronic transport coefficients, which are analyzed in detail. Our work provides a complete overview of what one may expect if nanoribbons are cut out from Pmmn sheets with a view to potential technological applications.

Molecular dynamics simulation of the behaviour of water in nano-confined ionic liquid-water mixtures.

This work describes the behaviour of water molecules in 1-butyl-3-methylimidazolium tetrafluoroborate ionic liquid under nanoconfinement, between graphene sheets. By means of molecular dynamics simulations, the adsorption of water molecules at the graphene surface is studied. A depletion of water molecules in the vicinity of the neutral and negatively charged graphene surfaces, and their adsorption at the positively charged surface are observed in line with the preferential hydration of the ionic liquid anions. The findings are appropriately described using a two-level statistical model. The confinement effect on the structure and dynamics of the mixtures is thoroughly analyzed using the density and the potential of mean force profiles, as well as by the vibrational densities of the states of water molecules near the graphene surface. The orientation of water molecules and the water-induced structural transitions in the layer closest to the graphene surface are also discussed.

Spin-polarized transport in hydrogen-passivated graphene and silicene nanoribbons with magnetic transition-metal substituents.

We present a systematic theoretical study of the electronic transport in hydrogen passivated zigzag graphene and silicene nanoribbons with between zero and four neighboring H atoms on one edge replaced by magnetic transition metals (Fe, Co, and Ni). The calculations were performed using equilibrium transport and density-functional theory with the generalized gradient approximation to exchange and correlation. We considered the magnetic moments of the two edges aligned both ferromagnetically (Ferro-F form) and antiferromagnetically (Ferro-A form). The Ferro-A graphene-based ribbons were all semiconducting and would support moderate spin-polarized currents of either sign by applying positive or negative gate voltages. The Ferro-F graphene-based ribbons were all metallic; the most interesting for possible spintronic applications being that with a single Ni atom, in which strong spin-filtering at low bias resulted from a deep trough in the transmission of one spin component around the Fermi level. By contrast, in the Si-based analog this trough was split, partially eliminating the polarization of the current. This splitting was found to be related to the buckled structure of the Si-based nanoribbon, which has its origin in its preference for sp(3)-like hybridization.

Structural and electronic properties of zigzag InP nanoribbons with Stone-Wales type defects.

By means of density-functional-theoretic calculations, we investigate the structural and electronic properties of a hexagonal InP sheet and of hydrogen-passivated zigzag InP nanoribbons (ZInPNRs) with Stone-Wales (SW)-type defects. Our results show that the influence of this kind of defect is not limited to the defected region but it leads to the formation of ripples that extend across the systems, in keeping with the results obtained recently for graphene and silicene sheets. The presence of SW defects in ZInPNRs causes an appreciable broadening of the band gap and transforms the indirect-bandgap perfect ZInPNR into a direct-bandgap semiconductor. An external transverse electric field, regardless of its direction, reduces the gap in both the perfect and defective ZInPNRs.

Spin currents and filtering behavior in zigzag graphene nanoribbons with adsorbed molybdenum chains.

By means of density-functional-theoretic calculations, we investigated the structural, electronic and transport properties of hydrogen-passivated zigzag graphene nanoribbons (ZGNRs) on which a one-atom-thick Mo chain was adsorbed (with or without one or two missing atoms), or in which the passivating hydrogen atoms were replaced by Mo atoms. Mo-passivated ZGNRs proved to be nonmagnetic. ZGNRs with an adsorbed defect-free Mo chain were most stable with the Mo atoms forming dimers above edge bay sites, which suppressed the magnetic moments of the C atoms in that half of the ribbon; around the Fermi level of these systems, each spin component had a transmission channel via the Mo spz band and one had an additional channel created by polarization of the ZGNR π* band, leading to a net spin current. The absence of an Mo dimer from an Mo chain adsorbed at the ZGNR edge made the system a perfect spin filter at low voltage bias by suppressing the Mo spz band channels. Thus this last kind of hybrid system is a potential spin valve.

Spin-dependent electronic conduction along zigzag graphene nanoribbons bearing adsorbed Ni and Fe nanostructures.

Using SMEAGOL, an ab initio computational method that combines the non-equilibrium Green's function formalism with density-functional theory, we calculated spin-specific electronic conduction in systems consisting of single Fen and Nin nanostructures (n = 1-4) adsorbed on a hydrogen-passivated zigzag graphene nanoribbon. For each cluster we considered both ferromagnetically and antiferromagnetically coupled ribbon edges (Ferro-F and Ferro-A systems, respectively). Adstructures located laterally on Ferro-A ribbons caused significant transmittance loss at energies 0.6-0.25 eV below the Fermi level for one spin and 0.2-0.4 eV above the Fermi level for the other, allowing the potential use of these systems in transistors to create a moderately spin-polarized current of one or the other sign depending on the gate voltage. Ni3 and Ni4 clusters located at the centre of Ferro-F ribbons exhibited a strong spin-filtering effect in a narrow energy window around the Fermi level.

Electronic properties of pure and p-type doped hexagonal sheets and zigzag nanoribbons of InP.

Unlike graphene, a hexagonal InP sheet (HInPS) cannot be obtained by mechanical exfoliation from the native bulk InP, which crystallizes in the zinc blende structure under ambient conditions. However, by ab initio density functional theory calculations we found that a slightly buckled HInPS is stable both in pristine form and when doped with Zn atoms; the same occurred for hydrogen-passivated zigzag InP nanoribbons (ZInPNRs), quasi-one-dimensional versions of the quasi-two-dimensional material. We investigated the electronic properties of both nanostructures, in the latter case also in the presence of an external transverse electric field, and the results are compared with those of hypothetical planar HInPS and ZInPNRs. The band gaps of planar ZInPNRs were found to be tunable by the choice of strength of this field, and to show an asymmetric behavior under weak electric fields, by which the gap can either be increased or decreased depending on their direction; however, this effect is absent from slightly buckled ZInPNRs. The binding energies of the acceptor impurity states of Zn-doped HInPS and ZInPNRs were found to be similar and much larger than that of Zn-doped bulk InP. These latter findings show that the reduction of the dimensionality of these materials limits the presence of free carriers.

Electronic structure and transport properties of monatomic Fe chains in a vacuum and anchored to a graphene nanoribbon.

The electronic structure and transport properties of monatomic Fe wires of different characteristics are studied within the density functional theory. In both equidistant and dimerized (more stable) isolated wires, magnetism plays an important role since it leads to different shapes of the transmission coefficients for each spin component. In equidistant wires, electron localization around the Fermi level leads to symmetry breaking between d(xy) and d(x(2)-y(2)) bands. The main effect of the structural dimerization is to decrease the number of channels available for the minority spin component. When anchored to the edges of a graphene nanoribbon, the dimerization of the chain is preserved, despite the hybridization of the d states of Fe with the C atoms which gives way to a reduction in the number of d channels around the Fermi level. Most conduction is then led by an electronic channel from the ribbon and the sp(z) bands from the Fe wires. Suggestions to improve the spintronic ability of Fe wires are proposed.

Magnetic cooperative effects in small Ni-Ru clusters.

We report ab initio calculations of the structures, magnetic moments, and electronic properties of Ni(7)-xRu(x) clusters (x = 0-7) using a density-functional method that employs linear combinations of pseudoatomic orbitals as basis sets, nonlocal norm-conserving pseudopotentials, and the generalized gradient approximation to exchange and correlation. The pure clusters Ni(7) and Ru(7) were predicted to have octahedral and cubic structures, respectively, and the binary clusters were found to be either decahedral (Ni(6)Ru, Ni(5)Ru(2), and Ni(4)Ru(3)) or cubic (Ni(3)Ru(4), Ni(2)Ru(5), and NiRu(6)). For Ni(5)Ru(2) and Ni(4)Ru(3) we found a magnetic cooperative phenomenon, which is due to both geometrical effects and electronic contributions through Ni-Ru hybridization.

Molecular dynamics simulations of the structural and thermodynamic properties of imidazolium-based ionic liquid mixtures.

In this work, extensive molecular dynamics simulations of mixtures of alcohols of several chain lengths (methanol and ethanol) with the ionic liquids (ILs) composed of the cation 1-hexyl-3-methylimidazolium and several anions of different hydrophobicity degrees (Cl(-), BF(4)(-), PF(6)(-)) are reported. We analyze the influence of the nature of the anion, the length of the molecular chain of the alcohol, and the alcohol concentration on the thermodynamic and structural properties of the mixtures. Densities, excess molar volumes, total and partial radial distribution functions, coordination numbers, and hydrogen bond degrees are reported and analyzed for mixtures of the ILs with methanol and ethanol. The aggregation process is shown to be highly dependent on the nature of the anion and the size of the alcohol, since alcohol molecules tend to interact predominantly with the anionic part of the IL, especially in mixtures of the halogenated IL with methanol. Particularly, our results suggest that the formation of an apolar network similar to that previously reported in mixtures of ILs with water does not take place in mixtures with alcohol when the chloride anion is present, the alcohol molecules being instead homogeneously distributed in the polar network of IL. Moreover, the alcohol clusters formed in mixtures of [HMIM][PF(6)] with alcohol were found to have a smaller size than in mixtures with water. Additionally, we provide a semiquantitative analysis of the dependence of the hydrogen bonding degree of the mixtures on the alcohol concentration.

Prediction of phonon thermal transport in thin GaAs, InAs and InP nanowires by molecular dynamics simulations: influence of the interatomic potential.

A number of different potentials are currently being used in molecular dynamics simulations of semiconductor nanostructures. Confusion can arise if an inappropriate potential is used. To illustrate this point, we performed direct molecular dynamics simulations to predict the room temperature lattice thermal conductivity λ of thin GaAs, InAs and InP nanowires. In each case, simulations performed using the classical Harrison potential afforded values of λ about an order of magnitude smaller than those obtained using more elaborate potentials (an Abell-Tersoff, as parameterized by Hammerschmidt et al for GaAs and InAs, and a potential of Vashishta type for InP). These results will be a warning to those wishing to use computer simulations to orient the development of quasi-one-dimensional systems as heat sinks or thermoelectric devices.

Magnetism of substitutional Fe impurities in graphene nanoribbons.

Using the generalized gradient approximation to exchange and correlation, we perform density functional calculations on an Fe atom at a single vacancy of graphene nanoribbons. Our results show that, after relaxation, the Fe atom is magnetic, in contrast to the behavior recently found for Fe at a single vacancy of the graphene sheet.

A density-functional study of the vertical ionization potentials of the cluster Mn13.

Using the density-functional method SIESTA, nonlocal norm-conserving pseudopotentials and the generalized gradient approximation to exchange and correlation, we show that the two ionization energies of the photoionization spectrum of Mn(13) can be attributed to the presence of the ground state and an excited spin isomer of this species.

A density-functional study of the structures and electronic properties of neutral, anionic, and endohedrally doped In(x)P(x) clusters.

We report extensive ab initio calculations of the structures, binding energies, and magnetic moments of In(x)P(x) and In(x)P(x) (-) clusters (x=1-15) using a density-functional method that employs linear combinations of pseudoatomic orbitals as basis sets, nonlocal norm-conserving pseudopotentials, and the generalized gradient approximation for exchange and correlation. Our results, which are compared with those obtained previously for some of these clusters by means of all-electron calculations, show that hollow cages with alternating In-P bonds are energetically preferred over other structures for both the neutral and anionic species within the range x=6-15. We also consider the endohedrally doped X@In(10)P(10) (X=Cr,Mn,Fe,Co) and Ti@In(x)P(x) (x=7-12) clusters. Our results show that, except for Ti@In(7)P(7) and Ti@In(8)P(8), the transition metal atoms preserve their atomic spin magnetic moments when encapsulated in the InP cages, instead of suffering either a spin crossover or a spin quenching due to hybridization effects. We also show that the stabilities of some empty and doped InP cages can be explained on the basis of the jellium model.

Nonequilibrium nanothermodynamics.

Entropy production for a system outside the thermodynamic limit is formulated using Hill's nanothermodynamics, in which a macroscopic ensemble of such systems is considered. The external influence of the environment on the average nanosystem is connected to irreversible work with an explicit formula based on the Jarzynski equality. The entropy production retains its usual form as a sum of products of fluxes and forces and Onsager's symmetry principle is proven to hold for the average nanosystem, if it is assumed to be valid for the macroscopic ensemble, by two methods. The first one provides expressions that relate the coefficients of the two systems. The second gives a general condition for a system under an external force to preserve Onsager's symmetry.

A density-functional study of the possibility of noncollinear magnetism in small Mn clusters using SIESTA and the generalized gradient approximation to exchange and correlation.

We investigated the possibility of noncollinear magnetism in small Mn(n) clusters (n=2-6) using the density-functional method SIESTA with the generalized gradient approximation (GGA) to exchange and correlation. The lowest-energy states identified were collinear, with the atomic spin magnetic moments pointing in the same direction, for Mn(2) and Mn(3), and noncollinear for Mn(4), Mn(5) and, most decidedly, Mn(6). These SIESTA/GGA results, which are compared with those of an earlier SIESTA study that used the local spin density approximation, are qualitatively in keeping with the result obtained by VASP/GGA calculations.

A density-functional study of the structures, binding energies and magnetic moments of the clusters Mo(N) (N = 2-13), Mo(12)Fe, Mo(12)Co and Mo(12)Ni.

We report the results of density-functional calculations of the structures, binding energies and magnetic moments of the clusters Mo(N) (N = 2-13), Mo(12)Fe, Mo(12)Co and Mo(12)Ni that were performed using the SIESTA method within the generalized gradient approximation for exchange and correlation. For pure Mo(N) clusters, we obtain collinear magnetic structures in all cases, even when the self-consistent calculations were started from non-collinear inputs. Our results for these clusters show that both linear, planar and three-dimensional clusters have a strong tendency to form dimers. In general, even-numbered clusters are more stable than their neighbouring odd-numbered clusters because they can accommodate an integer number of tightly bound dimers. As a consequence, the binding energies of pure Mo(N) clusters, in their lowest-energy states, exhibit an odd-even effect in all dimensionalities. Odd-even effects are less noticeable in the magnetic moments than in the binding energies. When comparing our results for pure Mo clusters with those obtained recently by other authors, we observe similarities in some cases, but striking differences in others. In particular, the odd-even effect in three-dimensional Mo clusters was not observed before, and our results for some clusters (e.g. for planar Mo(3) and Mo(7) and for three-dimensional Mo(7) and Mo(13)) differ from those reported by other authors. For Mo(12)Fe and Mo(12)Ni, we obtain that the icosahedral configuration with the impurity atom at the cluster surface is more stable than the configuration with the impurity at the central site, while the opposite occurs in the case of Mo(12)Co. In Mo(12)Co and Mo(12)Ni, the impurities exhibit a weak magnetic moment parallely coupled to the total magnetic moment of the Mo atoms, whereas in Mo(12)Fe the impurity shows a high moment with antiparallel coupling.

Engineering the magnetic structure of Fe clusters by Mn alloying.

We propose to tailor the magnetic structure of atomic clusters by suitable doping, which produces the nanometric equivalent to alloying. As a proof of principle, we perform a theoretical analysis of Fe(6-x)Mn(x) clusters (x = 0-5), which shows a modulation of the magnetic moment of the clusters as a function of Mn doping and, more importantly, a collinear to noncollinear transition at x = 4.