# Role of Polymer Interactions with Clays and Modifiers on Nanomechanical Properties and Crystallinity in Polymer Clay Nanocomposites

Authors: Dinesh R. Katti; Kalpana S. Katti; M. Raviprasad; Chunju Gu
Journal: Journal of Nanomaterials
 (2012)

---

## Abstract

We present the effect of organic modifier on crystallinity and nanomechanical properties of polymer clay nanocomposites (PCNs) using two different polymers while maintaining the same nanoclay and organic modifier. Experimental results and interaction energy maps of Polybutylene-Terephthalate- (PBT-) PCN system indicate that the underlying mechanisms of change in crystallinity and improvement in mechanical properties as proposed in altered phase theory are valid. Experimental and molecular simulation studies of PBT-PCN and Nylon6-PCN reveal that a higher crystallinity polymer could require significantly higher attractive and repulsive interaction energies between polymer and organic modifiers to change the crystallinity of the polymer in the PCN significantly and thus improve mechanical properties of the PCN.


---

## Body

## 1. Introduction

Polymer clay nanocomposites (PCNs) are widely researched material systems for wide variety of applications in industry because of enhanced mechanical, thermal, optical, and barrier properties. PCNs are categorized as intercalated or exfoliated and are typically synthesized by adding small amount (about 1–9% of weight polymer) [1, 2] of organically modified nanosized expansive clay. Addition of organically modified clay enhances the mechanical properties [3], thermal properties [4–6], and barrier properties [4, 6–9] of polymer clay nanocomposites. Most importantly, the amount of clay in PCN influences the crystallinity of polymer and the mechanical properties [10]. Studies show that the properties of PCNs depend on individual phases and interaction between individual constituents [10–12]. Our previous work on PCNs has shown that the nonbonded interactions between different constituents alter the polymer phase [10–12] and thus significantly influence the mechanical properties of PCNs [1]. The organic modifier, which is used to modify the clay, has a significant influence on crystallinity of polymer and elastic modulus of PCNs [11]. Effect of different organic modifiers, which contain different functional groups and different chain lengths, is studied with Nylon6 in detail [11] and the results indicate that the functional group and backbone chain length of organic modifier have a significant influence on d-spacing of PCNs, crystallinity of polymer, and the elastic modulus of PCNs. In this paper, we present the results of experimental and quantitative modeling studies on PCN with two different polymers, Nylon6 and Polybutylene Terephthalate (PBT), and same organic modifier, and nano clay.PBT is a semicrystalline thermoplastic polymer which shows excellent properties such as rate of crystallization, thermal stability and processability [13]. PBT has been used widely in nanocomposites [14–17] and shows improved tensile modulus and strength [17], enhanced thermo mechanical properties and elastic modulus [6]. High viscosity PBT nanocomposites show improved mechanical, dynamic mechanical properties [16]. Several experimental studies conducted on PBT clay nanocomposites can be found in the literature [6, 13–23]. In this paper, PCN containing Nylon6 is referred to as Nylon6-PCN and the PCN containing PBT is referred to as PBT-PCN. Properties of Nylon6-PCN are studied at the molecular scale using molecular dynamics and experiments [10–12, 24–28]. Experimental results of Nylon6-PCN used in the discussion is from our previous work [10]. Molecular dynamics simulations provide an insight into mechanisms responsible for the enhancement of PCN properties. Molecular models are designed to be close to experimental condition and models are validated with experimental d-spacing. Quantitative interaction energies between PBT-PCN constituents are calculated and the results are compared quantitatively with Nylon6-PCN. Crystallinity and elastic modulus tests are conducted using differential scanning calorimetry (DSC) and nanoindentation technique, respectively. This paper describes a comprehensive quantitative study of interaction energies and the effect of interactions on properties of polymer clay nanocomposites.


## 2. Materials and Methods

### 2.1. Materials

Na-montmorillonite with cationic exchange capacity of 76.4 mequiv/100 g was obtained from the Clay Minerals Repository at the University of Missouri, Columbia, MO, USA. 12-Aminolauric acid [NH2(CH2)11COOH] was supplied by the TCI America. Poly(1, 4-butylene terephthalate) (Mv=38,000) was obtained from the Sigma-Aldrich, Inc., USA. Polyamide 6 (PA6) [(CH2)5CONH]n (Mv=16000) was supplied by the Polysciences, Inc., USA. Trifluoroacetic acid (99.9%) and chloroform (99.97%) were purchased from the EMD chemicals Inc., Germany. Methanol anhydrous (99.9%), formic acid (97%), silver nitrate, standard volumetric solution, and hydrochloric acid (HCl) were purchased from the Alfa Aesar, MA, USA.


### 2.2. Preparation of Organically Modified Montmorillonite (OMMT)

Ten grams of Na-montmorillonite (Na-MMT) are crushed into fine powder and screened through #325 sieve (45μm) and then dispersed into deionized (DI) preheated (333 K) water (400 mL for PBT, 200 mL for PA6). The clay solution is stirred vigorously for 1 hour to disperse the clay uniformly and swell sufficiently. 3.312 g (15.4 mmol) of 12-aminolauric acid and 2.5 mL of 0.1 N HCl are dissolved in DI water at 333 K and stirred vigorously to form a uniform solution. The clay suspension is gradually added into the acid solution and is stirred vigorously at 333 K for 1 hour to obtain the organically modified montmorillonite (OMMT). Prepared OMMT is washed several times with DI water to completely remove Cl- and then the solution is tested by 0.1 N AgNO3 to verify the absence of Cl- ions. Finally, the OMMT cake is dried in an oven at 333 K for 24 hours and then ground and screened through #325 sieve for synthesis of PCN.


### 2.3. Preparation of PBT Film and PBT-Clay Nanocomposite (PBT-PCN) Film

Five grams of PBT are added to 42 ml chloroform and 7 mL trifluoroacetic acid (TFA) and stirred for 2 hours to make a uniform solution. Subsequently, 0.45 g (9 wt% of OMMT with respect to pure PBT) OMMT is added to the PBT-CHCl3-TFA solution and then stirred vigorously for 3 hours to obtain PBT-PCN. PBT and PBT-PCN films are obtained by casting the corresponding solutions on a glass slide using a drawdown bar (Paul N. Gardner Company, Inc.) set at 25 mils thickness.


### 2.4. Preparation of Nylon6-Clay Nanocomposite (Nylon6-PCN)

Ten grams of Nylon6 are added to 155.73 ml formic acid (FA) and stirred vigorously for 40 minutes to prepare a uniform solution. Next, 0.9 g (9 wt% of OMMT with respect to pure nylon6) of OMMT is dispersed into the nylon6-FA solution and stirred vigorously for 6 hours at room temperature. DI water is added to the resulting solution to form PCN and to remove the FA from PCN. Finally, the remaining DI water is removed by adding methanol anhydrous and the resulting solution is placed in a vacuum oven set at 313 K temperature for a period of 24 hours to obtain the PCN.


### 2.5. Characterization

#### 2.5.1. XRD Characterization

The X-ray analysis is performed by X-ray diffractometer (Philips Analatical X’pert MPD, Almelo, Netherlands) with secondary monochromator and Cu-tube by CuKα radiation wavelength of 1.5406 Å. X-ray diffraction of MMT and OMMT is conducted at a scan rate of 0.01°/s in each case and with scan range of 2θ = 2.07–35.02° for MMT and 2θ = 2.11–35.06° for OMMT. The same OMMT was used to prepare PBT-PCN and nylon6-PCN.  For the PBT-PCN film a scan rate of 0.02°/s and the scan range of 2θ = 2.03–65.01° are used. The Nylon6-PCN is compressed into thin sheets by compression testing machine (SATEC, model 22 EMF, PA, USA) and the X-ray diffraction of Nylon6-PCN was done at a scan rate of 4°/min with the scan range of 2θ = 2.01–15°.


#### 2.5.2. Photoacoustic- (PA-)FTIR Characterization

PA-FTIR experiments are conducted in FTIR bench (model Nexus, Thermo Nicolet, USA) with a photoacoustic cell (MTEC Photoacoustics model no. 300). Sample is prepared by pressing the OMMT powder into a pellet. The nylon6 sample is heated gradually from 303 K to 493 K and then pressed into a pellet. The PBT and PBT-PCN samples are prepared in smaller size and placed in the sample holder. The PA-FTIR is conducted in the linear photoacoustic mode for 500 scans at a resolution of4cm-1, with moving mirror velocity of 0.1581 cm/s. FTIR data are acquired from 400 to 4000 cm-1.


#### 2.5.3. Differential Scanning Calorimetry (DSC)

DSC experiments are conducted using differential scanning calorimeter (model Q1000 DSC) of TA instruments, New Castle, Delaware, USA. PBT and PBT-PCN samples are heated from 25°C to 250°C with nitrogen atmosphere with the heating rate of 10°C/min for heating/cooling/heating cycle with retention time of 3 minutes at the end of each treatment for temperature equilibration. Nylon6 and Nylon6-PCN samples are heated from 30°C to 100°C in 3 minutes and the temperature is maintained at 100°C for 10 minutes. A heating rate of 10°C/min is used for the first heating and cooling cycles with retention time of 10 minutes at 250°C and the heating rate of 20°C/min [11] is used in the final heating cycle. We have taken the net weight of the respective polymers into consideration. The enthalpy of fusion was determined by integrating the area under the DSC endotherm peak which was divided by the net weight of the material. The percentage of crystallinity Xc is determined from DSC measurements by
(1)Xc=ΔHmΔHm°100%,
where ΔH is melting enthalpy of sample [J/g]; ΔHm° represents theoretical value of enthalpy for 100% crystalline polymer. For PBT homopolymer, ΔHm°=141 J/g [29]; for nylon6, ΔHm°=203 J/g [11]. This is generally considered a standard method to measure crystallinity of polymer which has been described earlier [30] and extensively used in literature [10, 29, 31, 32].


#### 2.5.4. Nanoindentation

Nanoindentation experiments are conducted using triboscope nanoindenter (Hysitron Inc., Minneapolis, MN). The triboscope is operated in association with a multimode AFM controlled by the nanoscope IIIa controller and with the J-type piezo scanner from Veeco (Veeco Metrology Group, Santa Barbara, CA). For PBT and PBT-PCN, load controlled nanoindentation tests are performed using a 100–200 nm Berkovich (three-sided pyramid) diamond indenter tip. Experiments are conducted for maximum load of 2900μN. The reduced elastic modulus (Er) and hardness (H) values are calculated from initial part of the unloading curve using Oliver-Pharr method [33]. The average Er and H are obtained from 30 numbers of indentations for each specimen. For nylon6 and Nylon6-PCN, load controlled nanoindentation tests are performed [11]. The elastic modulus is calculated using the following equation:(2)1Er=1-v2E+1-vi2Ei,
where E and ν are the elastic modulus and Poisson’s ratio of the sample respectively. Ei and νi are the elastic modulus and Poisson’s ratio for the diamond indenter. In the case of polymeric materials and its composites, due to the fact that Ei≫E, (2) can be simply given as [34]                                                       (3)E=Er(1-ν2).


### 2.6. Simulations Using Molecular Dynamics

Molecular dynamics (MD) is a computational method used to study the position and interaction of atoms and molecules. MD is used extensively to study proteins, bio molecules, polymers, and minerals. MD provides a good understanding of molecules with respect to their interactions and molecular structure. This method involves solving the equation of motion in small time steps. The equation of motion isFi=miai, where Fi is the force acting on an atom i, mi is the mass of an atom i, and ai is the acceleration of an atom i. Change in energy of an atom i with respect to its position (r) is equal to the force acting on an atom i. By solving these equations of motion, MD provides time-dependent response for the molecular system. This response contains position of every atom in the system and their velocities. In this work, MD has been used to study the interactions between atoms and group of atoms in PBT-PCN and Nylon6-PCN.


#### 2.6.1. Force Field Parameters for PBT, Nylon6, and Na-MMT

Force field describes the potential energy of a molecular system. Potential energy of any molecular system consists of bonded energy and nonbonded energy. Bonded energy consists of bond, angle, dihedral, and improper energy of the molecular system. Nonbonded energy consists of van der Waals and electro static energy of the molecular system. Total potential energy of the molecular system can be described by the following expressionEtotal=[Ebond+Eangle+Edihedral+Eimproper]+[Evan+Eelec]. These individual potential energy terms can be described by mathematical functions and the functions depend on force constants and variables. Some of the general force fields used by researchers are CHARMM, CFF, AMBER, and COMPASS. In our previous work, CHARMM [41] (Chemistry at HARvard Macromolecular Mechanics) force field parameters were used for Polyamide6 and 12-Aminolauric acid and the same set of parameters were used in the present work. We have transformed CFF force field parameters of sodium Montmorillonite to CHARMM force field parameters [42–44]. Force field parameters for PBT were obtained by considering number of methyl group and model components which are similar to polyethylene terephthalate [45].


#### 2.6.2. Molecular Model Construction and Simulations

