Next Article in Journal
Two-Step Dopamine-to-Polydopamine Modification of Polyethersulfone Ultrafiltration Membrane for Enhancing Anti-Fouling and Ultraviolet Resistant Properties
Previous Article in Journal
Performance Evaluation of Red Clay Binder with Epoxy Emulsion for Autonomous Rammed Earth Construction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High-Barrier Polyimide Containing Carbazole Moiety: Synthesis, Gas Barrier Properties, and Molecular Simulations

1
National and Local Joint Engineering Center of Advanced Packaging Materials R & D Technology, Key Laboratory of Advanced Packaging Materials and Technology of Hunan Province, School of Packaging and Materials Engineering, Hunan University of Technology, Zhuzhou 412007, China
2
Key Laboratory of Polymeric Materials and Application Technology of Hunan Province, Key Laboratory of Advanced Functional Polymer Materials of Colleges, Universities of Hunan Province, College of Chemistry, Xiangtan University, Xiangtan 411105, China
*
Author to whom correspondence should be addressed.
Polymers 2020, 12(9), 2048; https://doi.org/10.3390/polym12092048
Submission received: 2 August 2020 / Revised: 22 August 2020 / Accepted: 28 August 2020 / Published: 8 September 2020
(This article belongs to the Section Polymer Physics and Theory)

Abstract

:
A high-barrier polyimide (2,7-CPI) was synthesized through the polymerization of pyromellitic dianhydride (PMDA) and a novel diamine (2,7-CDA) containing carbazole moiety. The synthesized diamine and polyimide were fully characterized by elemental analyses, FTIR and NMR. The 2,7-CPI displays very attractive barrier performances, with oxygen transmission rate (OTR) and water vapor transmission rate (WVTR) low to 0.14 cm3·m−2·day−1 and 0.05 g·m−2·day−1, respectively. Meanwhile, 2,7-CPI also exhibits exceptional thermal stability with a glass transition temperature (Tg) of 467 °C, 5% weight-loss temperature (Td5%) of 550 °C under N2 and coefficient of thermal expansion (CTE) of 3.4 ppm/K. The barrier performances of 2,7-CPI are compared with those of a structural analogue (2,7-CPPI) and a typical polyimide (Kapton). Their barrier performances with respect to microstructure were investigated by molecular simulations, wide angle X-ray diffraction (WAXD), and positron annihilation lifetime spectroscopy (PALS). The results show that 2,7-CPI possesses better coplanar structure and more number of intermolecular hydrogen bonds among the three PIs, which result in tight chain packing and thereby high crystallinity, low free volume, and decreased chains mobility. That is, the high crystallinity and low free volume of 2,7-CPI reduce the diffusion and solubility of gases. Meanwhile, the poor chains mobility further decreases the gases diffusion. The reduced diffusion and solubility of gases consequently promote the improvement of barrier properties for 2,7-CPI. The polyimide has a wide application prospect in the field of flexible electronic packaging industries.

Graphical Abstract

1. Introduction

In recent years, flexible displays have attracted considerable attention because of the low weight, high mobility, and their portable, wearable, and even foldable applications [1,2]. Flexible active matrix organic light emitting diode (AMOLED) display is widely accepted as the most promising development trend for flexible displays technology, because they show low power consumption, faster response, and are more power efficient [3,4]. Nowadays, plastic films are widely used as substrates for flexible AMOLED, owing to their light weight, good flexibility, and excellent robustness [5,6]. With the wide use of top-emitting architectures in flexible AMOLED, transparency is not a prerequisite for the substrates [7]. However, to realize the application of plastic substrates in AMOLED, two key requirements should first be satisfied for the plastic films. One is the high gas barrier properties (WVTR < 10−6 g·m−2·day−1 and OTR < 10−5 cm3·m−2·day−1), which can prevent the degradation of organic materials in the devices caused by the penetration of moisture and oxygen through the substrates [8,9]. The other is the high thermal and dimensional stability. As the driven device, high-performance thin-film transistors (TFTs) are crucial for AMOLED displays. High processing temperature is significantly important for the performance improvement of TFTs [10]. The widely used low-temperature-poly-silicon (LTPS) TFTs require the plastic substrates to withstand process temperature of higher than 400 °C [11]. In addition, high dimensional stability can prevent the substrates to shrink during the temperature cycles associated with the TFT fabrication, and guarantee good overlay alignment in the TFT array [12]. Hence, the plastic substrates are required to possess excellent thermal and dimensional stability (Tg ≥ 400 °C, CTE ≤ 8 ppm/K) [9,11]. However, the demands for thermal and dimensional stability far exceed the tolerance range of conventional plastic substrates such as polyethylene terephthalate (PET) and polyethylene naphthalate (PEN) (Tg ˂ 200 °C, CTE of 20~80 ppm/K) [13].
In this regard, polyimides (PIs) have been considered as one of the best choices for flexible AMOLED substrates because of their excellent thermal stability and mechanical properties [14,15]. However, the insufficient barrier performances of PIs affect their application in AMOLED substrates. To improve the barrier performances of PIs, inorganic barrier layers such as aluminum oxide and silicon dioxide or inorganic/organic alternating layers are usually fabricated on PIs by sputtering or chemical vapor deposition (CVD) [16,17]. However, the deposition methods are usually complicated, and require large-scale equipment, thus lead to fairly high fabrication costs. Apparently, these deposition processes are on the basis of intrinsic PI matrix. Hence, enhancing the barrier performances of intrinsic PIs can reduce the thickness of surface barrier layers [18], simplify the deposition processes and decrease costs.
Many researches reported that the gas permeability of polymers was greatly determined by their chemical structures, aggregation structures, and interchain cohesion [19,20]. Polyimides possess the advantage of high designability of molecular structures. In this case, structural modifications are considered as effective routes to enhance the barrier performances of PIs. In our previous work, a diamine (2,7-CPDA, see Figure 1) containing rigid planar carbazole moiety was synthesized and polymerized with PMDA to obtain polyimide (2,7-CPPI) [21]. The introduction of carbazole moieties increases the regularity and rigidity of polymer chains and the interchain cohesion, which improves the chains packing efficiency and thereby the barrier performances of 2,7-CPPI.
In this work, a novel diamine (2,7-CDA) containing rigid planar carbazole moiety was designed and synthesized. Then, the 2,7-CDA was polymerized with PMDA to form polyimide (2,7-CPI). The chemical structure of 2,7-CDA is shown in Figure 1. Compared with 2,7-CPDA, the two amine groups in 2,7-CDA are directly attached to the 2,7-positions of the carbazole ring without the bridging of benzene rings. The stereo view of the lowest energy conformations of repeat units for 2,7-CPI and 2,7-CPPI are presented in Figure 1. It can be observed that the diamine residues in 2,7-CPI exhibit a better coplanarity than those in 2,7-CPPI, which is conducive to the close chains packing. In addition, the direct connection of two large rigid planar structures (carbazole ring and dianhydride residue) in 2,7-CPI will increase the chains rigidity and reduce the mobility of polymer chains. These above-mentioned structural characteristics of 2,7-CPI are beneficial to the improvement of barrier performances.
Here, the synthesis, characterizations, as well as barrier and thermal properties of 2,7-CPI were reported. Additionally, the barrier performances of 2,7-CPI are compared with those of a typical polyimide (Kapton) based on diamine 4,4′-oxydianiline (ODA) and PMDA and our previous reported structural analogue 2,7-CPPI. Molecular simulation is an increasingly vital tool to provide detailed insights into the structural characteristics and gas transport performances of polymers at a microscopic level, which is hard to obtain by experiments [22,23]. In order to elucidate the effect of different diamine structures on the barrier properties of PIs, molecular simulations, WAXD and PALS were utilized to study the chains packing, intermolecular hydrogen bonds, free volume, chains mobility, penetrant trajectories, diffusivities, and sorption behavior of PIs. This will be helpful in revealing the barrier mechanism and developing high-barrier polymer membranes.

2. Experimental Section

2.1. Materials

O-aminobiphenyl, decalin, palladium 10% on activated carbon (10% Pd/C) were purchased from Energy-chemical Company (Shanghai, China). Hydrazine hydrate and sodium nitrite (NaNO2) were obtained from Aladdin Company (Shanghai, China). Sodium azide and concentrated sulfuric acid were purchased from Wuhan changcheng Chemical Technology Co., Ltd. (Wuhan, China). Potassium nitrate (KNO3), sodium hydroxide (NaOH), sodium chloride (NaCl), absolute ethanol (EtOH), N-hexane, methylene chloride, ethyl acetate, N,N-dimethylformamide (DMF), and silica gel were obtained from Sinopharm Group Chemical Reagent Co., Ltd. (Shanghai, China). ODA and PMDA were purchased from Alfa-Aesar company (Shanghai, China). Before use, the PMDA was dried at 110 °C under vacuum for 6 h and DMF was purified by distillation under inert nitrogen atmosphere.

