Search tips
Search criteria 


Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
Sci Rep. 2017; 7: 11460.
Published online 2017 September 13. doi:  10.1038/s41598-017-11602-5
PMCID: PMC5597609

Extremely high electrical conductance of microporous 3D graphene-like zeolite-templated carbon framework


We report the remarkably high electrical conductance of microporous 3D graphene-like carbons that were formed using lanthanum (La)-catalyzed synthesis in a Y zeolite (LaY) template investigated using conductive atomic force microscopy (C-AFM) and theoretical calculations. To uncover the relation between local electrical conductance and the microporous structures, we tuned the crystallographic ordering of LaY-templated carbon systems by changing the heating temperature. The structure of the LaY-templated carbon prepared at the higher temperature has graphene-like sp 2 hybridized bonds, which was confirmed using high-resolution transmission electron microscopy and X-ray diffraction measurements. C-AFM current–voltage spectroscopy revealed that the local current flow in the LaY-templated carbon depends on the quantity of C–C bonds within the narrow neck between the closed supercages (i.e. there are three types of carbon: carbon with heat treatment, carbon without heat treatment, and carbon synthesized at low temperature). The difference in electrical conductance on the LaY-templated carbon was also confirmed via theoretical computation using the Boltzmann transport theory and the deformation potential theory based on the density functional theory. These results suggest that the degree of order of the pores in the 3D zeolite-templated carbon structures is directly related to electrical conductance.


Ordered microporous carbons (i.e. pore size < 2 nm) synthesized using the nanocasting method with a zeolite as the sacrificial template are of scientific and technological interest owing to their regular and three-dimensional (3D) interconnected pore structure and high surface area greater than 2000 m2 g−1 110. The synthesis requires a carbon deposition process within the zeolite micropores that uses small hydrocarbon molecules (e.g. acetylene, ethylene, propylene)11. This process is often troublesome because pyrolytic carbonization of the hydrocarbons can also occur outside the zeolite template. Carbon that forms on the external surfaces of the zeolite hinders further diffusion of the hydrocarbon molecules into the micropores, and consequently causes insufficient formation of the carbon framework to replicate the zeolite pore structure. To mitigate this problem, many researchers have developed carbon deposition techniques, including a two-step carbon infiltration method involving the impregnation of liquid furfuryl alcohol, chemical vapor deposition (CVD) of propylene2, 3, 12, a CVD process at very low pressure, and a pulsed CVD method5, 11. Very recently, we reported a facile method to selectively form a carbon structure inside zeolite micropores. Certain cations (e.g. La3+, Y3+, Ca2+) catalytically affect carbon formation in the zeolite pores at sub-nanometer scale through d–π coordination bonding with small olefinic hydrocarbons. That is, using a zeolite containing these cations, the carbon structure can be formed at a lower temperature where pyrolytic carbonization of the hydrocarbon normally does not occur. In this way, a 3D graphene-like sp 2 carbon framework was formed only inside the zeolite pores without carbon deposition on the external zeolite surfaces13, 14.

Zeolite-templated carbons (ZTCs) are now considered to be highly promising as an electrode material for lithium batteries, supercapacitors, and fuel cells. In battery and supercapacitor applications, ZTC electrodes can exhibit high specific capacitance even under high current density conditions because of the 3D connected and ordered pore network, which promotes the fast transport of electrolyte ions. In fuel cell applications, the novel structure of ZTC is highly effective for loading a large amount of Pt catalyst with high dispersion, leading to high activity for the electrocatalytic oxygen reduction reaction. Distinct from the structural characteristics of ZTCs, electrical conductance is another critical physical property for high performance in such electrical applications. The ability to conduct an electric current influences the internal resistance of the devices.

The electrical properties of ordered microporous carbon formed using a zeolite Y template were previously investigated. T. Kyotani et al. reported that for a porous carbon formed using a zeolite Y with the CVD process at high reaction temperature (i.e. about 800 °C), the electronic conductivity is lower than that of carbon from bulk polyacrylonitrile because a graphitized structure with high conductivity is not formed because of the spatial limitation in the zeolite Y channels12. In a further study, K. Nueangnoraj et al. reported the formation of a crosslinked fullerene-like framework that is close to carbon Schwarzite (i.e. a 3D graphene-based framework) using a zeolite Y replica with heat treatment at 900 °C after the CVD process. They mentioned that the ordered framework is structurally diverse (i.e. is both semiconductive with a bandgap and metallic with zero bandgap)11. To better understand and develop the electrochemical performance of ZTCs, investigation of local electrical conductance is very important.

