Interdiffusion of small molecules into a glassy polymer film via coarse-grained molecular dynamics simulations

All-atom molecular dynamics simulation of solvent diffusion in an unentangled polystyrene film

The diffusion behavior of low molecular weight solvents within an unentangled polystyrene film below and above its glass transition temperature is investigated. The diffusion behavior in the glassy state exhibits a distinct behavior known as case II or class II diffusion, noticeably diverging from conventional Fickian diffusion observed above the glass transition temperature of the polymer film. In the context of case II diffusion, the primary experimental observation entails the emergence of a well-defined concentration front moving at a constant speed, delineating a swollen, rubbery region from a glassy region within the polymer system. Despite the prevalence of this phenomenon in experimental settings, simulating case II diffusion has posed a significant challenge, primarily due to the computationally intensive nature of the diffusion process. To address this, we have developed an all-atom molecular dynamics simulation approach for the observation of case II diffusion in glassy polymers. This method aims to unravel the intricacies of the diffusion process, providing valuable insights into the dynamic interactions between solvents and the polymer matrix.

Diffusion Coefficients of Variable-Size Amphiphilic Additives in a Glass-Forming Polyethylene Matrix

Diffusion of additives in polymers is an important issue in the plastics industry since migratory-type molecules are widely used to tune the properties of polymeric composites. Predicting the diffusional behavior of new additives can minimize the need for repetitive experiments. This work presents molecular dynamics simulations at the microsecond time scale and uses the MARTINI force field to estimate self-diffusion coefficients, D, of six monounsaturated amides and their analogs carboxylic acids in polyethylene matrices (PE, MW = 5600 Da). The results are strongly influenced by the glass-forming properties of the PE matrix, which we characterize by three distinct temperatures. The metastability region (T < 325 K), the glass transition temperature (Tg = 256-260 K), and the end of the transition (T ≅ 200 K). Self-diffusion mechanisms are inferred from the results of the dependence of D on the molecular mass of the additive, observing a Rouse-like behavior at high temperatures and deviations from it within the metastability region of the matrix. Interestingly, D values are nonsensitive to the nature of the considered polar head for additives of similar size. The temperature-dependent behavior of D follows, at fixed additive size, a linear Arrhenius pattern at high temperatures and a super Arrhenius trend at lower temperatures, which is well represented with a power law equation as predicted by the Mode Coupling Theory (MCT). We offer a conceptual explanation for the observed super-Arrhenius behavior. This explanation draws on Truhlar and Kohen's interpretation of the available energies at both the initial and the transition states along the diffusion pathway. The matrix's mobility significantly affects solute self-diffusion, yielding equal activation enthalpies for the Arrhenius region or the same power law parameters for the super-Arrhenius regime. Finally, we establish a one-to-one time-equivalence of the self-diffusion processes between CG and all-atom systems for the largest additives and the PE matrix in the high-temperature regime.

The Martini Model in Materials Science

The Martini model, a coarse-grained force field initially developed with biomolecular simulations in mind, has found an increasing number of applications in the field of soft materials science. The model's underlying building block principle does not pose restrictions on its application beyond biomolecular systems. Here, the main applications to date of the Martini model in materials science are highlighted, and a perspective for the future developments in this field is given, particularly in light of recent developments such as the new version of the model, Martini 3.

Polyelectrolyte stiffness on gold nanorods mediates cell membrane damage

Charge and surface chemistry of gold nanorods (AuNRs) are often considered the predictive factors for cell membrane damage. Unfortunately, extensive research on AuNR passivated with polyelectrolyte (PE) ligand shell (AuNR-PE) has hitherto left a vital knowledge gap between the mechanical stability of the ligand shell and the cytotoxicity of AuNR-PEs. Here, the agreement between unbiased coarse-grained molecular dynamics (CGMD) simulation and empirical outcomes on hemolysis of red blood cells by AuNR-PEs demonstrates for the first time, a direct impact of the mechanical stability of the PE shell passivating the AuNRs on the lipid membrane rupture. Such mechanical stability is ultimately modulated by the rigidity of the PE components. The CGMD simulation results also reveal the mechanism where the PE chain adsorbs near the surface of the lipid bilayer without penetrating the hydrophobic core of the bilayer, which allows the hydrophobic AuNR core to be in direct contact with the hydrophobic interior of the lipid bilayer, thereby perforating the lipid membrane to induce membrane damage.

Diffusion Behavior of Drug Molecules in Acrylic Pressure-Sensitive Adhesive

Acrylic pressure-sensitive adhesive (PSA) is widely used in transdermal drug delivery systems, while the diffusion behavior of drug molecules in PSA is of great importance. In this paper, PSAs with different cross-link densities were prepared by adjusting the ratio of cross-linkers. The effects of cross-link density and temperature on the diffusion of drugs in PSA were investigated by Fourier transform infrared attenuated total reflectance and molecular dynamics simulation. The consistency between the experimental and simulation results demonstrated that molecular dynamics simulation could be used to predict the diffusion behavior of drugs in PSA. The results showed that free volume and the wriggling of polymer chains are positively related to the diffusion coefficient of drug molecules, while hydrogen bonds hinder drug diffusion.