Search

This computational study illustrates how modifying the end groups of polymer chains offers a viable strategy to control polymer crystal nucleation in entangled melts. More specifically, altering ch...

Spatial distributions are presented that quantitatively capture how polymer properties (e.g., segment alignment, density, and potential energy) vary with distance from nascent polymer crystals (nuclei) in prototypical polyethylene melts. It is revealed that the spatial extent of nuclei and their interfaces is metric-dependent as is the extent to which nucleus interiors are solid-like. As distance from a nucleus increases, some properties, such as density, decay to melt-like behavior more rapidly than polymer segment alignment, indicating that a polymer nucleus resides in a nematic-like droplet. This nematic-like droplet region coincides with enhanced formation of ordered polymer segments that are not part of the nucleus. It is more favorable to find nonconstituent ordered polymer segments near a nucleus than in the surrounding metastable melt, pointing to the possibility of one nucleus inducing the formation of other nuclei. In this vein, there is also a second region of enhanced ordering that lies along the nematic director of a nucleus, but beyond its nematic droplet and fold regions. These results indicate that crystal stacking, a key characteristic of lamellae in semicrystalline polymeric materials, begins to emerge during the earliest stages of polymer crystallization (i.e., crystal nucleation). More generally, the findings of this study provide a conceptual bridge between polymer crystal nucleation under nonflow and flow conditions and are used to rationalize previous results.

The high anisotropy of polymer-based fibres confers them high tensile strength, but also makes them more vulnerable against non-uniform mechanical loads. This is even more important for Kevlar® fib...

Halogen bonding has emerged as an important noncovalent interaction in a myriad of applications, including drug design, supramolecular assembly, and catalysis. Current understanding of the halogen bond is informed by electronic structure calculations on isolated molecules and/or crystal structures that are not readily transferable to liquids and disordered phases. To address this issue, we present a first-principles simulation-based approach for quantifying halogen bonds in molecular systems rooted in an understanding of nuclei-nuclei and electron-nuclei spatial correlations. We then demonstrate how this approach can be used to quantify the structure and dynamics of halogen bonds in condensed phases, using solid and liquid molecular chlorine as prototypical examples with high concentrations of halogen bonds. We close with a discussion of how the knowledge generated by our first-principles approach may inform the development of classical empirical models, with a consistent representation of halogen bonding., Comment: 8 pages, 8 figures

Alumina is one of the most abundant minerals and has a wide range of industrial applications, with catalysis as one of the most important. Of particular relevance for catalysis is the structure of the mineral/water interface. In this work, water structure and sodium halide adsorption at the neutral α-alumina(0001)/water interface are investigated using molecular dynamics simulations. This work demonstrates the accuracy of the chosen model of the alumina/water interface and shows that high charge density monovalent ions, such as Na+ and F–, have a strong affinity for the interface due to the specific pattern of alumina surface OH groups, such that the adsorbed ions displace waters that are hydrogen-bonded to the surface in their absence. A significant portion of the driving force for anion adsorption arises from surface bound Na+, which reverse the intrinsic surface dipole field and drive the accumulation of halides at the interface. The resulting electrolytic interfacial structure reorients water molecule...

The intercalation of Co in NiFe layered double hydroxides (LDH) significantly improves the electrocatalytic performance of the LDH for the oxygen evolution reaction. The mechanism behind the improv...

Oxide–water interfaces are ubiquitous, with many applications in industry and the environment, yet there is a great deal of controversy over their properties and microscopic structure. This controversy stems, in part, from the unique H-bond networks formed at different surface terminations and mineral compositions. Density functional theory simulations of these interfaces require an accurate description of both the oxide mineral and water in diverse H-bond environments. Thus, herein we simulate the Al2O3(001)–H2O interface using the PBE, PBE-TS, RPBE, SCAN, and HSE06-TS functionals to determine how calculated interfacial properties depend on the choice of functional. We find that the structure of the first few layers of water at the surface is determined by electron correlation in a way that cannot be approximated using semiemipirical van der Waals corrections. Of the functionals investigated, we find that SCAN yields the most accurate interfacial structure, dynamics, and sum frequency generation spectrum...

Development of inexpensive and efficient photo- and electro-catalysts is vital for clean energy applications. Electronic and structural properties can be tuned by the introduction of defects to achieve the desirable electrocatalytic activity. Using first-principles molecular dynamics simulations, the structural, dynamical, and electronic properties of 2D borocarbonitride (h-BCN) sheets have been investigated, highlighting how anti-site defects in B and N doped graphene significantly influence the bandgap, and thereby open up new avenues to tune the chemical behavior of the 2D sheets. In the present work, all of the monolayers investigated display direct bandgaps, which reduce from 0.99 eV to 0.24 eV with increasing number of anti-site defects. The present results for the electronic structure and findings for bandgap engineering open up applications of BCN monolayers in optoelectronic devices and solar cells. The influence of the anti-site distribution of B and N atoms on the ultra-high hole/electron mobility and conductivity is discussed based on density functional theory coupled with the Boltzmann transport equation. The BCN defect monolayer is predicted to have carrier mobilities three times higher than that of the pristine sheet. The present results demonstrate that BN doped graphene monolayers are likely to be useful in the next-generation 2D field-effect transistors.