Materials Studio software, developed by Accelrys Software Inc., is used to build the molecular models. This software is widely used to construct molecular models. NAMD2.5b2 [46] and VMD1.8.7 [47] are used to conduct simulations of OMMT, PBT, Nylon6, and PCNs and visualize the results from simulation, respectively. These programs are developed by Theoretical Biophysics Group at Illinois Beckman Institute and widely used by researchers to study complex molecular system.Organically Modified Montmorillonite (OMMT) Model
OMMT model consists of two sheets of montmorillonite (MMT) and eighteen numbers of 12-aminolauric acid molecules. Size of the OMMT model is 31.68 Å inX direction, 27.44 Å in Y direction, and 24.16 Å in Z direction. Each montmorillonite sheet consists of eighteen unit cells where six unit cells are in X direction and three unit cells are in Y direction. Each of these MMT sheets has nine negative net charges and sheets are placed at 17.60 Å initial d-spacing. In our previous work [26], we have studied different initial d-spacing and we found that 17.60 Å initial d-spacing model provides closest experimental d-spacing at the end of simulation. In order to make OMMT model charge neutral, eighteen numbers of 12-aminolauric acid chains are inserted in the model where nine are inserted in between two MMT sheets and nine are inserted above the MMT sheet as shown in Figure 1. Previous work [25] showed that these 12-aminolauric acid molecules should be placed flat and parallel to the MMT sheet. Partial charges and structure of 12-aminolauric acid are shown in Figure 2(a).Figure 1
Initial OMMT model.Distribution of partial charges (a) 12-aminolauric acid, (b) PBT, (c) Nylon6.
(a)
(b)
(c)Polybutylene Terephthalate (PBT) Model
PBT monomers are used to build polymers. A representative molecular model of PBT is shown in Figure3. Six different lengths of PBT polymers were built to study the length effect of polymers on d-spacing of clay in PCN. These polymer chains contain two, four, six, eight, ten, and twelve monomers. All PBT polymers were annealed at 700 K under atmospheric pressure. Annealing process was conducted three times to obtain globally minimized confirmation of PBT chains and minimization was confirmed from the potential energy.Figure 3
Single chain molecular model of PBT.


#### 2.6.3. Simulation of OMMT

Initial OMMT model is minimized for 25 ps at a temperature of zero Kelvin in vacuum using molecular dynamics program NAMD. Next, the temperature of the OMMT system is raised to 300 K in three equal steps by Langevin temperature parameter. Once the temperature is increased to 300 K, pressure of the OMMT system is increased to 1.01325 bars in four equal steps by Langevin piston pressure control. Next, the system temperature is increased by 33 K and subsequently reduced to atmospheric temperature. The increase and reduction of temperature is done to mimic the experimental conditions as mentioned in Section2.2. A 0.5 fs time step is used for all OMMT simulations. Finally, OMMT model is simulated for 100 ps at room temperature and pressure with periodic boundary condition. During OMMT simulation, constraints were applied to MMT atoms in X and Y direction but MMT atoms are allowed to move in Z direction. There are no constraints applied to organic modifiers thus organic modifier atoms are allowed to move in all three directions. Simulation of OMMT is done under constant number of particles, constant pressure, and constant temperature (NPT). 16 Å cut-off distance and 14 Å switch distance are used for van der waals and electrostatic energy calculation in all simulations.


#### 2.6.4. Polymer Clay Nanocomposites (PCN) Models

Selected OMMT model, as shown in Figure4, and minimized PBT models are used to construct PCN models. Six different lengths of PBT chains are used in PCN to select the suitable PCN system. One of the initial PCN models is shown in Figure 5. The suitable PCN system is selected based on experimental d-spacing as shown in Figure 6. The estimated mass of PBT chains in PCN model with respect to clay is shown in Table 1. Two PBT chains are inserted into the OMMT model and one is placed in between clay sheets and the other is placed above top clay sheet to ensure the periodic boundary condition during the simulation. Physical dimensions of PBT are measured and compared with the simulated interlayer spacing of OMMT to make sure that the PBT chain fits in the model. The same conditions described above were implemented to construct the nylon6 PCN model. The number of Nylon6 chains in PCN models is selected based on our previous work [11]. Initial and final Nylon6-PCN model are shown in Figures 7, and 8 respectively. Final model is selected based on the experimental d-spacing.Table 1
Estimated mass percentage of PBT with respect to clay.
Model
Intercalated polymer wt% with clay
PCN(1)
1.96
PCN(2)
4.19
PCN(3)
6.41
PCN(4)
8.64
PCN(5)
10.87
PCN(6)
13.10Figure 4
Final OMMT model.Figure 5
Initial PBT-PCN model.Figure 6
Final PBT-PCN model.Figure 7
Initial Nylon6-PCN model.Figure 8
Final Nylon6-PCN model.


#### 2.6.5. Simulation of PCNs

Molecular simulation of PCN is conducted for PBT-PCN and Nylon6-PCN system. Procedure for minimization and dynamics of PCNs is done as described in Section2.6.3 Cell basis vector for PCNs is adjusted according to the physical size of different models. Periodic boundary conditions are applied to calculate the electrostatic interactions between PCN constituents. Initially, PCN models are simulated for 25 ps. Based on d-spacing, final PCN model is selected and then the final simulation is conducted for 200 ps. Constraints are applied to MMT as mentioned in Section 2.6.3 and polymers are allowed to move freely in all three directions. Constant number, pressure, and temperature (NPT) are used for all PCN simulation. Constant pressure is maintained by Nose’-Hoover Langevin piston pressure control method. This method is used with Langevin dynamics, such as temperature control, in order to perform NPT ensemble. Equilibrium of PCN system is confirmed from the plot of total energy of the system versus simulation time and it is found that 200 ps simulation time is sufficient for the PCN simulations.


## 2.1. Materials

Na-montmorillonite with cationic exchange capacity of 76.4 mequiv/100 g was obtained from the Clay Minerals Repository at the University of Missouri, Columbia, MO, USA. 12-Aminolauric acid [NH2(CH2)11COOH] was supplied by the TCI America. Poly(1, 4-butylene terephthalate) (Mv=38,000) was obtained from the Sigma-Aldrich, Inc., USA. Polyamide 6 (PA6) [(CH2)5CONH]n (Mv=16000) was supplied by the Polysciences, Inc., USA. Trifluoroacetic acid (99.9%) and chloroform (99.97%) were purchased from the EMD chemicals Inc., Germany. Methanol anhydrous (99.9%), formic acid (97%), silver nitrate, standard volumetric solution, and hydrochloric acid (HCl) were purchased from the Alfa Aesar, MA, USA.


## 2.2. Preparation of Organically Modified Montmorillonite (OMMT)

Ten grams of Na-montmorillonite (Na-MMT) are crushed into fine powder and screened through #325 sieve (45μm) and then dispersed into deionized (DI) preheated (333 K) water (400 mL for PBT, 200 mL for PA6). The clay solution is stirred vigorously for 1 hour to disperse the clay uniformly and swell sufficiently. 3.312 g (15.4 mmol) of 12-aminolauric acid and 2.5 mL of 0.1 N HCl are dissolved in DI water at 333 K and stirred vigorously to form a uniform solution. The clay suspension is gradually added into the acid solution and is stirred vigorously at 333 K for 1 hour to obtain the organically modified montmorillonite (OMMT). Prepared OMMT is washed several times with DI water to completely remove Cl- and then the solution is tested by 0.1 N AgNO3 to verify the absence of Cl- ions. Finally, the OMMT cake is dried in an oven at 333 K for 24 hours and then ground and screened through #325 sieve for synthesis of PCN.


## 2.3. Preparation of PBT Film and PBT-Clay Nanocomposite (PBT-PCN) Film

Five grams of PBT are added to 42 ml chloroform and 7 mL trifluoroacetic acid (TFA) and stirred for 2 hours to make a uniform solution. Subsequently, 0.45 g (9 wt% of OMMT with respect to pure PBT) OMMT is added to the PBT-CHCl3-TFA solution and then stirred vigorously for 3 hours to obtain PBT-PCN. PBT and PBT-PCN films are obtained by casting the corresponding solutions on a glass slide using a drawdown bar (Paul N. Gardner Company, Inc.) set at 25 mils thickness.


## 2.4. Preparation of Nylon6-Clay Nanocomposite (Nylon6-PCN)

Ten grams of Nylon6 are added to 155.73 ml formic acid (FA) and stirred vigorously for 40 minutes to prepare a uniform solution. Next, 0.9 g (9 wt% of OMMT with respect to pure nylon6) of OMMT is dispersed into the nylon6-FA solution and stirred vigorously for 6 hours at room temperature. DI water is added to the resulting solution to form PCN and to remove the FA from PCN. Finally, the remaining DI water is removed by adding methanol anhydrous and the resulting solution is placed in a vacuum oven set at 313 K temperature for a period of 24 hours to obtain the PCN.


## 2.5. Characterization

### 2.5.1. XRD Characterization

The X-ray analysis is performed by X-ray diffractometer (Philips Analatical X’pert MPD, Almelo, Netherlands) with secondary monochromator and Cu-tube by CuKα radiation wavelength of 1.5406 Å. X-ray diffraction of MMT and OMMT is conducted at a scan rate of 0.01°/s in each case and with scan range of 2θ = 2.07–35.02° for MMT and 2θ = 2.11–35.06° for OMMT. The same OMMT was used to prepare PBT-PCN and nylon6-PCN.  For the PBT-PCN film a scan rate of 0.02°/s and the scan range of 2θ = 2.03–65.01° are used. The Nylon6-PCN is compressed into thin sheets by compression testing machine (SATEC, model 22 EMF, PA, USA) and the X-ray diffraction of Nylon6-PCN was done at a scan rate of 4°/min with the scan range of 2θ = 2.01–15°.


### 2.5.2. Photoacoustic- (PA-)FTIR Characterization

PA-FTIR experiments are conducted in FTIR bench (model Nexus, Thermo Nicolet, USA) with a photoacoustic cell (MTEC Photoacoustics model no. 300). Sample is prepared by pressing the OMMT powder into a pellet. The nylon6 sample is heated gradually from 303 K to 493 K and then pressed into a pellet. The PBT and PBT-PCN samples are prepared in smaller size and placed in the sample holder. The PA-FTIR is conducted in the linear photoacoustic mode for 500 scans at a resolution of4cm-1, with moving mirror velocity of 0.1581 cm/s. FTIR data are acquired from 400 to 4000 cm-1.


### 2.5.3. Differential Scanning Calorimetry (DSC)

DSC experiments are conducted using differential scanning calorimeter (model Q1000 DSC) of TA instruments, New Castle, Delaware, USA. PBT and PBT-PCN samples are heated from 25°C to 250°C with nitrogen atmosphere with the heating rate of 10°C/min for heating/cooling/heating cycle with retention time of 3 minutes at the end of each treatment for temperature equilibration. Nylon6 and Nylon6-PCN samples are heated from 30°C to 100°C in 3 minutes and the temperature is maintained at 100°C for 10 minutes. A heating rate of 10°C/min is used for the first heating and cooling cycles with retention time of 10 minutes at 250°C and the heating rate of 20°C/min [11] is used in the final heating cycle. We have taken the net weight of the respective polymers into consideration. The enthalpy of fusion was determined by integrating the area under the DSC endotherm peak which was divided by the net weight of the material. The percentage of crystallinity Xc is determined from DSC measurements by
(1)Xc=ΔHmΔHm°100%,
where ΔH is melting enthalpy of sample [J/g]; ΔHm° represents theoretical value of enthalpy for 100% crystalline polymer. For PBT homopolymer, ΔHm°=141 J/g [29]; for nylon6, ΔHm°=203 J/g [11]. This is generally considered a standard method to measure crystallinity of polymer which has been described earlier [30] and extensively used in literature [10, 29, 31, 32].


### 2.5.4. Nanoindentation

Nanoindentation experiments are conducted using triboscope nanoindenter (Hysitron Inc., Minneapolis, MN). The triboscope is operated in association with a multimode AFM controlled by the nanoscope IIIa controller and with the J-type piezo scanner from Veeco (Veeco Metrology Group, Santa Barbara, CA). For PBT and PBT-PCN, load controlled nanoindentation tests are performed using a 100–200 nm Berkovich (three-sided pyramid) diamond indenter tip. Experiments are conducted for maximum load of 2900μN. The reduced elastic modulus (Er) and hardness (H) values are calculated from initial part of the unloading curve using Oliver-Pharr method [33]. The average Er and H are obtained from 30 numbers of indentations for each specimen. For nylon6 and Nylon6-PCN, load controlled nanoindentation tests are performed [11]. The elastic modulus is calculated using the following equation:(2)1Er=1-v2E+1-vi2Ei,
where E and ν are the elastic modulus and Poisson’s ratio of the sample respectively. Ei and νi are the elastic modulus and Poisson’s ratio for the diamond indenter. In the case of polymeric materials and its composites, due to the fact that Ei≫E, (2) can be simply given as [34]                                                       (3)E=Er(1-ν2).


## 2.5.1. XRD Characterization

The X-ray analysis is performed by X-ray diffractometer (Philips Analatical X’pert MPD, Almelo, Netherlands) with secondary monochromator and Cu-tube by CuKα radiation wavelength of 1.5406 Å. X-ray diffraction of MMT and OMMT is conducted at a scan rate of 0.01°/s in each case and with scan range of 2θ = 2.07–35.02° for MMT and 2θ = 2.11–35.06° for OMMT. The same OMMT was used to prepare PBT-PCN and nylon6-PCN.  For the PBT-PCN film a scan rate of 0.02°/s and the scan range of 2θ = 2.03–65.01° are used. The Nylon6-PCN is compressed into thin sheets by compression testing machine (SATEC, model 22 EMF, PA, USA) and the X-ray diffraction of Nylon6-PCN was done at a scan rate of 4°/min with the scan range of 2θ = 2.01–15°.


## 2.5.2. Photoacoustic- (PA-)FTIR Characterization

PA-FTIR experiments are conducted in FTIR bench (model Nexus, Thermo Nicolet, USA) with a photoacoustic cell (MTEC Photoacoustics model no. 300). Sample is prepared by pressing the OMMT powder into a pellet. The nylon6 sample is heated gradually from 303 K to 493 K and then pressed into a pellet. The PBT and PBT-PCN samples are prepared in smaller size and placed in the sample holder. The PA-FTIR is conducted in the linear photoacoustic mode for 500 scans at a resolution of4cm-1, with moving mirror velocity of 0.1581 cm/s. FTIR data are acquired from 400 to 4000 cm-1.


## 2.5.3. Differential Scanning Calorimetry (DSC)

