Insights into the nature of the coupling interactions between uracil corrosion inhibitors and copper : a DFT and molecular dynamics study

In this study, the nature of the coupling interactions between copper and uracil as well as its several derivatives has been systematically investigated employing the atoms in molecules (AIM) theory and energy decomposition analyses. The whole interaction process has been investigated through the analyses of the radial distribution functions of the Cu⋯X (X = S and O) contact on the basis of the ab initio molecular dynamics. No direct relationship between the adsorption strengths and inhibition efficiencies of the inhibitors has been observed. Additionally, the possibility of the methyl-substituted dithiouracil species to act as copper corrosion inhibitors has been tested.

Vacancy mediated desorption of hydrogen from a sodium alanate surface : An ab initio spin-polarized study

Ab initio spin-polarized density functional theory calculations are performed to explore the effect of single Na vacancy on NaAlH4(001) surface on the initial dehydrogenation kinetics. The authors found that two Al–H bond lengths become elongated and weakened due to the presence of a Na vacancy on the NaAlH4(001) surface. Spontaneous recombination from the surface to form molecular hydrogen is observed in the spin-polarized ab initio molecular dynamics simulation. The authors’ results indicate that surface Na vacancies play a critical role in accelerating the dehydrogenation kinetics in sodium alanate. The understanding gained here will aid in the rational design and development of complex hydride materials for hydrogen storage

Study of the initial stage of solid electrolyte interphase formation upon chemical reaction of lithium metal and N-Methyl-N-Propyl-Pyrrolidinium-Bis (Fluorosulfonyl) Imide

Chemical reaction studies of N-methyl-N-propyl-pyrrolidinium-bis(fluorosulfonyl)imide-based ionic liquid with the lithium metal surface were performed using ab initio molecular dynamics (aMD) simulations and X-ray Photoelectron Spectroscopy (XPS). The molecular dynamics simulations showed rapid and spontaneous decomposition of the ionic liquid anion, with subsequent formation of long-lived species such as lithium fluoride. The simulations also revealed the cation to retain its structure by generally moving away from the lithium surface. The XPS experiments showed evidence of decomposition of the anion, consistent with the aMD simulations and also of cation decomposition and it is envisaged that this is due to the longer time scale for the XPS experiment compared to the time scale of the aMD simulation. Overall experimental results confirm the majority of species suggested by the simulation. The rapid chemical decomposition of the ionic liquid was shown to form a solid electrolyte interphase composed of the breakdown products of the ionic liquid components in the absence of an applied voltage.

Graphene-covered perovskites: An effective strategy to enhance light absorption and resist moisture degradation

The long-term stability of methylammonium lead triiodide (MAPbI3) perovskite in moist environments is a paramount challenge to realise the commercialization of perovskite solar cells. In an attempt to address this concern, we have carried out systematic first-principles studies on the MAPbI3 perovskite with a hydrophobic graphene layer interfaced as a water barrier. We find there is a charge transfer at the graphene/MAPbI3 interface and electrons can be excited from graphene into the perovskite surface, leading to well separated electron–hole pairs, i.e. reduced recombination. By studying the optical properties, we find the hybrid graphene/MAPbI3 nanocomposite displays enhanced light absorption compared with the pristine MAPbI3. Furthermore, from an ab initio molecular dynamics simulation, the graphene/MAPbI3 nanocomposite is confirmed to be able to resist the reaction with water molecules, highlighting a great advantage of this nanocomposite in promoting long-term photovoltaic performance.

Graphene-like two-dimensional ionic boron with double dirac cones at ambient condition

Recently, partially ionic boron (γ-B28) has been predicted and observed in pure boron, in bulk phase and controlled by pressure [Nature, 457 (2009) 863]. By using ab initio evolutionary structure search, we report the prediction of ionic boron at a reduced dimension and ambient pressure, namely, the two-dimensional (2D) ionic boron. This 2D boron structure consists of graphene-like plane and B2 atom pairs, with the P6/mmm space group and 6 atoms in the unit cell, and has lower energy than the previously reported α-sheet structure and its analogues. Its dynamical and thermal stability are confirmed by the phonon-spectrum and ab initio molecular dynamics simulation. In addition, this phase exhibits double Dirac cones with massless Dirac fermions due to the significant charge transfer between the graphene-like plane and B2 pair that enhances the energetic stability of the P6/mmm boron. A Fermi velocity (vf) as high as 2.3 x 106 m/s, which is even higher than that of graphene (0.82 x 106 m/s), is predicted for the P6/mmm boron. The present work is the first report of the 2D ionic boron at atmospheric pressure. The unique electronic structure renders the 2D ionic boron a promising 2D material for applications in nanoelectronics.

