CN104658941A - Method for measuring parameters of chip lamination structure - Google Patents

Method for measuring parameters of chip lamination structure Download PDF

Info

Publication number
CN104658941A
CN104658941A CN201510089321.8A CN201510089321A CN104658941A CN 104658941 A CN104658941 A CN 104658941A CN 201510089321 A CN201510089321 A CN 201510089321A CN 104658941 A CN104658941 A CN 104658941A
Authority
CN
China
Prior art keywords
beta
semiconductor layer
formula
chip
insulating barrier
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201510089321.8A
Other languages
Chinese (zh)
Other versions
CN104658941B (en
Inventor
郭宇锋
李曼
张长春
吉新村
夏晓娟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Posts and Telecommunications Nantong Institute Limited
Original Assignee
Nanjing Post and Telecommunication University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing Post and Telecommunication University filed Critical Nanjing Post and Telecommunication University
Priority to CN201510089321.8A priority Critical patent/CN104658941B/en
Publication of CN104658941A publication Critical patent/CN104658941A/en
Application granted granted Critical
Publication of CN104658941B publication Critical patent/CN104658941B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The invention provides a method for testing parameters of a chip lamination structure. According to the method, a high-frequency capacitance-voltage curve and a low-frequency capacitance-voltage curve of the chip lamination structure are tested by the aid of a capacitance-voltage tester, the insulation layer thickness and the semiconductor doping concentration of the structure are measured respectively according to a low-frequency saturation value and a high-frequency saturation value, and the insulation layer fixed charge density of the structure is measured according to a point with zero externally-applied voltage in the low-frequency capacitance-voltage curve as well as a characteristic extreme value point. A simple and non-destructive characterization method is provided for evaluation of the reliability of the chip lamination structure in three-dimensional integration.

Description

A kind of method of measurement of chip-stack structural parameters
Technical field
The present invention relates to field of semiconductor device test, particularly relate to the method for testing of thickness of insulating layer in chip-stack structure, doping content of semiconductor and insulating barrier fixed charge density.
Background technology
Three-dimensional is integrated is the final solution overcoming " More Moore " application, raising packaging density and circuit working speed and realize the multi-functional challenge of integrated circuit.Chip-stacked technology realizes one of three-dimensional integrated key technology, and to realize chip-stacked basic skills be characteristics of Direct Wafer Bonded.In order to assess bonding quality, although the cross-sectional analysis that often adopts of current people and bond strength detection method directly perceived, namely the chip after detecting goes to pot, and be a kind of destructively to detect, its application is restricted.Although infrared thermography is non-damaged data, only can detects bonding cavity, and more efficiently information can not be obtained.Therefore, one can be found neither to destroy three-dimensional stacking structure, the non-damaged data mode of effective evaluation can be carried out in para-linkage interface again, significant.
It is three-dimensional that integrated basic structure---chip-stack structure is divided into five layers, as shown in Figure 1, and the first semiconductor layer 1, second semiconductor layer 5, first insulating barrier 2, second insulating barrier 4, metal level 3.Wherein, the doping type of the first semiconductor layer 1, second semiconductor layer 5 is identical, namely can be P type, also can be N-type.
Current, non-damaged data mode conventional in semiconductor technology is C-V characteristic method, and low frequency and High Frequency C-V curve are widely used and extract the isostructural physical parameter of MIS, MOS.Tang Yi, Chinese patent, 200710046681.5, utilize charge pump method of testing to extract metal-oxide-semiconductor interfacial state.Duan little Jin, Chinese patent, ZL200710120481.X, is integrated with photocarrier radiometric technique and free-carrier Absorption measuring technique, provides a kind of method measuring doping content of semiconductor.
The invention provides a kind of method of nondestructive characterization chip-stack structure, i.e. low-and high-frequency C-V method, the thickness of insulating layer of chip-stack structure, doping content of semiconductor and insulating barrier fixed charge density can be recorded based on the method.
Summary of the invention
Goal of the invention: the object of the invention is to overcome the deficiencies in the prior art, provides a kind of method of testing of chip-stack structural parameters, can the thickness of insulating layer of this structure of Quick Measurement, doping content of semiconductor and insulating barrier fixed charge density.
Technical scheme: the method for testing of extraction chip-stack structural parameters of the present invention, it comprises the following steps:
Step 1: the High Frequency C-V curve and the low frequency C-V curve that record chip-stack structure with electric capacity-voltage tester and C-V tester;
Step 2: read electric capacity saturation value from low frequency C-V curve, utilize formula 1 to record the thickness of insulating layer of chip-stack structure first insulating barrier (2), the second insulating barrier (4);
t i > > e i 2 C LF max Formula 1
Wherein, t ifor thickness of insulating layer, ε ithe dielectric constant of insulating barrier, C lFmaxit is low frequency capacitive saturation value;
Step 3: bear semiaxis from High Frequency C-V curve and read electric capacity saturation value, utilize formula 2 to record the doping content of semiconductor of chip-stack structure first semiconductor layer (1);
C HF min > > e i 2 t i + e i e s 2 qe s kT q 2 N A ln ( N A n i ) Formula 2
Wherein, C hFminhigh frequency capacitance saturation value, ε sbe the dielectric constant of semiconductor layer, k is Boltzmann constant, and T is absolute temperature, and q is electron charge.θ is transoid modifying factor, and its value is 2.125, n iintrinsic carrier concentration, N asemiconductor layer doped concentration;
Step 4: read electric capacity saturation value from High Frequency C-V curve positive axis, utilize formula 2 to record the doping content of semiconductor of chip-stack structure second semiconductor layer (5);
Step 5: capacitance when being zero from high frequency or low frequency C-V curve reading applied voltage, sets up formula 3 and formula 4;
1 C 0 = 2 C i + 1 C s 10 + 1 C s 20 Formula 3
( V s 10 + V s 20 ) - 2 V fermi = Q s 10 + Q s 20 C i + Q f 1 + Q f 2 C i Formula 4
Wherein, C 0for capacitance when applied voltage is zero, C ibe capacitive dielectric layer, its value is: C i>>e i/ t i.C s10with C s20be respectively lower semiconductor layer differential capacitance, they are about V s10, V s20function, V s10with V s20for applied voltage be zero time, upper and lower semiconductor layer surface gesture;
C s 10 = ϵ s L D - exp ( - βV s 10 ) + 1 + n p 0 p p 0 ( exp ( βV s 10 ) - 1 ) F ( βV s 10 , n p 0 / p p 0 )
C s 20 = ϵ s L D - exp ( - βV s 20 ) + 1 + n p 0 p p 0 ( exp ( βV s 20 ) - 1 ) F ( βV s 20 , n p 0 / p p 0 )
Wherein, n p0and p p0equilibrium electron concentration and the hole concentration of P-type silicon respectively, L dbe the Debye length of semiconductor layer, the expression formula of F (a, b) is: q f1, Q f2be the fixed charge density of the first insulating barrier (2), the second insulating barrier (4), V fermifor Fermi potential.Q s10with Q s20be respectively the charge density of the first semiconductor layer (1), the second semiconductor layer (5), they are about V s10, V s20function;
Q s 10 = - sign ( V s 10 ) 2 ϵ s βL D F ( βV s 10 , n p 0 / p p 0 ) / q
Q s 20 = - sign ( V s 20 ) 2 ϵ s βL D F ( βV s 20 , n p 0 / p p 0 ) / q
Step 6: read feature extreme point (V from high frequency or low frequency C-V curve g *, C *), set up formula 5 and formula 6;
Q f 1 + Q f 2 = - 4 ϵ s βL D F ( βV s * , n p 0 / p p 0 ) - Q m 0 Formula 5
Q f 1 - Q f 2 = V g * C i Formula 6
Wherein and Q m0tried to achieve by following formula:
ϵ s L D - exp ( - βV s * ) + 1 + n p 0 p p 0 ( exp ( βV s * ) - 1 ) F ( βV s * , n p 0 / p p 0 ) = 2 C * C i C i - 2 C *
Q m0=C i(2V fermi-V s01-V s02);
Step 7: simultaneous formula 3,4,5,6, utilizes Program to record the fixed charge density of the first insulating barrier (2), the second insulating barrier (4).
The dopant material of the first semiconductor layer (1), the second semiconductor layer (5) is Si, Ge, GaAs or InP.
The material of the first insulating barrier (2), the second insulating barrier (4) is SiO 2, SiN or Si 3n 4dielectric material.
The material of metal level (3) is Cu, TiN or Ti.
The doping type of the first semiconductor layer (1), the second semiconductor layer (5) is identical, is namely P type or N-type.
The doping content of the first semiconductor layer (1), the second semiconductor layer (5) is identical or not identical.
Fixed charge density in first insulating barrier (2), the second insulating barrier (4) is identical or not identical.
Beneficial effect:
1) the present invention is a kind of method of nondestructive characterization chip-stack structure, i.e. low-and high-frequency C-V method.
2) the present invention can utilize the low frequency of chip-stack structure and High Frequency C-V curve to measure the many kinds of parameters of this structure, low frequency saturation value and high frequency saturation value can be utilized to measure thickness of insulating layer and the doping content of semiconductor of this structure, utilize low frequency C-V extra curvature making alive be zero point and feature extreme point measure the insulating barrier fixed charge density of this structure.
3) the present invention goes for dissimilar C-V curve.When applied voltage is zero, according to the surface state of semiconductor layer 1 and 5, low frequency and the High Frequency C-V curve of chip-stack structure are divided into Four types: I, II, III, IV, respectively as shown in Figure 2-5.The surface state of its correspondence is as shown in table 1.
Table 1 zero partially time the surface state of layer 1 and 5 and chip-stack structure C-V curve type relation table
4) the auto-adaptive parameter method of measurement not relying on concrete device model of the present invention is when carrying out parameter measurement for new device model, the relative actual value error of measured value is all below 4%, both consistency are good, and therefore the present invention effectively can measure the structural parameters of this structure.Method of measurement of the present invention decreases carries out the time spent by computing to chip-stack structural model formula in earlier stage, substantially reduces the cycle that this model parameter is measured, simple and quick.
Accompanying drawing explanation
Fig. 1 is chip-stack structural representation;
Fig. 2 is chip-stack structure low frequency provided by the invention and High Frequency C-V curve type I;
Fig. 3 is chip-stack structure low frequency provided by the invention and High Frequency C-V curve type II;
Fig. 4 is chip-stack structure low frequency provided by the invention and High Frequency C-V curve type III;
Fig. 5 is chip-stack structure low frequency provided by the invention and High Frequency C-V curve type IV.
Wherein have: the first semiconductor layer 1, second semiconductor layer 5, first insulating barrier 2, second insulating barrier 4, metal level 3.
Embodiment
Further illustrate the present invention below by embodiment, but and limit the scope of the invention never in any form.
In the present embodiment, the substrate silicon layer of sample is the doping of P type, and its thickness is 725 μm, insulating barrier SiO 2thickness be 500nm, metal level adopts TiN+Ti+TiN laminated construction, and wherein, the thickness of two-layer TiN is the thickness of 60nm, Ti is 80nm.The work function of TiN and Ti is respectively 4.7eV and 4.33eV.
Be positioned over by sample on test platform, chassis material is copper, and top pressure probe material is tungsten, is calculated the impedance of sample, and then obtain low frequency capacitive value and high-frequency electrical capacitance by the alternating voltage between record pressure probe and chassis and alternating current.Substrate ground during experiment, the DC offset voltage on pressure probe scans 40V from-40V, scanning 0.1V per second.The frequency of low frequency small-signal is 20HZ, and amplitude is 20mV.The frequency of high frequency small-signal is 2MHz, and amplitude is 20mV.
In order to the change making the low frequency of chip-stack structure and High Frequency C-V curve present Fig. 2 ~ 5, the device size that the present invention selects comprises (thickness of insulating layer t i, doping content of semiconductor N a, insulating barrier fixed charge density Q f1, Q f2) 0.3 μm, 3 × 10 14cm -3,-9 × 10 11cm -3v -1,-1 × 10 11cm -3v -1; 0.5 μm, 5 × 10 14cm -3, 2 × 10 11cm -3v -1, 8 × 10 11cm -3v -1; 0.7 μm, 7 × 10 14cm -3, 4 × 10 11cm -3v -1, 10 × 10 11cm -3v -1with 0.9 μm, 9 × 10 14cm -3, 7 × 10 11cm -3v -1, 14 × 10 11cm -3v -1deng yardstick.
Measure all devices at different applied voltage V gunder low frequency and high frequency C-V characteristic, calculate thickness of insulating layer, doping content of semiconductor and insulating barrier fixed charge density according to technique scheme computing formula.Errors is as shown in table 2 ~ 4.
Table 2 thickness of insulating layer actual value and the relative error with measured value thereof
Table 3 doping content of semiconductor actual value and the relative error with measured value thereof
Table 4 fixed charge density actual value and the relative error with measured value thereof

