Atomistic molecular dynamics simulation of the temperature and pressure dependences of local and terminal relaxations in cis-1,4-polybutadiene

The effect of functional groups on the glass transition temperature of atmospheric organic compounds: a molecular dynamics study

Organic compounds constitute a substantial part of atmospheric particulate matter not only in terms of mass concentration but also in terms of distinct functional groups. The glass transition temperature provides an indirect way to investigate the phase state of the organic compounds, playing a crucial role in understanding their behavior and influence on aerosol processes. Molecular dynamics (MD) simulations were implemented here to predict the glass transition temperature (Tg) of atmospherically relevant organic compounds as well as the influence of their functional groups and length of their carbon chain. The cooling step used in the simulations was chosen to be neither too low (to supress crystallization) nor too high (to avoid Tg overprediction). According to the MD simulations, the predicted Tg is sensitive to the functional groups as follows: carboxylic acid (-COOH) > hydroxyl (-OH) and (-COOH) > carbonyls (-CO). Increasing the number of carbon atoms leads to higher Tg for the linearly structured compounds. Linear compounds with lower molecular weight were found to exhibit a lower Tg. No clear correlation between O : C and Tg was observed. The architecture of the carbon chain (linear, or branched, or ring) was also found to impact the glass transition temperature. Compounds containing a non-aromatic carbon ring are characterized by a higher Tg compared to linear and branched ones with the same number of carbon atoms.

Amorphous cis-1,4-polybutadiene P-V-T properties from atomistic simulations

CONTEXT

As a result of the diversity of microstructures encountered in cis-1,4-polybutadiene and the variety of measurement methods used, experimental values of variation of glass transition temperature (Tg) with pressure are relatively dispersed. However, atomistic simulations enable access to valuable information for very well-controlled chemistry and structures with a well-defined and systematic acquisition protocol. By varying the temperature and pressure, the specific volume of the melt was computed, yielding results that deviated by only 2% from experimental data. A linear relationship between Tg and pressure was observed, with Tg predicted to be 162 K at zero pressure and a rate of change of Tg with respect to pressure (dTg/dP) of 0.24 K/MPa.

METHOD

The atomistic dilatometry experiments were conducted on a model of amorphous cis-1,4 polybutadiene with an approximate molecular weight of 5400 g/mol using the LAMMPS code and the all-atom forcefield pcff + . The dilatometry process involved cooling and heating at a rate of 9 × 10¹² K/min. The specific volume was calculated by averaging over seven independent configurations for each temperature. The Tait equation was employed to fit the specific volume evolution within the temperature range of 10 to 700 K under different pressures of 0, 60, and 100 MPa.

Simulational Tests of the Rouse Model

An extensive review of literature simulations of quiescent polymer melts is given, considering results that test aspects of the Rouse model in the melt. We focus on Rouse model predictions for the mean-square amplitudes ⟨(Xp(0))2⟩ and time correlation functions ⟨Xp(0)Xp(t)⟩ of the Rouse mode Xp(t). The simulations conclusively demonstrate that the Rouse model is invalid in polymer melts. In particular, and contrary to the Rouse model, (i) mean-square Rouse mode amplitudes ⟨(Xp(0))2⟩ do not scale as sin-2(pπ/2N), N being the number of beads in the polymer. For small p (say, p≤3) ⟨(Xp(0))2⟩ scales with p as p-2; for larger p, it scales as p-3. (ii) Rouse mode time correlation functions ⟨Xp(t)Xp(0)⟩ do not decay with time as exponentials; they instead decay as stretched exponentials exp(-αtβ). β depends on p, typically with a minimum near N/2 or N/4. (iii) Polymer bead displacements are not described by independent Gaussian random processes. (iv) For p≠q, ⟨Xp(t)Xq(0)⟩ is sometimes non-zero. (v) The response of a polymer coil to a shear flow is a rotation, not the affine deformation predicted by Rouse. We also briefly consider the Kirkwood-Riseman polymer model.

Arriving at the most plausible interpretation of the dielectric spectra of glycerol with help from quasielastic γ-ray scattering time-domain interferometry