Conductive atomic force microscopy (C-AFM) is a technique that shows the local electrical properties of electrically heterogeneous surfaces at nanometer scale because of the precise and stable electrical contact between the sharp AFM tip and the surface15. C-AFM has been widely used to reveal the local electrical conductance of various nanostructures, including nanodots16, nanorods17, tetrapods18, and nanodumbells19. In a recent report, we showed the synthesis of highly ordered microporous carbons using zeolites functionalized by carbonization catalysts (e.g. La3+, Y3+ or Ca2+) and their electrical conductance measured using C-AFM13. The carbons obtained from faujasite (FAU) zeolite templates exhibited electrical conductance two orders of magnitude greater than that of amorphous mesoporous carbon, which is explained by the graphene-like framework of the ZTC consisting of sp 2-bonded carbons. In this study, we shed further light on the electrical properties of the ZTC (i.e. lanthanum-catalyzed synthesis in a FAU zeolite template carbon) by means of different combinations of synthesis, measurement, and theoretical calculation. The crystallinity of the ZTC was controlled by changing the heating conditions, which produced several model systems for studying their electrical properties. Variation of the conductance, dependent upon the pore ordering and topology of the ZTC, is systematically investigated by using samples with tailored microporous structures (e.g. carbon with and without heat treatment, or synthesized at a lower temperature). We compare the electrical conductance of ZTCs in this study to theoretical values that are calculated using a computational model structure that utilizes the Boltzmann transport theory and deformation potential theory based on density functional theory (DFT).


LaY zeolite-templated carbon

Three types of carbon materials were prepared through carbonization of ethylene within a La3+-exchanged Y (LaY) zeolite. For two carbon samples (C500-HT and C650-HT), carbon deposition was performed at 500 and 650 °C, respectively, followed by thermal treatment at 900 °C under dry N2 flow. The other sample (C650) was synthesized using ethylene CVD at 650 °C without the subsequent thermal treatment. Figure 1 shows X-ray diffraction (XRD) patterns and transmission electron microscopy (TEM) images for the final carbon products after template removal. Both the C500-HT and C650-HT samples exhibit a well-resolved XRD peak at 2θ = 6.4° (Fig. 1a). This XRD peak position is almost the same as that of the (111) diffraction of the LaY zeolite. Accordingly, the XRD patterns of C650-HT and C500-HT indicate that the zeolite pore structure is successfully replicated by the carbon framework. Faithful carbon replication is supported by TEM images (Fig. 1b,c), showing a highly ordered microporous structure with lattice fringes of 1.4 nm. The fringe spacing is consistent with the XRD peak at 2θ = 6.4°. Previous reports in the literature show that heat treatment is required to densify the deposited carbon framework in the zeolite and to improve pore regularity2022. On the other hand, the C650 sample obtained without thermal treatment exhibits a broad and weak XRD peak at around 2θ = 6.0° (Fig. 1a) and poorly-ordered pore arrays in the TEM image (Fig. 1d), which are probably caused by insufficient connectivity within the carbon frameworks without thermal treatment to retain the ordered porous structure after removing the template13, 5, 12. It is worth mentioning that none of the carbon samples have dense graphitic carbon layers outside the microporous carbon domains. The inhibition of carbon deposition on the external surfaces of the zeolite template can be attributed to the lower carbonization temperature13.

Figure 1
Structure of the LaY zeolite-templated carbon. (a) XRD patterns and (bd) TEM images of the carbon samples obtained using LaY zeolite as the template under different synthesis conditions: (b) C650-HT, (c) C500-HT, and (d) C650.

We investigated the carbon bonding in the samples using solid-state 13C magic-angle spinning NMR spectroscopy (Supplementary Fig. 1) and elemental analysis (Supplementary Table 1) (See Supplementary Information). All the NMR spectra are almost the same, showing predominant sp 2 carbon peaks at chemical shifts of 123 and 130 ppm. These peaks can be assigned to a 6-membered carbon ring and a distorted carbon ring (e.g. a 5- or 7-membered ring), respectively13. Thus, all the carbon samples are mainly composed of sp 2-hybridized carbon bonds. However, they have different hydrogen atom content, as evaluated by elemental analysis. The C/H atomic ratio of C650-HT is 8.9, while that of C500-HT is 5.7. The hydrogen content indicates the presence of edge-plane sites or vacancy defects in the graphene-like frameworks. That is, carbon deposition at temperatures that are too low can lead to formation of many fractures or vacant holes in the carbon framework even if the overall pore structure is well-ordered. On the other hand, the C650 sample exhibits a C/H ratio of 3.2, which can be related to a lack of framework connectivity, as expected from poor pore ordering. In addition, the Raman spectra exhibited no significant difference in the intensity ratio of the D- and G-bands or the band broadness for the carbon samples obtained using the different synthesis conditions because zeolite-templated carbons normally show a very broad intrinsic D-band that is attributed to the presence of five- or seven-membered carbon rings in the curved carbon structure (Supplementary Fig. 2)3, 13, 23, 24.