2.2. Instrumentation

All nuclear magnetic resonance spectra (NMR) were recorded on a Bruker ARX400 MHz spectrometer (Bruker Corporation, Fallanden, Switzerland). Samples were prepared as solution of 5~10 mg of compound in 0.5 mL of deuterated dimethyl sulfoxide (DMSO) using tetramethylsilane (TMS) as the internal reference. Mass spectra were measured on a AcquityUPLC/UPC2/Xevo G2-XS QTOFMS (Waters Corporation, Milford, MA, USA). Elemental analysis was carried out on a Vario EL cube Elemental Analyzer (Elementar Corporation, Langenselbold, Germany). Infrared spectra were recorded on a Nicolet iS10 Fourier-transform infrared (FT-IR) spectrometer (Thermo Fisher Scientific, Waltham, MA, USA). The molecular mass of the poly(amic acid) (PAA) was estimated by gel permeation chromatography on multi-angle laser light scattering (GPC-MALLS) system (Wyatt Technology Corporation, Santa Barbara, CA, USA). DMF is the eluent with flow rate of 1 mL/min and test temperature of 50 °C. Wide angle X-ray diffractograms (WAXD) were recorded by a Rigaku, Ultema III X-ray diffractometer (Rigaku Corporation, Tokyo, Japan) using a Cu Kα radiation. Density was obtained with an ALFA MIRAGE SD-200L electronic density balance (Alfa Mirage Corporation, Osaka, Japan).
Thermogravimetric analyses (TGA) were performed with a TA thermal analyzer (TGA55, TA instruments, New Castle, DE, USA) under N2 with a heating rate of 20 °C/min from 50 to 800 °C, and heated under flowing nitrogen (40 mL/min). The dynamic mechanical (DMA) spectra of the samples were obtained by using TA thermal analyzer (DMA Q850, TA instruments, New Castle, DE, USA). The specimens were analyzed in tensile mode at a constant frequency of 1 Hz, amplitude of 20 μm, and a temperature range from 25 to 500 °C at a heating rate of 5 °C/min. Thermal mechanical analysis (TMA) was used to study the coefficient of thermal expansion (CTE) of the film with a heating rate of 5 °C/min from 25 °C to 350 °C by TMA Q400 instrument (TA instruments, New Castle, DE, USA) under nitrogen.
Tensile test was performed on samples cut from 35~50 μm thick sheet and tested using SANS CMT6103 instrument (New Sansi Material Testing Co., Ltd. Shenzhen, China) according to GB/T16421-1996. The specimen size is 10 × 100 mm, Jaw separation is 50 mm. Jaw speed was first set to 2 mm/min, when elongation reached 1 mm, the Jaw speed was changed to 20 mm/min.
The oxygen permeability of a PI film was measured by Mocon (Mocon Corporation, Minneapolis, MN, USA), in accordance with ASTM-D 3985, using an OX-TRAN 2/21 ML instrument at 23 °C and 0% RH. The moisture vapor permeability was measured using the model PERMATRAN-W® 3/33 of the Mocon Corporation (Indio, CA, USA) at 90% RH and 37.8 °C according to ASTM F-1249. The specimen of 75 μm was fixed on an aluminum foil with the film testing area of 5 cm2.
Positron lifetime measurements were performed as follows: two identical samples with dimension of 1.5 × 10 × 10 mm3 were sandwiched with a 10 μCi 22Na positron source. The 22Na nucleus emits a 1.28 MeV γ-ray simultaneously (within a few ps) with the positron. The positron lifetime is determined from the time delay between the emission of the birth gamma (1.28 MeV) and one of the 0.511 MeV annihilation photons. Lifetime measurements were carried out using a fast-fast coincidence system with a time resolution of about 210 ps and a channel width of 12.6 ps. We have analyzed the lifetime spectra with the help of the data processing programs PATFIT. Before analyzing each spectrum, the positron source components (377 ps/11.75%, 1.04 ns/0.17%) were subtracted. The variance of the fits was around 1.1.

2.3. Synthesis of 4,4′-Dinitro-[1,1′-Biphenyl]-2-Amine (DPNA)

O-aminobiphenyl (4.231 g, 25 mmol) and concentrated sulfuric acid (50 mL) were added to a 100-mL three-necked flask under argon atmosphere. Then, KNO3 (5.303 g, 52.5 mmol) was added with stirring at 0 °C for 4 h. After that, the reaction mixture was poured into ice water. NaOH was added to adjust the solution to neutral. Subsequently, the orange intermediate (DPNA) was obtained through filtration, which was vacuum-dried at 80 °C for 24 h. The yield of the product was about 75 % (4.860 g). IR (KBr, v, cm−1): 3380 (-NH2 stretching), 1507 (−NO2 stretching), 1260 (C−N stretching), 1100~700 (Ar–H stretching). 1H NMR (400 MHz, DMSO-d6, δ, ppm): 8.31 (t, J = 12.8 Hz, 2H), 7.85 – 7.61 (m, 3H), 7.45 (d, J = 7.9 Hz, 1H), 7.28 (d, J = 8.1 Hz, 1H), 5.75 (s, 2H). 13C NMR (100 MHz, DMSO-d6, δ, ppm): 148.79, 147.29, 147.22, 145.18, 131.76, 130.48, 129.60, 124.49, 110.95, 109.66. MS (EI, m/z): 259(100) ([M+], calcd for C12H9N3O4, 259.06). Anal. Calcd for C12H9N3O4: C, 55.60; H, 3.50 and N, 16.21; found: C, 55.53; H, 3.52 and N, 16.24.

2.4. Synthesis of 2-Azido-4,4′-Dinitro-1,1′-Biphenyl (DPNN3)

Sulfuric acid (10 mL, 6 mol/L) was added into a three-necked flask, and then cooled to 0 °C to remove oxygen. After that, DPNA (0.518 g, 2 mmol) was added with stirring at 0 °C for 0.5 h. Subsequently, NaNO2 (0.207 g, 3 mmol) and NaN3 (0.195 g, 3 mmol) aqueous solutions with a concentration of 1 mol/L were slowly added with stirring at 0 °C for 3 h. Then, the intermediate (DPNN3) was extracted by 30 mL of dichloromethane, and washed by saturated salt water to remove impurities. The yield of the product was about 72% (0.411 g). IR (KBr, v, cm−1): 2127 (R-N=N=N stretching), 1509 (−NO2 stretching), 1347 (C−N stretching), 1100~700 (Ar–H stretching). 1H NMR (400 MHz, DMSO-d6, δ, ppm): 8.33 (d, J = 8.7 Hz, 2H), 8.16 (s, 1H), 8.14 – 8.06 (m, 1H), 7.83 (d, J = 8.7 Hz, 2H), 7.73 (d, J = 8.4 Hz, 1H). 13C NMR (100 MHz, DMSO-d6, δ, ppm): 148.51, 147.77, 142.87, 139.20, 136.83, 132.65, 131.29, 123.86, 120.25, 115.22. MS (EI, m/z): 285(10) ([M+], calcd for C12H7N5O4, 285.05). Anal. Calcd for C12H7N5O4: C, 50.53; H, 2.47 and N, 24.55; found: C, 50.49; H, 2.50 and N, 24.59.

2.5. Synthesis of 2,7-Dinitro-9H-Carbazole (2,7-CDN)

The mixture of DPNN3 (0.285 g, 1 mmol) and decalin (15 mL) was placed in a three-necked flask under argon atmosphere, and then heated to 155 °C with stirring for 6 h. After cooling down to room temperature, the product (2,7-CDN) was concentrated and purified by silica-gel column chromatography using ethyl acetate/n-hexane (v/v = 1/3), and then vacuum-dried at 80 °C for 24 h. The yield of the product was about 65 % (0.167 g). IR (KBr, v, cm−1): 1510 (−NO2 stretching), 1345 (C−N stretching), 3379 (N-H stretching), 1100~700 (Ar–H stretching). 1H NMR (400 MHz, DMSO-d6, δ, ppm): 12.23 (s, 1H), 8.47 (d, J = 8.8 Hz, 4H), 8.07 (dd, J = 8.7, 1.8 Hz, 2H). 13C NMR (100 MHz, DMSO-d6, δ, ppm): 147.08, 141.28, 126.47, 123.08, 115.02, 108.40. MS (EI, m/z): 257(100) ([M+], calcd for C12H7N3O4, 257.04). Anal. Calcd for C12H7N3O4: C, 56.04; H, 2.74 and N, 16.34; found: C, 56.00; H, 2.78 and N, 16.36.

2.6. Synthesis of 9H-Carbazole-2,7-Diamine (2,7-CDA)