DSC experiments are conducted using differential scanning calorimeter (model Q1000 DSC) of TA instruments, New Castle, Delaware, USA. PBT and PBT-PCN samples are heated from 25°C to 250°C with nitrogen atmosphere with the heating rate of 10°C/min for heating/cooling/heating cycle with retention time of 3 minutes at the end of each treatment for temperature equilibration. Nylon6 and Nylon6-PCN samples are heated from 30°C to 100°C in 3 minutes and the temperature is maintained at 100°C for 10 minutes. A heating rate of 10°C/min is used for the first heating and cooling cycles with retention time of 10 minutes at 250°C and the heating rate of 20°C/min [11] is used in the final heating cycle. We have taken the net weight of the respective polymers into consideration. The enthalpy of fusion was determined by integrating the area under the DSC endotherm peak which was divided by the net weight of the material. The percentage of crystallinity Xc is determined from DSC measurements by
(1)Xc=ΔHmΔHm°100%,
where ΔH is melting enthalpy of sample [J/g]; ΔHm° represents theoretical value of enthalpy for 100% crystalline polymer. For PBT homopolymer, ΔHm°=141 J/g [29]; for nylon6, ΔHm°=203 J/g [11]. This is generally considered a standard method to measure crystallinity of polymer which has been described earlier [30] and extensively used in literature [10, 29, 31, 32].


## 2.5.4. Nanoindentation

Nanoindentation experiments are conducted using triboscope nanoindenter (Hysitron Inc., Minneapolis, MN). The triboscope is operated in association with a multimode AFM controlled by the nanoscope IIIa controller and with the J-type piezo scanner from Veeco (Veeco Metrology Group, Santa Barbara, CA). For PBT and PBT-PCN, load controlled nanoindentation tests are performed using a 100–200 nm Berkovich (three-sided pyramid) diamond indenter tip. Experiments are conducted for maximum load of 2900μN. The reduced elastic modulus (Er) and hardness (H) values are calculated from initial part of the unloading curve using Oliver-Pharr method [33]. The average Er and H are obtained from 30 numbers of indentations for each specimen. For nylon6 and Nylon6-PCN, load controlled nanoindentation tests are performed [11]. The elastic modulus is calculated using the following equation:(2)1Er=1-v2E+1-vi2Ei,
where E and ν are the elastic modulus and Poisson’s ratio of the sample respectively. Ei and νi are the elastic modulus and Poisson’s ratio for the diamond indenter. In the case of polymeric materials and its composites, due to the fact that Ei≫E, (2) can be simply given as [34]                                                       (3)E=Er(1-ν2).


## 2.6. Simulations Using Molecular Dynamics

Molecular dynamics (MD) is a computational method used to study the position and interaction of atoms and molecules. MD is used extensively to study proteins, bio molecules, polymers, and minerals. MD provides a good understanding of molecules with respect to their interactions and molecular structure. This method involves solving the equation of motion in small time steps. The equation of motion isFi=miai, where Fi is the force acting on an atom i, mi is the mass of an atom i, and ai is the acceleration of an atom i. Change in energy of an atom i with respect to its position (r) is equal to the force acting on an atom i. By solving these equations of motion, MD provides time-dependent response for the molecular system. This response contains position of every atom in the system and their velocities. In this work, MD has been used to study the interactions between atoms and group of atoms in PBT-PCN and Nylon6-PCN.


### 2.6.1. Force Field Parameters for PBT, Nylon6, and Na-MMT

Force field describes the potential energy of a molecular system. Potential energy of any molecular system consists of bonded energy and nonbonded energy. Bonded energy consists of bond, angle, dihedral, and improper energy of the molecular system. Nonbonded energy consists of van der Waals and electro static energy of the molecular system. Total potential energy of the molecular system can be described by the following expressionEtotal=[Ebond+Eangle+Edihedral+Eimproper]+[Evan+Eelec]. These individual potential energy terms can be described by mathematical functions and the functions depend on force constants and variables. Some of the general force fields used by researchers are CHARMM, CFF, AMBER, and COMPASS. In our previous work, CHARMM [41] (Chemistry at HARvard Macromolecular Mechanics) force field parameters were used for Polyamide6 and 12-Aminolauric acid and the same set of parameters were used in the present work. We have transformed CFF force field parameters of sodium Montmorillonite to CHARMM force field parameters [42–44]. Force field parameters for PBT were obtained by considering number of methyl group and model components which are similar to polyethylene terephthalate [45].


### 2.6.2. Molecular Model Construction and Simulations

Materials Studio software, developed by Accelrys Software Inc., is used to build the molecular models. This software is widely used to construct molecular models. NAMD2.5b2 [46] and VMD1.8.7 [47] are used to conduct simulations of OMMT, PBT, Nylon6, and PCNs and visualize the results from simulation, respectively. These programs are developed by Theoretical Biophysics Group at Illinois Beckman Institute and widely used by researchers to study complex molecular system.Organically Modified Montmorillonite (OMMT) Model
OMMT model consists of two sheets of montmorillonite (MMT) and eighteen numbers of 12-aminolauric acid molecules. Size of the OMMT model is 31.68 Å inX direction, 27.44 Å in Y direction, and 24.16 Å in Z direction. Each montmorillonite sheet consists of eighteen unit cells where six unit cells are in X direction and three unit cells are in Y direction. Each of these MMT sheets has nine negative net charges and sheets are placed at 17.60 Å initial d-spacing. In our previous work [26], we have studied different initial d-spacing and we found that 17.60 Å initial d-spacing model provides closest experimental d-spacing at the end of simulation. In order to make OMMT model charge neutral, eighteen numbers of 12-aminolauric acid chains are inserted in the model where nine are inserted in between two MMT sheets and nine are inserted above the MMT sheet as shown in Figure 1. Previous work [25] showed that these 12-aminolauric acid molecules should be placed flat and parallel to the MMT sheet. Partial charges and structure of 12-aminolauric acid are shown in Figure 2(a).Figure 1
Initial OMMT model.Distribution of partial charges (a) 12-aminolauric acid, (b) PBT, (c) Nylon6.
(a)
(b)
(c)Polybutylene Terephthalate (PBT) Model
PBT monomers are used to build polymers. A representative molecular model of PBT is shown in Figure3. Six different lengths of PBT polymers were built to study the length effect of polymers on d-spacing of clay in PCN. These polymer chains contain two, four, six, eight, ten, and twelve monomers. All PBT polymers were annealed at 700 K under atmospheric pressure. Annealing process was conducted three times to obtain globally minimized confirmation of PBT chains and minimization was confirmed from the potential energy.Figure 3
Single chain molecular model of PBT.


### 2.6.3. Simulation of OMMT

Initial OMMT model is minimized for 25 ps at a temperature of zero Kelvin in vacuum using molecular dynamics program NAMD. Next, the temperature of the OMMT system is raised to 300 K in three equal steps by Langevin temperature parameter. Once the temperature is increased to 300 K, pressure of the OMMT system is increased to 1.01325 bars in four equal steps by Langevin piston pressure control. Next, the system temperature is increased by 33 K and subsequently reduced to atmospheric temperature. The increase and reduction of temperature is done to mimic the experimental conditions as mentioned in Section2.2. A 0.5 fs time step is used for all OMMT simulations. Finally, OMMT model is simulated for 100 ps at room temperature and pressure with periodic boundary condition. During OMMT simulation, constraints were applied to MMT atoms in X and Y direction but MMT atoms are allowed to move in Z direction. There are no constraints applied to organic modifiers thus organic modifier atoms are allowed to move in all three directions. Simulation of OMMT is done under constant number of particles, constant pressure, and constant temperature (NPT). 16 Å cut-off distance and 14 Å switch distance are used for van der waals and electrostatic energy calculation in all simulations.


### 2.6.4. Polymer Clay Nanocomposites (PCN) Models

Selected OMMT model, as shown in Figure4, and minimized PBT models are used to construct PCN models. Six different lengths of PBT chains are used in PCN to select the suitable PCN system. One of the initial PCN models is shown in Figure 5. The suitable PCN system is selected based on experimental d-spacing as shown in Figure 6. The estimated mass of PBT chains in PCN model with respect to clay is shown in Table 1. Two PBT chains are inserted into the OMMT model and one is placed in between clay sheets and the other is placed above top clay sheet to ensure the periodic boundary condition during the simulation. Physical dimensions of PBT are measured and compared with the simulated interlayer spacing of OMMT to make sure that the PBT chain fits in the model. The same conditions described above were implemented to construct the nylon6 PCN model. The number of Nylon6 chains in PCN models is selected based on our previous work [11]. Initial and final Nylon6-PCN model are shown in Figures 7, and 8 respectively. Final model is selected based on the experimental d-spacing.Table 1
Estimated mass percentage of PBT with respect to clay.
Model
Intercalated polymer wt% with clay
PCN(1)
1.96
PCN(2)
4.19
PCN(3)
6.41
PCN(4)
8.64
PCN(5)
10.87
PCN(6)
13.10Figure 4
Final OMMT model.Figure 5
Initial PBT-PCN model.Figure 6
Final PBT-PCN model.Figure 7
Initial Nylon6-PCN model.Figure 8
Final Nylon6-PCN model.


### 2.6.5. Simulation of PCNs

Molecular simulation of PCN is conducted for PBT-PCN and Nylon6-PCN system. Procedure for minimization and dynamics of PCNs is done as described in Section2.6.3 Cell basis vector for PCNs is adjusted according to the physical size of different models. Periodic boundary conditions are applied to calculate the electrostatic interactions between PCN constituents. Initially, PCN models are simulated for 25 ps. Based on d-spacing, final PCN model is selected and then the final simulation is conducted for 200 ps. Constraints are applied to MMT as mentioned in Section 2.6.3 and polymers are allowed to move freely in all three directions. Constant number, pressure, and temperature (NPT) are used for all PCN simulation. Constant pressure is maintained by Nose’-Hoover Langevin piston pressure control method. This method is used with Langevin dynamics, such as temperature control, in order to perform NPT ensemble. Equilibrium of PCN system is confirmed from the plot of total energy of the system versus simulation time and it is found that 200 ps simulation time is sufficient for the PCN simulations.


## 2.6.1. Force Field Parameters for PBT, Nylon6, and Na-MMT

Force field describes the potential energy of a molecular system. Potential energy of any molecular system consists of bonded energy and nonbonded energy. Bonded energy consists of bond, angle, dihedral, and improper energy of the molecular system. Nonbonded energy consists of van der Waals and electro static energy of the molecular system. Total potential energy of the molecular system can be described by the following expressionEtotal=[Ebond+Eangle+Edihedral+Eimproper]+[Evan+Eelec]. These individual potential energy terms can be described by mathematical functions and the functions depend on force constants and variables. Some of the general force fields used by researchers are CHARMM, CFF, AMBER, and COMPASS. In our previous work, CHARMM [41] (Chemistry at HARvard Macromolecular Mechanics) force field parameters were used for Polyamide6 and 12-Aminolauric acid and the same set of parameters were used in the present work. We have transformed CFF force field parameters of sodium Montmorillonite to CHARMM force field parameters [42–44]. Force field parameters for PBT were obtained by considering number of methyl group and model components which are similar to polyethylene terephthalate [45].


## 2.6.2. Molecular Model Construction and Simulations

Materials Studio software, developed by Accelrys Software Inc., is used to build the molecular models. This software is widely used to construct molecular models. NAMD2.5b2 [46] and VMD1.8.7 [47] are used to conduct simulations of OMMT, PBT, Nylon6, and PCNs and visualize the results from simulation, respectively. These programs are developed by Theoretical Biophysics Group at Illinois Beckman Institute and widely used by researchers to study complex molecular system.Organically Modified Montmorillonite (OMMT) Model
OMMT model consists of two sheets of montmorillonite (MMT) and eighteen numbers of 12-aminolauric acid molecules. Size of the OMMT model is 31.68 Å inX direction, 27.44 Å in Y direction, and 24.16 Å in Z direction. Each montmorillonite sheet consists of eighteen unit cells where six unit cells are in X direction and three unit cells are in Y direction. Each of these MMT sheets has nine negative net charges and sheets are placed at 17.60 Å initial d-spacing. In our previous work [26], we have studied different initial d-spacing and we found that 17.60 Å initial d-spacing model provides closest experimental d-spacing at the end of simulation. In order to make OMMT model charge neutral, eighteen numbers of 12-aminolauric acid chains are inserted in the model where nine are inserted in between two MMT sheets and nine are inserted above the MMT sheet as shown in Figure 1. Previous work [25] showed that these 12-aminolauric acid molecules should be placed flat and parallel to the MMT sheet. Partial charges and structure of 12-aminolauric acid are shown in Figure 2(a).Figure 1
Initial OMMT model.Distribution of partial charges (a) 12-aminolauric acid, (b) PBT, (c) Nylon6.
(a)
(b)
(c)Polybutylene Terephthalate (PBT) Model
PBT monomers are used to build polymers. A representative molecular model of PBT is shown in Figure3. Six different lengths of PBT polymers were built to study the length effect of polymers on d-spacing of clay in PCN. These polymer chains contain two, four, six, eight, ten, and twelve monomers. All PBT polymers were annealed at 700 K under atmospheric pressure. Annealing process was conducted three times to obtain globally minimized confirmation of PBT chains and minimization was confirmed from the potential energy.Figure 3
Single chain molecular model of PBT.


## 2.6.3. Simulation of OMMT

Initial OMMT model is minimized for 25 ps at a temperature of zero Kelvin in vacuum using molecular dynamics program NAMD. Next, the temperature of the OMMT system is raised to 300 K in three equal steps by Langevin temperature parameter. Once the temperature is increased to 300 K, pressure of the OMMT system is increased to 1.01325 bars in four equal steps by Langevin piston pressure control. Next, the system temperature is increased by 33 K and subsequently reduced to atmospheric temperature. The increase and reduction of temperature is done to mimic the experimental conditions as mentioned in Section2.2. A 0.5 fs time step is used for all OMMT simulations. Finally, OMMT model is simulated for 100 ps at room temperature and pressure with periodic boundary condition. During OMMT simulation, constraints were applied to MMT atoms in X and Y direction but MMT atoms are allowed to move in Z direction. There are no constraints applied to organic modifiers thus organic modifier atoms are allowed to move in all three directions. Simulation of OMMT is done under constant number of particles, constant pressure, and constant temperature (NPT). 16 Å cut-off distance and 14 Å switch distance are used for van der waals and electrostatic energy calculation in all simulations.