C-BN Single-Walled Nanotubes from Hybrid Connection of BN/C Nanoribbons: Prediction by ab initio Density Functional Calculations

We demonstrated for the first time by ab initio density functional calculation and molecular dynamics simulation that C0.5(BN)0.5 armchair single-walled nanotubes (NT) are gapless semiconductors and can be spontaneously formed via the hybrid connection of graphene/BN Nanoribbons (GNR/BNNR) at room temperature. The direct synthesis of armchair C0.5(BN)0.5 via the hybrid connection of GNR/BNNR is predicted to be both thermodynamically and dynamically stable. Such novel armchair C0.5(BN)0.5 NTs possess enhanced conductance as that observed in GNRs. Additionally, the zigzag C0.5(BN)0.5 SWNTs are narrow band gap semiconductors, which may have potential application for light emission. In light of recent experimental progress and the enhanced degree of control in the synthesis of GNRs and BNNR, our results highlight an interesting avenue for synthesizing a novel specific type of C0.5(BN)0.5 nanotube (gapless or narrow direct gap semiconductor), with potentially important applications in BNC-based nanodevices.

Hydrogen spillover mechanism on a pd-doped mg surface as revealed by ab initio density functional calculation

The hydrogenation kinetics of Mg is slow, impeding its application for mobile hydrogen storage. We demonstrate by ab initio density functional theory (DFT) calculations that the reaction path can be greatly modified by adding transition metal catalysts. Contrasting with Ti doping, a Pd dopant will result in a very small activation barrier for both dissociation of molecular hydrogen and diffusion of atomic H on the Mg surface. This new computational finding supports for the first time by ab initio simulationthe proposed hydrogen spillover mechanism for rationalizing experimentally observed fast hydrogenation kinetics for Pd-capped Mg materials.

Structural and electronic properties of diazonium functionalized (4, 4) single walled carbon nanotube : an ab initio study

In this work, ab initio density functional theory (DFT) calculations are performed to study the structural and electronic properties of diazonium reagent functionalized (4, 4) single-walled carbon nanotube (SWCNT). We find the aryl group covalently bonds with SWCNT and prefers to be perpendicular to the side wall of nanotube. It has a rotational barrier of 0.35 eV around the formed aryl-tube bond axis and should be thermodynamically stable at room temperature. Additionally, new peaks appeared around the Fermi energy in the density of state (DOS) due to the weak band dispersion. Increasing of the coverage of the functional group will result in significant upshift of the Fermi level.

The role of V2O5 on the dehydrogenation and hydrogenation in magnesium hydride : An ab initio study

Ab initio density functional theory calculations are performed to study the experimentally observed catalytic role of V2O5 in the recycling of hydrogen in magnesium hydride. We find that the Mg–H bond length becomes elongated when MgH2 clusters are positioned on single, two, and three coordinated oxygen sites (O1, O2, and O3) on the V2O5(001) surface. Molecular hydrogen is predicted to spontaneously form at the hole site on the V2O5(001) surface. Additionally, the activation barrier for the dissociation of hydrogen on V-doped Mg(0001) surface is 0.20 eV, which is only 1/5 of that on pure Mg(0001) surface. Our results indicate that oxygen sites on the V2O5(001)surface and the V dopant in Mg may be important facilitators for dehydrogenation and rehydrogenation, respectively. The understanding gained here will aid in the rational design and development of Mg-based hydrogen storage materials.

A coupled meshless technique/molecular dynamics approach for deformation characterization of monocrystalline metal

