influence of specific intermolecular interactions on the thermal and dielectric properties of bulk polymers atomistic molecular dynamics simulations of nylon 6

Specific intermolecular interactions, in particular H-bonding, have a strong influence on the structural, thermal and relaxation characteristics of polymers. We report here the results of molecular dynamics simulations of Nylon 6 which provides an excellent example for the investigation of such an influence. To demonstrate the effect of proper accounting for H-bonding on bulk polymer properties, the AMBER99sb force field is used with two different parametrization approaches leading to two different sets of partial atomic charges. The simulations allowed the study of the thermal and dielectric properties in a wide range of temperatures and cooling rates. The feasibility of the use of the three methods for the estimation of the glass transition temperature not only from the temperature dependence of structural characteristics such as density, but also by using the electrostatic energy and dielectric constant is demonstrated. The values of glass transition temperatures obtained at different cooling rates are practically the same for the three methods. By proper accounting for partial charges in the simulations, a reasonable agreement between the results of our simulations and experimental data for the density, thermal expansion coefficient, static dielectric constant and activation energy of γ and β relaxations is obtained demonstrating the validity of the modeling approach reported

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

Welding dynamics in an atomistic model of an amorphous polymer blend with polymer-polymer interface

We consider an atomistic model of thermal welding at the polymer-polymer interface of a polyetherimide/polycarbonate blend, motivated by applications to 3D manufacturing in space. We follow diffusion of semiflexible chains at the interface and analyze strengthening of the samples as a function of the welding time tw by simulating the strain-stress and shear viscosity curves. The time scales for initial wetting, and for fast and slow diffusion, are revealed. It is shown that each component of the polymer blend has its own characteristic time of slow diffusion at the interface. Analysis of strainstress demonstrates saturation of the Young’s modulus at tw = 240 ns, while the tensile strength continues to increase. The shear viscosity is found to have a very weak dependence on the welding time for tw > 60 ns. It is shown that both strain-stress and shear viscosity curves agree with experimental data

The Transport Properties of Semi-Crystalline Polyetherimide BPDA-P3 in Amorphous and Ordered States: Computer Simulations

The effect of polymer chain ordering on the transport properties of the polymer membrane was examined for the semi-crystalline heterocyclic polyetherimide (PEI) BPDA-P3 based on 3,3′,4,4′-biphenyltetracarboxylic dianhydride (BPDA) and diamine 1,4-bis [4-(4-aminophenoxy)phenoxy]benzene (P3). All-atom Molecular Dynamics (MD) simulations were used to investigate the gas diffusion process carried through the pores of a free volume several nanometers in size. The long-term (~30 μs) MD simulations of BPDA-P3 were performed at T = 600 K, close to the experimental value of the melting temperature (Tm ≈ 577 K). It was found during the simulations that the transition of the PEI from an amorphous state to an ordered one occurred. We determined a decrease in solubility for both the gases examined (CO2 and CH4), caused by the redistribution of free volume elements occurring during the structural ordering of the polymer chains in glassy state (Tg ≈ 481 K). By analyzing the diffusion coefficients in the ordered state, the presence of gas diffusion anisotropy was found. However, the averaged values of the diffusion coefficients did not differ from each other in the amorphous and ordered states. Thus, permeability in the observed system is primarily determined by gas solubility, rather than by gas diffusion

Structural ordering in SWCNT-polyimide nanocomposites and its influence on their mechanical properties

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

ALGORITHM FOR OPTIMIZING THE COMBUSTION PROCESS OF ASSOCIATED PETROLEUM GAS IN THERMAL POWER PLANTS, TAKING INTO ACCOUNT THE VARIABILITY OF ITS COMPOSITION

Problems of natural resources rational using, energy efficiency and energy saving are topical at present. As part of this problem is actual task of rational and efficient use of associated petroleum gas. The composition of associated petroleum gas due to the nature of its production variable over time. These changes lead to the violation of the combustion optimal mode in thermal power plants. This article presents the algorithm of combustion optimization of hydrocarbon fuels variable composition, in particular of associated petroleum gas. There were developed the algorithms: the determination of initial ratio "fuel-air" in thermal power plants and to bring this ratio to the optimum; optimization of the combustion process during change of the specific heat of combustion. The algorithms are based on registration of change coolant temperature at the outlet of thermal power plant

Simulating local mobility and mechanical properties of thermostable polyimides with different dianhydride fragments

\u3cp\u3eThe dynamic and mechanical properties of the three thermoplastic polyimides—crystallizable polyimide BPDA-P3 and amorphous polyimides ODPA-P3 and aBPDA-P3—have been simulated using the atomistic molecular dynamics technique. The three simulated polyimides differ in the chemical structure of their corresponding dianhydride fragments. Analyzing the local orientational mobility of different phenylene rings, it has been established that the increase of the glass-transition temperature (T \u3csub\u3eg\u3c/sub\u3e) in the ordered set (Formula presented.) is caused by the slowing down of the phthalimide rings relaxation in the corresponding dianhydride fragments. It has been observed that rather poor mechanical characteristics upon aBPDA-P3 stretching in the strain-hardening regime are also due to the low orientational mobility of the phthalimide rings. The correlation between the dynamic fragility and the polyimides strain-hardening moduli has been observed; the increase of the dynamic fragility leads to the increase of the strain-hardening modulus.\u3c/p\u3

Molecular dynamics simulations of uniaxial deformation of thermoplastic polyimides

\u3cp\u3eThe 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.\u3c/p\u3

Structural ordering in SWCNT-polyimide nanocomposites and its influence on their mechanical properties

\u3cp\u3eUsing 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.\u3c/p\u3

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