Electron capture of tetracyanoethylene oxide in the gas phase. Rearrangement of the parent radical anion to form [(NC)(3)C](-). A joint experimental and ab initio study

Capture of an electron by tetracyanoethylene oxide can initiate a number of decomposition pathways. One of these decompositions yields [(NC)3C]− as the ionic product. Ab initio calculations (at the B3LYP/6-31+G∗ level of theory) indicate that the formation of [(NC)3C]− is initiated by capture of an electron into the LUMO of tetracyanoethylene oxide to yield the anion radical [(NC)2C–O–C(CN)2]−· that undergoes internal nucleophilic substitution to form intermediate [(NC)3C–OCCN]−·. This intermediate dissociates to form [(NC)3C]− (m/z 90) as the ionic product. The radical (NC)3C· has an electron affinity of 4.0 eV (385 kJ mol−1). Ab initio calculations show that [(NC)3C]− is trigonal planar with the negative charge mainly on the nitrogens. A pictorial representation of this structure is the resonance structure formed from three degenerate contributing structures (NC)2–CCN−. The other product of the reaction is nominally (NCCO)·, but there is no definitive experimental evidence to indicate whether this radical survives intact, or decomposes to NC· and CO. The overall process [(NC)2C–O–C(CN)2]−· → [(NC)3C]− + (NCCO)· is calculated to be endothermic by 21 kJ mol−1 with an overall barrier of 268 kJ mol−1.

Resonance of graphene nanoribbons doped with nitrogen and boron : a molecular dynamics study

Based on its enticing properties, graphene has been envisioned with applications in the area of electronics, photonics, sensors, bioapplications and others. To facilitate various applications, doping has been frequently used to manipulate the properties of graphene. Despite a number of studies conducted on doped graphene regarding its electrical and chemical properties, the impact of doping on the mechanical properties of graphene has been rarely discussed. A systematic study of the vibrational properties of graphene doped with nitrogen and boron is performed by means of a molecular dynamics simulation. The influence from different density or species of dopants has been assessed. It is found that the impacts on the quality factor, Q, resulting from different densities of dopants vary greatly, while the influence on the resonance frequency is insignificant. The reduction of the resonance frequency caused by doping with boron only is larger than the reduction caused by doping with both boron and nitrogen. This study gives a fundamental understanding of the resonance of graphene with different dopants, which may benefit their application as resonators.

Intercalation of dodecylamine into kaolinite and its layering structure investigated by molecular dynamics simulation

Dodecylamine was successfully intercalated into the layer space of kaolinite by utilizing the methanol treated kaolinite–dimethyl sulfoxide (DMSO) intercalation complex as an intermediate. The basal spacing of kaolinite, measured by X-ray diffraction (XRD), increased from 0.72 nm to 4.29 nm after the intercalation of dodecylamine. Also, the significant variation observed in the Fourier Transform Infrared Spectroscopy (FTIR) spectra of kaolinite when intercalated with dodecylamine verified the feasibility of intercalation of dodecylamine into kaolinite. Isothermal-isobaric (NPT) molecular dynamics simulation with the use of Dreiding force field was performed to probe into the layering behavior and structure of nanoconfined dodecylamine in the kaolinite gallery. The concentration profiles of the nitrogen atom, methyl group and methylene group of intercalated dodecylamine molecules in the direction perpendicular to the kaolinite basal surface indicated that the alkyl chains within the interlayer space of kaolinite exhibited an obvious layering structure. However, the unified bilayer, pseudo-trilayer, or paraffin-type arrangements of alkyl chains deduced based on their chain length combined with the measured basal spacing of organoclays were not found in this study. The alkyl chains aggregated to a mixture of ordered paraffin-type-like structure and disordered gauche conformation in the middle interlayer space of kaolinite, and some alkyl chains arranged in two bilayer structures, in which one was close to the silica tetrahedron surface, and the other was close to the alumina octahedron surface with their alkyl chains parallel to the kaolinite basal surface.

Prediction of thermal expansion properties of carbon nanotubes using molecular dynamics simulations

