Atomic level picture of stress relaxation in polymer melts

We have been developing a physical picture on the atomic level of stress relaxation in polymer melts by means of computer simulation of the process in model systems. In this article we treat a melt of freely jointed chains, each with N = 200 bonds and with excluded-volume interactions between all nonbonded atoms, that has been subjected to an initial constant-volume uniaxial extension. We consider both the stress relaxation history σ(t) based on atomic interactions, and the stress history σ_e(t; N_R) based on subdividing the chain into segments with N_R bonds each, with each segment regarded as an entropic spring. It is found that at early times σ(t) > σ_e(t; N_R) for all N_R, and that, for the remainder of the simulation, there is no value of N_R for which σ(t) = σ_e(t; N_R) for an extended period; by the end of the simulation σ(t) has fallen just below the value σ_e(t; 50). The decay of segment orientation, 〈P₂(t; N_R)〉, and of bond orientation 〈P₂(t; 1)〉, is computed during the simulation. It is found that the decay of the atom-based stress σ(t) is closely related to that of 〈P₂(t; 1)〉. This result may be understood through the concept of steric shielding. The change in local structure of the polymer melt during relaxation is also studied. © 1998 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 36 : 143–154, 1998     

Molecular model of phenolic polymer dissolution in photolithography

The resolution of photolithographic processes has advanced to the point that difficulties, such as line‐edge roughness, associated with phenomena occurring at molecular length scales are becoming important. In order to control these phenomena, it is necessary to understand them. To that end, a numerical model has been used to simulate the dissolution of phenolic polymers in aqueous base. The simulation applies the Critical Ionization Model to a rectangular‐lattice representation of the polymer matrix. The model has been adapted to describe the dissolution process that is responsible for photoresist function. Both continuum and molecular versions of the model are presented. The Continuum Model yields dissolution profiles that approximate the contours of the calculated spatial variations in chemical blocking (blocking profile). An algorithm has been developed to connect individual cells to form polymer chains, and to fill the lattice in a way that produces a Gaussian chain length distribution. The model employs only a single adjustable parameter, the time‐step correction factor (assuming one can measure the probability of ionization once a site encounters the developer). The Molecular Model predicts a dissolution rate that decreases non‐linearly with respect to degree of chemical blocking, as is observed experimentally. Dissolution profiles can be generated with the Molecular Model based either on this calculated dependence of the dissolution rate on blocking fraction or from direct application of the model to a blocking profile. The probabilistic nature of the model introduces edge roughness of the same degree as that observed experimentally. © 1999 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 37: 2103–2113, 1999     

Salt concentration dependence of the mechanical properties of LiPF6/poly(propylene glycol) acrylate electrolyte at a graphitic carbon interface: A reactive molecular dynamics study

This reactive molecular dynamics study explores the salt concentration dependence of the viscoelastic and mechanical failure properties of a poly(propylene glycol)/LiPF₆‐based solid polymer electrolyte (SPE) at a graphitic carbon electrode interface. To account for the finite‐size effect of interface‐confined SPE films, the properties of two distinct film thicknesses are compared with the respective bulk properties. Additionally, the effect of uniaxial compression in the interface‐normal direction on free energy profiles of Li‐ion SPE‐desolvation is studied. © 2018 The Authors. Journal of Polymer Science Part B: Polymer Physics Published by Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2018 , 56, 718–730     

Synthesis and characterization of a new network polymer electrolyte containing polyether in the main chains and side chains

A new network polymer electrolyte matrix with polyether in the side chains and main chains was synthesized by the azo-macroinitiator method and urethane reaction. The macroinitiator, polymer and network polymer were confirmed by Fourier-transform infrared (FT-IR) spectroscopy and ¹H NMR. FT-IR was also used to study the environment of lithium ions doped in these network polymer electrolytes. Three important groups are considered: N-H, carbonyl, and ether groups. The thermal properties of the polymer electrolytes were measured by differential scanning calorimetry and thermogravimetric analysis. The T_g value of this polymer is less than that of a general comb-like polymer. Added lithium ions interact with the oxygen atoms on ether groups, causing the T_g of the polymer electrolyte to increase. Moreover, the interaction between lithium ions and ether groups decreases the decomposition temperature of the polymer. The conductivity measured by AC impedance reached a maximum of 10⁻⁴ S cm⁻¹. A plot of conductivity vs. temperature fit the Vogel-Tamman-Fulcher equation, indicating that ionic mobility in this network polymer electrolyte is coupled to segmental chain movements.     

