RESUMEN

In this study, we have conducted a comparative analysis of the structural ordering of short oligoetherimide chains (dimers) near the bounding surface, depending on the structure of that surface. In order to clarify the possibility of oligoetherimide ordering along the symmetry axes of graphene, two types of bounding surfaces were considered: graphene, with a regular discrete position of interaction centers (carbon atoms), and a smooth, structureless impermeable wall. The chemical structures of the considered dimers consist of two repeating units of BPDA-P3, ODPA-P3, or aBPDA-P3 thermoplastic polyetherimides. Using all-atom molecular dynamics simulations, the process of structural ordering of the dimers near the surface of the graphene or wall was established. The ODPA-P3 and BPDA-P3 dimers form an ordered state near the graphene surface, while the aBPDA-P3 dimers do not demonstrate structural ordering. The simulation results confirmed that the ordering direction of the BPDA-P3 and ODPA-P3 dimers near the graphene surface is chosen randomly. Comparison of the oligoetherimide structure formed near the attracting wall without a symmetrical location of the interaction centers shows the similarity of the ordering of dimers near the graphene surface and the wall. As in the case of the graphene surface, the ordering of oligoetherimide molecules near the structureless wall demonstrates one direction of ordering. Therefore, we confirmed that the key factor for the onset of ordering is the presence of a confining surface, rather than the symmetrical arrangement of interaction centers in the substrate structure.

RESUMEN

Recent experiments and atomistic computer simulations have shown that asphaltene byproducts of oil refineries can serve as thermal conductivity enhancers for organic phase-change materials such as paraffin and therefore have the potential to improve the performance of paraffin-based heat storage devices. In this work, we explore how the size of the polycyclic aromatic cores of asphaltenes affects the properties of paraffin-asphaltene systems by means of atomistic molecular dynamics simulations. We show that increasing the size of the asphaltene core from 7-8 aromatic rings to ∼20 rings drastically changes the aggregation behavior of asphaltenes. Instead of relatively small, compact aggregates formed by small-core asphaltene molecules, enlarged cores promote the formation of extended single-column structures stabilized in paraffin by asphaltene's aliphatic periphery. Chemical modification of the asphaltenes by removing the periphery leads to the formation of bundles of columns. In contrast to small-core molecules, asphaltenes with enlarged cores do not suppress paraffin crystallization even at high filler concentrations. Remarkably, asphaltenes with enlarged aromatic cores are able to increase the thermal conductivity of liquid paraffin to a greater extent compared to their small-core counterparts. This effect becomes even more pronounced for modified asphaltenes without the aliphatic side groups. Overall, our computational findings suggest that asphaltenes with enlarged aromatic cores can significantly improve the performance of heat storage devices based on organic phase change materials.

RESUMEN

Asphaltenes represent a novel class of carbon nanofillers that are of potential interest for many applications, including polymer nanocomposites, solar cells, and domestic heat storage devices. In this work, we developed a realistic coarse-grained Martini model that was refined against the thermodynamic data extracted from atomistic simulations. This allowed us to explore the aggregation behavior of thousands of asphaltene molecules in liquid paraffin on a microsecond time scale. Our computational findings show that native asphaltenes with aliphatic side groups form small clusters that are uniformly distributed in paraffin. The chemical modification of asphaltenes via cutting off their aliphatic periphery changes their aggregation behavior: modified asphaltenes form extended stacks whose size increases with asphaltene concentration. At a certain large concentration (44 mol. %), the stacks of modified asphaltenes partly overlap, leading to the formation of large, disordered super-aggregates. Importantly, the size of such super-aggregates increases with the simulation box due to phase separation in the paraffin-asphaltene system. The mobility of native asphaltenes is systematically lower than that of their modified counterparts since the aliphatic side groups mix with paraffin chains, slowing down the diffusion of native asphaltenes. We also show that diffusion coefficients of asphaltenes are not very sensitive to the system size: enlarging the simulation box results in some increase in diffusion coefficients, with the effect being less pronounced at high asphaltene concentrations. Overall, our findings provide valuable insight into the aggregation behavior of asphaltenes on spatial and time scales that are normally beyond the scales accessible for atomistic simulations.

RESUMEN

