Toward Predictive Molecular Dynamics Simulations of Asphaltenes in Toluene and Heptane

The 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). Systems with standard GAFF partial charges calculated by the AM1-BCC and HF/6-31G*(RESP) methods were simulated alongside systems without partial charges. The partial charges implemented differ in terms of the resulting electrical negativity of the asphaltene polyaromatic core, with the AM1-BCC method giving the greatest magnitude of the total core charge. Based on our analysis of the molecular relaxation and orientation, and on the aggregation behavior of asphaltenes in toluene and heptane, we proposed to use the partial charges obtained by the AM1-BCC method for the study of asphaltene aggregates. A good agreement with available experimental data was observed on the sizes of the aggregates, their fractal dimensions, and the solvent entrainment for the model asphaltenes in toluene and heptane. From the results obtained, we conclude that for a better predictive ability, simulation parameters must be carefully chosen, with particular attention paid to the partial charges owing to their influence on the electrical negativity of the asphaltene core and on the asphaltenes aggregation.

Transport Properties of Thermoplastic R-BAPB Polyimide: Molecular Dynamics Simulations and Experiment

The present work evaluates the transport properties of thermoplastic R-BAPB polyimide based on 1,3-bis(3,3',4,4'-dicarboxyphenoxy)benzene (dianhydride R) and 4,4'-bis(4-aminophenoxy)biphenyl (diamine BAPB). Both experimental studies and molecular dynamics simulations were applied to estimate the diffusion coefficients and solubilities of various gases, such as helium (He), oxygen (O2), nitrogen (N2), and methane (CH4). The validity of the results obtained was confirmed by studying the correlation of the experimental solubilities and diffusion coefficients of He, O2, and N2 in R-BAPB, with their critical temperatures and the effective sizes of the gas molecules, respectively. The solubilities obtained in the molecular dynamics simulations are in good quantitative agreement with the experimental data. A good qualitative relationship between the simulation results and the experimental data is also observed when comparing the diffusion coefficients of the gases. Analysis of the Robeson plots shows that R-BAPB has high selectivity for He, N2, and CO2 separation from CH4, which makes it a promising polymer for developing gas-separation membranes. From this point of view, the simulation models developed and validated in the present work may be put to effective use for further investigations into the transport properties of R-BAPB polyimide and nanocomposites based on it.

Atomistic Molecular Dynamics Simulations of the Initial Crystallization Stage in an SWCNT-Polyetherimide Nanocomposite

Crystallization 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. We selected two amorphous thermoplastic polyetherimides ODPA-P3 and aBPDA-P3 based on 3,3′,4,4′-oxydiphthalic dianhydride (ODPA), 2,3′,3,4′-biphenyltetracarboxylic dianhydride (aBPDA) and diamine 1,4-bis[4-(4-aminophenoxy)phenoxy]benzene (P3) and simulated the onset of crystallization in the presence of SWCNTs using atomistic molecular dynamics. For ODPA-P3, we found that the planar phthalimide and phenylene moieties show pronounced ordering near the CNT (carbon nanotube) surface, which can be regarded as the initial stage of crystallization. We will discuss two possible mechanisms for ODPA-P3 crystallization in the presence of SWCNTs: the spatial confinement caused by the CNTs and π–π interactions at the CNT-polymer matrix interface. Based on our simulation results, we propose that ODPA-P3 crystallization is most likely initiated by favorable π–π interactions between the carbon nanofiller surface and the planar ODPA-P3 phthalimide and phenylene moieties.

Molecular dynamics simulations of uniaxial deformation of thermoplastic polyimides

The results of atomistic molecular-dynamics simulations of mechanical properties of heterocyclic polymer subjected to uniaxial deformation are reported. A new amorphous thermoplastic polyimide R-BAPO with a repeat unit consisting of dianhydride 1,3-bis-(3',4,-dicarboxyphenoxy)diphenyl (dianhydride R) and diamine 4,4'-bis-(4''-aminophenoxy)diphenyloxide (diamine BAPO) was chosen for the simulations. Our primary goal was to establish the impact of various factors (sample preparation method, molecular mass, and cooling and deformation rates) on the elasticity modulus. In particular, we found that the elasticity modulus was only slightly affected by the degree of equilibration, the molecular mass and the size of the simulation box. This is most likely due to the fact that the main contribution to the elasticity modulus is from processes on scales smaller than the entanglement length. Essentially, our simulations reproduce the logarithmic dependence of the elasticity modulus on cooling and deformation rates, which is normally observed in experiments. With the use of the temperature dependence analysis of the elasticity modulus we determined the flow temperature of R-BAPO to be 580 K in line with the experimental data available. Furthermore, we found that the application of high external pressure to the polymer sample during uniaxial deformation can improve the mechanical properties of the polyimide. Overall, the results of our simulations clearly demonstrate that atomistic molecular-dynamics simulations represent a powerful and accurate tool for studying the mechanical properties of heterocyclic polymers and can therefore be useful for the virtual design of new materials, thereby supporting cost-effective synthesis and experimental research.

Influence of the carbon nanotube surface modification on the microstructure of thermoplastic binders

The structural properties of polymer nanocomposites based on thermoplastic polyimides filled with surface-modified carbon nanotubes (CNT) have been studied by means of fully-atomistic molecular-dynamics simulations.

Thermal properties of bulk polyimides: insights from computer modeling versus experiment

Due to the great importance for many industrial applications it is crucial from the point of view of theoretical description to reproduce thermal properties of thermoplastic polyimides as accurate as possible in order to establish "chemical structure-physical properties" relationships of new materials. In this paper we employ differential scanning calorimetry, dilatometry, and atomistic molecular dynamics (MD) simulations to explore whether the state-of-the-art computer modeling can serve as a precise tool for probing thermal properties of polyimides with highly polar groups. For this purpose the polyimide R-BAPS based on dianhydride 1,3-bis(3',4-dicarboxyphenoxy)benzene (dianhydride R) and diamine 4,4'-bis(4''-aminophenoxy)biphenyl sulphone) (diamine BAPS) was synthesized and extensively studied. Overall, our findings show that the widely used glass-transition temperature Tg evaluated from MD simulations should be employed with great caution for verification of the polyimide computational models against experimental data: in addition to the well-known impact of the cooling rate on the glass-transition temperature, correct definition of Tg requires cooling that starts from very high temperatures (no less than 800 K for considered polyimides) and accurate evaluation of the appropriate cooling rate, otherwise the errors in the measured values of Tg become undefined. In contrast to the glass-transition temperature, the volumetric coefficient of thermal expansion (CTE) does not depend on the cooling rate in the low-temperature domain (T < Tg) so that comparison of computational and experimental values of CTE provides a much safer way for proper validation of the theoretical model when electrostatic interactions are taken into account explicitly. Remarkably, this conclusion is most likely of generic nature: we show that it also holds for the commercial polyimide EXTEM, another polyimide with a similar chemical structure.

Influence of the carbon nanofiller surface curvature on the initiation of crystallization in thermoplastic polymers

The segments of crystallizable polyimide tend to lay parallel to the graphene nanofiller surface and this tendency is stronger than for carbon nanotubes.