Electrical conductivity

To investigate the electrical conductance of the LaY-templated carbon, C-AFM was utilized because of its capability for local probing and for obtaining morphological and electrical properties under ambient conditions. Figure 2a shows the scheme for measuring the local electrical conductance of C650-HT on Au (111) (100 nm thickness) on a mica substrate using C-AFM with a Pt/Ir metal coated tip. The topography indicates randomly stacked LaY-templated carbons (bright parts) on the Au (111) surface after drop casting of carbon-dispersed isopropanol (Fig. 2b). The height and width of the stacked C650-HT are about 300 nm and 2 μm, respectively (Fig. 2c). Figure 2d shows representative current–voltage (I–V) curves measured on C650-HT with Au (111) as the reference. The Au (111) surface (3.5 × 10−4 Ω−1) shows a current flow only two times higher than that of C650-HT (1.7 × 10−4 Ω−1). Furthermore, we measured local I–V curves at the white square dots, moving in the direction of the red dotted arrow from left (first point) to right (fifteenth point), as marked in Fig. 2b. Positions 1–4 and 13–15 show the Au (111) surface and the other positions (i.e. 5–12) show C650-HT. The electrical conductance of C650-HT and the Au (111) surface seems indistinguishable because surface contamination of the Au (111) from the wet chemical processes and air can lead to varied and reduced local electrical conductance, as shown in Fig. 2e. Nevertheless, C650-HT consistently shows high electrical conductance.

