Ab-initio calculations for a realistic sensor: A study of CO sensors based on nitrogen-rich carbon nanotubes

The use of nanoscale low-dimensional systems could boost the sensitivity of gas sensors. In this work we simulate a nanoscopic sensor based on carbon nanotubes with a large number of binding sites using ab initio density functional electronic structure calculations coupled to the Non-Equilibrium Green's Function formalism. We present a recipe where the adsorption process is studied followed by conductance calculations of a single defect system and of more realistic disordered system considering different coverages of molecules as one would expect experimentally. We found that the sensitivity of the disordered system is enhanced by a factor of 5 when compared to the single defect one. Finally, our results from the atomistic electronic transport are used as input to a simple model that connects them to experimental parameters such as temperature and partial gas pressure, providing a procedure for simulating a realistic nanoscopic gas sensor. Using this methodology we show that nitrogen-rich carbon nanotubes could work at room temperature with extremely high sensitivity. Copyright 2012 Author(s). This article is distributed under a Creative Commons Attribution 3.0 Unported License. [http://dx.doi.org/10.1063/1.4739280]

Photodynamic Efficiency of Cationic meso-Porphyrins at Lipid Bilayers: Insights from Molecular Dynamics Simulations

Porphyrin derivatives have applications as photoactive drugs in photodynamic therapy. However, little is known about their interactions with phospholipid membranes at the molecular level. We employed molecular dynamics simulations to model the binding between a series of cationic meso-(N-methyl-4-pyridinium)phenylporphyrins and anionic phosphatidylglycerol lipid bilayers. This was done in the presence of molecular oxygen within the membrane. The ability of various porphyrins to cause photodamage was quantified in terms of their immersion depth and degree of exposition to a higher oxygen concentration inside the membrane. Simulations showed that the photodynamic efficiency could be improved as the number of hydrophobic phenyl substituents attached to the porphyrinic ring increased. In the specific case of porphyrins containing two hydrophobic and two charged substituents, the cis isomer was significantly more efficient than the trans. These results correlate well with previous experimental observations. They highlight the importance of both the total charge and amphiphilicity of the photosensitizer for its performance in photodynamic therapy.

Influence of Molecular Dynamics on the Dielectric Properties of Poly(9,9-di-n-octylfluorene-alt-benzothiadiazole) -Based Devices

This paper uses Nuclear Magnetic Resonance (NMR) and Differential Scanning Calorimetry (DSC) techniques to study the molecular relaxations and phase transitions in poly(9,9-di-n-octylfluorene-alt-benzothiadiazole) (F8BT), which has been extensively studied as the active thin film in organic devices. Besides the identification of the glass transition, beta relaxation and crystal-to-crystal phase transition, we correlate such phenomena with dielectric and transport mechanisms in diodes with F8BT as the active layer. The beta relaxation has been assigned to a transition at about 210 K measured by H-1 and C-13 solid state NMR, and can be attributed to local motions in the side chains. The glass transition has been detected by DSC and H-1 NMR. Dielectric spectroscopy (DS) carried out at low frequencies on diodes made from F8BT show two peaks which are coincident with the above transitions. This allowed us to correlate the electrical changes in the film with the onset of specific molecular motions. In addition, DS indicates a third peak related with a crystal-to-crystal phase transition. Finally, these transitions were correlated with changes in the carrier mobility recorded in thin films and published recently.

Effect of the thermostat in the molecular dynamics simulation on the folding of the model protein chignolin

Molecular dynamics simulations of the model protein chignolin with explicit solvent were carried out, in order to analyze the influence of the Berendsen thermostat on the evolution and folding of the peptide. The dependence of the peptide behavior on temperature was tested with the commonly employed thermostat scheme consisting of one thermostat for the protein and another for the solvent. The thermostat coupling time of the protein was increased to infinity, when the protein is not in direct contact with the thermal bath, a situation known as minimally invasive thermostat. In agreement with other works, it was observed that only in the last situation the instantaneous temperature of the model protein obeys a canonical distribution. As for the folding studies, it was shown that, in the applications of the commonly utilized thermostat schemes, the systems are trapped in local minima regions from which it has difficulty escaping. With the minimally invasive thermostat the time that the protein needs to fold was reduced by two to three times. These results show that the obstacles to the evolution of the extended peptide to the folded structure can be overcome when the temperature of the peptide is not directly controlled.