2,7-CDN (2.572 g, 10 mmol) and absolute ethanol (450 mL) were added into a three-necked flask with stirring for 0.5 h at room temperature under argon atmosphere to remove oxygen. Then, the mixture was heated to 80 °C, and Pd/C (0.02 g) and hydrazine hydrate (5 mL) were added with stirring at 80 °C for 24 h under argon atmosphere. After cooling to room temperature, the Pd/C was removed by filtration and the grey white diamine (2,7-CDA) was obtained by lower temperature crystallization. Yield: 93 % (1.834 g). IR (KBr, v, cm−1): 3390 (N-H stretching), 3318 (-NH2 stretching), 1617 (δ N-H), 1273 (C-N stretching), 1100~700 (Ar–H stretching). 1H NMR (400 MHz, DMSO-d6, δ, ppm): 10.24 (s, 1H), 7.43 (d, J = 8.2 Hz, 2H), 6.49 (d, J = 1.7 Hz, 2H), 6.34 (dd, J = 8.2, 1.8 Hz, 2H), 4.87 (s, 4H). 13C NMR (100 MHz, DMSO-d6, δ, ppm): 146.09, 141.63, 119.10, 114.89, 107.69, 95.39. MS (EI, m/z): 197(100) ([M+], calcd for C12H11N3, 197.10). Anal. Calcd for C12H11N3: C, 73.07; H, 5.62 and N, 21.30; found: C, 73.02; H, 5.65 and N, 21.32.

2.7. Synthesis of Polyimides

In clean rooms, the polyimide (2,7-CPI) was synthesized following the procedure: 2,7-CDA (0.3379 g, 1.713 mmol) was dissolved in 10 mL of DMF in a 50-mL flask under argon atmosphere. Then, PMDA (0.3736 g, 1.713 mmol) was added, reaching a solid content of about 7 wt%. The reaction mixture was stirred at 0 °C for 8 h to produce a viscous poly(amic acid) (PAA) solution. Then, the PAA is defoamed and coated uniformly on a glass plate. Subsequently, the thermal imidization was carried out in a vacuum oven with the temperature program of 100 °C (1 h) / 100~200 °C (1 h) / 200~300 °C (1 h)/300~400 °C (1 h) to form 2,7-CPI film. The 2,7-CPI film was then released from the glass plate after the oven cooled down to room temperature. IR (KBr, v, cm−1): 1775 and 1713 (C=O stretching), 1359 (C–N stretching), 1100~700 (Ar–H stretching). Anal. Calcd for C22H9N3O4: C, 69.66; H, 2.39 and N, 11.08; found: C, 69.25; H, 2.55 and N, 11.41.

2.8. Molecular Simulation

2.8.1. Construction of Polymer Microstructures

Biovia Materials Studio software was employed to perform simulations in this study using the COMPASS (Condensed-phase Optimized Molecular Potentials for Atomistic Simulation Studies) forcefield [24,25]. Initially, polymer chains containing 25 repeat units were built and geometry optimization was performed. Then, a periodic model of the PI, comprising of five polymer chains in cubic unit cell, was constructed and the total energy of the system was minimized using smart minimizing method. Annealing was performed by the NPT (constant number of particles (N), pressure (P), and temperature (T)) dynamics procedure through heating and cooling the system at 1 atm in the temperature range of 300 to 1000 K in steps of 50 K. Interactions of non-bond, van der Waals, and electrostatic forces, were calculated using an atom-based summation method and an Ewald summation method, respectively. The annealed cell was later put through a stage-wise equilibration procedure. First the cell was heated to 1000 K, and then the temperature was decreased in several stages to 600 K in steps of 100 K. After this, the temperature was decreased in several stages to 400 K in steps of 50 K and then decreased to 300 K in steps of 25 K. Each stage consisted of two consecutive NVT (constant number of particles (N), volume (V) and temperature (T)) and NPT runs at 1 atm and a specific temperature. The aim of the procedure was to obtain a refined system that would relax at the experimental density of the amorphous polymer at 1 atm and 300 K. Finally, the cell was relaxed by consecutive NVT (at 300 K) and NPT dynamics (at 1 atm and 298 K) to ensure that a constant density has been reached. Two criteria were used to determine the equilibrium of the system: (1) The density of the system remained stable for a long time; (2) the fluctuation of energy was lower than 10 % [26]. The plots of density and energy versus simulation time in the last NPT for PI systems are shown in supporting information Figures S1 and S2. It can be seen that the PI systems have reached equilibrium states. The equilibrated densities of Kapton, 2,7-CPPI and 2,7-CPI systems are 1.41, 1.56, and 1.57 g/cm3, respectively, which are fairly close to the experimental densities (1.42, 1.57 and 1.58 g/cm3). In all runs, Nosé method was used for temperature control [27]. In NPT runs, the pressure was controlled by Berendsen’s method. During these simulations, the cutoff for the nonbonded interactions was taken as 15.5 Å.

2.8.2. Free Volume

The free volume was determined by a grid scanning method using the Connolly task. The void distribution was estimated by a method previously used for micro-crystalline materials [28,29]. Specifically, the simulation cell was divided into three-dimensional fine grids with a size of approximately 0.25 Å. The void size at a grid was determined as the diameter of the maximum cavity that encloses the grid and additionally has no overlap with any polymer atom.

2.8.3. Radial Distribution Functions

The radial distribution functions (RDF) refer to a measure of the probability that, given the presence of an atom at the origin of an arbitrary reference frame, there will be an atom with its center located in a spherical shell of infinitesimal thickness at a distance r from the reference atom. The RDF was calculated by the average of the static relationship of every given pair of particles A B using the following equation [30]:
g A B ( r ) = n A B ( r ) 4 π r 2 Δ ρ A B
where n A B ( r ) is the average number of atom pairs between r and r + Δ r , and Δ ρ A B is the density of atom pairs of type A B .

2.8.4. Diffusion Coefficients

The diffusivity of gas molecules through PI was estimated by first inserting ten molecules of each gas into the equilibrated simulation box. Then, minimization of the potential energy was performed using “smart minimizing method” run. After this, the cell was put through annealing and stage-wise equilibration procedures using the same parameters as described before. The resulting structure was then equilibrated by NVT and NPT simulations at 298 K in order to ensure that its minimized total energy remained approximately constant with respect to the simulation time. An NVE simulation of the system was performed for 10,000 ps. The diffusion coefficients can be calculated by means of the Einstein relation [31,32]:
D = 1 6 N lim t d d t i N | r i ( t ) r i ( 0 ) | 2
where N is the number of penetrants, r i ( 0 ) and r i ( t ) are the initial and final positions of the center of mass of penetrant i over the time interval t , and | r i ( t ) r i ( 0 ) | 2 is the averaged mean-square displacement (MSD) of the penetrant. The diffusion coefficient was determined from the slope of MSD versus time data. In this work, MSD of H2O and O2 were calculated from the trajectories of ten penetrant molecules in the PI microstructures.

2.8.5. Determination of Local Mobility of Polymer Chains

The backbone mean square displacements (MSD) of PIs were determined to investigate the local mobility of PI chains. An NVE simulation of the equilibrated system was performed for 5000 ps to study the backbone MSD.

2.8.6. Dihedral Angle Profiles

Dihedral angle profiles are used to examine the mobility of linking groups of polymer chains. Four consecutive atoms in the backbone are used to calculate the dihedral angle θ. The dihedral angle is defined as the angle between the two planes formed by the first three and last three atoms. For the most extended planar conformation, the dihedral angle θ is ±180°. Figure S3 illustrates the schematic diagrams of the analyzed atomic segments for the dihedral angle analysis.

2.8.7. Sorption Isotherm

The equilibrated cell was used for grand canonical Monte Carlo (GCMC) simulations employing the standard Metropolis algorithm using the “Sorption Isotherm” module [33]. Both the polymer framework and the penetrant molecules were treated as rigid bodies. The degrees of freedom of the system were accordingly specified by the center-of-mass position and orientation of the molecules. Metropolis sampling was used for inserting or deleting permeant molecules as well as accepting or rejecting their translational and rotational configurational moves. The COMPASS force field and force field assigned partial charges on atoms were used. A VTμ simulation was performed at each fixed pressure and 298 K. The pressure of the penetrant gas was varied from 10 to 3000 kPa. For each pressure value, 105 equilibration steps were first performed to ensure proper relaxation of the polymer chains in response to the insertion of the penetrant molecule, following which 106 steps of production run were carried out. The sorption isotherm can be obtained in the form of a plot of the concentration of sorbed gas, C , as a function of pressure at constant temperature. The solubility coefficient, S , is then obtained from the limiting slope of the sorption isotherm at zero pressure as [34]:
S = lim p 0 ( C p )
where C is in units of cm3(STP)/cm3(polymer) and p is pressure.

3. Results and Discussion

3.1. Synthesis and Characterization of Monomers and Polyimides

