Tr. .mv.cvuar u.- rK-nw m- Ilillillilllllll’l‘lllllmill» 3 1293 01025 0920 This is to certify that the thesis entitled SORPTION OF ORGANIC PENETRANTS BY AMORPHOUS POLYAMIDE presented by MYUNGHOON LEE has been accepted towards fulfillment of the requirements for MASTER degree inPACKAGING Major professor DateélliIL 1. 1994 0-7639 MS U is an Affirmative Action/Equal Opportunity Institution LIBRARY Mlchlgan State Unlverslty PLACE IN RETURN BOX to romovo this checkout from your record. TO AVOID FIN IS rotum on or before date duo. DATE ‘DU E DATE DU E DATE DU E r ' 45?”: FEB 0 “ 21.059 . “.991 FE! 2 MAGIC 2 W911i. W“ L‘flll MSU!.A." .u. .L A - '1 ”.rr SORPTION OF ORGANIC PENETRANTS BY AMORPHOUS POLYAMIDE By Myunghoon Lee A THESIS Submitted to Michigan State University in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE School of Packaging 1994 ABSTRACT SORPTION OF ORGANIC PENETRANTS BY AMORPHOUS POLYAMIDE By Myunghoon Lee Sorption studies involving the sorption of n-propanol by an amorphous nylon ( Nylon 61 /6T ) were carried out as a function of sorbate vapor activity at 23 0 C. Vapor activity levels from 0.035 to 0.91 were investigated to evaluate the concentration dependency of sorption mechanism. Sorption behavior of propanol by Nylon 6I/6T showed distinctive two mode sorption phenomena as a function of vapor activity . At vapor activity levels below a=0.11, equilibrium sorption was achievediwithin a short period of time ( less than 20 hours ) , which can be interpreted as fellowing a Fickian diffusion model. i A Langmuir- Flory-Huggins dual mode sorption model can also be applied at these concentration levels. However , for vapor activities above a=0.11, the sorption process appeared to be non-Fickian and resulted in a lack of equilibrium being attained. ACKNOWLEDGEMENTS I'd like to express my heartfelt thanks to my major advisor Dr. Jack . R . Giacin. His devoted and endless supports enabled me to complete this research. I would also like to thank Dr. Ruben Hernandez and Dr. Eric. A. Grulke for their precious advice and discussion throughout the course of research. I also thank to Dr. Heidi and Mr. Bob Hurwitz for helping me to design and assemble the test system and to provide comfortable circumstance during the research. I am deeply indebted to the Korea Institute of Industrial Design and Packaging , which sent me to study in the School of Packaging. Finally , my lovely wife and two daughters deserve to have my deep thanks and love for their devoted supports during the research. ii TABLE OF CONTENTS LIST OF TABLES LIST OF FIGURES INTRODUCTION LITERATURE REVIEW Characteristics of Amorphous Polymers Properties of Amorphous Polyamide Properties of Penetrant Sorption Mechanisms Ideal sorption and diffusion Non ideal sorption Modified dual mode sorption Factors affecting the sorption process Temperature Crystallinity Orientation Sorption measurement iii Page vi 13 14 14 16 20 22 22 25 27 30 Diffusion and relaxation in glassy polymers Correlation of sorption with penetrant size MATERIALS AND METHODS Polymer film n-Propanol Sorption measurement Density experiment RESULTS AND DISCUSSION Data preparation and prelimenary test Equilibrium sorption Analysis of sorption at vapor activity a < 0.11 Analysis of sorption at vapor activity a > 0.11 CONCLUSION APPENDIX Appendix A Procedure for standard calibration curve constraction Appendix B Electrobalance calibration procedure Appendix C Vapor activity fluctuation BIBLIOGRAPHY iv 31 33 35 35 35 36 39 47 49 53 67 68 70 71 76 LIST OF TABLES Table Title Page 1. Comparison of the properties between Amorphous 12 Nylons and Nylon 66 2. Relationship between rotameter readings and 45 vapor activities 3. Volume fraction vs. vapor activity from Flory- 51 Huggins equation 4. Average vapor activity and standard deviation 52 in each step of vapor activity Figure l. 2. 10. ll. 12. 13. 14. 15. 16. 17. LIST OF FIGURES Title Structure of Nylon 61 / 6T Structural formulas of various alcohols Schematic diagram of electrobalance test apparatus Saturated vapor pressure of n-propanol Standard calibration graph of n-propanol at 23 0 C Illustration of sample flushing by nitrogen at 23 0 C Evaluation of the reproducability of the sorption system at a = 0.06 Interval sorption of n-propanol at a = O - 0.035 Interval sorption of n-propanol at a = 0 .035- 0.050 Interval sorption of n-propanol at a = 0.050 - 0. 80 Interval sorption of n-propanol at a = 0.080 - 0.110 Interval sorption of n—propanol at a = 0.110 - 0.295 Interval sorption of n-propanol at a = 0.295 - 0.390 Interval sorption of n-propanol at a = 0.390 - 0.490 Interval sorption of n-propanol at a = 0.490 - 0.840 Interval sorption of n-propanol at a = 0.840 - 0.910 Successive interval sorption of n-propanol into Nylon 61/6T at 23 OC(a=0-0.9l) vi Page 10 13 41 42 43 56 56 57 57 58 58 59 59 60 6O 61 18. 19. 20. 21. 22. 23. 24. 25. 26. 27. 28. Successive interval sorption of n- propanol into Nylon 61/6T at 23 OC(a=0-0.110) Experimental sorption isotherm for n-propanol in Nylon 6I/6T at the low vapor activities ( a< 0.11 ) Flory - Huggins plot when x = 1.6 from the equation Ina1=1nV1+ v2 + xv22,v,=1- V3 Plot of Mt/Moo vs. t 1/2 for n-propanol sorption by Nylon 6I/6T, at 23 0C ( a = 0 - 0.035) Plot of Mt/Moo vs. t 1/2 for n-propanol sorption by Nylon 6l/6T, at 23 0C ( a = 0.035 - 0.050) Plot of Mt/Moo vs. t 1/2 for n-propanol sorption by Nylon 61/6T, at 23 OC ( a = 0.050 - 0.080) Plot of Mt/Moo vs. t 1/2 for n-propanol sorption by Nylon 6I/6T, at 23 0C ( a = 0.080 - 0.110) Comparison between experimental and calculated data at vapor activity a =0.035 Deviation from Equation(21) due to non - Fickian sorption for sorption of n-propanol by Nylon 6I/6T at 23 0C( a = 0.11 - 0.295): Graphic estimation of Fickian Moo value Plot of Mt/Mco vs. t 1/2 for n-propanol sorption by Nylon 6I/6T, at 23 0C; a = 0.11 - 0.295, Equation (21) for Moo = 0.275 and t 1/2 = 84 sec. Plot of Mt/Moo vs. t 1/2 for n-propanol sorption by Nylon 6I/6T, at 23 0C; a = 0.11 - 0.295 ; points observed curve,Equation(21) for Moo = 0.275 and t 1/2 = 84 sec. vii 61 62 62 63 63 65 65 66 66 29. 30. 32. 33. 35. 36. 37. Distribution of vapor activity at Distribution of vapor activity at Distribution of vapor activity at Distribution of vapor activity at Distribution of vapor activity at Distribution of vapor activity at Distribution of vapor activity at Distribution of vapor activity at Distribution of vapor activity at a = 0.035 a = 0.050 a = 0.080 a = 0.1 10 a = 0.295 a = 0.390 a = 0.490 a = 0.840 a = 0.910 viii 71 71 72 72 73 74 74 75 INTRODUCTION Over the past decade , the development and use of high barrier plastic materials has increasingly replaced traditional glass and metal containers for food and beverage packaging. The advantages of plastic packaging are numerous and include : low cost . light weight , a wide range of mechanical properties ,transparency. flexibility , direct food contact , and general consumeripreference because of convenience and microwaveability ( Salame , 1986). However , with the use of plastic packaging there are concomitant concerns related to product / package interactions. This is a broad based topic that includes transport of gases and organic vapor of low molecular weight : (i) from the product through the package , as well as (ii) from the environment through the package to the product. The specific mass transport ( permeability ) process may be described as a function of penetrant concentration , temperature , and time to equilibrium, as well as the composition of the penetrant / polymer system. Therefore , polymer structure , free volume , chain stiffness or segmental mobility , and the availability of specific sites of interaction in the polymer,plus the physicochemical characteristics of the penetrant (gas, water vapor , or organic vapor) determine the mode and mechanism of sorption 1 2 and thus . the resultant transport and mechanical properties of the polymer ( Blatz , 1989 ). The dual mode sorption model has been widely used to describe the solubility of gases in glassy polymers . and in glassy , polar polymers as well. This model assumes that the solute molecules in the glassy polymer consist of non-specifically absorbed and specifically adsorbed species , which are in dynamic equilibrium in the medium ( Michaels et al., 1963 ) The solubility of the absorbed species is represented by Henry's law , and the solubility of the adsorbed species is described by a Langmuir type adsorption isotherm. The Langmuir sorption is believed to occur at specific sites , usually considered to be units of stable free volume present in the amorphous polymer structure. Local equilibrium between the absorbed and adsorbed populations is maintained throughout the polymer matrix. The total amount of solute sorbed by both mechanisms is : CH' b C=CD+CH=KDP+——— (1) l + bp where C is the total concentration of sorbed solute in the polymer , CD and CH are the solubilities due to absorption ( Henry's Law ) and Langmuir type adsorption respectively , K1) is the Henry's law 3 constant . b is the hole affinity constant . CH' is the hole saturation constant and p is the pressure. Recently , Hernandez et a1. ( 1991 ) have studied the sorption of water vapor into an amorphous polyamide ( Nylon 61 / 6T ). The authors have applied a dual mode sorption model , based on Langmuir and Flory - Huggins equations , to describe the sorption process. The total amount of water sorbed by the polymer was described by the summation of a Langmuir type association , and by a solution component mode given by the Flory-Huggins model , which is given by: V1 = V1L + V1FH (2) Where V1 is the total volume fraction of water within the polymer , and the superscript L and FH refer to Langmuir and Flory - Huggins water sorption contributions , respectively. These contributions are expressed as : K a V1L= (3) 1 + Ba and a=exp[1nV1FH+(l~V1FH)+x(1-V1FH)2] (4) Where a is water activity , x is the Flory - Huggins interaction parameter , and K and B are parameters of the Langmuir equation ( Hernandez et a1., 1991). The behavior of oxygen solubility within the polymer / water system 4 was related to V114 , according to the following expression. V0:2 = V* - FV]L (5) where V02 is the solubility of oxygen in the polymer as a function of a , V* is the solubility of oxygen in the dry polymer , and F is a factor that relates sorption values of oxygen at dry conditions and the fraction volume of water described by the Langmuir sorption mode , when water activity equals one. The fact that 80 % of the total oxygen dissolved by the polymer at dry conditions was displaced by molecules of water associated with active sites of the polymer matrix , indicated the importance of these active sites in the mechanism of the solubility of oxygen within the Nylon 61 / 6T bulk phase. These results also suggested that molecular size may be an important factor in determining the final equilibrium sorption values in the three component system of polymer , water and oxygen ( Hernandez et al. , 1991 ). Based on the results described above , it can be proposed that a similar model can be applied to the binary system of polymer and organic penetrants ,which has not yet been reported. A quantitative relationship between molecular weight of sorbates and active sites concentration in polyethylene was reported by Gedraitite et a1. ( 1989) , who studied the sorption by polyethylene , of selected low molecular weight compounds dissolved in n - hexane , and in isopropanol. D In addition to the limited amount of data describing the effect of moisture content on the mass transport and solubility of oxygen in hydrophilic polymers , no references have been found in the literature which considered the importance of size distribution of sorption sites in determining the final equilibrium sorption values for the sorption of organic penetrants by hydrophilic polymers , such as the amorphous polyamide. There is also a paucity of data describing the effect of organic vapor concentration on the solubility and diffusion of organic penetrants in hydrophilic polymers , and the importance of active binding sites associated with the polymer matrix in determining the final equilibrium sorption value for the binary system of polymer / organic vapor. The proposed studies will provide for a better understanding of the parameters contributing to the diffusion and sorption of organic penetrants in hydrophilic polymer membranes . In terms of practical application , the data obtained describing the relationship between sorption parameters and the molecular weight of the sorbate molecule may allow calculation of the solubility of a penetrant in the polymer matrix from knowledge of the penetrants solubility in a contact phase , and its molecular weight. This could further provide a means of designing a barrier structure for a specific end use application. For example , in the packaging of a juice product , when there is concern for flavor loss due to sorption _. or sorption of ingredients from a 6 pharmaceutical preparation . which could result in loss of product efficacy. Objectives 1. Determine the equilibrium sorption isotherm for n-propanol in an amorphous polyamide ( Nylon 61 / 6T ) at water activity 2 0. 2. Provide an empirical framework to study the sorption of organic vapor by a glassy polymer , and the relationship of Langmuir type association ( active site binding) , to the solubility of organic penetrants in hydrophilic polymers. 3. Evaluate the appropriateness of the Langmuir-Flory-Huggins dual mode sorption model to describe the sorption isotherm of n-propanol in the amorphous polyamide( Nylon 61 / 6T ). LITERATURE REVIEW Characteristics of Amorphous Polymers A vast majority of well-known polymers are characterized by their structural differences. The conformation of flexible chain molecules in the melt or glassy state is determined both by the local position and orientational distribution of chain segments belonging to the same chain ( intramolecular correlations ) or belonging to different chains (intermolecular correlations ) (Wendorff , 1987 ). Polymers differ from other low molecular weight organic or inorganic materials , in that they are composed of long chain molecules , usually containing a huge number of atoms. This gives rise to a very large number of internal degrees of freedom for flexible chain molecules ( Blumstein , 1978 ). The chain is represented by a line in space , the curvature of which depends upon the persistence length A1. The parameter(Ai) is thus a measure of the chain stiffness. It corresponds to the distance along the chain over which correlations in the orientation of successive chain units extend. The mean square end-to-end distance < h2 > is given by

