Simulating startup shear of entangled polymer melts

Start-up shear rheology is a standard experiment used for characterizing polymer flow, and to test various models of polymer dynamics. A rich phenomenology is developed for behavior of entangled monodisperse linear polymers in such tests, documenting shear stress overshoots as a function of shear rates and molecular weights. A tube theory does a reasonable qualitative job at describing these phenomena, although it involves several drastic approximations and the agreement can be fortuitous. Recently, Lu and coworkers published several papers [e.g. Lu {\it et al.} {\it ACS Macro Lett}. 2014, 3, 569-573] reporting results from molecular dynamics simulations of linear entangled polymers, which contradict both theory and experiment. Based on these observations, they made very serious conclusions about the tube theory, which seem to be premature. In this letter, we repeat simulations of Lu {\it et al.} and systematically show that neither their simulation results, nor their comparison with theory are confirmed.

Selectivity in the mechanism of action of antimicrobial mastoparan peptide Polybia-MP1

Many potent antimicrobial peptides also present hemolytic activity, an undesired collateral effect for the therapeutic application. Unlike other mastoparan peptides, Polybia-MP1 (IDWKKLLDAAKQIL), obtained from the venom of the social wasp Polybia paulista, is highly selective of bacterial cells. The study of its mechanism of action demonstrated that it permeates vesicles at a greater rate of leakage on the anionic over the zwitterionic, impaired by the presence of cholesterol or cardiolipin; its lytic activity is characterized by a threshold peptide to lipid molar ratio that depends on the phospholipid composition of the vesicles. At these particular threshold concentrations, the apparent average pore number is distinctive between anionic and zwitterionic vesicles, suggesting that pores are similarly formed depending on the ionic character of the bilayer. To prospect the molecular reasons for the strengthened selectivity in Polybia-MP1 and its absence in Mastoparan-X, MD simulations were carried out. Both peptides presented amphipathic alpha-helical structures, as previously observed in Circular Dichroism spectra, with important differences in the extension and stability of the helix; their backbone solvation analysis also indicate a different profile, suggesting that the selectivity of Polybia-MP1 is a consequence of the distribution of the charged and polar residues along the peptide helix, and on how the solvent molecules orient themselves according to these electrostatic interactions. We suggest that the lack of hemolytic activity of Polybia-MP1 is due to the presence and position of Asp residues that enable the equilibrium of electrostatic interactions and favor the preference for the more hydrophilic environment.

Study of the mechanism of action of anoplin, a helical antimicrobial decapeptide with ion channel-like activity, and the role of the amidated C-terminus