The axial coefficients of thermal expansion (CTE) of various carbon nanotubes (CNTs), i.e., single-wall carbon nanotubes (SWCNTs), and some multi-wall carbon nanotubes (MWCNTs), were predicted using molecular dynamics (MDs) simulations. The effects of two parameters, i.e., temperature and the CNT diameter, on CTE were investigated extensively. For all SWCNTs and MWCNTs, the obtained results clearly revealed that within a wide low temperature range, their axial CTEs are negative. As the diameter of CNTs decreases, this temperature range for negative axial CTEs becomes narrow, and positive axial CTEs appear in high temperature range. It was found that the axial CTEs vary nonlinearly with the temperature, however, they decrease linearly as the CNT diameter increases. Moreover, within a wide temperature range, a set of empirical formulations was proposed for evaluating the axial CTEs of armchair and zigzag SWCNTs using the above two parameters. Finally, it was found that the absolute value of the negative axial CTE of any MWCNT is much smaller than those of its constituent SWCNTs, and the average value of the CTEs of its constituent SWCNTs. The present fundamental study is very important for understanding the thermal behaviors of CNTs in such as nanocomposite temperature sensors, or nanoelectronics devices using CNTs.

Formation of carbon nanoscrolls from graphene nanoribbons : a molecular dynamics study

Carbon nanoscrolls (CNSs) are one of the carbon-based nanomaterials similar to carbon nanotubes (CNTs) but are not widely studied in spite of their great potential applications. Their practical applications are hindered by the challenging fabrication of the CNSs. A physical approach has been proposed recently to fabricate the CNS by rolling up a monolayer graphene nanoribbon (GNR) around a CNT driven by the interaction energy between them. In this study, we perform extensive molecular dynamics (MD) simulations to investigate the various factors that impact the formation of the CNS from GNR. Our simulation results show that the formation of the CNS is sensitive to the length of the CNT and temperature. When the GNR is functionalized with hydrogen, the formation of the CNS is determined by the density and distribution of the hydrogen atoms. Graphyne, the allotrope of graphene, is inferior to graphene in the formation of the CNS due to the weaker bonds and the associated smaller atom density. The mechanism behind the rolling of GNR into CNS lies in the balance between the GNR–CNT van der Waals (vdW) interactions and the strain energy of GNR. The present work reveals new important insights and provides useful guidelines for the fabrication of the CNS.

Thermal behavior of kaolinite–urea intercalation complex and molecular dynamics simulation for urea molecule orientation

The thermal behavior of kaolinite–urea intercalation complex was investigated by thermogravimetry–differential scanning calorimetry (TG–DSC), X-ray diffraction (XRD), and fourier transform infrared spectroscopy (FTIR). In addition, the interaction mode of urea molecules intercalated into the kaolinite gallery was studied by means of molecular dynamics simulation. Three main mass losses were observed at 136 °C, in the range of 210–270 °C, and at 500 °C in the TG–DSC curves, which were, respectively, attributed to (1) melting of the surface-adsorbed urea, (2) removal of the intercalated urea, and (3) dehydroxylation of the deintercalated kaolinite. The three DSC endothermic peaks at 218, 250, and 261 °C were related to the successive removals of intercalated urea with three different distribution structures. Based on the angle between the dipole moment vector of urea and the basal surface of kaolinite, the three urea models could be described as follows: (1) Type A, the dipole moment vector is nearly parallel to the basal surface of kaolinite; (2) Type B, the dipole moment vector points to the silica tetrahedron with the angle between it and the basal surface of kaolinite ranging from 20°to 40°; and (3) Type C, the dipole moment vector is nearly perpendicular to the basal surface of kaolinite. The three distribution structures of urea molecules were validated by the results of the molecular dynamics simulation. Furthermore, the thermal behavior of the kaolinite–urea intercalation complex investigated by TG–DSC was also supported by FTIR and XRD analyses.

Fractal and complex network analyses of protein molecular dynamics

Based on protein molecular dynamics, we investigate the fractal properties of energy, pressure and volume time series using the multifractal detrended fluctuation analysis (MF-DFA) and the topological and fractal properties of their converted horizontal visibility graphs (HVGs). The energy parameters of protein dynamics we considered are bonded potential, angle potential, dihedral potential, improper potential, kinetic energy, Van der Waals potential, electrostatic potential, total energy and potential energy. The shape of the h(q)h(q) curves from MF-DFA indicates that these time series are multifractal. The numerical values of the exponent h(2)h(2) of MF-DFA show that the series of total energy and potential energy are non-stationary and anti-persistent; the other time series are stationary and persistent apart from series of pressure (with H≈0.5H≈0.5 indicating the absence of long-range correlation). The degree distributions of their converted HVGs show that these networks are exponential. The results of fractal analysis show that fractality exists in these converted HVGs. For each energy, pressure or volume parameter, it is found that the values of h(2)h(2) of MF-DFA on the time series, exponent λλ of the exponential degree distribution and fractal dimension dBdB of their converted HVGs do not change much for different proteins (indicating some universality). We also found that after taking average over all proteins, there is a linear relationship between 〈h(2)〉〈h(2)〉 (from MF-DFA on time series) and 〈dB〉〈dB〉 of the converted HVGs for different energy, pressure and volume.