/2AiL=[1-(l/x)][1-exp(-x)] (6) where L is the contour length of the chain and x = L /Ai. A worm- like chain becomes a flexible Gaussian chain for x —> 00 , and a rigid rod- like 7 8 chain for x —> 0 ( Wendorff . I987 )- Flexible chain molecules have a persistence length of the order of 1 or 2 nm. A flexible chain in a random configuration occupies only a small portion of the space it pervades. The density of an isolated chain molecule is typically of the order of l % of that of the condensed state , and it decreases with increasing chain length ( Wendorff , 1987 ). Intermolecular forces require , however , that the empty space be filled in the condensed , amorphous state. This may be achieved in various ways as described below. i (l) The chains may collapse , a process which is known to occur in solutions with poor solvents ( Williams et al. , 1981) ; (2) The chains may aggregate in bundles in which neighboring chains or chain segments are parallel to each other as in a nematic liquid crystalline phase ( Pechhold et al. , 1979 ) ; (3) The chains may remain in their random configuration and many different chains are able to pervade the space taken up by the reference chain. This leads to a highly entangled system ( Gennes , 1979 ). The distribution of the centers of the chains as well as their orientation is as close to random , as in a real gas. Macroscopic properties of amorphous polymers are known to depend heavily on the local spatial and orientational distribution of chain segments , belonging to the same or to a different chain (Wendorff, 1987). 9 Properties of Amorphous Polyamide Polyamides containing aromatic residues are normally very high melting, highly insoluble semi-crystalline materials , which are difficult to process. In recent years , a number of more easily processed amorphous polyamides , based on iso— and terephthalic acids have been claimed in the patent literature. In general , amorphous polyamides can possess a wide range of good mechanical properties , including tensile and impact strength , and a high specific modulus , which are desirable for engineering applications ( Dolden , 1976 ). The amorphous polyamide ( Nylon 61 / 6T ) which was developed by the E.I.Du pont De Nemours & Co., Ltd., has a random structure based on hexamethylenediamine and a 70/ 30 mixture of isophthalic and terephthalic acids ( Fig. l ), and when dry has a glass transition temperature (Tg ) of 125 0 C ( Blatz ,1989). It has good processability by primary extrusion and molding operations , as well as secondary processing such as thermoforrning. Film has been produced by both the blown film and flat cast processes. Coextruded film and sheet can be produced using three to five layers , at thicknesses from 2 mils to 60 mils. IO NHz CH2 - (CH2)4 - CH2 NHz -> Hexamethylene diamine + COzH C02” | + | \ COzH | COzH lsophthalic acid Terephthalic acid (70 % ) ( 30 % ) Figure 1. Structure of the Nylon 61 / 6T Nylon 61 / 6T also has good tensile strength , elongation and tensile modulus. These properties are similar to those for the other resins used for quality packaging applications. The tensile strength of amorphous polyamide is higher than nylon 66 dry , but decreases in a manner similar to nylon 66 at 50 % RH. The tensile modulus , a measure of stiffness , shows an atypical trend with increasing humidity. As expected , the modulus of nylon 66 decreases substantially as the humidity increases , however the modulus of amorphous polyamide 11 increases slightly. Thus . it maintains its stiffness at high humidities ( Blatz , 1989 ). The barrier properties of amorphous polyamide approaches the levels of polyvinylidene chloride , ethylene vinyl alcohol copolymer and the extrudable acrylonitrile copolymers. As was previously indicated, water vapor has an unexpected effect on the permeability of amorphous polyamide. Whereas the oxygen permeability of nylon 66 film increases on exposure to water vapor , that for the amorphous polyamide film decreases. Thus , a film of amorphous polyamide has a significantly better barrier as the humidity increases ( Blatz , 1989 ). The chemical resistance is usually not as good as the semi-crystalline resins , such as PET or 6 and 66 Nylon. The resin is unaffected by dilute acids and bases , but is attacked by acetic acid ( Blatz, 1989 ). It is also attacked by the common low molecular weight alcohols such as ethyl alcohol and solutions containing more than 10 % ethyl alcohol. The amorphous nylon is resistant to both aliphatic and aromatic hydrocarbons , to chlorinated solvents , ketones , and esters. However , methylene chloride will cause swelling ( Blatz, 1989 ). As with all glassy polymers , as the melt is cooled down and reaches the glass transition temperature ( Tg ) , the molecular motions are frozen at that temperature , which results in a void content that is higher 12 than that in a polymer with low Tg . For an amorphous glassy polymer . the larger the difference in temperature between the Tg and the test temperature the greater the amount of unrelaxed volume ( sometimes called the free volume ). Table 1 summarizes selected mechanical and barrier properties of amorphous nylon ( Nylon 61 / 6T ) and nylon 66. Items Unit Amorphous Nylon* Nylon 66 Tensile strength K Psi 10.0 8.0 Elongation % 50.0 70.0 Tensile modulus K Psi 330 340 Heat deflection T. @ 66 psi 26] 0 F (128 0 C) 365 0 F (185 0 C) of packaging resin @ 264 psi 248 0 F (120 0 C) 151 0 F (66 0 C) Oxygen Dry 3.3 2.5 Permeability (1) 50 % RH 1.4 5.0 CO2 8.0 10.0( 50 % RH) Permeability (2) Water vapor 5.0 5.7 Permeability (3) * SELAR® PA of Du Pont Co., (1) CC-mil / 100 SQ. IN. / 24 Hours-Atm. at 30 0 C (2) CC-mil / 100 SQ. IN. / 24 Hours- Atm. at 80 % RH, 30 0 C (3) g-mil / 100 SQ. IN. / 24 Hours- Atm. at 90 % RH, 23 0 C Table 1. Comparison of mechanical and barrier properties of amorphous nylon* and nylon 66 13 Properties of penetrant : l - Propanol Alcohols may be viewed as alkyl derivatives of water and are indicated by the formula R - OH. They are structurally similar to water , but have one of the hydrogens replaced by an alkyl group. The reactions of alcohols are those of the -OH functional group , known as the hydroxyl group . Figure 2 lists both the IUPAC and common names for a variety of simple alcohols , which include the 1 - propanol used in this study. CH3OH CH3CH20H CH3CH2CH20H methanol ethanol l-propanol methyl alcohol ethyl alcohol n-propyl alcohol CH3CHCH3 CH3CHCH20H CH3 | | OH CH3 CH3COH l CH3 2-propanol 2-methyl-l-propanol 2-methyl-2-propanol isopropyl alcohol isobutyl alcohol tert-butyl alcohol Figure 2. Structural formulas of alcohols. 14 Sorption Mechanisms Organic vapors can exhibit concentration dependent mass transport and sorption processes. The permeant vapor pressure and the type of vapors that come in contact with the package will also determine the magnitude of sorption and permeation into and out of polymeric packaging systems. The study of the solubility of gases in polymers has a twofold objective : (1) it permits the establishment of correlations of gas solubilities , with readily available gas molecular parameters ; and (2) it also supplies information on the morphology of polymers ( Vieth ,1966 ). Ideal sorption and diffusion In ideal gas - polymer systems , both the solubility and diffusion coefficients are constant at any given temperature , and so the permeability is also a constant. Permeation or transport through a polymer film can be described in terms of its component parts by Equation (7): P = D S (7) where S is the solubility coefficient , and D is the diffusion coefficient. S characterizes the amount of permeant that can be dissolved into the 15 polymer under the given vapor pressure . and D describes the rate at which the permeant molecules are advancing through the barrier film. For a simple permeation process , the sorption and desorption steps are described by the assumption of Henry's law , which relates the concentration of the penetrant in the polymer , to the vapor pressure in equilibrium with the polymer. The partial pressure of the penetrant is further related to the penetrant concentration in the gas phase through the ideal gas law. Generally , the application of the ideal gas law is justified since the concentration of the diffusant in the gas phase is very low. The diffusion step is described by Fick's first (8) and second laws (9) of diffusion ( Hopfenberg , 1973 ). dc F = - D(c) (8) 6x dc 6] 6c\ = —| D(C)- | (9) dt 6x \ 6x) where F is the flux or the rate of transfer of penetrant per unit area , expressed as mass of diffusant per unit area per time ;c is the concentration of the penetrant in the film , expressed in the same unit of mass of diffusant per unit of volume or mass of the polymer. D is the mutual diffusion coefficient , in ( length )3/ time ; t is time ; and x is the length in the direction in which transport of the penetrant 16 molecules occurs. To obtain the flux (F) or the diffusion coefficient (D) from Equations (8) or (9) , initial and boundary conditions associated with the experimental method are needed , and the expressions solved to give the desired values. It should be recognized that when the diffusion coefficient is calculated using Equations (8) and (9) , only approximate values will be obtained. More accurate estimations of this parameter can be carried out by using , for example , a nonlinear maximum likelihood sequential method , based on the Gauss linearization method ( Beck and Arnold, 1977). The partial pressure of the penetrant is further related to the penetrant concentration in the gas phase through the ideal gas law. Application of the ideal gas law is justified since the concentration of the diffusant in the gas phase is , in general , very low ( Hernandez et al. , 1986 ). Non-ideal sorption The sorption of a low molecular weight compound proceeds by two steps. Firstly the sorbate penetrates into the relatively ordered polymeric substance , forming a true solution. The concentration of the homogeneously dissolved compound is related to the concentrations outside the polymer by Henry's Law. Secondly the dissolved compound is 1 7 reversibly sorbed by certain centers . which are the zones of destruction of the short range order. The total sorbate concentration in the polymer will be the sum of the concentrations of the homogeneously dissolved compound. and of that sorbed by the sorption centers ( Shlyapnikov and Mar'in , 1987 ). For systems in which the solubility doesn't conform to Henrys law , both the diffusion and the solubility parameters are concentration dependent. The solubility coefficient often is essentially constant at low vapor activities for the more volatile vapors, and only the diffusion coefficient exhibits significant concentration dependence ( Rogers et a1. , 1960 ). This concentration dependence can be understood in the context of the free volume theory ( Cohen and Tumbull , 1959 ; Fujita , 1961; Kwei and Wang, 1972 ; Peterlin , 1975 ) , which gives a relationship between the diffusion coefficient D and the fractional free volume f : D=Aexp(-B/f) (10) where A is a frequency factor and B is a measure of the minimum hole size for the jump process ( Choy et al. , 1984 ). Fujita ( 1961 ) further assumed that f = fo + B Cs (11) where Cs is the mass fraction of the penetrant and 8 denotes the effectiveness of the penetrant molecule for increasing the free volume of the polymer. Substituting Equation ( l 1) into Equation (10) and assuming BCs << fo , Equation (10) becomes 18 D=Doexp(st) (12) where y = 86/ to2 (13) Peterlin ( 1975 ) has obtained an alternative expression for y in terms of the fractional free volume of the penetrant , but the physical meaning is essentially the same. Equation (13) shows that the concentration coefficient 7 depends on fo and f3 ; i.e.,y is a function not only of the fractional free volume but also of the interaction between the penetrant and the polymer. For the two penetrant-polymer systems , methylene chloride-HDPE and methylene chloride-LDPE , B does not seem to vary significantly with drawing , so the change in y is determined largely by the I / fo2 factor , which gives rise to a substantial increase in y with increasing draw ratio ( Choy et al. , 1984 ). The idea that the molecules of the penetrant present in a polymer are of two different types (those mobile , homogeneously dissolved , and those immobile , sorbed in certain centers ) , was called the dual mode sorption model. This model has been widely applied to describe nonideal sorption , which can explain the irregularities of diffusion and sorption of low molecular weight compounds ( i.e. , additives ) in glassy polymers ( Paul and Koros , 1976 ). The magnitude of the negative enthalpies of sorption reported by 19 Meats ( 1954 ) for neon . nitrogen . oxygen and argon in glassy polyvinyl acetate and by Barrer et al. ( 1957 ) for organic vapors in ethyl cellulose were inconsistent with the sorption theories of rubbery systems . and led Barter et al. to suggest a two-mode , concurrent sorption mechanism for glassy polymers , namely ordinary dissolution and 'hole' filling. In glassy polymers , substantial deviations from linearity are often observed in sorption isotherms at penetrant pressures above 1 atm. Barrer et al. ( 1957 ) first suggested that two distinct mechanisms may be operating in the sorption of low molecular mass species by glassy polymers. They observed isotherms for the sorption of C4 and C5 isomers in glassy ethyl cellulose that curved towards the pressure axis , and appeared to become asymptotically linear at high penetrant pressures. These results were interpreted in terms of a "dual mode" sorption model , which assumes a combination of Langmuir - type trapping within preexisting vacancies , plus "true" Henry's law solution. Quantitatively , this may be written : C'H P C = 1(1) p + -——— (14) l + b p where C is the total sorbate concentration , K1) is the effective Henry's law solubility constant , Cu is the Langmuir "capacity factor", b is the Langmuir site affinity parameter , and p is the gas pressure. 20 Michaels et al. ( 1963 ) also reported the sorption of helium . nitrogen . oxygen , argon and methane in glassy amorphous , and glassy semi-crystalline polyethylene terephthalate. For penetrant partial pressures up to 1.0 Henry's law was followed. However , carbon dioxide at 25 0 C and 40 0 C , and ethane at 25 0 C, in the same pressure range , deviated from Henry's law. The carbon dioxide isotherms prompted Michaels et al. ( 1963 ) to propose a two - mode sorption model of ordinary dissolution and adsorption in microvoids ( 'holes' ) , for gas sorption in glassy amorphous polymers. Michaels and Bixler( I961 ) attempted to adapt the Flory — Huggins theory of polymer solution to the sorption of simple gases, picturing the process as first condensation of the gas , followed by mixing of the polymer and liquid penetrant. Modified dual mode sorption model Recently , Hernandez et al. ( 1991 ) also proposed a modified dual mode sorption model , which described the sorption of water vapor by an amorphous polyamide at 23 0 C. They modified Eqn. (14) by using the Flory-Huggins equation to describe non-specific solution , rather than Henry's Law. This modification allowed the model to fit over the 21 activity range . 0 < a < l . Although other choices for the solution model are possible , they applied the Flory - Huggins Equation (15) because of its simplicity. Ina =an1+ V2 + )(V22 (15) where x is the interaction parameter. The non-specific sorption term is nonlinear , the complete model can have an inflection point and should predict clustering somewhere over the range of activity values ( Hernandez et al., 1990 ). The modified dual mode sorption model is expressed in terms of volume fractions and solute activity: V1 = VlFH + V1L (2) where V 1FH refers to the Flory - Huggins contribution to the solute volume fraction , and V1L is Langmuir volume fraction contribution. Since Equation (2) is nonlinear , it is convenient to determine the value for V 1FH by numerical methods , such as Newton - Raphson technique and Equation (2) becomes: K a V1: FH(a,x)+ (16) 1+ Ba Factors affecting sorption behavior Temperature The temperature dependence of the sorption parameters have been studied for a number of penetrant / polymer systems. The K1) and b terms are generally decreasing as a function of T , and appear to be well correlated at a given temperature with the penetrant Lennard - Jones energy parameter E/K values ( Pace and Datyner , 1980). The data for C‘ H are more ambiguous. Some results suggest that OH is independent , or varies only weakly with temperature , and others suggest that OH is a strongly decreasing function of T(Pace, 1980). However, it is reasonable to assume that the immobile species is contained in microvoids formed by local segmental motions in the rubbery state which have become ”frozen" in the glass. This frozen void fraction must be , to a first approximation , at least independent of T to explain the low thermal expansivity of the glass , which is comparable to that of the crystalline state. Hence C'H should also be essentially independent of T ( Pace and Datyner , 1980 ). The glass transition temperature ( T8) of any amorphous substance , whether polymeric or not , is defined as the point where the thermal expansion coefficient undergoes a discontinuity. In polymers , there 22 23 expansion coefficient undergoes a discontinuity. In polymers . there may be more than one discontinuity associated with the thermal expansion coefficient term. The largest discontinuity is usually associated with the loss of the molecular mobility . which permits configurational rearrangements of the chain backbones ( Ferry , 1970 ). Below the Tg of a polymer , not enough energy is provided to produce the Micro-Brownian motion of chain segments 20-50 carbon atoms in length , and the chains are locked into the position created from the processing history. Sorption and diffusion in polymers at temperatures below the polymer's Tg consist of more complex behaviors , and are not well understood at this time ( Hernandez , 1984). Above the Tg of a polymer , it seems likely that hole filling is an important sorption mechanism , which is supported by the fact that the low-pressure ( Henry's law ) solubility for simple gases varies smoothly on going through Tg , although the apparent heats of solution are different and constant in the glassy and rubbery regions. Meares ( 1957 ) proposed the following idealised concept of the gaseous sorption and diffusion process in polymers at temperatures above and below Tg , based on the relative activation energies for diffusion and heats of solution for polyvinyl acetate in the two temperature regimes , and the zone 24 of activation theory. Above Tg , gas molecules which dissolve in the polymer must create their own "holes" by separating the interchain polymer contacts. The penetrant then diffuses through the polymer matrix along cylindrical voids created by the synchronized rotation of polymer segments about the C - C bonds ( Meares , 1957). Below Tg , the polymer consists of regions of densely packed and arranged chains which have limited freedom for rotation , separated by less dense regions of disordered chains that form the 'holes' into which the gas sorbs. The gas molecules diffuse between these 'holes' by the slight compressing of localized chains in the dense regions enabling the gas molecules to pass through. This compression in the glassy state does not create the long cavities common to the rubbery state, indicating that the zone of chain activation is much larger in the rubbery state. This zone of activation was found to be essentially independent of the size of the penetrating molecule ( Meares , 1954 ). Above the Tg the polymer is assumed to be in thermodynamic equilibrium , while below Tg the hole concentration is that which was frozen in at Tg , provided the glass is formed by slow cooling from the rubbery region. One therefore expects C'H=Coexp(—Hh/RT), T>Tg (17a) C'H:Coexp(’Hh/RTg), T 0.11 ). The rate of sorption in the initial stage is principally controlled by the Fickian process and therefore, may be used to estimate a true diffusion coefficient, even when the total sorption does not follow the Fickian model ( Berens , 1977) . Furthermore , since the second sorption stage is almost entirely relaxation - controlled , the kinetic and equilibrium parameters describing the slow relaxation process can be obtained from long time sorption data. To determine if the initial portion of the sorption profile , which exhibited "two stage" sorption behavior , followed Fickian behavior , and to estimate the diffusion coefficient the following treatment was applied with the aid of a computer - assisted fitting routine. Initially , values of M.» were intuitively selected by inspection of a plot 0f the experimental Mt vs. 1 data. This is illustrated for sorption studies discribed in Figure 26, as shown, the Mac value for a Fickian diffusion process is assigned , the t 1/ 2 value can be obtained from a plot of Mt / Moo versus t“2 ( Figure 27 ) and the diffusion coefficient ( Ds ) calculated by solution of equation (22). By substituting the Ds value into Equation (21) , calculated values of Mt are obtained. The calculated and experimental values of Mt / Moo versus t“2 are then plotted and the quality of the fit is evaluated by graphical analysis. If the experimental and calculated values show poor agreement , this 55 procedure is repeated with a new value of Moo until a satisfactory f it between the calculated and experimental values is obtained ( Figure 26 ). As shown in Figure 28, an assigned value of Moo = 0.275 , gave good agreement ( Figure 28 ) between the experimental and calculated data for the sorption profile plot at a = 0.295. The estimated diffusion coefficient for the initial stage of the sorption profile ( i.e, Fickian ) was 1.264 13—10 cm2/sec. Wt. Gain [m glm 9] Weight Gainlm g] 56 0.07 _ 0.06 0.05 002: —-o-- WLGng “n.6mn1 I k I l l 1 4 . I 00 50 100 150 200 250 Time(hours) Figure 7. Evaluation of the reproducability of the sorption system at a = O - 0.06 , 001 l 000 000 000 000 000 ‘ ' ‘ 4‘ ‘ ' 0 10 20 30 fimaHm) Figure 8. Interval sorption of n- propanol at a = O - 0.035 Wt. Gain[m 9]] 0.03 0.02 ’ 0.01 ' 0.00 . 1 1 1 1 1 . . 1 :13 35 37 39 Times(Hrs) Figure 9. Interval sorption of n- propanol at a =0035 - 0.050 Weight Gain [m g] 0.10 0.05 0.00 0 2 4 6 8 10 12 Times(Hrs) Figure 10. Interval sorption of n- propanol at a = 0.050- 0.080 Weight Gain [mg] wt. gainImg] 0.0 A L A I 1 L l l 0.0 5.0 10.0 15.0 20.0 Times(Hrs) Figurei 1. Interval sorption of n- propanol at a=0.080-0.110 80 Times(Hrs) Figure12. Interval sorption of n- propanol at a=0.11-0.295 Wt. Gejn[m 9] WI. Gain[mg] 59 Times(Hrs) Figure13. Interval sorption of n- propanol at a=0.295-O.390 I l 1 l l L l 4 J '0 so 100 Times(Hrs) Figure14. Interval sorption of n- propanol at a=0.39-0.490 Wt. GainIm 9] WI. Gain[m g] 60 0.0 1 l . l . 50 100 150 Times(Hrs) Figure15. Interval sorption of n- propanol at a=0.49-0.840 l 0.0 A a 4 1 a a 1 a n A 50 1 00 Times(Hrs) Figure16. Interval sorption of n- propanol at a=0.84-0.910 61 oooooooooooooo 1.5 - ,.-- 9 .8 ......... I ......... :I' ’3 1 O _ "7 ........... " E ' f..- g 6. =3 5 1: a = 0.035 6: a = 0.390 E o s _ / """"" 2:a=0.050 7:a=0.490 ' 3: a=0.080 8: a=0.840 3‘3"; 4:a=0.110 923:0.910 1 2 .---: 5: a = 0.295 0.0 ("""" 1 I 1 1 I 1 1 I L 1 I 1 1 I 1 1 I L 4 0 100 200 300 400 500 600 700 Time(hrs.) Figure 17. Successive interval sorption of n- propanol into Nylon 61/6T at 23 0C (a = 0 - 0.91) P 0.20 - g a=0.05 - 0.08 I. ‘E’ '8 Q g 0.10 _ a=0.035 - 0.05 a = 0 - 0.035 0.00 1 1 I 1 1 I 0 50 100 Time(hrs.) Figure 18. Successive interval sorption of n- propanol into Nylon 6I/6T at2300(a=0-0.110) Volume Fraction Weight Fraction[glgx100) 0.30 0.25 0.20 0.05 0.00 62 1 I 1 J I 4 1 I 1 1 l #4 I J 4 0.02 0.04 0.06 0.08 0.10 0.12 Vapor Activity Figure 19. Experimental sorption isotherm for propanol in 0.011 0.009 0.007 0.005 0.003 0.001 Nylon 61/6T at low vapor activities (a < 0.11) I I ' Y T T V I 4 j A L M J l A j l l L I I 1 0.06 0.08 0.10 0.12 0.14 Vapor activity I 0.04 Figure 20. Flory-Huggins plot when I =1.6,from the equation Ina=InV1+V2+V2 63 8 2 ‘2" 1 I 1 I 1 l 1 l 1 0 40 80 120 160 200 t1/2,Sec1/2 Figure 21. Plot of Mthoo vs t1/2 for n- propanol sorption by Nylon 6|I6T, at 23 00 (a: 0.0 - 0.035) 110 _ 100 [- 90 fi- 80 :- 70 :- 6 60 E 5 so 3 5 40 Z- 30 Z- 20 E- 10 r 0 g I 1 L 1 l 1 0 4o 80 120 160 11/2,sec1/2 Figure 22. Plot of Mt/Moo vs t1/2 for n- propanol sorption by Nylon 6|/6T,at 23 00 (a = 0.035 - 0.050) MIIMoo Mtllvloo 110 100 64 I L I m I 1 I 40 80 120 160 200 t1l2,sec1l2 Figure 23. Plot of Mt/Moo vs t1l2 for n- propanol ’sorption 100' 90 80 70 60 50 40 30 20 10 by Nylon 6l/6T,at 23 0C (a = 0.05 - 0.08) I 1 J 40 80 120 160 200 t1l2,sec1/2 Figure 24. Plot of Mthoo vs t1l2 for n- propanol sorption by Nylon 6|/6T,at 23 00 (a = 0.08 - 0.11) 65 110 100' ‘1 J 90 .. 80':- 9 7O " 60 50 40 30 A (Mt/Moo)experi x100 28 r - (Mt/Moo)ca|cu. x100 1 O . I 1 I 1 I 0 100 200 300 11/25601/2 'U' [MtIMoolxIUU Y'U Figure 25. Comparison between experimental and calculated data at vapor activity a = 0.035 0.45 1 0.40 - 0.35 . 0.30 0.25 ._. 0.20 CD 3; 0.15 -§ '1 :Moo 0.275 mg 0 0'10 2 :Moo 0.284mg is; 0.05 3 :Moo 0.298 mg g 0.0011I11111I14L11I1PIL41.J 0 10 20 30 4O 50 60 70 80 t, Hrs Figure 26. Deviation from Eqn.(21) due to non-Fickian sorption for sorption of propanol by Nylon 6l/6T at 230C(a=0.11- 0.295):Graphical Estimation of Fickian Moo value 66 Ivlt I Moo [gig] OJ 0.0 1 I 1 I 1 I 1 I 1 I 1 I4 [*1 1 I 1 O 20 4O 60 100 120 140 160 180 200 11/2 = 84 sec. t1l2.sec1/2 Figure 27. Plot of Mt/Moo vs. t1l2 for n- propanol sorption by Nylon 6I/6T ,23 00, a = 0.11 - 0.295,Eqn.(21) for Moo = 0.275 and 11/2 = 84 sec. MI I Moo [gig] 0.0 1 41 J 100 200 tI/2,sec1l2 Figure 28. Plot of Mt/Moo vs. t 1/2 for n- propanol sorption by Nylon 6II6T, 23 0C (a = 0.11 - 0.295 ): point observed, curve, Eqn.(21) for Moo = 0.275 and t 1/2 = 84 sec. Conclusion The sorption of propanol by amorphous polyamide, evaluated by an incremental sorption technique, showed that simple Fickian kinetics were obeyed and equilibrium was achieved within a short time at very low vapor activities ( a1 < 0.11 ) . However , at high activity levels ( a1 > 0.11 ) , the initial,rapid sorption was followed by a slower further uptake of propanol,which was qualitatively attributed to a relaxation-controlled swelling process. The initial rapid stage of sorption may be related to the Fickian diffusion of penetrant into pre-existing and available vacancies or sites in the glassy polymer. The relaxation processes , which occur over a long time scale may be related to a structural reordering or redistribution of free volume elements to provide additional sites of suitable size and accessability to accommodate more penetrant molecules. 67 Appendix A Procedure of standard calibration curve construction In all cases , a standard curve of response vs. penetrant concentration was constructed from standard solutions of known concentration. Calibration solutions were prepared by dissolution of known quantities of n- propanol in dichlorobenzene. The detailed procedure follows: 1. Bake out vials and syringes in oven prior to use to remove any residual solvent or permeant. Cool to room temperature. 2. Using the gas chromatograph , check the retention time of each solvent by head space technique to ensure there is no interferring peaks at the propanol retention times. 3. Prepare dilute solutions for the permeant standards by the following procedure: - fill up 100 ml volumetric flask with 10 ul of n-propanol and dilute to make with dichlorobenzene. This becomes 100 ppm standard solution. - make up similar standard solution with dichlorobenzene , which becomes a standard solution of 50 ppm concentration. - Make 2 or 3 more different concentration solutions by using the same method 4. From the lowest concentration solution , inject 0.5 111 sample solution directly into the gas chromatograph and the area response recorded. 68 69 Replicate the run and calculate the average Plot the gas chromatograph area unit response vs. the number of grams injected per sample. The slope of this curve equals the calibration factor. Injected quantity can be calculated as follows W .—: concentration(v/v) x 0.5 111 x 1 ml/ 1000 111 x propanol density(g/ml) Setting conditions of gas chromatograph are as follows : o Oven temp. : 50 0 C 0 Initial time : 5 min. 0 Rate : 3 0 / min. 0 Final time : 0 min. 0 Final temp. : 165 0 C o Inject temp.: 200 0 C o Detect temp.: 250 0 C o Oven max. : 250 0 C 0 He carrier gas : 33 ml/min. 0 Sample injected : 0.5 ul 0 Head space test : 100 1.11 inject 0 Concentration used : 100 , 50 , 25 , 10 , 5 ppm ( 5 kinds ) Appendix B Electrobalance Calibration Procedure Cahn D-200 Electrobalance can be calibrated by the following procedure. . Suspend the sample container on the hangdown wire . Press the down arrow key three times and the NET WEIGHT box will display the weight that the balance is actually seeing. . Enter "Y" into the REZERO box. Next, press the up arrow key once to go back to the Balance Calibrate menu. . Enter into Tare Balance box "Y" and the balance will tare. . Press the down arrow key to input the calibration weight . Remove the hangdown tube so that the pan is exposed. Then put the 100 mg calibration weight onto the sample weight pan. Replace the hangdown tube for stable reading. . Place the cursur into the Calibrate Balance box and enter "Y" and press ENTER when the pan with the calibration weight is steady. . Lower the hangdown tube and remove the calibration weight. . Suspend the sample and set up run conditions. 70 Yap or Activity trap or Actin’ty 7 1 Appendix C. Vapor activity fluctuation 0.10 0.08 - 0.06 - 0.04 "' Bun B D— 002 .. Mean: 0.035 Standard deviation: 0.003 000 I 1 I 1 l - 1 0 10 20 30 40 Time(Hrs) Figure 29. Distribution of vapor activityat a = 0.035 0.10 0.08 0.06 0.04 0.02 0.00 . .._B___n—-u—F T Mean:0.050 " Standard deviation:0.0016 1 l 1 I 1 J 1 I L I 1 I 1 1 2 3 4 5 6 7 8 ' Times(Hrs) Figure 30. Distribution of vapor activity at a = 0.050 i.iap or Activity trap or Activity \1 N 0.12 _ 0.10 :- E a E, 0.08 _ 43—43— r 0.06 _- 0.04 :- I Mean: 0.08 0,02 :. Standard devration: 0.0042 0.00 : J— 1 J l L l 1 L 4 6 8 10 12 14 Times(Hrs) Figure 31. Distribution of vapor activity at a = 0.080 0.30 0.25 1 0.20 - 0.15 1 0.10 fir—”J; F + O 05 E Mean: 0.110 ' - Standard deviation: 0.003 O’COI- l 1 l 1 l 1 I 1 O 10 15 20 25 3O Times(Hrs) Figure 32. Distribution of vapor activity at a = 0.110 Yap or Activity Yap or Activity 0.50 0.40 - 0.35 0.45 - h 030E. a 3334 W 5' Mean: 0.295 0‘25 0 Standard deviation:0.0069 0.2011I11I11I11J11l11I1LI11 0 10 20 30 40 50 60 70 80 Times(Hrs) Figure 33. Distribution of vapor activity at a = 0.295 0.50 0.45 i- 0.40 - fl 4L .35 0.35 i' Mean: 0.390 Standard deviation: 0.0084 0.30 LIAAIAIIJIIILIAALAIII; 0 10 20 30 40 so 60 70 80 Times(Hrs) Figure 34. Distribution of vapor activity at a = 0.390 Yap or Activity Yap or Activity 0.60 0.58 0.56 0.54 0.52 0.50 in a 0.48 0.46 0.44 0.42 0'40 4 1 1 1 1 L 1 1 1 1 1 I 1 1 1 1 0 50 100 150 Times(Hrs) U I I I V ljl l E) ' I mean: 0.490 Standard deviation: 0.015 I V ‘l I J Figure 35. Distribution of vapor activity at a = 0.490 1.00 0.95 '- O.90 '- a L In 0.85 . a a . a a a 0.80 . : Mean: 0.84 0-75 _” Standard deviation: 0.027 0.70 h 1 l 1 l 1 1 1 l 1 L 1 I 1 l l_ J 1 l 1 0 20 40 60 80 100 120 140 160 180 200 Times(Hrs) Figure 36. Distribution of vapor activity at a = 0.840 Yap or Activity 1.00 0.95 0.85 0.80 0.90 Mean: 0.91 Standard deviation: 0.0046 l 1 1 l 1 1 l 1 1 l l 1 1 I 1 1 1 1 20 40 ' 60 80 100 i 20 140 Times(Hrs) Figure 37. Distribution of vapor activity at a = 0.910 BIBLIOGRAPHY Bagley , F and F. A. Long, 1955. J. Am. Chem. Soc., 77 , 2172. Barrer, R., J .Barrie and J. Slater, 1957. J. Polymer Sci.,23, 315 Beck,J . V and Arnold. K. J , 1977. "Parameter estimation in engineering science" , John Wiley and Son , New York , NY. Berens, A. R. 1975. Angew. Macromol. Chem., 47, 97. Blatz, P. S. 1989. "Innovations in polyamide technology" AICHE Meeting, Houston, TX, April 2 - 6. Blumstein, A., 1978. "Liquid Crystalline Order in Polymers", Academic Press, New York. Choy, C. L., W. P. Leung and T. L. Ma, 1984. J. Polym. Sci: Polym. Physics Edition, Vol. 22, 707-719 Cohen, M. H. and D. Tumbull , 1959. J. Chem. Physics , 31 , 1164. Dolden, J G., 1976. "Structure property relationships in amorphous polyamides" Polym. Vol. 17 , pp 875 Gedraitite, G. B., A. P. Martin and V. A. Shlyapnikov, 1989. Eur. Polym. J ., 25:39. Gennes,P. G., 1979. "Scaling Concepts in Polymer Physics", Cornell Univ. Press, Ithaca, NY. Hartley, G. S. , 1946. Trans. Faraday Soc., 423, 6. Hernandez, R. J ., J. R. Giacin and A. L.Baner, 1986. J. Plastic Film and Sheeting, 2 (3) : 187. Hernandez, R. J ., J. R. Giacin and E. A. Glulke, 1990. J. Membrane Science 76 77 Hopfenberg,H.B. and V.Stannett, 1973. "The physics of glassy polymers", R. N. Haward(ed), chapter 9, The diffusion and sorption of gases and vapors in glassy polymers, Applied science publishers Ltd, London, pp 504—506. Kwei, T. K. and T. T. Wang , 1972. Macromolecules , 5 , 128. Meares,P. , 1954. J. Am. Chem. Soc., 76, 3415 - 22. Meares,P., 1957. Trans. Faraday Soc., 53 , 101 - 6. Meares,P., 1958. Trans. Faraday Soc., 54 , 40 - 6. l Michaels. A. S.,W.Vieth and J.Barrie, 1963. J. Appl. Phys.,34(1),1,l3 : Michaels. A. S. and H. J. Bixler, 1961. J. Polym. Sci., 50, 393 Ohashi , K., R. J. Hernandez , J . R. Giacin and E. A. Gruke , 1991. " Modeling the sorption of water vapor by a semicrystalline polyamide ", Thesis , school of packaging , Michigan state univ. Pace, R. J. and A. Datyner, 1980. "Model of sorption of simple molecules in polymer", J. Polymer Sci., Vol. 18, 1103-1124 Pechhold, W. and H. P. Grossmann, 1979. Disc. Faraday Soc., 68, 58. Peterlin , A. 1975. J. macromol. sci. physics , B 11 , 57 Rogers, CE, 1965. " physics and Chemistry of the Organic solid state" , (Eds Fox , Labes , ans Weissbergcr) , Interscience , NewYork , Ch 6 Salame, M., 1986. J. Plastic Film and Sheeting, 2 (4) : 321 Vieth, W. R., H. H.Alcalay and A. J. Frabe, 1961. J. Appl. Polymer Sci. 8, 2125. Wendorff, J. H., 1987. "Studies on the nature of order in the amorphous polymers", International Symposium. Williams, C., F. Brochard and H. L. Frisch, 1981. Ann. Rev. Phys. Chem. 32, 433-451. "‘11111111111111“