Anoplin, an antimicrobial, helical decapeptide from wasp venom, looses its biological activities by mere deamidation of its C-terminus. Secondary structure determination, by circular dichroism spectroscopy in amphipathic environments, and lytic activity in zwitterionic and anionic vesicles showed quite similar results for the amidated and the carboxylated forms of the peptide. The deamidation of the C-terminus introduced a negative charge at an all-positive charged peptide, causing a loss of amphipathicity, as indicated by molecular dynamics simulations in TFE/water mixtures and this subtle modification in a peptide`s primary structure disturbed the interaction with bilayers and biological membranes. Although being poorly lytic, the amidated form, but not the carboxylated, presented ion channel-like activity on anionic bilayers with a well-defined conductance step; at approximately the same concentration it showed antimicrobial activity. The pores remain open at trans-negative potentials, preferentially conducting cations, and this situation is equivalent to the interaction of the peptide with bacterial membranes that also maintain a high negative potential inside. Copyright (C) 2007 European Peptide Society and John Wiley & Sons, Ltd.

Micelles forming biaxial lyotropic nematic phases

The paper by Yu and Saupe on the first biaxial nematic phase created excitement for a number of reasons. Some theories of biaxial phases already existed, but experimental observation was still lacking. The phase was discovered in a lyotropic system with three components, which in theory is difficult. Lyotropic liquid crystals are composed of supramolecular assemblies of amphiphilic molecules, which may change shape and size as a function of concentration and temperature. The experimental phase diagram of the lyotropic biaxial phase was rather complex, with the biaxial region inserted between nematic cylindrical and nematic discotic phases via second-order transitions. In addition, re-entrant behaviour was evident. Saupe investigated further systems experimentally, observing that the biaxial phase might be absent in cases where a direct transition between the cylindrical and discotic phases occurred. He provided a range of theoretical and experimental contributions on the properties of these lyotropics, but was very cautious regarding the detailed amphiphilic assemblies involved. The present paper reviews this area, focusing on proposals for the structure of the micellar assemblies. Emphasis is placed on recent papers which indicate a transformation of the two uniaxial shapes, in mixing conditions, both from the theoretical and the experimental point of view, and to questions still requiring further study.

Mixture of changing uniaxial micellar forms in lyotropic biaxial nematics

Lyotropic nematics consisting of surfactant-cosurfactant water solutions may present a biaxial phase or direct U(+) <-> U(-) transitions, in different regions of the temperature-relative concentration phase diagram, for different systems and compositions. We propose that these may be related to changes of uniaxial micellar form, which may occur either smoothly or abruptly. Smooth change of cylinder-like into disc-like shapes requires a distribution of Maier-Saupe interaction constants and we consider two limiting cases for the distribution of forms: a single Gaussian and a double Gaussian. Alternatively, an abrupt change of form is described by a discontinuous distribution of interaction constants. Our results show that the dispersive distributions yield a biaxial phase, while an abrupt change of shape leads to coexistence of uniaxial phases. Fitting the theory to the experiment for the ternary system KL/decanol/D2O leads to transition lines in very good agreement with experimental results. In order to rationalise the results of the comparison, we analyse temperature and concentration form dependence, which connects micellar and experimental macroscopic parameters. Physically consistent variations of micellar asymmetry, amphiphile partitioning and volume are obtained. To the best of the authors` knowledge, this is the first truly statistical microscopic approach that is able to model experimentally observed lyotropic biaxial nematic phases.

Thermodynamic Stabillity of Hydrogen-Bonded Systems in Polar and Nonpolar Environments

The thermodynamic properties of a selected set of benchmark hydrogen-bonded systems (acetic acid dimer and the complexes of acetic acid with acetamide and methanol) was studied with the goal of obtaining detailed information on solvent effects on the hydrogen-bonded interactions using water, chloroform, and n-heptane as representatives for a wide range in the dielectric constant. Solvent effects were investigated using both explicit and implicit solvation models. For the explicit description of the solvent, molecular dynamics and Monte Carlo simulations in the isothermal isobaric (NpT) ensemble combined with the free energy perturbation technique were performed to determine solvation free energies. Within the implicit solvation approach, the polarizable continuum model and the conductor-like screening model were applied. Combination of gas phase results with the results obtained from the different solvation models through an appropriate thermodynamic cycle allows estimation of complexation free energies, enthalpies, and the respective entropic contributions in solution. Owing to the strong solvation effects of water the cyclic acetic acid dimer is not stable in aqueous solution. In less polar solvents the double hydrogen bond structure of the acetic acid dimer remains stable. This finding is in agreement with previous theoretical and experimental results. A similar trend as for the acetic acid dimer is also observed for the acetamide complex. The methanol complex was found to be thermodynamically unstable in gas phase as well as in any of the three solvents. (C) 2010 Wiley Periodicals, Inc. J Comput Chem 31: 2046-2055, 2010

Surfactant-nanotube interactions in water and nanotube separation by diameter: atomistic simulations

A non-destructive sorting method to separate single-walled carbon nanotubes (SWNTs) by diameter was recently proposed. By this method, SWNTs are suspended in water by surfactant encapsulation and the separation is carried out by ultracentrifugation in a density gradient. SWNTs of different diameters are distributed according to their densities along the centrifuge tube. A mixture of two anionic surfactants, namely sodium dodecylsulfate (SDS) and sodium cholate (SC), presented the best performance in discriminating nanotubes by diameter. Unexpectedly, small diameter nanotubes are found at the low density part of the centrifuge tube. We present molecular dynamics studies of the water-surfactant-SWNT system to investigate the role of surfactants in the sorting process. We found that surfactants can actually be attracted towards the interior of the nanotube cage, depending on the relationship between the surfactant radius of gyration and the nanotube diameter. The dynamics at room temperature showed that, as the amphiphile moves to the hollow cage, water molecules are dragged together, thereby promoting the nanotube filling. The resulting densities of filled SWNT are in agreement with measured densities.