Claims (7)

1. a method of measurement for chip-stack structural parameters, is characterized in that: the method comprises the following steps:
Step 1: the High Frequency C-V curve and the low frequency C-V curve that record chip-stack structure with electric capacity-voltage tester and C-V tester;
Step 2: read electric capacity saturation value from low frequency C-V curve, utilize formula 1 to record the thickness of insulating layer of chip-stack structure first insulating barrier (2), the second insulating barrier (4);
t i > > e i 2 C LF max Formula 1
Wherein, t ifor thickness of insulating layer, ε ithe dielectric constant of insulating barrier, C lFmaxit is low frequency capacitive saturation value;
Step 3: bear semiaxis from High Frequency C-V curve and read electric capacity saturation value, utilize formula 2 to record the doping content of semiconductor of chip-stack structure first semiconductor layer (1);
C HF min > > e i 2 t i + e i e s 2 qe s kT q 2 N A ln ( N A n i ) Formula 2
Wherein, C hFminhigh frequency capacitance saturation value, ε sbe the dielectric constant of semiconductor layer, k is Boltzmann constant, and T is absolute temperature, and q is electron charge.θ is transoid modifying factor, and its value is 2.125, n iintrinsic carrier concentration, N asemiconductor layer doped concentration;
Step 4: read electric capacity saturation value from High Frequency C-V curve positive axis, utilize formula 2 to record the doping content of semiconductor of chip-stack structure second semiconductor layer (5);
Step 5: capacitance when being zero from high frequency or low frequency C-V curve reading applied voltage, sets up formula 3 and formula 4;
1 C 0 = 2 C i + 1 C s 10 + 1 C s 20 Formula 3
( V s 10 + V s 20 ) - 2 V fermi = Q s 10 + Q s 20 C i + Q f 1 + Q f 2 C i Formula 4
Wherein, C 0for capacitance when applied voltage is zero, C ibe capacitive dielectric layer, its value is: C i> > e i/ t i.C s10with C s20be respectively lower semiconductor layer differential capacitance, they are about V s10, V s20function, V s10with V s20for applied voltage be zero time, upper and lower semiconductor layer surface gesture;
C s 10 = ϵ s L D - exp ( - βV s 10 ) + 1 + n p 0 p p 0 ( exp ( βV s 10 ) - 1 ) F ( βV s 10 , n p 0 / p p 0 )
C s 20 = ϵ s L D - exp ( - βV s 20 ) + 1 + n p 0 p p 0 ( exp ( βV s 20 ) - 1 ) F ( βV s 20 , n p 0 / p p 0 )
Wherein, n p0and p p0equilibrium electron concentration and the hole concentration of P-type silicon respectively, L dbe the Debye length of semiconductor layer, the expression formula of F (a, b) is: q f1, Q f2be the fixed charge density of the first insulating barrier (2), the second insulating barrier (4), V fermifor Fermi potential.Q s10with Q s20be respectively the charge density of the first semiconductor layer (1), the second semiconductor layer (5), they are about V s10, V s20function;
Q s 10 = - sign ( V s 10 ) 2 ϵ s β L D F ( βV s 10 , n p 0 / p p 0 ) / q
Q s 20 = - sign ( V s 20 ) 2 ϵ s β L D F ( βV s 20 , n p 0 / p p 0 ) / q
Step 6: read feature extreme point (V from high frequency or low frequency C-V curve g *, C *), set up formula 5 and formula 6;
Q f 1 + Q f 2 = 4 ϵ s βL D F ( βV s * , n p 0 / p p 0 ) - Q m 0 Formula 5
Q f 1 - Q f 2 = V g * C i Formula 6
Wherein and Q m0tried to achieve by following formula:
ϵ s L D - exp ( - βV s * ) + 1 + n p 0 p p 0 ( exp ( βV s * ) - 1 ) F ( βV s * , n p 0 / p p 0 ) = 2 C * C i C i - 2 C *
Q m0=C i(2V fermi-V s01-V s02);
Step 7: simultaneous formula 3,4,5,6, utilizes Program to record the fixed charge density of the first insulating barrier (2), the second insulating barrier (4).
2. chip-stack structural parameters method of testing as claimed in claim 1, is characterized in that: the dopant material of the first semiconductor layer (1), the second semiconductor layer (5) is Si, Ge, GaAs or InP.
3. chip-stack structural parameters method of testing as claimed in claim 1, is characterized in that: the material of the first insulating barrier (2), the second insulating barrier (4) is SiO 2, SiN or Si 3n 4dielectric material.
4. chip-stack structural parameters method of testing as claimed in claim 1, is characterized in that: the material of metal level (3) is Cu, TiN or Ti.
5. chip-stack structural parameters method of testing as claimed in claim 1, is characterized in that: the doping type of the first semiconductor layer (1), the second semiconductor layer (5) is identical, is namely P type or N-type.
6. chip-stack structural parameters method of testing as claimed in claim 1, is characterized in that: the doping content of the first semiconductor layer (1), the second semiconductor layer (5) is identical or not identical.
7. chip-stack structural parameters method of testing as claimed in claim 1, is characterized in that: the fixed charge density in the first insulating barrier (2), the second insulating barrier (4) is identical or not identical.
CN201510089321.8A 2015-02-26 2015-02-26 A kind of measuring method of chip-stack structural parameters Active CN104658941B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510089321.8A CN104658941B (en) 2015-02-26 2015-02-26 A kind of measuring method of chip-stack structural parameters

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510089321.8A CN104658941B (en) 2015-02-26 2015-02-26 A kind of measuring method of chip-stack structural parameters