## 2.6.4. Polymer Clay Nanocomposites (PCN) Models

Selected OMMT model, as shown in Figure4, and minimized PBT models are used to construct PCN models. Six different lengths of PBT chains are used in PCN to select the suitable PCN system. One of the initial PCN models is shown in Figure 5. The suitable PCN system is selected based on experimental d-spacing as shown in Figure 6. The estimated mass of PBT chains in PCN model with respect to clay is shown in Table 1. Two PBT chains are inserted into the OMMT model and one is placed in between clay sheets and the other is placed above top clay sheet to ensure the periodic boundary condition during the simulation. Physical dimensions of PBT are measured and compared with the simulated interlayer spacing of OMMT to make sure that the PBT chain fits in the model. The same conditions described above were implemented to construct the nylon6 PCN model. The number of Nylon6 chains in PCN models is selected based on our previous work [11]. Initial and final Nylon6-PCN model are shown in Figures 7, and 8 respectively. Final model is selected based on the experimental d-spacing.Table 1
Estimated mass percentage of PBT with respect to clay.
Model
Intercalated polymer wt% with clay
PCN(1)
1.96
PCN(2)
4.19
PCN(3)
6.41
PCN(4)
8.64
PCN(5)
10.87
PCN(6)
13.10Figure 4
Final OMMT model.Figure 5
Initial PBT-PCN model.Figure 6
Final PBT-PCN model.Figure 7
Initial Nylon6-PCN model.Figure 8
Final Nylon6-PCN model.


## 2.6.5. Simulation of PCNs

Molecular simulation of PCN is conducted for PBT-PCN and Nylon6-PCN system. Procedure for minimization and dynamics of PCNs is done as described in Section2.6.3 Cell basis vector for PCNs is adjusted according to the physical size of different models. Periodic boundary conditions are applied to calculate the electrostatic interactions between PCN constituents. Initially, PCN models are simulated for 25 ps. Based on d-spacing, final PCN model is selected and then the final simulation is conducted for 200 ps. Constraints are applied to MMT as mentioned in Section 2.6.3 and polymers are allowed to move freely in all three directions. Constant number, pressure, and temperature (NPT) are used for all PCN simulation. Constant pressure is maintained by Nose’-Hoover Langevin piston pressure control method. This method is used with Langevin dynamics, such as temperature control, in order to perform NPT ensemble. Equilibrium of PCN system is confirmed from the plot of total energy of the system versus simulation time and it is found that 200 ps simulation time is sufficient for the PCN simulations.


## 3. Results and Discussion

### 3.1. Experiments

#### 3.1.1. XRD Results

d-spacing of MMT, OMMT, and PCN from several samples are summarized in Table 2. MMT d-spacing is less than the OMMT d-spacing thus indicating that the organic modifier has intercalated the clay. d-spacing of both PCNs lie in between pure MMT d-spacing and OMMT d-spacing. Intercalation of nylon6 with clay was verified from the XRD, spectroscopic, and simulation studies [26, 27]. PBT-PCN has intercalated structure which has higher d-spacing than nylon6-PCN. The reduction of d-spacing in PCN compared to OMMT might be caused by dissolution of the modifier in the solvents.Table 2
Experimentald-spacing.
Sample
d-spacing range/Å
Pure MMT
11.56–12.25
OMMT
15.60–17.43
PBT-PCN
13.90–14.53
Nylon6-PCN
13.50


#### 3.1.2. FTIR Results of PBT, OMMT, and PCN

In our previous work, we have observed that the organic modifier, 12-aminolauric acid, is intercalated between clay sheets with hydrogen bonds and ionic interactions generated by end functional groups of organic modifier with surface oxygen of interlayer clay sheets. Nanocomposites are prepared by intercalation of polymer and MMT with the aid of organic modifier [26]. The spectroscopic analysis on nylon6-PCN system has been reported and revealed the hydrogen formation between intercalated polymer and organic modifier and surface oxygen of clay sheet lying in the interlayer clay gallery in the polymer clay nanocomposite [26]. In the present work, the same OMMT is used in conjunction with PBT and the spectroscopic results of PBT, OMMT, and PBT-PCN are discussed here. The PA FTIR spectra of OMMT, PBT, and PCN are shown in Figure 9 and the corresponding band assignments are shown in Table 3. As indicated in Figure 9 and Table 3, the band contribution in PCN from PBT is much more prominent than that from OMMT due to low weight fraction of OMMT in PCN as well as the strong bands from PBT. FTIR spectra of PCN, PBT, and OMMT are shown in Figure 10 in the range of 3800–2800 cm-1. It is observed that the sharp O–H stretching band at 3645–3628 cm-1 has affected the spectrum of PCN and created a broad band in the range of 3633–3679 cm-1 and the broad band of ionically bonded N–H stretching at 3256–3035 cm-1 generating a bulge at 3260–3047 cm-1 in PCN associated with aromatic C–H stretching from PBT. FTIR spectra of PCN, PBT, and OMMT are shown in Figure 11 in between 2000 and 1000 cm-1. In this energy-region, bands from PBT are very strong and thus the bands in both PBT and PCN are almost identical. Most bands from OMMT are overlapped by PBT, and Si–O stretching band is not visible in PCN. The PA FTIR spectra of PCN, PBT, and OMMT in the range of 1000–400 cm-1 are shown in Figure 12. The in-plane ring deformation vibration of PBT at 633 cm-1 has shifted to the lower energy side by 9 cm-1 in the PCN. The tripled Al–O, Si–O vibration and in-plane ring deformation vibration shows broader band than that of pure PBT. From Figure 12, it can be seen that the sharp out-of-plane ring deformation vibration of PBT at 497 cm-1 is suppressed and the out-of-plane ring deformation vibration of PBT at 476 cm-1 is absent in PCN. These phenomena indicate that ring deformations in PBT are restricted by OMMT due to intercalation. It is evident that Si–O–Si deformation from OMMT appears in PCN at 463 cm-1 and shifted to the lower energy side by 15 cm-1 and this shift in PCN is probably due to the formation of bonding of PBT with silica tetrahedral oxygen through the hydrogen bond between the functional groups of PBT and 12-aminolauric acid. Due to this effect the Si–O stretching band in PCN disappears. Wang et al. [36] found that during melt processing when PBT/PET blends are filled with silica nanoparticles, hydroxyl end groups of PET and PBT react with the surface hydroxyl groups of SiO2 which is shown by X-ray photoelectron Spectroscopy (XPS) tests. Fermeglia et al. [48] investigated the binding energy between poly(Butylene terephtalate-co-thiodiethylene terephtalate), clay and different organic modifiers using molecular simulations and have found that –NH2, –COOH, or –SH groups in organic modifiers generally lead to a good interaction between modifiers and polymers.Table 3
Band assignments of photoacoustic FTIR spectra of OMMT, PBT, and PCN.
Sample
Wave number range (cm-1)
Bands
Reference
3646–3628
O–H stretching of structural hydroxyl group
[26]
3256–3035
Ionic bonded N–H stretching
2932
C–H asymmetric stretching
2853
C–H symmetric stretching
1713–1701
C=O stretching
1623
Combination of O–H deformation and N–H bending
1511
Asymmetric R–COO- stretching
1466
CO–H bending
1396
Symmetric R–COO- stretching
1112
Si–O in-plane stretching
1060
Si–O out-of-plane stretching
OMMT
1022,997
Si–O in-plane stretching
963
Si–OH vibrations
921
Al–O/Al–OH stretching vibration
884
Al–FeOH deformation
847
Al–MgOH deformation
799
Si–O stretching of quartz and silica
727
Si–O deformation perpendicular to optical axis
699
Si–O deformation parallel to optical axis
623
Coupled Al–O and Si–O out-of-plane vibration
531
Al–OSi deformation
478
Si–OSi vibration
3417
–OH stretching
3052
Aromatic C–H stretching
[35]
2962
C–H2 stretching
[36]
2895
Aliphatic C–H asymmetric stretching
[35]
2863
Aliphatic C–H symmetric stretching
[35]
1722
C=O stretching
[37]
1583
Aromatic C=C stretching
[38]
1504
Aliphatic C–H2 bending
[37]
1453
Aliphatic C–H2 bending
[39]
1407
Aliphatic aromatic C–H in-of-plane bending
[36]
PBT
1387
CH2 wagging
[39]
1281
(C=O–)–O stretching
[35]
1110
CH2–O stretching
[35]
1017
Aromatic C–H bending
[35]
873
Aromatic C–H bending
[40]
808
Aromatic C–H bending
[37]
729
Aromatic C–H out-of-plane bending
[35, 36]
685
Aromatic =C–H out-of-plane deformation
[38]
633
Aromatic in-plane deformation
[38]
497
Aromatic out-of-plane deformation
[38]
476
Aromatic out-of-plane deformation
[38]
3633–3679
O–H stretching of  structural hydroxyl group-OMMT
[26]
3413
–OH stretching
3260–3047
Coupled N–H stretching and aromatic C–H stretching
[35]
2963
C–H2 stretching
[36]
2895
Aliphatic C–H asymmetric stretching
[35]
2858
Aliphatic C–H symmetric stretching
[35]
1719
C=O stretching
[37]
1576
Aromatic C=C stretching
[38]
1503
Aliphatic C–H2 bending
[37]
1452
Aliphatic C–H2 bending
[39]
1406
Aromatic C–H in-of-plane bending
[36]
PCN
1387
CH2 wagging
[39]
1278
(C=O–)–O stretching
[35]
1105
CH2–O stretching
[35]
1020
Aromatic C–H bending
[35]
873
Aromatic C–H bending
[40]
810
Aromatic C–H bending
[37]
729
Aromatic C–H out-of-plane bending
[35, 36]
685
Aromatic =C–H out-of-plane deformation
[38]
652
Aromatic in-plane deformation
[38]
624
Tripled Al–O, Si–O out-of-plane and aromatic in-plane deformation
[26]
497
Aromatic out-of-plane deformation
[38]
463
Si–OSi deformation
[26]Figure 9
Photo acoustic FTIR spectra of OMMT, PBT, and PCN in the energy range of 4000–400cm-1 at scan speed of 0.15 cm-1.Figure 10
Photo acoustic FTIR spectra of OMMT, PBT, and PCN in the energy range of 3800–2800cm-1 at scan speed of 0.15 cm-1.Figure 11
Photo acoustic FTIR spectra of OMMT, PBT, and PCN in the energy range of 2000–1000cm-1 at scan speed of 0.15 cm-1.Figure 12
Photo acoustic FTIR spectra of OMMT, PBT, and PCN in the energy range of 1000–400cm-1 at scan speed of 0.15 cm-1.


#### 3.1.3. Crystallinity of PBT and PCN

DSC tests are conducted on the pure PBT and PCN and results are shown in Table4; it is observed that the pure PBT shows 0.5% higher crystallinity than that of PCN. This observation of reduced crystallinity of PCN with respect to a polymer has been observed earlier in Nylon6-PCN [10]. As shown in Table 4, Nylon6-PCN with 9% of clay loading shows 4.56% reduction in crystallinity.Table 4
DSC results of PBT, Nylon6, and PCNs.
Materials
Enthalpy (J/g)
Crystallinity (%)
PBT
50.21
35.6
PBT-PCN
49.56
35.1
Nylon6
56.06
27.61
Nylon6-PCN
46.81
23.05


#### 3.1.4. Nanomechanical Properties of PBT and PCN

In the present work, nanoindentation tests for PBT and PBT-PCN are performed in load control mode. Figure13 is the representative L-D curves for PBT and PBT-PCN films. Maximum load was maintained for period of 5 s between loading and unloading sequence in order to reduce the time-dependent deformation. The elastic modulus and hardness values for PBT and PBT-PCN are presented in Table 5. Nylon6 and Nylon6-PCN nanomechanical tests were performed in load control mode [10] and the results are presented in Table 5.Table 5
Nanomechanical properties of PCNs.
Parameters
PBT
PBT-PCN
Nylon6
Nylon6-PCN
Crystallinity (%)
35.6
35.1
27.6
23.1
Elastic modulus (GPa)
1.40
2.20
3.35
5.46
Hardness (GPa)
0.09
0.14
0.12
0.17Figure 13
Representative load displacement (L-D) curves for PBT and PCN.


### 3.2. Interactions in PCNs System

#### 3.2.1. Nonbonded Interactions between Organic Modifier and Clay