Figure 2
Electrical conductance of LaY zeolite-templated carbon with heat treatment. (a) Scheme for probing the electrical conductance of the LaY-templated carbon on Au (111)/mica using conductive AFM with a Pt/Ir metal-coated tip. (b) Topography (2.5 × 6.8 µm ...

Figure 3a show the topography of C650 on a Au (111) substrate. The morphology and dispersity of the LaY-templated carbon on the Au (111) surface are also confirmed by scanning electron microscope (SEM) imaging, as shown in Supplementary Fig. 3 (See Supplementary Information). The line profile in Fig. 3b shows a height and width of ~600 nm and ~1.5 μm, respectively. The representative I–V curves of C650 and Au (111) are shown in Fig. 3d. The inset in Fig. 3d shows the magnified I–V curve of C650. C650 (2.4 × 10−7 Ω−1) exhibits an electrical conductance three orders of magnitude lower than that of the Au (111) surface (1.2 × 10−4 Ω−1). Moreover, Fig. 3c shows the local electrical conductance of C650 and Au (111), which were measured at the white square dots, moving in the direction of the red dotted arrow from left to right, as marked in Fig. 3a. Positions 1–7 and 13 are on the Au (111) surface and positions 8–12 are on C650. The difference in electrical conductance is quite clear between the C650 and Au (111). The topography of C650 on Au (111) was taken before and after the I–V curve measurements to compare the morphology and position of the carbon (See Supplementary Fig. 4 in the Supplementary Information).

Figure 3
Electrical conductance of the LaY zeolite-templated carbon without heat treatment. (a) Topography (4.2 × 4.2 µm2) of C650 on Au (111). (b) Height line profile along the white solid line in (a). (c) Local electrical ...

Figure 4 shows a comparison of the electrical conductance of the three different LaY-templated carbon samples and Au (111) as the reference. C650-HT clearly shows higher conductance than that of C650 and C500-HT, while C500-HT has a relatively high conductance range. The contact area between the probe and the surface can be obtained by using the Derjaguin–Muller–Toporov (DMT) model25, 26, given by A = π(R 2/3/K 2/3)(L + L c)2/3, where R is the curvature radius of the Pt/Ir-coated tip, K is the combined elastic modulus of the carbon (graphite) and the tip, L is the applied load, and L c is the adhesion. The calculated contact area A is 6.6 ± 0.6 nm2 using K = 530 GPa26, 27; R = 50 nm ± 10; L = 13 nN; and L c = 15.6 nN (C650-HT), 24.3 nN (C650), and 18.2 nN (C500-HT). The conductivity σ is 8.4 × 106 S m−1 and 2.0 × 104 S m−1 for C650-HT and C650, respectively.

Figure 4
Electrical conductance of the LaY-templated carbon as a function of synthesis temperature and heat treatment. For the measurement, an applied load of 13 nN was used. The adhesion (effective load) is 15.6 nN (28.6 nN), 24.3 nN (37.3 nN), 18.2 nN (31.2 ...


The main question here is: “What is the origin of the difference in electrical conductance between C650-HT and C650?”. The most likely possibility is the degree of pore order in the carbon product due to the role of heat treatment. If a poorly ordered, synthesized carbon/zeolite composite is heat treated at about 900 °C, it becomes a carbon product with highly ordered pores. TEM data also show highly well-ordered carbon structures after heat treatment, as shown in Fig. 1. The lack of pore order is due to inadequate C–C bonds within the narrow necks between the closed supercages. Lack of pore order can also be caused by using a lower temperature during the carbon synthesis process even though heat treatment is performed. Thus, we also confirmed the electrical conductance of C500-HT, as shown in Supplementary Fig. 5 (See Supplementary Information). Additionally, X-ray photoelectron spectroscopy (XPS) was employed to analyze the chemical components of the three types of LaY-templated carbon (i.e. C650-HT, C650, and C500-HT), as shown in Supplementary Fig. 6 (See Supplementary Information). The C1s core-level spectrum revealed three major peaks, which are assigned to C–C/C–H bonds at 284.4 eV; the C–OH alcohol group at 286.1 eV; and the O = C–OH/O–C = O carbonate group at 288.6 eV. Comparing C650 and C500-HT, the portion of C–C/C–H was lower and C–OH was higher in C500-HT. The O1 s spectrum revealed two major peaks, which are assigned to C–O bonds at 531.7 eV and C=O bonds (i.e. of the alcohol group) at 532.9 eV.

To computationally evaluate the electrical conductivity of the well-ordered porous 3D graphene-like zeolite-templated carbon framework, we used a Schwarzite structure28 as the model shown in Fig. 5a. The unit cell contains a core with O h symmetry and negative local Gaussian curvature; six (4,4) carbon nanotube segments connect the cores located in the neighboring unit cells. Our electronic structure calculation based on DFT revealed that the Schwarzite structure is metallic28, as shown in Fig. 5b. With the electronic band structure, its conductivity was calculated using the Boltzmann Transport (BT)29 and deformation potential (DP) theories30. The conductivity tensor σαβ() for a given energy can be calculated using29


where N is the number of k-points sampled, n, k is the summation over the band index n in the momentum k-space, and δ(x) is Dirac delta function. The conductivity tensor for a given n and k can be obtained using σαβ(nk) = e2τ(nk) vα(nk)vβ(nk). Here, τ is the relaxation time and vα is the α-component of the group velocity v(nk) = ∇k(nk), where (nk) is the electronic band energy shown in Fig. 5b, and k is the gradient in the momentum space.

Figure 5
Electronic structure calculation. (a) Equilibrium geometry of a Schwarzite structure with 200 carbon atoms per cubic unit cell with a lattice constant a; model of the well-ordered pore 3D graphene-like zeolite-templated carbon framework. (b) Electronic ...

To determine the electric conductivity σ, it is important to calculate the relaxation time τ, which is difficult. Usually, because of such difficulty, the BT theory is used under a constant relaxation time approximation to get σ/τ rather than σ. However, in this study, we calculated the relaxation time τ using the DP theory30 along with the effective mass approximation, which provides the following relation


where μ is the carrier mobility, m* is the effective mass at the band edge, and e is the electronic charge. The elastic constant was obtained using C = [∂2Et/∂(ΔL/L0)2]/V0, that is, the second derivative of the total energy E t of the system with respect to the strain ΔL/l0 along the transport direction. Here, l 0 is the equilibrium length and V 0 is the cell volume. From DP, E DP was determined using Δ/(Δl/l0) or the ratio of the strain-induced energy level shifts Δ at the band edges to the strain. Δ was evaluated with respect to the 1 s core level obtained using an all-electron calculation. In this study, only one transport direction was considered because of the simple cubic symmetry of our model system. Note that the DP theory was originally developed for hole and electron carriers located at the valence band maximum and at the conduction band minimum, respectively, in semiconductors. To apply the theory to our metallic system, two doping situations were considered. In each doping case, the Fermi level is located at either of the energy levels at the band edges specified with a blue or red dashed line in Fig. 5b, which can be achieved using either a positive or negative doping process, respectively. The elastic constant C, the effective mass m*, and the DP value E DP were estimated by fitting various data obtained by applying the DP theory, as shown in Supplementary Fig. 7 (See Supplementary Information). The relaxation time τ and the mobility μ were then evaluated using Eq. 2 with these estimated parameters.

Figure 5c,d shows the electrical conductivity of our model system (i.e. Schwarzite) as a function of the doping levels that we calculated with all the necessary quantities evaluated, as mentioned above. In the positive (negative) doping case that resulted in the Fermi level of the system shifting down (up), the conductivity was calculated to be σ = 2~3 × 106 S m− 1 (~1 × 106 S m− 1), as shown in Fig. 5c (Fig. 5d). To verify our method, we calculated the relaxation time τ of graphene and thus its electrical conductivity σ using the same method. They were calculated to be τ = ~ 24 ps, and σ = 5 × 107 S m− 1, which are in good agreement with experimental results31, 32. Our calculated conductivity for the model system is about an order of magnitude smaller than that of graphene, which is reasonably acceptable because of the structural differences between the two. Interestingly, the conductivity of the model system is comparable to that of normal metals. Even more interestingly, it is several orders of magnitude larger than that of amorphous carbon (σ = ~ 2 × 103 S m− 1)33. These theoretical results agree well with our experimental observations, though there are deviations of conductivity from the multipath through the carbon chains. Additionally, we measured the sheet resistance Rs of the zeolite-carbon composites that were deposited using the drop casting method on a printed circuit board substrate with a stripe patterned metal electrode (interval distance of 1 mm) using a probe station. Relatively high Rs of 67.5 kΩ/sq and infinite (higher than 100 MΩ/sq) were measured on the carbon with and without heat treatment (at 650 °C), respectively. We found that it was difficult to form a uniform film when depositing zeolite-carbon composites because particle formation occurs, rather than the geometrically uniform thin film that forms after drop casting. In addition, we found that it was difficult to obtain the absolute sheet resistance for the powder samples because there are many factors that influence the resistance (e.g. high contact resistance, low density and uniformity, and packing density), which makes the reliable measurement of electrical conductance on carbon samples with a macroscopic four-probe configuration quite challenging. In the case of C-AFM, we can make a reliable electrical contact with carbon particles with the AFM tip. Therefore, we think that this study shows a unique capability for nanoscale electrical measurements on carbon particles with C-AFM.

In conclusion, we demonstrate that the electrical conductance of LaY-templated carbon depends on the heat treatment of the carbon structure, which increases the C–C bonds within the narrow necks between neighboring supercages. The higher conductance of the LaY-templated carbon with heat treatment than that of the LaY-templated carbon without heat treatment or low synthesis temperature is revealed because of the loss of pore order and differing hydrogen content. The crystallinity of these microporous carbon products is confirmed using HR-TEM and XRD measurements, and the results are well matched with the conductance trends. In addition, the electrical conductivity of the LaY-templated carbon was computed using the Boltzmann transport theory and the deformation potential theory based on density functional theory. These results suggest that the degree of order for the pores in the 3D ZTC structures is related to the electrical conductivity and the LaY-templated carbon is potentially applicable in supercapacitors and rechargeable batteries.


Synthesis of LaY-templated carbon

Y zeolite was synthesized according to methods found in the literature34, 35. Ion exchange was performed with a 0.5 M aqueous solution of lanthanum nitrate hexahydrate (98%, Sigma-Aldrich) at room temperature. For carbon deposition within the La3+ form of the Y zeolite (i.e. LaY zeolite), 0.5 g of LaY zeolite was placed in a quartz plug-flow reactor equipped with a 15-mm diameter fritted disk. The reactor was then heated to different temperatures (i.e. 500 °C and 650 °C) with a heating rate of 5 °C min−1 under dry N2 flow. After stabilizing the temperature at the desired value, a mixture of ethylene, N2, and water vapor (10/87/3) was passed through the sample bed. The flow rate of the gas mixture was 130 cm3 min−1. Injection of the ethylene gas and water vapor was stopped when carbon deposition was complete. The carbon deposition process gave a carbon yield of approximately 0.34 g per 1 g zeolite, which was confirmed using thermogravimetric analysis. After carbon deposition, the reactor was heated to 900 °C under N2 flow, and this temperature was maintained for 2 h to rigidify the carbon frameworks formed in the LaY zeolite. After cooling the reactor, the resultant sample was repeatedly treated in a 0.3 M HF/0.15 M HCl solution to remove the LaY zeolite. The template-free carbon products were collected using filtration, washed with deionized water, and dried at 100 °C. Complete removal of the zeolite template was confirmed by thermogravimetry using a TGA Q50 (TA Instruments) (Supplementary Fig. 8). The final carbon samples were denoted as C500-HT and C650-HT, according to the carbon deposition temperature. For comparison, a carbon sample was also synthesized through carbon deposition at 650 °C without the subsequent thermal treatment at 900 °C. The carbon product, thus obtained, is denoted as C650. The synthesized LaY-templated carbon was dispersed in isopropanol for drop casting on a Au (111) substrate. After drop casting, the LaY-templated carbon was randomly dispersed and stacked on the Au (111) substrate. The Au (111) film was prepared using flame annealing after e-beam evaporation of a 100-nm thick Au layer on a mica substrate.

Characterization of LaY-templated carbon

Characterization of the LaY-templated carbon was performed using high-resolution transmission electron microscopy (HR-TEM), XRD, SEM, and XPS. The crystallinity of the LaY-templated carbon was determined using HR-TEM (FEI Titan E-TEM G2) with 300 kV acceleration using a holey carbon grid (300 mesh) after ethanol dispersion on the support. SEM images were obtained using a Magellan 400 (FEI company) at a bias voltage of 2 kV and work distance of 4.1 mm. XPS data for the LaY-templated carbon on Au (111) were taken using a Thermo VG Scientific Sigma Probe spectrometer equipped with an Al-Kα X-ray source (1486.3 eV). Powder XRD patterns were recorded on a Rigaku Multiplex instrument using Cu Kα radiation (30 kV, 40 mA).

AFM measurements

The electrical conductance and morphology of the LaY-templated carbon on Au (111) substrate were characterized using an environmental AFM (5500, Agilent) in combined conductance sensing mode with a Pt/Ir-coated AFM probe (PPP-EFM-50, Nanosensors) in air. A crushed ZTC sample was measured using C-AFM. The AFM tip has a typical force constant of 2.8 N/m and resonance frequency of 75 kHz. The topographical images were taken using non-contact mode AFM, then the I–V curves of the local positions were obtained in I–V spectroscopy mode. We confirm that the tip–sample contact is in the elastic regime by using a series of reproducible measurements of the I–V curves that indicate the absence of plastic deformation of the tip and the sample.

Theoretical calculations

A first-principles calculation based on DFT was performed using the DMol3 package36 to study the electronic structures and electrical conductance of the schwarzite structure28 selected as the model for the well-ordered pore 3D graphene-like zeolite-templated carbon framework. We used the generalized gradient approximation with the Perdew–Burke–Ernzerhof 37 functional to describe the exchange-correlations functional. A double numerical basis set was chosen with a real-space cutoff of 3.7 Å to expand the electronic wave functions in the all-electron calculation. The octupole scheme for the multipolar fitting procedure and a fine grid scheme for numerical integration were employed to accurately evaluate the charge density. We performed structure relaxation using the Broyden–Fletcher–Goldfarb–Shanno algorithm3843 without any symmetry constraints until all the atomic forces became smaller than 4 × 10−4 Ha Å−1. For such structural relaxation, only the Г point was used for the Brillouin zone integration since the model schwarzite structure is very large. The electronic transport properties were calculated using the semi-classical BT theory29 with the relaxation time obtained using the DP theory30. To evaluate reliable transport coefficients, a denser 2 × 2 × 2 k-grid was used for the Monkhorst–Pack scheme44.

Electronic supplementary material


This work was supported by the Institute for Basic Science (IBS) [IBS-R004-A2-2017-a00 and IBS-R004-D1]. Y.-K.K. and S.-H.K. were supported by the Ministry of Trade, Industry & Energy (MOTIE) of Korea (Project No. 10045360) and the National Research Foundation of Korea (NRF-2015R1A2A2A01006204). Some portions of our computational work were carried out using resources at the KISTI Supercomputing Center (KSC-2016-C3-0034).

Author Contributions

Author Contributions

H.L., K.K., R.R. and J.Y.P. designed the experiments. H.L. and J.H.K. performed the AFM experiments, K.K. and Y.K. performed the synthesis and characterization experiments, S.-H.K. and Y.-K.K. performed the theoretical analysis. All authors discussed the results and wrote the manuscript.


Competing Interests

The authors declare that they have no competing interests.


Hyunsoo Lee and Kyoungsoo Kim contributed equally to this work.

Electronic supplementary material

Supplementary information accompanies this paper at doi:10.1038/s41598-017-11602-5

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Contributor Information

Young-Kyun Kwon,

Ryong Ryoo,

Jeong Young Park,


1. Ma Z, Kyotani T, Tomita A. Preparation of a high surface area microporous carbon having the structural regularity of Y zeolite. Chem. Commun. 2000;2:2365–2366. doi: 10.1039/b006295m. [Cross Ref]
2. Ma Z, Kyotani T, Liu Z, Terasaki O, Tomita A. Very high surface area microporous carbon with a three-dimensional nano-array structure: synthesis and its molecular structure. Chem. Mater. 2001;13:4413–4415. doi: 10.1021/cm010730l. [Cross Ref]
3. Nishihara H, et al. A possible buckybowl-like structure of zeolite templated carbon. Carbon. 2009;47:1220–1230. doi: 10.1016/j.carbon.2008.12.040. [Cross Ref]
4. Parmentier J, Gaslain FOM, Ersen O, Centeno TA, Solovyov LA. Structure and sorption properties of a zeolite-templated carbon with the EMT structure type. Langmuir. 2014;30:297–307. doi: 10.1021/la402762v. [PubMed] [Cross Ref]
5. Kyotani T, Ma Z, Tomita A. Template synthesis of novel porous carbons using various types of zeolites. Carbon. 2003;41:1451–1459. doi: 10.1016/S0008-6223(03)00090-3. [Cross Ref]
6. Yang Z, Xia Y, Mokaya R. Enhanced hydrogen storage capacity of high surface area zeolite-like carbon materials. J. Am. Chem. Soc. 2007;129:1673–1679. doi: 10.1021/ja067149g. [PubMed] [Cross Ref]
7. Yang Z, Xia Y, Sun X, Mokaya R. Preparation and hydrogen storage properties of zeolite-templated carbon materials nanocast via chemical vapor deposition: effect of the zeolite template and nitrogen doping. J. Phys. Chem. B. 2006;110:18424–18431. doi: 10.1021/jp0639849. [PubMed] [Cross Ref]
8. Zhou J, et al. Effect of cation nature of zeolite on carbon replicas and their electrochemical capacitance. Electrochim. Acta. 2013;89:763–770. doi: 10.1016/j.electacta.2012.11.068. [Cross Ref]
9. Johnson SA, Brigham ES, Ollivier PJ, Mallouk TE. Effect of micropore topology on the structure and properties of zeolite polymer replicas. Chem. Mater. 1997;9:2448–2458. doi: 10.1021/cm9703278. [Cross Ref]
10. Corma A. From microporous to mesoporous molecular sieve materials and their use in catalysis. Chem. Rev. 1997;97:2373–2420. doi: 10.1021/cr960406n. [PubMed] [Cross Ref]
11. Nueangnoraj K, et al. Formation of crosslinked-fullerene-like framework as negative replica of zeolite Y. Carbon. 2013;62:455–464. doi: 10.1016/j.carbon.2013.06.033. [Cross Ref]
12. Kyotani T, Nagai T, Inoue S, Tomita A. Formation of new type of porous carbon by carbonization in zeolite nanochannels. Chem. Mater. 1997;9:609–615. doi: 10.1021/cm960430h. [Cross Ref]
13. Kim K, et al. Lanthanum-catalysed synthesis of microporous 3D graphene-like carbons in a zeolite template. Nature. 2016;535:131–135. doi: 10.1038/nature18284. [PubMed] [Cross Ref]
14. Kim K, Kwon Y, Lee T, Cho SJ, Ryoo R. Facile large-scale synthesis of three-dimensional graphene-like ordered microporous carbon via ethylene carbonization in CaX zeolite template. Carbon. 2017;118:517–523. doi: 10.1016/j.carbon.2017.03.082. [Cross Ref]
15. Park JY, Maier S, Hendriksen B, Salmeron M. Sensing current and forces with SPM. Mater. Today. 2010;13:38–45. doi: 10.1016/S1369-7021(10)70185-1. [Cross Ref]
16. Lee H, Lee YK, Van TN, Park JY. Nanoscale Schottky behavior of Au islands on TiO2 probed with conductive atomic force microscopy. Appl. Phys. Lett. 2013;103:173103. doi: 10.1063/1.4826140. [Cross Ref]
17. Park WI, Yi GC, Kim J-W, Park S-M. Schottky nanocontacts on ZnO nanorod arrays. Appl. Phys. Lett. 2003;82:4358–4360. doi: 10.1063/1.1584089. [Cross Ref]
18. Fang L, et al. Mechanical and electrical properties of CdTe tetrapods studied by atomic force microscopy. J. Chem. Phys. 2007;127:184704. doi: 10.1063/1.2786993. [PubMed] [Cross Ref]
19. Kwon S, et al. Probing the nanoscale Schottky barrier of metal/semiconductor interfaces of Pt/CdSe/Pt nanodumbbells by conductive-probe atomic force microscopy. Nanoscale. 2015;7:12297–12301. doi: 10.1039/C5NR02285A. [PubMed] [Cross Ref]
20. Ma Z, Kyotani T, Tomita A. Synthesis methods for preparing microporous carbons with a structural regularity of zeolite Y. Carbon. 2002;40:2367–2374. doi: 10.1016/S0008-6223(02)00120-3. [Cross Ref]
21. Hou P-X, Yamazaki T, Orikasa H, Kyotani T. An easy method for the synthesis of ordered microporous carbons by the template technique. Carbon. 2005;43:2624–2627. doi: 10.1016/j.carbon.2005.05.001. [Cross Ref]
22. Choi S, et al. Large-scale synthesis of high-quality zeolite-templated carbons without depositing external carbon layers. Chem. Eng. J. 2015;280:597–605. doi: 10.1016/j.cej.2015.06.055. [Cross Ref]
23. Leyva-García S, et al. Characterization of a zeolite-templated carbon by electrochemical quartz crystal microbalance and in situ Raman spectroscopy. Carbon. 2015;89:63–73. doi: 10.1016/j.carbon.2015.03.016. [Cross Ref]
24. Biener J, et al. Macroscopic 3D nanographene with dynamically tunable bulk properties. Adv. Mater. 2012;24:5083–5087. doi: 10.1002/adma.201202289. [PubMed] [Cross Ref]
25. Derjaguin BV, Muller VM, Toporov YP. Effect of contact deformations on the adhesion of particles. J. Colloid Interface Sci. 1975;53:314–326. doi: 10.1016/0021-9797(75)90018-1. [Cross Ref]
26. Park JY, et al. Tribological properties of quasicrystals: effect of aperiodic versus periodic surface order. Phys. Rev. B. 2006;74:024203. doi: 10.1103/PhysRevB.74.024203. [Cross Ref]
27. Blakslee OL, Proctor DG, Seldin EJ, Spence GB, Weng T. Elastic constants of compression-annealed pyrolytic graphite. J. Appl. Phys. 1970;41:3373–3382. doi: 10.1063/1.1659428. [Cross Ref]
28. Park S, Kittimanapun K, Ahn JS, Kwon Y-K, Tománek D. Designing rigid carbon foams. J. Phys. Condens. Matter. 2010;22:334220. doi: 10.1088/0953-8984/22/33/334220. [PubMed] [Cross Ref]
29. Madsen GKH, Singh DJ. BoltzTraP. A code for calculating band-structure dependent quantities. Comput. Phys. Commun. 2006;175:67–71. doi: 10.1016/j.cpc.2006.03.007. [Cross Ref]
30. Bardeen J, Shockley W. Deformation potentials and mobilities in non-polar crystals. Phys. Rev. 1950;80:72–80. doi: 10.1103/PhysRev.80.72. [Cross Ref]
31. Neugebauer P, Orlita M, Faugeras C, Barra A-L, Potemski M. How perfect can graphene be? Phys. Rev. Lett. 2009;103:136403. doi: 10.1103/PhysRevLett.103.136403. [PubMed] [Cross Ref]
32. Novoselov KS, et al. Two-dimensional gas of massless Dirac fermions in graphene. Nature. 2005;438:197–200. doi: 10.1038/nature04233. [PubMed] [Cross Ref]
33. Pauleau, Y. & Barna, P. B. Protective coatings and thin films: synthesis, characterization, and applications, Springer, p. 215 (1997).
34. Delprato F, Delmotte L, Guth JL, Huve L. Synthesis of new silica-rich cubic and hexagonal faujasites using crown-etherbased supramolecules as templates. Zeolites. 1990;10:546–552. doi: 10.1016/S0144-2449(05)80310-0. [Cross Ref]
35. Ferchiche S, Warzywoda JAS., Jr. Direct synthesis of zeolite Y with large particle size. Int. J. Inorg. Mater. 2001;3:773–780. doi: 10.1016/S1466-6049(01)00046-0. [Cross Ref]
36. Delley B. An all‐electron numerical method for solving the local density functional for polyatomic molecules. J. Chem. Phys. 1990;92:508–517. doi: 10.1063/1.458452. [Cross Ref]
37. Perdew JP, Burke K, Ernzerhof M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996;77:3865–3868. doi: 10.1103/PhysRevLett.77.3865. [PubMed] [Cross Ref]
38. Broyden CG. The convergence of a class of double-rank Minimization Algorithms 1. General Considerations. IMA J. Appl. Math. 1970;6:76–90. doi: 10.1093/imamat/6.1.76. [Cross Ref]
39. Fletcher R. A new approach to variable metric algorithms. Comput. J. 1970;13:317–322. doi: 10.1093/comjnl/13.3.317. [Cross Ref]
40. Fletcher, R. Practical methods of optimization; (2nd ed.). Wiley-Interscience, (1987).
41. Goldfarb D. A family of variable-metric methods derived by variational means. Math. Comp. 1970;24:23–26. doi: 10.1090/S0025-5718-1970-0258249-6. [Cross Ref]
42. Shanno DF. Conditioning of quasi-newton methods for function minimization. Math. Comp. 1970;24:647–656. doi: 10.1090/S0025-5718-1970-0274029-X. [Cross Ref]
43. Shanno DF, Kettler PC. Optimal conditioning of quasi-newton methods. Math. Comp. 1970;24:657–664. doi: 10.1090/S0025-5718-1970-0274030-6. [Cross Ref]
44. Monkhorst HJ, Pack JD. Special points for Brillouin-zone integrations. Phys. Rev. B. 1976;13:5188–5192. doi: 10.1103/PhysRevB.13.5188. [Cross Ref]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group