Publications (2)

Publication Number Publication Date
CN104658941A true CN104658941A (en) 2015-05-27
CN104658941B CN104658941B (en) 2018-03-30

Family

ID=53249913

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510089321.8A Active CN104658941B (en) 2015-02-26 2015-02-26 A kind of measuring method of chip-stack structural parameters

Country Status (1)

Country Link
CN (1) CN104658941B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110531243A (en) * 2019-10-16 2019-12-03 启发(天津)电子科技有限公司 A kind of semiconducter device testing method and system
CN112768367A (en) * 2021-02-04 2021-05-07 微龛(广州)半导体有限公司 SOI wafer bonding quality detection method and system
CN112951735A (en) * 2021-01-28 2021-06-11 微龛(广州)半导体有限公司 Wafer bonding quality detection method and system

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5485091A (en) * 1995-05-12 1996-01-16 International Business Machines Corporation Contactless electrical thin oxide measurements
US20040251927A1 (en) * 2003-06-10 2004-12-16 Alex Salnik Characterization of ultra shallow junctions in semiconductor wafers
CN101136347A (en) * 2007-09-29 2008-03-05 上海集成电路研发中心有限公司 MOS pipe interface state testing method
CN101159243A (en) * 2007-08-20 2008-04-09 中国科学院光电技术研究所 Method for measuring semiconductor doping concentration

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5485091A (en) * 1995-05-12 1996-01-16 International Business Machines Corporation Contactless electrical thin oxide measurements
US20040251927A1 (en) * 2003-06-10 2004-12-16 Alex Salnik Characterization of ultra shallow junctions in semiconductor wafers
CN101159243A (en) * 2007-08-20 2008-04-09 中国科学院光电技术研究所 Method for measuring semiconductor doping concentration
CN101136347A (en) * 2007-09-29 2008-03-05 上海集成电路研发中心有限公司 MOS pipe interface state testing method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
MAN LI等: "Capacitance-Voltage Characteristics of Symmetric Stacked Bonding Structure in 3D Integration", 《ECS TRANSACTIONS》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110531243A (en) * 2019-10-16 2019-12-03 启发(天津)电子科技有限公司 A kind of semiconducter device testing method and system
CN112951735A (en) * 2021-01-28 2021-06-11 微龛(广州)半导体有限公司 Wafer bonding quality detection method and system
CN112951735B (en) * 2021-01-28 2022-11-08 微龛(广州)半导体有限公司 Wafer bonding quality detection method and system
CN112768367A (en) * 2021-02-04 2021-05-07 微龛(广州)半导体有限公司 SOI wafer bonding quality detection method and system