Interaction energies between different constituents of OMMT and PCNs are calculated using NAMDEnergy module of NAMD software. Nonbonded interactions, van der waals and electrostatic, are calculated in terms of attractive and repulsive interactions. Among different initiald-spacing of OMMT, d-spacing of 17.60 Å OMMT model predicts final d-spacing of 15.50 Å which is very close to experimental d-spacing of 15.60 Å. Based on this comparison of d-spacing, the OMMT model is chosen for simulation. Among six different PCNs models, the PCN model which contains 6.41 wt% of PBT predicts the simulated d-spacing of 14.80 Å which is close to experimental d-spacing of 14.53 Å. The same OMMT model with Nylon6 is simulated and it predicts d-spacing of 15.35 Å. Non-bonded interactions between organic modifier and clay are calculated and the interaction energies are shown in Table 6. Positive values indicate repulsive interactions and negative values indicate attractive interactions between constituents and higher magnitude denotes greater interactions. Interaction energy values of organic modifiers are the sum of backbone and functional group interaction energies. Among these two components, functional group of organic modifier shows higher interaction with clay than the backbone of the organic modifier. This is because of higher positive net charge in the functional group compared to the backbone. Functional group interactions are composed of −1096 kcal/mol in electrostatic interaction energy and −80 kcal/mol in van der waals interaction energy and backbone interactions are composed of −549 kcal/mol electrostatic interaction energy and −345 kcal/mol van der waals interaction energy. It is observed that the electrostatic interactions are the major components of both modifier functional group and modifier backbone interaction energies with the clay. More detailed quantitative results are presented in Table 7. Here, it is observed that the major component of attractive interaction with clay is from the hydrogen atoms of modifier. Even though the partial charge of the backbone hydrogen smaller is compared to functional group hydrogen as shown in Figure 2(a), backbone hydrogen shows higher interaction because of its presence in significant numbers.Table 6
Nonbonded interaction energies of OMMT constituents.
OMMT constituents
Nonbonded energy (kcal/mol)
Clay-modifier backbone
−894
Clay-modifier functional group
−1176Table 7
Nonbonded energies of individual atoms of organic modifier with clay.
Component of modifier and clay
Electrostatic energy (kcal/mol) [Col-A]
Van der waals energy (kcal/mol) [Col-B ]
Total nonbonded energy (kcal/mol) [Col-A + Col-B]
Clay-modifier backbone hydrogen
−2503
−87
−2590
Clay-modifier backbone carbon
+1954
−257
+1697
Clay-modifier functional hydrogen
−2392
−8
−2400
Clay-modifier functional nitrogen
+497
−0.5
+497
Clay-modifier functional oxygen
+1689
−48
+1641
Clay-modifier functional carbon
−890
−25
−915


#### 3.2.2. Nonbonded Interactions of PBT-PCN and Nylon6-PCN System

The nonbonded interaction energies between constituents of PBT-PCN and Nylon6-PCN are presented in Table8. Interaction energy between the clay and modifier is attractive and higher than interaction energy between clay and polymer. This is because of higher partial charges of individual atoms in modifiers than the polymers. Interaction energies between clay and polymer are attractive and higher than the interaction energy between polymer and modifier. Detailed results are shown in Table 9 and it is observed that the interaction energy between clay and modifier functional group is attractive and higher than among other constituents in both PCNs. This higher interaction between clay and modifier functional groups could be attributed to higher positive partial charges in functional groups atoms as shown in Figure 2(a). Similarly, modifier backbone also has higher attractive interactions with clay compared to the polymer functional group and backbone in both PCNs.Table 8
Nonbonded energies PCNs.
Component of PCN
Electrostatic energy (kcal/mol)
Van der waals energy (kcal/mol)
Total nonbonded energy (kcal/mol)
[Col-A]
[Col-B ]
[Col-A + Col-B]
PBT-PCN
Clay-polymer
−36
−165
−201
Clay-modifier
−1836
−316
−2152
Polymer-modifier
−93
−43
−136
Nylon6-PCN
Clay-polymer
12
−182
−170
Clay-modifier
−1464
−375
−1839
Polymer-modifier
−209
−30
−239Table 9
Nonbonded energies of modifier and polymer with clay.
Component of modifier, polymer, and clay in PCN
Electrostatic energy (kcal/mol)
Van der waals energy (kcal/mol)
Total nonbonded energy (kcal/mol)
[Col-A]
[Col-B ]
[Col-A + Col-B]
PBT-PCN
Clay-modifier backbone
−582
−282
−864
Clay-modifier functional group
−1255
−34
−1289
Clay-polymer backbone
−250
−108
−358
Clay-polymer functional group
+213
−56
+157
Nylon6-PCN
Clay-modifier backbone
−485
−297
−782
Clay-modifier functional group
−979
−77
−1056
Clay-polymer backbone
−183
−123
−306
Clay-polymer functional group
+155
−62
+93Polymer functional group has negative partial charges as shown in Figure2(b) thus it induces repulsive interaction with clay in both PCNs systems. Interaction energies between clay and individual atoms of organic modifier are shown in Table 10. It is observed that the organic modifier hydrogen atoms show higher magnitude of attractive interactions with clay as observed in the OMMT system. Interaction energies of individual polymer atoms with clay are shown in Table 11. Here, hydrogen and carbon atoms show attractive interaction with clay while the other atoms show repulsive interactions. Functional group carbon atoms contribute towards more attractive interactions and greater negative partial charge of oxygen atoms contribute to more repulsive interactions. As shown in Table 12, polymer backbone has higher repulsive interaction with modifier functional group and lesser repulsive interaction with modifier backbone. Polymer functional group has higher attractive interaction with modifier functional group and lesser attractive interaction with backbone. Similar nature of interactions have been reported in our previous work [11].Table 10
Nonbonded energies between individual atoms of organic modifier and clay.
Individual atoms of organic modifiers and clay
Electrostatic energy (kcal/mol)
Van der waals energy (kcal/mol)
Total nonbonded energy (kcal/mol)
[Col-A]
[Col-B ]
[Col-A + Col-B]
PBT-PCN
Clay-modifier backbone hydrogen
−2512
−67
−2579
Clay-modifier backbone carbon
+1930
−215
+1715
Clay-modifier functional hydrogen
−2564
−7
−2571
Clay-modifier functional nitrogen
+546
+37
+583
Clay-modifier functional oxygen
+1628
−42
+1586
Clay-modifier functional carbon
−865
−22
−887
Nylon6-PCN
Clay-modifier backbone hydrogen
−2085
−72
−2157
Clay-modifier backbone carbon
+1600
−226
+1374
Clay-modifier functional hydrogen
−2071
−7
−2078
Clay-modifier functional nitrogen
+444
−2
+442
Clay-modifier functional oxygen
+1359
−44
+1315
Clay-modifier functional carbon
−711
−24
−735Table 11
Nonbonded energies between individual atoms of polymer and clay.
Component of polymer and clay in PCN
Electrostatic energy (kcal/mol) [Col-A]
Van der waals energy (kcal/mol) [Col-B ]
Total nonbonded energy (kcal/mol) [Col-A + Col-B]
PBT-PCN
Clay-polymer aromatic  hydrogen
−301
−10
−311
Clay-polymer backbone hydrogen
−380
−10
−390
Clay-polymer aromatic  carbon
+188
−59
+129
Clay-polymer backbone carbon
+243
−29
+214
Clay-polymer functional hydrogen
−90
−0.2
−90.2
Clay-polymer functional oxygen
+1127
−36
1091
Clay-polymer functional carbon
−824
−20
−844
Nylon6-PCN
Clay-polymer backbone hydrogen
−1056
−19
−1075
Clay-polymer backbone carbon
+872
−104
+768
Clay-polymer functional nitrogen
+593
−27
+566
Clay-polymer functional hydrogen
−424
−2
−426
Clay-polymer functional oxygen
+589
−13
+576
Clay-polymer functional carbon
−600
−20
−620
Clay-polymer backbone hydrogen
−1056
−19
−1075Table 12
Nonbonded energies between polymer and organic modifier.
Component of modifier and polymer in PCN
Electrostatic energy (kcal/mol)[Col-A]
Van der waals energy (kcal/mol)[Col-B ]
Total nonbonded energy (kcal/mol)       [Col-A + Col-B]
PBT-PCN
Polymer backbone—modifier backbone
+115
−29
+86
Polymer backbone—modifier functional
+209
−10
+199
Polymer functional—modifier backbone
−132
−12
−144
Polymer functional—modifier functional
−285
+7
−278
Nylon6-PCN
Polymer backbone—modifier backbone
+58
−20
+38
Polymer backbone—modifier functional
+91
−11
+80
Polymer functional—modifier backbone
−94
−11
−105
Polymer functional—modifier functional
−264
+11
−253


#### 3.2.3. Summary of Interaction of OMMT, PBT-PCN, and Nylon6-PCN System

As mentioned in Section1, nonbonded interactions play a major role in PCNs and they affect the crystallinity and mechanical properties of PCNs. Based on our results, more than 80%, 75%, 70% of interaction energies are a result of electrostatic interaction for OMMT, PBT-PCN and Nylon6-PCN respectively. Figures 14 and 15 show overall interaction energy map for PBT-PCN, and Nylon6-PCN, respectively. The two energy maps have similar pattern in terms of attractive and repulsive interactions between different constituents of PCNs. These energy maps indicate that the PBT-PCN system has greater magnitude of interaction energies between different constituents than the Nylon6-PCN system. In all three systems, OMMT, PBT-PCN, and Nylon6-PCN, modifier shows attractive interactions with the clay. Further, in PBT-PCN system, modifiers have 17% higher attractive interactions with clay than in Nylon6-PCN system. Polymer has attractive interaction with clay and modifier in both systems. In terms of interaction energy magnitudes, almost equal interaction of polymer with clay is observed in both systems and 75 % higher interaction with organic modifier is observed in Nylon6-PCN compared to PBT-PCN.Figure 14
Energy map of PBT-PCN system.Figure 15
Energy map of nylon6-PCN system.


### 3.3. Discussion

Our previous work has shown that the molecular interactions between clay, modifier, and polymer results in significant alteration of polymer to a very large extent in PCN [10]. This is the basis of the altered phase theory to explain the significant improvement in mechanical properties of PCN as a result of addition of small amounts of well-dispersed clay nanoparticles in the polymer. The attractive interactions between functional groups of the polymer with the modifier and concurrent repulsive interaction between modifier and polymer backbone result in reduction of crystallinity of the surrounding polymer.In the current work with PBT-PCN system we have observed reduction in crystallinity and improved mechanical properties of the PCN with respect to the pristine polymer. We have seen similar phenomena of reduction in crystallinity and improved mechanical properties of Nylon6-PCN [10]; however, the magnitudes of reduction in crystallinity for the two PCNs are different and the improvement in mechanical properties is different. Based on the interaction energy maps for the two systems shown in Figures 14 and 15, the mechanism of attractive interaction between functional group of the polymer and modifier and repulsive interaction between backbone of the polymer and the modifier leading to changes in polymer crystallinity is present, although the magnitudes of the interaction energies are different. In the case of Nylon6-PCN the reduction in crystallinity is 16.5 % while the reduction in crystallinity is 0.5% in PBT-PCN. The improvement in elastic modulus and hardness values for Nylon6-PCN is 63% and 42%, respectively, with respect to pristine Nylon6 while the improvement in elastic modulus and hardness values for PBT-PCN is only 10% and 23%, respectively, with respect to pristine PBT. It appears that in the case of Nylon6-PCN, the clay-modifier system has significantly altered crystallinity and improved the mechanical properties, while the same clay-modifier system has impacted only moderately the crystallinity and mechanical properties of the PBT-PCN system.The comparison of interaction energy maps for the two PCNs shown in Figures14 and  15 reveals that the interaction energies between polymer and modifier and polymer and clay are higher in the PBT-PCN system than the Nylon6-PCN system. The attractive interaction energy and repulsive interaction energy between polymer and modifier are 422 kcal/mol and 285 kcal/mol for PBT-PCN while the attractive and repulsive interaction energies for Nylon6-PCN are 358 kcal/mol and 118 kcal/mol, respectively. The attractive and repulsive interaction energies between clay and polymer are 157 kcal/mol and 358 kcal/mol for PBT-PCN and 93 kcal/mol and 306 kcal/mol for Nylon6-PCN. Inspite of the higher interaction energies between polymer and modifier in PBT-PCN system than the Nylon6-PCN system, the alteration of crystallinity and improvement in mechanical properties is lower for PBT-PCN system. It appears that the initial crystallinity of the polymer may have affected the magnitude of interaction energy needed to alter the crystallinity and improve the mechanical properties of the PCN. The crystallinity of pure PBT is 35.6% while the crystallinity of Nylon6 is 27.6%. Thus, it appears that significantly higher interaction energies between clay-modifier and polymer may be required to alter the crystallinity and improve mechanical properties of a polymer that is more crystalline to begin with. Modest differences are observed in the attractive interaction energies between the same modifier and clay in the presence of different polymers.


## 3.1. Experiments

### 3.1.1. XRD Results

d-spacing of MMT, OMMT, and PCN from several samples are summarized in Table 2. MMT d-spacing is less than the OMMT d-spacing thus indicating that the organic modifier has intercalated the clay. d-spacing of both PCNs lie in between pure MMT d-spacing and OMMT d-spacing. Intercalation of nylon6 with clay was verified from the XRD, spectroscopic, and simulation studies [26, 27]. PBT-PCN has intercalated structure which has higher d-spacing than nylon6-PCN. The reduction of d-spacing in PCN compared to OMMT might be caused by dissolution of the modifier in the solvents.Table 2
Experimentald-spacing.
Sample
d-spacing range/Å
Pure MMT
11.56–12.25
OMMT
15.60–17.43
PBT-PCN
13.90–14.53
Nylon6-PCN
13.50


### 3.1.2. FTIR Results of PBT, OMMT, and PCN

