Study of the Mechanical and Thermal Properties of Poly(Lactic Acid) and Poly(Ethylene Glycol) Block Copolymer with Molecular Dynamics

Micromechanical Dilution of PLA/PETG-Glass/Iron Nanocomposites: A More Efficient Molecular Dynamics Approach

Polylactic acid (PLA) and poly(ethylene terephthalate glycol) (PETG) are popular thermoplastics used in additive manufacturing applications. The mechanical properties of PLA and PETG can be significantly improved by introducing fillers, such as glass and iron nanoparticles (NPs), into the polymer matrix. Molecular dynamics (MD) simulations with the reactive INTERFACE force field were used to predict the mechanical responses of neat PLA/PETG and PLA-glass/iron and PETG-glass/iron nanocomposites with relatively high loadings of glass/iron NPs. We found that the iron and glass NPs significantly increased the elastic moduli of the PLA matrix, while the PETG matrix exhibited modest increases in elastic moduli. This difference in reinforcement ability may be due to the slightly greater attraction between the glass/iron NP and PLA matrix. The NASA Multiscale Analysis Tool was used to predict the mechanical response across a range of volume percent glass/iron filler by using only the neat and highly loaded MD predictions as input. This provides a faster and more efficient approach than creating multiple MD models per volume percent per polymer/filler combination. To validate the micromechanics predictions, experimental samples incorporating hollow glass microspheres (MS) and carbonyl iron particles (CIP) into PLA/PETG were developed and tested for elastic modulus. The CIP produced a larger reinforcement in elastic modulus than the MS, with similar increases in elastic modulus between PLA/CIP and PETG/CIP at 7.77 vol % CIP. The micromechanics-based mechanical predictions compare excellently with the experimental values, validating the integrated micromechanical/MD simulation-based approach.

Coarse-Grained Molecular Dynamics Simulations of Paclitaxel-Loaded Polymeric Micelles

A continuous manufacturing technology based on coaxial turbulent jet in coflow was previously developed to produce paclitaxel-loaded polymeric micelles. Herein, coarse-grained molecular dynamics (CG-MD) simulations were implemented to better understand the effect of the material attributes (i.e., the drug-polymer ratio and the ethanol concentration) and process parameters (i.e., temperature) on the self-assembly process of polymeric micelles as well as to provide molecular details on micelle instability. An all-atom (AA) poly (ethylene glycol)-poly (lactic acid) (PEG-PLA) polymer model was developed as the reference for parameterizing a coarse-grained (CG) model, and the AA polymer model was further validated with experimental glass transition temperature (T_g). The model transferability was verified by comparing structural properties between the AA and CG models. The CG model was further validated with experimental data, including micelle particle size measurements and drug encapsulation efficiency. Furthermore, the encapsulation of paclitaxel into the polymeric micelles was included in the simulations, taking into consideration the interactions between the paclitaxel and the polymers. The results from various points of view demonstrated a strong dependence of the shape of the micelles on the drug encapsulation, with micelles transitioning from spherical to ellipsoidal structures with an increasing paclitaxel amount. Simulation data were also used to identify the critical aggregation number (i.e., the number of polymer and drug molecules required for transition from one shape to another). Improved micellar structural stability was found with a larger micellar size and less solvent accessibility. Lastly, an evaluation was performed on the micellar dissociation free energy using a steered molecular dynamics simulation over a range of temperatures and ethanol concentrations. These simulations revealed that at higher ethanol and temperature conditions, micelles become destabilized, resulting in greater paclitaxel release. The increased drug release was determined to originate from the solvation of the hydrophobic core, which promoted micellar swelling and an associated reduction in hydrophobic interactions, leading to a loosely packed micellar structure.

Improved Core Viscosity Achieved by PDLLA10kCo-Incorporation Promoted Drug Loading and Stability of mPEG2k-b-PDLLA2.4k Micelles

PURPOSE

This study aims to investigate the effect of poly(D, L-lactic acid)10K (PDLLA10K) incorporation on the drug loading and stability of poly(ethylene glycol)2K-block-poly(D, L-lactide)2.4K (mPEG2k-b-PDLLA2.4k) micelles. In addition, a suitable lyophilization protector was screened for this micelle to obtain favorable lyophilized products.

METHODS

The incorporation ratios of PDLLA10k were screened based on the particle size and drug loading. The dynamic stability, core viscosity, drug release, stability in albumin, and in vivo pharmacokinetic characteristics of PDLLA10k incorporated micelles were compared with the original micelles. In addition, the particle size variation was used as an indicator to screen the most suitable lyophilization protectant for the micelles. DSC, FTIR, XRD were used to illustrate the mechanism of the lyophilized protectants.