Carbon nanodot decorated graphitic carbon nitride: New insights into the enhanced photocatalytic water splitting from ab initio studies

Interfacing carbon nanodots (C-dots) with graphitic carbon nitride (g-C3N4) produces a metal-free system that has recently demonstrated significant enhancement of photo-catalytic performance for water splitting into hydrogen [Science, 2015, 347, 970–974]. However, the underlying photo-catalytic mechanism is not fully established. Herein, we have carried out density functional theory (DFT) calculations to study the interactions between g-C3N4 and trigonal/hexagonal shaped C-dots. We find that hybrid C-dots/g-C3N4 can form a type-II van der Waals heterojunction, leading to significant reduction of band gap. The C-dot decorated g-C3N4 enhances the separation of photogenerated electron and hole pairs and the composite's visible light response. Interestingly, the band alignment of C-dots and g-C3N4 calculated by the hybrid functional method indicates that C-dots act as a spectral sensitizer in hybrid C-dots/g-C3N4 for water splitting. Our results offer new theoretical insights into this metal-free photocatalyst for water splitting.

Can current force fields reproduce ring puckering in 2-O-sulfo-α-l-iduronic acid? A molecular dynamics simulation study

The monosaccharide 2-O-sulfo-α-l-iduronic acid (IdoA2S) is one of the major components of glycosaminoglycans. The ability of molecular mechanics force fields to reproduce ring-puckering conformational equilibrium is important for the successful prediction of the free energies of interaction of these carbohydrates with proteins. Here we report unconstrained molecular dynamics simulations of IdoA2S monosaccharide that were carried out to investigate the ability of commonly used force fields to reproduce its ring conformational flexibility in aqueous solution. In particular, the distribution of ring conformer populations of IdoA2S was determined. The GROMOS96 force field with the SPC/E water potential can predict successfully the dominant skew-boat to chair conformational transition of the IdoA2S monosaccharide in aqueous solution. On the other hand, the GLYCAM06 force field with the TIP3P water potential sampled transitional conformations between the boat and chair forms. Simulations using the GROMOS96 force field showed no pseudorotational equilibrium fluctuations and hence no inter-conversion between the boat and twist boat ring conformers. Calculations of theoretical proton NMR coupling constants showed that the GROMOS96 force field can predict the skew-boat to chair conformational ratio in good agreement with the experiment, whereas GLYCAM06 shows worse agreement. The omega rotamer distribution about the C5–C6 bond was predicted by both force fields to have torsions around 10°, 190°, and 360°.

Molecular dynamics simulations of CXCL-8 and its interactions with a receptor peptide, heparin fragments, and sulfated linked cyclitols

CXCL-8 (Interleukin 8) is a CXC chemokine with a central role in the human immune response. We have undertaken extensive in silico analyses to elucidate the interactions of CXCL-8 with its various binding partners, which are crucial for its biological function. Sequence and structure analyses showed that residues in the thirdq β-sheet and basic residues in the heparin binding site are highly variable, while residues in the second β-sheet are highly conserved. Molecular dynamics simulations in aqueous solution of dimeric CXCL-8 have been performed with starting geometries from both X-ray and NMR structures showed shearing movements between the two antiparallel C-terminal helices. Dynamic conservation analyses of these simulations agreed with experimental data indicating that structural differences between the two structures at quaternary level arise from changes in the secondary structure of the N-terminal loop, the 310-helix, the 30s, 40s, and 50s loops and the third β-sheet, resulting in a different interhelical separation. Nevertheless, the observation of these different states indicates that CXCL-8 has the potential to undergo conformational changes, and it seems likely that this feature is relevant to the mode of binding of glycosaminoglycan (GAG) mimetics such as cyclitols. Simulations of the receptor peptide fragment−CXCL-8 complex identified several specific interactions of the receptor peptide with CXCL-8 that could be exploited in the structure-based design of competitive peptides and nonpeptidic molecules targeting CXCL-8 for combating inflammatory diseases. Simulations of the CXCL-8 dimer complexed with a 24-mer heparin fragment and of the CXCL-8−receptor peptide complex revealed that Arg60, Lys64, and Arg68 in the dimer bind to cyclitols in a horseshoe pattern, defining a region which is spatially distinct from the receptor binding site. There appears to be an optimum number of sulfates and an optimum length of alkyl spacers required for the interaction of cyclitol inhibitors with the dimeric form of CXCL-8. Calculation of the binding affinities of cyclitol inhibitors reflected satisfactorily the ranking of experimentally determined inhibitory potencies. The findings of these molecular modeling studies will help in the search for inhibitors which can modulate various CXCL-8 biological activities and serve as an excellent model system to study CXC-inhibitor interactions.