In our previous work, we have observed that the organic modifier, 12-aminolauric acid, is intercalated between clay sheets with hydrogen bonds and ionic interactions generated by end functional groups of organic modifier with surface oxygen of interlayer clay sheets. Nanocomposites are prepared by intercalation of polymer and MMT with the aid of organic modifier [26]. The spectroscopic analysis on nylon6-PCN system has been reported and revealed the hydrogen formation between intercalated polymer and organic modifier and surface oxygen of clay sheet lying in the interlayer clay gallery in the polymer clay nanocomposite [26]. In the present work, the same OMMT is used in conjunction with PBT and the spectroscopic results of PBT, OMMT, and PBT-PCN are discussed here. The PA FTIR spectra of OMMT, PBT, and PCN are shown in Figure 9 and the corresponding band assignments are shown in Table 3. As indicated in Figure 9 and Table 3, the band contribution in PCN from PBT is much more prominent than that from OMMT due to low weight fraction of OMMT in PCN as well as the strong bands from PBT. FTIR spectra of PCN, PBT, and OMMT are shown in Figure 10 in the range of 3800–2800 cm-1. It is observed that the sharp O–H stretching band at 3645–3628 cm-1 has affected the spectrum of PCN and created a broad band in the range of 3633–3679 cm-1 and the broad band of ionically bonded N–H stretching at 3256–3035 cm-1 generating a bulge at 3260–3047 cm-1 in PCN associated with aromatic C–H stretching from PBT. FTIR spectra of PCN, PBT, and OMMT are shown in Figure 11 in between 2000 and 1000 cm-1. In this energy-region, bands from PBT are very strong and thus the bands in both PBT and PCN are almost identical. Most bands from OMMT are overlapped by PBT, and Si–O stretching band is not visible in PCN. The PA FTIR spectra of PCN, PBT, and OMMT in the range of 1000–400 cm-1 are shown in Figure 12. The in-plane ring deformation vibration of PBT at 633 cm-1 has shifted to the lower energy side by 9 cm-1 in the PCN. The tripled Al–O, Si–O vibration and in-plane ring deformation vibration shows broader band than that of pure PBT. From Figure 12, it can be seen that the sharp out-of-plane ring deformation vibration of PBT at 497 cm-1 is suppressed and the out-of-plane ring deformation vibration of PBT at 476 cm-1 is absent in PCN. These phenomena indicate that ring deformations in PBT are restricted by OMMT due to intercalation. It is evident that Si–O–Si deformation from OMMT appears in PCN at 463 cm-1 and shifted to the lower energy side by 15 cm-1 and this shift in PCN is probably due to the formation of bonding of PBT with silica tetrahedral oxygen through the hydrogen bond between the functional groups of PBT and 12-aminolauric acid. Due to this effect the Si–O stretching band in PCN disappears. Wang et al. [36] found that during melt processing when PBT/PET blends are filled with silica nanoparticles, hydroxyl end groups of PET and PBT react with the surface hydroxyl groups of SiO2 which is shown by X-ray photoelectron Spectroscopy (XPS) tests. Fermeglia et al. [48] investigated the binding energy between poly(Butylene terephtalate-co-thiodiethylene terephtalate), clay and different organic modifiers using molecular simulations and have found that –NH2, –COOH, or –SH groups in organic modifiers generally lead to a good interaction between modifiers and polymers.Table 3
Band assignments of photoacoustic FTIR spectra of OMMT, PBT, and PCN.
Sample
Wave number range (cm-1)
Bands
Reference
3646–3628
O–H stretching of structural hydroxyl group
[26]
3256–3035
Ionic bonded N–H stretching
2932
C–H asymmetric stretching
2853
C–H symmetric stretching
1713–1701
C=O stretching
1623
Combination of O–H deformation and N–H bending
1511
Asymmetric R–COO- stretching
1466
CO–H bending
1396
Symmetric R–COO- stretching
1112
Si–O in-plane stretching
1060
Si–O out-of-plane stretching
OMMT
1022,997
Si–O in-plane stretching
963
Si–OH vibrations
921
Al–O/Al–OH stretching vibration
884
Al–FeOH deformation
847
Al–MgOH deformation
799
Si–O stretching of quartz and silica
727
Si–O deformation perpendicular to optical axis
699
Si–O deformation parallel to optical axis
623
Coupled Al–O and Si–O out-of-plane vibration
531
Al–OSi deformation
478
Si–OSi vibration
3417
–OH stretching
3052
Aromatic C–H stretching
[35]
2962
C–H2 stretching
[36]
2895
Aliphatic C–H asymmetric stretching
[35]
2863
Aliphatic C–H symmetric stretching
[35]
1722
C=O stretching
[37]
1583
Aromatic C=C stretching
[38]
1504
Aliphatic C–H2 bending
[37]
1453
Aliphatic C–H2 bending
[39]
1407
Aliphatic aromatic C–H in-of-plane bending
[36]
PBT
1387
CH2 wagging
[39]
1281
(C=O–)–O stretching
[35]
1110
CH2–O stretching
[35]
1017
Aromatic C–H bending
[35]
873
Aromatic C–H bending
[40]
808
Aromatic C–H bending
[37]
729
Aromatic C–H out-of-plane bending
[35, 36]
685
Aromatic =C–H out-of-plane deformation
[38]
633
Aromatic in-plane deformation
[38]
497
Aromatic out-of-plane deformation
[38]
476
Aromatic out-of-plane deformation
[38]
3633–3679
O–H stretching of  structural hydroxyl group-OMMT
[26]
3413
–OH stretching
3260–3047
Coupled N–H stretching and aromatic C–H stretching
[35]
2963
C–H2 stretching
[36]
2895
Aliphatic C–H asymmetric stretching
[35]
2858
Aliphatic C–H symmetric stretching
[35]
1719
C=O stretching
[37]
1576
Aromatic C=C stretching
[38]
1503
Aliphatic C–H2 bending
[37]
1452
Aliphatic C–H2 bending
[39]
1406
Aromatic C–H in-of-plane bending
[36]
PCN
1387
CH2 wagging
[39]
1278
(C=O–)–O stretching
[35]
1105
CH2–O stretching
[35]
1020
Aromatic C–H bending
[35]
873
Aromatic C–H bending
[40]
810
Aromatic C–H bending
[37]
729
Aromatic C–H out-of-plane bending
[35, 36]
685
Aromatic =C–H out-of-plane deformation
[38]
652
Aromatic in-plane deformation
[38]
624
Tripled Al–O, Si–O out-of-plane and aromatic in-plane deformation
[26]
497
Aromatic out-of-plane deformation
[38]
463
Si–OSi deformation
[26]Figure 9
Photo acoustic FTIR spectra of OMMT, PBT, and PCN in the energy range of 4000–400cm-1 at scan speed of 0.15 cm-1.Figure 10
Photo acoustic FTIR spectra of OMMT, PBT, and PCN in the energy range of 3800–2800cm-1 at scan speed of 0.15 cm-1.Figure 11
Photo acoustic FTIR spectra of OMMT, PBT, and PCN in the energy range of 2000–1000cm-1 at scan speed of 0.15 cm-1.Figure 12
Photo acoustic FTIR spectra of OMMT, PBT, and PCN in the energy range of 1000–400cm-1 at scan speed of 0.15 cm-1.


### 3.1.3. Crystallinity of PBT and PCN

DSC tests are conducted on the pure PBT and PCN and results are shown in Table4; it is observed that the pure PBT shows 0.5% higher crystallinity than that of PCN. This observation of reduced crystallinity of PCN with respect to a polymer has been observed earlier in Nylon6-PCN [10]. As shown in Table 4, Nylon6-PCN with 9% of clay loading shows 4.56% reduction in crystallinity.Table 4
DSC results of PBT, Nylon6, and PCNs.
Materials
Enthalpy (J/g)
Crystallinity (%)
PBT
50.21
35.6
PBT-PCN
49.56
35.1
Nylon6
56.06
27.61
Nylon6-PCN
46.81
23.05


### 3.1.4. Nanomechanical Properties of PBT and PCN

In the present work, nanoindentation tests for PBT and PBT-PCN are performed in load control mode. Figure13 is the representative L-D curves for PBT and PBT-PCN films. Maximum load was maintained for period of 5 s between loading and unloading sequence in order to reduce the time-dependent deformation. The elastic modulus and hardness values for PBT and PBT-PCN are presented in Table 5. Nylon6 and Nylon6-PCN nanomechanical tests were performed in load control mode [10] and the results are presented in Table 5.Table 5
Nanomechanical properties of PCNs.
Parameters
PBT
PBT-PCN
Nylon6
Nylon6-PCN
Crystallinity (%)
35.6
35.1
27.6
23.1
Elastic modulus (GPa)
1.40
2.20
3.35
5.46
Hardness (GPa)
0.09
0.14
0.12
0.17Figure 13
Representative load displacement (L-D) curves for PBT and PCN.


## 3.1.1. XRD Results

d-spacing of MMT, OMMT, and PCN from several samples are summarized in Table 2. MMT d-spacing is less than the OMMT d-spacing thus indicating that the organic modifier has intercalated the clay. d-spacing of both PCNs lie in between pure MMT d-spacing and OMMT d-spacing. Intercalation of nylon6 with clay was verified from the XRD, spectroscopic, and simulation studies [26, 27]. PBT-PCN has intercalated structure which has higher d-spacing than nylon6-PCN. The reduction of d-spacing in PCN compared to OMMT might be caused by dissolution of the modifier in the solvents.Table 2
Experimentald-spacing.
Sample
d-spacing range/Å
Pure MMT
11.56–12.25
OMMT
15.60–17.43
PBT-PCN
13.90–14.53
Nylon6-PCN
13.50


## 3.1.2. FTIR Results of PBT, OMMT, and PCN

In our previous work, we have observed that the organic modifier, 12-aminolauric acid, is intercalated between clay sheets with hydrogen bonds and ionic interactions generated by end functional groups of organic modifier with surface oxygen of interlayer clay sheets. Nanocomposites are prepared by intercalation of polymer and MMT with the aid of organic modifier [26]. The spectroscopic analysis on nylon6-PCN system has been reported and revealed the hydrogen formation between intercalated polymer and organic modifier and surface oxygen of clay sheet lying in the interlayer clay gallery in the polymer clay nanocomposite [26]. In the present work, the same OMMT is used in conjunction with PBT and the spectroscopic results of PBT, OMMT, and PBT-PCN are discussed here. The PA FTIR spectra of OMMT, PBT, and PCN are shown in Figure 9 and the corresponding band assignments are shown in Table 3. As indicated in Figure 9 and Table 3, the band contribution in PCN from PBT is much more prominent than that from OMMT due to low weight fraction of OMMT in PCN as well as the strong bands from PBT. FTIR spectra of PCN, PBT, and OMMT are shown in Figure 10 in the range of 3800–2800 cm-1. It is observed that the sharp O–H stretching band at 3645–3628 cm-1 has affected the spectrum of PCN and created a broad band in the range of 3633–3679 cm-1 and the broad band of ionically bonded N–H stretching at 3256–3035 cm-1 generating a bulge at 3260–3047 cm-1 in PCN associated with aromatic C–H stretching from PBT. FTIR spectra of PCN, PBT, and OMMT are shown in Figure 11 in between 2000 and 1000 cm-1. In this energy-region, bands from PBT are very strong and thus the bands in both PBT and PCN are almost identical. Most bands from OMMT are overlapped by PBT, and Si–O stretching band is not visible in PCN. The PA FTIR spectra of PCN, PBT, and OMMT in the range of 1000–400 cm-1 are shown in Figure 12. The in-plane ring deformation vibration of PBT at 633 cm-1 has shifted to the lower energy side by 9 cm-1 in the PCN. The tripled Al–O, Si–O vibration and in-plane ring deformation vibration shows broader band than that of pure PBT. From Figure 12, it can be seen that the sharp out-of-plane ring deformation vibration of PBT at 497 cm-1 is suppressed and the out-of-plane ring deformation vibration of PBT at 476 cm-1 is absent in PCN. These phenomena indicate that ring deformations in PBT are restricted by OMMT due to intercalation. It is evident that Si–O–Si deformation from OMMT appears in PCN at 463 cm-1 and shifted to the lower energy side by 15 cm-1 and this shift in PCN is probably due to the formation of bonding of PBT with silica tetrahedral oxygen through the hydrogen bond between the functional groups of PBT and 12-aminolauric acid. Due to this effect the Si–O stretching band in PCN disappears. Wang et al. [36] found that during melt processing when PBT/PET blends are filled with silica nanoparticles, hydroxyl end groups of PET and PBT react with the surface hydroxyl groups of SiO2 which is shown by X-ray photoelectron Spectroscopy (XPS) tests. Fermeglia et al. [48] investigated the binding energy between poly(Butylene terephtalate-co-thiodiethylene terephtalate), clay and different organic modifiers using molecular simulations and have found that –NH2, –COOH, or –SH groups in organic modifiers generally lead to a good interaction between modifiers and polymers.Table 3
Band assignments of photoacoustic FTIR spectra of OMMT, PBT, and PCN.
Sample
Wave number range (cm-1)
Bands
Reference
3646–3628
O–H stretching of structural hydroxyl group
[26]
3256–3035
Ionic bonded N–H stretching
2932
C–H asymmetric stretching
2853
C–H symmetric stretching
1713–1701
C=O stretching
1623
Combination of O–H deformation and N–H bending
1511
Asymmetric R–COO- stretching
1466
CO–H bending
1396
Symmetric R–COO- stretching
1112
Si–O in-plane stretching
1060
Si–O out-of-plane stretching
OMMT
1022,997
Si–O in-plane stretching
963
Si–OH vibrations
921
Al–O/Al–OH stretching vibration
884
Al–FeOH deformation
847
Al–MgOH deformation
799
Si–O stretching of quartz and silica
727
Si–O deformation perpendicular to optical axis
699
Si–O deformation parallel to optical axis
623
Coupled Al–O and Si–O out-of-plane vibration
531
Al–OSi deformation
478
Si–OSi vibration
3417
–OH stretching
3052
Aromatic C–H stretching
[35]
2962
C–H2 stretching
[36]
2895
Aliphatic C–H asymmetric stretching
[35]
2863
Aliphatic C–H symmetric stretching
[35]
1722
C=O stretching
[37]
1583
Aromatic C=C stretching
[38]
1504
Aliphatic C–H2 bending
[37]
1453
Aliphatic C–H2 bending
[39]
1407
Aliphatic aromatic C–H in-of-plane bending
[36]
PBT
1387
CH2 wagging
[39]
1281
(C=O–)–O stretching
[35]
1110
CH2–O stretching
[35]
1017
Aromatic C–H bending
[35]
873
Aromatic C–H bending
[40]
808
Aromatic C–H bending
[37]
729
Aromatic C–H out-of-plane bending
[35, 36]
685
Aromatic =C–H out-of-plane deformation
[38]
633
Aromatic in-plane deformation
[38]
497
Aromatic out-of-plane deformation
[38]
476
Aromatic out-of-plane deformation
[38]
3633–3679
O–H stretching of  structural hydroxyl group-OMMT
[26]
3413
–OH stretching
3260–3047
Coupled N–H stretching and aromatic C–H stretching
[35]
2963
C–H2 stretching
[36]
2895
Aliphatic C–H asymmetric stretching
[35]
2858
Aliphatic C–H symmetric stretching
[35]
1719
C=O stretching
[37]
1576
Aromatic C=C stretching
[38]
1503
Aliphatic C–H2 bending
[37]
1452
Aliphatic C–H2 bending
[39]
1406
Aromatic C–H in-of-plane bending
[36]
PCN
1387
CH2 wagging
[39]
1278
(C=O–)–O stretching
[35]
1105
CH2–O stretching
[35]
1020
Aromatic C–H bending
[35]
873
Aromatic C–H bending
[40]
810
Aromatic C–H bending
[37]
729
Aromatic C–H out-of-plane bending
[35, 36]
685
Aromatic =C–H out-of-plane deformation
[38]
652
Aromatic in-plane deformation
[38]
624
Tripled Al–O, Si–O out-of-plane and aromatic in-plane deformation
[26]
497
Aromatic out-of-plane deformation
[38]
463
Si–OSi deformation
[26]Figure 9
Photo acoustic FTIR spectra of OMMT, PBT, and PCN in the energy range of 4000–400cm-1 at scan speed of 0.15 cm-1.Figure 10
Photo acoustic FTIR spectra of OMMT, PBT, and PCN in the energy range of 3800–2800cm-1 at scan speed of 0.15 cm-1.Figure 11
Photo acoustic FTIR spectra of OMMT, PBT, and PCN in the energy range of 2000–1000cm-1 at scan speed of 0.15 cm-1.Figure 12
Photo acoustic FTIR spectra of OMMT, PBT, and PCN in the energy range of 1000–400cm-1 at scan speed of 0.15 cm-1.