Glycerol is one of the glass-forming liquids selected by Robert H. Cole in 1950 to start his study of molecular dynamics by dielectric spectroscopy. Seventy-one years have gone by and remarkably no consensus has been reached on the nature and identity of the relaxation processes observed in the dielectric spectra. The macroscopic dielectric relaxation data allow different interpretations to yield contrasting results, and it is not possible to determine which one is most plausible. Coming to the rescue is the application of the nuclear γ-resonance time-domain interferometry (TDI) to glycerol by Saito et al. [Phys. Rev. E 105, L012605 (2022)10.1103/PhysRevE.105.L012605]. Their microscopic TDI data potentially can decide which interpretation of the dielectric spectra of glycerol is most plausible. The attempt was made by Saito et al., but there is a problem in their analysis of the dielectric data of glycerol and hence their conclusion is untenable. In this paper, we critically compare four major interpretations with the TDI data in an effort to identify the most plausible interpretation of the relaxation processes constituting the dielectric spectra of glycerol.

Isochronal superpositioning of the caged dynamics, the α, and the Johari-Goldstein β relaxations in metallic glasses

The superposition of the frequency dispersions of the structural α relaxation determined at different combinations of temperature T and pressure P while maintaining its relaxation time τ_α(T, P) constant (i.e., isochronal superpositioning) has been well established in molecular and polymeric glass-formers. Not known is whether the frequency dispersion or time dependence of the faster processes including the caged molecule dynamics and the Johari-Goldstein (JG) β relaxation possesses the same property. Experimental investigation of this issue is hindered by the lack of an instrument that can cover all three processes. Herein, we report the results from the study of the problem utilizing molecular dynamics simulations of two different glass-forming metallic alloys. The mean square displacement 〈Δr²t〉, the non-Gaussian parameter α₂t, and the self-intermediate scattering function F_sq,t at various combinations of T and P were obtained over broad time range covering the three processes. Isochronal superpositioning of 〈Δr²t〉, α₂t, and F_sq,t was observed over the entire time range, verifying that the property holds not only for the α relaxation but also for the caged dynamics and the JG β relaxation. Moreover, we successfully performed density ρ scaling of the time τ_α2,maxT,P at the peak of α₂t and the diffusion coefficient D(T, P) to show both are functions of ρ^γ/T with the same γ. It follows that the JG β relaxation time τ_β(T, P) is also a function of ρ^γ/T since τ_α2,maxT,P corresponds to τ_β(T, P).

Why is the change of the Johari–Goldstein β-relaxation time by densification in ultrastable glass minor?

Coupling-Model-based theoretical explanation of the minor change of JG β-relaxation achieved by ultrastability in contrast to the dramatic change in α-relaxation.

Scaling of the dynamics of flexible Lennard-Jones chains: Effects of harmonic bonds

The previous paper [A. A. Veldhorst et al., J. Chem. Phys. 141, 054904 (2014)] demonstrated that the isomorph theory explains the scaling properties of a liquid of flexible chains consisting of ten Lennard-Jones particles connected by rigid bonds. We here investigate the same model with harmonic bonds. The introduction of harmonic bonds almost completely destroys the correlations in the equilibrium fluctuations of the potential energy and the virial. According to the isomorph theory, if these correlations are strong a system has isomorphs, curves in the phase diagram along which structure, dynamics, and the excess entropy are invariant. The Lennard-Jones chain liquid with harmonic bonds does have curves in the phase diagram along which the structure and dynamics are invariant. The excess entropy is not invariant on these curves, which we refer to as "pseudoisomorphs." In particular, this means that Rosenfeld's excess-entropy scaling (the dynamics being a function of excess entropy only) does not apply for the Lennard-Jones chain with harmonic bonds.

Scaling of the dynamics of flexible Lennard-Jones chains

The isomorph theory provides an explanation for the so-called power law density scaling which has been observed in many molecular and polymeric glass formers, both experimentally and in simulations. Power law density scaling (relaxation times and transport coefficients being functions of ρ(γ(S)), where ρ is density, T is temperature, and γ(S) is a material specific scaling exponent) is an approximation to a more general scaling predicted by the isomorph theory. Furthermore, the isomorph theory provides an explanation for Rosenfeld scaling (relaxation times and transport coefficients being functions of excess entropy) which has been observed in simulations of both molecular and polymeric systems. Doing molecular dynamics simulations of flexible Lennard-Jones chains (LJC) with rigid bonds, we here provide the first detailed test of the isomorph theory applied to flexible chain molecules. We confirm the existence of isomorphs, which are curves in the phase diagram along which the dynamics is invariant in the appropriate reduced units. This holds not only for the relaxation times but also for the full time dependence of the dynamics, including chain specific dynamics such as the end-to-end vector autocorrelation function and the relaxation of the Rouse modes. As predicted by the isomorph theory, jumps between different state points on the same isomorph happen instantaneously without any slow relaxation. Since the LJC is a simple coarse-grained model for alkanes and polymers, our results provide a possible explanation for why power-law density scaling is observed experimentally in alkanes and many polymeric systems. The theory provides an independent method of determining the scaling exponent, which is usually treated as an empirical scaling parameter.