Molecular Dynamics Investigations of PRODAN in a DLPC Bilayer

Molecular dynamics computer simulations have been performed to identify preferred positions of the fluorescent probe PRODAN in a fully hydrated DLPC bilayer in the fluid phase. In addition to the intramolecular charge-transfer first vertical excited state, we considered different charge distributions for the electronic ground state of the PRODAN molecule by distinct atomic charge models corresponding to the probe molecule in vacuum as well as polarized in a weak and a strong dielectric solvent (cyclohexane and water). Independent on the charge distribution model of PRODAN, we observed a preferential orientation of this molecule in the bilayer with the dimethylamino group pointing toward the membrane's center and the carbonyl oxygen toward the membrane's interface. However, changing the charge distribution model of PRODAN, independent of its initial position in the equilibrated DLPC membrane, we observed different preferential positions. For the ground state representation without polarization and the in-cyclohexane polarization, the probe maintains its position close to the membrane's center. Considering the in-water polarization model, the probe approaches more of the polar headgroup region of the bilayer, with a strong structural correlation with the choline group, exposing its oxygen atom to water molecules. PRODAN's representation of the first vertical excited state with the in-water polarization also approaches the polar region of the membrane with the oxygen atom exposed to the bilayer's hydration shell. However, this model presents a stronger structural correlation with the phosphate groups than the ground state. Therefore, we conclude that the orientation of the PRODAN molecule inside the DLPC membrane is well-defined, but its position is very sensitive to the effect of the medium polarization included here by different models for the atomic charge distribution of the probe.

The T-2 phase in the Nb-Si-B system studied by ab initio calculations and synchrotron X-ray diffraction

The solid solution based on Nb5Si3 (Cr5B3 structure type, D8(l), tl32, 14/mcm, No140, a=6.5767 angstrom, c=11.8967 angstrom) in the Nb-Si-B system was studied from the structural and thermodynamic point of view both experimentally and by ab initio calculations. Rietveld refinement of powder X-ray synchrotron data allowed to determine the boron to silicon substitution mechanism and the structural parameters. Ab initio calculations of different ordered compounds and selected disordered alloys allowed to obtain in addition to the enthalpy of formation of the solution, substitution mechanism and structural parameters which are in excellent agreement with the experimental data. The stability of the phase is discussed. (C) 2012 Elsevier Inc. All rights reserved.

Selectivity and Mechanisms Driven by Reaction Dynamics: The Case of the Gas-Phase OH-+CH3ONO2 Reaction

Well-established statistical approaches such as transition-state theory based on high-level calculated potential energy profiles are unable to account for the selectivity observed in the gas-phase OH- + CH3ONO2 reaction. This reaction can undergo bimolecular nucleophilic displacement at either the carbon center (S(N)2@C) or the nitrogen center (S(N)2@N) as well as a proton abstraction followed by dissociation (E(CO)2) pathway. Direct dynamics simulations yield an S(N)2:E(CO)2 product ratio in close agreement with experiment and show that the lack of reactivity at the nitrogen atom is due to the highly negative electrostatic potential generated by the oxygen atoms in the ONO2 group that scatters the incoming OH-. In addition to these dynamical effects, the nonstatistical behavior of these reactions is attributed to the absence of equilibrated reactant complexes and to the large number of recrossings, which might be present in several ion-molecule gas-phase reactions.

Characterization of suramin binding sites on the human group IIA secreted phospholipase A(2) by site-directed mutagenesis and molecular dynamics simulation

Suramin is a polysulphonated naphthylurea with inhibitory activity against the human secreted group IIA phospholipase A(2) (hsPLA2GIIA), and we have investigated suramin binding to recombinant hsPLA2GIIA using site-directed mutagenesis and molecular dynamics (MD) simulations. The changes in suramin binding affinity of 13 cationic residue mutants of the hsPLA2GIIA was strongly correlated with alterations in the inhibition of membrane damaging activity of the protein. Suramin binding to hsPLA2GIIA was also studied by MD simulations, which demonstrated that altered intermolecular potential energy of the suramin/mutant complexes was a reliable indicator of affinity change. Although residues in the C-terminal region play a major role in the stabilization of the hsPLA2GIIA/suramin complex, attractive and repulsive hydrophobic and electrostatic interactions with residues throughout the protein together with the adoption of a bent suramin conformation, all contribute to the stability of the complex. Analysis of the h5PLA2GIIA/suramin interactions allows the prediction of the properties of suramin analogues with improved binding and higher affinities which may be candidates for novel phospholipase A(2) inhibitors. (C) 2012 Elsevier Inc. All rights reserved.