Path-integral ab initio molecular dynamics (PI-AIMD) calculations have been employed to probe the nature of chloride ion solvation in aqueous solution. Nuclear quantum effects (NQEs) are shown to weaken hydrogen bonding between the chloride anion and the solvation shell of water molecules. As a consequence, the disruptive effect of the anion on the solvent water structure is significantly reduced compared to what is found in the absence of NQEs. The chloride hydration structure obtained from PI-AIMD agrees well with information extracted from neutron scattering data. Inparticular, the observed satellite peak in the hydrogen-chloride-hydrogen triple angular distribution serves as a clear signature of NQEs. The present results suggest that NQEs are likely to play acrucial role in determining the structure of saline solutions., 6 pages, 4figures

Poly($p$-phenylene terephthalamide) (PPTA) exhibits van der Waals (vdW) bonding along its $a$ axis, hydrogen bonding along its $b$ axis, and covalent bonding along its $c$ axis. We explore the structural and mechanical properties of PPTA using density functional theory with various functionals including LDA, PBE, PBE+rVV10L, SCAN, and SCAN+rVV10, compared with available experiments. The hierarchy of nonempirical semilocal functionals LDA, PBE, and SCAN (not fitted to any multicenter bonded system) includes differing amounts of intermediate-range vdW interaction. rVV10 is the long-range vdW correction. (rVV10L differs from rVV10 only in the value of a range parameter.) Among the tested functionals, SCAN shows the best performance for the lattice parameters of PPTA along the two crystal directions involving vdW or hydrogen-bond interaction. The equilibrium lattice constants obtained by SCAN and PBE+rVV10L are closest to experimental data, while SCAN+rVV10 slightly overbinds the system. We study the mechanical response of PPTA by applying strain along three lattice directions. Due to the inclusion of vdW interaction, SCAN, PBE+rVV10L, and SCAN+rVV10 all exhibit correct bonding strain-energy curves. On the contrary, PBE strongly underestimates the vdW interaction needed to resist uniaxial stretching along the $a$ axis. The Young's modulus and yield strength of PPTA are computed and compared with previous results. The experimental values are much smaller than the computed ones, mainly due to the fact that the PPTA fiber samples used for measurements are mechanically weaker than the perfect molecular crystals considered in the simulations. Interestingly, when a compressive uniaxial stress of 25 Kbar is applied along the $b$ axis, a structural phase transition, in which the hydrogen bonds reform along one diagonal of the $ab$ rectangle, is predicted by SCAN.

The specificity of aqueous halide solvation is fundamental to a wide range of bulk and interfacial phenomena spanning from biology to materials science. Halide polarizability is thought to drive the ion specificity, and if so, it is essential to have an accurate description of the electronic properties of halide ions in water. To this end, the solvation of the chloride anion, Cl− has been reinvestigated with state-of-the-art density functional theory. Specifically, the PBE-D3, PBE0-D3, and SCAN functionals have been employed to probe the impact of correlation and exchange approximations. Anticipating the findings, adding exact exchange improves the electronic structure, but simultaneously significantly reduces the Cl− polarizability, resulting in an over-structured Cl–O radial distribution function (RDF) and longer water H-bond lifetimes to Cl−. SCAN does not yield as much improvement in the energetics of Cl− relative to bulk water, but does result in a smaller reduction of the polarizability and thus a less structured Cl–O RDF, which agrees better with experiment. Special consideration is therefore warranted in assessing the impact of exchange on the energy, charge density, and the charge density response when designing and testing hybrid functionals for aqueous halide solvation.

This study demonstrates that monodisperse entangled polymer melts crystallize via the formation of nanoscale nascent polymer crystals (i.e., nuclei) that exhibit substantial variability in terms of their constituent crystalline polymer chain segments (stems). More specifically, large-scale coarse-grain molecular simulations are used to quantify the evolution of stem length distributions and their properties during the formation of polymer nuclei in supercooled prototypical polyethylene melts. Stems can adopt a range of lengths within an individual nucleus (e.g., &sim, 1&ndash, 10 nm) while two nuclei of comparable size can have markedly different stem distributions. As such, the attainment of chemically monodisperse polymer specimens is not sufficient to achieve physical uniformity and consistency. Furthermore, stem length distributions and their evolution indicate that polymer crystal nucleation (i.e., the initial emergence of a nascent crystal) is phenomenologically distinct from crystal growth. These results highlight that the tailoring of polymeric materials requires strategies for controlling polymer crystal nucleation and growth at the nanoscale.