Molecular dynamics simulation of the phosphorylation-induced conformational changes of a tau peptide fragment

Aggregation of the microtubule associated protein tau (MAPT) within neurons of the brain is the leading cause of tauopathies such as Alzheimer's disease. MAPT is a phospho-protein that is selectively phosphorylated by a number of kinases in vivo to perform its biological function. However, it may become pathogenically hyperphosphorylated, causing aggregation into paired helical filaments and neurofibrillary tangles. The phosphorylation induced conformational change on a peptide of MAPT (htau225−250) was investigated by performing molecular dynamics simulations with different phosphorylation patterns of the peptide (pThr231 and/or pSer235) in different simulation conditions to determine the effect of ionic strength and phosphate charge. All phosphorylation patterns were found to disrupt a nascent terminal β-sheet pattern (226VAVVR230 and 244QTAPVP249), replacing it with a range of structures. The double pThr231/pSer235 phosphorylation pattern at experimental ionic strength resulted in the best agreement with NMR structural characterization, with the observation of a transient α-helix (239AKSRLQT245). PPII helical conformations were only found sporadically throughout the simulations. Proteins 2014; 82:1907–1923. © 2014 Wiley Periodicals, Inc.

The synthesis and application of novel profluorescent nitroxides as probes for polymer degradation