In the present work, we address the problem of utilizing machine learning (ML) methods to predict the thermal properties of polymers by establishing "structure-property" relationships. Having focused on a particular class of heterocyclic polymers, namely polyimides (PIs), we developed a graph convolutional neural network (GCNN), being one of the most promising tools for working with big data, to predict the PI glass transition temperature T g as an example of the fundamental property of polymers. To train the GCNN, we propose an original methodology based on using a "transfer learning" approach with an enormous "synthetic" data set for pretraining and a small experimental data set for its fine-tuning. The "synthetic" data set contains more than 6 million combinatorically generated repeating units of PIs and theoretical values of their T g values calculated using the well-established Askadskii's quantitative structure-property relationship (QSPR) computational scheme. Additionally, an experimental data set for 214 PIs was also collected from the literature for training, fine-tuning, and validation of the GCNN. Both "synthetic" and experimental data sets are included into a PolyAskInG database (Polymer Askadskii's Intelligent Gateway). By using the PolyAskInG database, we developed GCNN which allows estimation of T g of PI with a mean absolute error (MAE) of about 20 K, which is 1.5 times lower than in the case of Askadskii QSPR analysis (33 K). To prove the efficiency and usability of the proposed GCNN architecture and training methodology for predicting polymer properties, we also employed "transfer learning" to develop alternative GCNN pretrained on proxy-characteristics taken from the popular quantum-chemical QM9 database for small compounds and fine-tuned on an experimental T g values data set from PolyAskInG database. The obtained results indicate that pretraining of GCNN on the "synthetic" polymer data set provides MAE which is almost twice as low as that in the case of using the QM9 data set in the pretraining stage (∼41 K). Furthermore, we address the questions associated with the influence of the differences in the size of the experimental and "synthetic" data sets (so-called "reality gap" problem), as well as their chemical composition on the training quality. Our results state the overall priority of using polymer data sets for developing deep neural networks, and GCNN in particular, for efficient prediction of polymer properties. Moreover, our work opens up a challenge for the theoretically supported generation of large "synthetic" data sets of polymer properties for the training of the complex ML models. The proposed methodology is rather versatile and may be generalized for predicting other properties of different polymers and copolymers synthesized through the polycondensation reaction.

RESUMEN

A 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. To this end, here we performed atomic-scale molecular dynamics simulations of organic phase-change materials (paraffins), in which the cooling rate was varied over four orders of magnitude. Our computational results clearly show that a certain threshold (1.2 × 1011 K/min) in the values of cooling rates exists. When cooling is slower than the threshold, the simulations qualitatively reproduce an experimentally observed abrupt change in the temperature dependence of the density, enthalpy, and thermal conductivity of paraffins upon crystallization. Beyond this threshold, when cooling is too fast, the paraffin's properties in simulations start to deviate considerably from experimental data: the faster the cooling, the larger part of the system is trapped in the supercooled liquid state. Thus, a proper choice of a cooling rate is of tremendous importance in computer simulations of organic phase-change materials, which are of great promise for use in domestic heat storage devices.

Asunto(s)

RESUMEN