Many atomic liquids can form transient covalent bonds reminiscent of those in the corresponding solid states. These directional interactions dictate many important properties of the liquid state, necessitating a quantitative, atomic-scale understanding of bonding in these complex systems. A prototypical example is liquid silicon, wherein transient covalent bonds give rise to local tetrahedral order and consequent non-trivial effects on liquid state thermodynamics and dynamics. To further understand covalent bonding in liquid silicon, and similar liquids, we present an ab initio simulation-based approach for quantifying the structure and dynamics of covalent bonds in condensed phases. Through the examination of structural correlations among silicon nuclei and maximally localized Wannier function centers, we develop a geometric criterion for covalent bonds in liquid Si. We use this to monitor the dynamics of transient covalent bonding in the liquid state and estimate a covalent bond lifetime. We compare covalent bond dynamics to other processes in liquid Si and similar liquids and suggest experiments to measure the covalent bond lifetime., Comment: 6 pages, 4 figures

We use density functional theory molecular dynamics simulations to investigate the structure, dynamics, and vibrational sum frequency generation (vSFG) spectra at the Al2O3(0001)–H2O and Al2O3(1120)–H2O interfaces. We find that the differences in the vSFG spectra between the two interfaces can be explained by significantly weaker surface–water interactions at the (0001) vs (1120) interface. The weaker interactions at the (0001) surface are caused by the flat surface plane and high density of OH groups, leading to a decoupling of the vibrational modes of the surface OH groups and H2O molecules. The (0001) vSFG spectrum thus displays two well-separated peaks at the near-neutral pH, in contrast to the vSFG spectrum of the corrugated (1120) interface, which has stronger surface–water interactions and thereby a narrower band in the vSFG spectrum with closely spaced peaks. By simulating the interfaces with both the Perdew–Burke–Ernzerhof (PBE)–Tkatchenko–Scheffler and revised PBE (RPBE) functionals, we find ...

Solid, liquid and alloyed phases of gallium play a role in a variety of important technological applications. While many of the gallium phases involved in these applications are metallic, some have...

We report the observation of a sizable photostrictive effect of 5.7% with fast, submillisecond response times, arising from a light-induced lattice dilation of a molecular nanosheet, composed of the molecular charge-transfer compound dibenzotetrathiafulvalene (DBTTF) and C60. An interfacial self-assembly approach is introduced for the thickness-controlled growth of the thin films. From photoabsorption measurements, molecular simulations, and electronic structure calculations, we suggest that photostriction within these films arises from a transformation in the molecular structure of constituent molecules upon photoinduced charge transfer, as well as the accommodation of free charge carriers within the material. Additionally, we find that the photostrictive properties of the nanosheets are thickness-dependent, a phenomenon that we suggest arises from surface-induced conformational disorder in the molecular components of the film. Moreover, because of the molecular structure in the films, which results largely from interactions between the constituent π-systems and the sulfur atoms of DBTTF, the optoelectronic properties are found to be anisotropic. This work enables the fabrication of 2D molecular charge-transfer nanosheets with tunable thicknesses and properties, suitable for a wide range of applications in flexible electronic technologies.

Supramolecular orientational memory (SOM) provides a route to otherwise inaccessible nanoscale architectures for certain molecules. In these privileged cases, columnar domains organized from self-assembling dendrons undergo reorientation during heating to, and subsequent cooling from, a 3D phase composed of “spheres”, such as a body-centered cubic phase or a Pmn cubic phase, known also as Frank-Kasper A15. The directions of the reoriented columns preserve key interactions from the preceding cubic phase. However, SOM was observed so far in a very limited number of assemblies. The molecular determinants enabling SOM, and its generality, remain poorly understood. Here we report the synthesis and structural and retrostructural analysis of a perylene bisimide (PBI) with two self-assembling benzyl ether dendrons, 3,5-G2-PBI, and compare its assemblies with those of a previously reported PBI, 3,4,5-G2-PBI, which exhibits SOM and has an additional minidendritic building block in its dendrons. The removal of this minidendron in 3,5-G2-PBI eliminates its ability to self-assemble into supramolecular spheres and organize into a cubic phase, thereby precluding 3,5-G2-PBI from exhibiting SOM. This finding demonstrates hierarchical transfer of structural information from primary structure to material function, analogous to the misfolding of proteins into toxic structures such as those implicated in Alzheimer's and Prion diseases. The concepts exemplified here provide new insights into the hierarchical basis for SOM and will aid in the translation of the SOM concept to a broader diversity of soft matter such as block copolymers and surfactants.

Hexagonal borocarbonitrides (BCN) are a class of 2D materials, which display excellent catalytic activity for water splitting. Here, we report analysis of thermal stability, phonons and thermal conductivity of BCN monolayers over a wide range of temperatures using classical molecular dynamics simulations. Our results show that in contrast to the case of graphene and boron nitride monolayers, the out-of-plane phonons in BCN monolayers induce an asymmetry in the phonon density of states at all temperatures. Despite possessing lower thermal conducting properties compared to graphene and BN monolayers, the BCN nanosheets do not lose thermal conductivity as much as graphene and BN in the studied temperature range of 200-1000 K, and thus, the BCN nanosheets are suitable for thermal interface device applications over a wide range of temperatures. Besides their promising role in water splitting, the above results highlight the possibility of expanding the use of BCN 2D materials in thermal management applications and thermoelectrics., Comment: 6 pages, 4 figures