This PhD project has expanded the knowledge in the area of profluorescent nitroxides with regard to the synthesis and characterisations of novel profluorescent nitroxide probes as well as physical characterisation of the probe molecules in various polymer/physical environments. The synthesis of the first example of an azaphenalene-based fused aromatic nitroxide TMAO, [1,1,3,3-tetramethyl-2,3-dihydro-2-azaphenalen-2-yloxyl, was described. This novel nitroxide possesses some of the structural rigidity of the isoindoline class of nitroxides, as well as some properties akin to TEMPO nitroxides. Additionally, the integral aromatic ring imparts fluorescence that is switched on by radical scavenging reactions of the nitroxide, which makes it a sensitive probe for polymer degradation. In addition to the parent TMAO, 5 other azaphenalene derivatives were successfully synthesised. This new class of nitroxide was expected to have interesting redox properties when the structure was investigated by high-level ab initio molecular orbitals theory. This was expected to have implications with biological relevance as the calculated redox potentials for the azaphenalene ring class would make them potent antioxidant compounds. The redox potentials of 25 cyclic nitroxides from four different structural classes (pyrroline, piperidine, isoindoline and azaphenalene) were determined by cyclic voltammetry in acetonitrile. It was shown that potentials related to the one electron processes of the nitroxide were influenced by the type of ring system, ring substituents or groups surrounding the moiety. Favourable comparisons were found between theoretical and experimental potentials for pyrroline, piperidine and isoindoline ring classes. Substitution of these ring classes, were correctly calculated to have a small yet predictable effect on the potentials. The redox potentials of the azaphenalene ring class were underestimated by the calculations in all cases by at least a factor of two. This is believed to be due to another process influencing the redox potentials of the azaphenalene ring class which is not taken into account by the theoretical model. It was also possible to demonstrate the use of both azaphenalene and isoindoline nitroxides as additives for monitoring radical mediated damage that occurs in polypropylene as well as in more commercially relevant polyester resins. Polymer sample doped with nitroxide were exposed to both thermo-and photo-oxidative conditions with all nitroxides showing a protective effect. It was found that isoindoline nitroxides were able to indicate radical formation in polypropylene aged at elevated temperatures via fluorescence build-up. The azaphenalene nitroxide TMAO showed no such build-up of fluorescence. This was believed to be due to the more labile bond between the nitroxide and macromolecule and the protection may occur through a classical Denisov cycle, as is expected for commercially available HAS units. Finally, A new profluorescent dinitroxide, BTMIOA (9,10-bis(1,1,3,3- tetramethylisoindolin-2-yloxyl-5-yl)anthracene), was synthesised and shown to be a powerful probe for detecting changes during the initial stages of thermo-oxidative degradation of polypropylene. This probe, which contains a 9,10-diphenylanthracene core linked to two nitroxides, possesses strongly suppressed fluorescence due to quenching by the two nitroxide groups. This molecule also showed the greatest protective effect on thermo-oxidativly aged polypropylene. Most importantly, BTMIOA was found to be a valuable tool for imaging and mapping free-radical generation in polypropylene using fluorescence microscopy.

Experimental and theoretical studies of the redox potentials of cyclic nitroxides

The redox potentials of 25 cyclic nitroxides from four different structural classes (pyrrolidine, piperidine, isoindoline, and azaphenalene) were determined experimentally by cyclic voltammetry in acetonitrile, and also via high-level ab initio molecular orbital calculations. It is shown that the potentials are influenced by the type of ring system, ring substituents and/or groups surrounding the radical moiety. For the pyrrolidine, piperidine, and isoindolines there is excellent agreement (mean absolute deviation of 0.05 V) between the calculated and experimental oxidation potentials; for the azaphenalenes, however, there is an extraordinary discrepancy (mean absolute deviation of 0.60 V), implying that their one-electron oxidation might involve additional processes not considered in the theoretical calculations. This recently developed azaphenalene class of nitroxide represents a new variant of a nitroxide ring fused to an aromatic system and details of the synthesis of five derivatives involving differing aryl substitution are also presented.

Desorption electrospray ionisation mass spectrometry of stabilised polyesters reveals activation of hindered amine light stabilisers

The use of hindered amine light stabilizers (HALS) to retard thermo- and photo-degradation of polymers has become increasingly common. Proposed mechanisms of polymer stabilisation involve significant changes to the HALS chemical structure; however, reports of the characterisation of these modified chemical species are limited. To better understand the fate of HALS and determine their in situ modifications, desorption electrospray ionisation mass spectrometry (DESI-MS) was employed to characterise ten commercially available HALS present in polyester-based coil coatings. TINUVIN® 770, 292, 144, 123, 152, and NOR371; HOSTAVIN® 3052, 3055, 3050, and 3058 were separately formulated with a pigmented, thermosetting polyester resin, cured on metal at 262 C and analysed directly by DESI-MS. High-level ab initio molecular orbital theory calculations were also undertaken to aid the mechanistic interpretation of the results. For HALS containing N-substituted piperidines (i.e., N-CH3, N-C(O)CH3, and N-OR) a secondary piperidine (N-H) analogue was detected in all cases. The formation of these intermediates can be explained either through hydrogen abstraction based mechanisms or direct N-OR homolysis with the former dominant under normal service temperatures (ca. 25-80 C), and the latter potentially becoming competitive under the high temperatures associated with curing (ca. 230-260 C). © 2013 Elsevier Ltd. All rights reserved.

Gas-phase reactions of aryl radicals with 2-butyne : an experimental and theoretical investigation employing the N-methyl-pyridinium-4-yl radical cation

Aromatic radicals form in a variety of reacting gas-phase systems, where their molecular weight growth reactions with unsaturated hydrocarbons are of considerable importance. We have investigated the ion-molecule reaction of the aromatic distonic N-methyl-pyridinium-4-yl (NMP) radical cation with 2-butyne (CH3C CCH3) using ion trap mass spectrometry. Comparison is made to high-level ab initio energy surfaces for the reaction of NMP and for the neutral phenyl radical system. The NMP radical cation reacts rapidly with 2-butyne at ambient temperature, due to the apparent absence of any barrier. The activated vinyl radical adduct predominantly dissociates via loss of a H atom, with lesser amounts of CH3 loss. High-resolution Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometry allows us to identify small quantities of the collisionally deactivated reaction adduct. Statistical reaction rate theory calculations (master equation/RRKM theory) on the NMP + 2-butyne system support our experimental findings, and indicate a mechanism that predominantly involves an allylic resonance-stabilized radical formed via H atom shuttling between the aromatic ring and the C-4 side-chain, followed by cyclization and/or low-energy H atom beta-scission reactions. A similar mechanism is demonstrated for the neutral phenyl radical (Ph center dot)+2-butyne reaction, forming products that include 3-methylindene. The collisionally deactivated reaction adduct is predicted to be quenched in the form of a resonance-stabilized methylphenylallyl radical. Experiments using a 2,5-dichloro substituted methyl-pyridiniumyl radical cation revealed that in this case CH3 loss from the 2-butyne adduct is favoured over H atom loss, verifying the key role of ortho H atoms, and the shuttling mechanism, in the reactions of aromatic radicals with alkynes. As well as being useful phenyl radical analogues, pyridiniumyl radical cations may form in the ionosphere of Titan, where they could undergo rapid molecular weight growth reactions to yield polycyclic aromatic nitrogen hydrocarbons (PANHs).