Scheme 1 presents the synthesis route of 2,7-CDA. First, DPNA was prepared from o-aminobiphenyl by nitrification reaction. Second, DPNN3 was synthesized from DPNA by azide reaction. Subsequently, the as-obtained DPNN3 was reacted with decalin to form 2,7-CDN by cyclization reaction, which was then reduced to diamine 2,7-CDA. NMR, FT-IR, mass spectra, and elemental analysis were used to characterize the structures of the diamine and intermediate. The NMR, MS spectra of the intermediates DPNA, DPNN3, and 2,7-CDN and MS spectrum of 2,7-CDA are presented in supporting information Figures S4–S7. The NMR results of 2,7-CDA are given in Figure 2. The FT-IR spectra of diamine and intermediate are shown in Figure S8. The elemental analysis results are shown in experimental section. All these results were consistent with the proposed molecular structure of diamine and intermediate, confirming the successful synthesis of diamine by the mentioned four-step reactions.
The conventional two‒step thermal imidization method was used to prepare polyimide, as shown in Scheme 2. The 2,7-CDA was reacted with PMDA in DMF to form a precursor PAA, which was then converted to 2,7-CPI by thermal cyclodehydration. GPC was employed to measure the molecular weight of PAA. The weight-average molecular weight (Mw) of the PAA was 3.67 × 105, with a polydispersity index (Mw/Mn) of 1.96. The FT-IR was used to confirm the formation of 2,7-CPI in Figure S8. Compared with 2,7-CDA, the absorption peaks at 3318–3390 cm−1 (N‒H stretching) and 1617 cm−1 (δ N‒H) disappeared in 2,7-CPI, meanwhile the characteristic peaks of the imide group at 1359 cm−1 (stretching of C‒N) and 1775, 1713 cm−1 (stretching of carbonyl) were observed, suggesting the successful reaction between 2,7-CDA and PMDA, and the complete imidization reaction of the PAA.

3.2. Thermal and Mechanical Properties

The thermal performances of 2,7-CPI were studied by DMA, TMA, and TGA, which are also compared with those of a typical polyimide (Kapton) and our previous reported structural analogue 2,7-CPPI. The results are listed in Table 1. The 2,7-CPI displayed the highest Tg of 467 °C among the three PIs (Figure 3a). Additionally, the 2,7-CPI possessed outstanding dimensional stability, with a CTE low to 3.4 ppm/K in the range of 100~300 °C (Figure 3b). Furthermore, the 2,7-CPI exhibited exceptional thermal stability, with Td5% and Td10% of 550 °C and 590 °C, respectively (Figure S9). The tensile strength and modulus of 2,7-CPI were 134 MPa and 5.7 GPa, respectively. The dimensional and thermal stability and mechanical properties of 2,7-CPI were comparable to those of 2,7-CPPI, but better than Kapton. The excellent thermal and mechanical performances of 2,7-CPI and 2,7-CPPI were mainly attributed to the aromatic rigid carbazole moieties in the PI main chains and the strong interchain cohesion. Also, the 2,7-CPI had good flexibility, which was indicated from the rolled film in Figure S10. The excellent comprehensive performances of 2,7-CPI and 2,7-CPPI make them promising materials in the micro-electronics packaging industry.

3.3. Barrier Properties

The barrier performances of 2,7-CPI film are compared with Kapton and 2,7-CPPI, as shown in Table 2. Herein the diamine ODA used for the synthesis of Kapton has an ether bond between the two benzene rings (see Figure 1). The 2,7-CPI showed better barrier performances than Kapton and 2,7-CPPI. Its WVTR and OTR were low to 0.05 g·m−2·day−1 and 0.14 cm3·m−2·day−1, respectively, which reduced by three orders of magnitude in comparison with Kapton. Compared with 2,7-CPPI, the WVTR and OTR of 2,7-CPI showed reductions of 50% and 30%, respectively. The 2,7-CPI also showed superior barrier properties in comparison with most reported intrinsic polyimides [36,37] and other barrier polymers such as PET [38], ethylene-vinyl alcohol copolymer [39], nylon-6 [40], and polyvinylidene chloride [41].

3.4. Aggregation Structures Analysis

In order to reveal the better barrier properties of 2,7-CPI in comparison with Kapton and 2,7-CPPI, WAXD was carried out to study the packing degree and order of polymer chains. The experimental WAXD scattering spectra of the three PIs are presented in Figure 4. The 2,7-CPPI and 2,7-CPI presented strong diffraction peaks, suggesting their higher crystallinity, while the broad bands of the Kapton demonstrated the amorphous nature of the polymer. The crystals in 2,7-CPPI and 2,7-CPI matrices will hinder the permeation of H2O and O2 by forcing them to follow a more tortuous diffusion pathway, and thereby improving their barrier properties [42]. The peak in WAXD patterns is ascribed to the chain-to-chain distance (d-spacing) between the polymer backbones [43]. The 2,7-CPI showed an intensive reflection peak at 22.11° corresponding to the d-spacing of 4.01 Å, which was lower than those of Kapton and 2,7-CPPI (Table 3). In addition, 2,7-CPI presented higher density than Kapton and 2,7-CPPI, as shown in Table 3. The low d-spacing and high crystallinity and density suggested the close chains packing of 2,7-CPI. The Kapton displayed lowest density and largest d-spacing among the three PIs, demonstrating its loose chains packing. This was mainly due to the nonplanar structure of ODA caused by the ether bond.

3.5. Hydrogen Bonds Analysis

Intermolecular interactions play an important role in chains packing. The radial distribution functions (RDFs) were calculated to study the intermolecular H-bonding interactions in PI matrices. In theory, the hydrogen bonds between ‒NH‒ of carbazole and O=C‒ of imide rings could be formed in the 2,7-CPPI and 2,7-CPI matrices, while no hydrogen bonds could be formed in Kapton. The radial distribution functions (gAB(r)) of 2,7-CPPI and 2,7-CPI were calculated to analyze the hydrogen atoms of ‒HN‒ in carbazole and the oxygen atoms of O=C‒ in imide rings, and the relevant calculated results are shown in Figure 5. Intermolecular interactions can be classified into H-bonding, strong and weak van der Waals’ forces, corresponding to the distances between atoms of 2.6~3.1, 3.1~5.0, and above 5.0 Å, respectively [44]. In Figure 5, there were peaks appearing in the range of 2.6~3.1 Å, suggesting the formation of hydrogen bonds in 2,7-CPPI and 2,7-CPI. Notably, 2,7-CPI exhibited higher g(r) value than 2,7-CPPI in that range, indicating the larger number of hydrogen bonds in 2,7-CPI compared to that in 2,7-CPPI. Figure 6 presents the hydrogen bonds formed in the simulation cell of 2,7-CPI, confirming the existence of the above-mentioned hydrogen bonds. Statistical results showed that 49 and 60 hydrogen bonds were formed in the 2,7-CPPI and 2,7-CPI systems, respectively (Table 3), while no hydrogen bonds were found in the Kapton system. The cohesive energy density (CED) of the three PIs systems were calculated by molecular simulations and shown in Table 3. The 2,7-CPI showed the highest CED among the three PIs, probably owing to the more number of hydrogen bonds formed in 2,7-CPI matrix. The strong intermolecular hydrogen bonding of 2,7-CPI contributed to its close chain packing. The low CED value of Kapton was also a reason for its loose chains packing.

3.6. Free Volume and Cavity Size Distribution Analysis

The above-discussed chain packing and interaction will affect the free volumes formed in PI matrix, which play a vital role in the gas transport performances. PALS was utilized to study the free volume of the PI films. The positron lifetime curves of Kapton, 2,7-CPPI and 2,7-CPI film are shown in Figure 7. The slope of the decay curve for Kapton was lower than those of 2,7-CPPI and 2,7-CPI, indicating the larger size of the free volume in Kapton [45]. The raw data obtained were resolved into two finite lifetime components using the PATFIT program which assumes a Gaussian distribution for the logarithm of the lifetime for each component. The obtained data of positron lifetime are listed in Table S1. Jean et al. have reported a relation between the second lifetime component τ2 and the mean free-volume radius (R) of polymeric materials [46]. According to such correlation, the R of PI are calculated and shown in Table 3. In Table 3, the relative fractional free volume (FFV) is also calculated according to our previous work [21]. The FFV of the three PIs decreased in the order: Kapton ˃ 2,7-CPPI ˃ 2,7-CPI.
Molecular simulations were used to study the characteristics of free volumes including size distribution, three-dimensional arrangement, and connectivity. The distributions of void size in Kapton, 2,7-CPPI and 2,7-CPI are shown in Figure 8a. Compared to Kapton, the 2,7-CPPI and 2,7-CPI had fewer voids with radius bigger than 1.2 Å. The kinetic radii of O2 and H2O are 1.73 Å and 1.325 Å, respectively, which are larger than 1.2 Å. This suggested that 2,7-CPPI and 2,7-CPI had fewer voids available for the diffusion of H2O and O2. In addition, 2,7-CPI possessed fewer voids with radius bigger than 0.8 Å in comparison with 2,7-CPPI, indicating the fewer number of large-size free volumes in 2,7-CPI.
In order to understand the available free volume for various penetrants, the FFV was analyzed using spherical probes with different radii. Figure 8b illustrates the FFV as a function of probe radius for Kapton, 2,7-CPPI and 2,7-CPI. The total fractional free volume (FFV0) for the three PIs was obtained when the probe radius was set as 0 Å. As expected, the FFV decreased with an increase in the probe radius. It is observed that 2,7-CPI showed the lowest FFV among the three PIs in all the related probe radius range. For reference, the kinetic radii of O2 and H2O considered in this study are indicated on the figure. The 2,7-CPI presented the lowest FFV0, H2O accessible FFV and O2 accessible FFV, as shown in Table 3. This was consistent with the FFV results of PALS. Furthermore, the morphological representation of the O2 and H2O accessible volumes in Kapton, 2,7-CPPI and 2,7-CPI are shown in Figure 9. The voids in 2,7-CPI were smaller and discontinuous compared with those of Kapton and 2,7-CPPI. This was due to the close chain packing caused by the better coplanar structures and strong intermolecular hydrogen bonding of 2,7-CPI. The less connected and low accessible free volume of 2,7-CPI resulted in low gas permeability and thereby high gas barrier. The Kapton showed larger and well-connected voids resulting from the nonplanar structure and low interchain cohesion, which led to high gas permeability.