## 3.1.3. Crystallinity of PBT and PCN

DSC tests are conducted on the pure PBT and PCN and results are shown in Table4; it is observed that the pure PBT shows 0.5% higher crystallinity than that of PCN. This observation of reduced crystallinity of PCN with respect to a polymer has been observed earlier in Nylon6-PCN [10]. As shown in Table 4, Nylon6-PCN with 9% of clay loading shows 4.56% reduction in crystallinity.Table 4
DSC results of PBT, Nylon6, and PCNs.
Materials
Enthalpy (J/g)
Crystallinity (%)
PBT
50.21
35.6
PBT-PCN
49.56
35.1
Nylon6
56.06
27.61
Nylon6-PCN
46.81
23.05


## 3.1.4. Nanomechanical Properties of PBT and PCN

In the present work, nanoindentation tests for PBT and PBT-PCN are performed in load control mode. Figure13 is the representative L-D curves for PBT and PBT-PCN films. Maximum load was maintained for period of 5 s between loading and unloading sequence in order to reduce the time-dependent deformation. The elastic modulus and hardness values for PBT and PBT-PCN are presented in Table 5. Nylon6 and Nylon6-PCN nanomechanical tests were performed in load control mode [10] and the results are presented in Table 5.Table 5
Nanomechanical properties of PCNs.
Parameters
PBT
PBT-PCN
Nylon6
Nylon6-PCN
Crystallinity (%)
35.6
35.1
27.6
23.1
Elastic modulus (GPa)
1.40
2.20
3.35
5.46
Hardness (GPa)
0.09
0.14
0.12
0.17Figure 13
Representative load displacement (L-D) curves for PBT and PCN.


## 3.2. Interactions in PCNs System

### 3.2.1. Nonbonded Interactions between Organic Modifier and Clay

Interaction energies between different constituents of OMMT and PCNs are calculated using NAMDEnergy module of NAMD software. Nonbonded interactions, van der waals and electrostatic, are calculated in terms of attractive and repulsive interactions. Among different initiald-spacing of OMMT, d-spacing of 17.60 Å OMMT model predicts final d-spacing of 15.50 Å which is very close to experimental d-spacing of 15.60 Å. Based on this comparison of d-spacing, the OMMT model is chosen for simulation. Among six different PCNs models, the PCN model which contains 6.41 wt% of PBT predicts the simulated d-spacing of 14.80 Å which is close to experimental d-spacing of 14.53 Å. The same OMMT model with Nylon6 is simulated and it predicts d-spacing of 15.35 Å. Non-bonded interactions between organic modifier and clay are calculated and the interaction energies are shown in Table 6. Positive values indicate repulsive interactions and negative values indicate attractive interactions between constituents and higher magnitude denotes greater interactions. Interaction energy values of organic modifiers are the sum of backbone and functional group interaction energies. Among these two components, functional group of organic modifier shows higher interaction with clay than the backbone of the organic modifier. This is because of higher positive net charge in the functional group compared to the backbone. Functional group interactions are composed of −1096 kcal/mol in electrostatic interaction energy and −80 kcal/mol in van der waals interaction energy and backbone interactions are composed of −549 kcal/mol electrostatic interaction energy and −345 kcal/mol van der waals interaction energy. It is observed that the electrostatic interactions are the major components of both modifier functional group and modifier backbone interaction energies with the clay. More detailed quantitative results are presented in Table 7. Here, it is observed that the major component of attractive interaction with clay is from the hydrogen atoms of modifier. Even though the partial charge of the backbone hydrogen smaller is compared to functional group hydrogen as shown in Figure 2(a), backbone hydrogen shows higher interaction because of its presence in significant numbers.Table 6
Nonbonded interaction energies of OMMT constituents.
OMMT constituents
Nonbonded energy (kcal/mol)
Clay-modifier backbone
−894
Clay-modifier functional group
−1176Table 7
Nonbonded energies of individual atoms of organic modifier with clay.
Component of modifier and clay
Electrostatic energy (kcal/mol) [Col-A]
Van der waals energy (kcal/mol) [Col-B ]
Total nonbonded energy (kcal/mol) [Col-A + Col-B]
Clay-modifier backbone hydrogen
−2503
−87
−2590
Clay-modifier backbone carbon
+1954
−257
+1697
Clay-modifier functional hydrogen
−2392
−8
−2400
Clay-modifier functional nitrogen
+497
−0.5
+497
Clay-modifier functional oxygen
+1689
−48
+1641
Clay-modifier functional carbon
−890
−25
−915


### 3.2.2. Nonbonded Interactions of PBT-PCN and Nylon6-PCN System

The nonbonded interaction energies between constituents of PBT-PCN and Nylon6-PCN are presented in Table8. Interaction energy between the clay and modifier is attractive and higher than interaction energy between clay and polymer. This is because of higher partial charges of individual atoms in modifiers than the polymers. Interaction energies between clay and polymer are attractive and higher than the interaction energy between polymer and modifier. Detailed results are shown in Table 9 and it is observed that the interaction energy between clay and modifier functional group is attractive and higher than among other constituents in both PCNs. This higher interaction between clay and modifier functional groups could be attributed to higher positive partial charges in functional groups atoms as shown in Figure 2(a). Similarly, modifier backbone also has higher attractive interactions with clay compared to the polymer functional group and backbone in both PCNs.Table 8
Nonbonded energies PCNs.
Component of PCN
Electrostatic energy (kcal/mol)
Van der waals energy (kcal/mol)
Total nonbonded energy (kcal/mol)
[Col-A]
[Col-B ]
[Col-A + Col-B]
PBT-PCN
Clay-polymer
−36
−165
−201
Clay-modifier
−1836
−316
−2152
Polymer-modifier
−93
−43
−136
Nylon6-PCN
Clay-polymer
12
−182
−170
Clay-modifier
−1464
−375
−1839
Polymer-modifier
−209
−30
−239Table 9
Nonbonded energies of modifier and polymer with clay.
Component of modifier, polymer, and clay in PCN
Electrostatic energy (kcal/mol)
Van der waals energy (kcal/mol)
Total nonbonded energy (kcal/mol)
[Col-A]
[Col-B ]
[Col-A + Col-B]
PBT-PCN
Clay-modifier backbone
−582
−282
−864
Clay-modifier functional group
−1255
−34
−1289
Clay-polymer backbone
−250
−108
−358
Clay-polymer functional group
+213
−56
+157
Nylon6-PCN
Clay-modifier backbone
−485
−297
−782
Clay-modifier functional group
−979
−77
−1056
Clay-polymer backbone
−183
−123
−306
Clay-polymer functional group
+155
−62
+93Polymer functional group has negative partial charges as shown in Figure2(b) thus it induces repulsive interaction with clay in both PCNs systems. Interaction energies between clay and individual atoms of organic modifier are shown in Table 10. It is observed that the organic modifier hydrogen atoms show higher magnitude of attractive interactions with clay as observed in the OMMT system. Interaction energies of individual polymer atoms with clay are shown in Table 11. Here, hydrogen and carbon atoms show attractive interaction with clay while the other atoms show repulsive interactions. Functional group carbon atoms contribute towards more attractive interactions and greater negative partial charge of oxygen atoms contribute to more repulsive interactions. As shown in Table 12, polymer backbone has higher repulsive interaction with modifier functional group and lesser repulsive interaction with modifier backbone. Polymer functional group has higher attractive interaction with modifier functional group and lesser attractive interaction with backbone. Similar nature of interactions have been reported in our previous work [11].Table 10
Nonbonded energies between individual atoms of organic modifier and clay.
Individual atoms of organic modifiers and clay
Electrostatic energy (kcal/mol)
Van der waals energy (kcal/mol)
Total nonbonded energy (kcal/mol)
[Col-A]
[Col-B ]
[Col-A + Col-B]
PBT-PCN
Clay-modifier backbone hydrogen
−2512
−67
−2579
Clay-modifier backbone carbon
+1930
−215
+1715
Clay-modifier functional hydrogen
−2564
−7
−2571
Clay-modifier functional nitrogen
+546
+37
+583
Clay-modifier functional oxygen
+1628
−42
+1586
Clay-modifier functional carbon
−865
−22
−887
Nylon6-PCN
Clay-modifier backbone hydrogen
−2085
−72
−2157
Clay-modifier backbone carbon
+1600
−226
+1374
Clay-modifier functional hydrogen
−2071
−7
−2078
Clay-modifier functional nitrogen
+444
−2
+442
Clay-modifier functional oxygen
+1359
−44
+1315
Clay-modifier functional carbon
−711
−24
−735Table 11
Nonbonded energies between individual atoms of polymer and clay.
Component of polymer and clay in PCN
Electrostatic energy (kcal/mol) [Col-A]
Van der waals energy (kcal/mol) [Col-B ]
Total nonbonded energy (kcal/mol) [Col-A + Col-B]
PBT-PCN
Clay-polymer aromatic  hydrogen
−301
−10
−311
Clay-polymer backbone hydrogen
−380
−10
−390
Clay-polymer aromatic  carbon
+188
−59
+129
Clay-polymer backbone carbon
+243
−29
+214
Clay-polymer functional hydrogen
−90
−0.2
−90.2
Clay-polymer functional oxygen
+1127
−36
1091
Clay-polymer functional carbon
−824
−20
−844
Nylon6-PCN
Clay-polymer backbone hydrogen
−1056
−19
−1075
Clay-polymer backbone carbon
+872
−104
+768
Clay-polymer functional nitrogen
+593
−27
+566
Clay-polymer functional hydrogen
−424
−2
−426
Clay-polymer functional oxygen
+589
−13
+576
Clay-polymer functional carbon
−600
−20
−620
Clay-polymer backbone hydrogen
−1056
−19
−1075Table 12
Nonbonded energies between polymer and organic modifier.
Component of modifier and polymer in PCN
Electrostatic energy (kcal/mol)[Col-A]
Van der waals energy (kcal/mol)[Col-B ]
Total nonbonded energy (kcal/mol)       [Col-A + Col-B]
PBT-PCN
Polymer backbone—modifier backbone
+115
−29
+86
Polymer backbone—modifier functional
+209
−10
+199
Polymer functional—modifier backbone
−132
−12
−144
Polymer functional—modifier functional
−285
+7
−278
Nylon6-PCN
Polymer backbone—modifier backbone
+58
−20
+38
Polymer backbone—modifier functional
+91
−11
+80
Polymer functional—modifier backbone
−94
−11
−105
Polymer functional—modifier functional
−264
+11
−253


### 3.2.3. Summary of Interaction of OMMT, PBT-PCN, and Nylon6-PCN System

As mentioned in Section1, nonbonded interactions play a major role in PCNs and they affect the crystallinity and mechanical properties of PCNs. Based on our results, more than 80%, 75%, 70% of interaction energies are a result of electrostatic interaction for OMMT, PBT-PCN and Nylon6-PCN respectively. Figures 14 and 15 show overall interaction energy map for PBT-PCN, and Nylon6-PCN, respectively. The two energy maps have similar pattern in terms of attractive and repulsive interactions between different constituents of PCNs. These energy maps indicate that the PBT-PCN system has greater magnitude of interaction energies between different constituents than the Nylon6-PCN system. In all three systems, OMMT, PBT-PCN, and Nylon6-PCN, modifier shows attractive interactions with the clay. Further, in PBT-PCN system, modifiers have 17% higher attractive interactions with clay than in Nylon6-PCN system. Polymer has attractive interaction with clay and modifier in both systems. In terms of interaction energy magnitudes, almost equal interaction of polymer with clay is observed in both systems and 75 % higher interaction with organic modifier is observed in Nylon6-PCN compared to PBT-PCN.Figure 14
Energy map of PBT-PCN system.Figure 15
Energy map of nylon6-PCN system.