Collapse‐expansion transition of elastic shell induced by grafted polymer chains

A coarse‐grained model for an elastic shell grafted with polymer chains is investigated by molecular dynamics methods. With increasing the number of grafted polymer chains (GPCs), it is found that the conformation of the shell undergoes from expansion to collapse and back to the expansion. By varying the density of the GPCs, the phase transition of the elastic shell can be successfully controlled at moderate bending energy of the shell and at moderate binding energy between the shell and GPCs. Furthermore, the self‐assembly structures of the GPCs are also affected by the elastic shell in certain conditions. In the case of a few GPCs on the shell, the chains tend to be adsorbed on the shell surface unfolded at high value of bending energy. However, when the bending energy is small, the chains can be folded several times easily. This may be an important step toward a deeper understanding of how to control the microstructure in the production of biocomposites. © 2012 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys, 2012     

Polyaromatic-terminated iron(ii) clathrochelates as electrocatalysts for efficient hydrogen production in water electrolysis cells with polymer electrolyte membrane

1. Download : Download high-res image (182KB)
2. Download : Download full-size image

     

Molecular simulation of structure formation and rheological properties of mixtures of telechelic and monofunctional associating polymer

We study the phase behavior, the characteristic times, and the rheological properties under the steady shear flow of the mixtures consisting of telechelic and monofunctional associating polymers by a coarse‐grained molecular dynamics simulation. The mixtures form the transient networks, the closely packed spherical micelles, and the wormlike micelles. We confirm the molecular origins of the several characteristic times of the mixtures. The dependencies of the characteristic times on the composition ratio between telechelic and monofunctional associating polymers show good agreement with reported experimental results. Under the steady shear flow, the mixtures show the shear thinning induced by the change of the spatial configuration of the micelles. The telechelic associating polymers especially play an important role in connecting the micelles at the shear thinning regime and enhance the steady shear viscosity. Furthermore, at the wormlike micellar region, the mixtures show the second shear thinning initiated by the transformation of the association conformation of the telechelic associating polymers.     

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 t_w 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 strain–stress demonstrates saturation of the Young's modulus at t_w = 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 t_w > 60 ns. It is shown that both strain–stress and shear viscosity curves agree with experimental data.     

Effects of ethylene glycol dimethacrylate as cross-linker in ionic liquid gel polymer electrolyte based on poly(glycidyl methacrylate)

In this study, modified poly(glycidyl methacrylate)-based films for gel electrolyte were prepared by an in situ UV photopolymerization technique. The effects of adding ethylene glycol dimethacrylate (EGDMA) to the polymer host were studied through X-ray diffraction analysis and differential scanning calorimetry. The results from Fourier transform-infrared spectra indicate complete polymerization among the monomers. The addition of EGDMA to the formulation of gel polymer electrolyte increased the loading of 1-butyl-3-methylimidazolium bis(fluoromethylsulfonyl)imide up to 200 wt.% with the highest value of 8.2 × 10^?4 S cm^?1. All the gel polymer electrolyte membranes obeyed the Arrhenius law.     

Sorption and diffusion of methanol and water in PVDF‐g‐PSSA and Nafion® 117 polymer electrolyte membranes