3.7. Local Mobility of Polymer Chains

In addition to free volume, the local mobility of polymer chains also has a vital impact on the gas transport performances of polymers. The backbone mean square displacements (MSD) of Kapton, 2,7-CPPI and 2,7-CPI systems were calculated and shown in Figure 10. It is evident that the MSD of the three PIs decreased in the order: Kapton ˃ 2,7-CPPI ˃ 2,7-CPI. The introduction of ether bonds into the main chains of Kapton facilitated its chains mobility, while the rigid planar carbazoles in 2,7-CPPI and 2,7-CPI hindered the chains mobility. The 2,7-CPI showed lower chains mobility than 2,7-CPPI. The dihedral angle profiles of the polymer chain segments were employed to analyze the mobility of linking groups in 2,7-CPPI and 2,7-CPI, as shown in Figure 11. Compared to 2,7-CPPI, the dihedral angle peaks of the 2,7-CPI were relatively higher and narrower, demonstrating that the rotational freedoms of imide rings and carbazoles along the axis of C‒N bond were restricted [47]. This was a reason for the lower chains mobility of 2,7-CPI. In addition, the more number of hydrogen bonds formed in 2,7-CPI also contributed to its lower chains mobility. The lower chains mobility was another reason for the high barrier properties observed in 2,7-CPI. The mobility of polymer chains has a crucial effect on the channel formation for the gas molecules to permeate [48]. The poor chain mobility of 2,7-CPI reduced the degree of channel formation for gas diffusion, thus led to the high gas barrier. The high chain mobility of Kapton contributed to its high gas permeability.

3.8. Gas Transport Behaviour

In order to better understand the barrier performances of the three PIs, gas transport behavior including gas diffusion and solubility are investigated.

3.8.1. Gas Diffusion

The diffusion behaviors of O2 and H2O through the PI matrices were analyzed. The representative trajectories of O2 and H2O molecules in the three PI systems during 10 ns simulation run are shown in Figure 12. Additionally, the displacements |ri (t) − ri (0)| of O2 and H2O are presented in Figure S11. The curves for 2,7-CPPI and 2,7-CPI have been shifted vertically for better visualization in Figure S11. As shown in Figure 12 and Figure S11, the diffusions of O2 and H2O molecules in the polymer matrix involved random oscillations within the available free volume and occasional jumps from one voids to another [49]. The length of the trajectories for O2 and H2O molecules in the three PI matrices decreased in the order: Kapton ˃ 2,7-CPPI ˃ 2,7-CPI. The jumps frequency and lengths of O2 and H2O in the three PI matrices followed the same trend (Figure S11). It also can be seen that H2O showed lower mobility than O2 in a particular PI matrix.
The log(MSD) of O2 and H2O in the three PI systems was plotted as a function of time in Figure 13. The diffusion coefficient was then calculated by the Einstein’s equation in the normal diffusion regime, where the slope of the log(MSD) vs. log(t) curve was unity [50]. The obtained diffusion coefficients are listed in Table 4. The diffusion coefficients of O2 and H2O in the three PI matrices decreased in the order: Kapton ˃ 2,7-CPPI ˃ 2,7-CPI, which was consistent with the free volumes and chain mobility results. The 2,7-CPI had small FFV, poor cavities connectivity, and low chains mobility, which restricted the diffusion of H2O and O2 and resulted in the lowest diffusion coefficients. The high FFV, good cavities connectivity, and high chains mobility of Kapton facilitated the gas diffusion. As for a particular PI matrix, H2O presented lower diffusion coefficient compared with O2. Although the kinetic diameter of H2O was smaller, its diffusivity was lower than that of O2. This could be due to the strong interaction between polar H2O and PI matrix [51], which hindered the diffusion of H2O.

3.8.2. Gas Solubility

To understand the gas sorption behavior, the sorption isotherms of O2 and H2O in the three PI matrices were investigated by GCMC simulations. The results are illustrated in Figure 14. The initially rapid uptake at low pressures followed by a slower gain in loading at high pressures was observed. At low pressures, the Langmuir-mode adsorption at the micro-voids was dominant, while the Henry-mode adsorption at the inter-polymer chain free volume was primary at high pressures [52]. The sorption data of O2 and H2O fitted well to dual mode sorption model and the obtained solubility coefficients are listed in Table 4. The solubility coefficients of O2 and H2O in the three PI matrices showed the trend of Kapton ˃ 2,7-CPPI ˃ 2,7-CPI. The smaller size and fewer number of cavities in 2,7-CPI resulted in fewer gas sorption sites and thereby the lowest solubility coefficients. The large size and more number of cavities in Kapton led to its highest solubility coefficients. Additionally, H2O presented larger solubility coefficient than O2 in a particular PI matrix, which was contrary to the trend of their diffusion coefficients. This was mainly due to the strong affinities between H2O and PI matrix [38], high critical temperature (H2O: 647 K, O2: 155 K) and small size of H2O, which favored the adsorption processes of H2O in the polymer matrix.

3.8.3. Gas Permeability

According to solution-diffusion mechanism, permeability coefficient (P) can be obtained by the product of the diffusion coefficient and solubility coefficient. Table 4 shows the permeability coefficients of O2 and H2O in the three PIs. The simulated permeabilities of O2 and H2O in the three PIs showed similar trends as the experimental results (Table 2), validating the feasibility of simulations for gas transport analysis. The permeability coefficients of O2 and H2O in the three PIs decreased in the following the order: Kapton ˃ 2,7-CPPI ˃ 2,7-CPI, which was in agreement with the trends of diffusion and solubility coefficients. It can be seen that the reductions of diffusion and solubility coefficients of O2 and H2O in 2,7-CPI gave rise to the decreases of permeability coefficients. Consequently, the barrier performances of 2,7-CPI were enhanced. The high gas barrier of 2,7-CPI was mainly ascribed to its high crystallinity, low free volume, and poor chains mobility.

4. Conclusions

A intrinsic high-barrier polyimide (2,7-CPI) based on PMDA and a novel diamine (2,7-CDA) containing carbazole moiety was prepared. The 2,7-CPI exhibited high gas barrier, exceptional thermal stability, and good mechanical properties. Using molecular simulations, WAXD and PALS, the barrier performances of 2,7-CPI are compared with those of a structural analogue (2,7-CPPI) and a typical polyimide (Kapton) to reveal the barrier improvement of 2,7-CPI in relation to microstructure. The results showed that 2,7-CPI exhibited compact chain packing caused by the better coplanar structure and more number of intermolecular hydrogen bonds, which resulted in high crystallinity and low FFV. The Kapton showed loose chain packing resulted from the nonplanar structure and low interchain cohesion, which led to high FFV. Molecular simulations revealed that the chain mobility of the three PIs decreased in the order: Kapton ˃ 2,7-CPPI ˃ 2,7-CPI. The high crystallinity, low free volume, and decreased chain mobility reduced the diffusion and solubility of the gases in the polymer matrix, thus leading to the improvement of barrier performances for 2,7-CPI. The high FFV and chain mobility of Kapton increased the gases diffusion and solubility, thus leading to high gas permeability. The excellent comprehensive performances of 2,7-CPI suggest its promising application in micro-electronics encapsulation and high grade packaging industry.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4360/12/9/2048/s1. Figures S1–S2: Plots of density and energy versus simulation time in the NPT simulation for PIs. Figure S3: Schematic diagrams of the analyzed atomic segments for the dihedral angle analysis. Figures S4–S6: 1H NMR, 13C NMR spectra and mass spectrum of DPNA, DPNN3, 2,7-CDN. Figure S7: Mass spectrum of 2,7-CDA. Figure S8: FT-IR spectra of DPNA, DPNN3, 2,7-CDN, 2,7-CDA, and 2,7-CPI. Figures S9–S10: TGA curve and photo image of 2,7-CPI film. Figure S11 and Table S1: The positron lifetime data and displacement of O2 and H2O for Kapton, 2,7-CPPI, and 2,7-CPI.