The Shinoda-DeVane-Klein (SDK) model is herein demonstrated to be a viable coarse-grain model for performing molecular simulations of polyethylene (PE), affording new opportunities to advance molecular-level, scientific understanding of PE materials and processes. Both structural and dynamical properties of entangled PE melts are captured by the SDK model, which also recovers important aspects of PE crystallization phenomenology. Importantly, the SDK model can be used to represent a variety of materials beyond PE and has a simple functional form, making it unique among coarse-grain PE models. This study expands the suite of tools for studying PE in silico and paves the way for future work probing PE and PE-based composites at the molecular level.

Few polymers are as well-known as PPTA, the main constituent of Kevlar\textregistered{} fibers. To achieve high mechanical strength, PPTA chains must be treated with sulfuric acid, which is removed after fibers are formed. However, simulations show that tiny clusters of sulfuric acid remain embedded deeply within the fibers. Their presence is likely to go undetected, but can severely affect the material's strength under the high-strain conditions of its intended use. If the process behind the strength of PPTA fibers is also directly responsible for their main weakness, a solvent-free process is a promising route toward stronger materials

We demonstrate that nascent polymer crystals (i.e., nuclei) are anisotropic entities, with neither spherical nor cylindrical geometry, in contrast to previous assumptions. In fact, cylindrical, spherical, and other high symmetry geometries are thermodynamically unfavorable. Moreover, post-critical transitions are necessary to achieve the lamellae that ultimately arise during the crystallization of semicrystalline polymers. We also highlight how inaccurate treatments of polymer nucleation can lead to substantial errors (e.g., orders of magnitude discrepancies in predicted nucleation rates). These insights are based on quantitative analysis of over four million crystal clusters from the crystallization of prototypical entangled polyethylene melts. New comprehensive bottom-up models are needed to capture polymer nucleation., Comment: 19 pages long with 5 figures

Memory effects, including shape, chirality, and liquid-crystallinity, have enabled macroscopic materials with novel functions. However, the generation of complex supramolecular nanosystems via memory effects has not yet been investigated. Here, we report a cyclotriveratrylene-crown (CTV) compound that self-assembles into supramolecular columns and spheres forming, respectively, hexagonal and cubic mesophases. Upon transition from one phase to the other, an epitaxial relationship holds, via an unprecedented supramolecular orientational memory effect. Specifically, the molecular orientation and columnar character of supramolecular packing is preserved in the cubic phase, providing an otherwise inaccessible structure comprising orthogonally oriented domains of supramolecular columns. The continuous columnar character of tetrahedrally distorted supramolecular spheres self-organized from the CTV derivative in the faces of the Pm3̅n lattice is the basis of this supramolecular orientational memory, which holds throughout cycling in temperature between the two phases. This concept is expected to be general for other combinations of periodic and quasiperiodic arrays generated from supramolecular spheres upon transition to supramolecular columns.

Halide perovskites form the foundation of an emerging class of materials for broad application in renewable and sustainable applications, including photocatalysis and solar energy harvesting. These materials exhibit beneficial photophysical properties, including bandgaps suitable for solar energy harvesting and efficient charge screening that underlies efficient charge carrier separation and resistance to defects. For organic–inorganic hybrid perovskites, these benefits are thought to arise, in part, from dipolar molecular cations that can reorient in response to charged particles and defects. In this work, we provide a similar perspective for inorganic metal halide perovskites, which do not contain molecular species with permanent dipoles. We discuss how lone pair electrons lead to dipolar ions that exhibit dynamics in analogy with traditional molecular plastic crystals and hybrid perovskites. We argue that further understanding these electronic plastic crystal motions with first principles simulations and synchrotron scattering can help create a basic understanding of photophysical properties of metal halide perovskites and inform the design of advanced functional materials.

Ultrahigh-molecular-weight polyethylene (UHMWPE) is of great interest as a next-generation body armor material because of its superior mechanical properties. However, such unique properties depend critically on its microscopic structure characteristics, including the degree of crystallinity, chain alignment, and morphology. Here, we present a highly aligned UHMWPE and its composite sheets containing uniformly dispersed boron nitride (BN) nanosheets. The dispersion of BN nanosheets into the UHMWPE matrix increases its mechanical properties over a broad temperature range. Experiments and simulation confirm that the alignment of chain segments in the composite matrix increases with temperature, leading to an improvement in mechanical properties at high temperature. Together with the large thermal conductivity of UHMWPE and BN, our findings serve to expand the application spectrum of highly aligned polymer nanocomposite materials for ballistic panels and body armor over a broad range of temperatures.