RESULTS

After the incorporation of 5 wt% PDLLA10K, the maximum loading of mPEG2k-b-PDLLA2.4k micelles for TM-2 was increased from 26 wt% to 32 wt%, and the in vivo half-life was increased by 2.25-fold. Various stability of micelles was improved. Also, the micelles with hydroxypropyl-β-cyclodextrin (HP-β-CD) as lyophilization protectants had minimal variation in particle size.

CONCLUSIONS

PDLLA10k incorporation can be employed as a strategy to increase the stability of mPEG2k-b-PDLLA2.4k micelles, which can be attributed to the viscosity building effect. HP-β-CD can be used as an effective lyophilization protectant since mPEG and HP-β-CD form the pseudopolyrotaxanesque inclusion complexes.

Effects of Graphene Oxidation on Interaction Energy and Interfacial Thermal Conductivity of Polymer Nanocomposite: A Molecular Dynamics Approach

Interfacial characteristics of polymer nanocomposites represent a crucial aspect to understand their global properties and to evaluate the interaction between nanofillers and matrix. In this work we used a molecular dynamics (MD) approach to characterize the interfacial region at the atomistic scale of graphene-based polymer nanocomposites. Three different polymer matrixes were considered, polylactic acid (PLA), polypropylene (PP) and epoxy resin (EPO), which were reinforced with three types of graphene fillers: pristine graphene (G), graphene oxide (GO) and reduced graphene oxide (rGO). In particular, the compatibility of the nanofillers in polymer matrixes were evaluated in terms of the interaction energy, while the interfacial thermal resistance (Kapitza resistance) between matrices and fillers was calculated with a nonequilibrium molecular dynamics (NEMD) method. Results showed that the oxidation degree plays an important role on the studied properties of the interfacial region. In particular, it was observed that the Kapitza resistance is decreased in the oxidized graphene (GO and rGO), while interaction energy depended on the polarity of the polymer matrix molecules and the contribution of the Coulombic component.

Investigation of nanoparticle-polymer interaction in bio-based nanosilica-filled PLA/NR nanocomposites: molecular dynamics simulation

Molecular dynamics (MD) simulation, by employing the COMPASS force field, was utilized to investigate structural and thermal characteristics as well as interfacial interactions between components of nanocomposite consisting of poly(lactic acid) (PLA)/natural rubber (NR)/nanosilica, abbreviated as PSxN, where 1 ≤ x ≤ 7 and it represents the parts of SiO₂ nanoparticles added to the PLA/NR (PN) blend. Analysis of the obtained results including density (ρ), fractional free volume (FFV), glass transition temperature (T_g), interaction energy (E_interaction), and radial distribution function (RDF) of these nanocomposites was performed. Comparing E_interaction of nanocomposites with that of the PN blend showed that the interactions between the chains of the two polymers are highly dependent on the added amounts of silica nanoparticles, so that by adding silica to the PN blend to obtain PS1N and PS3N nanocomposites, the amount of E_interaction was reduced to a smaller positive value, which indicates the tendency of the nanocomposite's components to interact with each other. By further addition of silica nanoparticles to have PS5N and PS7N nanocomposites and then by analysis of the RDF results, it was found that the nanoparticles were not well dispersed in these two nanocomposites and they were accumulated in the NR rubbery phase. Therefore, the percolation threshold for silica loading on the PN blend is at most 3 parts (x = 3). These results as well as the other obtained simulation results were compared with the available experimental data, and the agreement observed between them approved the simulation procedure and validated the obtained results.

Discovery and design of soft polymeric bio-inspired materials with multiscale simulations and artificial intelligence

Establishing the “Materials 4.0” paradigm requires intimate knowledge of the virtual space in materials design.

Structural characterization of poly-l-lactic acid (P(L)LA) and poly(glycolic acid)(PGA) oligomers

Structural characterization of poly-l-lactic acid (P(L)LA) and poly(glycolic acid) (PGA) oligomers containing three units was carried out with an atomistic approach. Oligomer structures were first optimized through quantum chemical calculations, using density functional theory (DFT); rotational barriers concerning dihedral angles along the chain were then investigated. Diffusion coefficients of l-lactic acid and glycolic acid in pure water were estimated through molecular dynamic (MD) simulations. Monomer structures were obtained with quantum chemical computation in implicit water using DFT method; atomic charges were fitted with Restrained Electrostatic Potentials (RESP) formalism, starting from electrostatic potentials calculated with quantum chemistry. MD simulations were carried out in explicit water, in order to take into account solvent presence.