Sorption and diffusion properties of poly(vinylidene fluoride)‐graft‐poly(styrene sulfonic acid) (PVDF‐g‐PSSA) and Nafion® 117 polymer electrolyte membranes were studied in water/methanol mixtures. The two types of membranes were found to have different sorption properties. The Nafion 117 membrane was found to have a maximum in‐solvent uptake around 0.4 to 0.6 mole fraction of methanol, while the PVDF‐g‐PSSA membranes took up less solvent with increasing methanol concentration. The proton NMR spectra were recorded for membranes immersed in deuterated water/methanol mixtures. The spectra showed that the hydroxyl protons inside the membrane exhibit resonance lines different from the resonance lines of hydroxyl protons in the external solvent. The spectral features of the lines of these internal hydroxyl groups in the membranes were different in the Nafion membrane compared with the PVDF‐g‐PSSA membranes. Diffusion measurements with the pulsed field gradient NMR (PFG‐NMR) method showed that the diffusion coefficient of the internal hydroxyl groups in the solvent immersed Nafion membrane mirrors the changes in the diffusion coefficients of hydroxyl and methyl protons in the external solvent. For the PVDF‐g‐PSSA membranes, a decrease in the diffusion coefficient of the internal hydroxyl protons was seen with increasing methanol concentration. These results indicate that the morphology and chemical structure of the membranes have an effect on their solvent sorption and diffusion characteristics. © 2000 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 38: 3277–3284, 2000     

Active layer of the cathode of a fuel cell with a solid polymer electrolyte: The effect of the Nafion concentration on the overall characteristics

A computer-aided simulation of the structure of the active layer of the cathode of a fuel cell with a solid polymer electrolyte (Nafion) is performed under the assumption about equidimensionalness of dimensions of grains of the substrate (with platinum crystallites in them) and grains (agglomerates of molecules) of Nafion. It is analyzed how the Nafion concentration affects principal parameters, which include the specific surface area, in the vicinity of which electrochemical process goes on; the effective ionic electroconductivity, and the effective diffusion coefficient of a gas. It is demonstrated how one can determine the Nafion concentration at which the overall current takes on a maximum value. Dependences of the optimum value of the overall current and the thickness of the active layer and the weight of platinum, which correspond to it, on the Nafion concentration are calculated. It is demonstrated that there in principle cannot exist one individual optimum concentration of Nafion, which is suitable for all techniques used for the preparation of the active layer. The mutual relationship between values of the effective diffusion coefficient of a gas and the effective ionic electroconductivity of Nafion determines the value of the optimum of the Nafion concentration.     

Molecular dynamics simulation study of association in trifluoroethanol/water mixtures

Mixtures of Trifluoroethanol (TFE) and water with different proportions are studied using molecular dynamics simulations. The radial and spatial distribution functions, as well as the size distribution of TFE clusters are obtained from the trajectories. The variation of radial and spatial distribution functions with composition show that the addition of TFE enhances the water structure, but the hydrogen bonds between TFE molecules are broken as TFE is diluted with water. The TFE‐rich solutions have stronger TFE–water hydrogen bonds. The clustering of TFE molecules in low concentration region is attributed to the hydrophobic interactions between CF₃ groups. The distribution of cluster sizes in solution supports these conclusions. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Cathode of a hydrogen-oxygen fuel cell with a solid polymer electrolyte: The effect of the flooding of pores by water on the characteristics of the active layer

A computer model of the active layer of the cathode of a hydrogen-oxygen fuel cell with a solid polymer electrolyte is studied. The active mass of the electrode consists of equidimensional grains of the substrate (agglomerates of carbon particles with platinum particles embedded in them) and a solid polymer electrolyte (Nafion). The flooding by water can be experienced by both the pores in the substrate grains, which facilitate the oxygen penetration into the active layer of the electrode, and the voids between the grains. All possible versions of the flooding of these pores by water are considered. A calculation of the optimum, at a given polarization of the electrode, value of electrochemical activity, the thickness of the active layer, and the weight of platinum is performed. The major parameters of the system are the concentrations of grains of the substrate and solid polymer electrolyte, the size of these grains, the platinum concentration in the substrate grains, the average diameter of pores in the substrate grains, and the polarization of electrodes. The ultimate aim of the work is to estimate how the flooding of pores of the active layer of the cathode by water affects the magnitude of the optimum current, the effective thickness of the active layer, and the weight of platinum.Translated from Elektrokhimiya, Vol. 41, No. 1, 2005, pp. 35–47.Original Russian Text Copyright © 2005 by Chirkov, Rostokin.     

Mass transfer in porous systems, flooded in part with water, of gas diffusion and catalytic layers of the cathode of a fuel cell with a solid polymer electrolyte

