Fardin Khabaz

Effect of Carbon Nanotube Functionalization on Mechanical and Thermal Properties of Cross-Linked Epoxy–Carbon Nanotube Nanocomposites: Role of Strengthening the Interfacial Interactions

We have used amido-amine functionalized carbon nanotubes (CNTs) that form covalent bonds with cross-linked epoxy matrices to elucidate the role of the matrix-filler interphase in the enhancement of mechanical and thermal properties in these nanocomposites. For the base case of nanocomposites of cross-linked epoxy and pristine single-walled CNTs, our previous work (Khare, K. S.; Khare, R. J. Phys. Chem. B 2013, 117, 7444-7454) has shown that weak matrix-filler interactions cause the interphase region in the nanocomposite to be more compressible. Furthermore, because of the weak matrix-filler interactions, the nanocomposite containing dispersed pristine CNTs has a glass transition temperature (Tg) that is ∼66 K lower than the neat polymer. In this work, we demonstrate that in spite of the presence of stiff CNTs in the nanocomposite, the Youngs modulus of the nanocomposite containing dispersed pristine CNTs is virtually unchanged compared to the neat cross-linked epoxy. This observation suggests that the compressibility of the matrix-filler interphase interferes with the ability of the CNTs to reinforce the matrix. Furthermore, when the compressibility of the interphase is reduced by the use of amido-amine functionalized CNTs, the mechanical reinforcement due to the filler is more effective, resulting in a ∼50% increase in the Youngs modulus compared to the neat cross-linked epoxy. Correspondingly, the functionalization of the CNTs also led to a recovery in the Tg making it effectively the same as the neat polymer and also resulted in a ∼12% increase in the thermal conductivity of the nanocomposite containing functionalized CNTs compared to that containing pristine CNTs. These results demonstrate that the functionalization of the CNTs facilitates the transfer of both mechanical load and thermal energy across the matrix-filler interface.

Molecular Topology and Local Dynamics Govern the Viscosity of Imidazolium-Based Ionic Liquids.

A series of branched ionic liquids (ILs) based on the 1-(iso-alkyl)-3-methylimidazolium cation from 1-(1-methylethyl)-3-methylimidazolium bistriflimide to 1-(5-methylhexyl)-3-methylimidazolium bistriflimide and linear ILs based on the 1-(n-alkyl)-3-methylimidazolium cation from 1-propyl-3-methylimidazolium bistriflimide to 1-heptyl-3-methylimidazolum bistriflimide were recently synthesized and their physicochemical properties characterized. For the ILs with the same number of carbons in the alkyl chain, the branched IL was found to have the same density but higher viscosity than the linear one. In addition, the branched IL 1-(2-methylpropyl)-3-methylimidazolium bistriflimide ([2mC3C1Im][NTf2]) was found to have an abnormally high viscosity. Motivated by these experimental observations, the same ILs were studied using molecular dynamics (MD) simulations in the current work. The viscosities of each IL were calculated using the equilibrium MD method at 400 K and the nonequilibrium MD method at 298 K. The results agree with the experimental trend. The ion pair (IP) lifetime, spatial distribution function, and associated potential of mean force, cation size and shape, and interaction energy components were calculated from MD simulations. A quantitative correlation between the liquid structure and the viscosity was observed. Analysis shows that the higher viscosities in the branched ILs are due to the relatively more stable packing between the cations and anions indicated by the lower minima in the potential of mean force (PMF) surface. The abnormal viscosity of [2mC3C1Im][NTf2] was found to be the result of the specific side chain length and molecular structure.

Glass Transition and Molecular Mobility in Styrene-Butadiene Rubber Modified Asphalt.

Asphalt, a soft matter consisting of more than a thousand chemical species, is of vital importance for the transportation infrastructure, yet it poses significant challenges for microscopic theory and modeling approaches due to its multicomponent nature. Polymeric additives can potentially enhance the thermo-mechanical properties of asphalt, thus helping reduce the road repair costs; rational design of such systems requires knowledge of the molecular structure and dynamics of these systems. We have used molecular dynamics (MD) simulations to investigate the volumetric, structural, and dynamic properties of the neat asphalt as well as styrene-butadiene rubber (SBR) modified asphalt systems. The volume-temperature behavior of the asphalt systems exhibited a glass transition phenomenon, akin to that observed in experiments. The glass transition temperature, room temperature density, and coefficient of volume thermal expansion of the neat asphalt systems so evaluated were in agreement with experimental data when the effect of the high cooling rate used in simulations was accounted for. While the volumetric properties of SBR modified asphalt were found to be insensitive to the presence of the SBR additive, the addition of SBR led to an increase in the aggregation of asphaltene molecules. Furthermore, addition of SBR caused a reduction in the mobility of the constituent molecules of asphalt, with the reduction being more significant for the larger constituent molecules. Similar to other glass forming liquids, the reciprocal of the diffusion coefficient of the selected molecules was observed to follow the Vogel-Fulcher-Tammann (VFT) behavior as a function of temperature. These results suggest the potential for using polymeric additives for enhancing the dynamic mechanical properties of asphalt without affecting its volumetric properties.

