Electonic transport properties of nitrate-doped carbon nanotube networks

The conductivity of carbon nanotube (CNT) networks can be improved markedly by doping with nitric acid. In the present work, CNTs and junctions of CNTs functionalized with NO$_3$ molecules are investigated to understand the microscopic mechanism of nitric acid doping. According to our density functional theory band structure calculations, there is charge transfer from the CNT to adsorbed molecules indicating p-type doping. The average doping efficiency of the NO$_3$ molecules is higher if the NO$_3$ molecules form complexes with water molecules. In addition to electron transport along individual CNTs, we have also studied electron transport between different types (metallic, semiconducting) of CNTs. Reflecting the differences in the electronic structures of semiconducting and metallic CNTs, we have found that besides turning semiconducting CNTs metallic, doping further increases electron transport most efficiently along semiconducting CNTs as well as through a junction between them.Comment: 13 pages, 12 figure

Equilibrium shape and dislocation nucleation in strained epitaxial nanoislands

We study numerically the equilibrium shapes, shape transitions and dislocation nucleation of small strained epitaxial islands with a two-dimensional atomistic model, using simple interatomic pair potentials. We first map out the phase diagram for the equilibrium island shapes as a function of island size (up to N = 105 atoms) and lattice misfit with the substrate and show that nanoscopic islands have four generic equilibrium shapes, in contrast with predictions from the continuum theory of elasticity. For increasing substrate-adsorbate attraction, we find islands that form on top of a finite wetting layer as observed in Stranski-Krastanow growth. We also investigate energy barriers and transition paths for transitions between different shapes of the islands and for dislocation nucleation in initially coherent islands. In particular, we find that dislocations nucleate spontaneously at the edges of the adsorbate-substrate interface above a critical size or lattice misfit.Comment: 4 pages, 3 figures, uses wrapfig.sty and epsfig.st

Stress release mechanisms for Cu on Pd(111) in the submonolayer and monolayer regimes

We study the strain relaxation mechanisms of Cu on Pd(111) up to the monolayer regime using two different computational methodologies, basin-hopping global optimization and energy minimization with a repulsive bias potential. Our numerical results are consistent with experimentally observed layer-by-layer growth mode. However, we find that the structure of the Cu layer is not fully pseudomorphic even at low coverages. Instead, the Cu adsorbates forms fcc and hcp stacking domains, separated by partial misfit dislocations. We also estimate the minimum energy path and energy barriers for transitions from the ideal epitaxial state to the fcc-hcp domain pattern.Comment: 4 pages, 4 figure

A chain-of-states acceleration method for the efficient location of minimum energy paths

We describe a robust and efficient chain-of-states method for computing Minimum Energy Paths (MEPs) associated to barrier-crossing events in poly-atomic systems, which we call the acceleration method. The path is parametrized in terms of a continuous variable t ∈ [0, 1] that plays the role of time. In contrast to previous chain-of-states algorithms such as the nudged elastic band or string methods, where the positions of the states in the chain are taken as variational parameters in the search for the MEP, our strategy is to formulate the problem in terms of the second derivatives of the coordinates with respect to t, i.e., the state accelerations. We show this to result in a very simple and efficient method for determining the MEP. We describe the application of the method to a series of test cases, including two low-dimensional problems and the Stone-Wales transformation in C60This work has been supported by the Spanish Research and Innovation Office through Project Nos. FIS2012-31713 and FIS2012-3754

Meteorological effects of the solar eclipse of 20 March 2015: analysis of UK Met Office automatic weather station data and comparison with automatic weather station data from the Faroes and Iceland.

Here, we analyse high-frequency (1 min) surface air temperature, mean sea-level pressure (MSLP), wind speed and direction and cloud-cover data acquired during the solar eclipse of 20 March 2015 from 76 UK Met Office weather stations, and compare the results with those from 30 weather stations in the Faroe Islands and 148 stations in Iceland. There was a statistically significant mean UK temperature drop of 0.83±0.63°C, which occurred over 39 min on average, and the minimum temperature lagged the peak of the eclipse by about 10 min. For a subset of 14 (16) relatively clear (cloudy) stations, the mean temperature drop was 0.91±0.78 (0.31±0.40)°C but the mean temperature drops for relatively calm and windy stations were almost identical. Mean wind speed dropped significantly by 9% on average during the first half of the eclipse. There was no discernible effect of the eclipse on the wind-direction or MSLP time series, and therefore we can discount any localized eclipse cyclone effect over Britain during this event. Similar changes in air temperature and wind speed are observed for Iceland, where conditions were generally clearer, but here too there was no evidence of an eclipse cyclone; in the Faroes, there was a much more muted meteorological signature.This article is part of the themed issue 'Atmospheric effects of solar eclipses stimulated by the 2015 UK eclipse'