Author Contributions

Y.L. (Yiwu Liu) and J.T.: writing—original draft, writing—review and editing, conceptualization, project administration, funding acquisition; A.T., X.Z., Y.L. (Yuhui Li), and P.H.: project administration, resources. Investigation, methodology; C.C. and D.W.: formal analysis, investigation, methodology; H.Z.: writing—review and editing, supervision, project administration. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China (No. 51603066), NSAF (U1730142, U1930203), the Natural Science Foundation of Hunan Province (No. 2019JJ40071), the Science Research Project of Hunan Provincial Department of Education (No. 18A257) and the Student Innovation and Entrepreneurship Training Program of China (No. 201911535005).

Acknowledgments

We would like to thank Zhiquan Chen from Hubei Nuclear Solid Physics Key Laboratory, Department of Physics, Wuhan University for their help in the measurement of the polymer free volume.

Conflicts of Interest

The authors declare no competing financial interest.

References

  1. Huang, Y.; Hsiang, E.-L.; Deng, M.-Y.; Wu, S.-T. Mini-LED, Micro-LED and OLED displays: Present status and future perspectives. Light. Sci. Appl. 2020, 9, 1–16. [Google Scholar] [CrossRef] [PubMed]
  2. Zou, S.-J.; Shen, Y.; Xie, F.-M.; Chen, J.-D.; Li, Y.-Q.; Tang, J.-X. Recent advances in organic light-emitting diodes: Toward smart lighting and displays. Mater. Chem. Front. 2020, 4, 788–820. [Google Scholar] [CrossRef]
  3. Choi, M.; Park, Y.J.; Sharma, B.K.; Bae, S.-R.; Kim, S.Y.; Ahn, J.-H. Flexible active-matrix organic light-emitting diode display enabled by MoS2thin-film transistor. Sci. Adv. 2018, 4, eaas8721. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Lee, H.-Y.; Su, P.-S.; Lee, C.-T. Micromesh-structured flexible polymer white organic light-emitting diodes using single emissive layer of blended polymer and quantum dots. Org. Electron. 2020, 82, 105722. [Google Scholar] [CrossRef]
  5. Mamo, M.D.; Noh, Y.-Y.; Shin, E.S. Self-aligned patterning of conductive films on plastic substrates for electrodes of flexible electronics. J. Mater. Chem. C 2017, 5, 10900–10906. [Google Scholar] [CrossRef]
  6. Macdonald, W.A.; Looney, M.K.; MacKerron, D.; Eveson, R.; Adam, R.; Hashimoto, K.; Rakos, K. Latest advances in substrates for flexible electronics. J. Soc. Inf. Disp. 2007, 15, 1075–1083. [Google Scholar] [CrossRef]
  7. Thomschke, M.; Reineke, S.; Lüssem, B.; Leo, K. Highly Efficient White Top-Emitting Organic Light-Emitting Diodes Comprising Laminated Microlens Films. Nano Lett. 2011, 12, 424–428. [Google Scholar] [CrossRef]
  8. Gardonio, S.; Gregoratti, L.; Melpignano, P.; Aballe, L.; Biondo, V.; Zamboni, R.; Murgia, M.; Caria, S.; Kiskinova, M. Degradation of organic light-emitting diodes under different environment at high drive conditions. Org. Electron. 2007, 8, 37–43. [Google Scholar] [CrossRef]
  9. Choi, M.-C.; Kim, Y.; Ha, C.-S. Polymers for flexible displays: From material selection to device applications. Prog. Polym. Sci. 2008, 33, 581–630. [Google Scholar] [CrossRef]
  10. Gao, X.; Lin, L.; Liu, Y.; Huang, X. LTPS TFT Process on Polyimide Substrate for Flexible AMOLED. J. Disp. Technol. 2015, 11, 666–669. [Google Scholar] [CrossRef]
  11. Salleo, A.; Street, R.A. Publisher’s Note: Light-induced bias stress reversal in polyfluorene thin-film transistors. J. Appl. Phys. 2003, 94, 4231. [Google Scholar] [CrossRef] [Green Version]
  12. Kao, S.-C.; Li, L.-J.; Hsieh, M.-C.; Zhang, S.; Tsai, P.-M.; Sun, Z.-Y.; Wang, D. 71-1: Invited Paper: The Challenges of Flexible OLED Display Development. SID Symp. Dig. Tech. Pap. 2017, 48, 1034–1037. [Google Scholar] [CrossRef]
  13. Tetsuka, H.; Ebina, T.; Tsunoda, T.; Nanjo, H.; Mizukami, F. Flexible organic electroluminescent devices based on transparent clay films. Nanotechnology 2007, 18, 355701. [Google Scholar] [CrossRef]
  14. Lee, J.; Lee, Y.; Kang, T.; Chu, H.; Kwag, J. Alleviation of abnormal NBTI phenomenon in LTPS TFTs on polyimide substrate for flexible AMOLED. J. Soc. Inf. Disp. 2020, 28, 333–341. [Google Scholar] [CrossRef]
  15. Park, J.-S.; Kim, T.; Stryakhilev, D.; Lee, J.-S.; An, S.-G.; Pyo, Y.-S.; Lee, D.-B.; Mo, Y.G.; Jin, D.U.; Chung, H.K. Flexible full color organic light-emitting diode display on polyimide plastic substrate driven by amorphous indium gallium zinc oxide thin-film transistors. Appl. Phys. Lett. 2009, 95, 13503. [Google Scholar] [CrossRef]
  16. Seo, S.-W.; Chung, H.K.; Chae, H.; Seo, S.J.; Cho, S.M. Flexible organic/inorganic moisture barrier using plasma-polymerized layer. Nano 2013, 8, 1350041. [Google Scholar] [CrossRef]
  17. Kwon, J.H.; Jeong, E.G.; Jeon, Y.; Kim, D.-G.; Lee, S.; Choi, K.C. Design of Highly Water Resistant, Impermeable, and Flexible Thin-Film Encapsulation Based on Inorganic/Organic Hybrid Layers. ACS Appl. Mater. Interfaces 2018, 11, 3251–3261. [Google Scholar] [CrossRef]
  18. Tsai, M.-H.; Chang, C.-J.; Lu, H.-H.; Liao, Y.-F.; Tseng, I.-H. Properties of magnetron-sputtered moisture barrier layer on transparent polyimide/graphene nanocomposite film. Thin Solid Films 2013, 544, 324–330. [Google Scholar] [CrossRef]
  19. Sanaeepur, H.; Amooghin, A.E.; Bandehali, S.; Moghadassi, A.; Matsuura, T.; Van Der Bruggen, B. Polyimides in membrane gas separation: Monomer’s molecular design and structural engineering. Prog. Polym. Sci. 2019, 91, 80–125. [Google Scholar] [CrossRef]
  20. Madkour, T.; Madkour, T.M. Development of the molecular design rules of ultra-permeable poly[1-(trimethylsilyl)-1-propyne] membranes. Polymer 2000, 41, 7489–7497. [Google Scholar] [CrossRef]
  21. Liu, Y.; Huang, J.; Tan, J.; Zeng, Y.; Liu, J.; Zhang, H.; Pei, Y.; Xiang, X.; Liu, Y. Intrinsic high-barrier polyimide with low free volume derived from a novel diamine monomer containing rigid planar moiety. Polymer 2017, 114, 289–297. [Google Scholar] [CrossRef]
  22. Yang, Y.; Nair, A.K.N.; Sun, S. Sorption and Diffusion of Methane and Carbon Dioxide in Amorphous Poly(alkyl acrylates): A Molecular Simulation Study. J. Phys. Chem. B 2020, 124, 1301–1310. [Google Scholar] [CrossRef] [PubMed]
  23. Velioğlu, S.; Ahunbay, M.G.; Tantekin-Ersolmaz, S.B. Propylene/propane plasticization in polyimide membranes. J. Membr. Sci. 2016, 501, 179–190. [Google Scholar] [CrossRef]
  24. Sun, H. COMPASS: An ab Initio Force-Field Optimized for Condensed-Phase ApplicationsOverview with Details on Alkane and Benzene Compounds. J. Phys. Chem. B 1998, 102, 7338–7364. [Google Scholar] [CrossRef]
  25. McQuaid, M.J.; Sun, H.; Rigby, D. Development and validation of COMPASS force field parameters for molecules with aliphatic azide chains. J. Comput. Chem. 2003, 25, 61–71. [Google Scholar] [CrossRef]
  26. Zhu, J.; Zhao, X.; Liu, L.; Song, M.; Wu, S. Quantitative relationships between intermolecular interaction and damping parameters of irganox-1035/NBR hybrids: A combination of experiments, molecular dynamics simulations, and linear regression analyses. J. Appl. Polym. Sci. 2018, 135, 46202. [Google Scholar] [CrossRef]
  27. Nosé, S. A unified formulation of the constant temperature molecular dynamics methods. J. Chem. Phys. 1984, 81, 511. [Google Scholar] [CrossRef] [Green Version]
  28. Ban, S.; Vlugt, T.J. Zeolite microporosity studied by molecular simulation. Mol. Simul. 2009, 35, 1105–1115. [Google Scholar] [CrossRef]
  29. Bhattacharya, S.; Gubbins, K.E. Fast Method for Computing Pore Size Distributions of Model Materials. Langmuir 2006, 22, 7726–7731. [Google Scholar] [CrossRef]
  30. Zeng, J.; Zhang, Y.; Dai, Y.; Chen, S. Molecular dynamics simulation of nitrobenzene in heterocyclic ionic liquids. J. Mol. Liq. 2014, 198, 274–279. [Google Scholar] [CrossRef]
  31. Pant, P.V.K.; Boyd, R.H. Molecular-dynamics simulation of diffusion of small penetrants in polymers. Macromolecules 1993, 26, 679–686. [Google Scholar] [CrossRef]
  32. Müller-Plathe, F. Molecular dynamics simulation of gas transport in amorphous polypropylene. J. Chem. Phys. 1992, 96, 3200–3205. [Google Scholar] [CrossRef]
  33. Akkermans, R.L.; Spenley, N.A.; Robertson, S.H. Monte Carlo methods in Materials Studio. Mol. Simul. 2013, 39, 1153–1164. [Google Scholar] [CrossRef]
  34. Fried, J.; Sadat-Akhavi, M.; Mark, J. Molecular simulation of gas permeability: Poly(2,6-dimethyl-1,4-phenylene oxide). J. Membr. Sci. 1998, 149, 115–126. [Google Scholar] [CrossRef]
  35. Liu, J.; Tan, J.-H.; Zeng, Y.; Liu, Y.; Zeng, K.-J.; Liu, Y.-J.; Wu, R.-M.; Chen, H. Synthesis and characterization of high-barrier polyimide containing rigid planar moieties and amide groups. Polym. Test. 2017, 61, 83–92. [Google Scholar] [CrossRef]
  36. Tseng, I.-H.; Liao, Y.-F.; Chiang, J.-C.; Tsai, M.-H. Transparent polyimide/graphene oxide nanocomposite with improved moisture barrier property. Mater. Chem. Phys. 2012, 136, 247–253. [Google Scholar] [CrossRef]
  37. Sykes, G.F.; Clair, A.K.S. The effect of molecular structure on the gas transmission rates of aromatic polyimides. J. Appl. Polym. Sci. 1986, 32, 3725–3735. [Google Scholar] [CrossRef]
  38. Inagaki, N.; Cech, V.; Narushima, K.; Takechi, Y. Oxygen and water vapor gas barrier poly(ethylene naphthalate) films by deposition of SiOx plasma polymers from mixture of tetramethoxysilane and oxygen. J. Appl. Polym. Sci. 2007, 104, 915–925. [Google Scholar] [CrossRef]
  39. Kim, S.-W.; Cha, S.-H. Thermal, mechanical, and gas barrier properties of ethylene-vinyl alcohol copolymer-based nanocomposites for food packaging films: Effects of nanoclay loading. J. Appl. Polym. Sci. 2013, 131, 40289. [Google Scholar] [CrossRef]
  40. Gargalaka, J.G., Jr.; Couto, R.A.A.; Constantino, V.R.L.; Toma, H.E.; Araki, K. Influence of the relative amounts of crystalline and amorphous phases on the mechanical properties of polyamide-6 nanocomposites. J. Appl. Polym. Sci. 2012, 125, 3239–3249. [Google Scholar] [CrossRef]
  41. Ebina, T.; Mizukami, S. Flexible Transparent Clay Films with Heat-Resistant and High Gas-Barrier Properties. Adv. Mater. 2007, 19, 2450–2453. [Google Scholar] [CrossRef]
  42. Bai, H.; Huang, C.; Xiu, H.; Zhang, Q.; Deng, H.; Wang, K.; Chen, F.; Fu, Q. Significantly Improving Oxygen Barrier Properties of Polylactide via Constructing Parallel-Aligned Shish-Kebab-Like Crystals with Well-Interlocked Boundaries. Biomacromolecules 2014, 15, 1507–1514. [Google Scholar] [CrossRef] [PubMed]
  43. La Cruz, D.S.-D.; Green, M.D.; Ye, Y.; Elabd, Y.A.; Long, T.E.; Winey, K.I. Correlating backbone-to-backbone distance to ionic conductivity in amorphous polymerized ionic liquids. J. Polym. Sci. Part. B Polym. Phys. 2011, 50, 338–346. [Google Scholar] [CrossRef]
  44. Mattozzi, A.; Hedenqvist, M.; Gedde, U. Diffusivity of n-hexane in poly(ethylene-stat-octene)s assessed by molecular dynamics simulation. Polymer 2007, 48, 5174–5180. [Google Scholar] [CrossRef]
  45. Fong, C.; Dong, A.W.; Hill, A.J.; Boyd, B.; Drummond, C.J. Positron annihilation lifetime spectroscopy (PALS): A probe for molecular organisation in self-assembled biomimetic systems. Phys. Chem. Chem. Phys. 2015, 17, 17527–17540. [Google Scholar] [CrossRef] [PubMed]
  46. Liao, K.-S.; Chen, H.; Awad, S.; Yuan, J.-P.; Hung, W.-S.; Lee, K.-R.; Lai, J.-Y.; Hu, C.-C.; Jean, Y.C. Determination of Free-Volume Properties in Polymers Without Orthopositronium Components in Positron Annihilation Lifetime Spectroscopy. Macromolecules 2011, 44, 6818–6826. [Google Scholar] [CrossRef]
  47. Chen, Y.-R.; Chen, L.-H.; Chang, K.-S.; Chen, T.-H.; Lin, Y.-F.; Tung, K.-L. Structural characteristics and transport behavior of triptycene-based PIMs membranes: A combination study using ab initio calculation and molecular simulations. J. Membr. Sci. 2016, 514, 114–124. [Google Scholar] [CrossRef]
  48. Park, S.; Lee, A.S.; Do, Y.S.; Kim, J.F.; Hwang, S.S.; Lee, Y.M.; Lee, J.-H.; Lee, J.S. Side-chain engineering of ladder-structured polysilsesquioxane membranes for gas separations. J. Membr. Sci. 2016, 516, 202–214. [Google Scholar] [CrossRef]
  49. Neyertz, S.; Brown, D.; Pandiyan, S.; Van Der Vegt, N.F.A. Carbon Dioxide Diffusion and Plasticization in Fluorinated Polyimides. Macromol. 2010, 43, 7813–7827. [Google Scholar] [CrossRef]
  50. Mozaffari, F.; Eslami, H.; Moghadasi, J. Molecular dynamics simulation of diffusion and permeation of gases in polystyrene. Polymer 2010, 51, 300–307. [Google Scholar] [CrossRef]
  51. Zhang, K.; Yu, Q.; Zhu, L.; Liu, S.; Chi, Z.; Chen, X.; Zhang, Y.; Xu, J. The Preparations and Water Vapor Barrier Properties of Polyimide Films Containing Amide Moieties. Polymer 2017, 9, 677. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Petropoulos, J.H. A generalized, topologically consistent, dual-mode transport model for glassy polymer-gas systems. J. Polym. Sci. Part B 1989, 27, 603–620. [Google Scholar] [CrossRef]