Structure and Hydrogen Bonding of Water in Polyacrylate Gels: Effects of Polymer Hydrophilicity and Water Concentration.

The ability to tune the hydrophilicity of polyacrylate copolymers by altering their composition makes these materials attractive candidates for membranes used to separate alcohol-water mixtures. The separation behavior of these polyacrylate membranes is governed by a complex interplay of factors such as water and alcohol concentrations, water structure in the membrane, polymer hydrophilicity, and temperature. We use molecular dynamics simulations to investigate the effect of polymer hydrophilicity and water concentration on the structure and dynamics of water molecules in the polymer matrix. Samples of poly(n-butyl acrylate) (PBA), poly(2-hydroxyethyl acrylate) (PHEA), and a 50/50 copolymer of BA and HEA were synthesized in laboratory, and their properties were measured. Model structures of these systems were validated by comparing the simulated values of their volumetric properties with the experimental values. Molecular simulations of polyacrylate gels swollen in water and ethanol mixtures showed that water exhibits very different affinities toward the different (carbonyl, alkoxy, and hydroxyl) functional groups of the polymers. Water molecules are well dispersed in the system at low concentrations and predominantly form hydrogen bonds with the polymer. However, water forms large clusters at high concentrations along with the predominant formation of water-water hydrogen bonds and the acceleration of hydrogen bond dynamics.

Effect of chain architecture on the size, shape, and intrinsic viscosity of chains in polymer solutions: A molecular simulation study

Effect of chain architecture on the chain size, shape, and intrinsic viscosity was investigated by performing molecular dynamics simulations of polymer solutions in a good solvent. Four types of chains—linear, comb shaped, H-shaped, and star—were studied for this purpose using a model in which the solvent particles were considered explicitly. Results indicated that the chain length (N) dependence of the mean squared radius of gyration of the chains followed a power-law behavior ⟨R(g)(2)⟩(1/2)∼N(υ) with scaling exponents of υ = 0.605, 0.642, 0.602, and 0.608, for the linear, comb shaped, H-shaped, and star shaped chains, respectively. The simulation results for the geometrical shrinking factor were higher than the prior theoretical predictions for comb shaped chains. Analysis of chain shape demonstrated that the star chains were significantly smaller and more spherical than the others, while the comb and H-shaped polymer chains showed a more cylindrical shape. It is shown that the intrinsic viscosity of the chains can be calculated by plotting the specific viscosity determined from simulations against the solution concentration. The intrinsic viscosity exhibited linear behavior with the reciprocal of the overlap concentration for all chain architectures studied. The molecular weight dependence of the intrinsic viscosity followed the Mark-Houwink relation, [η] = KM(a), for all chain architectures. When comparing the calculated values of exponent a with the literature experimental values, agreement was found only for the H and star chains, and a disagreement for the linear and comb chains. The viscosity shrinking factor of the branched chains was compared with the available experimental data and the theoretical predictions and a general agreement was found.

Molecular simulations of asphalt rheology: Application of time–temperature superposition principle