Also Published As

Publication number Publication date
CN104658941B (en) 2018-03-30

Similar Documents

Publication Publication Date Title
CN101136347A (en) MOS pipe interface state testing method
CN108318796B (en) Three-port silicon carbide-based power device interface state testing method
CN104658941A (en) Method for measuring parameters of chip lamination structure
CN102522386B (en) Gate-oxidizing-layer interface-trap density-testing structure and testing method
CN103471740B (en) A kind of capacitive temperature sensor
CN102832203A (en) Structure and method for testing trap density of gate oxide interface
Arefin et al. A MEMS capacitive pH sensor for high acidic and basic solutions
Recart et al. Application of junction capacitance measurements to the characterization of solar cells
CN106247920B (en) A kind of surface strain detection device based on the sandwich interdigital capacitor of elastic substrates
CN105136352A (en) Capacitive pressure sensor and preparation method thereof
JP5881037B2 (en) Method and apparatus for evaluating insulator / semiconductor interface
CN101635292B (en) Contact pad for measuring electrical thickness of gate dielectric layer and measurement structure thereof
CN103943606B (en) Detection structure and detection method for defect in semiconductor device
Tikhov et al. Capacitors with nonlinear characteristics based on stabilized zirconia with built-in gold nanoparticles
CN104716065B (en) Capacitance-voltage characteristic correction method for metal oxide semiconductor field-effect transistor
CN103822948A (en) Testing method for semiconductor
Guenther et al. Extracting accurate capacitance voltage curves from impedance spectroscopy
CN104457548B (en) A kind of non-contact type silicon wafer thickness thinning measurement apparatus
CN106021667B (en) Method for extracting structure parameters of silicon on insulating substrate and bulk silicon transverse power diode
CN100356182C (en) Quick characterization method for charge density of buried oxide layer of SOI
CN1687800A (en) Electrical parameter characterization method of silicon on insulation body
CN103868952B (en) Ion implanted layer carrier concentration method of testing
CN109884122A (en) Organic gas detection chip based on sulfuration rhenium nano-device
CN102338767A (en) Electrolyte-insulator-semiconductor (EIS) structure based on H<+> sensitive Ta2O5 material, preparation thereof and pH online detection method
CN103928364B (en) A kind of for detecting the structure of alloying level in eutectic bonding

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20181212

Address after: 226000 Room 8319, Building 11, Happy New Town, Gangzha District, Nantong City, Jiangsu Province

Patentee after: Nanjing University of Posts and Telecommunications Nantong Institute Limited

Address before: 210023 9 Wen Yuan Road, Ya Dong new town, Qixia District, Nanjing, Jiangsu.

Patentee before: Nanjing Post & Telecommunication Univ.