Figure 1. The chemical structure of 2,7-CPDA, 2,7-CDA, ODA, and PMDA and the stereo view of the lowest energy conformation of repeat units for 2,7-CPPI and 2,7-CPI.
Figure 1. The chemical structure of 2,7-CPDA, 2,7-CDA, ODA, and PMDA and the stereo view of the lowest energy conformation of repeat units for 2,7-CPPI and 2,7-CPI.
Polymers 12 02048 g001
Figure 2. (a) 1H NMR (above: the spectra in the range of 0~11; below: the partial enlargement in the range of 4.5~11), (b) 13C NMR, (c) H-H COSY, (d) C-H QC and (e) C-H BC spectra of 2,7-CPA in DMSO-d6.
Figure 2. (a) 1H NMR (above: the spectra in the range of 0~11; below: the partial enlargement in the range of 4.5~11), (b) 13C NMR, (c) H-H COSY, (d) C-H QC and (e) C-H BC spectra of 2,7-CPA in DMSO-d6.
Polymers 12 02048 g002
Figure 3. (a) DMA and (b) TMA curves of the 2,7-CPI film.
Figure 3. (a) DMA and (b) TMA curves of the 2,7-CPI film.
Polymers 12 02048 g003
Figure 4. WAXD curve of Kapton, 2,7-CPPI and 2,7-CPI films.
Figure 4. WAXD curve of Kapton, 2,7-CPPI and 2,7-CPI films.
Polymers 12 02048 g004
Figure 5. The RDFs of the hydrogen atoms of ‒HN‒ in carbazole and oxygen atoms of O=C‒ in imide rings for 2,7-CPPI and 2,7-CPI.
Figure 5. The RDFs of the hydrogen atoms of ‒HN‒ in carbazole and oxygen atoms of O=C‒ in imide rings for 2,7-CPPI and 2,7-CPI.
Polymers 12 02048 g005
Figure 6. H-bonds between ‒HN‒ and O=C‒ in imide rings in simulation cell of 2,7-CPI. (red: oxygen atoms; blue: nitrogen atoms; white: hydrogen atoms; grey: carbon atoms).
Figure 6. H-bonds between ‒HN‒ and O=C‒ in imide rings in simulation cell of 2,7-CPI. (red: oxygen atoms; blue: nitrogen atoms; white: hydrogen atoms; grey: carbon atoms).
Polymers 12 02048 g006
Figure 7. Positron lifetime curve measured for Kapton, 2,7-CPPI and 2,7-CPI films.
Figure 7. Positron lifetime curve measured for Kapton, 2,7-CPPI and 2,7-CPI films.
Polymers 12 02048 g007
Figure 8. Distributions of void radius (a) and fractional free volume (b) as a function of probe radius in Kapton, 2,7-CPPI and 2,7-CPI. The kinetic radii of O2 and H2O are indicated by vertical lines.
Figure 8. Distributions of void radius (a) and fractional free volume (b) as a function of probe radius in Kapton, 2,7-CPPI and 2,7-CPI. The kinetic radii of O2 and H2O are indicated by vertical lines.
Polymers 12 02048 g008
Figure 9. Three-dimensional representation of H2O and O2 accessible volume for (a1) (a2) 2,7-CPI, (b1) (b2) 2,7-CPPI and (c1) (c2) Kapton. (grey: Van der Waals surface; blue: Connolly surface; a1, b1 and c1: probe radius of 1.325 Å; a2, b2 and c2: probe radius of 1.73 Å).
Figure 9. Three-dimensional representation of H2O and O2 accessible volume for (a1) (a2) 2,7-CPI, (b1) (b2) 2,7-CPPI and (c1) (c2) Kapton. (grey: Van der Waals surface; blue: Connolly surface; a1, b1 and c1: probe radius of 1.325 Å; a2, b2 and c2: probe radius of 1.73 Å).
Polymers 12 02048 g009
Figure 10. MSD of polymer chains as a function of time in Kapton, 2,7-CPPI and 2,7-CPI.
Figure 10. MSD of polymer chains as a function of time in Kapton, 2,7-CPPI and 2,7-CPI.
Polymers 12 02048 g010
Figure 11. Dihedral angle profiles of the 2,7-CPPI and 2,7-CPI. (Schematic diagrams of the analyzed atomic segments are shown in Figure S3).
Figure 11. Dihedral angle profiles of the 2,7-CPPI and 2,7-CPI. (Schematic diagrams of the analyzed atomic segments are shown in Figure S3).
Polymers 12 02048 g011
Figure 12. The representative trajectories of O2 and H2O in (a) Kapton, (b) 2,7-CPPI and (c) 2,7-CPI.
Figure 12. The representative trajectories of O2 and H2O in (a) Kapton, (b) 2,7-CPPI and (c) 2,7-CPI.
Polymers 12 02048 g012
Figure 13. log(MSD) of O2 and H2O in Kapton, 2,7-CPPI and 2,7-CPI as a function of time.
Figure 13. log(MSD) of O2 and H2O in Kapton, 2,7-CPPI and 2,7-CPI as a function of time.
Polymers 12 02048 g013
Figure 14. The adsorption isotherms of O2 and H2O in Kapton, 2,7-CPPI and 2,7-CPI.
Figure 14. The adsorption isotherms of O2 and H2O in Kapton, 2,7-CPPI and 2,7-CPI.
Polymers 12 02048 g014
Scheme 1. Synthesis route of 2,7-CDA.
Scheme 1. Synthesis route of 2,7-CDA.
Polymers 12 02048 sch001
Scheme 2. Synthesis route of polyimide.
Scheme 2. Synthesis route of polyimide.
Polymers 12 02048 sch002
Table 1. Thermal and mechanical properties of the 2,7-CPI film.
Table 1. Thermal and mechanical properties of the 2,7-CPI film.
PITga
(°C)
Td5%
(°C)
Td10%
(°C)
CTE b
(ppm∙K−1)
Tensile Strength
(MPa)
Tensile Modulus
(GPa)
Kapton c37455057138.0121 ± 3.02.3 ± 0.2
2,7-CPPI d4375565802.89143.8 ± 3.54.5 ± 0.2
2,7-CPI4675505903.40134 ± 3.05.7 ± 0.2
a Measured by DMA; b CTE within the range of 100~300 °C; c The values are obtained from Ref. [35]; d The values are obtained from Ref. [21].
Table 2. Barrier properties of Kapton, 2,7-CPPI and 2,7-CPI films.
Table 2. Barrier properties of Kapton, 2,7-CPPI and 2,7-CPI films.
PIsWVP
(g·mil·m−2·day−1)
OP
(cm3·mil·m−2·day−1)
WVTR
(g·m−2·day−1)
OTR
(cm3·m−2·day−1)
Kapton a114.5 ± 2.12330.7 ± 2.5038.8 ± 0.72112.0 ± 0.80
2,7-CPPI a0.3 ± 0.020.5 ± 0.030.1 ± 0.010.2 ± 0.01
2,7-CPI0.15 ± 0.010.41 ± 0.020.05 ± 0.0030.14 ± 0.007
a The barrier properties values are obtained from Ref. [21].
Table 3. Physical properties and FFV of Kapton, 2,7-CPPI and 2,7-CPI.
Table 3. Physical properties and FFV of Kapton, 2,7-CPPI and 2,7-CPI.
PIsDensity (g.cm−3)
(°)
d-Spacing (Å)CED (J/cm3)N a
(H-bonds)
Rb
(Å)
FFVb
(%)
FFVc
(O2,%)
FFVc
(H2O,%)
FFV0c
(%)
Kapton1.4218.644.7645402.6011.5010.4815.8038.51
2,7-CPPI1.5721.914.05570492.146.826.2911.5135.33
2,7-CPI1.5822.114.01613602.146.095.4810.3334.16
a The number of hydrogen bonds in the simulation cells. b The free-volume parameters tested by PALS. c Free volume parameters calculated by MD simulations; FFV(O2), FFV(H2O) and FFV0 based on probe radii of 1.73 Å, 1.325 Å and 0 Å, respectively.
Table 4. Simulated diffusion coefficients, solubility coefficients, and permeability for O2 and H2O in Kapton, 2,7-CPPI and 2,7-CPI.
Table 4. Simulated diffusion coefficients, solubility coefficients, and permeability for O2 and H2O in Kapton, 2,7-CPPI and 2,7-CPI.
PIsD aS bP c
H2OO2H2OO2H2OO2
Kapton12.889.43.920.05650.185.00
2,7-CPPI3.06.130.640.0131.920.08
2,7-CPI2.15.320.270.0080.570.04
a Units of (10−8 cm2/s). b Units of (cm3(STP)·cm−3·cmHg−1). c Units of (108 cm2·cm3(STP) s−1·cm−3·cmHg−1).