Molecular dynamics of poly(Ethylene Glycol) intercalated in clay, studied using 13C solid-state NMR

In this study, Cross-Polarization Magic-angle Spinning CP/MAS, 2D Exchange, Centerband-Only Detection of Exchange (CODEX), and Separated-Local-Field (SLF) NMR experiments were used to study the molecular dynamics of poly(ethylene glycol) (PEG) inside Hectorite/PEG intercalation compounds in both single- and double-layer configurations. The results revealed that the overall amplitude of the motions of the PEG chain in the single-layer configuration is considerably smaller than that observed for the double-layer intercalation compound. This result indicates that the effect of having the polymer chain interacting with both clay platelets is to produce a substantial decrease in the motional amplitudes of those chains. The presence of these dynamically restricted segments might be explained by the presence of anchoring points between the clay platelets and the PEG oxygen atoms, which was induced by the Na+ cations. By comparing the PEG motional amplitudes of the double-layered nanocomposites composed of polymers with different molecular weights, a decrease in the motional amplitude for the smaller PEG chain was observed, which might also be understood using the presence of anchoring points.

A molecular dynamics study of the influence of chain branching on the properties of polymer systems

Die vorliegende Arbeit beschäftigt sich mit dem Einfluß von Kettenverzweigungen unterschiedlicher Topologien auf die statischen Eigenschaften von Polymeren. Diese Untersuchungen werden mit Hilfe von Monte-Carlo- und Molekular-Dynamik-Simulationen durchgeführt.Zunächst werden einige theoretische Konzepte und Modelle eingeführt, welche die Beschreibung von Polymerketten auf mesoskopischen Längenskalen gestatten. Es werden wichtige Bestimmungsgrößen eingeführt und erläutert, welche zur quantitativen Charakterisierung von Verzweigungsstrukturen bei Polymeren geeignet sind. Es wird ebenso auf die verwendeten Optimierungstechniken eingegangen, die bei der Implementierung des Computerprogrammes Verwendung fanden. Untersucht werden neben linearen Polymerketten unterschiedliche Topolgien -Sternpolymere mit variabler Armzahl, Übergang von Sternpolymeren zu linearen Polymeren, Ketten mit variabler Zahl von Seitenketten, reguläre Dendrimere und hyperverzweigte Strukturen - in Abhängigkeit von der Lösungsmittelqualität. Es wird zunächst eine gründliche Analyse des verwendeten Simulationsmodells an sehr langen linearen Einzelketten vorgenommen. Die Skalierungseigenschaften der linearen Ketten werden untersucht in dem gesamten Lösungsmittelbereich vom guten Lösungsmittel bis hin zu weitgehend kollabierten Ketten im schlechten Lösungsmittel. Ein wichtiges Ergebnis dieser Arbeit ist die Bestätigung der Korrekturen zum Skalenverhalten des hydrodynamischen Radius Rh. Dieses Ergebnis war möglich aufgrund der großen gewählten Kettenlängen und der hohen Qualität der erhaltenen Daten in dieser Arbeit, insbesondere bei den linearen ketten, und es steht im Widerspruch zu vielen bisherigen Simulations-Studien und experimentellen Arbeiten. Diese Korrekturen zum Skalenverhalten wurden nicht nur für die linearen Ketten, sondern auch für Sternpolymere mit unterchiedlicher Armzahl gezeigt. Für lineare Ketten wird der Einfluß von Polydispersität untersucht.Es wird gezeigt, daß eine eindeutige Abbildung von Längenskalen zwischen Simulationsmodell und Experiment nicht möglich ist, da die zu diesem Zweck verwendete dimensionslose Größe eine zu schwache Abhängigkeit von der Polymerisation der Ketten besitzt. Ein Vergleich von Simulationsdaten mit industriellem Low-Density-Polyäthylen(LDPE) zeigt, daß LDPE in Form von stark verzweigten Ketten vorliegt.Für reguläre Dendrimere konnte ein hochgradiges Zurückfalten der Arme in die innere Kernregion nachgewiesen werden.

Comparison of model potentials for molecular dynamics simulation of crystalline silica

Simulationen von SiO2 mit dem von van Beest, Kramer und vanSanten (BKS) entwickelten Paarpotenzial erzeugen vielezufriedenstellende Ergebnisse, aber auch charakteristischeSchwachstellen. In dieser Arbeit wird das BKS-Potenzial mitzwei kürzlich vorgeschlagenen Potenzialen verglichen, dieeffektiv Mehrteilchen-Wechselwirkungen beinhalten. Der ersteAnsatz erlaubt dazu fluktuierende Ladungen, der zweiteinduzierbare Polarisierungen auf den Sauerstoffatomen. Die untersuchten Schwachstellen des BKS Potenzialsbeinhalten das Verhältnis der zwei Gitterkonstanten a und cim Quarzübergang, das von BKS falsch beschrieben wird.Cristobalit und Tridymit erscheinen instabil mit BKS.Weiterhin zeigt die BKS-Zustandsdichte charakteristischeAbweichungen von der wahren Zustandsdichte. DerÜbergangsdruck für den Stishovit I-II Übergang wird deutlichüberschätzt. Das Fluktuierende-Ladungs-Modell verbesserteinige der genannten Punkte, reproduziert aber viele andereEigenschaften schlechter als BKS. DasFluktierende-Dipol-Modell dagegen behebt alle genanntenArtefakte. Zusätzlich wird der druckinduzierte Phasenübergang imalpha-Quarz untersucht. Alle Potentiale finden die selbeStruktur für Quarz II. Bei anschliessender Dekompressionerzeugt BKS eine weitere Phase, während die beiden anderenPotentiale wieder zum alpha-Quarz zurückkehren. Weiterhinwerden zwei Methoden entwickelt, um die piezoelektrischenKonstanten bei konstantem Druck zu bestimmen. Die Ergebnissegeben Hinweise auf eine möglicherweisenicht-elektrostatische Natur der Polarisierungen imFluktuierende-Dipole-Modell. Mit dieser Interpretation scheint das Fluktuierende-DipolPotential alle verfügbaren experimentellen Daten am bestenvon allen drei untersuchten Ansätzen zu reproduzieren.

Polymer-liquid crystal interface: a molecular dynamics study

Liquid crystals (LCs) are an interesting class of soft condensed matter systems characterized by an unusual combination of fluidity and long-range order, mainly known for their applications in displays (LCDs). However, the interest in LC continues to grow pushed by their application in new technologies in medicine, optical imaging, micro and nano technologies etc. In LCDs uniaxial alignment of LCs is mainly achieved by a rubbing process. During this treatment, the surfaces of polymer coated display substrates are rubbed in one direction by a rotating cylinder covered with a rubbing cloth. Basically, LC alignment involves two possible aligning directions: uniaxial planar (homogeneous) and vertical (homeotropic) to the display substrate. An interesting unresolved question concerning LCs regards the origin of their alignment on rubbed surfaces, and in particular on the polymeric ones used in the display industry. Most studies have shown that LCs on the surface of the rubbed polymer film layer are lying parallel to the rubbing direction. In these systems, micrometric grooves are generated on the film surface along the rubbing direction and also the polymer chains are stretched in this direction. Both the parallel aligned microgrooves and the polymer chains at the film surface may play a role in the LC alignment and it is not easy to quantify the effect of each contribution. The work described in this thesis is an attempt to find new microscopic evidences on the origin of LC alignment on polymeric surfaces through molecular dynamics (MD) simulations, which allow the investigation of the phenomenon with atomic detail. The importance of the arrangement of the polymeric chains in LCs alignment was studied by performing MD simulations of a thin film of a typical nematic LC, 4-cyano-4’-pentylbiphenyl (5CB), in contact with two different polymers: poly(methyl methacrylate)(PMMA) and polystyrene (PS). At least four factors are believed to influence the LC alignment: 1. the interactions of LCs with the backbone vinyl chains; 2. the interactions of LCs with the oriented side groups; 3. the anisotropic interactions of LCs with nanometric grooves; 4. the presence of static surface charges. Here we exclude the effect of microgrooves and of static surface charges from our virtual experiment, by using flat and neutral polymer surfaces, with the aim of isolating the chemical driving factors influencing the alignment of LC phases on polymeric surfaces.