Abstract: We present results of Monte Carlo simulations of chains of C-60 and chains of C-70 molecules encapsulated in a single-walled carbon nanotube (SWCNT). We observe the changes in the configuration of the fullerene molecules when varying tube radius and temperature. In particular, the evolution of the pair correlation functions reveal a transition from linear harmonic chain behavior to a hard-sphere liquid upon heating, demonstrating the possibility of tuning properties of C-60- and C-70@SWCNT peapods with radius and temperature.

Abstract: We performed a detailed investigation of the oxidation of the technologically relevant GaAs(001)-beta 2(2x4) surface via density functional calculations. The purpose is to gain insights on the atomistic mechanisms and local bondings that underlie the degradation of the surface properties once exposed to oxygen. The study comprises the adsorption of single O atoms, through the sampling of several adsorption sites, and the subsequent formation of the O adsorbate at increasing coverage by taking into account multiple-atom adsorption. Based on the evaluation of the energetics and the structural properties of the atomistic models generated, the results here reported delineate a consistent picture of the initial stage of the surface oxidation: (i) at low coverage, in the limit of single O insertions, oxygen is incorporated on the surface forming a twofold-bridging Ga-O-As bond; (ii) at increasing coverage, as multiple O atoms are involved, this is accompanied by the formation of a threefold-coordinated bond (with two Ga and one As atoms); (iii) the latter has important implications regarding the electronic properties of the adsorbate since this O bonding may result in the formation of As dangling bonds. Moreover, a clear trend of increased energy gain for the incorporation of neighboring O atoms compared to single O insertions indicates that the formation of oxide clusters is favored over a regime of uniform oxidation. Our findings provide a detailed description of the O bonding and stress the importance of modeling the adsorption of multiple O atoms for an accurate description of the surface oxidation.

Abstract: The treatment of sp2 carbon materials, including micrographite, nanographite, HOPG, onion-like-carbon, and single-walled carbon nanotubes, in a 3:1 sulfuric to nitric acid mixture produced photoluminescent reaction solutions. These colloidal, aqueous solutions appeared photoluminescently stable under a UV lamp and ranged in color from red to blue. The photoluminescent wavelength shifted to shorter wavelength with increasing reaction time or increasing reaction temperature. Raman spectroscopy showed evidence of defect structures in graphitic residue, and transmission electron microscopy showed unusual structures present in the supernatant including graphitic balls.

Abstract: Spin rotations induced by an ac electric field in a two-electron double quantum dot are studied by an exact numerical solution of the time-dependent Schrodinger equation in the context of recent electric-dipole spin resonance experiments on gated nanowires. We demonstrate that the splitting of the main resonance line by the spin exchange coupling is accompanied by the appearance of fractional resonances and that both these effects are triggered by interdot tunnel coupling. We find that the ac-driven system generates residual but distinct harmonics of the driving frequency, which are amplified when tuned to the main transition frequency. The mechanism is universal for electron systems in electrically driven potentials and works also in the absence of electron-electron interaction or spin-orbit coupling.

Abstract: RuCp2 (ruthenocene) and RuCpPy (cyclopentadienyl pyrrolyl ruthenium) complexes are used in ruthenium (Ru) atomic layer deposition (ALD) but exhibit a markedly different reactivity with respect to the substrate and co-reactant. In search of an explanation, we report here the results of a comparative study of the heterolytic and homolytic dissociation enthalpy of these two ruthenium complexes, making use of either density functional theory (DFT) or multiconfigurational perturbation theory (CASPT2). While both methods predict distinctly different absolute dissociation enthalpies, they agree on the relative values between both molecules. A reduced heterolytic dissociation enthalpy is obtained for RuCpPy compared to RuCp2, although the difference obtained from CASPT2 (19.9 kcal/mol) is slightly larger than the one obtained with any of the DFT functionals (around 17 kcal/mol). Both methods also agree on the more pronounced stability of the Cp- ligand in RuCpPy than in RuCp2 (by around 9 kcal/mol with DFT and by 6 kcal/mol with CASPT2).

Abstract: We present a comprehensive study of the crystal structure, magnetic structure, and microscopic magnetic model of (CuBr)LaNb2O7, the Br analog of the spin-gap quantum magnet (CuCl) LaNb2O7. Despite similar crystal structures and spin lattices, the magnetic behavior and even peculiarities of the atomic arrangement in the Cl and Br compounds are very different. The high- resolution x-ray and neutron data reveal a split position of Br atoms in (CuBr) LaNb2O7. This splitting originates from two possible configurations developed by [CuBr] zigzag ribbons. While the Br atoms are locally ordered in the ab plane, their arrangement along the c direction remains partially disordered. The predominant and energetically more favorable configuration features an additional doubling of the c lattice parameter that was not observed in (CuCl) LaNb2O7. (CuBr) LaNb2O7 undergoes long-range antiferromagnetic ordering at T-N = 32 K, which is nearly 70% of the leading exchange coupling J4 similar or equal to 48 K. The Br compound does not show any experimental signatures of low-dimensional magnetism because the underlying spin lattice is three-dimensional. The coupling along the c direction is comparable to the couplings in the ab plane, even though the shortest Cu-Cu distance along c (11.69 angstrom) is three times larger than nearest-neighbor distances in the ab plane (3.55 angstrom). The stripe antiferromagnetic long-range order featuring columns of parallel spins in the ab plane and antiparallel spins along c is verified experimentally and confirmed by the microscopic analysis.

Abstract: Iron carbides play a crucial role in steel manufacturing and processing and to a large extent determine the physical properties of steel products. The modified embedded atom method (MEAM) in combination with Lee's Fe-C potential is a good candidate for molecular dynamics simulations on larger Fe-C systems. Here, we investigate the stability and crystal structures of pure iron and binary iron carbides using MEAM and compare them with the experimental data and quantum-mechanical density functional theory calculations. The analysis shows that the Fe-C potential gives reasonable results for the relative stability of iron and iron carbides. The performance of MEAM for the prediction of the potential energy and the calculated lattice parameters at elevated temperature for pure iron phases and cementite are investigated as well. The conclusion is that Lee's MEAM Fe-C potential provides a promising basis for further molecular dynamics simulations of Fe-C alloys and steels at lower temperatures (up to 800 K).

Abstract: Using an effective one-dimensional cluster expansion in combination with first-principles electronic structure calculations we have studied the energetics and electronic properties of Sb-Te layered systems. For a Te concentration between 0 and 60 at. % an almost continuous series of metastable structures is obtained consisting of consecutive Sb bilayers next to consecutive Sb2Te3 units, with the general formula (Sb-2)(n)(Sb2Te3)(m) (n, m = 1,2, ... ). Between 60 and 100 at.% no stable structures are found. We account explicitly for the weak van derWaals bonding between Sb bilayers and Sb2Te3 units by using a recently developed functional, which strongly improves the interlayer bonding distances. At T = 0 K, no evidence is found for the existence of two separate single-phase regions delta and gamma and a two-phase region delta + gamma. Metastable compounds with a Te concentration between 0 and 40 at. % are semimetallic, whereas compounds with a Te concentration between 50 and 60 at. % are semiconducting. Compounds with an odd number of Sb layers are metallic and have a much higher formation energy than those with an even number of consecutive Sb layers, thereby favoring the formation of Sb bilayers.

Abstract: The strain-induced pseudomagnetic field in supported graphene deposited on top of a nanostructured substrate is investigated by using atomistic simulations. A step, an elongated trench, a one-dimensional barrier, a spherical bubble, a Gaussian bump, and a Gaussian depression are considered as support structures for graphene. From the obtained optimum configurations we found very strong induced pseudomagnetic fields which can reach up to similar to 1000 T due to the strain-induced deformations in the supported graphene. Different magnetic confinements with controllable geometries are found by tuning the pattern of the substrate. The resulting induced magnetic fields for graphene on top of a step, barrier, and trench are calculated. In contrast to the step and trench the middle part of graphene on top of a barrier has zero pseudomagnetic field. This study provides a theoretical background for designing magnetic structures in graphene by nanostructuring substrates. We found that altering the radial symmetry of the deformation changes the sixfold symmetry of the induced pseudomagnetic field.

Abstract: The electronic properties of two-dimensional honeycomb structures of molybdenum disulfide (MoS(2)) subjected to biaxial strain have been investigated using first-principles calculations based on density functional theory. On applying compressive or tensile bi-axial strain on bi-layer and mono-layer MoS(2), the electronic properties are predicted to change from semiconducting to metallic. These changes present very interesting possibilities for engineering the electronic properties of two-dimensional structures of MoS(2).

Abstract: Unidirectional chiral states are predicted in single layer graphene which originate from the breaking of the sublattice symmetry due to an asymmetric mass potential. The latter can be created experimentally using boron-nitride (BN) substrates with a line defect (B-B or N-N) that changes the induced mass potential in graphene. Solving the Dirac-Weyl equation, the obtained energy spectrum is compared with the one calculated using ab initio density functional calculations. We found that these one-dimensional chiral states are very robust and they can even exist in the presence of a small gap between the mass regions. In the latter case additional bound states are found that are topologically different from those chiral states.

Abstract: We apply neutron diffraction, high-resolution synchrotron x-ray diffraction, magnetization measurements, electronic structure calculations, and quantum Monte-Carlo simulations to unravel the structure and magnetism of (CuCl)LaTa2O7. Despite the pseudo-tetragonal crystallographic unit cell, this compound features an orthorhombic superstructure, similar to the Nb-containing (CuX)LaNb2O7 with X = Cl and Br. The spin lattice entails dimers formed by the antiferromagnetic fourth-neighbor coupling J(4), as well as a large number of nonequivalent interdimer couplings quantified by an effective exchange parameter J(eff). In (CuCl)LaTa2O7, the interdimer couplings are sufficiently strong to induce the long-range magnetic order with the Neel temperature T-N similar or equal to 7 K and the ordered magnetic moment of 0.53 mu(B), as measured with neutron diffraction. This magnetic behavior can be accounted for by J(eff)/J(4) similar or equal to 1.6 and J(4) similar or equal to 16 K. We further propose a general magnetic phase diagram for the (CuCl)LaNb2O7-type compounds, and explain the transition from the gapped spin-singlet (dimer) ground state in (CuCl)LaNb2O7 to the long-range antiferromagnetic order in (CuCl)LaTa2O7 and (CuBr)LaNb2O7 by an increase in the magnitude of the interdimer couplings J(eff)/J(4), with the (CuCl)LaM2O7 (M = Nb, Ta) compounds lying on different sides of the quantum critical point that separates the singlet and long-range-ordered magnetic ground states.

Abstract: Using highly efficient simulations of the tight-binding Bogoliubov-de-Gennes model, we solved self-consistently for the pair correlation and the Josephson current in a superconducting-bilayer graphene-superconducting Josephson junction. Different doping levels for the non-superconducting link are considered in the short- and long-junction regimes. Self-consistent results for the pair correlation and superconducting current resemble those reported previously for single-layer graphene except at the Dirac point, where remarkable differences in the proximity effect are found, as well as a suppression of the superconducting current in the long-junction regime. Inversion symmetry is broken by considering a potential difference between the layers and we found that the supercurrent can be switched if the junction length is larger than the Fermi length.

Abstract: We investigate the optimal orientations and positions of tubular fullerene molecules C-70, C-80 and C-90 encapsulated in single-walled carbon nanotubes (SWCNTs). We find that increasing the tube radius leads to the following succession of energetically stable regimes: 1) lying molecules positioned on the tube's long axis, 2) tilted molecules on the tube's long axis and 3) lying molecules shifted away from the tube's long axis. In the case of C-70 and C-80 molecules, standing on-axis configurations also occur. Our findings are relevant for the possible application of molecular-orientation-dependent electronic properties of fullerene nanopeapods.

Abstract: Monte Carlo (MC) methods have a long-standing history as partners of molecular dynamics (MD) to simulate the evolution of materials at the atomic scale. Among these techniques, the uniform-acceptance force-bias Monte Carlo (UFMC) method [ G. Dereli Mol. Simul. 8 351 (1992)] has recently attracted attention [ M. Timonova et al. Phys. Rev. B 81 144107 (2010)] thanks to its apparent capacity of being able to simulate physical processes in a reduced number of iterations compared to classical MD methods. The origin of this efficiency remains, however, unclear. In this work we derive a UFMC method starting from basic thermodynamic principles, which leads to an intuitive and unambiguous formalism. The approach includes a statistically relevant time step per Monte Carlo iteration, showing a significant speed-up compared to MD simulations. This time-stamped force-bias Monte Carlo (tfMC) formalism is tested on both simple one-dimensional and three-dimensional systems. Both test-cases give excellent results in agreement with analytical solutions and literature reports. The inclusion of a time scale, the simplicity of the method, and the enhancement of the time step compared to classical MD methods make this method very appealing for studying the dynamics of many-particle systems.