Force-matched embedded-atom method potential for niobium

Large-scale simulations of plastic deformation and phase transformations in alloys require reliable classical interatomic potentials. We construct an embedded-atom method potential for niobium as the first step in alloy potential development. Optimization of the potential parameters to a well-converged set of density-functional theory (DFT) forces, energies, and stresses produces a reliable and transferable potential for molecular dynamics simulations. The potential accurately describes properties related to the fitting data, and also produces excellent results for quantities outside the fitting range. Structural and elastic properties, defect energetics, and thermal behavior compare well with DFT results and experimental data, e.g., DFT surface energies are reproduced with less than 4% error, generalized stacking-fault energies differ from DFT values by less than 15%, and the melting temperature is within 2% of the experimental value.Comment: 17 pages, 13 figures, 7 table

QMCube (QM3): An all‐purpose suite for multiscale QM/MM calculations

QMCube (QM3) is a suite written in the Python programming language, initially focused on multiscale QM/MM simulations of biological systems, but open enough to address other kinds of problems. It allows the user to combine highly efficient QM and MM programs, providing unified access to a wide range of computational methods. The suite also supplies additional modules with extra functionalities. These modules facilitate common tasks such as performing the setup of the models or process the data generated during the simulations. The design of QM3 has been carried out considering the least number of external dependencies (only an algebra library, already included in the distribution), which makes it extremely portable. Also, the modular structure of the suite should help to expand and develop new computational methods

Maximum Flux Transition Paths of Conformational Change

Given two metastable states A and B of a biomolecular system, the problem is to calculate the likely paths of the transition from A to B. Such a calculation is more informative and more manageable if done for a reduced set of collective variables chosen so that paths cluster in collective variable space. The computational task becomes that of computing the "center" of such a cluster. A good way to define the center employs the concept of a committor, whose value at a point in collective variable space is the probability that a trajectory at that point will reach B before A. The committor "foliates" the transition region into a set of isocommittors. The maximum flux transition path is defined as a path that crosses each isocommittor at a point which (locally) has the highest crossing rate of distinct reactive trajectories. (This path is different from that of the MaxFlux method of Huo and Straub.) It is argued that such a path is nearer to an ideal path than others that have been proposed with the possible exception of the finite-temperature string method path. To make the calculation tractable, three approximations are introduced, yielding a path that is the solution of a nonsingular two-point boundary-value problem. For such a problem, one can construct a simple and robust algorithm. One such algorithm and its performance is discussed.Comment: 7 figure

The International Surface Pressure Databank version 2

The International Surface Pressure Databank (ISPD) is the world's largest collection of global surface and sea-level pressure observations. It was developed by extracting observations from established international archives, through international cooperation with data recovery facilitated by the Atmospheric Circulation Reconstructions over the Earth (ACRE) initiative, and directly by contributing universities, organizations, and countries. The dataset period is currently 1768–2012 and consists of three data components: observations from land stations, marine observing systems, and tropical cyclone best track pressure reports. Version 2 of the ISPD (ISPDv2) was created to be observational input for the Twentieth Century Reanalysis Project (20CR) and contains the quality control and assimilation feedback metadata from the 20CR. Since then, it has been used for various general climate and weather studies, and an updated version 3 (ISPDv3) has been used in the ERA-20C reanalysis in connection with the European Reanalysis of Global Climate Observations project (ERA-CLIM). The focus of this paper is on the ISPDv2 and the inclusion of the 20CR feedback metadata. The Research Data Archive at the National Center for Atmospheric Research provides data collection and access for the ISPDv2, and will provide access to future versions

Changing a semantics: opportunism or courage?

The generalized models for higher-order logics introduced by Leon Henkin, and their multiple offspring over the years, have become a standard tool in many areas of logic. Even so, discussion has persisted about their technical status, and perhaps even their conceptual legitimacy. This paper gives a systematic view of generalized model techniques, discusses what they mean in mathematical and philosophical terms, and presents a few technical themes and results about their role in algebraic representation, calibrating provability, lowering complexity, understanding fixed-point logics, and achieving set-theoretic absoluteness. We also show how thinking about Henkin's approach to semantics of logical systems in this generality can yield new results, dispelling the impression of adhocness. This paper is dedicated to Leon Henkin, a deep logician who has changed the way we all work, while also being an always open, modest, and encouraging colleague and friend.Comment: 27 pages. To appear in: The life and work of Leon Henkin: Essays on his contributions (Studies in Universal Logic) eds: Manzano, M., Sain, I. and Alonso, E., 201