The non-selective cation channel TRPV1 is responsible for transducing noxious stimuli into action potentials propagating through peripheral nerves. It is activated by temperatures greater than 43 °C, while remaining completely non-conductive at temperatures lower than this threshold. The origin of this sharp response, which makes TRPV1 a biological temperature sensor, is not understood. Here we used molecular dynamics simulations and free energy calculations to characterize the molecular determinants of the transition between non-conductive and conductive states. We found that hydration of the pore and thus ion permeation depends critically on the polar character of its molecular surface: in this narrow hydrophobic enclosure, the motion of a polar side-chain is sufficient to stabilize either the dry or wet state. The conformation of this side-chain is in turn coupled to the hydration state of four peripheral cavities, which undergo a dewetting transition at the activation temperature.

In this work, the solvation and electronic structure of the aqueous chloride ion solution was investigated using Density Functional Theory (DFT) based \textit{ab initio} molecular dynamics (AIMD). From an analysis of radial distribution functions, coordination numbers, and solvation structures, we found that exact exchange ($E_{\rm xx}$) and non-local van der Waals (vdW) interactions effectively \textit{weaken} the interactions between the Cl$^-$ ion and the first solvation shell. With a Cl-O coordination number in excellent agreement with experiment, we found that most configurations generated with vdW-inclusive hybrid DFT exhibit 6-fold coordinated distorted trigonal prism structures, which is indicative of a significantly disordered first solvation shell. By performing a series of band structure calculations on configurations generated from AIMD simulations with varying DFT potentials, we found that the solvated ion orbital energy levels (unlike the band structure of liquid water) strongly depend on the underlying molecular structures. In addition, these orbital energy levels were also significantly affected by the DFT functional employed for the electronic structure; as the fraction of $E_{\rm xx}$ was increased, the gap between the highest occupied molecular orbital of Cl$^-$ and the valence band maximum of liquid water steadily increased towards the experimental value., Comment: 13 pages, 8 figures

Oxygen K-edge X-ray absorption spectra of liquid water are computed based on the configurations from advanced ab initio molecular dynamics simulations, as well as an electron excitation theory from the GW method. One one hand, the molecular structures of liquid water are accurately predicted by including both van der Waals interactions and hybrid functional (PBE0). On the other hand, the dynamic screening effects on electron excitation are approximately described by the recently developed enhanced static Coulomb hole and screened exchange approximation by Kang and Hybertsen [Phys. Rev. B 82, 195108 (2010)]. The resulting spectra of liquid water are in better quantitative agreement with the experimental spectra due to the softened hydrogen bonds and the slightly broadened spectra originating from the better screening model., 10 pages, 5 figures, accepted by Phys. Rev. B

A three-component system of Janus dendrimers (JDs) including hydrogenated, fluorinated, and hybrid hydrogenated–fluorinated JDs are reported to coassemble by film hydration at specific ratios into an unprecedented class of supramolecular Janus particles (JPs) denoted Janus dendrimersomes (JDSs). They consist of a dumbbell-shaped structure composed of an onion-like hydrogenated vesicle and an onion-like fluorinated vesicle tethered together. The synthesis of dye-tagged analogs of each JD component enabled characterization of JDS architectures with confocal fluorescence microscopy. Additionally, a simple injection method was used to prepare submicron JDSs, which were imaged with cryogenic transmission electron microscopy (cryo-TEM). As reported previously, different ratios of the same three-component system yielded a variety of structures including homogenous onion-like vesicles, core-shell structures, and completely self-sorted hydrogenated and fluorinated vesicles. Taken together with the JDSs reported herein, a self-sorting pathway is revealed as a function of the relative concentration of the hybrid JD, which may serve to stabilize the interface between hydrogenated and fluorinated bilayers. The fission-like pathway suggests the possibility of fusion and fission processes in biological systems that do not require the assistance of proteins but instead may result from alterations in the ratios of membrane composition.

Two-dimensional van der Waals heterostructures are of considerable interest for the next generation nanoelectronics because of their unique interlayer coupling and optoelectronic properties. Here, we report a modified Langmuir–Blodgett method to organize two-dimensional molecular charge transfer crystals into arbitrarily and vertically stacked heterostructures, consisting of bis(ethylenedithio)tetrathiafulvalene (BEDT–TTF)/C60 and poly(3-dodecylthiophene-2,5-diyl) (P3DDT)/C60 nanosheets. A strong and anisotropic interfacial coupling between the charge transfer pairs is demonstrated. The van der Waals heterostructures exhibit pressure dependent sensitivity with a high piezoresistance coefficient of −4.4 × 10−6 Pa−1, and conductance and capacitance tunable by external stimuli (ferroelectric field and magnetic field). Density functional theory calculations confirm charge transfer between the n-orbitals of the S atoms in BEDT–TTF of the BEDT–TTF/C60 layer and the π* orbitals of C atoms in C60 of the P3DDT/C60 layer contribute to the inter-complex CT. The two-dimensional molecular van der Waals heterostructures with tunable optical–electronic–magnetic coupling properties are promising for flexible electronic applications., Two-dimensional van der Waals heterostructures are of interest due to their unique interlayer coupling and optoelectronic properties. Here authors develop a Langmuir-Blodgett method to organize charge transfer molecular heterostructures with externally tunable conductance and capacitance and broadband photoresponse.

In addition to its technological relevance, silicon poses a challenge for first principles simulations because it undergoes a semiconductor-to-metal transition upon melting. Moreover, the resulting metallic liquid contains a mixture of metallic and covalent bonding. This coexistence of fundamentally different interactions is difficult to describe within approximate density functional methods, which oftentimes cannot accurately describe these two extremes simultaneously. We report an investigation of the structure, dynamics, and thermodynamics of liquid silicon using ab initio molecular dynamics simulations with three density functional approximations: the local density approximation, the Perdew-Burke-Ernzerhof generalized gradient approximation, and the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation. We demonstrate that SCAN describes this liquid with better accuracy than the other often-used functionals because it can simultaneously capture covalent and metallic bonding with similar high accuracy.

A constitutional isomeric library synthesized by a modular approach has been used to discover six amphiphilic Janus dendrimer primary structures, which self-assemble into uniform onion-like vesicles with predictable dimensions and number of internal bilayers. These vesicles, denoted onion-like dendrimersomes, are assembled by simple injection of a solution of Janus dendrimer in a water-miscible solvent into water or buffer. These dendrimersomes provide mimics of double-bilayer and multibilayer biological membranes with dimensions and number of bilayers predicted by the Janus compound concentration in water. The simple injection method of preparation is accessible without any special equipment, generating uniform vesicles, and thus provides a promising tool for fundamental studies as well as technological applications in nanomedicine and other fields.

This study reveals important features of polymer crystal formation at high-driving forces in entangled polymer melts based on simulations of polyethylene. First and in contrast to small-molecule crystallization, the heat released during polymer crystallization does not appreciably influence structural details of early-stage, crystalline clusters (crystal nuclei). Second, early-stage polymer crystallization (crystal nucleation) can occur without substantial chain-level relaxation and conformational changes. This study's results indicate that local structures and environments guide crystal nucleation in entangled polymer melts under high-driving force conditions. Given that such conditions are often used to process polyethylene, local structures and the separation of time scales associated with crystallization and chain-level processes are anticipated to be of substantial importance to processing strategies. This study highlights new research directions for understanding polymer crystallization.

Polymer micelles and vesicles form upon hydration of amphiphilic block copolymers in dilute aqueous solution. Present challenges with these self-assemblies include understanding how molecular structure and polydispersity determine nano- and mesoscopic shape and properties such as flexibility. With charged copolymers, divalent ionic ligands can rigidify polymer vesicle membranes and also induce microphase separation into domains of weak polyelectrolyte gels. In this review, we focus on the underlying physical and molecular questions concerning copolymer assembly and associated challenges and implications for nano-delivery materials. We also highlight molecular simulation techniques that can be used to investigate properties of assemblies, such as curvature, patterning, and other ionic effects in functional polymeric membranes and micelles.

Nanoparticle shape can improve drug delivery, based in part on recent fi ndings that fl exible, worm-like nanocarriers (Worms) increase the amount of drug delivered to tumors and shrink the tumors more effectively than spherical micelles (Spheres). Here, all-atom molecular dynamics (MD) simulations are used to build a rational coarse grain (rCG) model that helps clarify shape-dependent effects in delivery of the widely used anticancer drug Taxol by block copolymer micelles. Potentials for rCG-MD were developed to examine the partitioning of this hydrophobic-aromatic drug into Worms and Spheres that selfassemble in water from poly(ethyleneglycol)-poly(caprolactone) (PEG-PCL), a weakly segregating amphiphile. PCL is a biodegradable, hydrophobic polymer widely used in biomaterials and accurately modeled here. Thermodynamic integration of the force to pull a single Taxol molecule from the micelles into solvent shows that twice as much drug loads into Worms than Spheres, fully consistent with experiments. Diffusivity of drug in the hydrated PEG corona is surprisingly slow compared to that in the core, indicative of strong but transient drug-polymer interactions. The distinctly distended corona of the Worms enhances such interactions and refl ects the same balance of molecular forces that underlie an experimentally-validated phase diagram for simulated Spheres, Worms, and Bilayers. Moreover, with realistic drug loadings in micro-second simulations, Taxol is seen to draw PEG chains into the PCL core, dispersing the drug while localizing it near the interface—thus providing a molecular explanation for a measurable burst release of drug as well as the enhanced delivery seen with Worms. Cancer is diagnosed today in about 50% of the population over the average lifetime of a person, with nearly all therapies pushed from bench to bedside based solely on trial-and-error experimentation, particularly when it comes to formulation. Current treatments generally include chemotherapy, and many of the top anticancer drugs in the clinic have a hydrophobic and/or aromatic character that allows them to permeate cell

We present the results of coarse grained molecular dynamics simulation using a charge free model that is able to capture different regions of the morphological phase diagram of charged diblock copolymers. Specifically, we were able to reproduce many phases of the poly(acrylic acid)-(1,4)-polybutadiene (PAA-PBA) diblock copolymer, Ca(2+) and water systems as a function of pH and calcium concentration with short-range LJ type potentials. The morphologies observed range from bilayers to cylinders to spherical micelles. Such polyanionic/cationic amphiphiles in water typically present multiple challenges for molecular simulations, particularly due to the many charge interactions that are long ranged and computationally costly. Further, it is precisely these interactions that are thought to modulate large amphiphile assemblies of interest such as lipid rafts. However, our model is able to reproduce different morphologies due to pH and with or without the addition of Ca(2+) as well as the lateral phase segregation and the domain registration observed in neutral and charged diblock copolymer mixtures. The results suggest that the overall effect of charges is a local structural rearrangement that renormalizes the steric repulsion between the headgroups. This simple model, which is devoid of charges, is able to reproduce the complex phase diagram and can be used to investigate collective phenomena in these charged systems such as domain formation and registration or colocalization of lipid rafts across bilayer leaflets.

The transfer of various nano-scale fullerenes into lipid bilayers has been studied using all-atom (AA) and coarse-grained (CG) molecular dynamics (MD) simulations. The free energy change, when C60, C180, and C540 fullerenes are transferred from water to the interior of a lipid [dioleoylphosphatidylcholine (DOPC)] bilayer, has been calculated. Upon entering the lipid bilayer, the largest (2.4 nm diameter) fullerene causes local distortions in the bilayer surface, which were previously observed in carbon nanotube simulations. These local distortions, however, do not lead to any free energy barriers to bilayer entry. The free energy profiles confirm spontaneous absorption of all three fullerenes. Qualitative agreement was observed when comparing fullerene partitioning in water/bilayer systems to water–hexane systems. In contrast to these nonspecific single fullerene properties, extensive CG-MD simulations of fullerene rich lipid bilayers reveal substantial impact of fullerene-size on the bilayer stability. While previous CG-MD simulations indicated that bilayer bound C60 aggregates have little effect on the bilayer structure, the present MD simulations indicate that C540 aggregation has drastic effects. Specifically, the observed destabilization likely has implications for understanding the cytotoxic mechanisms of nano-carbon particles upon uptake by biological cells.

Janus Drug Delivery Vehicle Efficient drug delivery vehicles need to be produced in a limited size range and with uniform size distribution. The self-assembly of traditional small-molecule and polymeric amphiphiles has led to the production of micelles, liposomes, polymeric micelles, and polymersomes for use in drug delivery applications. Now, Percec et al. (p. 1009 ) describe the self-assembly of Janus-type (i.e., two-headed) dendrimers to produce monodisperse supramolecular constructs, termed “dendrimersomes,” and other complex architectures. The structures, which showed long-term stability as well as very narrow size distributions, were easily produced by the injection of an ethanolic solution of the dendrimer into water. The dendrimersomes could be loaded with the anticancer drug doxorubicin and exhibited controlled drug release with changing pH.

Combining single-walled carbon nanotubes (CNT) with biological molecules provides a route to novel nanoscale materials with many promising applications in nanotechnology and nanomedicine. Recent experiments show that CNTs covalently functionalized with the coxsackie-adenovirus receptor (CAR) serve as biosensors capable of specifically recognizing Knob proteins from the adenovirus capsid. These experiments suggest that CAR retains its biologically active form when bound to CNT, but a detailed understanding of the structural changes that occur within CAR after CNT attachment is lacking. To address this, we have performed all-atom classical molecular dynamics (MD) simulations of CAR and the CAR-Knob complex in aqueous solution alone and also when covalently linked to CNT. The MD results show that the CNT damps structural fluctuations in CAR and reduces the internal mobility of the protein. However, CNT induces very little structural deformation and does not affect CAR's ability to specifically bind Knob. This MD study verifies that CAR retains its biological functionality when attached to CNT and provides a computational approach to rationalize nanobiosensing devices.

We demonstrate a versatile class of nanoscale chemical sensors based on single-stranded DNA (ssDNA) for chemical recognition and single-walled carbon nanotube field effect transistors (SWNT FETs) for electronic read-out. SWNT FETs with a nanoscale coating of ssDNA respond to vapours that cause no detectable conductivity change in bare devices. The gases tested are methanol, trimethylamine, propionic acid, dimethylmethylphosphonate and dinitrotoluene. Sensor responses differ in sign and magnitude for different gases and can be tuned by choice of the ssDNA base sequence. Sensors respond and recover rapidly (seconds), and the sensor surface is self-regenerating. Preliminary results of all-atom molecular dynamics simulations are consistent with experiment.

A DPD model of PEO-based block copolymer vesicles in water is developed by introducing a new density based coarse graining and by using experimental data for interfacial tension. Simulated as a membrane patch, the DPD model is in excellent agreement with experimental data for both the area expansion modulus and the scaling of hydrophobic core thickness with molecular weight. Rupture simulations of polymer vesicles, or "polymersomes", are presented to illustrate the system sizes feasible with DPD. The results should provide guidance for theoretical derivations of scaling laws and also illustrate how spherical polymer vesicles might be studied in simulation.