## 3.2.1. Nonbonded Interactions between Organic Modifier and Clay

Interaction energies between different constituents of OMMT and PCNs are calculated using NAMDEnergy module of NAMD software. Nonbonded interactions, van der waals and electrostatic, are calculated in terms of attractive and repulsive interactions. Among different initiald-spacing of OMMT, d-spacing of 17.60 Å OMMT model predicts final d-spacing of 15.50 Å which is very close to experimental d-spacing of 15.60 Å. Based on this comparison of d-spacing, the OMMT model is chosen for simulation. Among six different PCNs models, the PCN model which contains 6.41 wt% of PBT predicts the simulated d-spacing of 14.80 Å which is close to experimental d-spacing of 14.53 Å. The same OMMT model with Nylon6 is simulated and it predicts d-spacing of 15.35 Å. Non-bonded interactions between organic modifier and clay are calculated and the interaction energies are shown in Table 6. Positive values indicate repulsive interactions and negative values indicate attractive interactions between constituents and higher magnitude denotes greater interactions. Interaction energy values of organic modifiers are the sum of backbone and functional group interaction energies. Among these two components, functional group of organic modifier shows higher interaction with clay than the backbone of the organic modifier. This is because of higher positive net charge in the functional group compared to the backbone. Functional group interactions are composed of −1096 kcal/mol in electrostatic interaction energy and −80 kcal/mol in van der waals interaction energy and backbone interactions are composed of −549 kcal/mol electrostatic interaction energy and −345 kcal/mol van der waals interaction energy. It is observed that the electrostatic interactions are the major components of both modifier functional group and modifier backbone interaction energies with the clay. More detailed quantitative results are presented in Table 7. Here, it is observed that the major component of attractive interaction with clay is from the hydrogen atoms of modifier. Even though the partial charge of the backbone hydrogen smaller is compared to functional group hydrogen as shown in Figure 2(a), backbone hydrogen shows higher interaction because of its presence in significant numbers.Table 6
Nonbonded interaction energies of OMMT constituents.
OMMT constituents
Nonbonded energy (kcal/mol)
Clay-modifier backbone
−894
Clay-modifier functional group
−1176Table 7
Nonbonded energies of individual atoms of organic modifier with clay.
Component of modifier and clay
Electrostatic energy (kcal/mol) [Col-A]
Van der waals energy (kcal/mol) [Col-B ]
Total nonbonded energy (kcal/mol) [Col-A + Col-B]
Clay-modifier backbone hydrogen
−2503
−87
−2590
Clay-modifier backbone carbon
+1954
−257
+1697
Clay-modifier functional hydrogen
−2392
−8
−2400
Clay-modifier functional nitrogen
+497
−0.5
+497
Clay-modifier functional oxygen
+1689
−48
+1641
Clay-modifier functional carbon
−890
−25
−915


## 3.2.2. Nonbonded Interactions of PBT-PCN and Nylon6-PCN System

The nonbonded interaction energies between constituents of PBT-PCN and Nylon6-PCN are presented in Table8. Interaction energy between the clay and modifier is attractive and higher than interaction energy between clay and polymer. This is because of higher partial charges of individual atoms in modifiers than the polymers. Interaction energies between clay and polymer are attractive and higher than the interaction energy between polymer and modifier. Detailed results are shown in Table 9 and it is observed that the interaction energy between clay and modifier functional group is attractive and higher than among other constituents in both PCNs. This higher interaction between clay and modifier functional groups could be attributed to higher positive partial charges in functional groups atoms as shown in Figure 2(a). Similarly, modifier backbone also has higher attractive interactions with clay compared to the polymer functional group and backbone in both PCNs.Table 8
Nonbonded energies PCNs.
Component of PCN
Electrostatic energy (kcal/mol)
Van der waals energy (kcal/mol)
Total nonbonded energy (kcal/mol)
[Col-A]
[Col-B ]
[Col-A + Col-B]
PBT-PCN
Clay-polymer
−36
−165
−201
Clay-modifier
−1836
−316
−2152
Polymer-modifier
−93
−43
−136
Nylon6-PCN
Clay-polymer
12
−182
−170
Clay-modifier
−1464
−375
−1839
Polymer-modifier
−209
−30
−239Table 9
Nonbonded energies of modifier and polymer with clay.
Component of modifier, polymer, and clay in PCN
Electrostatic energy (kcal/mol)
Van der waals energy (kcal/mol)
Total nonbonded energy (kcal/mol)
[Col-A]
[Col-B ]
[Col-A + Col-B]
PBT-PCN
Clay-modifier backbone
−582
−282
−864
Clay-modifier functional group
−1255
−34
−1289
Clay-polymer backbone
−250
−108
−358
Clay-polymer functional group
+213
−56
+157
Nylon6-PCN
Clay-modifier backbone
−485
−297
−782
Clay-modifier functional group
−979
−77
−1056
Clay-polymer backbone
−183
−123
−306
Clay-polymer functional group
+155
−62
+93Polymer functional group has negative partial charges as shown in Figure2(b) thus it induces repulsive interaction with clay in both PCNs systems. Interaction energies between clay and individual atoms of organic modifier are shown in Table 10. It is observed that the organic modifier hydrogen atoms show higher magnitude of attractive interactions with clay as observed in the OMMT system. Interaction energies of individual polymer atoms with clay are shown in Table 11. Here, hydrogen and carbon atoms show attractive interaction with clay while the other atoms show repulsive interactions. Functional group carbon atoms contribute towards more attractive interactions and greater negative partial charge of oxygen atoms contribute to more repulsive interactions. As shown in Table 12, polymer backbone has higher repulsive interaction with modifier functional group and lesser repulsive interaction with modifier backbone. Polymer functional group has higher attractive interaction with modifier functional group and lesser attractive interaction with backbone. Similar nature of interactions have been reported in our previous work [11].Table 10
Nonbonded energies between individual atoms of organic modifier and clay.
Individual atoms of organic modifiers and clay
Electrostatic energy (kcal/mol)
Van der waals energy (kcal/mol)
Total nonbonded energy (kcal/mol)
[Col-A]
[Col-B ]
[Col-A + Col-B]
PBT-PCN
Clay-modifier backbone hydrogen
−2512
−67
−2579
Clay-modifier backbone carbon
+1930
−215
+1715
Clay-modifier functional hydrogen
−2564
−7
−2571
Clay-modifier functional nitrogen
+546
+37
+583
Clay-modifier functional oxygen
+1628
−42
+1586
Clay-modifier functional carbon
−865
−22
−887
Nylon6-PCN
Clay-modifier backbone hydrogen
−2085
−72
−2157
Clay-modifier backbone carbon
+1600
−226
+1374
Clay-modifier functional hydrogen
−2071
−7
−2078
Clay-modifier functional nitrogen
+444
−2
+442
Clay-modifier functional oxygen
+1359
−44
+1315
Clay-modifier functional carbon
−711
−24
−735Table 11
Nonbonded energies between individual atoms of polymer and clay.
Component of polymer and clay in PCN
Electrostatic energy (kcal/mol) [Col-A]
Van der waals energy (kcal/mol) [Col-B ]
Total nonbonded energy (kcal/mol) [Col-A + Col-B]
PBT-PCN
Clay-polymer aromatic  hydrogen
−301
−10
−311
Clay-polymer backbone hydrogen
−380
−10
−390
Clay-polymer aromatic  carbon
+188
−59
+129
Clay-polymer backbone carbon
+243
−29
+214
Clay-polymer functional hydrogen
−90
−0.2
−90.2
Clay-polymer functional oxygen
+1127
−36
1091
Clay-polymer functional carbon
−824
−20
−844
Nylon6-PCN
Clay-polymer backbone hydrogen
−1056
−19
−1075
Clay-polymer backbone carbon
+872
−104
+768
Clay-polymer functional nitrogen
+593
−27
+566
Clay-polymer functional hydrogen
−424
−2
−426
Clay-polymer functional oxygen
+589
−13
+576
Clay-polymer functional carbon
−600
−20
−620
Clay-polymer backbone hydrogen
−1056
−19
−1075Table 12
Nonbonded energies between polymer and organic modifier.
Component of modifier and polymer in PCN
Electrostatic energy (kcal/mol)[Col-A]
Van der waals energy (kcal/mol)[Col-B ]
Total nonbonded energy (kcal/mol)       [Col-A + Col-B]
PBT-PCN
Polymer backbone—modifier backbone
+115
−29
+86
Polymer backbone—modifier functional
+209
−10
+199
Polymer functional—modifier backbone
−132
−12
−144
Polymer functional—modifier functional
−285
+7
−278
Nylon6-PCN
Polymer backbone—modifier backbone
+58
−20
+38
Polymer backbone—modifier functional
+91
−11
+80
Polymer functional—modifier backbone
−94
−11
−105
Polymer functional—modifier functional
−264
+11
−253


## 3.2.3. Summary of Interaction of OMMT, PBT-PCN, and Nylon6-PCN System

As mentioned in Section1, nonbonded interactions play a major role in PCNs and they affect the crystallinity and mechanical properties of PCNs. Based on our results, more than 80%, 75%, 70% of interaction energies are a result of electrostatic interaction for OMMT, PBT-PCN and Nylon6-PCN respectively. Figures 14 and 15 show overall interaction energy map for PBT-PCN, and Nylon6-PCN, respectively. The two energy maps have similar pattern in terms of attractive and repulsive interactions between different constituents of PCNs. These energy maps indicate that the PBT-PCN system has greater magnitude of interaction energies between different constituents than the Nylon6-PCN system. In all three systems, OMMT, PBT-PCN, and Nylon6-PCN, modifier shows attractive interactions with the clay. Further, in PBT-PCN system, modifiers have 17% higher attractive interactions with clay than in Nylon6-PCN system. Polymer has attractive interaction with clay and modifier in both systems. In terms of interaction energy magnitudes, almost equal interaction of polymer with clay is observed in both systems and 75 % higher interaction with organic modifier is observed in Nylon6-PCN compared to PBT-PCN.Figure 14
Energy map of PBT-PCN system.Figure 15
Energy map of nylon6-PCN system.


## 3.3. Discussion

Our previous work has shown that the molecular interactions between clay, modifier, and polymer results in significant alteration of polymer to a very large extent in PCN [10]. This is the basis of the altered phase theory to explain the significant improvement in mechanical properties of PCN as a result of addition of small amounts of well-dispersed clay nanoparticles in the polymer. The attractive interactions between functional groups of the polymer with the modifier and concurrent repulsive interaction between modifier and polymer backbone result in reduction of crystallinity of the surrounding polymer.In the current work with PBT-PCN system we have observed reduction in crystallinity and improved mechanical properties of the PCN with respect to the pristine polymer. We have seen similar phenomena of reduction in crystallinity and improved mechanical properties of Nylon6-PCN [10]; however, the magnitudes of reduction in crystallinity for the two PCNs are different and the improvement in mechanical properties is different. Based on the interaction energy maps for the two systems shown in Figures 14 and 15, the mechanism of attractive interaction between functional group of the polymer and modifier and repulsive interaction between backbone of the polymer and the modifier leading to changes in polymer crystallinity is present, although the magnitudes of the interaction energies are different. In the case of Nylon6-PCN the reduction in crystallinity is 16.5 % while the reduction in crystallinity is 0.5% in PBT-PCN. The improvement in elastic modulus and hardness values for Nylon6-PCN is 63% and 42%, respectively, with respect to pristine Nylon6 while the improvement in elastic modulus and hardness values for PBT-PCN is only 10% and 23%, respectively, with respect to pristine PBT. It appears that in the case of Nylon6-PCN, the clay-modifier system has significantly altered crystallinity and improved the mechanical properties, while the same clay-modifier system has impacted only moderately the crystallinity and mechanical properties of the PBT-PCN system.The comparison of interaction energy maps for the two PCNs shown in Figures14 and  15 reveals that the interaction energies between polymer and modifier and polymer and clay are higher in the PBT-PCN system than the Nylon6-PCN system. The attractive interaction energy and repulsive interaction energy between polymer and modifier are 422 kcal/mol and 285 kcal/mol for PBT-PCN while the attractive and repulsive interaction energies for Nylon6-PCN are 358 kcal/mol and 118 kcal/mol, respectively. The attractive and repulsive interaction energies between clay and polymer are 157 kcal/mol and 358 kcal/mol for PBT-PCN and 93 kcal/mol and 306 kcal/mol for Nylon6-PCN. Inspite of the higher interaction energies between polymer and modifier in PBT-PCN system than the Nylon6-PCN system, the alteration of crystallinity and improvement in mechanical properties is lower for PBT-PCN system. It appears that the initial crystallinity of the polymer may have affected the magnitude of interaction energy needed to alter the crystallinity and improve the mechanical properties of the PCN. The crystallinity of pure PBT is 35.6% while the crystallinity of Nylon6 is 27.6%. Thus, it appears that significantly higher interaction energies between clay-modifier and polymer may be required to alter the crystallinity and improve mechanical properties of a polymer that is more crystalline to begin with. Modest differences are observed in the attractive interaction energies between the same modifier and clay in the presence of different polymers.


## 4. Conclusions

A comprehensive assessment of PBT-PCN and Nylon6-PCN systems using molecular dynamics and experimental techniques has been conducted. Experimental results and interaction energy maps of PBT-PCN system indicate that the underlying mechanisms of change in crystallinity and improvement in mechanical properties as proposed by altered phase theory are valid. Although the same amount of modified clay is used in the preparation of PCN with PBT and Nylon6, a dramatic difference in change in crystallinity and improvement in mechanical properties is observed. It appears that a polymer with higher crystallinity could require significantly higher attractive and repulsive interaction energies between the polymer and modifier compared to a polymer with lower crystallinity to achieve similar magnitude of percent change in crystallinity and improvement in mechanical properties of PCN.



---

*Source: 341056-2012-06-20.xml*