We show that the several decades gap between the time scales of experiments and atomistic simulations can be significantly reduced by applying time-temperature superposition (TTS) to the simulation results for various rheological properties of asphalt. Molecular dynamics simulations with atomistically detailed models were performed to characterize the temperature dependence of shear viscosity, dynamic modulus, and tensile creep compliance of AAM–1 asphalt. The TTS principle was successfully applied to the data to construct master curves of these rheological properties, and a comparison of simulation data with experiments showed good agreement. Application of TTS resulted in a six decade increase in the frequency range available for moduli determination and a two decade increase in both the shear rate range available for viscosity determination and the time scale available for characterization of creep compliance. The horizontal shift factors (aT) used in the construction of the master curves for the different rheological properties were quantitatively consistent with each other. The corresponding temperature dependence of the scaled relaxation times was compared with the expectations from three literature models, namely, the Vogel–Fulcher–Tammann, double exponential, and parabolic–Arrhenius models. A distinct change from super-Arrhenius to Arrhenius behavior was observed for the scaled relaxation times below the glass transition temperature (Tg). Only the parabolic-Arrhenius model was able to capture the temperature dependence of the data over the entire temperature range that covers both the glassy and the rubbery states.We show that the several decades gap between the time scales of experiments and atomistic simulations can be significantly reduced by applying time-temperature superposition (TTS) to the simulation results for various rheological properties of asphalt. Molecular dynamics simulations with atomistically detailed models were performed to characterize the temperature dependence of shear viscosity, dynamic modulus, and tensile creep compliance of AAM–1 asphalt. The TTS principle was successfully applied to the data to construct master curves of these rheological properties, and a comparison of simulation data with experiments showed good agreement. Application of TTS resulted in a six decade increase in the frequency range available for moduli determination and a two decade increase in both the shear rate range available for viscosity determination and the time scale available for characterization of creep compliance. The horizontal shift factors (aT) used in the construction of the master curves for the diffe...

Temperature Dependence of Volumetric and Dynamic Properties of Imidazolium-Based Ionic Liquids

Atomistically detailed molecular dynamics simulations were used to investigate the temperature dependence of the specific volume, dynamic properties, and viscosity of linear alkyl chain ([CnC1Im][NTf2], n = 3-7) and branched alkyl chain ([(n - 2)mCn-1C1Im][NTf2]) ionic liquids (ILs). The trend of the glass transition temperature (Tg) values obtained in the simulations as a function of the alkyl chain length of cations was similar to the trend seen in experiments. In addition, the system relaxation behavior as determined from the temperature dependence of the diffusion coefficient, rotational relaxation time, and viscosity close to Tg was observed to follow the Vogel-Fulcher-Tammann expression. Furthermore, the reciprocal of the diffusion coefficient of the anion and cation in both linear and branched IL systems showed a linear correlation with viscosity, thus confirming the validity of the Stokes-Einstein relationship for these systems. Similarly, the average rotational relaxation time of the ions was also found to correlate linearly with the viscosity of the ILs over a wide range of temperatures, thereby validating the Debye-Stokes-Einstein relationship for the ILs. These simulation findings suggest that the temperature dependence of the relaxation time of ILs is very similar to that of other glass-forming liquids.

Swelling of Random Copolymer Networks in Pure and Mixed Solvents: Multi-Component Flory–Rehner Theory

A generalized extension of Flory-Rehner (FR) theory is derived to describe equilibrium swelling of polymer networks, including copolymers with two or more chemically distinct repeat units, in either pure or mixed solvents. The model is derived by equating the chemical potential of each solvent in the liquid and gel phases at equilibrium, while assuming the deformation of the network chains is affine. Simplifications of the model are derived for specific cases involving homopolymer networks, copolymer networks, pure solvents, and binary solvent mixtures. With reasonable assumptions, the number of polymer-solvent interaction parameters that must be determined by experiments can be reduced to two effective parameters (θ1 and θ2), which describe the net interactions between water/copolymer (θ1) and ethanol/copolymer (θ2), respectively. Experimental measurements of the swelling of random copolymer networks of n-butyl acrylate and 2-hydroxyethyl acrylate in water, ethanol, and a 100 g/L ethanol/water mixture are utilized to validate the model. For a random copolymer network, θ1 and θ2 can be obtained by fitting the three-component FR model to equilibrium swelling data obtained in the pure solvents. Predicted solvent volume fractions for swelling in water-ethanol mixtures obtained by inserting fitted values of θ1 and θ2 into the four-component FR model are in reasonable agreement with experimental measurements.

Temperature dependence of creep compliance of highly cross-linked epoxy: A molecular simulation study

We have used molecular dynamics (MD) simulations to study the effect of temperature on the creep compliance of neat cross-linked epoxy. Experimental studies of mechanical behavior of cross-linked epoxy in literature commonly report creep compliance values, whereas molecular simulations of these systems have primarily focused on the Young’s modulus. In this work, in order to obtain a more direct comparison between experiments and simulations, atomistically detailed models of the cross-linked epoxy are used to study their creep compliance as a function of temperature using MD simulations. The creep tests are performed by applying a constant tensile stress and monitoring the resulting strain in the system. Our results show that simulated values of creep compliance increase with an increase in both time and temperature. We believe that such calculations of the creep compliance, along with the use of time temperature superposition, hold great promise in connecting the molecular insight obtained from molecular ...