Recently, a strong structural ordering of thermoplastic semi-crystalline polyimides near single-walled carbon nanotubes (SWCNTs) was found that can enhance their mechanical properties. In this study, a comparative analysis of the results of microsecond-scale all-atom computer simulations and experimental measurements of thermoplastic semi-crystalline polyimide R-BAPB synthesized on the basis of dianhydride R (1,3-bis-(3',4-dicarboxyphenoxy) benzene) and diamine BAPB (4,4'-bis-(4″-aminophenoxy) biphenyl) near the SWCNTs on the rheological properties of nanocomposites was performed. We observe the viscosity increase in the SWCNT-filled R-BAPB in the melt state both in computer simulations and experiments. For the first time, it is proven by computer simulation that this viscosity change is related to the structural ordering of the R-BAPB in the vicinity of SWCNT but not to the formation of interchain linkage. Additionally, strong anisotropy of the rheological properties of the R-BAPB near the SWCNT surface was detected due to the polyimide chain orientation. The increase in the viscosity of the polymer in the viscous-flow state and an increase in the values of the mechanical characteristics (Young's modulus and yield peak) of the SWCNT-R-BAPB nanocomposites in the glassy state are stronger in the directions along the ordering of polymer chains close to the carbon nanofiller surface. Thus, the new experimental data obtained on the R-BAPB-based nanocomposites filled with SWCNT, being extensively compared with simulation results, confirm the idea of the influence of macromolecular ordering near the carbon nanotube on the mechanical characteristics of the composite material.

RESUMEN

We studied the structure of brushes consisting of branched oligolactide (OLA) chains grafted onto the surface of cellulose nanoparticles (CNPs) in polylactide (PLA) and compared the outcomes to the case of grafting linear OLA chains using atomistic molecular dynamics simulations. The systems were considered in a melt state. The branched model OLA chains comprised one branching point and three branches, while the linear OLA chains examined had a molecular weight similar to the branched chains. It was shown that free branches of the branched OLA chains tend to fold back toward the CNPs due to dipole-dipole interactions within the grafted layer, in contrast to the well-established behavior of the grafted uncharged branched chains. This result, however, is in qualitative agreement with the conformational behavior known for linear OLA chains. At the same time, no significant difference in the effectiveness of covering the filler surface with grafted branched or linear OLA chains was found. In terms of the expelling ability of the grafted chains and the interaction between PLA and CNP or OLA, the linear chains were broadly similar (sparse grafting) or better (intermediate or dense grafting) compared to the branched ones. Thus, the grafted lactide chains with a linear architecture, rather than their branched counterpart, may be preferable for the covalent modification of cellulose nanoparticles.

RESUMEN

[This corrects the article DOI: 10.1039/C9RA07325F.].

RESUMEN

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.

RESUMEN

The structure of a grafted layer of lactide chains in the "dry brush" regime immersed in a melt of chemically similar polymer was examined while varying graft lengths. To this end, microsecond atomistic molecular dynamics simulations were performed. Almost no influence of graft length on the fraction of the grafted chains backfolded to the grafting surface was found. However, a structural ordering was unexpectedly observed in the system when the length of the grafted lactide chains was close to approximately 10 Kuhn segments. This ordering of the grafts is characterized by the formation of helical fragments whose structure is in good agreement with the experimental data for the α crystal of the lactide chains. Both the backfolding and the structural ordering may be viewed as the initial stage of the crystallization of the layer of grafted lactide chains. In contrast to the known behavior for conventional polymer brushes in the "dry brush" regime, the structure of the grafted lactide chains can be either amorphous or ordered, depending on the graft length N and the grafting density σ when their product Nσ is fixed.

RESUMEN

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.

RESUMEN

Paraffin-based composites represent a promising class of materials with numerous practical applications such as e.g. heat storage. Computer modeling of these complex multicomponent systems requires a proper theoretical description of both the n-alkane matrix and the non-alkane filler molecules. The latter can be modeled with the use of a state-of-the-art general-purpose force field such as GAFF, CHARMM, OPLS-AA and GROMOS, while the paraffin matrix is traditionally described in the frame of relatively old, alkane-specific force fields (TraPPE, NERD, and PYS). In this paper we link these two types of models and evaluate the performance of several general-purpose force fields in computer modeling of paraffin by their systematic comparison with earlier alkane-specific models as well as with experimental data. To this end, we have performed molecular dynamics simulations of n-eicosane bulk samples with the use of 10 different force fields: TraPPE, NERD, PYS, OPLS-UA, GROMOS, GAFF, GAFF2, OPLS-AA, L-OPLS-AA, and CHARMM36. For each force field we calculated several thermal, structural and dynamic characteristics of n-eicosane over a wide temperature range. Overall, our findings show that the general-purpose force fields such as CHARMM36, L-OPLS-AA and GAFF/GAFF2 are able to provide a realistic description of n-eicosane samples. While alkane-specific models outperform most general-purpose force fields as far as the temperature dependence of mass density, the coefficient of volumetric thermal expansion in the liquid state, and the crystallization temperature are concerned, L-OPLS-AA, CHARMM36 and GAFF2 force fields provide a better match with experiment for the shear viscosity and the diffusion coefficient in melt. Furthermore, we show that most general-purpose force fields are able to reproduce qualitatively the experimental triclinic crystal structure of n-eicosane at low temperatures.

RESUMEN

Using 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. The structural properties of the nanocomposites at temperature close to the bulk polymer melting point were studied. The mechanical properties were determined for the nanocomposites cooled down to the glassy state. It was found that the SWCNT nanofiller initiates' structural ordering not only in the crystallizable BPDA-P3 but also in the amorphous ODPA-P3 samples were in agreement with previously obtained experimental results. Two stages of the structural ordering were detected in the presence of SWCNTs, namely the orientation of the planar moieties followed by the elongation of whole polymer chains. The first type of local ordering was observed on the microsecond time scale and did not lead to the change of the mechanical properties of a polymer binder in considered nanocomposites. At the end of the second stage, both BPDA-P3 and ODPA-P3 PI chains extended completely along the SWCNT surface, which in turn led to enhanced mechanical characteristics in their glassy state.

RESUMEN

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.

RESUMEN

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.

RESUMEN

The complexes formed by two dendrimers with charged terminal groups and oppositely charged long linear polyelectrolyte (LPE) have been studied using Brownian dynamics simulations. The structural properties of the complexes and their dependence on the LPE chain length were investigated. It was observed that dendrimers in the considered complexes are sufficiently overcharged; i.e., the number of adsorbed LPE monomers is larger than required for the neutralization. The degree of overcharging increases with the increase of the LPE length and is accompanied by the linker appearance until saturation in overcharging is reached. Nonmonotonic dependence of the linker size on the LPE length was observed. To describe the structural properties of the complexes formed by two macroions and a polyelectrolyte chain, the correlation theory has been developed.

Asunto(s)

RESUMEN

Molecular dynamics simulations are used to probe the structural organization of nonstoichiometric interpolyelectrolyte complexes (IPECs) formed by oppositely charged starlike and linear polyelectrolytes (PEs) in dilute aqueous solution. We demonstrate that undercompensated star-IPEC consists of a denser coacervate core and a charged starlike corona. Two distinctive populations of star branches completely embedded in a coacervate core and stretched in a lyophilizing corona are found. The scaling arguments support the stability of IPEC with partitioned star branches.

Asunto(s)