Mass transfer in porous gas diffusion and catalytic layers of the cathode of a hydrogen-air fuel cell with a solid polymer electrolyte is considered. The transport processes are considered with allowance made for the partial flooding of porous systems of these layers with water, which forms during the fuel cell operation. The consideration also allows for the influence of the diluent gas present when air oxygen is used as the oxidant. The fraction of water-flooded pores is calculated within percolation theory as a function of structural parameters of the porous system. Conditions leading to the beginning of the gas diffusion layer flooding are presented.     

Nonadiabatic relaxation dynamics of water anion cluster and its isotope effects by ring‐polymer molecular dynamics simulation

We have applied a recently developed hybrid quantum ring‐polymer molecular dynamics method to the nonadiabatic p → s relaxation dynamics in water anion clusters to understand the isotope effects observed in previous experiments. The average relaxation times for (H₂O)₅₀^? and (D₂O)₅₀^? were calculated at 120 and 207 fs, respectively, and are comparable to the experimental results. Therefore, we conclude that nuclear quantum effects play an essential role in understanding the observed isotope effects for water anion cluster nonadiabatic dynamics. The nonadiabatic relaxation mechanisms are also discussed in detail. © 2014 Wiley Periodicals, Inc.     

Effects of polymer adsorption on the effective viscosity in microchannel flows: phenomenological slip layer model from molecular simulations

Molecular simulations (Dissipative Particle Dynamics - DPD) were used to quantify the effect of polymer adsorption on the effective shear viscosity of a semi-dilute polymer solution in microchannel Poseuille flow. It is well known that polymer depletion layers develop adjacent to solid walls due to hydrodynamic forces, causing an apparent wall slip and reduced effective viscosity (increased total flow rate). We found that depletion layers also developed in the presence of hydrodynamically rough adsorbed layers on the wall. Polymer-polymer (steric) repulsion between flowing and adsorbed polymer expanded the depletion layer compared to no-adsorption cases, and the effective viscosity was reduced further. Desorption occurred for higher shear rates, reducing the repulsion effect and shrinking the depletion layers. A phenomenological algebraic model for the depletion layer thickness, including a shear modified adsorption isotherm, was developed based on the simulation data. The depletion layer model can be used together with the effective viscosity model we developed earlier.     

Simulation of phase separation in alcohol/water mixtures using two-body force field and standard molecular dynamics

Standard molecular dynamics simulations have been carried out on pure alcohols and alcohol/water mixtures. A simple atom-atom force field consisting of Lennard-Jones potentials plus coulombic terms over atomic point charges, but without explicit polarization terms, has been specifically fitted to reproduce several experimental properties of the pure alcohols, and has been used for mixtures by developing combination rules with the TIP3P water model. Densities, enthalpies of vaporization, radial distribution functions, self-diffusion coefficients, and rotational correlation functions of the pure alcohols are well reproduced and compare favorably with those from more sophisticated force fields. Some key aspects of the phase behaviour are correctly reproduced by the molecular dynamics simulation, showing a distinct demixing process for the n-butanol/water mixture as opposed to the stability of the t-butanol/water mixtures. The results demonstrate the ability of a molecular dynamics simulation, even in its standard form and with easily accessible time ranges, but with a carefully optimized force field, to simulate and, to a certain extent, predict the properties of binary mixtures.     

电解质水溶液传递性质的布朗动力学模拟研究

在传统布朗动力学的基础上, 考虑流体力学的影响, 并且引入Smart Monte Carlo方法的接受概率, 对电解质溶液进行布朗动力学模拟, 得到不同浓度和温度下KCl溶液中离子间的径向分布函数, 并且与超网链积分方程理论计算结果进行了比较, 同时, 模拟了KCl和NaCl溶液的摩尔电导率. 模拟过程基于电解质溶液的原始模型, 溶剂被看作连续介质, 溶质分子之间的相互作用采用软核加静电的势能函数模型, 长程静电力的处理采用Ewald加和方法. 结果显示, 流体力学的作用对于电解质溶液的结构性质没有明显的影响, 但是对于传递性质的影响显著; 考虑流体力学作用的布朗动力学模拟结果与实验数据吻合良好.