Share and Cite

MDPI and ACS Style

Liu, Y.; Tang, A.; Tan, J.; Zhao, X.; Chen, C.; Wu, D.; Li, Y.; He, P.; Zhang, H. High-Barrier Polyimide Containing Carbazole Moiety: Synthesis, Gas Barrier Properties, and Molecular Simulations. Polymers 2020, 12, 2048. https://doi.org/10.3390/polym12092048

AMA Style

Liu Y, Tang A, Tan J, Zhao X, Chen C, Wu D, Li Y, He P, Zhang H. High-Barrier Polyimide Containing Carbazole Moiety: Synthesis, Gas Barrier Properties, and Molecular Simulations. Polymers. 2020; 12(9):2048. https://doi.org/10.3390/polym12092048

Chicago/Turabian Style

Liu, Yiwu, Ao Tang, Jinghua Tan, Xianqing Zhao, Chengliang Chen, Ding Wu, Yuhui Li, Pan He, and Hailiang Zhang. 2020. "High-Barrier Polyimide Containing Carbazole Moiety: Synthesis, Gas Barrier Properties, and Molecular Simulations" Polymers 12, no. 9: 2048. https://doi.org/10.3390/polym12092048

APA Style

Liu, Y., Tang, A., Tan, J., Zhao, X., Chen, C., Wu, D., Li, Y., He, P., & Zhang, H. (2020). High-Barrier Polyimide Containing Carbazole Moiety: Synthesis, Gas Barrier Properties, and Molecular Simulations. Polymers, 12(9), 2048. https://doi.org/10.3390/polym12092048

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop