J Phys Chem B
September 2024
The test particle insertion method is used to study the solubility of oxygen in two commonly used polymers: polyethylene (PE) and isotactic polypropylene (iPP). Amorphous samples for both polymers were prepared by means of Monte Carlo and molecular dynamics simulations, and the oxygen solubility was measured across different temperatures. The solubility-temperature dependence for iPP proved to be nonmonotonic due to the interplay between binding and reorganizational enthalpy, while for PE, it appeared to be monotonic based on the available data in the studied temperature range.
View Article and Find Full Text PDFAmphiphilic block copolymer and lipids can be assembled into hybrid vesicles (HVs), which are an alternative to liposomes and polymersomes. Block copolymers that have either poly(sitostryl methacrylate) or statistical copolymers of sitosteryl methacrylate and butyl methacrylate as the hydrophobic part and a poly(carboxyethyl acrylate) hydrophilic segment are synthesized and characterized. These block copolymers assemble into small HVs with soybean L-α-phosphatidylcholine (soyPC), confirmed by electron microscopy and small-angle X-ray scattering.
View Article and Find Full Text PDFThis research focuses on modeling heat transfer in heterogeneous media composed of stacked spheres of paraffin as a perspective polymeric phase-change material. The main goal is to study the requirements of the numerical scheme to correctly predict the thermal conductivity in a periodic system composed of an indefinitely repeated configuration of spherical particles subjected to a temperature gradient. Based on OpenFOAM, a simulation platform is created with which the resolution requirements for accurate heat transfer predictions were inferred systematically.
View Article and Find Full Text PDFParaffin waxes are a promising material for heat storage with high energy density. Their low thermal conductivity, which limits the speed of charging and discharging in heat buffers, was previously shown to be improved by adding graphene nanofillers. In the present study, using molecular dynamics simulations, the segregation by molecular weight of polydisperse paraffin near graphene flakes is investigated.
View Article and Find Full Text PDFA molecular-level insight into phase transformations is in great demand for many molecular systems. It can be gained through computer simulations in which cooling is applied to a system at a constant rate. However, the impact of the cooling rate on the crystallization process is largely unknown.
View Article and Find Full Text PDFWe study the effect of solvent-free annealing and explicit solvent evaporation protocols in classical molecular dynamics simulations on the interface properties of a blend of a diketopyrrolopyrrole (DPP) polymer with conjugated substituents (DPP2Py2T) and PCBM[60]. We specifically analyze the intramolecular segmental mobility of the different polymer building blocks as well as intermolecular radial and angular distribution functions between donor and acceptor. The annealing simulations reveal an increase of the glass-transition temperature of 45 K in the polymer-fullerene blend compared to that of pure DPP2Py2T.
View Article and Find Full Text PDFOne of the important parameters in water management of proton exchange membranes is the electro-osmotic drag (EOD) coefficient of water. The value of the EOD coefficient is difficult to justify, and available literature data on this for Nafion membranes show scattering from in experiments and simulations. Here, we use a classical all-atom model to compute the EOD coefficient and thermodynamic properties of water from molecular dynamics simulations for temperatures between 330 and 420 K, and for different water contents between λ = 5 and λ = 20.
View Article and Find Full Text PDFUnderstanding the process-property relations of helical polymers using molecular simulations has been an attractive research field over the years. Specifically, isotactic polypropylene still remains a challenge for current computational experimentation, as it exhibits phenomena such as crystallization that emerge on large spatial and temporal scales. Coarse-graining is an efficient technique for approaching such phenomena, although previous coarse-grained models lack in preserving important atomistic and structural details.
View Article and Find Full Text PDFWe develop an all-atom force field for a series of diketopyrrolopyrrole polymers with two aromatic pyridine substituents and a variable number of π-conjugated thiophene units in the backbone (DPP2PyT), used as donor materials in organic photovoltaic devices. Available intrafragment parameterizations of the individual fragment building blocks are combined with interfragment bonded and nonbonded parameters explicitly derived from density functional theory calculations. To validate the force field, we perform classical molecular dynamics simulations of single polymer chains with = 1, 2, 3 in good and bad solvents and of melts.
View Article and Find Full Text PDFPrevious experiments of water vapor penetration into polyelectrolyte membrane (PEM) thin films have indicated the influence of the water concentration gradient and polymer chemistry on the interface evolution, which will eventually affect the efficiency of the fuel cell operation. Moreover, PEMs of different side chains have shown differences in water cluster structure and diffusion. The evolution of the interface between water and polyelectrolyte membranes (PEMs), which are used in fuel cells and flow batteries, of three different side-chain lengths has been studied using dissipative particle dynamics (DPD) simulations.
View Article and Find Full Text PDFThrough coarse-grained molecular dynamics simulation of polymer-grafted nanoparticles (NPs) in a lamellar-forming diblock copolymer (BCP), we systematically study the effects of the grafting density (), the compatibility between the grafted chains and the A-block of BCPs (ε), and the NP number () on the distance () of the NPs from the interface by proposing novel characterization parameters of the orientation and distribution of the grafted chains. The NP gradually migrates away from the interface and into the A-block region with the increase of ε for all studied , while slightly returning toward the interface at high ε and great , which is the first observation of nonmonotonic migration at the molecular level. We ascribe the reason of this to the behavior of the grafted chains that are near the interface.
View Article and Find Full Text PDFThe conventional definition of asphaltenes is based on their solubility in toluene and their insolubility in heptane. We have utilized this definition to study the influence of partial charge parametrization on the aggregation behavior of asphaltenes using classical atomistic molecular dynamics simulations performed on the microsecond time scale. Under consideration here are toluene- and heptane-based systems with different partial charges parametrized using the general AMBER force field (GAFF).
View Article and Find Full Text PDFThe improvement of mechanical properties of polymer nanocomposites (PNCs) has been studied for many years, with the main focus on the structure of the nanofillers. Much less effort has been devoted to unraveling the factors controlling the structure of the grafted chains. Herein, through coarse-grained molecular-dynamics simulations, we have successfully fabricated an ideal, mechanically-interlocked composite structure composed of end-functionalized chains grafted to the nanoparticle surface forming rings and making the matrix chains thread through these rings.
View Article and Find Full Text PDFWe present the results of the atomistic molecular dynamics modeling of different protonation states of Nafion at varying hydration levels. Previous experiments have shown that the degree of deprotonation (DDP) of the sulfonic acid groups in a Nafion membrane varies significantly upon hydration. Our goal is to provide insights into the effects of variable protonation states and water content on the internal structure and vehicular transport inside the Nafion membrane.
View Article and Find Full Text PDFUsing fully-atomistic models, tens-microseconds-long molecular-dynamic modelling was carried out for the first time to simulate the kinetics of polyimides ordering induced by the presence of single-walled carbon nanotube (SWCNT) nanofillers. Three polyimides (PI) were considered with different dianhydride fragments, namely 3,3',4,4'-biphenyltetracarboxylic dianhydride (BPDA), 2,3',3,4'-biphenyltetracarboxylic dianhydride (aBPDA), and 3,3',4,4'-oxidiphthalic dianhydride (ODPA) and same diamine 1,4-bis[4-(4-aminophenoxy)phenoxy]benzene (diamine P3). Both crystallizable PI BPDA-P3 and two amorphous polyimides ODPA-P3 and aBPDA-P3 reinforced by SWCNTs were studied.
View Article and Find Full Text PDFThis reactive molecular dynamics study explores the salt concentration dependence of the viscoelastic and mechanical failure properties of a poly(propylene glycol)/LiPF-based solid polymer electrolyte (SPE) at a graphitic carbon electrode interface. To account for the finite-size effect of interface-confined SPE films, the properties of two distinct film thicknesses are compared with the respective bulk properties. Additionally, the effect of uniaxial compression in the interface-normal direction on free energy profiles of Li-ion SPE-desolvation is studied.
View Article and Find Full Text PDFNafion nanocomposites for energy-related applications are being used extensively because of the attractive properties such as enhanced water retention, low unwanted crossover of electrolytes, and high proton conductivity. We present the results of the molecular dynamics modeling of Nafion films confined between two walls (substrates) of different polymer-wall interaction strengths and of different separation distances to model Nafion nanocomposites. Our goal is to provide insights into the effects of varying hydrophilicity and volume fraction of fillers/nanoparticles on the internal structure and water transport inside the Nafion membrane.
View Article and Find Full Text PDFBound rubber plays a key role in the mechanical reinforcement of elastomer nanocomposites. In the present work, we reveal the formation mechanism of bound rubber in elastomer nanocomposites, using the coarse-grained molecular dynamics simulations. For the polymer-nanoparticle system, the "chain bridge" connected with neighboring nanoparticles forms, once the gap between two neighboring nanoparticles is less than the polymer size.
View Article and Find Full Text PDFThis paper deals with molecular-dynamics simulations of the mechanical properties of prestretched double network filled elastomers. To this end, we firstly validated the accuracy of this method, and affirmed that the produced stress-strain characteristics were qualitatively consistent with Lesser's experimental results on the prestretched tri-block copolymers with a competitive double network. Secondly, we investigated the effect of the crosslinking network ratio on the mechanical properties of the prestretched double network homopolymers under uniaxial tension.
View Article and Find Full Text PDFCrystallization of all-aromatic heterocyclic polymers typically results in an improvement of their thermo-mechanical properties. Nucleation agents may be used to promote crystallization, and it is well known that the incorporation of nanoparticles, and in particular carbon-based nanofillers, may induce or accelerate crystallization through nucleation. The present study addresses the structural properties of polyetherimide-based nanocomposites and the initial stages of polyetherimide crystallization as a result of single-walled carbon nanotube (SWCNT) incorporation.
View Article and Find Full Text PDFIn this paper we adopt molecular dynamics simulations to study the amphiphilic AB block copolymer (BCP) mediated nanoparticle (NP) dispersion in polymer nanocomposites (PNCs), with the A-block being compatible with the NPs and the B-block being miscible with the polymer matrix. The effects of the number and components of BCP, as well as the interaction strength between A-block and NPs on the spatial organization of NPs, are explored. We find that the increase of the fraction of the A-block brings different dispersion effect to NPs than that of B-block.
View Article and Find Full Text PDFUnited-atom molecular-dynamics computer simulations of atactic polystyrene (PS) were performed for the bulk and free-standing films of 2 nm-20 nm thickness, for both linear and cyclic polymers comprised of 80 monomers. Simulated volumetric glass-transition temperatures (T) show a strong dependence on the film thickness below 10 nm. The glass-transition temperature of linear PS is 13% lower than that of the bulk for 2.
View Article and Find Full Text PDFUsing short polymer chains and through molecular dynamics simulation, we designed a well-dispersed nanoparticle (NP) network, which was then incorporated with the polymer matrix. First, we examined the effects of the dual-end grafted chains flexibility and density on the spatial distribution of this particular polymer nanocomposites system. By changing the interaction strength between the matrix polymer chains and the dual-end grafted chains in the semi-interpenetrating network system (NP network), we analyzed the interpenetration state between the linear polymer matrix and the NP network via calculating the total interfacial interaction energy.
View Article and Find Full Text PDFIn this work we perform molecular-dynamics simulations, both on the coarse-grained and the chemistry-specific levels, to study the influence of morphology on the mechanical properties of polymer nanocomposites (PNCs) filled with uniform spherical nanoparticles (which means without chemical modification) and patchy spherical nanoparticles (with discrete, attractive interaction sites at prescribed locations on the particle surface). Through the coarse-grained model, the nonlinear decrease of the elastic modulus (G') and the maximum of the viscous modulus (G″) around the shear strain of 10% is clearly reproduced. By turning to the polybutadiene model, we examine the effect of the shear amplitude and the interaction strength among uniform NPs on the aggregation kinetics.
View Article and Find Full Text PDFConstant temperature-constant pressure (NpT) molecular-dynamics computer simulations have been carried out for the united-atom model of a non-crosslinked (1,4) cis-polyisoprene (PI) melt confined between two amorphous, fully coordinated silica surfaces. The Lennard-Jones 12-6 potential was implemented to describe the polymer-silica interactions. The thickness H of the produced PI-silica film has been varied in a wide range, 1 < H/R(g) < 8, where R(g) is the individual PI chain radius of gyration measured under the imposed confinement.
View Article and Find Full Text PDF