Dielectric α-relaxation of 1,4-polybutadiene confined between graphite walls : Molecular dynamics investigations through numerical simulations of polymer molecules relaxation

We present results of Molecular Dynamics (MD) simulations of a chemically realistic model of 1,4-polybutadiene confined by crystalline graphite walls. The simulations cover a large range of temperatures from T ≈ 2T g to T ≈ 1.15T g, where relevant time scales are accessible using such computational methods. We investigate the dielectric relaxation close to the walls in comparison to the one in the center of the film, and study the latter as a function of the film thickness from the walls. The segmental dynamics in the film is slowed down close to the walls, in comparison to the bulk. In addition to the α-process, the relaxation exhibits an additional long time decay, the so-called wall desorption process. We focus here on the α-process and find no significant shift of the dielectric T g as a function of layer thickness, in agreement with recent dielectric experiments. These findings can be correlated with the importance of the dihedral dynamics for all relaxation processes in polymers, which is unaltered except for the first nanometer next to the walls.

Effect of chain length on fragility and thermodynamic scaling of the local segmental dynamics in poly(methylmethacrylate)

Local segmental relaxation properties of poly(methylmethacrylate) (PMMA) of varying molecular weight are measured by dielectric spectroscopy and analyzed in combination with the equation of state obtained from PVT measurements. Significant variations of glass transition temperature and fragility with molecular weight are observed. In accord with the general properties of glass-forming materials, single molecular weight dependent scaling exponent gamma is sufficient to define the mean segmental relaxation time taualpha and its distribution. This exponent can be connected to the Gruneisen parameter and related thermodynamic quantities, thus demonstrating the interrelationship between dynamics and thermodynamics in PMMA. Changes in the relaxation properties ("dynamic crossover") are observed as a function of both temperature and pressure, with taualpha serving as the control parameter for the crossover. At longer taualpha another change in the dynamics is apparent, associated with a decoupling of the local segmental process from ionic conductivity.

Thermodynamic scaling of the viscosity of van der Waals, H-bonded, and ionic liquids

Viscosities eta and their temperature T and volume V dependences are reported for seven molecular liquids and polymers. In combination with literature viscosity data for five other liquids, we show that the superpositioning of relaxation times for various glass-forming materials when expressed as a function of TV(gamma), where the exponent gamma is a material constant, can be extended to the viscosity. The latter is usually measured to higher temperatures than the corresponding relaxation times, demonstrating the validity of the thermodynamic scaling throughout the supercooled and higher T regimes. The value of gamma for a given liquid principally reflects the magnitude of the intermolecular forces (e.g., steepness of the repulsive potential); thus, we find decreasing gamma in going from van der Waals fluids to ionic liquids. For some strongly H-bonded materials, such as low molecular weight polypropylene glycol and water, the superpositioning fails, due to the nontrivial change of chemical structure (degree of H bonding) with thermodynamic conditions.

Glass transition in 1,4-polybutadiene: Mode-coupling theory analysis of molecular dynamics simulations using a chemically realistic model

We present molecular dynamics simulations of the glass transition in a chemically realistic model of 1,4-polybutadiene (PBD). Around 40 K above the calorimetric glass transition of this polymer the simulations reveal a well-developed two-stage relaxation of all correlation functions. We have analyzed the time-scale separation between vibrational degrees of freedom (subpicosecond dynamics) and the alpha relaxation behavior (nanosecond to microsecond dynamics) using the predictions of mode-coupling theory (MCT). Our value for the mode-coupling critical temperature Tc agrees perfectly with prior experimental estimates for PBD. The predictions of MCT for the scaling behavior of the so-called beta relaxation, the plateau regime separating vibrational dynamics and the alpha relaxation, are well fulfilled. Furthermore, we are able to derive a consistent set of MCT exponents, completely characterizing the scaling behavior of relaxation processes in the vicinity of Tc. For the temperature dependence of the alpha relaxation we find deviations from MCT predictions which we trace to the influence of torsional barriers on the atomic motions.