Block-copolymer amphiphiles have been observed to assemble into vesicles and other morphologies long known for lipids but with remarkably different properties. Coarse-grain molecular dynamics (CG-MD) is used herein to elaborate the structures and properties of diblock copolymer assemblies in water. By varying the hydrophilic/hydrophobic ratio of the copolymer in line with experiment, bilayer, cylindrical and spherical micelle morphologies spontaneously assemble. Varying the molecular weight (MW) with hydrophilic/hydrophobic ratio appropriate to a bilayer yields a hydrophobic core thickness that scales for large MW as a random coil polymer, in agreement with experiment. The extent of hydrophobic-segment overlap in the core increases nonlinearly with MW, indicative of chain entanglements and consistent with the dramatic decrease reported for lateral mobility in polymer vesicles. Calculated trends with MW as well as hydrophilic/hydrophobic ratio thus agree with experiment, demonstrating that CG-MD simulations provide a rational design tool for diblock copolymer assemblies.

This paper illustrates the use of computer simulation in understanding the interaction of nano-scale synthetic molecules with cellular membranes, which are themselves a few nanometres thick. Specifically, the interaction of a membrane-bound synthetic molecule with a model membrane has been studied using so-called coarse-grain (CG) molecular dynamics (MD). The reported CG-MD simulations have been carried out for a pore-promoting hydraphile molecule in a dimyristoylphosphatidylcholine (DMPC) lipid bilayer. Analysis of the CG-MD simulation provides insights into the nature of the interaction of the hydraphile with the lipid molecules and how an initially fully extended trans-membrane hydraphile adjusts its end-to-end distance to match the bilayer thickness. The equilibrium membrane-bound hydraphile conformation observed in the CG-MD simulation agrees well with deductions based on experimental observations. The success of the CG-MD approach in the present example suggests that the methodology could be used as a design tool in related nano-science and engineering applications.

A coarse-grain model for amphiphilic diblock copolymers is developed by fitting the required parameters to properties taken from all-atom molecular dynamics simulations and experimental data. Computations with the present coarse-grain model yield spontaneous self-assembly of a random system into membrane bilayers when the amphiphilic diblock copolymers have a lipid-like hydrophilic/hydrophobic ratio. The model semiquantitatively reproduces a number of experimental results that were not explicitly used in the parametrization. In particular, diblock polymers with the appropriate ratio of hydrophobic−hydrophilic segment lengths self-assemble into membranes whose hydrophobic thickness (determined from mass density profiles) and scaling with molecular weight are found to be in good agreement with the experiment.

The interaction of surfactant diblock poly(ethylene oxide)–poly(ethylethylene) copolymers (PEO–PEE) with a lipid bilayer of dimyristoylphosphatidylcholine has been studied by means of coarse-grain molecular dynamics simulations. The effect of the surfactants on the lipid bilayer was studied over a wide range of diblock copolymer concentrations. The simulations show that the hydrophilic PEO chains adopt different structures at low and high concentrations. In particular, the computed density profiles reveal that the PEO chains extend over a longer range from the bilayer surface, with increasing copolymer concentration. The simulated density profiles are in agreement with the scaling law predictions.

A computationally efficient coarse grain (CG) model designed to mimic the lipid molecule, dimyristoylphosphatidylcholine (DMPC) is used to study the self-assembly of a lamellar bilayer starting from a disordered configuration. The utility of the CG model is illustrated by examining structural and dynamical properties of a hydrated bilayer system containing 1024 lipid molecules. Comparisons with results for an all-atom model of DMPC suggest that the CG model is about four orders of magnitude less demanding of CPU time.

The effect of stretching a polymer sample containing a single trefoil knot has been studied by computer simulation molecular dynamics calculations. Under axial load that approximates a fiber extrusion process, the knot is found to nucleate crystallization of the sample, which occurs on the ns time scale. The extension of the strain field associated with the knot has been quantified.

The rupture properties of a polyethylene knot are studied by first-principles molecular dynamics calculations. While earlier works were conducted in isolated-molecule approximation, we presently report a study focused on a more realistic system where the bulk environment is taken into account through the inclusion of nearest-neighboring chains and a proper set of constraints. Our work shows that the nature of the break in real samples containing knots is essentially intramolecular and that the multimolecules model reproduces results very similar to those of both bulk polyethylene and single-chain systems, thus confirming the relevance of earlier findings.

A coarse grain model for phospholipids was systematically parametrized to mimic structural properties obtained from an atomistic simulation of a dimyristoylphosphatidylcholine bilayer. The model semiquantitatively reproduces the cross-sectional structure of a preassembled phospholipid bilayer obtained from an atomistic simulation; a property that was not directly fit. The model is sufficiently fast to permit the simulation of the self-assembly of the bilayer starting from a random configuration.

The entanglement of neighboring chains is a common topological defect found in polymers. The limiting stress and mechanism of rupture of two entangled polyethylene chains, subjected to external loading, has been investigated using first-principles molecular dynamics calculations. Our results show that intermolecular entanglement significantly weakens each polymer strand. Interchain bond−bond or atom−bond friction is the major driving force in the rupture process.