Role of Surfactants in Carbon Nanotubes Density Gradient Separation

Several strategies aimed at sorting single-walled carbon nanotubes (SWNT) by diameter and/or electronic structure have been developed in recent years. A nondestructive sorting method was recently proposed in which nanotube bundles are dispersed in water-surfactant solutions and submitted to ultracentrifugation in a density gradient. By this method, SWNTs of different diameters are distributed according to their densities along the centrifuge tube. A mixture of two anionic amphiphiles, namely sodium dodecylsulfate (SIDS) and sodium cholate (SC), presented the best performance in discriminating nanotubes by diameter. We present molecular dynamics studies of the water-surfactant-SWNT system. The simulations revealed one aspect of the discriminating power of surfactants: they can actually be attracted toward the interior of the nanotube cage. The binding energies of SDS and SC on the outer nanotube surface are very similar and depend weakly on diameter. The binding inside the tubes, on the contrary, is strongly diameter dependent: SDS fits best inside tubes with diameters ranging from 8 to 9 angstrom, while SC is best accommodated in larger tubes, with diameters in the range 10.5-12 angstrom. The dynamics at room temperature showed that, as the amphiphile moves to the hollow cage, water molecules are dragged together, thereby promoting the nanotube filling. The resulting densities of filled SWNT are in agreement with measured densities.

Classical simulation of deposition of thiophene oligomers on TiO(2)-anatase: Relevance of long-range electrostatic interactions

We performed classical molecular dynamics simulations of the vapor-deposition of alpha-T4 oligomers on the TiO(2)-anatase (101) surface, comparing different sets of charges associated with the atoms of the model. The potential energy surfaces for alpha-T4 and TiO(2) were described by re-parametrizations of the Universal force field with charges given by the charge equilibration (QEq) scheme, or with fixed charges obtained by an ab initio method using the Hirshfeld partition. The two sets of charges lead to completely different results for the interface formation, and for the characteristics of the organic film, with a clearly defined alpha-T4 contact layer in the QEq case, and a more homogeneous molecular distribution when using Hirshfeld charges. The main reason for the discrepancy was found to be the incorrect charge assignment given by QEq to the sulfur and alpha-carbon atoms in thiophenes, and highlight the relevance of long-range interactions in the organization of molecular films. (C) 2009 Elsevier B.V. All rights reserved.

Spin current in the Mobius cyclacene belts

We report on density functional theory studies of the electronic structure and magnetic properties of Mobius-[n]cyclacenes. The geometry of Mobius bands presents a modulation of bond lengths that is needed to accommodate the twist. This modulation takes the form of bond alternation defects analogous to those of solitons in polyacetylene. The ground state of all Mobius bands is a triplet, with a spin density distribution that follows the bond length modulation. A molecular dynamics simulation of the Mobius cyclacene at 300 K shows that the twist travels around the belt inducing a magnetic current. (C) 2009 Elsevier B. V. All rights reserved.

Theoretical investigation of Hf and Zr defects in c-Ge

The introduction of high-permittivity gate dielectric materials into complementary metal oxide semiconductor technology has reopened the interest in Ge as a channel material mainly due to its high hole mobility. Since HfO(2) and ZrO(2) are two of the most promising dielectric candidates, it is important to investigate if Hf and Zr may diffuse into the Ge channel. Therefore, using ab initio density functional theory calculations, we have studied substitutional and interstitial Hf and Zr impurities in c-Ge, looking for neutral defects. We find that (i) substitutional Zr and Hf defects are energetically more favorable than interstitial defects; (ii) under oxygen-rich conditions, neither Zr nor Hf migration towards the channel is likely to occur; (iii) either under Hf- or Zr-rich conditions it is very likely, particularly for Zr, that defects will be incorporated in the channel.