This paper presents a multiscale study using the coupled Meshless technique/Molecular Dynamics (M2) for exploring the deformation mechanism of mono-crystalline metal (focus on copper) under uniaxial tension. In M2, an advanced transition algorithm using transition particles is employed to ensure the compatibility of both displacements and their gradients, and an effective local quasi-continuum approach is also applied to obtain the equivalent continuum strain energy density based on the atomistic poentials and Cauchy-Born rule. The key parameters used in M2 are firstly investigated using a benchmark problem. Then M2 is applied to the multiscale simulation for a mono-crystalline copper bar. It has found that the mono-crystalline copper has very good elongation property, and the ultimate strength and Young's modulus are much higher than those obtained in macro-scale.

Molecular dynamics study of dynamic buckling properties of nanowires with defects

Based on the embedded atom method (EAM) and molecular dynamics (MD) method, the deformation properties of Cu nanowires with different single defects under dynamic compression have been studied. The mechanical behaviours of the perfect nanowire are first studied, and the critical stress decreases with the increase of the nanowire’s length, which is well agreed with the modified Euler theory. We then consider the effects to the buckling phenomenon resulted from different defects. It is found that obvious decrease of the critical stress is resulted from different defects, and the largest decrease is found in nanowire with the surface vertical defect. Surface defects are found exerting larger influence than internal defects. The buckling duration is found shortened due to different defects except the nanowire with surface horizon defect, which is also found possessing the largest deflection. Different deflections are also observed for different defected nanowires. It is find that due to surface defects, only deflection in one direction is happened, but for internal defects, more complex deflection circumstances are observed.

Molecular dynamics studies of the effects of branching characteristics on the crystalline structure of polyethylene

Molecular dynamics simulations were carried out on single chain models of linear low-density polyethylene in vacuum to study the effects of branch length, branch content, and branch distribution on the polymer’s crystalline structure at 300 K. The trans/gauche (t/g) ratios of the backbones of the modeled molecules were calculated and utilized to characterize their degree of crystallinity. The results show that the t/g ratio decreases with increasing branch content regardless of branch length and branch distribution, indicating that branch content is the key molecular parameter that controls the degree of crystallinity. Although t/g ratios of the models with the same branch content vary, they are of secondary importance. However, our data suggests that branch distribution (regular or random) has a significant effect on the degree of crystallinity for models containing 10 hexyl branches/1,000 backbone carbons. The fractions of branches that resided in the equilibrium crystalline structures of the models were also calculated. On average, 9.8% and 2.5% of the branches were found in the crystallites of the molecules with ethyl and hexyl branches while C13 NMR experiments showed that the respective probabilities of branch inclusion for ethyl and hexyl branches are 10% and 6% [Hosoda et al., Polymer 1990, 31, 1999–2005]. However, the degree of branch inclusion seems to be insensitive to the branch content and branch distribution.

Molecular dynamics investigation on edge stress and shape transition in graphene nanoribbons

Graphene nanoribbon (GNR) with free edges demonstrates unique pre-existing edge energy and edge stress, leading to non-flat morphologies. Using molecular dynamics (MD) methods, we evaluated edge energies as well as edge stresses for four different edge types, including regular edges (armchair and zigzag), armchair edge terminated with hydrogen and reconstructed armchair. The results showed that compressive stress exists in the regular and hydrogen-terminated edges along the edge direction. In contrast, the reconstructed armchair edge is generally subject to tension. Furthermore, we also investigated shape transition between flat and rippled configurations of GNRs with different free edges. It was found that the pre-existing stress at free edges can greatly influence the initial energy state and the shape transition.

Formation of single-walled carbon nanotube via the interaction of graphene nanoribbons : ab initio density functional calculations

The interaction of bare graphene nanoribbons (GNRs) was investigated by ab initio density functional theory calculations with both the local density approximation (LDA) and the generalized gradient approximation (GGA). Remarkably, two bare 8-GNRs with zigzag-shaped edges are predicted to form an (8, 8) armchair single-wall carbon nanotube (SWCNT) without any obvious activation barrier. The formation of a (10, 0) zigzag SWCNT from two bare 10-GNRs with armchair-shaped edges has activation barriers of 0.23 and 0.61 eV for using the LDA and the revised PBE exchange correlation functional, respectively, Our results suggest a possible route to control the growth of specific types SWCNT via the interaction of GNRs.