Selective Wrapping and Supramolecular Structures of Polyfluorene-Carbon Nanotube Hybrids

We report on the photophysical properties of single-walled carbon nanotube (SWNT) suspensions In toluene solutions of poly[9,9-dioctylfluorenyl-2,7-diyl](PFO). Steady-state and time-resolved photoluminescence spectroscopy in the near-infrared and visible spectral regions are used to study the interaction of the dispersed SWNTs with the wrapped polymer. Molecular dynamics simulations of the PFO-SWNT hybrids in toluene were carried out to evaluate the energetics of different wrapping geometries. The simulated fluorescence spectra in the visible region were obtained by the quantum chemical ZINDO-CI method, using a sampling of structures obtained from the dynamics trajectories. The tested schemes consider polymer chains aligned along the nanotube axis, where chirality has a minimal effect, or forming helical structures, where a preference for high chiral angles is evidenced. Moreover, toluene affects the polymer structure favoring the helical conformation. Simulations show that the most stable hybrid system is the PFO-wrapped (8,6) nanotube, in agreement with the experimentally observed selectivity.

Computer simulations of copper and gold nanowires and single-wall nanowires

Copper and gold nanowires under tension evolve to form linear atomic chains (LACs), and the study and understanding of this evolution is an important subject for the development of nanocontacts. Here we study the differences and similarities between copper and gold nanowires (NWs) under stress along the [111] crystallographic direction until their rupture using tight-binding molecular dynamics. In both metals, the first significant rearrangement occurs due to one inside atom that goes to the NW` surface. In an attempt to better understand this effect, for both metals we also consider hollow NW`s where the inside atoms were excluded after the initial relaxation to create single-wall NW`s (SWNWs). The dynamical evolution of these SWNWs provides insight on the formation of the constriction that evolves to form LACs. Studying the calculated forces supported by the NW`s we show that SWNWs can sustain larger forces before the first major rearrangement in the copper and gold when compared to the original NW`s.

Ground-state energy of a classical artificial molecule

We study the ground-state energy of a classical artificial molecule formed by two-dimensional clusters (artificial atoms) of N/2 charged particles separated by a distance d. For the small molecules of N = 2 and 4, we obtain analytical expressions for this energy. For the larger ones, we calculate the ground-state energy using molecular dynamics simulation for N up to 128. From our numerical results, we are able to find out a function to approximate the ground-state energy of the molecules covering the range from atoms to molecules for any inter-atom distance d and for particle number from N = 8 to 128 within a difference less than one percent from the MD data.

Adsorption Behavior of 5-Fluorouracil on Au(111): An In Situ STM Study

The biological effects of chemical substitution of DNA bases triggered several investigations of their physicochemical properties This paper studies the adsorption behavior of a halogenated uracil, 5-fluorouracil (5FU). at the electrochemical interface of Au(111) and sulfuric acid solution. Upon modulation of the electric field across the interface, four distinct phases could be inferred by means of cyclic voltammetry (CV) At negative potentials relative to the SCE electrode, limited by the threshold of hydrogen evolution, no molecular species could be detected by scanning tunneling microscopy (STM) at the reconstructed Au(111)-(23 x root 3) surface, indicating that any physisorbed molecules are randomly distributed Incursion into more positive potentials increases the surface population but doer not form any two-dimensional (2D) physisorbed ordered structure Instead, we observed metastable structures that are only detectable. on surfaces with high defect density At sufficiently high positive potentials. limited by gold oxidation, the molecules are chemisorbed in a (3 x 2 root 3) ordered structure. with the aromatic ring perpendicular to the surface We report the densest chemisorbed monolayer for pyrimidine-derivative molecules (area per molecule 0 14 +/- 0 04 nm(2)). A comparison of the adsorption behavior of uracil derivatives has been made based on recent results of chemical substitution and solvent effects. We propose that pi-stacking is enhanced when halogens are incorporated in the uracil structure, in a similar fashion to what is observed in then crystal structure